Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Agriculture”

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 1,387 records · Page 77Linked to original sources

Management of conservation reserve program grasslands to meet wildlife habitat objectives

Numerous studies document environmental and social benefits of the Conservation Reserve Program (CRP). This report offers a synopsis of findings regarding effects of establishing CRP conservation practices on the quality and distribution of wildlife habitat in agricultural landscapes. On individual farms, year-round provision of wildlife habitat by the CRP may appear relatively insignificant. However, considered from multi-farm to National scales, such improvements in habitat and wildlife response have proven to be extensive and profound. Because CRP acres historically have been dominated by plantings of introduced and native grasses, this report focuses on issues pertaining to wildlife response to grass-dominated conservation practices. While the majority of CRP acres have been concentrated largely in the Great Plains and Corn Belt regions, 47 states (and Puerto Rico) have participated, resulting in measurable environmental benefits throughout the United States. Numerous investigations of habitat use by a wide range of wildlife species reveal that periodic management of CRP lands can enhance benefits through and beyond a typical 10 year general CRP contract. Over its 28-year existence, the CRP has evolved into an effective integration of conservation and agricultural policies targeting fragile and environmentally-valuable lands. Landowners with fields enrolled in the CRP often are the first to observe improvement in the landscape, greater numbers and kinds of wildlife, cleaner water and air, less erosion, and they have the satisfaction of seeing fragile lands serve better purposes. There is persistent concern that improvement seen in wildlife habitat and other environmental profits delivered by the CRP are ephemeral and last only as long as funding supports the existence of the program and its vegetative cover is properly managed. An involved American population will continue to expect governmental policies to enhance long-term protection of natural resources and public health. Recent investigations furnish evidence that the collective economic value of environmental benefits delivered by the CRP likely exceed program costs. The mounting significance placed on environmentally-responsible land management is based in part on public recognition that social, aesthetic, and recreational values enhance the traditional uses of agricultural land.

Scientific Investigations Report↗

Streamflow and estimated loads of phosphorus and dissolved and suspended solids from selected tributaries to Lake Ontario, New York, water years 2012–14

This report presents results of the evaluation and interpretation of hydrologic and water-quality data collected as part of a cooperative program between the U.S. Geological Survey and the U.S. Environmental Protection Agency. Streamflow, phosphorus, and solids dissolved and suspended in stream water were the focus of monitoring by the U.S. Geological Survey at 10 sites on 9 selected tributaries to Lake Ontario during the period from October 2011 through September 2014. Streamflow yields (flow per unit area) were the highest from the Salmon River Basin due to sustained yields from the Tug Hill aquifer. The Eighteenmile Creek streamflow yields also were high as a result of sustained base flow contributions from a dam just upstream of the U.S. Geological Survey monitoring station at Burt. The lowest streamflow yields were measured in the Honeoye Creek Basin, which reflects a decrease in flow because of withdrawals from Canadice and Hemlock Lakes for the water supply of the City of Rochester. The Eighteenmile Creek and Oak Orchard Creek Basins had relatively high yields due in part to groundwater contributions from the Niagara Escarpment and seasonal releases from the New York State Barge Canal. Annual constituent yields (load per unit area) of suspended solids, phosphorus, orthophosphate, and dissolved solids were computed to assess the relative contributions and allow direct comparison of loads among the monitored basins. High yields of total suspended solids were attributed to agricultural land use in highly erodible soils at all sites. The Genesee River, Irondequoit Creek, and Honeoye Creek had the highest concentrations and largest mean yields of total suspended solids (165 short tons per square mile [t/mi 2 ], 184 t/mi 2 , and 89.7 t/mi 2 , respectively) of the study sites. Samples from Eighteenmile Creek, Oak Orchard Creek at Kenyonville, and Irondequoit Creek had the highest concentrations and largest mean yields of phosphorus (0.27 t/mi 2 , 0.26 t/mi 2 , and 0.20 t/mi 2 , respectively) and orthophosphate (0.17 t/mi 2 , 0.13 t/mi 2 , and 0.04 t/mi 2 , respectively) of the study sites. These results were attributed to a combination of sources, including discharges from wastewater treatment plants, diversions from the New York State Barge Canal, and manure and fertilizers applied to agricultural land. Yields of phosphorus also were high in the Genesee River Basin (0.17 t/mi 2 ) and were presumably associated with nutrient and sediment transport from agricultural land and from streambank erosion. The Salmon and Black Rivers, which drain a substantial amount of forested land and are influenced by large groundwater discharges, had the lowest concentrations and yields of phosphorus and orthophosphate of the study sites. Mean annual yields of dissolved solids were the highest in Irondequoit Creek due to a high percentage of urbanized area in the basin and in Oak Orchard Creek at Kenyonville and in Eighteenmile Creek due to groundwater contributions from the Niagara Escarpment. High yields of dissolved solids of 840 t/mi 2 , 829 t/mi 2 , and 715 t/mi 2 , respectively, from these basins can be attributed to seasonal chloride yields associated with use of road deicing salts. The Niagara Escarpment can produce large amounts of dissolved solids from the dissolution of minerals (a continual process reflected in base flow samples). Groundwater inflows in the Salmon River have very low concentrations of dissolved solids due to minimal bedrock interaction along the Tug Hill Plateau and discharge from the Tug Hill sand and gravel aquifer, which has minimal mineralization.

New York↗

Multiple-source tracking: Investigating sources of pathogens, nutrients, and sediment in the Upper Little River Basin, Kentucky, water years 2013–14

The South Fork Little River (SFLR) and the North Fork Little River (NFLR) are two major headwater tributaries that flow into the Little River just south of Hopkinsville, Kentucky. Both tributaries are included in those water bodies in Kentucky and across the Nation that have been reported with declining water quality. Each tributary has been listed by the Kentucky Energy and Environment Cabinet—Kentucky Division of Water in the 303(d) List of Waters for Kentucky Report to Congress as impaired by nutrients, pathogens, and sediment for contact recreation from point and nonpoint sources since 2002. In 2009, the Kentucky Energy and Environment Cabinet—Kentucky Division of Water developed a pathogen total maximum daily load (TMDL) for the Little River Basin including the SFLR and NFLR Basins. Future nutrient and suspended-sediment TMDLs are planned once nutrient criteria and suspended-sediment protocols have been developed for Kentucky. In this study, different approaches were used to identify potential sources of fecal-indicator bacteria (FIB), nitrate, and suspended sediment; to inform the TMDL process; and to aid in the implementation of effective watershed-management activities. The main focus of source identification was in the SFLR Basin. To begin understanding the potential sources of fecal contamination, samples were collected at 19 sites for densities of FIB ( E. coli ) in water and fluvial sediment and at 11 sites for Bacteroidales genetic markers (General AllBac, human HF183, ruminant BoBac, canid BacCan, and waterfowl GFD) during the recreational season (May through October) in 2013 and 2014. Results indicated 34 percent of all E. coli water samples ( n =227 samples) did not meet the U.S. Environmental Protection Agency 2012 recommended national criteria for primary recreational waters. No criterion currently exists for E. coli in fluvial sediment. By use of the Spearman’s rank correlation test, densities of FIB in fluvial sediments were observed to have a statistically significant positive correlation with drainage area. As drainage area increased, so did the densities of FIB in the fluvial sediments. There was no statistically significant correlation between drainage area and FIB in water. The human-associated marker (HF183) was found above the detection limit in 26 percent of the samples ( n =120 samples); a higher proportion of positive samples was in the NFLR Basin. The ruminant-associated marker (BoBac) was above the detection limit in 65 percent of samples; a higher proportion of positive samples was in the headwaters of the SFLR Basin. Nutrient yields differed between the SFLR and NFLR Basins. Comparatively, the SFLR Basin produced the largest estimated mean yields of total nitrogen (16,000 pounds per year per square mile (lb/yr/mi 2 ) and nitrite plus nitrate nitrogen (12,500 lb/yr/mi 2 ), and the NFLR Basin produced the largest estimated mean yields of ammonia plus organic nitrogen (4,700 lb/yr/mi 2 ), total phosphorus (1,100 lb/yr/mi 2 ), and orthophosphorus (590 lb/yr/mi 2 ). Nitrate sources in surface water were assessed in both basins using dual-nitrate isotope (nitrogen and oxygen) ratios. Data from the different land uses in the SFLR Basin showed differences in nitrate concentrations and overlapping, but moderately distinct, isotopic signatures. Predominantly forested sites consistently had low nitrate concentrations (median = 0.233 milligrams per liter) with minimal variability, and agricultural sites had the highest nitrate concentrations (median = 7.55 milligrams per liter) with the greatest variability. The median nitrate concentration for sites with mixed land use was 2.66 milligrams per liter. Dual-isotope data for forested sites plotted within ranges characteristic of soil-derived nitrate with possible but minimal influence from recycled atmospheric nitrate. Ranges of dual-isotope data for sites with agricultural and mixed land uses were characteristic of possible mixtures of chemical fertilizer, soil-derived nitrate, and manure and septic wastes. In the NFLR Basin, a positive linear relation was observed between nitrate concentrations and nitrogen isotope ratios (δ 15 N NO3 ) (R 2 =0.56; p -value <0.001) that potentially suggests the NFLR Basin has a higher proportion of δ 15 N NO3 -enriched sources, such as manure and sewage. However, mixing of other nitrate-derived sources cannot be excluded, because many values of δ 15 N NO3 and concentrations of nitrate showed minimal variation and plotted within dual-nitrate isotope ranges characteristic of fertilizer and soil-derived nitrate sources. A sediment-fingerprinting approach was used to quantify the relative contribution of four upland sources in the SFLR Basin (agricultural, pasture, riparian/forest, and streambank) to understand how land management affects suspended-sediment concentration. Carbon isotope ratios (δ 13 C), together with calcium and carbon concentrations, were the best indicators of sediment source; the uncertainty was less than 11 percent. Fine-sediment samples collected at the SFLR Basin outlet indicated streambanks as the largest source of the fine sediment to the stream followed by cropland and riparian/forest-source areas, respectively; pasture was a minor contributing source. Streambanks and cropland were essentially equal contributors of fine sediment at the NFLR Basin outlet.

Kentucky↗

Hydrogeologic characterization, groundwater chemistry, and vulnerability assessment, Ute Mountain Ute Reservation, Colorado and Utah

The U.S. Geological Survey, in cooperation with the Ute Mountain Ute Tribe (UMUT), initiated a study in 2016 to increase understanding of the hydrogeology and chemistry of groundwater within select areas of the Ute Mountain Ute Reservation (UMUR) in Colorado and Utah, identify vulnerabilities to the system and other natural resources, and outline information needs to aid in the understanding and protection of groundwater resources. The results presented for this study can be used to support the UMUT’s goal of protecting their vital groundwater resources on the UMUR. Hydrogeologic conditions were characterized for the surficial aquifer contained in Quaternary-age unconsolidated surficial deposits and the Dakota aquifer contained in the Cretaceous-age Dakota Sandstone. In the surficial aquifer, median depth to water ranges from about 5.4 to 17.2 feet below land surface in the Farm and Ranch Enterprise area and 11 to 34 feet below land surface in the Towaoc area, and the water table slopes generally southwest or south. A map of depth to the top of the Dakota Sandstone was constructed from existing well data. Depths range from zero in outcrop areas to more than 3,000 feet below land surface on mesas in the southeastern part of the UMUR. Groundwater-chemistry data were collected by the UMUT from 13 springs and 31 wells from 1996 through 2017. Specific conductance was much lower for samples from springs than from wells; median values were 512 and 6,024 microsiemens per centimeter at 25 degrees Celsius, respectively. Spring samples were well oxygenated. A few well samples were anoxic (dissolved oxygen concentrations less than 0.5 milligrams per liter [mg/L]), indicating reducing conditions in the aquifer. About 75 percent of spring samples had fresh water (total dissolved solids concentrations less than 1,000 mg/L), and about 85 percent of well samples had brackish or highly saline water (total dissolved solids concentrations greater than 1,000 mg/L). Water type for springs on the Ute Mountains was calcium bicarbonate. Lower-altitude springs had a calcium-sulfate water type. Most well samples had sodium as the dominant cation, and sulfate, bicarbonate, and chloride as the dominant anions. Fluoride concentrations in about 45 percent of well samples were greater than an agricultural-use standard of 2 mg/L. Nitrate plus nitrite concentrations in most spring and well samples were less than about 1.6 mg/L per liter. Concentrations in samples from wells in the irrigated agricultural area were elevated; the maximum concentration was 78.5 mg/L. About one-half of the trace-element samples had concentrations that were less than laboratory reporting limits. Only aluminum, arsenic, and selenium in spring samples, and boron and selenium in well samples, were detected at concentrations greater than surface-water standards or water-quality standards for agricultural use of groundwater. Only three organic compounds, the pesticides alachlor and atrazine and the volatile organic compound di(2-ethylhexyl) phthalate, were detected in well samples. The Escherichia coli bacteria was detected in 47 and 23 percent of samples from wells and springs, respectively. The E. coli detections included samples from three culturally significant springs, which did not meet the UMUT cultural-use standard of total absence of E. coli. Tritium and carbon-14 were the primary environmental tracers used for interpreting groundwater ages for Lopez 2 Spring and five wells (AP–1, 5000 Block, Cottonwood Spring, Goodknight, and SE Toe). Water from the AP–1 well contained a mixture of pre- and post-1950s recharge. Tritium and carbon-14 recharge ages for Lopez 2 Spring (post-1950s in age), Goodknight and SE Toe wells (pre-1950s in age), and Cottonwood Spring well (primarily pre-1950s in age) are supported by helium-4 data. The helium-4 data for the 5000 Block well are inconsistent with the tritium and carbon-14 age of pre-1950s recharge because of interference caused by high methane concentrations in the water. Springs and surficial deposits are more vulnerable to contamination from anthropogenic chemicals than deeper bedrock wells. Bedrock aquifers are vulnerable in areas where the geologic formations containing the aquifers are exposed at the land surface. Groundwater in deep bedrock aquifers is likely thousands of years old and is not currently affected by present-day land uses. Both shallow and deep groundwater are vulnerable to naturally occurring salts and minerals, such as of total dissolved solids, major ions, nitrate, and trace elements. Effects of a changing climate on water resources and other ecological characteristics of the UMUR could include changes in evapotranspiration, a decrease in snowpack, decreased aquifer recharge and flow of springs, a decrease in soil moisture, and increased occurrence of wildfires and forest mortality. Of particular interest for the UMUT are possible effects of a changing climate on medicinal and culturally important plants and springs Several information needs were identified during this study that would aid in the understanding and protection of groundwater resources on the UMUR. These include well-completion information for bedrock wells, the collection of environmental tracer data at additional wells, the addition of methane and hydrocarbon analysis to well sampling plans, and the resampling of springs and wells that were last sampled in 2002 or earlier.

Colorado↗

Spatially referenced models of streamflow and nitrogen, phosphorus, and suspended-sediment loads in the southeastern United States

Spatially Referenced Regression On Watershed attributes (SPARROW) models were applied to describe and estimate mean-annual streamflow and transport of total nitrogen (TN), total phosphorus (TP), and suspended sediment (SS) in streams and delivered to coastal waters of the southeastern United States on the basis of inputs and management practices centered near 2012, the base year of the model. Previously published TN and TP models for 2002 served as a starting point and reference for comparison. The datasets developed for the 2012 models not only represent updates of previous conditions but also incorporate new approaches for characterizing sources and transport processes that were not available for previous models. Variability in streamflow across the southeastern United States was explained as a function of precipitation adjusted for evapotranspiration, spring discharge, and municipal and domestic wastewater discharges to streams. Results from the streamflow model were used as input to the water-quality SPARROW models, and areas with large streamflow prediction errors—urban areas and karst areas—were used to provide guidance on where additional data are needed to improve routing of flow. Variability in TN transport in Southeast streams was explained by the following five sources in order of decreasing mass contribution to streams: atmospheric deposition, agricultural fertilizer, municipal wastewater, manure from livestock, and urban land. Variable rates of TN delivery from source to stream were attributed to variation among catchments in climate, soil texture, and vegetative cover, including the extent of cover crops in the watershed. Variability in TP transport in Southeast streams was explained by the following six sources in order of decreasing mass contribution to streams: parent-rock minerals, urban land, manure from livestock, municipal wastewater, agricultural fertilizer, and phosphate mining. Varying rates of TP delivery were attributed to variation in climate, soil erodibility, depth to water table, and the extent of conservation tillage practices in the watershed. Variability in SS transport in Southeast streams was explained by variable sediment export rates for different combinations of land cover and geologic setting (for upland sources of sediment) and by gains in stream power caused by longitudinal changes in channel hydraulics (for channel sources of sediment). Sediment yields for the transitional land cover (shrub, scrub, herbaceous, and barren) varied widely depending on geologic setting and on agricultural land cover. Varying rates of SS delivery, like those for TP, were attributed to variation in climate, soil erodibility, and the extent of conservation tillage practices in the watershed, as well as to areal extent of canopy land cover in the 100-meter buffer along the channel. Relatively large uncertainty, compared to the other three models, for almost all the SS source coefficients indicates the need for caution when interpreting the results from the sediment model. TN, TP, and SS inputs to streams from sources were balanced in the models with losses from physical processes in streams and reservoirs and with water withdrawals. The losses in streams and reservoirs along with withdrawals removed 35, 44, and 65 percent of the TN, TP, and SS load, respectively, that entered streams before reaching coastal waters.

Scientific Investigations Report↗

Hydrogeology, numerical simulation of groundwater flow, and effects of future water use and drought for reach 1 of the Washita River alluvial aquifer, Roger Mills and Custer Counties, western Oklahoma, 1980–2015

The Washita River alluvial aquifer is a valley-fill and terrace alluvial aquifer along the valley of the Washita River in western Oklahoma that provides a productive source of groundwater for agricultural irrigation and water supply. The Oklahoma Water Resources Board (OWRB) has designated the westernmost section of the aquifer in Roger Mills and Custer Counties, Okla., as reach 1 of the Washita River alluvial aquifer; reach 1 is the focus of this report. The OWRB issued an order on November 13, 1990, that established the maximum annual yield (MAY; 120,320 acre-feet per year [acre-ft/yr]) and equal-proportionate-share (EPS) pumping rate (2.0 acre-feet per acre per year [(acre-ft/acre)/yr]) for reach 1 of the Washita River alluvial aquifer. The MAY and EPS were based on hydrologic investigations that evaluated the effects of potential groundwater withdrawals on groundwater availability in the Washita River alluvial aquifer. Every 20 years, the OWRB is statutorily required to update the hydrologic investigation on which the MAY and EPS were based. Because 30 years have elapsed since the last order was issued, the U.S. Geological Survey, in cooperation with the OWRB, conducted a new hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the Washita River alluvial aquifer. The Washita River is the primary source of inflow to Foss Reservoir, a Bureau of Reclamation reservoir constructed in 1961 for flood control, water supply, and recreation. Foss Reservoir provides water for Bessie, Clinton, New Cordell, and Hobart, Okla. Nearly 98 percent of the total groundwater use from the Washita River alluvial aquifer during 1967 to 2015 was for irrigation; other uses of groundwater in the study area include public supply, mining, and agriculture. A hydrogeologic framework was developed for the Washita River alluvial aquifer and included the physical characteristics of the aquifer, the geologic setting, the hydraulic properties of hydrogeologic units, the potentiometric surface (water table), and groundwater-flow directions at a scale that captures the regional controls on groundwater flow. The Washita River alluvial aquifer consists of alluvium and terrace deposits that were transported primarily by water and range from clay to gravel in size. The terrace includes windblown deposits of silt size and, in some cases, contains gravel laid down at several levels along former courses of present-day rivers. A conceptual flow model is a simplified description of the aquifer system that includes hydrologic boundaries, major inflow and outflow sources of the groundwater-flow system, and a conceptual water budget with the estimated mean flows between those hydrologic boundaries. During the study period 1980–2015, mean annual groundwater withdrawals, predominantly used for agricultural irrigation, totaled 5,502 acre-ft/yr, or 14 percent of aquifer outflows. When applied across the 132-square-mile aquifer area used for modeling purposes (84,366 acres), mean annual recharge of 3.15 inches per year corresponds to a mean annual recharge volume of 22,169 acre-ft/yr, or 56 percent of aquifer inflows. The annual saturated-zone evapotranspiration outflow was 11,828 acre-ft/yr for the Washita River alluvial aquifer, or about 30 percent of aquifer outflows. For the Washita River alluvial aquifer, lateral flow was 17,157 acre-ft/yr, or 44 percent of the aquifer inflows. The conceptual flow model and hydrogeologic framework were used to conceptualize, design, and build the numerical groundwater-flow model. A numerical groundwater-flow model of the Washita River alluvial aquifer was constructed by using MODFLOW-2005. The Washita River alluvial aquifer groundwater-model grid was spatially discretized into 350-foot (ft) cells and two layers. Layer 1 represented the undifferentiated alluvium and terrace deposits of Quaternary age, and layer 2 represented the bedrock of Permian age, which was given a uniform nominal thickness of 100 ft. The groundwater-simulation period was temporally discretized into 433 monthly transient stress periods, representing January 1980 to December 2015. An initial 365-day steady-state stress period was configured to represent mean annual inflows and outflows from the Washita River alluvial aquifer for the study period. The groundwater-flow model was calibrated manually and by automated adjustment of model inputs by using PEST++. Calibration targets for the Washita River alluvial aquifer model included groundwater-level observations and reservoir-stage observations, as well as base-flow and stream-seepage estimates. Three groundwater-availability scenarios were used in the calibrated groundwater model to (1) estimate the EPS pumping rate that retains the saturated thickness that meets the minimum 20-year life of the aquifer, (2) quantify the effects of projected pumping rates on groundwater storage over a 50-year period, and (3) evaluate how projected pumping rates extended 50 years into the future and sustained hypothetical drought conditions over a 10-year period affect base flow and groundwater in storage. The results of the groundwater-availability scenarios could be used by the OWRB to reevaluate the established MAY of groundwater from the Washita River alluvial aquifer. EPS scenarios for the Washita River alluvial aquifer were run for periods of 20, 40, and 50 years. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were 1.7, 1.6, and 1.6 (acre-ft/acre)/yr, respectively. Given the aquifer area used for modeling purposes (84,366 acres), these rates correspond to annual yields of 142,579, 134,986, and 134,986 acre-ft/yr, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 281,000 acre-feet (acre-ft), or about 306,000 acre-ft (52 percent) less than the starting storage. Considering the land-surface area of the Washita River alluvial aquifer and using a specific yield of 0.12, this decrease in storage was equivalent to a mean groundwater-level decline of about 30 ft. The Washita River downstream from Foss Reservoir and most of the streams in the study area were dry at the end of the 20-year EPS scenario. Foss Reservoir stage was below the dead-pool stage of 1,597 ft after about 7 years of pumping in the 20-year EPS scenario. Four projected 50-year groundwater-use scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage in the Washita River alluvial aquifer. These four scenarios used (1) no groundwater use, (2) groundwater use at the 2015 pumping rate, (3) mean groundwater use for the simulation period, and (4) increasing groundwater use. Groundwater storage after 50 years with no groundwater use was 545,249 acre-ft, or 693 acre-ft (0.1 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean groundwater-level increase of 0.1 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 543,831 acre-ft, or 723 acre-ft (0.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage after 50 years with the mean pumping rate for the study period was 543,202 acre-ft, or 1,349 acre-ft (0.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.1 ft. Groundwater storage at the end of the 50-year period with an increasing demand groundwater-pumping rate, which was 38 percent greater than the 2015 groundwater-pumping rate, was 542,584 acre-ft, or 1,967 acre-ft (0.4 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean groundwater-level decrease of 0.2 ft. A hypothetical 10-year-drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage in the Washita River alluvial aquifer and Foss Reservoir stage and storage. To simulate the hypothetical drought, recharge in the calibrated model was reduced by 50 percent during the simulated drought period (1983–1992). Groundwater storage at the end of the drought period in December 1992 was 562,000 acre-ft, or 36,000 acre-ft (6 percent) less than the groundwater storage of the calibrated groundwater model (598,000 acre-ft). At the end of the hypothetical drought, the largest changes in saturated thickness (as great as 43.5 ft) were in the area upgradient from Foss Reservoir, particularly in the terrace at the model boundary. Substantial decreases in the Foss Reservoir stage began during the fall of 1985 in conjunction with base-flow decreases of up to 100 percent at U.S. Geological Survey streamgage 07324200 Washita River near Hammon, Okla. These lake-stage declines outpaced groundwater-level declines in the surrounding aquifer. The minimum Foss Reservoir storage simulated during the drought period was 77,954 acre-ft, which was a decrease of 46 percent from the nondrought storage.

Oklahoma↗

Areas contributing recharge to selected production wells in unconfined and confined glacial valley-fill aquifers in Chenango River Basin, New York

In the Chenango River Basin of central New York, unconfined and confined glacial valley-fill aquifers are an important source of drinking-water supplies. The risk of contaminating water withdrawn by wells that tap these aquifers might be reduced if the areas contributing recharge to the wells are delineated and these areas protected from land uses that might affect the water quality. The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation and the New York State Department of Health, began an investigation in 2019 to improve understanding of groundwater flow and delineate areas contributing recharge to 16 production wells clustered in three study areas in the basin as part of an effort to protect the source of water to these wells. Areas contributing recharge were delineated on the basis of numerical steady-state groundwater-flow models representing long-term average hydrologic conditions. In the Cortland study area, four water suppliers operate 10 production wells that withdraw a total average rate of 2,480 gallons per minute from an unconfined aquifer consisting of well-sorted sand and gravel deposits. Simulated areas contributing recharge to these wells at their average pumping rates covered a total area of 6.93 square miles. Simulated areas contributing recharge extend upgradient from the wells to upland till deposits and to groundwater divides. Some simulated areas contributing recharge include isolated areas remote from the wells. Short simulated groundwater traveltimes from recharging locations to discharging wells indicated that the wells are vulnerable to contamination from land-surface activities; 50 percent of the traveltimes were 10 years or less. Land cover in some of the areas contributing recharge included a substantial amount of urban and agriculture land use. The groundwater-flow model of the Cortland study area was calibrated to available hydrologic data by inverse modeling using nonlinear regression. The parameter variance-covariance matrix from model calibration was used to create parameter sets that reflect the uncertainty of the parameter estimates and the correlation among parameters to evaluate the uncertainty associated with the single, predicted contributing areas to the wells. This analysis led to contributing areas expressed as a probability distribution. Because of the effects of parameter uncertainty, the size of the probabilistic contributing areas was larger than the size of the single, predicted contributing area for the wells. Thus, some areas not in the single, predicted contributing area might actually be in the contributing area, including additional areas of urban and agriculture land use that have the potential to contaminate groundwater. Additional areas that might be in the contributing area included recharge originating near the pumping wells that have relatively short groundwater-flow paths and traveltimes. In each of the Greene and Cincinnatus study areas, one water supplier operates three wells that are screened near the top of the bedrock surface in a confined aquifer consisting of poorly to well-sorted sand and gravel deposits. This confined aquifer is overlain by a lacustrine confining unit of very fine sand, silt, and clay, which in turn is overlain by a thin unconfined aquifer of sand and gravel. The groundwater-flow models for these two areas were manually calibrated because of the limited hydrologic data. Simulated areas contributing recharge to the Greene study area wells covered a total area of 0.35 square mile for the average pumping rate of 170 gallons per minute. The contributing areas extended southeastward of the wells to the groundwater divide in the till uplands. The contributing areas also included remote, isolated areas on the opposite side of the Chenango River from the wells primarily in the till uplands. For the Cincinnatus study area wells, which have a low average pumping rate (34 gallons per minute), the simulated contributing areas totaled 0.06 square mile and were on the same side of the river as the wells, but they are isolated areas remote from the wells primarily in the till-covered bedrock uplands. Land cover in these contributing areas for both study areas is primarily agriculture and forested, with the contributing areas to the Greene study area wells also including some urban land uses. Because the Greene and Cincinnatus study area wells are screened relatively deep and some flow paths to the wells partly travel through the confining unit, which impedes the connection with surface sources of recharge, overall groundwater traveltimes are greater than for wells in the Cortland study area. Fifty percent of Cortland study area wells, but only 9 and 44 percent of Greene and Cincinnatus study area wells, respectively, have groundwater traveltimes of 10 years or less.

New York↗

Status and understanding of groundwater quality in the Sacramento Metropolitan Domestic-Supply Aquifer study unit, 2017—California GAMA Priority Basin Project

Groundwater quality in the Sacramento Metropolitan Domestic-Supply Aquifer study unit (SacMetro-DSA) was studied from August to November 2017 as part of the second phase of the Priority Basin Project of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is in parts of Amador, Placer, Sacramento, and Sutter Counties, and the extent of the study unit was defined by the location of three California Department of Water Resources groundwater subbasins: the North American, the South American, and the Cosumnes. The SacMetro-DSA focused on groundwater resources used for domestic drinking-water supply, which generally correspond to shallower parts of aquifer systems than those of groundwater resources used for public drinking water supply in the same area. The assessments characterized the quality of untreated groundwater, not the quality of drinking water. This study included two components: (1) a status assessment, which characterized the status of the quality of the groundwater resources used for domestic supply and (2) an understanding assessment, which evaluated the natural and human factors potentially affecting water quality in those resources. The first component of this study—the status assessment—was based on water-quality data collected from 49 sites sampled by the U.S. Geological Survey for the GAMA Priority Basin Project in 2017. The samples were analyzed for volatile organic compounds, pesticides, and naturally present inorganic constituents, such as major ions and trace elements. To provide context, concentrations of constituents measured in groundwater were compared to U.S. Environmental Protection Agency and California State Water Resources Control Board Division of Drinking Water regulatory and non-regulatory benchmarks for drinking-water quality. The status assessment used a grid-based method to estimate the proportion of the groundwater resources that had concentrations of water-quality constituents approaching or above benchmark concentrations. This method provides statistically unbiased results at the study-area scale and permits comparisons to other GAMA Priority Basin Project study areas. The second component of this study—the understanding assessment—identified the natural and human factors that potentially affect groundwater quality by evaluating land-use characteristics, groundwater age, and geochemical and hydrologic conditions of the domestic-supply aquifer and related these data to constituents identified in the status assessment for further evaluation. In the SacMetro-DSA study unit, arsenic was the only inorganic constituent detected above health-based benchmarks and was detected in 10 percent of the domestic-supply aquifer system. Inorganic constituents were detected above the non-health-based California State Water Resources Control Board—Division of Drinking Water secondary maximum contaminant levels (SMCL-CA) in 16 percent of the system. The inorganic constituents detected above the SMCL-CA were chloride, iron, manganese, and total dissolved solids (TDS). Organic constituents (volatile organic compounds and pesticides) with health-based benchmarks were not detected above health-based benchmarks; however, chloroform was detected at concentrations higher than 10 percent of the health-based benchmark (80 micrograms per liter) in 2 percent of the domestic-supply aquifer system. Of the 310 organic constituents analyzed, 16 constituents were detected; however, only bentazon and chloroform had detection frequencies greater than 10 percent. Inorganic constituents with health-based benchmarks that were evaluated in the understanding assessment included arsenic and hexavalent chromium. Arsenic and hexavalent chromium are natural constituents of aquifer sediments in the study unit and did not appear to be influenced by anthropogenic processes; rather, the presence of arsenic and hexavalent chromium appeared to be related to geochemical conditions controlled by oxidation–reduction reactions in the aquifer system. Naturally occurring inorganic constituents with SMCL-CAs evaluated in the understanding assessment were the trace elements iron and manganese, the major ion chloride, and TDS. Like arsenic and hexavalent chromium, the presence of iron and manganese was most strongly related to geochemical conditions in the aquifer system, specifically reducing conditions, which were most common near the western edge of the study unit close to the Sacramento River. Concentrations of chloride and TDS are indicators of salinity and were correlated with variables related to well location and included redox, agricultural land use, and elevation. Chloride and TDS were positively correlated to reducing conditions, and agricultural land use was negatively correlated to elevation and well depth. Observed correlations among variables were likely driven by the characteristics of the western part of the study unit, such as its higher proportion of agricultural land use and its relatively low elevation. A large portion of the western edge of the study unit is located in the center of the Sacramento Valley, defined by the location of the Sacramento River. The special-interest constituent perchlorate, also included in the understanding assessment, has natural and anthropogenic sources. Perchlorate was detected frequently and at moderate relative concentrations. In some areas of the study unit, concentrations of perchlorate were higher than what might be expected in nature; therefore, anthropogenic introduction of perchlorate or anthropogenically induced migration of native perchlorate could be occurring.

California↗

Update and recalibration of the Rio Grande Transboundary Integrated Hydrologic Model, New Mexico and Texas, United States, and northern Chihuahua, Mexico

The Rio Grande Transboundary Integrated Hydrologic Model (RGTIHM) was developed through an interagency effort between the U.S. Geological Survey and the Bureau of Reclamation to provide a tool for analyzing the hydrologic system response to the historical evolution of water use and potential changes in water supplies and demands in the Hatch Valley (also known as Rincon Valley in the study area) and Mesilla Basin, New Mexico and Texas, United States, and northern Chihuahua, Mexico. Reclamation operates the Rio Grande Project (RGP) to store and deliver surface water for irrigation and municipal use within the study area and in the El Paso Valley south of the El Paso Narrows. Biases in the RGTIHM’s simulation of streamflow and aquifer storage depletion and the availability of new estimates of historical agricultural consumptive use in the study area initiated an update and recalibration of the RGTIHM. In addition to the new estimates of historical agricultural consumptive use, updates were made to more accurately represent the natural system and included adjustments to the initial groundwater levels; streamflow rating tables; Rio Grande, canal, and drain streambed elevations; tributary streambed elevations; surface-water inflows and diversions; RGP surface-water deliveries and canal waste; on-farm efficiency; the routing of surface-water runoff within the MODFLOW Farm Process; and general head boundaries used to simulate interbasin groundwater flow. Model settings, including the assignment of hydraulic conductivity and storage properties to model layers and the MODFLOW solver package, were adjusted to improve numerical stability, and the model was recalibrated to better simulate the natural system. The updated and recalibrated RGTIHM demonstrates a robust ability to simulate the spatially and temporally variable measurements, estimates, or reports of hydraulic head, surface-water flows, agricultural pumping, RGP surface-water deliveries and canal waste, and decadal aquifer storage changes, with improvements over the previous version of the model.

Chihuahua, New Mexico, Texas↗

Estimated effects of pumping on groundwater storage and Walker River stream efficiencies in Smith and Mason Valleys, west-central Nevada

The Walker River originates in the Sierra Nevada Mountains and flows nearly 160 miles to its terminus at Walker Lake in west-central Nevada. The river provides a source of irrigation water for tens of thousands of acres of agricultural lands in California and Nevada and is the principal source of inflow to Walker Lake. Extraction of groundwater for agricultural use became prevalent in the late 1950s and early 1960s to supplement irrigation demands not met by surface-water diversions during times of drought. There is growing concern that continued groundwater withdrawals within the Walker River Basin are likely contributing to depleted streamflow of the Walker River and the long-term depletion of groundwater storage in the basin. This report documents changes in groundwater storage-volume and trends in Walker River stream efficiency, a measure of change in flow due to gaining or losing conditions, in the two largest agricultural valleys in the Walker River Basin, Smith and Mason Valleys, for a multi-decade period. Groundwater-level maps from previous studies were used for the beginning (1970) and middle (2006) points of this study. Groundwater levels measured from 1991–95 and 2016–20 were used to construct median groundwater-level maps that represented conditions in 1995 and 2020. Valley wide groundwater-level change was calculated by comparing groundwater-level maps for the periods 1970–95, 1996–2006, and 2007–20 and by observing the overall change from 1970 to 2020. Groundwater storage-volume change was calculated using groundwater-level change and previously defined specific yield values. Between 1970 and 2020, groundwater storage-volume declined 287,600 acre-feet in Smith Valley and 269,000 acre-feet in Mason Valley. Using groundwater storage-volume decline and annual groundwater pumpage rates, a maximum groundwater pumpage rate can be computed to support management of water resources. In Smith Valley, groundwater pumping in excess of 22,300 acre-feet per year would likely result in groundwater storage decline. In Mason Valley, groundwater pumping in excess of 75,200 acre-feet per year would likely result in groundwater storage decline. Stream efficiency was calculated using continuous streamflow data and monthly diversion volumes on two reaches: (1) the West Walker River in Smith Valley, from 1948 to 2020 and (2) the Walker River in Mason Valley, from 1958 to 2020. Stream efficiency during non-irrigation season in Smith and Mason Valleys declined at a statistically significant rate of 1.1 and 0.6 percent per year, respectively. Trends in stream efficiency corresponded to occurrence of prolonged drought, deviation from average annual streamflows, and total groundwater pumpage. Long-term declines in groundwater storage-volume and stream efficiency demonstrate that the alluvial aquifer system is becoming increasingly depleted, such that the river can no longer replenish groundwater storage while simultaneously balancing groundwater and surface-water withdrawals. The introduction of supplemental groundwater pumpage was intended to offset surface-water deficits during dry years; however, pumpage occurs even in years when average or above average streamflows meet surface-water demands. Reliance on supplemental groundwater pumpage has resulted in widespread groundwater storage-volume decline and decreased stream efficiency. With each successive drought cycle, the ability of Walker River to sustain streamflows and convey water downstream has diminished. Above average wet periods have a marginal and short-lived effect on rebounding the groundwater levels outside of the river corridor. Moreover, if the trend continues, each future drought cycle may further reduce groundwater supplies and that may further decrease streamflow reliability.

Nevada↗

Arsenic, chromium, uranium, and vanadium in rock, alluvium, and groundwater, Mojave River and Morongo Areas, western Mojave Desert, southern California

Trace elements within groundwater that originate from aquifer materials and pose potential public-health hazards if consumed are known as geogenic contaminants. The geogenic contaminants arsenic, chromium, and vanadium can form negatively charged ions with oxygen known as oxyanions. Uranium complexes with bicarbonate and carbonate to form negatively charged ions having aqueous chemistry similar to oxyanions. The concentrations of arsenic, chromium, uranium, and vanadium in groundwater result from the combined effects of (1) geologic abundance within aquifer materials; (2) the fraction of these elements that have weathered from and sorbed to the surfaces of mineral grains and are potentially available to groundwater; and (3) the aqueous chemistry of dissolved oxyanions in groundwater during different redox conditions and pH, both of which are affected by hydrogeology, including the length of time groundwater has been in contact with aquifer materials. Concentrations of arsenic, chromium, uranium, and vanadium were measured in samples of (1) rock, surficial alluvium, and drill cuttings using portable (handheld) X-ray fluorescence (pXRF); (2) operationally defined fractions extractable from these materials; and (3) water from wells sampled between 2000 and 2018 within the 3,500 square mile Mojave River area and Morongo area of the western Mojave Desert, southern California. Regionally, rock and surficial alluvium in the Mojave River and Morongo areas are high in arsenic, low in chromium and uranium, and near the average bulk continental crust concentration for vanadium. Locally, high chromium concentrations are present in mafic rock within the San Gabriel Mountains; high uranium concentrations are present in felsic rock within the San Bernardino Mountains; and high arsenic, uranium, and vanadium concentrations are present in extrusive (volcanic) felsic rock within uplands surrounding groundwater basins along the Mojave River downstream from Barstow, California. Elemental assemblages identified using principal component analyses (PCA) of pXRF data were used to characterize felsic, mafic, and felsic volcanic source terranes in rock, surficial alluvium, and in geologic material penetrated by selected monitoring wells drilled between 1994 and 2018. Highly felsic alluvium associated with recent deposition from the Mojave River was identified along the 90-mile length of the floodplain aquifer along the river. The thickness of these highly felsic alluvial deposits ranged from 200 feet (ft) near Victorville and near Barstow to a thin veneer about 30 ft thick downstream from Victorville and downstream portions of the floodplain aquifer within the Mojave Valley. Groundwater in the Mojave River and Morongo areas was generally oxic and alkaline (pH≥7.5). Maximum concentrations of arsenic, hexavalent chromium [Cr(VI)], uranium, and vanadium in water from as many as 498 wells sampled between 2000 and 2018 were 360, 140, 1,470, and 690 micrograms per liter (μg/L), respectively. Water from 22 percent of sampled wells exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) for arsenic of 10 μg/L, with arsenic concentrations commonly exceeding the MCL in water from wells east of Barstow, deep wells in the Victorville fan, and in suboxic or reduced groundwater within the floodplain aquifer. Water from about 1 percent of sampled wells had Cr(VI) concentrations greater than the California MCL for total chromium of 50 μg/L, whereas 13 percent of sampled wells had Cr(VI) concentrations greater than the former California MCL of 10 μg/L. Hexavalent chromium concentrations were highest in water from wells in the Sheep Creek alluvial fan, eroded from mafic rock in the San Gabriel Mountains, although Cr(VI) concentrations greater than the former California MCL also were present elsewhere in the study area where mafic materials or older groundwater were present. Water from about 9 percent of sampled wells exceeded the EPA MCL for uranium of 30 μg/L, with concentrations exceeding the MCL commonly associated with irrigation return from agricultural land overlying the floodplain aquifer. Water from about 7 percent of sampled wells had vanadium concentrations greater than the California notification level of 50 μg/L; most of these wells were in the Victorville fan within the Mojave River area. In general, arsenic concentrations were higher in suboxic or reduced water; chromium concentrations were higher in oxic, alkaline (pH≥7.5) water; uranium concentrations were higher in circumneutral to slightly alkaline water (pH≤7.4); and vanadium concentrations were higher in highly alkaline (pH≥8.0) water, independent of redox status. Concentrations within geologic source terranes are not the sole factor controlling the concentrations of geogenic elements in groundwater. Differences in mineral weathering, pH-dependent sorption to surface-exchange sites on mineral grains, and aqueous geochemistry (especially redox status and pH) affect geogenic element concentrations in groundwater. Consequently, the relative abundances of arsenic, Cr(VI), uranium, and vanadium in groundwater differ from their relative abundances in the average bulk continental crust and their regional abundances in rock and surficial alluvium within groundwater basins of the western Mojave Desert. Processes that control the concentrations of arsenic, chromium, uranium, and vanadium in groundwater operate at the mineral-grain and aquifer scale. At the mineral-grain scale, sequential chemical extraction data show arsenic and uranium are more available to groundwater (under specific geochemical conditions) than chromium or vanadium, which largely are unavailable within unweathered mineral grains. Additionally, chromium and vanadium form few aqueous complexes and bind tightly with iron minerals within surface coatings on mineral grains making them less available to groundwater, whereas complexation with other dissolved ions enhances the solubility of uranium and, to a lesser extent, arsenic. Complexation also increases the valence (less negative charge) and increases the size of dissolved oxyanions and uranium complexes with bicarbonate and carbonate making them less readily sorbed to aquifer materials. At the aquifer scale, hydrogeology (including isolation of water in aquifers from surface sources of recharge, older groundwater age, and long contact times between groundwater and aquifer materials) combined with geochemical processes (such as silicate weathering) to produce alkaline groundwater. Desorption from sorption sites on the surfaces of mineral grains with increasing pH increases arsenic, chromium, and vanadium concentrations in water from wells and increases Cr(VI) concentrations as long as water remains oxic. Aqueous geochemistry and concentrations of geogenic contaminants also are affected by anthropogenic activities including (1) discharge of treated municipal wastewater, which may change the redox status of groundwater; (2) return from irrigated agriculture, which may alter the chemistry of groundwater and increase the solubility of trace elements such as uranium; and (3) groundwater pumping and subsequent water-level declines, which may change the source of water yielded by wells. The quality of water imported from northern California and infiltrated from ponds for groundwater recharge may be altered by naturally present trace elements, especially uranium in areas of agricultural land use or chromium within mafic alluvium.

California↗

Surrogate regression models estimating nitrate concentrations at six springs in Gooding County, south-central Idaho, 2018–22

Populations of endangered Banbury Springs limpet ( Idaholanx fresti ) and threatened Bliss Rapids snail ( Taylorconcha serpenticola ) are declining in springs north of the Snake River along the southern Gooding County boundary, in south-central Idaho. One hypothesis for the decline is that increased macrophyte growth, associated with elevated nitrate concentrations in the springs, is decreasing aquatic habitat for the limpet and snail populations. In support of U.S. Fish and Wildlife Service efforts to understand the population declines, the U.S. Geological Survey developed surrogate regression models to estimate nitrate concentrations at six springs influenced by upgradient agriculture, which results in an increase and decrease each year of streamflow, specific conductance, and nitrate concentrations. The surrogate regression models use continuous specific conductance data and streamflow data (available at two springs from existing U.S. Geological Survey streamgages). The spring surrogate regression models showed that specific conductance can be an effective surrogate for nitrate in springs affected by agriculture and that the model results improved when streamflow data were included. Four of the six springs had surrogate regression models (using specific conductance and day of the year as explanatory variables) that performed well based on model summary statistics, and these models improved further with the inclusion of streamflow as an explanatory variable. The surrogate regression models at four springs had coefficient of determination (R 2 ) values ranging from 0.79 to 0.94. The root mean squared error of the four models ranged from 0.07 to 0.11 milligrams per liter. Two of the six springs were not well modeled, with adjusted R 2 values of 0.15 and 0.80. The surrogate regression models for these two springs also did not meet the required assumption of linearity between explanatory and response variables for linear regression. The surrogate regression models show that specific conductance can be an effective surrogate for nitrate in springs affected by agriculture and that models are improved where streamflow data are included. These surrogates improve understanding of nitrate concentration variability in the springs.

Idaho↗

Nitrogen load estimates from six nonpoint sources on Long Island, New York, from 1900 to 2019

Estimates of nitrogen loading from nonpoint sources on Long Island, New York, at or just below the land surface, are essential for assessing the current and future effects of nitrogen on the island’s drinking water and fresh and marine surface receiving waters. Annual estimates of nitrogen loading for the 120 years from 1900 to 2019 for major nonpoint nitrogen sources—septic systems, residential fertilizer, agricultural fertilizer, livestock waste, pet waste, and atmospheric deposition—were made by using a geographic information system to analyze, visualize, and process data sources, and format output data. This analysis provided spatial and temporal estimates of nitrogen loading derived from each nonpoint source at a 500- by 500-foot gridded resolution and represents the total mass of nitrogen applied on, or just below, the land surface annually from 1900 to 2019. These mass estimates are considered unattenuated as they do not reflect the various mechanisms of nitrogen loss, such as plant uptake, overland runoff, and chemical transformations in the soil and unsaturated zones that likely reduce the amount of nitrogen that reaches the water table. Island-wide and countywide summaries of the estimated nitrogen loading were analyzed to describe the long-term average nitrogen totals and the contributions from the six nonpoint sources. The island-wide average annual nitrogen load from 1900 to 2019 was 14.92 million kilograms (Mkg) nitrogen, which represents the aggregate of individual contributions from septic systems (4.15 Mkg), residential fertilizer (3.28 Mkg), agricultural fertilizer (3.19 Mkg), livestock waste (1.22 Mkg), pet waste (0.98 Mkg), and atmospheric deposition (2.10 Mkg). The island-wide average annual nitrogen load, normalized by area, was 4,100 kilograms nitrogen per square kilometer (kg N/km 2 ), which represents the aggregate of individual contributions from septic systems (1,100 kg N/km 2 ), residential fertilizer (910 kg N/km 2 ), agricultural fertilizer (880 kg N/km 2 ), livestock waste (340 kg N/km 2 ), pet waste (270 kg N/km 2 ), and atmospheric deposition (580 kg N/km 2 ).

New York↗

Evaluating drought risk of the Red River of the North Basin using historical and stochastic streamflow upstream from Emerson, Manitoba

Drought and its effect on streamflow are important to understand because of the potential to adversely affect water supply, agricultural production, and ecological conditions. The Red River of the North Basin in north-central United States and central Canada is susceptible to dry conditions. During an extended drought, streamflow conditions in the Red River of the North may become inadequate to support existing water supply needs in the basin for agriculture, industry, human use, and aquatic life. To understand potential future low-streamflow conditions in the Red River of the North Basin, the U.S. Geological Survey, in cooperation with the International Joint Commission, North Dakota Department of Water Resources, Red River Joint Water Resource District, and Red River Watershed Management Board, developed a water-balance model of the Red River of the North Basin upstream from Emerson, Manitoba, Canada, and coupled the model with stochastic weather inputs to simulate possible future low-streamflow conditions. Historical changes in low-streamflow conditions were characterized across the Red River of the North Basin using multiple change-point analysis for 12 streamgages. Across these stations, significant change-point years in 1943 and 1994 marked increases in the magnitude of low-streamflow conditions. During 1920–2015, conversion of primary land (not affected by human use) to agricultural and secondary land was followed by a conversion from smalls grains to corn and soybeans as the dominant crop type. From land-use analysis, 1940–2000 was determined to have relatively stable land use and therefore was used as the calibration period for the water-balance model. A deterministic water-balance model was developed for the Red River of the North Basin upstream from Emerson, Manitoba. The water-balance model was calibrated with data from 37 U.S. Geological Survey streamgages for 1940–2000 and verified using data for 2001–15. The calibrated water-balance model simulated streamflow distributions that mirrored the seasonal patterns of the observed mean monthly streamflow and the standard deviation of the monthly streamflow data, especially during the fall and winter months when streamflow was lowest. For the verification period, during the low-streamflow months of December through January, the difference between simulated and observed data was similar to the calibration comparison and successfully reproduced seasonal trends in the distribution of streamflow, even when using weather data that were outside the calibration period. To determine the future risk of low-streamflow conditions in the Red River of the North Basin, a block-bootstrap method was used to generate multiple possible future climates. These stochastically generated weather time series were then input to a water-balance model to simulate a distribution of possible streamflows. Three sets of experiments were performed, with each experiment containing a set of scenarios. The first set of experiments from the stochastic streamflow model were designed to investigate how changes in reservoir management would affect the distribution of low streamflow. Relative to scenario 1 (present-day [2023] reservoir operation), scenario 2 (no reservoir operation) shifted the low-streamflow frequency curves downward, reducing the annual minimum monthly streamflow for the Emerson subbasin. Subbasins were defined by the contributing area upstream from a selected streamgage station. Relative to scenario 1, scenario 3 (regulated streamflow with an increased reservoir capacity of 10 percent) shifted the low-streamflow frequency curves upward for the Emerson subbasin. The magnitude of this upward shift, caused by increased reservoir capacity, was lower than the magnitude of the shift caused by the absence of the reservoirs, which indicates that the streamflow was most affected when the reservoirs were first constructed. The second set of experiments from the stochastic streamflow model included two scenarios that were performed to better understand how the Red River of the North Basin responds to long periods of low or high precipitation. The results indicate that the model consistently overestimated streamflow, but the relative change between a wet and dry climate state of simulated streamflow distribution reasonably matched the relative change of historical streamflow. Across the subbasins, the model was most accurate for low-streamflow conditions associated with nonexceedance probabilities between 20 and 40 percent. The third set of experiments from the stochastic streamflow model were done to investigate low-streamflow response across the basin to several drought events. Low-end streamflow was reduced when the basin was exposed to a drought, and the magnitude of the reduction increased with longer or more intense droughts. Compared to the low-intensity drought scenarios, the range of percent reductions (as indicated by the interquartile range) was larger for the high-intensity drought scenarios for all subbasins, and the subbasins of Grand Forks and Emerson had a smaller range of reductions compared to the other three subbasins. The larger drainage area—combined with the large contribution of the Red Lake River and several other Minnesota tributaries that generally experience wetter climate conditions—upstream from the Emerson and Grand Forks subbasins may contribute to the smaller range in reductions under the high intensity scenarios. Comparison of the percent reduction in low-end streamflow among subbasins also indicated that the effects of drought duration and intensity could be cumulative. Combining factors of time and intensity produced a larger reduction in streamflow than when each effect was isolated. The array of drought scenarios can be used to determine how a subbasin would respond to multiple possible future conditions. Based on climate predictions, the drought scenario that best matches a future anticipated drought scenario can be used to estimate a low streamflow response for a given subbasin.

Manitoba, Minnesota, North Dakota↗

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California↗

Fluvial sediment dynamics in the Shoshone River and tributaries around Willwood Dam, Park County, Wyoming

Sedimentation affects many of the aging reservoirs in the United States. Dams and water diversions from rivers have been central elements of infrastructure supporting agricultural irrigation in the arid and semiarid regions of the Western United States for more than a century. The Willwood Irrigation District diversion dam (hereafter referred to as “Willwood Dam”) in Park County, Wyoming, is approximately 12 miles northeast of Cody, Wyo.; has a structural height of 70 feet; and impounds the Shoshone River for diversion into the Willwood Canal. Willwood Dam is part of a larger irrigation scheme supported by water storage in the much larger Buffalo Bill Dam, which is approximately 20 miles upstream. In October 2016, renovation construction activities at Willwood Dam and the Willwood Canal caused an unplanned evacuation of nearly 96,000 cubic yards of fine sediment. The fine sediment release in 2016 raised concerns that ongoing sediment management at Willwood Dam could impose limits on the long-term health of the aquatic ecosystem and fish populations. The U.S. Geological Survey, in cooperation with Wyoming Department of Environmental Quality and Willwood Work Groups 2 and 3, initiated an investigation of the dynamics of sediment transport in the Shoshone River and selected tributaries between Buffalo Bill Dam and Willwood Dam. The goal of the study was to quantify sediment transport into and out of Willwood Dam on an annual, seasonal, and event basis to better understand the relative quantities of sediment coming from natural sources and human activities on the landscape. The study ran from March 2019 through October 2021 and used observations of streamflow, turbidity, and acoustic backscatter collected at streamgages upstream and downstream from Willwood Dam to quantify suspended-sediment loads into and out of the dam during irrigation and fallow seasons, precipitation-runoff events, and deliberate sediment releases. Each tributary’s relative contribution to the sediment load upstream from Willwood Dam was examined using discrete measurements of suspended-sediment concentration and bedload during irrigation and fallow seasons, precipitation events, and stable conditions. Analysis of daily precipitation and temperature data indicated that conditions in the study area during the 2019 agricultural year were wetter and colder than period of record normal, and drier and near normal temperatures for the 2020 and 2021 agricultural years. Not all sediment load records between 2019 and 2021 are complete because of rejected observations (outliers), instrument failures or fouling, and instrument removal for calibrations. Statistical modeling of suspended-sediment concentration using paired values of turbidity and acoustic backscatter produced four models that, after refinement, had coefficients of determination indicating that more than 84 percent of the variance was explained by either turbidity or acoustic backscatter. A system of rules was developed to select the model predictions based on the seasonal operations of Willwood Dam, assumptions about the grain sizes mobilized during these operations, and assumed accuracy of the models at the downstream streamgage (Shoshone River below Willwood Dam, near Ralston, Wyo. [streamgage 06284010]) under different operational conditions. The sediment budget between upstream and downstream estimates of loads was interpreted using the mean predicted values bound by their respective model prediction intervals. When mean predicted loads of one streamgage were contained in the prediction intervals of the other streamgage, and vice-versa, difference in the sediment budget were interpreted as “indeterminate.” Modeled sediment load balances demonstrated the depositional and erosional behaviors expected from the conceptual model of dam operations whereby sediment tends to accumulate during irrigation seasons when the dam is spilling over the top, and sediment tends to evacuate during the fallow seasons when it is flowing through the sluice gates at the base of the dam. The sediment load calculations using the rules-based model criteria indicated that between 14,200 and 380,000 tons of suspended sediment moved through the Shoshone River around Willwood Dam during the irrigation seasons of 2019, 2020, and 2021; 380,000 tons of suspended sediment were transported during the cool, wet year of 2019, and 14,200 tons of suspended sediment were transported in 2020, which was relatively dry. During fallow seasons 2019, 2020, and 2021, which had fewer complete records, between 1,140 and 106,000 tons of suspended sediment was estimated to have moved through the river. For all seasons except fallow season 2022, the models estimated that more sediment was released from the dam than entered the dam, but the modeled mean loads at each streamgage were nearly always within the prediction intervals of each other, making the sediment balance indeterminant. Examination of suspended-sediment loads during irrigation seasons indicated that between 65 and 85 percent of fine sediment was transported during annual high flows and storm events, with the remainder transported during steady, lower streamflows. Examination of suspended loads during fallow seasons indicated that deliberate sediment releases through Willwood Dam accounted for between 39 and 67 percent of the total sediment moved during the fallow seasons. Deliberate sediment releases from Willwood Dam had estimated net exports of between 1,360 and 22,400 tons. Between August 2017 and July 2023, suspended-sediment concentration and bedload sediment samples were collected from 9 tributaries to the Shoshone River during 137 sampling events, including stable and precipitation-runoff conditions. During irrigation season precipitation events, the mean total sediment yields ranged from 0.33 to 9.51 tons per day per square mile; during fallow season precipitation events, the mean total yields ranged from 0.04 to 0.95 ton per day per square mile. The mean total sediment yield per unit area across all samples at each tributary site ranged from 0.26 to 3.08 tons per day per square mile. Bedload was a minor fraction of the total load, constituting a mean of 4 percent across all samples; 3 and 6 percent for events and nonevents, respectively, during irrigation season; and 3 and 1 percent for events and nonevents, respectively, during the fallow season. With the exception of one tributary, Dry Creek, these mean yield values were within the range of watershed-scale background sediment yield values estimated from reservoir surveys and previous suspended-sediment studies. Imagery from irrigation seasons 2012, 2015, 2017, 2019, and 2022 was used to determine the planimetric backwater extent of the pool area in the Shoshone River behind Willwood Dam to identify any changes in sediment storage. Active river channel widths in the Shoshone River upstream from Willwood Dam were all similar between years except 2015, which was determined to be statistically different from all other years. Bathymetric data taken in the pool behind Willwood Dam during three different surveys between November 2017 and April 2022 indicated no statistically significant differences in bed elevations between the years. Results from the planimetric and bathymetric survey data provide multiple lines of evidence indicating that sediment did not accumulate behind the dam within the error of the methods used. Examination of how precipitation affects sediment transport in the Shoshone River upstream from Willwood Dam indicated that accumulated rainfall from the natural runoff events captured during the study period varied from a trace to as much as 4.26 inches, with associated predicted suspended-sediment loads varying from 112 to 232,000 tons of suspended sediment. The behavior of the sediment loads relative to accumulated precipitation did not appear to change depending on irrigation or fallow season. A model of suspended-sediment concentrations relative to the 2-day accumulated precipitation indicated that suspended-sediment concentrations in the Shoshone River upstream from Willwood Dam increased exponentially for accumulations of 0.3 inch or more; such storms accounted for 10 percent or less of precipitation events observed during the 1981 to 2018 period of record. The gaps in records, precision of the instrumentation, and large variation in grain sizes in suspended-sediment mixtures downstream from the dam made closing the sediment budgets for most seasons unattainable. The biggest recent change in sediment storage measured using the planimetric area of deposits behind Willwood Dam took place between 2015 and 2017. The main event between these two measurements was the installation of new Willwood Canal gates in October 2016, which resulted in the large unplanned sediment release. Because the sediment budgets were nearly always indeterminate and the planimetric and bathymetric data indicated little change in the bed and bank material, it is likely that the change in sediment storage behind the dam during the study period was small relative to the precision of the statistical models and other uncertainties. This body of evidence suggests that, averaged during the 3-year study period, no major changes in storage took place, and that the current operations may be keeping storage at near-equilibrium. This condition could have been initiated because the middle sluice gate has now been operational since 2014, and the sediment release in October 2016 evacuated a large amount of legacy sediment from storage. Although the uncertainties are large, sluicing events allow for controlled releases of sediment that contributed to the near equilibrium conditions observed over an annual basis during this study.

Wyoming↗

Multidecadal change in pesticide concentrations relative to human health benchmarks in the Nation’s groundwater

Groundwater-quality trend assessments identify aquifers that are responding to changes in pesticide use and the compounds that may pose a threat to water availability. The U.S. Geological Survey has been monitoring pesticide concentrations in groundwater for 25 principal aquifers across the conterminous United States since 1993. The groundwater well locations represent a range of soils, climate, and landforms. The wells are used to monitor groundwater underlying selected agricultural and urban settings and groundwater used for domestic supply. This study examined changes in relative concentrations, defined here as the percentage of wells with pesticide concentrations exceeding a human health benchmark (HHB). HHBs used in this report are legally enforceable drinking-water standards and nonenforceable drinking water levels. Relative pesticide concentration increases may lead to decreased water availability, as restrictions may be put in place for groundwater used as a drinking-water source. This study focused on concentration changes in 22 pesticides that were included in laboratory analysis from 1993 to 2023. The analysis and interpretation of these pesticide concentrations in groundwater have been separated into approximate decadal intervals (decade 1 (1993–2001), decade 2 (2002–12), and decade 3 (2013–22). For one pesticide, 1,2-dibromo-3-chloropropane (DBCP), concentration data were also collected in decade 4 (2023–onward). Atrazine, deethylatrazine, alachlor, prometon, and simazine were 5 pesticides detected at moderate concentrations (greater than 10 percent of the HHB but less than or equal to the HHB). The percentage of wells that had groundwater pesticide concentrations in the moderate concentration category decreased from 7 percent in decade 1 to 2 percent in decade 3. The agricultural networks had the highest percentages of wells with moderate concentrations, and these percentages decreased from 13 percent in decade 1 to 4 percent in decade 3. Moderate concentrations in the urban networks decreased between decades 1 and 2 from 4 percent to 0 percent. No moderate concentrations occurred in the urban networks in decade 3. The percentage of wells with moderate concentrations in the domestic supply networks (1 percent) was the lowest of all the network types and did not change across the three decades. Moderate atrazine or deethylatrazine concentrations occurred across all three decades in aggregated ecoregions representing similar soils, climate, and landforms in the Semiarid West, Midcontinent, and Northeastern United States. Moderate concentrations of prometon, alachlor, and simazine also occurred in the Midcontinent, Arid West, Northeast, South Atlantic Gulf, and Semiarid West regions, but the moderate concentrations did not persist across all three decades. DBCP was the only pesticide that exceeded its respective HHB, and the exceedances occurred across all four decades. In this report, the DBCP analysis was limited to one well network in the Central Valley, California. Agricultural use of DBCP was suspended in 1977. Forty-five years after being banned, DBCP concentrations were greater than the maximum contaminant level of 2 micrograms per liter (μg/L), but the number of exceedances decreased from 50 percent to 15 percent of the samples between 1993 and 2023. This assessment of decadal groundwater pesticide concentrations provides a characterization of changes in water availability because of pesticide contamination in areas where groundwater is used as a drinking-water source. The results highlight the importance of continued long-term monitoring and assessment of groundwater pesticides to identify locations and specific compounds that may pose a potential risk to human health.

continental United States↗

Decadal trends in the quality of groundwater used for public drinking-water supply in California, 2004–23, California Groundwater Ambient Monitoring and Assessment Program, Priority Basin Project

This study provides a comprehensive assessment of decadal changes in the quality of groundwater used for public drinking-water supply at 444 monitoring sites across California during 2004–23. We assessed decadal step trends in groundwater quality for 145 water-quality constituents and geochemical indicators statewide and across geographic and land-use based network groups. We evaluated the statistical significance of directional changes (predominant increase or decrease of constituent concentrations) and the magnitude of those changes across all network groups. Uranium showed the most widespread directional and high-magnitude increases of all constituents with regulatory benchmarks statewide, particularly in the agriculture-dominated Central Valley as well as urban-and desert-dominated regions of Southern California. Fluoride and perchlorate showed the most widespread directional and high-magnitude decreases of all constituents with regulatory benchmarks statewide, which were also most pronounced in Southern California. Although arsenic and nitrate did not often register significant directional changes across network groups, they showed widespread, high-magnitude changes in both directions (increase and decrease) at levels often exceeding 10 percent of respective regulatory benchmarks statewide. Triazine herbicides (atrazine and simazine) and the gasoline oxygenate methyl tert -butyl ether showed significant directional decreases statewide, but not at levels considered to be of high magnitude compared to respective regulatory benchmarks. We observed significant directional and high-magnitude increases of total dissolved solids statewide, which were most pronounced in agricultural areas. Analysis of explanatory geochemical indicators indicated that prevalent statewide increases of alkalinity and calcium were the predominant components of the observed statewide increases in total dissolved solids by mass. Widespread increases in groundwater alkalinity and calcium across agricultural and urban areas may be related, in part, to warm-season irrigation and other anthropogenic factors that have shifted soil weathering dynamics over the long term. Increasing alkalinity concentrations were related to increasing uranium concentrations, particularly in areas with aquifer materials derived from granitic rocks. Conversely, increasing calcium concentrations were related to decreasing fluoride concentrations, particularly in areas where fluoride occurred naturally at elevated concentrations. Decrease of perchlorate, triazine herbicides, and methyl tert -butyl ether are likely related to decreased anthropogenic source inputs over time and natural attenuation in aquifers.

California↗