USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Surface Techniques”

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

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

At least 865 records · Page 48Linked to original sources

Hyperspectral imaging of river bathymetry using an ensemble of regression trees

Remote sensing has emerged as an effective tool for characterizing river systems, and machine learning (ML) techniques could make this approach even more powerful. To explore this possibility, we developed an ML-based workflow for hyperspectral imaging of river bathymetry using an ensemble of regression trees (HIRBERT). This approach involves using paired observations of depth and reflectance to select wavelength bands as predictors and then train a depth retrieval model; applying the model to the image yields a spatially continuous bathymetric map. We used data from five rivers with diverse morphologies and optical characteristics to assess whether HIRBERT can (1) provide more accurate depth estimates than a band ratio-based algorithm and (2) extend the range of depths detectable via remote sensing. Relative to single band combinations identified via optimal band ratio analysis (OBRA), regression tree ensembles improved depth retrieval performance, with observed versus predicted (OP) regression R 2 values increasing for all five sites. Similarly, HIRBERT provided more reliable depth estimates than OBRA over the full range of depths present along each river. These results suggest that by incorporating additional spectral information from multiple wavelength bands, ML could enhance bathymetric mapping across a range of river environments. In addition, we show how graphical tools can facilitate interpretation of ML-based depth retrieval models and yield insight regarding relationships between depth and reflectance. The HIRBERT workflow is packaged in free, standalone software developed to support applications in river research and management. Although ML can enhance remote sensing of river bathymetry, the limitations of this approach must also be acknowledged: Field measurements of water depth are required to train a depth retrieval model and the resulting model should only be applied to the image from which the training data were derived. The inherently image-specific nature of this approach implies that developing generalized regression tree ensembles that could be applied at larger scales would require additional research.

California, Idaho, Nebraska, Oregon, Wyoming↗

Brine contamination to aquatic resources from oil and gas development in the Williston Basin, United States

The Williston Basin, which includes parts of Montana, North Dakota, and South Dakota in the United States and the provinces of Manitoba and Saskatchewan in Canada, has been a leading domestic oil and gas producing region for more than one-half a century. Currently, there are renewed efforts to develop oil and gas resources from deep geologic formations, spurred by advances in recovery technologies and economic incentives associated with the price of oil. Domestic oil and gas production has many economic benefits and provides a means for the United States to fulfill a part of domestic energy demands; however, environmental hazards can be associated with this type of energy production in the Williston Basin, particularly to aquatic resources (surface water and shallow groundwater) by extremely saline water, or brine, which is produced with oil and gas. The primary source of concern is the migration of brine from buried reserve pits that were used to store produced water during recovery operations; however, there also are considerable risks of brine release from pipeline failures, poor infrastructure construction, and flow-back water from hydraulic fracturing associated with modern oilfield operations. During 2008, a multidisciplinary (biology, geology, water) team of U.S. Geological Survey researchers was assembled to investigate potential energy production effects in the Williston Basin. Researchers from the U.S. Geological Survey participated in field tours and met with representatives from county, State, tribal, and Federal agencies to identify information needs and focus research objectives. Common questions from agency personnel, especially those from the U.S. Fish and Wildlife Service, were “are the brine plumes (plumes of brine-contaminated groundwater) from abandoned oil wells affecting wetlands on Waterfowl Production Areas and National Wildlife Refuges?” and “are newer wells related to Bakken and Three Forks development different than the older, abandoned wells (in terms of potential for affecting aquatic resources)?” Of special concern were the wetland habitats of the ecologically important Prairie Pothole Region, which overlays a part of the Williston Basin and is recognized for the production of a majority of North America’s migratory waterfowl. On the basis of the concerns raised by on-the-ground land managers, as well as findings from previous research, a comprehensive study was developed with the following goals: summarize existing information pertaining to oil and gas production and aquatic resources in the Williston Basin; assess brine plume migration from new and previously studied sites in the Prairie Pothole Region; perform a regional, spatial evaluation of oil and gas production activities and aquatic resources; assess the potential for brine contamination to wetlands and streams; and hold a decision analysis workshop with key stakeholders to discuss issues pertaining to oil and gas production and environmental effects and to identify information gaps and research needs. This report represents an initial, multidisciplinary evaluation of measured and potential environmental effects associated with oil and gas production in the Williston Basin and Prairie Pothole Region. Throughout this report there are reviews of current knowledge, and discussions relating to data gaps and research needs. On the basis of the information presented, future research needs include: regional geophysical and water-quality assessments to establish baselines for current conditions and estimate the extent of previous brine contamination, investigations into the direct effects of brine to biotic communities, and evaluations to identify the most effective techniques to mitigate brine contamination.

Montana, North Dakota↗

Simulation of streamflow and sediment transport in two surface-coal-mined basins in Fayette County, Pennsylvania

The Hydrological Simulation Program - Fortran (HSPF) was used to simulate streamflow and sediment transport in two surface-mined basins of Fayette County, Pa. Hydrologic data from the Stony Fork Basin (0.93 square miles) was used to calibrate HSPF parameters. The calibrated parameters were applied to an HSPF model of the Poplar Run Basin (8.83 square miles) to evaluate the transfer value of model parameters. The results of this investigation provide information to the Pennsylvania Department of Environmental Resources, Bureau of Mining and Reclamation, regarding the value of the simulated hydrologic data for use in cumulative hydrologic-impact assessments of surface-mined basins. The calibration period was October 1, 1985, through September 30, 1988 (water years 1986-88). The simulated data were representative of the observed data from the Stony Fork Basin. Mean simulated streamflow was 1.64 cubic feet per second compared to measured streamflow of 1.58 cubic feet per second for the 3-year period. The difference between the observed and simulated peak stormflow ranged from 4.0 to 59.7 percent for 12 storms. The simulated sediment load for the 1987 water year was 127.14 tons (0.21 ton per acre), which compares to a measured sediment load of 147.09 tons (0.25 ton per acre). The total simulated suspended-sediment load for the 3-year period was 538.2 tons (0.30 ton per acre per year), which compares to a measured sediment load of 467.61 tons (0.26 ton per acre per year). The model was verified by comparing observed and simulated data from October 1, 1988, through September 30, 1989. The results obtained were comparable to those from the calibration period. The simulated mean daily discharge was representative of the range of data observed from the basin and of the frequency with which specific discharges were equalled or exceeded. The calibrated and verified parameters from the Stony Fork model were applied to an HSPF model of the Poplar Run Basin. The two basins are in a similar physical setting. Data from October 1, 1987, through September 30, 1989, were used to evaluate the Poplar Run model. In general, the results from the Poplar Run model were comparable to those obtained from the Stony Fork model. The difference between observed and simulated total streamflow was 1.1 percent for the 2-year period. The mean annual streamflow simulated by the Poplar Run model was 18.3 cubic feet per second. This compares to an observed streamflow of 18.15 cubic feet per second. For the 2-year period, the simulated sediment load was 2,754 tons (0.24 ton per acre per year), which compares to a measured sediment load of 3,051.2 tons (0.27 ton per acre per year) for the Poplar Run Basin. Cumulative frequency-distribution curves of the observed and simulated streamflow compared well. The comparison between observed and simulated data improved as the time span increased. Simulated annual means and totals were more representative of the observed data than hourly data used in comparing storm events. The structure and organization of the HSPF model facilitated the simulation of a wide range of hydrologic processes. The simulation results from this investigation indicate that model parameters may be transferred to ungaged basins to generate representative hydrologic data through modeling techniques.

Pennsylvania↗

Simulation of groundwater flow and streamflow depletion in the Branch Brook, Merriland River, and parts of the Mousam River watersheds in southern Maine

Watersheds of three streams, the Mousam River, Branch Brook, and Merriland River in southeastern Maine were investigated from 2010 through 2013 under a cooperative project between the U.S. Geological Survey and the Maine Geological Survey. The Branch Brook watershed previously had been deemed “at risk” by the Maine Geological Survey because of the proportionally large water withdrawals compared to estimates of the in-stream flow requirements for habitat protection. The primary groundwater withdrawals in the study area include a water-supply well in the headwaters of the system and three water-supply wells in the coastal plain near the downstream end of the system. A steady-state groundwater flow model was used to understand the movement of water within the system, to evaluate the water budget and the effect of groundwater withdrawals on streamflows, and to understand streamflow depletion in relation to the State of Maine’s requirements to maintain in-stream flows for habitat protection. Delineation of the simulated groundwater divides compared to the surface-water divides suggests that the groundwater divides in the headwater areas do not exactly correspond to the surface-water divides. Under both pumping and non-pumping conditions, groundwater flows from the headwaters of the Branch Brook watershed into the Mousam River watershed. Pumping in the Mousam River watershed captures a small amount of groundwater from the Branch Brook basin. The cumulative effect of groundwater withdrawals on base flows in two rivers in the study area (Branch Brook and the Merriland River) was evaluated using the groundwater flow model. Streamflow depletion in the headwaters of Branch Brook was 0.12 cubic feet per second (ft 3 /s) for the steady-state simulation, or about 10 percent of the average base flow at that location. Downstream on Branch Brook, the total streamflow depletion from all the wells was 0.59 ft 3 /s, or 3 percent of the average base flow at that location. In the Merriland River downstream from the Merriland River well, the total amount of streamflow depletion was 0.6 ft 3 /s, or about 7 percent of the average base flow. The groundwater model was used to evaluate several different scenarios that could affect streamflow and groundwater discharging to the rivers and streams in the study area. The scenarios were (1) no pumping from the water-supply wells; (2) current pumping from the water-supply wells, but simulated drought conditions (25 percent reduction in recharge); (3) current recharge, but with increased pumping from the large water-supply wells; and (4) drought conditions and increased pumping combined. Simulations of increased pumping in the water-supply wells resulted in streamflow depletion in the headwaters of Branch Brook increasing to 16 percent of the headwater base flow. Simulated increases in the pumping in the coastal plain wells increased the amount of streamflow depletion to 6 percent of the flow in Branch Brook and to 8 percent of the flow in the Merriland River. The additional stress of a drought imposed on the model (25 percent less recharge) had a substantial impact on streamflows, as expected. If the simulated drought occurred simultaneously with an increase in pumping, the base flows would be reduced 48 percent in the headwaters of Branch Brook, compared to the no-pumping scenario. Downstream in Branch Brook, the total reduction in flow would be 29 percent of the simulated base flows in the no-pumping scenario, and in the Merriland River, the reduction would be 33 percent of the base flows in the no-pumping scenario. The study evaluated two different methods of calculating in-stream flow requirements for Branch Brook and the Merriland River—a set of statewide equations used to calculate monthly median flows and the MOVE.1 record-extension technique used on site-specific streamflow measurements. The August median in-stream flow requirement in the Merriland River was calculated as 7.18 ft 3 /s using the statewide equations but was 3.07 ft 3 /s using the MOVE.1 analysis. In Branch Brook, the August median in-stream flow requirements were calculated as 20.3 ft 3 /s using the statewide equations and 11.8 ft 3 /s using the MOVE.1 analysis. In each case, using site-specific data yields an estimate of in-stream flow that is much lower than an estimate the statewide equations provide.

Maine↗

Study and interpretation of the chemical characteristics of natural water

The chemical composition of natural water is derived from many different sources of solutes, including gases and aerosols from the atmosphere, weathering and erosion of rocks and soil, solution or precipitation reactions occurring below the land surface, and cultural effects resulting from activities of man. Some of the processes of solution or precipitation of minerals can be closely evaluated by means of principles of chemical equilibrium including the law of mass action and the Nernst equation. Other processes are irreversible and require consideration of reaction mechanisms and rates. The chemical composition of the crustal rocks of the earth and the composition of the ocean and the atmosphere are significant in evaluating sources of solutes in natural fresh water. The ways in which solutes are taken up or precipitated and the amounts present in solution are influenced by many environmental factors, especially climate, structure and position: of rock strata, and biochemical effects associated with life cycles of plants and animals, both microscopic and macroscopic. Taken all together and in application with the further influence of the general circulation of all water in the hydrologic cycle, the chemical principles and environmental factors form a basis for the developing science of natural-water chemistry. Fundamental data used in the determination of water quality are obtained by the chemical analysis of water samples in the laboratory or onsite sensing of chemical properties in the field. Sampling is complicated by changes in composition of moving water and the effects of particulate suspended material. Most of the constituents determined are reported in gravimetric units, usually milligrams per liter or milliequivalents per liter. More than 60 constituents and properties are included in water analyses frequently enough to provide a basis for consideration of the sources from which each is generally derived, most probable forms of elements and ions in solution, solubility controls, expected concentration ranges and other chemical factors. Concentrations of elements that are commonly present in amounts less than a few tens of micrograms per liter cannot always be easily explained, but present information suggests many are controlled by solubility of hydroxide or carbonate or by sorption on solid particles. Chemical analyses may be grouped and statistically evaluated by averages, frequency distributions, or ion correlations to summarize large volumes of data. Graphing of analyses or of groups of analyses aids in showing chemical relationships among waters, probable sources of solutes, areal water-quality regimen, and water-resources evaluation. Graphs may show water type based on chemical composition, relationships among ions, or groups of ions in individual waters or many waters considered simultaneously. The relationships of water quality to hydrologic parameters, such as stream discharge rate or ground-water flow patterns, can be shown by mathematical equations, graphs, and maps. About 75 water analyses selected from the literature are tabulated to illustrate the relationships described, and some of these, along with many others that are not tabulated, are also utilized in demonstrating graphing and mapping techniques. Relationships of water composition to source rock type are illustrated by graphs of some of the tabulated analyses. Activities of man may modify water composition extensively through direct effects of pollution and indirect results of water development, such as intrusion of sea water in ground-water aquifiers. Water-quality standards for domestic, agricultural, and industrial use have been published by various agencies. Irrigation project requirements for water quality are particularly intricate. Fundamental knowledge of processes that control natural water composition is required for rational management of water quality.

Water Supply Paper↗

Study and interpretation of the chemical characteristics of natural water

The chemical composition of natural water is derived from many different sources of solutes, including gases and aerosols from the atmosphere, weathering and erosion of rocks and soil, solution or precipitation reactions occurring below the land surface, and cultural effects resulting from activities of man. Some of the processes of solution or precipitation of minerals can be closely evaluated by means of principles of chemical equilibrium including the law of mass action and the Nernst equation. Other processes are irreversible and require consideration of reaction mechanisms and rates. The chemical composition of the crustal rocks of the earth and the composition of the ocean and the atmosphere are significant in evaluating sources of solutes in natural fresh water. The ways in which solutes are taken up or precipitated and the amounts present in solution are influenced by many environmental factors, especially climate, structure and position: of rock strata, and biochemical effects associated with life cycles of plants and animals, both microscopic and macroscopic. Taken all together and in application with the further influence of the general circulation of all water in the hydrologic cycle, the chemical principles and environmental factors form a basis for the developing science of natural-water chemistry. Fundamental data used in the determination of water quality are obtained by the chemical analysis of water samples in the laboratory or onsite sensing of chemical properties in the field. Sampling is complicated by changes in composition of moving water and the effects of particulate suspended material. Most of the constituents determined are reported in gravimetric units, usually milligrams per liter or milliequivalents per liter. More than 60 constituents and properties are included in water analyses frequently enough to provide a basis for consideration of the sources from which each is generally derived, most probable forms of elements and ions in solution, solubility controls, expected concentration ranges and other chemical factors. Concentrations of elements that are commonly present in amounts less than a few tens of micrograms per liter cannot always be easily explained, but present information suggests many are controlled by solubility of hydroxide or carbonate or by sorption on solid particles. Chemical analyses may be grouped and statistically evaluated by averages, frequency distributions, or ion correlations to summarize large volumes of data. Graphing of analyses or of groups of analyses aids in showing chemical relationships among waters, probable sources of solutes, areal water-quality regimen, and water-resources evaluation. Graphs may show water type based on chemical composition, relationships among ions, or groups of ions in individual waters or many waters considered simultaneously. The relationships of water quality to hydrologic parameters, such as stream discharge rate or ground-water flow patterns, can be shown by mathematical equations, graphs, and maps. About 75 water analyses selected from the literature are tabulated to illustrate the relationships described, and some of these, along with many others that are not tabulated, are also utilized in demonstrating graphing and mapping techniques. Relationships of water composition to source rock type are illustrated by graphs of some of the tabulated analyses. Activities of man may modify water composition extensively through direct effects of pollution and indirect results of water development, such as intrusion of sea water in ground-water aquifiers. Water-quality standards for domestic, agricultural, and industrial use have been published by various agencies. Irrigation project requirements for water quality are particularly intricate. Fundamental knowledge of processes that control natural water composition is required for rational management of water quality.

Water Supply Paper↗

Survival, transport, and sources of fecal bacteria in streams and survival in land-applied poultry litter in the upper Shoal Creek basin, southwestern Missouri, 2001–2002

Densities of fecal coliform bacteria along a 5.7-mi (mile) reach of Shoal Creek extending upstream from State Highway 97 (site 3) to State Highway W (site 2) and in two tributaries along this reach exceeded the Missouri Department of Natural Resources (MDNR) standard of 200 col/100 mL (colonies per 100 milliliters) for whole-body contact recreation. A combination of techniques was used in this report to provide information on the source, transport, and survival of fecal bacteria along this reach of Shoal Creek. Results of water-quality samples collected during dye-trace and seepage studies indicated that at summer low base-flow conditions, pastured cattle likely were a substantial source of fecal bacteria in Shoal Creek at the MDNR monitoring site (site 3) at State Highway 97. Using repeat element Polymerase Chain Reaction (rep-PCR), cattle were the presumptive source of about 50 percent of the Escherichia coli ( E. coli ) isolates in water samples from site 3. Cattle, horses, and humans were the most common presumptive source of E. coli isolates at sites further upstream. Poultry was identified by rep-PCR as a major source of E. coli in Pogue Creek, a tributary in the upper part of the study area. Results of the rep-PCR were in general agreement with the detection and distribution of trace concentrations of organic compounds commonly associated with human wastewater, such as caffeine, the antimicrobial agent triclosan, and the pharmaceutical compounds acetaminophen and thiabendazole (a common cattle anthelmintic). Significant inputs of fecal bacteria to Shoal Creek occurred along a 1.6-mi reach of Shoal Creek immediately upstream from site 3. During a 36-hour period in July 2001, average densities of fecal coliform and E. coli bacteria increased from less than or equal to 500 col/100 mL upstream from this stream reach (sample site 2c) to 2,100 and 1,400 col/100 mL, respectively, at the MDNR sampling site. Fecal bacteria densities exhibited diurnal variability at all five sampling sites along the 5.7-mi study reach of Shoal Creek, but the trends at successive downstream sites were out of phase and could not be explained by simple advection and dispersion. At base-flow conditions, the travel time of bacteria in Shoal Creek along the 5.7-mi reach between State Highway W (site 2) and the MDNR sampling site (site 3) was about 26 hours. Substantial dispersion and dilution occurs along the upper 4.1 mi of this reach because of inflows from a number of springs and tributaries and the presence of several long pools and channel meanders. Minimal dispersion and dilution occurs along the 1.6-mi reach immediately upstream from the MDNR sampling site. Measurements of fecal bacteria decay in Shoal Creek during July 2001 indicated that about 8 percent of fecal coliform and E. coli bacteria decay each hour with an average first-order decay constant of 0.084 h -1 (per hour). Results of field test plots indicated that substantial numbers of fecal bacteria present in poul try litter can survive in fields for as much as 8 weeks after the application of the litter to the land surface. Median densities of fecal coliform and E. coli in slurry-water samples collected from fields increased from less than 60 col/100 mL before the application of turkey and broiler litter, to as large as 420,000 and 290,000 col/100 mL after the application of litter. Bacteria densities in the test plots generally decreased in a exponential manner over time with decay rates ranging from 0.085 to 0.185 d -1 (per day) for fecal coliform to between 0.100 and 0.250 d -1 for E. coli . The apparent survival of significant numbers of fecal bacteria on fields where poultry litter has been applied indicates that runoff from these fields is a potential source of fecal bacteria to vicinity streams for many weeks following litter application.

Water-Resources Investigations Report↗

Study and interpretation of the chemical characteristics of natural water

The chemical composition of natural water is derived from many different sources of solutes, including gases and aerosols from the atmosphere, weathering and erosion of rocks and soil, solution or precipitation reactions occurring below the land surface, and cultural effects resulting from human activities. Broad interrelationships among these processes and their effects can be discerned by application of principles of chemical thermodynamics. Some of the processes of solution or precipitation of minerals can be closely evaluated by means of principles of chemical equilibrium, including the law of mass action and the Nernst equation. Other processes are irreversible and require consideration of reaction mechanisms and rates. The chemical composition of the crustal rocks of the Earth and the composition of the ocean and the atmosphere are significant in evaluating sources of solutes in natural freshwater. The ways in which solutes are taken up or precipitated and the amounts present in solution are influenced by many environmental factors, especially climate, structure and position of rock strata, and biochemical effects associated with life cycles of plants and animals, both microscopic and macroscopic. Taken together and in application with the further influence of the general circulation of all water in the hydrologic cycle, the chemical principles and environmental factors form a basis for the developing science of natural-water chemistry. Fundamental data used in the determination of water quality are obtained by the chemical analysis of water samples in the laboratory or onsite sensing of chemical properties in the field. Sampling is complicated by changes in the composition of moving water and by the effects of particulate suspended material. Some constituents are unstable and require onsite determination or sample preservation. Most of the constituents determined are reported in gravimetric units, usually milligrams per liter or milliequivalents per liter. More than 60 constituents and properties are included in water analyses frequently enough to provide a basis for consideration of the sources from which each is generally derived, the most probable forms of elements and ions in solution, solubilitycontrols, expected concentration ranges, and other chemical factors. Mechanisms that control concentrations of elements commonly present in amounts less than a few tens of micrograms per liter cannot always be easily ascertained, but present information suggests that many are controlled by solubility of their hydroxides or carbonates or by sorption on solid particles. Many dissolved organic compounds can now be specifically determined. Chemical analyses may be grouped and statistically evaluated by means, medians, frequency distributions, or ion correlations to summarize large volumes of data. Graphing of analyses or of groups of analyses aids in showing chemical relationships among water, probable sources of solutes, areal water-quality regimen, temporal and spatial variation, and water-resources evaluation. Graphs may show water type based on chemical composition, relationships among ions, or groups of ions in individual waters or many waters considered simultaneously. The relationships of water quality to hydrogeologic characteristics, such as stream discharge rate or ground-water flow patterns, can be shown by mathematical equations, graphs, and maps. About 80 water analyses selected from the literature are tabulated to illustrate the relationships described, and some of these, along with many others that are not tabulated, are also used in demonstrating graphing and mapping techniques. Relationships of water composition to source rock type are illustrated by graphs of some of the tabulated analyses. Human activities may modify water composition extensively through direct effects of pollution and indirect results of water development, such as intrusion of seawater in groundwater aquifers. Water-quality standards for domestic, agricultural, and industrial use have been published by various agencies. Irrigation project requirements for water quality are particularly intricate. Fundamental knowledge of processes that control natural-water composition is required for rational management of water quality.

Water Supply Paper↗

Santa Barbara and Foothill groundwater basins Geohydrology and optimal water resources management—Developed using density dependent solute transport and optimization models

Groundwater has been a part of the city of Santa Barbara’s water-supply portfolio since the 1800s; however, since the 1960s, the majority of the city’s water has come from local surface water, and the remainder has come from groundwater, State Water Project, recycled water, increased water conservation, and as needed, seawater desalination. Although groundwater from the Santa Barbara and Foothill groundwater basins only accounts for a small percentage of the long-term supply, it is an important source of supplemental water during times of surface-water shortages. During the late 1980s and early 1990s, production wells extracted additional groundwater to compensate for drought related water-delivery shortfalls from other sources; in response, water levels declined substantially in the Santa Barbara and Foothill groundwater basins (below sea level in the Santa Barbara groundwater basin). In coastal basins that have groundwater extraction near shore, seawater intrusion is often a problem. Seawater intrusion in the Santa Barbara groundwater basin is thought to be more limited than in other coastal basins because of an offshore fault that acts as a partial barrier to groundwater flow. During the late 1980s and early 1990s, seawater intrusion was observed in the Santa Barbara groundwater basin, as indicated by increased chloride concentrations at several monitoring wells that ranged from 200 ft to 1,300 ft from the ocean and as close as 2,900 ft to the nearest pumping well. This demonstrated that seawater can intrude into the Santa Barbara groundwater basin when groundwater levels fall below sea level near the coast. The city of Santa Barbara is interested in developing a better understanding of the sustainability of its groundwater supplies. In 2014, California adopted historic legislation to manage its groundwater: the Sustainable Groundwater Management Act (SGMA). The SGMA requires the development and implementation of “Groundwater Sustainability Plans” in 127 priority groundwater basins; although Santa Barbara was not a designated priority basin, the city is taking steps to achieve sustainability. Sustainability was defined in the SGMA in terms of avoiding undesirable results: significant and unreasonable groundwater-level declines, reduction in groundwater storage, seawater intrusion, water-quality degradation, land subsidence, and surface-water depletion. In this project, a cooperative study between the U.S. Geological Survey (USGS) and the city of Santa Barbara, sustainable yield is defined as the volume of groundwater that can be pumped from storage without causing water-level drawdowns and the associated increases in seawater intrusion (as indicated by increases in measured chloride concentrations) at selected wells. In order to estimate the sustainability of Santa Barbara’s groundwater basins, a three-dimensional density-dependent groundwater-flow and solute-transport model (the Santa Barbara Flow and Transport Model, or SBFTM) was developed on the basis of an existing groundwater-flow model. To simulate seawater intrusion to the Santa Barbara Basin under various management strategies, the SBFTM uses the USGS code SEAWAT to simulate salinity transport and variable-density flow. The completed SBFTM was coupled with a management optimization tool, in this case a multi-objective evolutionary algorithm, to determine optimal pumping strategies that maximize the sustainable yield and at the same time satisfy user-defined drawdown and chloride-concentration constraints. As part of this study, a three-dimensional hydrogeologic framework model was developed to quantify the extent and hydrogeologic characteristics of the Santa Barbara and Foothill groundwater basins and to help define the discretization and hydraulic properties used in the SBFTM. The development of the hydrogeologic framework model required the collection and reconciliation of geologic and geophysical data from existing maps, reports, and databases, along with geologic and hydrologic data from recently drilled wells. These data were integrated into a three-dimensional hydrogeologic framework model that defines the stratigraphy and geometry of the aquifer zones and the major geologic structures in the basin. The hydrogeologic framework model also quantifies the variation in sediment grain size within each aquifer zone as the percentage of coarse-grained sediment. Previous studies indicated that there are two principal water-producing zones in the Santa Barbara groundwater basin, the upper and lower producing zones; an additional thin, productive zone was identified as part of this study. This “middle producing zone” is not as areally extensive as the upper and lower producing zones and only exists in the coastal part of Storage Unit I. These producing zones are bounded at depth by less productive shallow, middle, and deep zones. Two versions of the SBFTM were constructed: an initial-condition model and a modern transient model. The initial-condition model is a long-term transient model that simulates flow and solute-transport conditions during a period with limited anthropogenic influences preceeding the modern transient model. The simulation-transient model simulates flow and transport conditions from 1929 through 2013; however, because of data availability, the focus of the model calibration was 1972–2013. The SBFTM was calibrated to measured groundwater levels and drawdown, as well as measured chloride concentrations and change in concentrations, using a combination of automated and trial-and-error parameter-estimation techniques. A sensitivity analysis indicated that, in general, the SBFTM was most sensitive to recharge- and pumping-distribution parameters, specifically those controlling the amount of small-catchment recharge and the distribution of water extraction by hydrogeologic layer for production wells. The model was also sensitive to parameters controlling stream-recharge rates, horizontal and vertical hydraulic conductivity, and porosity. From 1929 to 1971, most of the water entering the area represented by the SBFTM was from creek and small-catchment recharge, and the majority of water leaving the SBFTM area was from pumping, discharge to creeks, and drains. In addition, about 37 percent of the total pumpage came from a net reduction in groundwater storage. From 1972 to 2013, the amount of water entering and leaving the SBFTM was fairly similar as that from 1929 to 1971, except the reduction in pumpage added about 17,000 acre-ft of water to storage. During this later period, there were also times of storage loss. For example, during July 1990, a month when approximately 705 acre-ft of groundwater was pumped in the study area, the pumpage was much greater than all sources of recharge combined, and about 382 acre-ft of water was removed from groundwater storage. Simulated hydraulic heads replicated the observed data to an acceptable matching of the measured water-level, flow direction, and vertical gradients. Simulated hydrographs for selected wells were in good agreement with the measured data, with an average residual of -2.7 ft and a standard deviation of 14.5 ft, indicating that the simulated heads, on average, underestimated the observed water levels. An examination of the model fit indicated that most of the discrepancies were lower simulated heads at wells proximal to production well sites. The simulated chloride concentrations reasonably matched the rising limbs of the measured breakthrough curves in terms of timing and magnitude; however, the simulation overestimated the chloride concentrations on the falling limbs. The overestimation of low chloride concentrations was attributed to the model overestimating the advance of the chloride front during periods of heavy pumping and underestimating the retreat of the chloride front during periods of low pumping. These simulation errors would result in a conservative response by local water managers to seawater intrusion. The SBFTM was used to develop a collection of predictive simulations optimized to produce pumping schedules that maximize yield, subject to a set of constraints and competing objectives. The simulations were grouped as scenarios that differed in their time horizon, initial conditions for groundwater levels and chloride concentrations, as well as precipitation, which was incorporated into the model through simulated recharge. Overall, five scenarios were developed in a multi-objective framework to obtain optimal pumping rates for all of the wells managed by the city, while minimizing excessive drawdown and seawater intrusion. For the current study, complexities in the simulation model and the optimization formulation required additional considerations. Incorporating the solute-transport equations to simulate chloride transport added a highly nonlinear process that is solved iteratively in each time step of the groundwater-flow model. These nonlinearities, coupled with the highly refined grid in the current model, creates challenges for many traditional optimization methods. Therefore, an optimization method was needed that could address nonlinear relationships as well as a very large problem size. Lastly, the optimization problem was reformulated to include multiple objectives without requiring convergence to a single solution. This approach, guided by the city’s objectives, allowed the maximum extraction of information from the complex simulation. Borg, a multi-objective evolutionary algorithm, was chosen as the optimization algorithm for this study for several reasons: (1) it is very computationally efficient; (2) it can run in parallel; (3) it requires little user input; and (4) it can solve for multiple competing objectives. The first three points allow the algorithm to proceed toward the optimal solutions at the fastest possible rate. The fourth point is advantageous for large, complex optimization problems because it is difficult to formulate the optimization problem in a way that produces only one optimal solution. The problem formulation consisted of four competing objectives and a constraint set in accordance with the main concerns of the city. The objectives were maximizing total pumpage, minimizing seawater intrusion, minimizing total drawdown in production wells, and minimizing the maximum drawdown. The constraints were pump capacity, meeting drinking-water standards for chloride, maintaining a specified minimum flowrate to a groundwater treatment plant, and maintaining minimum water levels in pumping wells. The decision variables either were quarterly pumpage by well or total pumpage by basin. Five optimization scenarios were developed that allow the decision makers to evaluate a range of optimal solutions for a variety of water levels and chloride concentrations as well as potential future climatic conditions. Three scenarios (1, 2, and 5) were multi-objective optimization formulations that allowed for variations in management preferences and climatic conditions. The other two scenarios (3 and 4) were designed to examine the optimization results to answer specific questions. Scenario 1 described the best-case sustainable yield assuming a “full” basin (that is, high initial water levels) and typical climate conditions for 10 years. Scenario 2 also started with a “full” basin; however, this was followed by a 10-year drought. Scenario 3 determined if an “empty” basin (that is, low initial water levels) would recover to full conditions (1998 conditions) given climate assumptions and optimal pumping schedules from scenarios 1 and 2. Scenario 4 was designed to produce decision rules that can be used by water managers to help choose an optimal pumping schedule based on measured water-level or chloride data. Scenario 5 identified future pumping schedules based on short-term climate variations during a 2-year management horizon. The results from scenarios 1 and 2 described the differences in maximum pumpage in the basin under typical and dry long-term climate projections, respectively. The scenario 1 results indicated the maximum 10-year pumpage of the basin was about 31,300 acre-ft under typical conditions and controlling simulated seawater intrusion and drawdowns. For scenario 2, less recharge over the 10-year dry climate produced a maximum pumpage estimate of 30,000 acre-ft to control seawater intrusion and drawdowns. The larger pumpage for scenario 1 resulted in more seawater intrusion, but less total drawdown, compared to that of scenario 2. Results for scenarios 3 and 4 showed the basin’s response to management actions combined with climate projections. Both scenarios used the optimal pumping schedules and the 10-year climates from scenarios 1 and 2. The scenario 3 results showed that under minimal pumping, the basin did not fully recover to 1998 water levels within 10 years under either climate scenario. The relatively larger recharge from the typical climate resulted in less drawdown at coastal monitoring wells after the 10-year recovery period than that from the dry climate. The location of the seawater intrusion front was not appreciably different between the scenarios, however. Scenario 4 used the optimal results from scenarios 1 and 2 to produce decision-rule curves that illustrated the pumpage for each basin, given measured levels of chloride concentration or drawdown. This allowed the use of additional measurements at monitoring wells to assess future management decisions on the basis of the sensitivity of observations of drawdown and seawater intrusion to various pumping rates. Scenario 5 allowed managers to investigate the effects of short-term climate variations on optimal pumping schedules. Three specific 2-year simulations were optimized: typical-to-dry (scenario 5A), dry-to-typical (scenario 5B), and dry-to-dry (scenario 5C). The most noteable result from scenario 5 was the overall reduction in optimal pumpage for most schedules in scenario 5C, when the climate is simulated as dry-to-dry. There are also many optimal pumping schedules that produced an overall increase in waterlevels over the two-year simulation period, regardless of climatic condition. Similar to scenario 2, the scenario 5C results represents conservative yield estimates under a minimal-precipitation climatic condition.

California↗

Apparatus and techniques for measuring bedload

The need for accurate determinations of the total sediment discharge of particles of bedload size has prompted this investigation of available and possible measuring apparatus and procedures. The accuracy of measurements of sediment discharge made with trap-type samplers is affected by the variability of sampler efficiency, by the oscillatory variation of bedload discharge, and by sampler placement. Equations that were developed for determining total discharge from measured bedioad discharge and measured suspended-sediment discharge are simplest if the bedload apparatus measures only the true bedload. Early bedload samplers are generally unsatisfactory. Recently developed or suggested apparatus include various improved samplers of the pressure-difference type, a pumping sampler, a magnetic sampler, acoustical instruments that measure the magnitude of the sound of particle collisions, an ultrasonic bedload sampler designed to measure and integrate electronically the concentration and velocity, and a tiltmeter designed to measure the total sediment discharge from the ground tilt that results from the passage of flow. All the pressure-difference samplers are improvements over early samplers, but none are void of the inherent shortcomings of trap-type apparatus; probably the Sphinx (Dutch) and VUV (Hungarian) samplers are the most satisfactory. The acoustical instruments are capable of measuring only the relative discharge. The ultrasonic sampler and the tiltmeter are not adequate without further development. Some new possible apparatus and means for measuring or aiding in measuring bedload discharge are small pit samplers, ultrasonic sounders, pressure transducers, and photography. A small pit sampler for measuring bedload discharge was designed to provide self-placement and portability ; however, its practicability and efficiency are undetermined. Exploratory films show that by using slowmotion photography the discharge of particles larger than about pea size can be determined provided the flow is clear; however, photography generally is not practical. Ultrasonic sounders provide continuous and accurate data on bed configuration and dune movement for use in equations that were developed for computing the bedload discharge. Computations with the equations indicate that the interpretation of the sounding data needs further study. Pressure transducers placed beneath the bed surface possibly can be used to provide information on dune movement; however, their installation would be difficult. The time required for collecting data on bed configuration and dune movement throughout a cross section could be substantially reduced by using several transducers simultaneously in conjunction with an ultrasonic sounder. A modified ultrasonic sounder that provides information on the shape and velocity of large particles and a method for determining the discharge of such particles were proposed; the method seems most feasible for particles of high sphericity.

Water Supply Paper↗

Numerical simulation of groundwater flow, resource optimization, and potential effects of prolonged drought for the Citizen Potawatomi Nation Tribal Jurisdictional Area, central Oklahoma

A hydrogeological study including two numerical groundwater-flow models was completed for the Citizen Potawatomi Nation Tribal Jurisdictional Area of central Oklahoma. One numerical groundwater-flow model, the Citizen Potawatomi Nation model, encompassed the jurisdictional area and was based on the results of a regional-scale hydrogeological study and numerical groundwater flow model of the Central Oklahoma aquifer, which had a geographic extent that included the Citizen Potawatomi Nation Tribal Jurisdictional Area. The Citizen Potawatomi Nation numerical groundwater-flow model included alluvial aquifers not in the original model and improved calibration using automated parameter-estimation techniques. The Citizen Potawatomi Nation numerical groundwater-flow model was used to analyze the groundwater-flow system and the effects of drought on the volume of groundwater in storage and streamflow in the North Canadian River. A more detailed, local-scale inset model was constructed from the Citizen Potawatomi Nation model to estimate available groundwater resources for two Citizen Potawatomi Nation economic development zones near the North Canadian River, the geothermal supply area and the Iron Horse Industrial Park. Groundwater pumping rates at potential well locations were optimized using the most recent version of the U.S. Geological Survey Groundwater-Management Process for MODFLOW. The objectives of optimization were to determine if a total pumping rate of 500 gallons per minute could be pumped from 5 wells at the geothermal supply area and to maximize discharge from 16 wells at the Iron Horse Industrial Park without exceeding specified head drawdown constraints at the pumping wells and thus prevent groundwater depletion. The inset model was used to estimate North Canadian River streamflow depletion caused by optimized pumping at the Iron Horse Industrial Park because water quality was a concern, and the river may have degraded water quality compared to water in other parts of the alluvial aquifer. The fate of streamflow that infiltrates into groundwater because of pumping was not directly determined, but it was assumed that this water could end up in the well discharge, and was considered to be a maximum proportion of well discharge derived from the North Canadian River. The total optimized continuous pumping rate from five managed wells at the geothermal supply area was 638 gallons per minute, which exceeded the target pumping rate of 500 gallons per minute. The total continuous pumping rate from 16 wells at the Iron Horse Industrial Park was 1,472 gallons per minute, which induced stream infiltration of approximately 4.1 gallons per minute (approximately 0.3 percent of the total well discharge) from the North Canadian River. To estimate the effects of drought on water resources in the Citizen Potawatomi Nation Tribal Jurisdictional Area, a hypothetical 10-year drought during which precipitation would decrease by 50 percent was simulated by decreasing model groundwater recharge by the same proportion for the period 1990–2000 of the transient model. The effects of the drought were estimated by calculating the change in the volume of groundwater storage and groundwater flow to streams at the end of the drought period, and the change in simulated streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., during and after the drought. The hypothetical decrease in recharge during the simulated drought caused groundwater in storage over the entire model in the study area to decrease by 361,500 acre-feet (14,100 acre-feet in the North Canadian River alluvial aquifer and 347,400 acre-feet in the Central Oklahoma aquifer), or approximately 0.2 percent of the total groundwater in storage over the drought period. This small percentage of groundwater loss showed that the Central Oklahoma aquifer as a bedrock aquifer has relatively low rates of recharge from the surface relative to the approximate storage. The budget for base flow to the North Canadian River indicated that the change in groundwater flow to the North Canadian River decreased during the 10-year drought by 386,500 acre-feet, or 37 percent. In all other parts of the Citizen Potawatomi Nation Tribal Jurisdictional Area, base flow decreased by 292,000 acre-feet, or 28 percent. Streamflow in the North Canadian River at the streamflow-gaging station at Shawnee, Okla., decreased during the hypothetical drought by as much as 28 percent, and the mean change in streamflow decreased as much as 16 percent. Streamflow at the Shawnee streamflow-gaging station did not recover to nondrought conditions until about 3 years after the simulated drought ended, during the relatively wet year of 2007.

Oklahoma↗

Geology and water resources of Owens Valley, California

Owens Valley, a long, narrow valley located along the east flank of the Sierra Nevada in east-central California, is the main source of water for the city of Los Angeles. The city diverts most of the surface water in the valley into the Owens River-Los Angeles Aqueduct system, which transports the water more than 200 miles south to areas of distribution and use. Additionally, ground water is pumped or flows from wells to supplement the surface-water diversions to the river-aqueduct system. Pumpage from wells needed to supplement water export has increased since 1970, when a second aqueduct was put into service, and local concerns have been expressed that the increased pumpage may have had a detrimental effect on the environment and the indigenous alkaline scrub and meadow plant communities in the valley. The scrub and meadow communities depend on soil moisture derived from precipitation and the unconfined part of a multilayered aquifer system. This report, which describes the hydrogeology of the aquifer system and the water resources of the valley, is one in a series designed to (1) evaluate the effects that ground-water pumping has on scrub and meadow communities and (2) appraise alternative strategies to mitigate any adverse effects caused by pumping. Two principal topographic features are the surface expression of the geologic framework--the high, prominent mountains on the east and west sides of the valley and the long, narrow intermountain valley floor. The mountains are composed of sedimentary, granitic, and metamorphic rocks, mantled in part by volcanic rocks as well as by glacial, talus, and fluvial deposits. The valley floor is underlain by valley fill that consists of unconsolidated to moderately consolidated alluvial fan, glacial and talus, and fluvial and lacustrine deposits. The valley fill also includes interlayered recent volcanic flows and pyroclastic rocks. The bedrock surface beneath the valley fill is a narrow, steep-sided graben that is structurally separated into the Bishop Basin to the north and the Owens Lake Basin to the south. These two structural basins are separated by (1) a bedrock high that is the upper bedrock block of an east-west normal fault, (2) a horst block of bedrock (the Poverty Hills), and (3) Quaternary basalt flows and cinder cones that intercalate and intrude the sedimentary deposits of the valley fill. The resulting structural separation of the basins allowed separate development of fluvial and lacustrine depositional systems in each basin. Nearly all the ground water in Owens Valley flows through and is stored in the saturated valley fill. The bedrock, which surrounds and underlies the valley fill, is virtually impermeable. Three hydrogeologic units compose the valley-fill aquifer system, a defined subdivision of the ground-water system, and a fourth represents the valley fill below the aquifer system and above the bedrock. The aquifer system is divided into horizontal hydrogeologic units on the basis of either (1) uniform hydrologic characteristics of a specific lithologic layer or (2) distribution of the vertical hydraulic head. Hydrogeologic unit 1 is the upper unit and represents the unconfined part of the system, hydrogeologic unit 2 represents the confining unit (or units), and hydrogeologic unit 3 represents the confined part of the aquifer system. Hydrogeologic unit 4 represents the deep part of the ground-water system and lies below the aquifer system. Hydrogeologic unit 4 transmits or stores much less water than hydrogeologic unit 3 and represents either a moderately consolidated valley fill or a geologic unit in the valley fill defined on the basis of geophysical data. Nearly all the recharge to the aquifer system is from infiltration of runoff from snowmelt and rainfall on the Sierra Nevada. In contrast, little recharge occurs to the system by runoff from the White and Inyo Mountains or from direct precipitation on the valley floor. Ground water flows from the margins of the valley towards the center of the valley; the ground water then flows south to the terminus of the system at Owens (dry) Lake. Ground water flows south from Bishop Basin to Owens Lake Basin through the narrows that constrict the flow opposite Poverty Hills. The aquifer system in the northern half of Owens Lake Basin is divided into east and west halves by the barrier effect caused by the Owens Valley fault. Discharge from the aquifer system is primarily by pumpage and evapotranspiration, and to a lesser extent by flowing wells, springs, underflow, and leakage to the Owens River-Los Angeles aqueduct system. Withdrawals from pumped or flowing wells is the largest component of discharge and accounts for about 50 percent of the outflow from the system. Transpiration by scrub and meadow plant communities, and to a lesser extent by irrigated alfalfa pasture, accounts for about 40 percent of the system's discharge. Natural hydraulic conductivity ranges from less than 400 to about 12,000 feet per day in the basalt flows, the more permeable material in the aquifer system. Where the basalts are fractured by explosives and drilling techniques, actual transmissivities can be greater than 1,000,000 feet squared per day. Hydraulic conductivities in sedimentary deposits of the aquifer system range from less than a few feet per day in lacustrine clays to more than 300 feet per day in gravel stringers and beach deposits in the transition zone between alluvial fan deposits and fluvial and lacustrine deposits. Degree of confinement in the aquifer system generally increases to the south and east in both the Bishop and Owens Lake Basins. The vertical hydraulic gradient across hydrogeologic unit 2 and confining beds in hydrogeologic units 1 and 3 is a function of (1) the asymmetric recharge and hydraulic head created by the dominant recharge from Sierra Nevada runoff and (2) the areal extent and thickness of the confining beds. Although most of the pumpage is from hydrogeologic unit 3, some coincident drawdown has been recorded in nonpumped wells that tap unit 1. Drawdown in hydrogeologic unit 1 is a function of changes in (1) lateral flow through hydrogeologic unit 1, (2) upward flow of ground water through the confining beds, (3) downward leakage of water from hydrogeologic unit 1 to unit 3 through wells, (4) direct withdrawal from well intervals open to hydrogeologic unit 1, and (5) increased evapotranspiration. The water in the aquifer system is generally of excellent quality for public supply and irrigation, with the exception of water stored in thick sequences of lacustrine silts and clays near Owens Lake. The water is principally a calcium bicarbonate type and dissolved-solids concentrations range from approximately 104 to 325 milligrams per liter. Water in the lacustrine sediments of Owens (dry) Lake is a sodium bicarbonate type and dissolved-solids concentrations are about 5,400 milligrams per liter.

California↗

Using cosmogenic nuclides to contrast rates of erosion and sediment yield in a semi-arid, arroyo-dominated landscape, Rio Puerco Basin, New Mexico

Analysis of in-situ-produced 10Be and 26Al in 52 fluvial sediment samples shows that millennial-scale rates of erosion vary widely (7 to 366 m Ma-1) through the lithologically and topographically complex Rio Puerco Basin of northern New Mexico. Using isotopic analysis of both headwater and downstream samples, we determined that the semi-arid, Rio Puerco Basin is eroding, on average, about 100 m Ma-1. This rapid rate of erosion is consistent with estimates made using other techniques and is likely to result from a combination of easily eroded lithologies, sparse vegetation, and monsoon-dominated rainfall. Data from 331 stream water samples collected by the US Geological Survey between 1960 and 1995 are consistent with basin-wide, average chemical denudation rates of only about 1??4 m Ma-1; thus, the erosion rates we calculate may be considered rates of sediment generation because physical weathering accounts for almost 99 per cent of mass loss. The isotopic data reveal that sediment is generally well mixed downstream with the area-weighted average sediment generation rate for 16 headwater samples (234 ton km-2 a-1 for basin area 170 to 1169 km2) matching well that estimated from a single sample collected far downstream (238 ton km-2 a-1, basin area = 14 225 km2). A series of 15 samples, collected from an arroyo wall and representing deposition through the late Holocene, indicates that 10Be concentration in sediment delivered by the fluvial system has not changed appreciably over the last 1200 years despite at least two cycles of arroyo cutting and filling. Other samples (n = 21) were collected along the drainage network. Rio Puerco erosion rates scale directly with a variety of metrics describing vegetation, precipitation, and rock erodibility. Using the headwater basins for calibration, the erosion rates for both the downstream samples and also the data set as a whole, are best modelled by considering a combination of relief and vegetation metrics, both of which co-vary with precipitation and erodibility as inferred from lithology. On average, contemporary sediment yields, determined by monitoring suspended-sediment discharge, exceed cosmogenically determined millennial-scale erosion rates by nearly a factor of two. This discrepancy, between short-term rates of sediment yield and long-term rates of erosion, suggests that more sediment is currently being exported from the basin than is being produced. Because the failure of incised channel walls and the head cutting of arroyo complexes appear to be the main sources of channel sediment today, this incongruence between rates of sediment supply and sediment yield is likely to be transitory, reflecting the current states of the arroyo cycle and perhaps the influence of current or past land-use patterns. Copyright ?? 2005 John Wiley & Sons, Ltd.

Earth Surface Processes and Landforms↗

Design of a shallow ground-water network to monitor agricultural chemicals, Lake Wales Ridge, central Florida

Extensive agricultural land use and dynamic hydraulic connections between the land surface and ground-water resources render many of Florida's aquifers vulnerable to chemical contamination. In these areas, there is a need to monitor shallow ground water for agricultural chemicals to evaluate potential migration of the chemicals to the subsurface and to deeper aquifers, and to assess the effects of agricultural practices on ground-water quality. Historically, efforts to monitor water-table aquifers have been minimal compared with monitoring deeper aquifers used for municipal drinking-water supply. The Florida Department of Agriculture and Consumer Services has proposed the establishment of long-term, water-quality networks to monitor shallow ground water in agricultural areas of Florida. These networks will bridge the current information gap between local-scale short-duration field studies, which are required for pesticide licensing, and statewide monitoring of deeper ground-water resources. Laboratory (experimental) and short-term field evaluations of the transport of agricultural chemicals to the subsurface contain some degree of uncertainty due to widespread variations in environmental conditions. The proposed regional-scale shallow ground-water networks will serve as "early warning" networks to avoid contamination of ground-water resources. The objectives of these shallow ground-water networks are to provide information for early detection of pesticides and nitrate in the subsurface, and for evaluating temporal trends in concentrations in relation to changes in land use and agricultural practices. Information on trends is important to the U.S. Environmental Protection Agency, State agencies, and agricultural managers, and is an integral part of Florida's Pesticide Management Plan, which evaluates the effectiveness of implementing "best management practices" to minimize environmental impacts. Lake Wales Ridge in Polk and Highlands Counties has been identified by the Florida Department of Agriculture and Consumer Services as a pilot study area for implementing a regional shallow ground-water network. Information gained from designing the Lake Wales Ridge Network will be used to develop plans as similar networks are considered for other areas of Florida. Lake Wales Ridge is heavily utilized for citrus production, but the well-drained, clean sands underlying the region, coupled with a hydraulic connection with deeper karst formations, render the area extremely vulnerable to the transport and migration of agricultural chemicals within the subsurface. The network design for Lake Wales Ridge was developed on the basis of the objectives of monitoring, factors affecting aquifer vulnerability, probability (statistical) sampling theory, and specifications regarding well construction, sampling, and laboratory analysis. The network will provide estimates of regional pesticide and nitrate concentrations in targeted ground water, monitor for trends over time, and will provide baseline, regional-scale information. The areas targeted for sampling were citrus groves located on soils classified as vulnerable to leaching of agricultural chemicals, and ground water in close proximity to the water table. The duration of monitoring is planned to be long term. Chemicals targeted for monitoring include selected nutrients, trace elements, and pesticides used for citrus management. Standardization of field sampling and laboratory methods, and construction of a computerized data-base repository will be important for long-term utility of the data. The proposed network design consists of existing and newly drilled wells. Existing wells were evaluated as candidates for the network based on location and well-construction specifications for water-quality monitoring. Using probability selection techniques, new well locations were selected in areas that lacked adequate existing wells. Implementation of the network will occur in phases. Quarterly sampling of 13 Phase I wells commenced in April 1999. An additional 11 Phase II wells were added to the network in April 2000. Proposed locations have been identified for 8 Phase III wells. Future topics for study include a regional water-quality analysis, a network design evaluation, and consideration of local-scale studies needed to support regional monitoring efforts. An analysis of historical water-quality data for Lake Wales Ridge, including Phase I wells, is needed for the purpose of examining trends, providing a preliminary description of spatial and short-term variability in water quality, and evaluating the proposed network design for cost effectively meeting program objectives. Local-scale, process-oriented studies are needed to provide important information for evaluating regional water-quality data. Coordination between individuals conducting such local-scale studies and those interpreting regional sampling results will be considered as a component of the network design, and will require coordination between multiple agencies.

Florida↗

Protocol for Monitoring Fish Assemblages in Pacific Northwest National Parks

Rivers and streams that drain from Olympic, Mount Rainier, and North Cascades National Parks are among the most protected corridors in the lower 48 States, and represent some of the largest tracts of contiguous, undisturbed habitat throughout the range of several key fish species of the Pacific Northwest. These watersheds are of high regional importance as freshwater habitat sanctuaries for native fish, where habitat conditions are characterized as having little to no disturbance from development, channelization, impervious surfaces, roads, diversions, or hydroelectric projects. Fishery resources are of high ecological and cultural importance in Pacific Northwest National Parks, and significantly contribute to economically important recreational, commercial, and tribal fisheries. This protocol describes procedures to monitor trends in fish assemblages, fish abundance, and water temperature in eight rivers and five wadeable streams in Olympic National Park during summer months, and is based on 4 years of field testing. Fish assemblages link freshwater, marine, and terrestrial ecosystems. They also serve as focal resources of national parks and are excellent indicators of ecological conditions of rivers and streams. Despite the vital importance of native anadromous and resident fish populations, there is no existing monitoring program for fish assemblages in the North Coast and Cascades Network. Specific monitoring objectives of this protocol are to determine seasonal and annual trends in: (1) fish species composition, (2) timing of migration of adult fish, (3) relative abundance, (4) age and size structure, (5) extent of non-native and hatchery fish, and (6) water temperature. To detect seasonal and annual trends in fish assemblages in reference sites, we rely on repeated and consistent annual sampling at each monitoring site. The general rationale for the repeated sampling of reference sites is to ensure that we account for the high interannual variability in fish movements and abundances in rivers. One underlying assumption is that the monitoring program is designed in perpetuity, and consequently our capability to detect trends substantially increases with time. The protocol describes sampling designs, methods, training procedures, safety considerations, data management, data analysis, and reporting. The allocation of sampling effort represents a balance between ecological considerations, a sound monitoring approach, and practical limitations caused by logistical constraints and a limited annual budget of $55,000. The widespread declines of native fish species in western North America highlights the importance and urgency of understanding trends in fish assemblages from undisturbed habitats. Seasonal and annual trends in fish assemblages will provide insights at the individual, population, and assemblage level. This protocol will allow managers to detect increases and decreases in abundance of priority management species, and occurrence of non-native, hatchery, and federally listed fish. The detection of trends in fish assemblages will allow for specific management actions that may include: implementation of more appropriate fishing regulations, evaluation of existing hatchery releases, control of non-native fish species, and prioritization of habitat restoration projects. Dissemination and communication of scientific findings on North Coast and Cascades Network fish assemblages will be a core product of this protocol, which will have much relevance to decision makers, park visitors, researchers, and educators.

Techniques and Methods↗

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↗

Ground-water, surface-water, and water-chemistry data, Black Mesa Area, northeastern Arizona — 2000–2001, and performance and sensitivity of the 1988 USGS numerical model of the N aquifer

The N aquifer is the major source of water in the 5,400-square-mile area of Black Mesa in northeastern Arizona. Availability of water is an important issue in this area because of continued industrial and municipal use, a growing population, and precipitation of about 6 to 14 inches per year. The monitoring program in Black Mesa has been operating since 1971 and is designed to determine the long-term effects of ground-water withdrawals from the N aquifer for industrial and municipal uses. The monitoring program includes measurements of (1) ground-water pumping, (2) ground-water levels, (3) spring discharge, (4) surface-water discharge, and (5) ground-water chemistry. In 2000, total ground-water withdrawals were 7,740 acre-feet, industrial use was 4,490 acre-feet, and municipal use was 3,250 acre-feet. From 1999 to 2000, total withdrawals increased by 9 percent, industrial use increased by 7 percent, and municipal use increased by 12 percent. From 1999 to 2001, water levels declined in 10 of 15 wells in the unconfined part of the aquifer, and the median change was -0.4 foot. Water levels declined in 8 of 16 wells in the confined part of the aquifer, and the median change was -0.2 foot. From the prestress period (prior to 1965) to 2001, the median water-level change for 33 wells was -17.2 feet. Median water-level changes were -1.2 feet for 15 wells in the unconfined part of the aquifer and -31.0 feet for 18 wells in the confined part. Discharges were measured once in 1999 and once in 2001 at four springs. Discharges decreased by 5 percent and 33 percent at two springs and increased by 3 percent and 81 percent at two springs. For about the past 10 years, discharges did not significantly change in Burro Spring, the unnamed spring near Dennehotso, and Moenkopi School Spring. The record of discharge from a consistent measuring point for Pasture Canyon Spring is too short for statistical analysis of trends. Continuous records of surface-water discharge have been collected from July 1976 to 2000 at Moenkopi Wash, July 1996 to 2000 at Laguna Creek, June 1993 to 2000 at Dinnebito Wash, and April 1994 to 2000 at Polacca Wash. Median flows for November, December, January, and February of each water year were used as an index of ground-water discharge to those streams. There is no significant trend in the median winter flows for Moenkopi Wash from 1977 to 2000. The records for the other three streams are too short for a statistical analysis of trends. The median winter flows for Dinnebito Wash and Polacca Wash, however, appear to have decreased during the last 6 years. There is no apparent trend in the median winter flows for Laguna Creek since 1997. In 2001, water samples were collected from 12 wells and 4 springs and analyzed for selected chemical constituents. Dissolved-solids concentrations ranged from 102 to 628 milligrams per liter. Water samples from 9 of the wells and from the 4 springs had less than 350 milligrams per liter of dissolved solids. Water-chemistry data with sufficient years of record for a statistical analysis of trends over time are available from 7 wells and 4 springs. From about the mid-1980s or early 1990s to 2001 there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from 6 of the 7 wells. The concentration of one tested constituent (dissolved solids) in samples from Rocky Ridge PM3 significantly increased from 1990 to 2001. From the late 1980s to 2001, there are no significant trends in the concentrations of dissolved solids, chloride, and sulfate in water samples from Burro Spring, the unnamed spring near Dennehotso, and Pasture Canyon Spring. From 1987 to 2001, concentrations of chloride and sulfate significantly increased in water samples from Moenkopi School Spring and concentrations of dissolved solids did not significantly change. The performance and sensitivity of the 1988 USGS numerical model of the N aquifer were analyzed. The overall performance of the model in steady-state conditions is reasonable for residuals of heads (difference between observed and simulated steady-state heads); 80 percent of the absolute values of residuals are less than 38 feet. Simulated flows are about 40 percent different than estimated flows at two of three discharge areas; however, this comparison is only a rough approximation of performance because the accuracy of the estimated steady-state flows is uncertain. The overall performance of the model for transient conditions is fair for residuals of changes in head (difference between observed and simulated changes in head from steady state to 1999); 80 percent of the absolute values of residuals are less than 31 feet. The model is biased in two areas. In the Tuba City area, simulated changes in head are more negative than observed changes in head; all six residuals are positive, and three residuals are between 75 and 155 feet. In the confined area of the aquifer, observed changes in head are more negative than simulated changes in head; 12 of the 17 residuals are negative, and 8 residuals are between -57 and -20 feet. Analysis of model sensitivity indicates that recharge, transmissivity, and storage coefficient are the most important parameters for estimating heads, changes in heads, and flows. A strong correlation between recharge and transmissivity and a lack of independent and reliable estimates of recharge, transmissivity, and discharge create a uniqueness problem in model calibration. Several models could be constructed and calibrated with different values of recharge or transmissivity and still have similar fits to the observed data. Information from recent data and studies and more advanced modeling techniques could be used to develop a more representative and less uncertain model. Future data collection and studies should focus on obtaining a better definition of recharge, discharge, transmissivity, and storage coefficient.

Arizona↗

A comparison of runoff quantity and quality from two small basins undergoing implementation of conventional- and low-impact-development (LID) strategies: Cross Plains, Wisconsin, water years 1999–2005

Environmental managers are often faced with the task of designing strategies to accommodate development while minimizing adverse environmental impacts. Low-impact development (LID) is one such strategy that attempts to mitigate environmental degradation commonly associated with impervious surfaces. The U.S. Geological Survey, in cooperation with the Wisconsin Department of Natural Resources, studied two residential basins in Cross Plains, Wis., during water years 1999–2005. A paired-basin study design was used to compare runoff quantity and quality from the two basins, one of which was developed in a conventional way and the other was developed with LID. The conventional-developed basin (herein called “conventional basin”) consisted of curb and gutter, 40-foot street widths, and a fully connected stormwater-conveyance system. The LID basin consisted of grassed swales, reduced impervious area (32-foot street widths), street inlets draining to grass swales, a detention pond, and an infiltration basin. Data collected in the LID basin represented predevelopment through near-complete build-out conditions. Smaller, more frequent precipitation events that produced stormwater discharge from the conventional basin were retained in the LID basin. Only six events with precipitation depths less than or equal to 0.4 inch produced measurable discharge from the LID basin. Of these six events, five occurred during winter months when underlying soils are commonly frozen, and one was likely a result of saturated soil from a preceding storm. In the conventional basin, the number of discharge events, using the same threshold of precipitation depth, was 180, with nearly one-half of those resulting from precipitation depths less than 0.2 inch. Precipitation events capable of producing appreciable discharge in the LID basin were typically those of high intensity or precipitation depth or those that occurred after soils were already saturated. Total annual discharge volume measured from the conventional basin ranged from 1.3 to 9.2 times that from the LID basin. Development of the LID basin did not appreciably alter the hydrologic response to precipitation characterized during predevelopment conditions. Ninety-five percent or more of precipitation in the LID basin was retained during each year of construction from predevelopment through near-complete build-out, surpassing the 90-percent benchmark established for new development by the Wisconsin Department of Natural Resources. The amount of precipitation retained in the conventional basin did not exceed 94 percent and fell below the 90-percent standard 2 of the 6 years monitored. Much of the runoff in the LID basin was retained by an infiltration basin, the largest control structure used to mitigate storm-runoff quantity and quality. The infiltration basin also was the last best-management practice (BMP) used to treat runoff before it left the LID basin as discharge. From May 25, 2002, to September 30, 2005, only 24 of 155 precipitation events exceeded the retention/ infiltrative capacity of the infiltration basin. The overall reduction in runoff volume from these few events was 51 percent. The effectiveness of the infiltration basin decreased as precipitation intensities exceeded 0.5 inch per hour. Annual loads were estimated to characterize the overall effectiveness of low-impact design practices for mitigating delivery of total solids, total suspended solids, and total phosphorus. Annual loads of these three constituents were greater in the LID basin than in the conventional basin in 2000 and 2004. Seventy percent or more of all constituent annual loads were associated with two discharge events in 2000, and a single discharge event produced 50 percent or more of constituent annual loads in 2004. Each of these discharge events was associated with considerable precipitation depths and (or) intensities, ranging from 4.89 to 6.21 inches and from 1.13 to 1.2 inches per hour, respectively. These same storms did not contribute as much of the annual load in the conventional basin. With large storms and saturated soils, the ability of low-impact design techniques to reduce runoff, and thus constituent loads, can be greatly diminished. For both the LID and conventional basins, the temperature of runoff was largely affected by ambient air temperatures. However, the temperature of discharge from the LID basin increased upon runoff cessation. This increase is likely due to solar heating of water that is temporarily stored in the detention pond and infiltration basin.

Wisconsin↗