SEARCH · Search USGS

Results for “Agriculture”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,387 records · Page 77Linked to original sources

USGS2015

Areas contributing recharge to production wells and effects of climate change on the groundwater system in the Chipuxet River and Chickasheen Brook Basins, Rhode Island

The Chipuxet River and Chickasheen Brook Basins in southern Rhode Island are an important water resource for public and domestic supply, irrigation, recreation, and aquatic habitat. The U.S. Geological Survey, in cooperation with the Rhode Island Department of Health, began a study in 2012 as part of an effort to protect the source of water to six large-capacity production wells that supply drinking water and to increase understanding of how climate change might affect the water resources in the basins. Soil-water-balance and groundwater-flow models were developed to delineate the areas contributing recharge to the wells and to quantify the hydrologic response to climate change. Surficial deposits of glacial origin ranging from a few feet to more than 200 feet thick overlie bedrock in the 24.4-square mile study area. These deposits comprise a complex and productive aquifer system. Simulated areas contributing recharge to the production wells covered a total area of 0.63 square miles for average well withdrawal rates from 2007 through 2011 (total rate of 583 gallons per minute). Simulated areas contributing recharge for the maximum well pumping capacities (total rate of 3,700 gallons per minute) covered a total area of 2.55 square miles. Most simulated areas contributing recharge extend upgradient of the wells to morainal and upland till deposits and to groundwater divides. Some simulated areas contributing recharge include small, isolated areas remote from the wells. Relatively short groundwater traveltimes from recharging locations to discharging wells indicated that the wells are vulnerable to contamination from land-surface activities; median traveltimes ranged from 3.5 to 8.6 years for the production wells examined, and 57 to 91 percent of the traveltimes were 10 years or less. Land cover in the areas contributing recharge includes a substantial amount of urban and agriculture land use for five wells adjacent to the Chipuxet River; for one well adjacent to a tributary stream, land use is less developed. The calibrated groundwater-flow model provided a single, best representation of the areas contributing recharge to a production well. 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 predicted contributing areas to the wells. A Monte Carlo analysis led to contributing areas expressed as a probability distribution that differed from a single deterministic contributing area. Because of the effects of parameter uncertainty, the size of the probabilistic contributing areas for both average and maximum pumping rates was larger than the size of the deterministic contributing areas for the wells. Thus, some areas not in the deterministic contributing area might actually be in the contributing area, including additional areas of urban and agricultural land use that has 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. At the maximum pumping rates, areas associated with low probabilities extended long distances along groundwater divides in the uplands remote from the wells. Climate projections for the Chipuxet River and Chickasheen Brook Basins from downscaled output from general circulation models indicate that mean annual temperature might increase by 4.7 degrees Fahrenheit and 8.0 degrees Fahrenheit by the late 21st century (2070–99) compared with the late 20th century (1970–99) under scenarios of lower and higher emissions of greenhouse gases, respectively. By the late 21st century, winter and spring precipitation is projected to increase by 12 to 17 percent, summer precipitation to increase by about the same as mean annual precipitation (8 percent), and fall precipitation to decrease by 5 percent for both emission scenarios compared with the late 20th century. Soil-water-balance simulations indicate that, although precipitation is expected to increase in three seasons, only in winter do precipitation increases exceed actual evapotranspiration increases. Recharge is projected to decrease in fall and generally change little in spring and summer. By the late 21st century, winter recharge is expected to increase by 13 percent for the lower emissions scenario and by 15 percent for the higher emissions scenario. In fall, recharge is projected to diminish by 13 percent for the lower emissions scenario and by 24 percent for the higher emissions scenario. Although recharge is projected to change seasonally in the 21st century, mean annual recharge changes minimally. Soil moisture is projected to decrease in the 21st century from spring through fall because of increases in potential evapotranspiration, and in fall because of decreases in precipitation in addition to increases in potential evapotranspiration. By the late 21st century, soil moisture for the lower emissions scenario is expected to decrease by 11 percent in summer and 15 percent in fall, and for the higher emissions scenario, decrease by 23 percent for both seasons. These decreases in soil moisture during the growing season might have implications for agriculture in the study area. Predicted changes in the magnitude and seasonal distribution of recharge in the 21st century increase simulated base flows and groundwater levels in the winter months for both emission scenarios, but because of less recharge in the fall and less or about the same recharge in the preceding months of spring and summer, base flows and groundwater levels in the fall months decrease for both emission scenarios. October has the largest base flow and groundwater level decreases. By the late 21st century, base flows at the Chipuxet River in October are projected to decrease by 9 percent for the lower emissions scenario and 18 percent for the higher emissions scenario. For a headwater stream in the upland till with shorter groundwater-flow paths and lower storage properties in its drainage area, base flows in October are projected to diminish by 28 percent and 42 percent for the lower and higher emissions scenarios by the late 21st century. Groundwater level changes in the uplands show substantial decreases in fall, but because of the large storage capacity of stratified deposits, water levels change minimally in the valley. By the late 21st century, water levels in large areas of upland till deposits in October are projected to decrease by up to 2 feet for the lower emissions scenario, whereas large areas decrease by up to 5 feet, with small areas with decreases of as much as 10 feet, for the higher emissions scenario. For both emission scenarios, additional areas of till go dry in fall compared with the late 20th century. Thus projected changes in recharge in the 21st century might extend low flows and low water levels for the year later in fall and there might be more intermittent headwater streams compared with the late 20th century with corresponding implications to aquatic habitat. Finally, the size and location of the simulated areas contributing recharge to the production wells are minimally affected by climate change because mean annual recharge, which is used to determine the contributing areas to the production wells, is projected to change little in the 21st century.

Rhode Island↗

USGS2015

Characterization of streamflow, salinity, and selenium loading and land-use change in Montrose Arroyo, western Colorado, from 1992 to 2013

Salinity and dissolved selenium are known water-quality impairments in the lower Gunnison River watershed of western Colorado. Salinity is a concern because of its adverse effects on agricultural land and equipment, and on municipal and industrial users. The Montrose Arroyo watershed in Montrose, Colorado, contains agricultural and residential areas as well as undeveloped land and has undergone substantial land-use change since the early 1990s. Previous sampling efforts indicated salinity concentrations and loads have remained constant since land-use change began in the early 1990s; however, recent sampling also indicated that dissolved-selenium concentrations and loads have begun to increase. In response to the potential increasing dissolved-selenium concentrations and loads, the U.S. Geological Survey—in cooperation with the Bureau of Reclamation; Colorado River Basin Salinity Control Forum; and Colorado River Water Conservation District—continued to monitor salinity and dissolved-selenium concentrations and loads in the Montrose Arroyo watershed. This report characterizes salinity and dissolved-selenium loads in Montrose Arroyo from 1992 to 2013 at three sites: Montrose Arroyo at East Niagara Street (MA2, U.S. Geological Survey site identification number 382802107513301), Montrose Arroyo at 6700 Road (MA3, U.S. Geological Survey site identification number 382711107500501), and Montrose Arroyo at 6750 and Ogden Roads (MA4, U.S. Geological Survey site identification number 382702107493701). A detailed land-use change analysis was also characterized in the MA3 subwatershed. The three sites were used to monitor salinity and dissolved-selenium concentrations and loads in Montrose Arroyo. Over a period from December 2011 through September 2013, 46 water-quality samples and streamflow measurements were collected at sites MA2, MA3, and MA4. Streamflow, salinity concentrations and loads, and selenium concentrations and loads were characterized and compared between the pre-lateral (before April 1999) and post-growth periods (March 2008 through September 2013) and between post-lateral (April 1999 through October 2000) and post-growth periods. Results from a previous USGS study on the characterization of salinity and selenium loading and land-use change in Montrose Arroyo from 1992 to 2010 indicated that there was no change in salinity load at site MA2 and a significant increase in dissolved-selenium load from the pre-lateral period to the post-growth period. Data associated with this report indicate the selenium loads at site MA2 show no significant change from the pre-lateral period to the post-growth period. In addition, both salinity load and dissolved-selenium load at site MA3 show significant decreases in salinity and dissolved-selenium load for both periods that are potentially associated with land-use change. Land use was characterized for 1992, 2002, and 2009 for site MA3. The common land-use change in the MA3 subwatershed was a conversion from previously irrigated agricultural land to urban land use. The MA3 subwatershed had 124 acres of irrigated land use converted to urban land use and 27.1 acres of unirrigated desert converted to urban land use from 1992 to 2009. Consistent with findings in previous land-use change reports, salinity and dissolved-selenium loading at site MA3 showed significant decreases as irrigated land was converted to urban land use.

Colorado↗

USGS2015

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↗

USGS2016

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↗

USGS2017

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↗

USGS2020

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↗

USGS2020

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↗

USGS2020

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↗

USGS2022

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↗

USGS2022

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↗

USGS2022

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↗

USGS2022

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↗

USGS2023

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↗

USGS2023

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↗

USGS2024

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↗

USGS2025

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↗

USGS2025

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↗

USGS2025

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↗