SEARCH · Geology Reports

Results for “Processes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47Linked to original sources

USGS2010

Processes of Terrace Formation on the Piedmont of the Santa Cruz River Valley During Quaternary Time, Green Valley-Tubac Area, Southeastern Arizona

In this report we describe a series of stepped Quaternary terraces on some piedmont tributaries of the Santa Cruz River valley in southeastern Arizona. These terraces began to form in early Pleistocene time, after major basin-and-range faulting ceased, with lateral planation of basin fill and deposition of thin fans of alluvium. At the end of this cycle of erosion and deposition, tributaries of the Santa Cruz River began the process of dissection and terrace formation that continues to the present. Vertical cutting alternated with periods of equilibrium, during which streams cut laterally and left thin deposits of channel fill. The distribution of terraces was mapped and compiled with adjacent mapping to produce a regional picture of piedmont stream history in the middle part of the Santa Cruz River valley. For selected tributaries, the thickness of terrace fill was measured, particle size and lithology of gravel were determined, and sedimentary features were photographed and described. Mapping of terrace stratigraphy revealed that on two tributaries, Madera Canyon Wash and Montosa Canyon Wash, stream piracy has played an important role in piedmont landscape development. On two other tributaries, Cottonwood Canyon Wash and Josephine Canyon Wash, rapid downcutting preempted piracy. Two types of terraces are recognized: erosional and depositional. Gravel in thin erosional terraces has Trask sorting coefficients and sedimentary structures typical of streamflood deposits, replete with bar-and-swale surface topography on young terraces. Erosional-terrace fill represents the channel fill of the stream that cuts the terrace; the thickness of the fill indicates the depth of channel scour. In contrast to erosional terraces, depositional terraces show evidence of repeated deposition and net aggradation, as indicated by their thickness (as much as 20+ m) and weakly bedded structure. Depositional terraces are common below mountain-front canyon mouths where streams drop their load in response to abrupt flattening of gradients and expansion of channel banks, and they extend down the piedmont along Josephine Canyon Wash. Gravel in depositional terraces also has sorting coefficients typical of streamflood deposits. Sedimentary features in both types of terraces are consistent with deposition by flash floods in ephemeral streams, suggesting the climate was arid. Bedding and clast armor are weakly developed, clast clusters and imbrication are common, and crossbedding is generally absent. Debris-flow deposits, even near the mountain front, are surprisingly rare. On the tectonically stable piedmont of southeastern Arizona, stream piracy and climate change are the most likely agents of terrace formation. Both piracy and climate change can cause rapid changes in discharge and sediment supply, which initiate cycles of incision, lateral cutting, and aggradation. Increased stream discharge initiates downcutting, but increased sediment supply interrupts downcutting and causes streams to cut laterally and aggrade. At times, on Madera Canyon Wash and Montosa Canyon Wash, stream piracy affected stream discharge and sediment supply, but on Cottonwood Canyon Wash and Josephine Canyon Wash, only climate change could have initiated terrace cutting. Terraces probably formed during extended arid intervals when sparse vegetation and flashy stream discharge combined to increase sediment supply. In most cases, sediment supply was sufficient to promote lateral cutting but not long-term aggradation. Thus, most streams formed erosional terraces. The middle Pleistocene Josephine Canyon Wash formed a depositional terrace because it had a source of abundant unconsolidated sediment.

Scientific Investigations Report↗

USGS2011

Digital signal processing and interpretation of full waveform sonic log for well BP-3-USGS, Great Sand Dunes National Park and Preserve, Alamosa County, Colorado

Along the Great Sand Dunes National Park and Preserve boundary (fig. 1), 10 monitoring wells were drilled by the National Park Service in order to monitor water flow in an unconfined aquifer spanning the park boundary. Adjacent to the National Park Service monitoring well named Boundary Piezometer Well No. 3, or BP-3, the U.S. Geological Survey (USGS) drilled the BP-3-USGS well. This well was drilled from September 14 through 17, 2009, to a total depth of 99.4 meters (m) in order to acquire additional subsurface information. The BP-3-USGS well is located at lat 37 degrees 43'18.06' and long -105 degrees 43'39.30' at a surface elevation of 2,301 m. Approximately 23 m of core was recovered beginning at a depth of 18 m. Drill cuttings were also recovered. The wireline geophysical logs acquired in the well include natural gamma ray, borehole caliper, temperature, full waveform sonic, density, neutron, resistivity, and induction logs. The BP-3-USGS well is now plugged and abandoned. This report details the full waveform digital signal processing methodology and the formation compressional-wave velocities determined for the BP-3-USGS well. These velocity results are compared to several velocities that are commonly encountered in the subsurface. The density log is also discussed in context of these formation velocities.

Scientific Investigations Report↗

USGS2011

Vegetation of the Elwha River estuary: Chapter 8 in Coastal habitats of the Elwha River, Washington--biological and physical patterns and processes prior to dam removal

The Elwha River estuary supports one of the most diverse coastal wetland complexes yet described in the Salish Sea region, in terms of vegetation types and plant species richness. Using a combination of aerial imagery and vegetation plot sampling, we identified 6 primary vegetation types and 121 plant species in a 39.7 ha area. Most of the estuary is dominated by woody vegetation types, with mixed riparian forest being the most abundant (20 ha), followed by riparian shrub (6.3 ha) and willow-alder forest (3.9 ha). The shrub-emergent marsh transition vegetation type was fourth most abundant (2.2 ha), followed by minor amounts of dunegrass (1.75 ha) and emergent marsh (0.2 ha). This chapter documents the abundance, distribution, and floristics of these six vegetation types, including plant species richness, life form, species origin (native or introduced), and species wetland indicator status. These data will serve as a baseline to which future changes can be compared, following the impending removal of Glines Canyon and Elwha Dams upstream on the Elwha River. Dam removals may alter many of the processes, materials, and biotic interactions that influence the estuary plant communities, including hydrology, salinity, sediment and wood transport, nutrients, and plant-microbe interactions.

Washington↗

USGS2013

Streamflow and water-quality conditions including geologic sources and processes affecting selenium loading in the Toll Gate Creek watershed, Aurora, Arapahoe County, Colorado, 2007

Toll Gate Creek is a perennial stream draining a suburban area in Aurora, Colorado, where selenium concentrations have consistently exceeded the State of Colorado aquatic-life standard for selenium of 4.6 micrograms per liter since the early 2000s. In cooperation with the City of Aurora, Colorado, Utilities Department, a synoptic water-quality study was performed along an 18-kilometer reach of Toll Gate Creek extending from downstream from Quincy Reservoir to the confluence with Sand Creek to develop a detailed understanding of streamflow and concentrations and loads of selenium in Toll Gate Creek. Streamflow and surface-water quality were characterized for summer low-flow conditions (July–August 2007) using four spatially overlapping synoptic-sampling subreaches. Mass-balance methods were applied to the synoptic-sampling and tracer-injection results to estimate streamflow and develop spatial profiles of concentration and load for selenium and other chemical constituents in Toll Gate Creek surface water. Concurrent groundwater sampling determined concentrations of selenium and other chemical constituents in groundwater in areas surrounding the Toll Gate Creek study reaches. Multivariate principal-component analysis was used to group samples and to suggest common sources for dissolved selenium and major ions. Hydrogen and oxygen stable-isotope ratios, groundwater-age interpretations, and chemical analysis of water-soluble paste extractions from core samples are presented, and interpretation of the hydrologic and geochemical data support conclusions regarding geologic sources of selenium and the processes affecting selenium loading in the Toll Gate Creek watershed.

Colorado↗

USGS2013

Arsenic in New Jersey Coastal Plain streams, sediments, and shallow groundwater: effects from different geologic sources and anthropogenic inputs on biogeochemical and physical mobilization processes

Arsenic (As) concentrations in New Jersey Coastal Plain streams generally exceed the State Surface Water Quality Standard (0.017 micrograms per liter (µg/L)), but concentrations seldom exceed 1 µg/L in filtered stream-water samples, regardless of geologic contributions or anthropogenic inputs. Nevertheless, As concentrations in unfiltered stream water indicate substantial variation because of particle inputs from soils and sediments with differing As contents, and because of discharges from groundwater of widely varying chemistry. In the Inner Coastal Plain, streams draining to lower reaches of the Delaware River traverse As-rich glauconitic sediments of marine origin in which As contents typically are about 20 milligrams per kilogram (mg/kg) or greater. In some of these sedimentary units, As concentrations exceed the New Jersey drinking-water maximum contaminant level (5 µg/L) in shallow groundwater that discharges to streams. Microbes, fueled by organic carbon beneath the streambed, reduce iron (Fe) and As, releasing As and Fe into solution in the shallow groundwater from geologic materials that likely include (in addition to glauconite) other phyllosilicates, apatite, and siderite. When the groundwater discharges to the stream, the dissolved Fe and As are oxidized, the Fe precipitates as a hydroxide, and the As sorbs or co-precipitates with the Fe. Because of the oxidation/precipitation process, dissolved As concentrations measured in filtered stream waters of the Inner Coastal Plain are about 1 µg/L, but the total As concentrations (and loads) are greater, substantially amplified by As-bearing suspended sediment in stormflows. In the Outer Coastal Plain, streams draining to the Atlantic Ocean traverse quartz-rich sediments of mainly deltaic origin where the As content generally is low (<8 mg/kg). In unfiltered and filtered water samples, As concentrations typically are less than 1 µg/L in the acidic stream water and groundwater of the Outer Coastal Plain, but are greater in waters from urban areas. Despite the generally small geologic contributions to Outer Coastal Plain groundwater, where wastewater inputs were indicated, concentrations of As in unfiltered shallow groundwater discharging to small urban streams exceeded the maximum contaminant level. With a history of agriculture in the New Jersey Coastal Plain, anthropogenic inputs of As, such as residues from former pesticide applications in soils, can amplify any geogenic As in runoff. Such inputs contribute to an increased total As load to a stream at high stages of flow. As a result of yet another anthropogenic influence, microbes that reduce and mobilize As beneath the streambeds are stimulated by inputs of dissolved organic carbon (DOC). Although DOC is naturally occurring, anthropogenic contributions from wastewater inputs may deliver increased levels of DOC to subsurface soils and ultimately groundwater. Arsenic concentrations may increase with the increases in pH of groundwater and stream water in developed areas receiving wastewater inputs, as As mobilization caused by pH-controlled sorption and desorption reactions are likely to occur in waters of neutral or alkaline pH (for example, Nimick and others, 1998; Barringer and others, 2007b). Because of the difference in As content of the geologic materials in the two sub-provinces of the Coastal Plain, the amount of As that is mobile in groundwater and stream water is, potentially, substantially greater in the Inner Coastal Plain than in the Outer Coastal Plain. In turn, streams within the Inner and Outer Coastal Plain can receive substantially more As in groundwater discharge from developed areas than from environments where DOC appears to be of natural origin.

New Jersey↗

USGS2014

Assessment of the quality of groundwater and the Little Wind River in the area of a former uranium processing facility on the Wind River Reservation, Wyoming, 1987 through 2010

In 2010, the U.S Geological Survey (USGS), in cooperation with the Wind River Environmental Quality Commission (WREQC), began an assessment of the effectiveness of the existing monitoring network at the Riverton, Wyoming, Uranium Mill Tailings Remedial Action (UMTRA) site. The USGS used existing data supplied by the U.S. Department of Energy (DOE). The study was to determine (1) seasonal variations in the direction of groundwater flow in the area of the former uranium processing facility toward the Little Wind River, (2) the extent of contaminated groundwater among the aquifers and between the aquifers and the Little Wind River, (3) whether current monitoring is adequate to establish the effectiveness of natural attenuation for the contaminants of concern, and (4) the influence of groundwater discharged from the sulfuric-acid plant on water quality in the Little Wind River.

Wyoming↗

USGS2019

Field trip guide to Mount St. Helens, Washington—Recent and ancient volcaniclastic processes and deposits

This field guide explores volcanic effusions, sediments, and landforms at Mount St. Helens in Washington. A detailed synopsis outlines the eruptive history of Mount St. Helens from about 300,000 years ago through 1980 and beyond. The five days in the field include about 28 stops and 12 potential stops. Exposures in valleys surrounding Mount St. Helens reveal records of diverse Pleistocene and Holocene processes including debris avalanche, lahar, huge water wave on a nearby lake, pyroclastic density currents (surge and flow), tephra fall, lava flow, the growth of domes, and Pleistocene glaciation. Many of the stops explore effects of the several catastrophes that constituted the 18 May 1980 eruption and made Mount St. Helens famous.

Washington↗

USGS2017

A process to estimate net infiltration using a site-scale water-budget approach, Rainier Mesa, Nevada National Security Site, Nevada, 2002–05

This report documents a process used to estimate net infiltration from precipitation, evapotranspiration (ET), and soil data acquired at two sites on Rainier Mesa. Rainier Mesa is a groundwater recharge area within the Nevada National Security Site where recharged water flows through bedrock fractures to a deep (450 meters) water table. The U.S. Geological Survey operated two ET stations on Rainier Mesa from 2002 to 2005 at sites characterized by pinyon-juniper and scrub-brush vegetative cover. Precipitation and ET data were corrected to remove measurement biases and gap-filled to develop continuous datasets. Net infiltration (percolation below the root zone) and changes in root-zone water storage were estimated using a monthly water-balance model. Site-scale water-budget results indicate that the heavily-fractured welded-tuff bedrock underlying thin (<40 centimeters) topsoil is a critical water source for vegetation during dry periods. Annual precipitation during the study period ranged from fourth lowest (182 millimeters [mm]) to second highest (708 mm) on record (record = 55 years). Annual ET exceeded precipitation during dry years, indicating that the fractured-bedrock reservoir capacity is sufficient to meet atmospheric-evaporative demands and to sustain vegetation through extended dry periods. Net infiltration (82 mm) was simulated during the wet year after the reservoir was rapidly filled to capacity. These results support previous conclusions that preferential fracture flow was induced, resulting in an episodic recharge pulse that was detected in nearby monitoring wells. The occurrence of net infiltration only during the wet year is consistent with detections of water-level rises in nearby monitoring wells that occur only following wet years.

Nevada↗

USGS2017

Modern landscape processes affecting archaeological sites along the Colorado River corridor downstream of Glen Canyon Dam, Glen Canyon National Recreation Area, Arizona

The landscape of the Colorado River through Glen Canyon National Recreation Area formed over many thousands of years and was modified substantially after the completion of Glen Canyon Dam in 1963. Changes to river flow, sediment supply, channel base level, lateral extent of sedimentary terraces, and vegetation in the post-dam era have modified the river-corridor landscape and have altered the effects of geologic processes that continue to shape the landscape and its cultural resources. The Glen Canyon reach of the Colorado River downstream of Glen Canyon Dam hosts many archaeological sites that are prone to erosion in this changing landscape. This study uses field evaluations from 2016 and aerial photographs from 1952, 1973, 1984, and 1996 to characterize changes in potential windblown sand supply and drainage configuration that have occurred over more than six decades at 54 archaeological sites in Glen Canyon and uppermost Marble Canyon. To assess landscape change at these sites, we use two complementary geomorphic classification systems. The first evaluates the potential for aeolian (windblown) transport of river-derived sand from the active river channel to higher elevation archaeological sites. The second identifies whether rills, gullies, or arroyos (that is, overland drainages that erode the ground surface) exist at the archaeological sites as well as the geomorphic surface, and therefore the relative base level, to which those flow paths drain. Results of these assessments are intended to aid in the management of irreplaceable archaeological resources by the National Park Service and stakeholders of the Glen Canyon Dam Adaptive Management Program.

Arizona↗

USGS2018

Water-quality observations of the San Antonio segment of the Edwards aquifer, Texas, with an emphasis on processes influencing nutrient and pesticide geochemistry and factors affecting aquifer vulnerability, 2010–16

As questions regarding the influence of increasing urbanization on water quality in the Edwards aquifer are raised, a better understanding of the sources, fate, and transport of compounds of concern in the aquifer—in particular, nutrients and pesticides—is needed to improve water management decision-making capabilities. The U.S. Geological Survey, in cooperation with the San Antonio Water System, performed a study from 2010 to 2016 to better understand how water quality changes under a range of hydrologic conditions and in contrasting land-cover settings (rural and urban) in the Edwards aquifer. The study design included continuous hydrologic monitoring, continuous water-quality monitoring, and discrete sample collection for a detailed characterization of water quality at a network of sites throughout the aquifer system. The sites were selected to encompass a “source-to-sink” (that is, from aquifer recharge to aquifer discharge) approach. Network sites were selected to characterize rainfall, recharging surface water, and groundwater; groundwater sites included wells in the unconfined part of the aquifer (unconfined wells) and in the confined part of the aquifer (confined wells) and a major discharging spring. Storm-related samples—including rainfall samples, stormwater-runoff (surface-water) samples, and groundwater samples—were collected to characterize the aquifer response to recharge. Elevated nitrate concentrations relative to national background values and the widespread detection of pesticides indicate that the Edwards aquifer is vulnerable to contamination and that vulnerability is affected by factors such as land cover, aquifer hydrogeology, and changes in hydrologic conditions. Greater vulnerability of groundwater in urban areas relative to rural areas was evident from results for urban groundwater sites, which generally had higher nitrate concentrations, elevated δ 15 N-nitrate values, a greater diversity of pesticides, and higher pesticide concentrations. The continuum of water quality from unconfined rural groundwater sites (least affected by anthropogenic contamination) to unconfined urban groundwater sites (most affected by anthropogenic contamination) demonstrates enhanced vulnerability of urban versus rural land cover. Differences in contaminant occurrences and concentration among unconfined urban wells indicate that the urban parts of the aquifer are not uniformly vulnerable, but rather are affected by spatial differences in the sources of nutrients and pesticides. In urban areas, the shallow, unconfined groundwater sites showed greater temporal variability in both nutrient and pesticide concentrations, as well as a greater degree of contamination, than did deeper, confined groundwater sites. In comparison to that of the shallow, unconfined groundwater sites, the water quality of the deeper, confined groundwater sites was relatively invariant during this multiyear study. Although aquifer hydrogeology is an important factor related to aquifer vulnerability, land cover likely has a greater influence on pesticide contamination of groundwater. Temporal variability in hydrologic conditions for the Edwards aquifer is apparent in data for surface water as a source of groundwater recharge, water-level altitude in wells, spring discharge, and groundwater quality. This temporal variability affects recharge sources, recharge amounts, groundwater traveltimes, flow routing, water-rock interaction processes, dilution, mixing, and, in turn, water quality. Relations of land cover, aquifer hydrogeology, and changing hydrologic conditions to water quality are complex but provide insight into the vulnerability of Edwards aquifer groundwater—a vital drinking-water resource.

Texas↗

USGS2021

Procedures for field data collection, processing, quality assurance and quality control, and archiving of relative- and absolute-gravity surveys

Repeat microgravity surveys carried out using relative- and absolute-gravity meters are useful for identifying changes in subsurface mass, such as the volume of water stored in an aquifer. These surveys require careful field procedures to achieve the part-per-billion accuracy required to measure the small changes in gravity relevant for hydrologic studies. This chapter describes techniques and methods for carrying out gravity surveys, requirements for assuring high-quality survey results, and data processing and archival procedures. The focus is on acquiring and documenting repeat gravity surveys for monitoring changes in groundwater storage. Similar gravity surveys may be completed to evaluate other causes of mass change, such as those caused by magma movement below volcanoes. The methods are also useful for one-time surveys that map spatial gravity variations associated with geologic structures such as faults or sedimentary basins. Repeat microgravity surveys can be carried out using relative-gravity meters, absolute-gravity meters, or both. Specific locations, known as gravity stations, are visited during each survey. Most commonly, absolute- and relative-gravity are combined using the least-squares method of network adjustment, much like benchmark elevations and relative-height differences in a leveling network. This chapter primarily describes the use of the A-10 absolute-gravity meter manufactured by Micro-g LaCoste, Inc., and relative-gravity meters made by LaCoste & Romberg (no longer in production) and ZLS Corporation, Inc. Field and office procedures are similar for other instruments such as the FG-5 absolute-gravity meter and Scintrex relative-gravity meters, but some adaptation may be required. Quality control for absolute-gravity data focuses primarily on proper field procedures and maintaining the time and distance calibration of the instrument. Quality control for relative-gravity surveys requires careful field procedures, an understanding of how the meter is behaving while in the field, and appropriate postprocessing. The techniques and methods described in this chapter were developed over 30 years at the USGS Arizona Water Science Center and the Southwest Gravity Program and are the basis for many studies on groundwater-storage change and geologic structure. A description of the Program and complete bibliography is available at https://www.usgs.gov/centers/az-water/science/azwsc-capabilities-hydrologic-gravity-monitoring .

Techniques and Methods↗

USGS2006

Use of the Multi-Node Well (MNW) package when simulating solute transport with the MODFLOW ground-water transport process

This report describes modifications to a U.S. Geological Survey (USGS) three-dimensional solute-transport model (MODFLOW-GWT), which is incorporated into the USGS MODFLOW ground-water model as the Ground-Water Transport (GWT) Process. The modifications were made to create compatibility between the Multi-Node Well (MNW) Package for MODFLOW and the MODFLOW-GWT model. This compatibility improves the capability of MODFLOW-GWT to represent accurately solute transport in simulations that include multi-node wells because long-screen wells or long open boreholes that extend through multiple model layers can provide fast pathways for solutes to move from one location to another in a ground-water flow system. For nonpumping multi-node wells (used to simulate open boreholes or observation wells, for example), a simple routing and local mixing model was developed to calculate nodal concentrations within the borehole. A depth-averaged concentration is calculated for such nonpumping wells. For pumping multi-node wells (either withdrawal or injection) in which the flow between the well and the ground-water system is in the same direction at all nodes, the average concentration in the well is calculated as a flux-based mean assuming complete and instantaneous mixing in the wellbore of all inflows. For pumping multi-node wells (either withdrawal or injection) in which the flow between the well and the ground-water system is not unidirectional, the concentration distribution within the well is calculated using the same routing and local mixing model used for a nonpumping multi-node well, with the added assumption that the flux pumped in or out of the well is added or removed above the first well node.

Techniques and Methods↗

USGS2006

User guide for the farm process (FMP1) for the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model, MODFLOW-2000

There is a need to estimate dynamically integrated supply-and-demand components of irrigated agriculture as part of the simulation of surface-water and ground-water flow. To meet this need, a computer program called the Farm Process (FMP1) was developed for the U.S. Geological Survey three-dimensional finite-difference modular ground-water flow model, MODFLOW- 2000 (MF2K). The FMP1 allows MF2K users to simulate conjunctive use of surface- and ground water for irrigated agriculture for historical and future simulations, water-rights issues and operational decisions, nondrought and drought scenarios. By dynamically integrating farm delivery requirement, surface- and ground-water delivery, as well as irrigation-return flow, the FMP1 allows for the estimation of supplemental well pumpage. While farm delivery requirement and irrigation return flow are simulated by the FMP1, the surface-water delivery to the farm can be simulated optionally by coupling the FMP1 with the Streamflow Routing Package (SFR1) and the farm well pumping can be simulated optionally by coupling the FMP1 to the Multi-Node Well (MNW) Package. In addition, semi-routed deliveries can be specified that are associated with points of diversion in the SFR1 stream network. Nonrouted surface-water deliveries can be specified independently of any stream network. The FMP1 maintains a dual mass balance of a farm budget and as part of the ground-water budget. Irrigation demand, supply, and return flow are in part subject to head-dependent sources and sinks such as evapotranspiration from ground water and leakage between the conveyance system and the aquifer. Farm well discharge and farm net recharge are source/sink terms in the FMP1, which depend on transpiration uptake from ground water and other head dependent consumptive use components. For heads rising above the bottom of the root zone, the actual transpiration is taken to vary proportionally with the depth of the active root zone, which can be restricted by anoxia or wilting. Depths corresponding to anoxia- or wilting-related pressure heads within the root zone are found using analytical solutions of a vertical pseudo steady-state pressure- head distribution over the depth of the total root zone (Consumptive Use Concept 1). Alternatively, a simpler, conceptual model is available, which defines how consumptive use (CU) components vary with changing head (CU Concept 2). Subtracting the ground water and precipitation transpiration components from the total transpiration yields a transpiratory irrigation requirement for each cell. The total farm delivery requirement (TFDR) then is determined as cumulative transpiratory and evaporative irrigation requirements of all farm cells and increased sufficiently to compensate for inefficient use from irrigation with respect to plant consumption. The TFDR subsequently is satisfied with surface- and ground-water delivery, respectively constrained by allotments, water rights, or maximum capacities. Five economic and noneconomic drought response policies can be applied optionally, if the potential supply of surface water and ground water is insufficient to meet the crop demand: acreage-optimization with or without a water conservation pool, deficit irrigation with or without water-stacking, and zero policy.

Techniques and Methods↗

USGS2007

Documentation of a Conduit Flow Process (CFP) for MODFLOW-2005

This report documents the Conduit Flow Process (CFP) for the modular finite-difference ground-water flow model, MODFLOW-2005. The CFP has the ability to simulate turbulent ground-water flow conditions by: (1) coupling the traditional ground-water flow equation with formulations for a discrete network of cylindrical pipes (Mode 1), (2) inserting a high-conductivity flow layer that can switch between laminar and turbulent flow (Mode 2), or (3) simultaneously coupling a discrete pipe network while inserting a high-conductivity flow layer that can switch between laminar and turbulent flow (Mode 3). Conduit flow pipes (Mode 1) may represent dissolution or biological burrowing features in carbonate aquifers, voids in fractured rock, and (or) lava tubes in basaltic aquifers and can be fully or partially saturated under laminar or turbulent flow conditions. Preferential flow layers (Mode 2) may represent: (1) a porous media where turbulent flow is suspected to occur under the observed hydraulic gradients; (2) a single secondary porosity subsurface feature, such as a well-defined laterally extensive underground cave; or (3) a horizontal preferential flow layer consisting of many interconnected voids. In this second case, the input data are effective parameters, such as a very high hydraulic conductivity, representing multiple features. Data preparation is more complex for CFP Mode 1 (CFPM1) than for CFP Mode 2 (CFPM2). Specifically for CFPM1, conduit pipe locations, lengths, diameters, tortuosity, internal roughness, critical Reynolds numbers (NRe), and exchange conductances are required. CFPM1, however, solves the pipe network equations in a matrix that is independent of the porous media equation matrix, which may mitigate numerical instability associated with solution of dual flow components within the same matrix. CFPM2 requires less hydraulic information and knowledge about the specific location and hydraulic properties of conduits, and turbulent flow is approximated by modifying horizontal conductances assembled by the Block-Centered Flow (BCF), Layer-Property Flow (LPF), or Hydrogeologic-Unit Flow Packages (HUF) of MODFLOW-2005. For both conduit flow pipes (CFPM1) and preferential flow layers (CFPM2), critical Reynolds numbers are used to determine if flow is laminar or turbulent. Due to conservation of momentum, flow in a laminar state tends to remain laminar and flow in a turbulent state tends to remain turbulent. This delayed transition between laminar and turbulent flow is introduced in the CFP, which provides an additional benefit of facilitating convergence of the computer algorithm during iterations of transient simulations. Specifically, the user can specify a higher critical Reynolds number to determine when laminar flow within a pipe converts to turbulent flow, and a lower critical Reynolds number for determining when a pipe with turbulent flow switches to laminar flow. With CFPM1, the Hagen-Poiseuille equation is used for laminar flow conditions and the Darcy-Weisbach equation is applied to turbulent flow conditions. With CFPM2, turbulent flow is approximated by reducing the laminar hydraulic conductivity by a nonlinear function of the Reynolds number, once the critical head difference is exceeded. This adjustment approximates the reductions in mean velocity under turbulent ground-water flow conditions.

Techniques and Methods↗

USGS2009

The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1

The ability to dynamically simulate the integrated supply-and-demand components of irrigated agricultural is needed to thoroughly understand the interrelation between surface water and groundwater flow in areas where the water-use by vegetation is an important component of the water budget. To meet this need, the computer program Farm Process (FMP1) was updated and refined for use with the U.S. Geological Survey's MODFLOW-2005 groundwater-flow model, and is referred to as MF2005-FMP2. The updated program allows the simulation, analysis, and management of nearly all components of human and natural water use. MF2005-FMP2 represents a complete hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for water consumption of irrigated agriculture, but also of urban use, and of natural vegetation. Supply and demand components of water use are analyzed under demand-driven and supply-constrained conditions. From large- to small-scale settings, the MF2005-FMP2 has the unique set of capabilities to simulate and analyze historical, present, and future conditions. MF2005-FMP2 facilitates the analysis of agricultural water use where little data is available for pumpage, land use, or agricultural information. The features presented in this new version of FMP2 along with the linkages to the Streamflow Routing (SFR), Multi-Node Well (MNW), and Unsaturated Zone Flow (UZF) Packages prevents mass loss to an open system and helps to account for 'all of the water everywhere and all of the time'. The first version, FMP1 for MODFLOW-2000, is limited to (a) transpiration uptake from unsaturated root zones, (b) on-farm efficiency defined solely by farm and not by crop type, (c) a simulation of water use and returnflows related only to irrigated agriculture and not also to non-irrigated vegetation, (d) a definition of consumptive use as potential crop evapotranspiration, (e) percolation being instantly recharged to the uppermost active aquifer, (f) automatic routing of returnflow from runoff either to reaches of tributary stream segments adjacent to a farm or to one reach nearest to the farm's lowest elevation, (g) farm-well pumping from cell locations regardless of whether an irrigation requirement from these cells exists or not, and (h) specified non-routed water transfers from an undefined source outside the model domain. All of these limitations are overcome in MF2005-FMP2. The new features include (a) simulation of transpiration uptake from variably saturated, fully saturated, or ponded root zones (for example, for crops like rice or riparian vegetation), (b) definition of on-farm efficiency not only by farm but also by crop, (c) simulation of water use and returnflow from non-irrigated vegetation (for example, rain-fed agriculture or native vegetation), (d) use of crop coefficients and reference evapotranspiration, (e) simulation of the delay between percolation from farms through the unsaturated zone and recharge into the uppermost active aquifer by linking FMP2 to the UZF Package, (f) an option to manually control the routing of returnflow from farm runoff to streams, (g) an option to limit pumping to wells located only in cells where an irrigation requirement exists, and (h) simulation of water transfers to farms from a series of well fields (for example, recovery well field of an aquifer-storage-and-recovery system, ASR). In addition to the output of an economic budget for each farm between irrigation demand and supply ('Farm Demand and Supply Budget' in FMP1), a new output option called 'Farm Budget' was created for FMP2, which allows the user to track all physical flows into and out of a water accounting unit at all times. Such a unit can represent individual farms, farming districts, natural areas, or urban areas. The example model demonstrates the application of MF2005-FMP2 with delayed recharge through an unsaturated zone, rejected infiltration in a riparian area, changes in de

Techniques and Methods↗

USGS2009

GWM-2005 - A Groundwater-Management Process for MODFLOW-2005 with Local Grid Refinement (LGR) Capability

This report describes the Groundwater-Management (GWM) Process for MODFLOW-2005, the 2005 version of the U.S. Geological Survey modular three-dimensional groundwater model. GWM can solve a broad range of groundwater-management problems by combined use of simulation- and optimization-modeling techniques. These problems include limiting groundwater-level declines or streamflow depletions, managing groundwater withdrawals, and conjunctively using groundwater and surface-water resources. GWM was initially released for the 2000 version of MODFLOW. Several modifications and enhancements have been made to GWM since its initial release to increase the scope of the program's capabilities and to improve its operation and reporting of results. The new code, which is called GWM-2005, also was designed to support the local grid refinement capability of MODFLOW-2005. Local grid refinement allows for the simulation of one or more higher resolution local grids (referred to as child models) within a coarser grid parent model. Local grid refinement is often needed to improve simulation accuracy in regions where hydraulic gradients change substantially over short distances or in areas requiring detailed representation of aquifer heterogeneity. GWM-2005 can be used to formulate and solve groundwater-management problems that include components in both parent and child models. Although local grid refinement increases simulation accuracy, it can also substantially increase simulation run times.

Techniques and Methods↗

USGS2011

Documentation for the State Variables Package for the Groundwater-Management Process of MODFLOW-2005 (GWM-2005)

Many groundwater-management problems are concerned with the control of one or more variables that reflect the state of a groundwater-flow system or a coupled groundwater/surface-water system. These system state variables include the distribution of heads within an aquifer, streamflow rates within a hydraulically connected stream, and flow rates into or out of aquifer storage. This report documents the new State Variables Package for the Groundwater-Management Process of MODFLOW-2005 (GWM-2005). The new package provides a means to explicitly represent heads, streamflows, and changes in aquifer storage as state variables in a GWM-2005 simulation. The availability of these state variables makes it possible to include system state in the objective function and enhances existing capabilities for constructing constraint sets for a groundwater-management formulation. The new package can be used to address groundwater-management problems such as the determination of withdrawal strategies that meet water-supply demands while simultaneously maximizing heads or streamflows, or minimizing changes in aquifer storage. Four sample problems are provided to demonstrate use of the new package for typical groundwater-management applications.

Techniques and Methods↗

USGS2013

Advective transport observations with MODPATH-OBS--documentation of the MODPATH observation process

The MODPATH-OBS computer program described in this report is designed to calculate simulated equivalents for observations related to advective groundwater transport that can be represented in a quantitative way by using simulated particle-tracking data. The simulated equivalents supported by MODPATH-OBS are (1) distance from a source location at a defined time, or proximity to an observed location; (2) time of travel from an initial location to defined locations, areas, or volumes of the simulated system; (3) concentrations used to simulate groundwater age; and (4) percentages of water derived from contributing source areas. Although particle tracking only simulates the advective component of conservative transport, effects of non-conservative processes such as retardation can be approximated through manipulation of the effective-porosity value used to calculate velocity based on the properties of selected conservative tracers. This program can also account for simple decay or production, but it cannot account for diffusion. Dispersion can be represented through direct simulation of subsurface heterogeneity and the use of many particles. MODPATH-OBS acts as a postprocessor to MODPATH, so that the sequence of model runs generally required is MODFLOW, MODPATH, and MODPATH-OBS. The version of MODFLOW and MODPATH that support the version of MODPATH-OBS presented in this report are MODFLOW-2005 or MODFLOW-LGR, and MODPATH-LGR. MODFLOW-LGR is derived from MODFLOW-2005, MODPATH 5, and MODPATH 6 and supports local grid refinement. MODPATH-LGR is derived from MODPATH 5. It supports the forward and backward tracking of particles through locally refined grids and provides the output needed for MODPATH_OBS. For a single grid and no observations, MODPATH-LGR results are equivalent to MODPATH 5. MODPATH-LGR and MODPATH-OBS simulations can use nearly all of the capabilities of MODFLOW-2005 and MODFLOW-LGR; for example, simulations may be steady-state, transient, or a combination. Though the program name MODPATH-OBS specifically refers to observations, the program also can be used to calculate model prediction of observations. MODPATH-OBS is primarily intended for use with separate programs that conduct sensitivity analysis, data needs assessment, parameter estimation, and uncertainty analysis, such as UCODE_2005, and PEST. In many circumstances, refined grids in selected parts of a model are important to simulated hydraulics, detailed inflows and outflows, or other system characteristics. MODFLOW-LGR and MODPATH-LGR support accurate local grid refinement in which both mass (flows) and energy (head) are conserved across the local grid boundary. MODPATH-OBS is designed to take advantage of these capabilities. For example, particles tracked between a pumping well and a nearby stream, which are simulated poorly if a river and well are located in a single large grid cell, can be simulated with improved accuracy using a locally refined grid in MODFLOW-LGR, MODPATH-LGR, and MODPATH-OBS. The locally-refined-grid approach can provide more accurate simulated equivalents to observed transport between the well and the river. The documentation presented here includes a brief discussion of previous work, description of the methods, and detailed descriptions of the required input files and how the output files are typically used.

Techniques and Methods↗