USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

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

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

At least 145 records · Page 8Linked to original sources
USGS2022

Hydrologic budget of the Harney Basin groundwater system, southeastern Oregon

Groundwater-level declines and limited quantitative knowledge of the groundwater-flow system in the Harney Basin prompted a cooperative study between the U.S. Geological Survey and the Oregon Water Resources Department to evaluate the groundwater-flow system and budget. This report provides a hydrologic budget of the Harney Basin groundwater system that includes separate groundwater budgets for upland and lowland areas to avoid double counting water that recharges in the uplands, discharges to streams and springs in the uplands, flows downstream to the lowlands, and recharges the lowland groundwater system. Lowlands generally represent the conterminous valleys within the center of the basin, including floodplains of the major streams and uplands represent all other areas in the basin. The upland groundwater budget is minimally affected by groundwater development and generally represents the budget of the natural system. In upland areas during 1982–2016, mean-annual recharge totaled 288,000 acre-feet (acre-ft) and mean-annual discharge totaled 239,000 acre-ft, resulting in a net recharge of 49,000 acre-ft. Upland groundwater recharge occurs as infiltration of precipitation and snowmelt and was estimated using the USGS Soil-Water-Balance model calibrated to estimates of runoff, evapotranspiration (ET), base flow, and snow-water equivalent. Groundwater discharge to streams is the predominant discharge mechanism in upland areas and was estimated as 225,000 acre-feet per year (acre-ft/yr) during 1982–2016 using hydrograph separation and summer low-flow estimates in streamgaged watersheds and a linear relation between estimated streamflow and base flow in ungaged watersheds. The remaining upland discharge occurs through springs (14,000 acre-ft/yr) that either emerge downgradient of locations where groundwater discharge to streams was estimated or are routed to irrigated areas. Spring discharge was estimated as a compilation of current and historical measurements. The net upland recharge, which is 17 percent of total upland recharge, ultimately recharges lowland areas as groundwater flow from uplands to lowlands. The lowland groundwater budget for the Harney Basin represents a combination of natural conditions and human activity as more than 99 percent of groundwater development has occurred either inside or within 2 miles of the lowland boundary. In lowland areas during 1982–2016, mean annual groundwater recharge totaled 173,000 acre-ft and groundwater discharge totaled 283,000 acre-ft, indicating discharge exceeded recharge by more than 60 percent. Excluding groundwater pumping, the lowland groundwater budget is more in balance with a mean annual recharge of 165,000 acre-ft and a mean annual discharge of 131,000 acre-ft during 1982–2016. The 23-percent difference between non-pumping recharge and discharge mostly represents the cumulative uncertainty in the estimates of the various groundwater budget components but also likely includes a small reduction in natural groundwater discharge captured by pumping. Lowland groundwater is predominantly recharged by infiltration of surface water (116,000 acre-ft/yr) through streams, floodwater, and irrigation, with a lesser amount as groundwater inflow from uplands and minimal recharge beneath Malheur and Harney Lakes. Recharge from streams and floodwater (natural and irrigation) was estimated using a balance of measured and estimated surface-water inflow to and outflow from lowland areas including streamflow, springflow, and ET where a portion of surface-water inflow to lowland areas is comprised of upland discharge to streams and springs. Groundwater ET (119,000 acre-ft/yr) is the predominant natural discharge mechanism in lowland areas and was estimated as the mean from two remote-sensing based approaches incorporating groundwater ET measurements from other similar basins and 23 years (1987–2015) of Landsat imagery. Discharge of lowland groundwater into Malheur and Harney Lakes is about 700 acre-ft/yr and is represented in groundwater ET estimates. The remaining natural groundwater discharge from lowland areas issues from Sodhouse Spring (8,900 acre-ft/yr) and as groundwater flow to the Malheur River Basin through Virginia Valley (3,100 acre-ft/yr). The relatively large amount of groundwater discharged to springs in Warm Springs Valley (25,000 acre-ft/yr) is accounted for in groundwater ET estimates. Natural groundwater discharge in lowland areas of the Harney Basin has remained relatively constant during the last 80 years based on comparisons with estimates north of Malheur Lake and west of Harney Lake published in the 1930s. Annual net amount of groundwater pumped (pumpage) from the Harney Basin during 2017–18 averaged 144,000 acre-ft. The net value is the difference between pumpage (about 152,000 acre-ft/yr) and reinfiltration of groundwater pumped for irrigation and non-irrigation purposes (about 8,000 acre-ft/yr). Net pumpage was estimated in concurrent studies that compiled groundwater-use data and coupled reported groundwater pumpage data from wells with remote-sensing-based ET estimates from groundwater-irrigated fields. Total pumpage for irrigation has increased from about 54,000 acre-ft/yr during 1991–92 to 145,000 acre-ft/yr during 2017–18. Presently, pumpage is greatest in the lowland region north of Malheur Lake (81,000 acre-ft/yr), with lesser amounts to the north and northwest of Harney Lake (41,000 acre-ft/yr) and to the south and east of Malheur Lake (22,000 acre-ft/yr). During this study, mean annual lowland groundwater discharge (including pumpage) exceeded mean annual recharge, indicating that the lowland hydrologic budget is out of balance. Net groundwater pumpage during 2017–18 is similar to groundwater discharge from all other sources in the lowlands and is four times the imbalance between non-pumping lowland recharge and discharge (34,000 acre-ft/yr). Declining groundwater levels at depth across many parts of the Harney Basin lowlands indicate that pumpage is depleting aquifer storage and is likely capturing a small amount of natural groundwater discharge to springs and ET in some lowland areas. If pumping continues, aquifer storage depletion will continue until the capture rate of natural discharge to springs and ET is equal to the pumping rate. If groundwater development occurs in upland areas and reduces either the streamflow or groundwater inflow to lowland areas, the deficit in the lowland water budget will increase.

Oregon↗
USGS2012

Groundwater simulation and management models for the upper Klamath Basin, Oregon and California

The upper Klamath Basin encompasses about 8,000 square miles, extending from the Cascade Range east to the Basin and Range geologic province in south-central Oregon and northern California. The geography of the basin is dominated by forested volcanic uplands separated by broad interior basins. Most of the interior basins once held broad shallow lakes and extensive wetlands, but most of these areas have been drained or otherwise modified and are now cultivated. Major parts of the interior basins are managed as wildlife refuges, primarily for migratory waterfowl. The permeable volcanic bedrock of the upper Klamath Basin hosts a substantial regional groundwater system that provides much of the flow to major streams and lakes that, in turn, provide water for wildlife habitat and are the principal source of irrigation water for the basin's agricultural economy. Increased allocation of surface water for endangered species in the past decade has resulted in increased groundwater pumping and growing interest in the use of groundwater for irrigation. The potential effects of increased groundwater pumping on groundwater levels and discharge to springs and streams has caused concern among groundwater users, wildlife and Tribal interests, and State and Federal resource managers. To provide information on the potential impacts of increased groundwater development and to aid in the development of a groundwater management strategy, the U.S. Geological Survey, in collaboration with the Oregon Water Resources Department and the Bureau of Reclamation, has developed a groundwater model that can simulate the response of the hydrologic system to these new stresses. The groundwater model was developed using the U.S. Geological Survey MODFLOW finite-difference modeling code and calibrated using inverse methods to transient conditions from 1989 through 2004 with quarterly stress periods. Groundwater recharge and agricultural and municipal pumping are specified for each stress period. All major streams and most major tributaries for which a substantial part of the flow comes from groundwater discharge are included in the model. Groundwater discharge to agricultural drains, evapotranspiration from aquifers in areas of shallow groundwater, and groundwater flow to and from adjacent basins also are simulated in key areas. The model has the capability to calculate the effects of pumping and other external stresses on groundwater levels, discharge to streams, and other boundary fluxes, such as discharge to drains. Historical data indicate that the groundwater system in the upper Klamath Basin fluctuates in response to decadal climate cycles, with groundwater levels and spring flows rising and declining in response to wet and dry periods. Data also show that groundwater levels fluctuate seasonally and interannually in response to groundwater pumping. The most prominent response is to the marked increase in groundwater pumping starting in 2001. The calibrated model is able to simulate observed decadal-scale climate-driven fluctuations in the groundwater system as well as observed shorter-term pumping-related fluctuations. Example model simulations show that the timing and location of the effects of groundwater pumping vary markedly depending on the pumping location. Pumping from wells close (within a few miles) to groundwater discharge features, such as springs, drains, and certain streams, can affect those features within weeks or months of the onset of pumping, and the impacts can be essentially fully manifested in several years. Simulations indicate that seasonal variations in pumping rates are buffered by the groundwater system, and peak impacts are closer to mean annual pumping rates than to instantaneous rates. Thus, pumping effects are, to a large degree, spread out over the entire year. When pumping locations are distant (more than several miles) from discharge features, the effects take many years or decades to fully impact those features, and much of the pumped water comes from groundwater storage over a broad geographic area even after two decades. Moreover, because the effects are spread out over a broad area, the impacts to individual features are much smaller than in the case of nearby pumping. Simulations show that the discharge features most affected by pumping in the area of the Bureau of Reclamation's Klamath Irrigation Project are agricultural drains, and impacts to other surface-water features are small in comparison. A groundwater management model was developed that uses techniques of constrained optimization along with the groundwater flow model to identify the optimal strategy to meet water user needs while not violating defined constraints on impacts to groundwater levels and streamflows. The coupled groundwater simulation-optimization models were formulated to help identify strategies to meet water demand in the upper Klamath Basin. The models maximize groundwater pumping while simultaneously keeping the detrimental impacts of pumping on groundwater levels and groundwater discharge within prescribed limits. Total groundwater withdrawals were calculated under alternative constraints for drawdown, reductions in groundwater discharge to surface water, and water demand to understand the potential benefits and limitations for groundwater development in the upper Klamath Basin. The simulation-optimization model for the upper Klamath Basin provides an improved understanding of how the groundwater and surface-water system responds to sustained groundwater pumping within the Bureau of Reclamation's Klamath Project. Optimization model results demonstrate that a certain amount of supplemental groundwater pumping can occur without exceeding defined limits on drawdown and stream capture. The results of the different applications of the model demonstrate the importance of identifying constraint limits in order to better define the amount and distribution of groundwater withdrawal that is sustainable.

Oregon;California↗
USGS2016

Assessment of groundwater availability in the Northern Atlantic Coastal Plain aquifer system From Long Island, New York, to North Carolina

Executive Summary The U.S. Geological Survey began a multiyear regional assessment of groundwater availability in the Northern Atlantic Coastal Plain (NACP) aquifer system in 2010 as part of its ongoing regional assessments of groundwater availability of the principal aquifers of the Nation. The goals of this national assessment are to document effects of human activities on water levels and groundwater storage, explore climate variability effects on the regional water budget, and provide consistent and integrated information that is useful to those who use and manage the groundwater resource. As part of this nationwide assessment, the USGS evaluated available groundwater resources within the NACP aquifer system from Long Island, New York, to northeastern North Carolina. The northern Atlantic Coastal Plain physiographic province depends heavily on groundwater to meet agricultural, industrial, and municipal needs. The groundwater assessment of the NACP aquifer system included an evaluation of how water use has changed over time; this evaluation primarily used groundwater budgets and development of a numerical modeling tool to assess system responses to stresses from future human uses and climate trends. This assessment focused on multiple spatial and temporal scales to examine changes in groundwater pumping, storage, and water levels. The regional scale provides a broad view of the sources and demands on the system with time. The sub-regional scale provides an evaluation of the differing response of the aquifer system across geographic areas allowing for closer examination of the interaction between different aquifers and confining units and the changes in these interactions under pumping and recharge conditions in 2013 and hydrologic stresses as much as 45 years in the future. By focusing on multiple scales, water-resource managers may utilize this study to understand system response to changes as they affect the system as a whole. The NACP aquifer system extends from Long Island to northeastern North Carolina, and includes aquifers primarily within New York, New Jersey, Delaware, Maryland, Virginia, and North Carolina. The seaward-dipping sedimentary wedge that underlies the northern Atlantic Coastal Plain physiographic province forms a complex groundwater system. Although the NACP aquifer system is recognized by the U.S. Geological Survey as one of the smallest of the 66 principal aquifer systems in the Nation, it ranks 13th overall in terms of total groundwater withdrawals and is 7th in population served. Despite abundant precipitation [about 45 inches per year (in/yr)], the supply of fresh surface water in this region is limited because many of the surface waters in this area are brackish estuaries, contributing to why many communities in the northern Atlantic Coastal Plain physiographic province rely heavily on groundwater to meet their water needs. Increases in population and changes in land use during the past 100 years have resulted in diverse increased demands for freshwater throughout the northern Atlantic Coastal Plain physiographic province with groundwater serving as a vital source of drinking water for the nearly 20 million people who live in the region. Total groundwater withdrawal in 2013 was estimated to be about 1,300 million gallons per day (Mgal/d) and accounts for about 40 percent of the drinking water supply with the densely populated areas tending to have the highest rates of withdrawals and, therefore, being most susceptible to effects from these withdrawals over time. Water levels in many of the confined aquifers are decreasing by as much as 2 feet per year (ft/yr) in response to extensive development and subsequent increased withdrawals throughout the region. Total water-level decreases (drawdowns) are more than 100 feet (ft) in some aquifers from their predevelopment (before 1900) levels. These drawdowns extend across state lines and under the Chesapeake and Delaware Bays, creating the potential for interstate aquifer management issues. Regional water-resources managers in the northern Atlantic Coastal Plain physiographic province face challenges beyond competing local domestic, industrial, agricultural, and environmental demands for water. Large changes in regional water use have made the State-level management of aquifer resources more difficult because of hydrologic effects that extend beyond State boundaries. The northern Atlantic Coastal Plain physiographic province is underlain by a wedge of unconsolidated to partially consolidated sediments that are typically thousands of feet thick along the coastline with a maximum thickness of about 10,000 ft near the edge of the continental shelf. The NACP aquifer system consists of nine confined aquifers and nine confining units capped by an unconfined surficial aquifer that is bounded laterally from the west by the contact between Coastal Plain sediments and the upland Piedmont bedrock. This aquifer system extends to the east to the limit of the Continental Shelf, however, the boundary between fresh and saline groundwater is considered to be much closer to the shoreline and varies vertically by aquifer. Precipitation over the region for average conditions from 2005 to 2009 is about 61,800 Mgal/d, but about 70 percent of it is lost to evapotranspiration resulting in an inflow of about 19,600 Mgal/d entering the groundwater system as aquifer recharge. Most of this recharge enters the aquifer system and flows through the shallow unconfined aquifer and either discharges to streams or directly to coastal waters without reaching the deep, confined aquifer system. In addition to recharge from precipitation, other sources of water include the return of wastewater from domestic septic systems of about 240 Mgal/d, about 60 Mgal/d of water released from storage in the confined system, and about 30 Mgal/d of lateral inflow at the boundary between freshwater and saltwater in response to pumping for conditions in 2013. The outflow needed to balance the inflows was subdivided between streamflow, discharge to tidal portions of streams, and coastal discharge. The hydrologic budget developed for current [2013] conditions determined that 93 percent of the total outflow was to surface waters with about 70 percent divided evenly between streamflow and shallow coastal discharge and 23 percent as discharge to tidal waters. The remaining 7 percent of the total outflow components include withdrawals from both the surficial and confined aquifers of the groundwater system. The groundwater availability assessment of the NACP aquifer system highlights the importance of analyses at both the regional and local scales to understand how changes in land use, water use, and climate have affected groundwater resources and how these resources may change in the future. The investigation included assessments of the regional changes in water levels and budgets across State lines, the importance of considering storage change in the confining units, the response of the aquifer system to a continuation of current [2013] hydrologic stresses into the future, and the potential effects of climate change and sea-level rise on the aquifer system. The Potomac aquifer group includes two of the most widely used aquifers in the NACP aquifer system, the Potomac-Patapsco and Potomac-Patuxent regional aquifers, providing about 24 percent of the total groundwater used in the region. Withdrawals from large pumping centers in this deep, confined aquifer group have resulted in substantial decreases in water-levels across state lines, particularly between southern Virginia and northeastern North Carolina as well as between southern New Jersey and northern Delaware where water levels in the Potomac-Patapsco aquifer have decreased by as much as 200 ft and 50 ft, respectively from predevelopment to current [2013] conditions. This response in water levels also is reflected in changes in water budgets where, for example, about 20 percent of the total response to pumping in Virginia is met by inducing flow from adjacent States. Understanding and quantifying these hydrologic effects that extend beyond State boundaries is critical for the State- and regional-level management of aquifer resources. The cumulative storage loss from the intervening confining units throughout the entire NACP aquifer system was about 35 percent of the total storage loss from predevelopment to current [2013] conditions. In geographic areas such as Delmarva Peninsula, Maryland, and New Jersey, the water released from storage in the confining units makes up the majority of the total storage release from the groundwater system and is becoming proportionally more important over time as the surficial aquifer approaches equilibrium with respect to pumping and recharge stresses as of 2013. Storage loss from the confining units is of particular concern because, unlike in the sands that comprise the confined aquifers, water removed from the clayey confining unit sediments cannot be replenished as these units gradually compress. This non-recoverable storage loss, if great enough, can result in land subsidence where these units are thick and the release from storage is relatively large and contributes to increased concerns for sea-level rise in areas such as the lower portion of the Chesapeake Bay. Groundwater usage increased dramatically in the NACP aquifer system during post-World War II era from the mid-1940s to early the 1980s, with withdrawals increasing from about 400 Mgal/d to more than 1,300 Mgal/d. Although groundwater withdrawals have been relatively constant since the early 1980s, about half of the total groundwater withdrawn from the NACP aquifer system since 1900 was withdrawn in the past 30 years. An analysis of the response of the groundwater system to a continuation of the current [2013] pumping for an additional 30 years into the future shows that the flow system continues to adjust in terms of changes in water budget components, water levels, and the boundary between freshwater and saltwater as it approaches equilibrium. The largest change in water budget components is the reduction in the amount of water released from storage. Across the entire NACP aquifer system, the reduction of storage release from 7 to 4 percent of the total water budget change is accounted for by reductions in groundwater discharge to streams and coastal waters. Locally, a similar response is calculated for each of the geographic areas except for Virginia where the amount of water released from storage accounts for about 25 percent of the total change in water budget. This finding suggests that the groundwater flow system in Virginia is not approaching equilibrium under the current [2013] stresses and, therefore, water levels will continue to decrease even if the pumping remains constant. An analysis of the change in water levels in the Potomac-Patapsco aquifer as pumping is continued 30 years into the future reveals that the largest decreases in water levels throughout the NACP aquifer system will occur in the southern Virginia and northeastern North Carolina parts of the study area. It is these areas that also see the greatest potential for increased lateral movement of saline groundwater in the deep, confined portion of the groundwater flow system in response to a continuation of the current [2013] pumping rates. The potential effects of long-term climate change and variability on the hydrologic system and availability of water resources in the NACP aquifer system continue to be of serious societal concern. These concerns include the effects of changes in aquifer recharge and in sea-level rise on the groundwater flow system. An assessment of the potential effects of a prolonged drought during current [2013] pumping conditions indicated that the reductions in recharge associated with droughts, including additional irrigation withdrawals required to meet increased crop water demand, have the greatest effects on water levels and streamflows in the surficial aquifer, and changes in water levels in the confined aquifers primarily resulted from the increased withdrawals associated with increased irrigation pumping; this response was most apparent in the Delmarva Peninsula. These results suggest that water levels may not be susceptible to the effects of droughts in the confined aquifers of the NACP aquifer system not used for irrigation, unlike in the unconfined surficial aquifer. A second analysis also was conducted to assess the effects of sea-level rise on the groundwater system throughout the northern Atlantic Coastal Plain physiographic province because recent analyses of the relative rates of sea-level rise along the Atlantic coast indicate that the Mid-Atlantic region represents a hot spot with anomalously higher rates of sea-level rise than observed elsewhere in the United States. Groundwater levels rose from 0 to 3 ft in response to a 3-ft simulated change in sea-level position, with the largest response occurring along the shoreline and away from non-tidal streams. About 37 percent (or 10,000 square miles) of the area of the northern Atlantic Coastal Plain physiographic province may experience about a 0.5-ft or more increase in water levels with the 3-ft increase in sea-level position, whereas about 18 percent (almost 5,000 square miles) of land of the northern Atlantic Coastal Plain physiographic province may experience a 2-ft or more increase in water levels with the 3-ft increase in sea-level position. These increases in the water table are of particular concern in low-lying areas where the unsaturated (vadose) zone is already thin, thus creating concerns for groundwater inundation of subsurface infrastructure, such as basements, septic systems, and subway systems. This increase in the water table also will likely alter the distribution of groundwater discharge to surface-water bodies thus increasing groundwater flow to streams that would have otherwise discharged directly to coastal waters. Throughout the NACP aquifer system, this redistribution of groundwater discharge results in an additional 2 percent of base flow in streams. Although the increases in groundwater discharge to streams (and corresponding decreases in discharge to coastal waters) calculated for the entire NACP aquifer system and its geographic areas represent only a small increase compared with current [2013] conditions, this redistribution of groundwater discharge from the coast to streams locally can alter the delivery of freshwater input to coastal receiving waters and have ecohydrological implications on the sensitive ecosystems which rely on a balance of groundwater discharge and surface-water flow.

Delaware, Maryland, New Jersey, New York, North Ca↗
USGS2017

Geophysics- and geochemistry-based assessment of the geochemical characteristics and groundwater-flow system of the U.S. part of the Mesilla Basin/Conejos-Médanos aquifer system in Doña Ana County, New Mexico, and El Paso County, Texas, 2010–12

One of the largest rechargeable groundwater systems by total available volume in the Rio Grande/Río Bravo Basin (hereinafter referred to as the “Rio Grande”) region of the United States and Mexico, the Mesilla Basin/Conejos-Médanos aquifer system, supplies water for irrigation as well as for cities of El Paso, Texas; Las Cruces, New Mexico; and Ciudad Juárez, Chihuahua, Mexico. The U.S. Geological Survey in cooperation with the Bureau of Reclamation assessed the groundwater resources in the Mesilla Basin and surrounding areas in Doña Ana County, N. Mex., and El Paso County, Tex., by using a combination of geophysical and geochemical methods. The study area consists of approximately 1,400 square miles in Doña Ana County, N. Mex., and 100 square miles in El Paso County, Tex. The Mesilla Basin composes most of the study area and can be divided into three parts: the Mesilla Valley, the West Mesa, and the East Bench. The Mesilla Valley is the part of the Mesilla Basin that was incised by the Rio Grande between Selden Canyon to the north and by a narrow valley (about 4 miles wide) to the southeast near El Paso, Tex., named the Paso del Norte, which is sometimes referred to in the literature as the “El Paso Narrows.” Previously published geophysical data for the study area were compiled and these data were augmented by collecting additional geophysical and geochemical data. Geophysical resistivity measurements from previously published helicopter frequency domain electromagnetic data, previously published direct-current resistivity soundings, and newly collected (2012) time-domain electromagnetic soundings were used in the study to detect spatial changes in the electrical properties of the subsurface, which reflect changes that occur within the hydrogeology. The geochemistry of the groundwater system was evaluated by analyzing groundwater samples collected in November 2010 for physicochemical properties, major ions, trace elements, nutrients, pesticides (reported but not used in the assessment), and environmental tracers. The data obtained from these samples (with the exception of the pesticide data) were used to gain insights into processes controlling the groundwater movement through the groundwater system in the study area. Results from the geophysical and geochemical assessments facilitated the interpretation of the geochemical characteristics of the groundwater sources and geochemical groups within the groundwater system. The groundwater-flow system in the study area consists primarily of the Mesilla Basin aquifer system, which can be divided into four hydrogeologic units by using an informal classification scheme based on basin-fill stratigraphy and sedimentology with an emphasis on aquifer characteristics. The four hydrogeologic units are (1) the Rio Grande alluvium, which is the shallow aquifer of the Mesilla Basin within the confines of the Mesilla Valley, and the three hydrogeologic units that compose the Santa Fe Group: (2) the lower part of the Santa Fe Group, which is the least productive zone, (3) the middle part of the Santa Fe Group, which is the primary water-bearing hydrogeologic unit in the basin and is generally saturated, and (4) the upper part of the Santa Fe Group, which is the most productive water-bearing unit within the Santa Fe Group but is only partially saturated in the north and largely unsaturated in the south and western parts of the Mesilla Basin. The helicopter frequency domain electromagnetic survey results indicated that approximately half of the resistivity values were less than 10 ohm-meters at depths of 50 and 100 feet with a transition where the resistivity values changed from relatively high values (greater than 20 ohm-meters) to relatively low resistivity values (less than 10 ohm-meters) near Vado, New Mexico. Slightly more than 25 percent of the gridded resistivity values from the three-dimensional grid of the combined inverse modeling results of the direct-current resistivity and time-domain electromagnetic soundings were equal to or less than 10 ohm-meters with large regions of low resistivity becoming apparent in the southernmost part of the study area near the Paso Del Norte where these low resistivity features are spatially the widest at or below the top of the bedrock. These low resistivity values might represent clayey deposits, sediments composed largely of sand and gravel saturated with saline water, or both. Historical dissolved-solids-concentration data within the surface geophysical subset area of the study area were compiled and compared to the inverse modeling results of the combined direct-current resistivity and time-domain soundings; this comparison was done to strengthen the interpretation made from the combined inverse modeling results that the low resistivity features were representative of sand and gravel deposits saturated with saline water and not clayey deposits. Water-level altitudes within the Rio Grande alluvium generally decreased from north to south, with a west to east decrease in water-level altitudes near Las Cruces, New Mexico, as a result of groundwater pumping. Groundwater flow within the Santa Fe Group is more complex than the groundwater flow within the Rio Grande alluvium because of the larger lateral and vertical extent of the Santa Fe Group compared to the Rio Grande alluvium. Groundwater from the Organ Mountains flows directly south towards the Paso del Norte. Groundwater from the Robledo Mountains, the Rough and Ready Hills, and the Sleeping Lady Hills generally flows to the southeast. Groundwater flowing near the north end of the midbasin uplift generally continues east towards the Rio Grande and then flows south on the east side of the midbasin uplift. Groundwater flowing near the west side of the midbasin uplift generally continues south parallel to the faults that make up the midbasin uplift and then flows east towards the Paso del Norte when it reaches the south end of the midbasin uplift. Groundwater from the Aden Hills and the East and West Potrillo Mountains flows to the south end of the midbasin uplift and then continues east towards the Paso del Norte. Throughout most of the Mesilla Valley, the vertical hydraulic gradient was downward because the water-level altitude in the Rio Grande alluvium was higher than it was in the Santa Fe Group, but in some areas (typically in the middle and southern parts of the Mesilla Valley), the vertical hydraulic gradient was substantially reduced or even reversed to an upward hydraulic gradient. The geochemistry data indicate that there was a complex system of multiple geochemical endmembers and mixing between these endmembers with recharge to the Rio Grande alluvium and Santa Fe Group composed mostly of seepage from the Rio Grande, inflows from deeper or neighboring water systems, and mountain-front recharge. Five distinct geochemical groups were identified in the Mesilla Basin study area: (1) ancestral Rio Grande (pre-Pleistocene) geochemical group, (2) modern Rio Grande (Pleistocene to present) geochemical group, (3) mountain-front geochemical group, (4) deep groundwater upwelling geochemical group, and (5) unknown freshwater geochemical group. The ancestral Rio Grande groundwater was water that recharged into the system as seepage losses from the ancestral Rio Grande; this groundwater generally flows from north to south-southeast towards the Paso del Norte. Groundwater on the west side of the midbasin uplift generally flows south until it reaches the southern part of the study area; from the southern part of the study area, the groundwater flows east towards the Paso del Norte. Groundwater on the east side of the midbasin uplift flows south-southeast towards the Paso del Norte where it mixes with groundwater from the modern Rio Grande, uplifted areas in the west, and the deep saline source. The water type of the modern Rio Grande geochemical group ranged from calcium-sulfate water type in the northern part of the study area to sodium-chloride-sulfate water type in the southern part of the study area; from north to south there was a substantial increase in specific conductance, strontium-87/strontium-86 ratio, potassium, and the trace metals of iron and lithium, changing the water chemistry such that it became similar to the water chemistry of the deep groundwater upwelling geochemical group. From age-dating results, water in the modern Rio Grande geochemical group was recharged to the Rio Grande alluvium within the past 10 years. The mountain-front geochemical group was generally old water (apparent age was greater than 10,000 carbon-14 years before present) that was somewhat mineralized and has relatively high concentrations of fluoride and silica, which might indicate longer exposure to volcanic and siliciclastic rocks or aluminosilicate minerals. There were five different locations of recharge determined from the groundwater geochemistry within the mountain-front geochemical group, all having a slightly different geochemical signature: (1) the Rough and Ready Hills, Robledo Mountains, and the Sleeping Lady Hills, (2) the Doña Ana Mountains, (3) the Aden Hills and West Potrillo Mountains, (4) the East Potrillo Mountains, and (5) the Sierra Juárez in Mexico. The groundwater from the Rough and Ready Hills, Robledo Mountains, the Sleeping Lady Hills, and the Doña Ana Mountains generally flows toward the Rio Grande and eventually mixes together and with the modern Rio Grande groundwater. The groundwater originating from the Aden Hills and East and West Potrillo Mountains generally flows east to southeast at a slow rate and eventually mixes and continues east, where it mixes with groundwater from the ancestral Rio Grande geochemical group and with the groundwater from the Sierra Juárez. The groundwater from the Sierra Juárez flows north and then east towards the Paso del Norte where it mixes with groundwater from the uplifted areas in the west, ancestral and modern Rio Grande groundwater, and the upwelling groundwater from a deep saline source. The deep groundwater upwelling geochemical group had the highest concentrations of bicarbonate, potassium, silica, aluminum, iron, and lithium within the study area, indicating that it had been in contact with carbonate and siliciclastic rocks for a much longer period of time and at higher temperatures compared to the other geochemical groups, and was most likely ancient marine groundwater originating from the Paleozoic and Cretaceous carbonate rocks which was upwelling into the Mesilla Basin aquifer system in the southeastern part of the study area through the extensive fault systems. Direct-current resistivity and time-domain electromagnetic soundings support the interpretation of ancient marine groundwater upwelling into the Mesilla Basin aquifer system, as do the analytical results from wells, and the helicopter frequency domain electromagnetic data collected along the Rio Grande. The hydrogen-2/hydrogen-1 ratio and oxygen-18/oxygen-16 ratio isotopic results for samples in the unknown freshwater geochemical group did not plot on the Rio Grande evaporation line, indicating this group did not have a Rio Grande signature (that is, there was no isotopic evidence of a component of Rio Grande water) and it also had the lowest mineralized content of any geochemical group in the study area.

New Mexico, Texas↗
USGS2014

Groundwater-flow and land-subsidence model of Antelope Valley, California

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley groundwater basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, groundwater provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most groundwater pumping in the valley occurs in the Antelope Valley groundwater basin, which includes the rapidly growing cities of Lancaster and Palmdale. Groundwater-level declines of more than 270 feet in some parts of the groundwater basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may increase reliance on groundwater. In 2011, the Los Angeles County Superior Court of California ruled that the Antelope Valley groundwater basin is in overdraft—groundwater extractions are in excess of the Court-defined safe yield of the groundwater basin. The Court determined that the safe yield of the adjudicated area of the basin was 110,000 acre-feet per year (acre-ft/yr). Natural recharge is an important component of total groundwater recharge in Antelope Valley; however, the exact quantity and distribution of natural recharge, primarily in the form of mountain-front recharge, is uncertain, with total estimates ranging from 30,000 to 160,000 acre-ft/yr. Technical experts, retained by parties to the adjudication, used 60,000 acre-ft/yr to estimate the sustainable yield of the basin, and this value was used in this study. In order to better understand the uncertainty associated with natural recharge and to provide a tool to aid in groundwater management, a numerical model of groundwater flow and land subsidence in the Antelope Valley groundwater basin was developed using old and new geohydrologic information. The groundwater-flow system consists of three aquifers: the upper, middle, and lower aquifers. The three aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to groundwater development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Groundwater flowed from the recharge areas to discharge areas around the playas where it discharged from the aquifer system as either evapotranspiration or from springs. Partial barriers to horizontal groundwater flow, such as faults, have been identified in the groundwater basin. Water-level declines owing to groundwater development have eliminated the natural sources of discharge, and pumping for agricultural and urban uses have become the primary source of discharge from the groundwater system. Infiltration of return flow from agricultural irrigation has become an important source of recharge to the aquifer system. The groundwater-flow model of the basin was discretized horizontally into a grid of 130 rows and 118 columns of square cells 1 kilometer (0.621 mile) on a side, and vertically into four layers representing the upper (two layers), middle (one layer), and lower (one layer) aquifers. Faults that were thought to act as horizontal-flow barriers were simulated in the model. The model was calibrated to simulate steady-state conditions, represented by 1915 water levels and transient-state conditions during 1915–95, by using water-level and subsidence data. Initial estimates of the aquifer-system properties and stresses were obtained from a previously published numerical model of the Antelope Valley groundwater basin; estimates also were obtained from recently collected hydrologic data and from results of simulations of groundwater-flow and land-subsidence models of the Edwards Air Force Base area. Some of these initial estimates were modified during model calibration. Groundwater pumpage for agriculture was estimated on the basis of irrigated crop acreage and crop consumptive-use data. Pumpage for public supply, which is metered, was compiled and entered into a database used for this study. Estimated annual agricultural pumpage peaked at 395,000 acre-feet (acre-ft) in 1951 and then declined because of declining agricultural production. Recharge from irrigation return flows was assumed to be 30 percent of agricultural pumpage; delays associated with return flow moving through the unsaturated zone were also simulated. The annual quantity of mountain-front recharge initially was based on estimates from previous studies. The model was calibrated using the PEST software suite; prior information from the area was incorporated through the use of Tikhonov regularization. During model calibration, the estimated mountain-front recharge was reduced from the previous estimate of 30,300 acre-ft/yr to 29,150 acre-ft/yr. Results of the simulations using the calibrated model indicate that simulated groundwater pumpage exceeded recharge in most years, resulting in an estimated cumulative depletion in groundwater storage of 8,700,000 acre-ft during the transient-simulation period (1915–2005). About 15,000,000 acre-ft of cumulative groundwater pumpage was simulated during the transient-simulation period (1915–2005), reaching a maximum rate of about 400,000 acre-ft/yr in 1951. Groundwater pumpage resulted in simulated hydraulic heads declining by more than 150 feet (ft) compared to 1915 conditions in agricultural areas. The decline in hydraulic head in the groundwater basin is the result of this depletion of groundwater storage. In turn, the simulated decline in hydraulic head in the groundwater basin has resulted in the decrease in natural discharge from the basin and has caused compaction of aquitards, resulting in land subsidence. The areal distribution of total simulated land subsidence for 2005, after about 90 years of groundwater development, indicates that land subsidence occurred throughout almost the entire Lancaster subbasin, with a maximum of about 9.4 ft in the central and eastern parts of the subbasin. An important objective of this study was to systematically address the uncertainty in estimates of natural recharge and related aquifer parameters by using the groundwater-flow and land-subsidence model with observational data and expert knowledge. After the model was calibrated to the observations and a reasonable parameter set obtained, the parameter null space—parameter values that do not appreciably affect the model calibration but may have importance for prediction—was identified. The effect of parameter uncertainty on the estimation of mountain-front recharge was addressed using the Null-Space Monte Carlo method. The Pareto trade-off method of visualizing uncertainty was also used to portray the reasonableness of larger natural-recharge rates. Results indicate that the total mountain-front recharge likely ranges between 28,000 and 44,000 acre-ft/yr, which is appreciably less than published estimates of 60,000 acre-ft/yr. Additionally, expected errors associated with agricultural pumpage estimates used in this study were found to have relatively little effect on the estimates of mountain-front recharge, reflecting the difficulty in increasing recharge through manipulation of other components of the water budget. The calibrated model was used to simulate the response of the aquifer to potential future pumping scenarios: (1) no change in the distribution of pumpage, or status quo; (2) redistribution of pumpage; and (3) artificial recharge. All three of these scenarios specify a total pumpage throughout the Antelope Valley of 110,000 acre-ft/yr according to the safe yield value ruled by the Los Angeles County Superior Court of California. This reduction in groundwater pumpage is assumed uniform throughout the basin, based on a 10-percent reduction of the total pumpage in 2005 to achieve the 110,000 acre-ft/yr level. The calibrated Antelope Valley groundwater-flow and land-subsidence model was used to simulate the hydrologic effects of the three groundwater-management scenarios during a 50-year period by using the reduced, temporally constant, pumpage distribution. Results from the first scenario indicated that the total drawdown observed since predevelopment would continue, with values exceeding 325 ft near Palmdale; consequently, land subsidence would also continue, with additional subsidence (since 2005) exceeding 3 ft in the central part of the Lancaster subbasin. The second scenario evaluated redistributing pumpage from areas in the Lancaster subbasin where simulated hydraulic-head declines were the greatest to areas where declines were smallest. Neither a formal optimization algorithm nor water-rights allocations were considered when redistributing the pumpage. Results indicated that hydraulic heads near Palmdale, where the pumpage was reduced, would recover by about 200 ft compared to 2005 conditions, with only 30 ft of additional drawdown in the northwestern part of the Lancaster subbasin, where the pumpage was increased. The magnitude of the simulated additional land subsidence decreased slightly compared to the first, status quo, scenario but land subsidence continued to be simulated throughout most of the northern part of the Lancaster subbasin. The third scenario consisted of two artificial-recharge simulations along the Upper Amargosa Creek channel and at a site located north of Antelope Buttes. Results indicate that applying artificial recharge at these sites would yield continued drawdowns and associated land subsidence. However, the magnitudes of drawdown and subsidence would be smaller than those simulated in the status quo scenario, indicating that artificial-recharge operations in the Antelope Valley could be expected to reduce the magnitude and extent of continued water-level declines and associated land subsidence.

California↗
USGS2016

Groundwater and surface-water interaction, water quality, and processes affecting loads of dissolved solids, selenium, and uranium in Fountain Creek, near Pueblo, Colorado, 2012–2014

In 2012, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group, initiated a study of groundwater and surface-water interaction, water quality, and loading of dissolved solids, selenium, and uranium to Fountain Creek near Pueblo, Colorado, to improve understanding of sources and processes affecting loading of these constituents to streams in the Arkansas River Basin. Fourteen monitoring wells were installed in a series of three transects across Fountain Creek near Pueblo, and temporary streamgages were established at each transect to facilitate data collection for the study. Groundwater and surface-water interaction was characterized by using hydrogeologic mapping, groundwater and stream-surface levels, groundwater and stream temperatures, vertical hydraulic-head gradients and ratios of oxygen and hydrogen isotopes in the hyporheic zone, and streamflow mass-balance measurements. Water quality was characterized by collecting periodic samples from groundwater, surface water, and the hyporheic zone for analysis of dissolved solids, selenium, uranium, and other selected constituents and by evaluating the oxidation-reduction condition for each groundwater sample under different hydrologic conditions throughout the study period. Groundwater loads to Fountain Creek and in-stream loads were computed for the study area, and processes affecting loads of dissolved solids, selenium, and uranium were evaluated on the basis of geology, geochemical conditions, land and water use, and evapoconcentration. During the study period, the groundwater-flow system generally contributed flow to Fountain Creek and its hyporheic zone (as a single system) except for the reach between the north and middle transects. However, the direction of flow between the stream, the hyporheic zone, and the near-stream aquifer was variable in response to streamflow and stage. During periods of low streamflow, Fountain Creek generally gained flow from groundwater. However, during periods of high streamflow, the hydraulic gradient between groundwater and the stream temporarily reversed, causing the stream to lose flow to groundwater. Concentrations of dissolved solids, selenium, and uranium in groundwater generally had greater spatial variability than surface water or hyporheic-zone samples, and constituent concentrations in groundwater generally were greater than in surface water. Constituent concentrations in the hyporheic zone typically were similar to or intermediate between concentrations in groundwater and surface water. Concentrations of dissolved solids, selenium, uranium, and other constituents in groundwater samples collected from wells located on the east side of the north monitoring well transect were substantially greater than for other groundwater, surface-water, and hyporheic-zone samples. With one exception, groundwater samples collected from wells on the east side of the north transect exhibited oxic to mixed (oxic-anoxic) conditions, whereas most other groundwater samples exhibited anoxic to suboxic conditions. Concentrations of dissolved solids, selenium, and uranium in surface water generally increased in a downstream direction along Fountain Creek from the north transect to the south transect and exhibited an inverse relation to streamflow with highest concentration occurring during periods of low streamflow and lowest concentrations occurring during periods of high streamflow. Groundwater loads of dissolved solids, selenium, and uranium to Fountain Creek were small because of the small amount of groundwater flowing to the stream under typical low-streamflow conditions. In-stream loads of dissolved solids, selenium, and uranium in Fountain Creek varied by date, primarily in relation to streamflow at each transect and were much larger than computed constituent loads from groundwater. In-stream loads generally decreased with decreases in streamflow and increased as streamflow increased. In-stream loads of dissolved solids and selenium increased between the north and middle transects but generally decreased between the middle and south transects. By contrast, uranium loads generally decreased between the north and middle transects but increased between the middle and south transects. In-stream load differences between transects appear primarily to be related to differences in streamflow. However, because groundwater typically flows to Fountain Creek under low-flow conditions, and groundwater has greater concentrations of dissolved solids, selenium, and uranium than surface water in Fountain Creek, increases in loads between transects likely are affected by inflow of groundwater to the stream, which can account for a substantial proportion of the in-stream load difference between transects. When loads decreased between transects, the primary cause likely was decreased streamflow as a result of losses to groundwater and flow through the hyporheic zone. However, localized groundwater inflow likely attenuated the magnitude by which the in-stream loads decreased. The combination of localized soluble geologic sources and oxic conditions likely is the primary reason for the occurrence of high concentrations of dissolved solids, selenium, and uranium in groundwater on the east side of the north monitoring well transect. To evaluate conditions potentially responsible for differences in water quality and redox conditions, physical characteristics such as depth to water, saturated thickness, screen depth below the water table, screen height above bedrock, and aquifer hydraulic conductivity were compared by using Wilcoxon rank-sum tests. Results indicated no significant difference between depth to water, screen height above bedrock, and hydraulic conductivity for groundwater samples collected from wells on the east side of the north transect and groundwater samples from all other wells. However, saturated thickness and screen depth below the water table both were significantly smaller for groundwater samples collected from wells on the east side of the north transect than for groundwater samples from other wells, indicating that these characteristics might be related to the elevated constituent concentrations found at that location. Similarly, saturated thickness and screen depth below the water table were significantly smaller for groundwater samples under oxic or mixed (oxic-anoxic) conditions than for those under anoxic to suboxic conditions. The greater constituent concentrations at wells on the east side of the north transect also could, in part, be related to groundwater discharge from an unnamed alluvial drainage located directly upgradient from that location. Although the quantity and quality of water discharging from the drainage is not known, the drainage appears to collect water from a residential area located upgradient to the east of the wells, and groundwater could become concentrated in nitrate and other dissolved constituents before flowing through the drainage. High levels of nitrate, whether from anthropogenic or natural geologic sources, could promote more soluble forms of selenium and other constituents by affecting the redox condition of groundwater. Whether oxic conditions at wells on the east side of the north transect are the result of physical characteristics or of groundwater inflow from the alluvial drainage, the oxic conditions appear to cause increased dissolution of minerals from the shallow shale bedrock at that location. Because ratios of hydrogen and oxygen isotopes indicate evaporation likely has not had a substantial effect on groundwater, constituent concentrations at that location likely are not the result of evapoconcentration.

Colorado↗
USGS2017

Hydrogeology and simulated groundwater flow and availability in the North Fork Red River aquifer, southwest Oklahoma, 1980–2013

On September 8, 1981, the Oklahoma Water Resources Board established regulatory limits on the maximum annual yield of groundwater (343,042 acre-feet per year) and equal-proportionate-share (EPS) pumping rate (1.0 acre-foot per acre per year) for the North Fork Red River aquifer. The maximum annual yield and EPS were based on a hydrologic investigation that used a numerical groundwater-flow model to evaluate the effects of potential groundwater withdrawals on groundwater availability in the North Fork Red River aquifer. The Oklahoma Water Resources Board is statutorily required (every 20 years) to update the hydrologic investigation on which the maximum annual yield and EPS were based. Because 20 years have elapsed since the final order was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in the North Fork Red River aquifer in Oklahoma. This report describes a hydrologic investigation of the North Fork Red River aquifer that includes an updated summary of the aquifer hydrogeology. As part of this investigation, groundwater flow and availability were simulated by using a numerical groundwater-flow model. The North Fork Red River aquifer in Beckham, Greer, Jackson, Kiowa, and Roger Mills Counties in Oklahoma is composed of about 777 square miles (497,582 acres) of alluvium and terrace deposits along the North Fork Red River and tributaries, including Sweetwater Creek, Elk Creek, Otter Creek, and Elm Fork Red River. The North Fork Red River is the primary source of surface-water inflow to Lake Altus, which overlies the North Fork Red River aquifer. Lake Altus is a U.S. Bureau of Reclamation reservoir with the primary purpose of supplying irrigation water to the Lugert-Altus Irrigation District. A hydrogeologic framework was developed for the North Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of a numerical groundwater-flow model of the North Fork Red River aquifer described in this report. A conceptual model of aquifer inflows and outflows was developed for the North Fork Red River aquifer to constrain the construction and calibration of a numerical groundwater-flow model that reasonably represented the groundwater-flow system. The conceptual-model water budget estimated mean annual inflows to and outflows from the North Fork Red River aquifer for the period 1980–2013 and included a sub-accounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from Lake Altus. The numerical groundwater-flow model simulated the period 1980–2013 and was calibrated to water-table-altitude observations at selected wells, monthly base flow at selected streamgages, net streambed seepage as estimated for the conceptual model, and Lake Altus stage. Groundwater-availability scenarios were performed by using the calibrated numerical groundwater-flow model to (1) estimate the EPS pumping rate that guarantees a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical (10-year) drought on base flow and groundwater storage. The results of the groundwater-availability scenarios could be used by the Oklahoma Water Resources Board to reevaluate the maximum annual yield of groundwater from the North Fork Red River aquifer. EPS scenarios for the North Fork Red River 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 0.59, 0.52, and 0.52 acre-foot per acre per year, respectively. Given the 497,582-acre aquifer area, these rates correspond to annual yields of about 294,000, 259,000, and 259,000 acre-feet per year, respectively. Groundwater storage at the end of the 20-year EPS scenario was about 951,000 acre-feet, or about 1,317,000 acre-feet (58 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean water-level decline of about 22 feet. Most areas of the active alluvium near the North Fork Red River, Elk Creek, and Elm Fork Red River remained partially saturated through the end of the EPS scenario because of streambed seepage. Lake Altus storage was reduced to zero after 6–7 years of EPS pumping in each scenario. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the North Fork Red River aquifer and base flows in the North Fork Red River upstream from Lake Altus. The effects of well withdrawals were evaluated by comparing changes in groundwater storage and base flow between four 50-year scenarios using (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2013), (3) 2013 pumping rates, and (4) increasing demand pumping rates. The increasing demand pumping rates assumed a 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwest Oklahoma. Groundwater storage after 50 years with no pumping was about 2,606,000 acre-feet, or 137,000 acre-feet (5.5 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 2.3 feet. Groundwater storage after 50 years with the mean pumping rate for the study period (1980–2013) was about 2,476,000 acre-feet, or about 7,000 acre-feet (0.3 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-level increase of 0.1 foot. Groundwater storage at the end of the 50-year period with 2013 pumping rates was about 2,398,000 acre-feet, or about 70,000 acre-feet (2.8 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.2 feet. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was about 2,361,000 acre-feet, or about 107,000 acre-feet (4.3 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-level decline of 1.8 feet. Mean annual base flow simulated at the Carter streamgage (07301500) on North Fork Red River increased by about 4,000 acre-feet (10 percent) after 50 years with no pumping and decreased by about 5,400 acre-feet (13 percent) after 50 years with increasing demand pumping rates. Mean annual base flow simulated at the North Fork Red River inflow to Lake Altus increased by about 7,400 acre-feet (15 percent) after 50 years with no pumping and decreased by about 5,800 acre-feet (12 percent) after 50 years with increasing demand pumping rates. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage and Lake Altus stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1984–1993). Groundwater storage at the end of the drought period was about 2,271,000 acre-feet, or about 426,000 acre-feet (15.8 percent) less than the groundwater storage of the calibrated numerical model. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 7.1 feet. At the end of the 10-year hypothetical drought period, base flows at the Sweetwater (07301420), Carter (07301500), Headrick (07305000), and Snyder (07307010) streamgages had decreased by about 37, 61, 44, and 45 percent, respectively. The minimum Lake Altus storage simulated during the drought period was 403 acre-feet, which was a decline of 92 percent from the nondrought storage. Reduced base flows in the North Fork Red River were the primary cause of Lake Altus storage declines.

Oklahoma↗
USGS2019

Groundwater movement and interaction with surface water near the confluence of the Platte and Elkhorn rivers, Nebraska, 2016–18

The State of Nebraska requires a sustainable balance between long-term water supplies and uses of groundwater and surface water and requires Natural Resources Districts to include the effect of groundwater use on surface-water systems as part of their respective integrated management plans. Recent droughts in Nebraska (2000–6; 2012–13) have amplified concerns about the long-term sustainability of groundwater and surface-water resources in the state, and concerns about the effect of groundwater irrigation on both streamflow and the water supplies needed to meet wildlife, recreational, and municipal needs. The lower Platte River provides nearly 100 percent of drinking-water supplies to Lincoln, Nebraska, 40 to 60 percent of drinking-water supplies to Omaha, Nebr., and critical aquatic and riparian habitat for threatened and endangered species. The Lower Platte River Basin-wide Management Plan has been jointly developed by the Nebraska Department of Natural Resources and seven Natural Resources Districts to address some of these concerns by managing groundwater and surface-water resources conjunctively. To sustain flows in the lower Platte River that are needed for municipal water supplies, water managers have proposed projects aimed at temporary storage of surface water in upstream parts of the basin to mitigate periods of low flow in the lower Platte River. To increase scientific understanding and provide support for any potential future streamflow augmentation projects, the Papio-Missouri River Natural Resources District, the Lower Platte North Natural Resources District, and the Nebraska Department of Natural Resources, in cooperation with the U.S. Geological Survey, initiated this study to examine groundwater/surface-water interaction along the lower Platte and Elkhorn Rivers upstream from their confluence. The study design described herein focused on understanding seasonal characteristics of groundwater movement and interaction with surface water during periods of high groundwater demand (June through August) and low groundwater demand (all other months). Understanding how groundwater movement and interaction with surface water are affected by streamflow conditions and local groundwater demand is critical to the development of any streamflow augmentation project intended to sustain streamflow and mitigate periods of low flow in the lower Platte River. The characteristics of groundwater movement and interaction with surface water are affected by hydrologic and local climatic conditions. For the study area, 2016–18 conditions can be broadly characterized as above normal precipitation. The flows measured at the Elkhorn River at Waterloo, Nebr., streamflow-gaging station (U.S. Geological Survey station 06800500) were above the long-term median, and the streamflow of the Platte River near Leshara, Nebr., streamflow-gaging station (06796500) remained normal or slightly above normal for the duration of this study. Continuous streamflow and water-level data were interpreted to examine differences in groundwater movement and interaction with surface water between the Platte and Elkhorn Rivers during high and low groundwater demand periods. Although the streamflow for the Platte and Elkhorn Rivers and their tributaries was less during the high groundwater demand period, the hydraulic gradient along a transect of recorder wells was identical (0.0012 foot per foot) during the high and low groundwater demand synoptic water-level and streamflow surveys. The hydraulic gradient between the Platte and Elkhorn Rivers generally remained between 0.0011 and 0.0012 foot per foot. It can be inferred that the hydraulic gradient, which is the only temporally variable factor in Darcy’s Law, is consistent throughout the study period and that groundwater flow does not vary appreciably along this transect. The northern part of the study area (north of the transect of recorder wells) has consistent groundwater and tributary flow from Big Slough, Rawhide Creek (Old Channel), and Rawhide Creek for low and high groundwater demand periods. In the southern part of the study area (south of the transect of recorder wells), tributary flow is more variable and dependent on local groundwater demand and flow conditions of the Platte River. Small decreases (less than 2 feet) in the groundwater levels, such as those measured during the high groundwater demand period, can have substantial changes in the streamflow in an unnamed tributary to the Elkhorn River. The streamflow measured during the high groundwater demand synoptic water-level and streamflow survey had decreased by nearly a factor of 20 when compared to the low groundwater demand period. The volume of groundwater discharge received by the Elkhorn River was estimated by examining the changes in streamflow between measurement locations. Streamflow measurements indicate that the groundwater discharge received by the Elkhorn River in the southern part of the study area was seasonably variable, making it difficult if not impossible to estimate an annual value. In the Elkhorn River, between the Elkhorn River at Waterloo, Nebr., streamflow-gaging station and the Q Street Bridge, streamflow measurements collected during the low groundwater demand period indicated a gain of 80 cubic feet per second, which is comparable to the gain estimated using aerial thermal infrared imagery and water temperature data. Streamflow measurements collected during the high groundwater demand period indicate a loss of 80 cubic feet per second across this same reach. In assessing water supply conditions in the lower Platte River system, the term “loss” in reference to streamflow in the Elkhorn River should be used with caution. Most likely, flow from the Elkhorn River which is “lost” to the groundwater system will later discharge to surface water closer to the confluence of the Platte and Elkhorn Rivers as underflow. A calibrated groundwater flow model of the study area likely is required to predict the fate of this water and to quantify groundwater discharge during varying hydrologic conditions along this reach. Aerial thermal infrared imagery indicated that much of the groundwater discharge in the southern part of the study area is focused across a 3-mile reach where the Elkhorn River turns southwest, perpendicular to the regional groundwater flow direction. Points of focused groundwater discharge were not detected with aerial thermal infrared imagery, indicating that groundwater discharge is diffuse rather than concentrated at focused points. Temperature-based streambed flux estimates indicated that strong regional groundwater gradients are not driving groundwater discharge and hyporheic flow is the dominant groundwater/surface-water exchange process.

Nebraska↗
USGS2014

Evaluation of alternative groundwater-management strategies for the Bureau of Reclamation Klamath Project, Oregon and California

The water resources of the upper Klamath Basin, in southern Oregon and northern California, are managed to achieve various complex and interconnected purposes. Since 2001, irrigators in the Bureau of Reclamation Klamath Irrigation Project (Project) have been required to limit surface-water diversions to protect habitat for endangered freshwater and anadromous fishes. The reductions in irrigation diversions have led to an increased demand for groundwater by Project irrigators, particularly in drought years. The potential effects of sustained pumping on groundwater and surface-water resources have caused concern among Federal and state agencies, Indian tribes, wildlife groups, and groundwater users. To aid in the development of a viable groundwater-management strategy for the Project, the U.S. Geological Survey, in collaboration with the Klamath Water and Power Agency and the Oregon Water Resources Department, developed a groundwater-management model that links groundwater simulation with techniques of constrained optimization. The overall goal of the groundwater-management model is to determine the patterns of groundwater pumping that, to the extent possible, meet the supplemental groundwater demands of the Project. To ensure that groundwater development does not adversely affect groundwater and surface-water resources, the groundwater-management model includes constraints to (1) limit the effects of groundwater withdrawal on groundwater discharge to streams and lakes that support critical habitat for fish listed under the Endangered Species Act, (2) ensure that drawdowns do not exceed limits allowed by Oregon water law, and (3) ensure that groundwater withdrawal does not adversely affect agricultural drain flows that supply a substantial portion of water for irrigators and wildlife refuges in downslope areas of the Project. Groundwater-management alternatives were tested and designed within the framework of the Klamath Basin Restoration Agreement (currently [2013] awaiting authorizing Federal legislation), which would establish a permanent limit on the amount of surface water that can be diverted annually to the Project. Groundwater-management scenarios were evaluated for the period 1970•2004; supplemental groundwater demand by the Project was estimated as the part of irrigation demand that would not have been satisfied by the surface-water diversion allowed under the Klamath Basin Restoration Agreement. Over the 35-year management period, 22 years have supplemental groundwater demand, which ranges from a few thousand acre-feet (acre-ft) to about 100,000 acre-ft in the driest years. The results of the groundwater-management model indicate that supplemental groundwater pumping by the Project can be managed to avoid adverse effects to groundwater discharge that supports critical aquatic habitat. The existing configuration of wells in the Project would be able to meet groundwater-pumping goals in 14 of the 22 years with supplemental groundwater demand; however, substantial irrigation shortages can be expected during drought periods when the demand for supplemental groundwater is highest. The maximum irrigation-season withdrawal calculated by the groundwater-management model is about 60,000 acre-ft, the average withdrawal in drought years is about 54,000 acre-ft, and the amount of unmet groundwater demand reaches a maximum of about 45,000 acre-ft. A comparison of optimized groundwater withdrawals by geographic region shows that the highest annual withdrawals are associated with wells in the Tule Lake and Klamath Valley regions of the Project. The patterns of groundwater withdrawal also show that a substantial amount of the available pumping capacity is unused due to the restrictions imposed by drawdown constraints. Subsequent model applications were used to evaluate the sensitivity of optimization results to various factors. A sensitivity analysis quantified the changes in optimized groundwater withdrawals that result from changes in drawdown-constraint limits. The analysis showed the potential for substantial increases in withdrawals of groundwater with less restrictive drawdown limits at drawdown-control sites in the California part of the model. Systematic variation of the drains-constraint limit yielded a trade-off curve between optimized groundwater withdrawals and the allowable reduction in groundwater discharge to the Project drain system. Additional model applications were used to assess the value of increasing the pumping capacity of the network of wells serving the Project, and the relation between reduced off-Project groundwater pumping and increased pumping by Project irrigators.

California;Oregon↗
USGS2022

Hydrogeology and groundwater quality in the San Agustin Basin, New Mexico, 1975–2019

This report describes the findings of a U.S. Geological Survey study, completed in cooperation with the Bureau of Land Management, focused on better understanding the present-day (1975–2019) hydrogeology and groundwater quality of the San Agustin Basin in west-central New Mexico to support sustainable groundwater resource management. The basin hosts a relatively undeveloped basin-fill and alluvium aquifer system and is topographically divided into east and west subbasins by the McClure Hills. Groundwater chemistry and groundwater elevation data were compiled, collected, and interpreted in the context of groundwater flow and quality. The analyses presented in this report consider groundwater chemistry data collected within the last decade (2010–19) and groundwater elevation data collected from 1975 through 2019 to provide insight into present-day conditions. Groundwater elevations show that groundwater typically moves from the highlands to the lowlands, with a prominent east to west regional trend. Groundwater elevations were lowest in the southwestern portion of the west subbasin, where estimated flow directions suggest underflow through the local highlands into the northern East Fork Gila River watershed, which is further supported by historical groundwater elevation data from the northern East Fork Gila River watershed. Gradual groundwater elevation gradients (about 2 feet per mile) near the east and west subbasin divide suggest that groundwater slowly flows from the east subbasin to the west subbasin. Quantitative analyses of groundwater chemistry data show that groundwater in both subbasins has similar chemical characteristics. A systematic east to west groundwater evolution in water chemistry was not observed despite evidenced subbasin connectivity. The absence of this pattern suggests that groundwater mixing is regionally prevalent, sediment reactivity is low and variable, and (or) recharge conditions are comparable in both subbasins. Groundwater chemistry was generally independent of aquifer type, suggesting that the aquifers are hydrologically well connected. Corrected carbon-14 groundwater age estimates in the basin ranged from 232 to 13,916 years before present with a median of 5,409 years. A wide range of groundwater ages is therefore present in the basin, with waters commonly being thousands of years old, thereby supporting generally slow regional groundwater movement. A component of relatively young groundwater, for which estimated ages could not be accurately computed, is also present in the basin, and it may commonly mix with older waters. The spatial distribution of categorical and quantitative groundwater ages indicates that most recharge likely occurs in the highlands through mountain-block recharge and as focused recharge within arroyos, although evidence of modern (1953 and after) groundwater was minimal at sampled sites. Median annual gradients (groundwater elevation change over time) indicate that most groundwater elevations in the lowlands changed little (−0.2 to 0.2 foot per year) from 1975 through 2019. Groundwater elevations in the highlands varied more annually, which is likely due to recharge from precipitation events. These more variable groundwater elevations in the highlands compared with the lowlands, along with groundwater ages, provide further evidence that most groundwater recharge takes place in the highlands, with minimal recharge in the lowlands. Median groundwater elevation change for all sites was −0.05 foot per year. Temporal consistency of lowland groundwater elevations suggests that regional groundwater dynamics have been more or less stable through time under current climate and development conditions, although median annual gradients indicate that groundwater elevations may have slightly declined on average between 1975 and 2019.

New Mexico↗
USGS2015

Geochemical conditions and the occurrence of selected trace elements in groundwater basins used for public drinking-water supply, Desert and Basin and Range hydrogeologic provinces, 2006-11: California GAMA Priority Basin Project

The geochemical conditions, occurrence of selected trace elements, and processes controlling the occurrence of selected trace elements in groundwater were investigated in groundwater basins of the Desert and Basin and Range (DBR) hydrogeologic provinces in southeastern California as part of the Priority Basin Project (PBP) of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The GAMA PBP is designed to provide an assessment of the quality of untreated (raw) groundwater in the aquifer systems that are used for public drinking-water supply. The GAMA PBP is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey and the Lawrence Livermore National Laboratory. The DBR hydrogeologic provinces consist of 141 defined groundwater basins separated by mountain ranges, faults, and other features. This report presents analyses of data collected from nine study areas within the DBR hydrogeologic provinces: Antelope Valley, Borrego Valley, the Central Desert area, Coachella Valley, Colorado River, Indian Wells Valley, Low-Use Basins of the Mojave and Sonoran Deserts, the Mojave, and Owens Valley. Collectively, these nine study areas are referred to as the DBR study unit. The study unit covers approximately 7,000 square miles and includes the 63 groundwater basins in the DBR hydrogeologic provinces in which groundwater is used for public drinking-water supply. The vast majority of the 223 wells sampled for this study were long-screened production wells used primarily for public supply. Uncorrected carbon-14 ( 14 C) groundwater ages for samples collected in the DBR study unit ranged from less than (<) 100 to 33,700 years before present (BP). Sixty-six percent of sample ages were greater than (>) 100 years BP, and 40 percent were >3,800 years BP. Samples collected from wells located adjacent to mountain-front recharge areas or major surface-water features generally had younger groundwater ages than did samples collected from wells located away from mountain fronts or towards the distal ends of basin groundwater flow paths. Most groundwater sampled in the DBR study unit had alkaline pH: 89 percent of sample pH values ranged from 7.1 to 9.8, with 37 percent greater than or equal to (&ge;) 7.9. Groundwater age was significantly correlated (positively) with pH, likely because silicate weathering is a primary control on groundwater pH and is a slow process. The oxidation-reduction (redox) condition of the groundwater sampled in the DBR study unit was predominantly oxic (71 percent), except in the Colorado River study area where organic-rich fluvial aquifers provide the electron donors necessary to support iron-reducing (anoxic-Fe) redox processes. The cation type of 78 percent of the samples was either sodium- or mixed-type, and the anion type of 83 percent of the samples was either bicarbonate- or mixed-type. Sodium-type groundwaters generally were older and more alkaline than calcium-type groundwaters, consistent with the change in water chemistry expected from cation exchange between groundwater and aquifer sediments over long periods of time. Because of the correlation with young groundwater, calcium-type groundwater was predominantly from wells located adjacent to mountain-front recharge areas. Arsenic (As), boron (B), fluoride (F), molybdenum (Mo), strontium (Sr), uranium (U), and vanadium (V) were selected for assessment in this study because they occurred at concentrations greater than California Department of Public Health or U.S. Environmental Protection Agency regulatory or non-regulatory drinking-water-quality benchmarks in more than 2 percent of the 223 samples collected in the DBR study unit. As and F were detected most commonly (18 and 13 percent, respectively) at concentrations above associated water-quality benchmarks and Sr and V least frequently (both at 3 percent). Given that 14 C groundwater ages are predominantly >100 years BP, land use in the study unit is primarily undeveloped, and chemicals derived from anthropogenic sources, such as volatile organic compounds, were infrequently detected, high concentrations of these trace elements in groundwater were most likely the result of natural factors and not anthropogenic factors. As, F, Mo, and V concentrations showed significant positive correlations to groundwater age and to pH. This relation is partly due to the sources of trace elements likely being the weathering of primary minerals, such as silicate minerals, which is a slow process that takes place over hundreds to thousands of years. This relation also reflects the positive correlation between groundwater age and pH. Geochemical modeling predicted that the dominant species of As, Mo, and V in solution were oxyanions (HAsO 4 2&ndash; , MoO 4 2&ndash; , and H 2 VO 4&ndash; ), which are likely to be mobile in alkaline groundwater because mineral surfaces composing aquifer matrices have a predominantly negative surface charge under alkaline conditions. F also exists predominantly as a negatively charged ion (F &ndash; ). At pH values >7.5, saturation indices generated by the geochemical modeling program PHREEQC indicated that F solubility may be somewhat limited by the precipitation of the mineral fluorapatite [Ca 5 (PO 4 ) 3 F]. Speciation modeling of As in anoxic-Fe groundwater (iron-reducing conditions) showed that samples were supersaturated with orpiment (As 2 S 3 ), indicating that mineral precipitation may be responsible for low As concentrations observed in reducing groundwater. In contrast, U concentrations showed significant negative correlations to groundwater age and to pH. Higher U concentrations generally occurred in samples for which geochemical modeling indicated that the uncharged ternary complex Ca 2 UO 2 (CO 3 ) 3 was the dominant aqueous U species. This uncharged complex is not attracted to the charged surfaces of minerals and thus increases U solubility. Formation of Ca 2 UO 2 (CO 3 ) 3 was greater in younger groundwaters because calcium and uranium concentrations generally were lower in older groundwaters, likely due to cation-exchange processes and precipitation of the mineral calcite as groundwater pH increased. Co-precipitation of U with the calcite (CaCO 3 ) may remove U from the aqueous phase. Saturation indices indicated that the anoxic-Fe groundwaters from the Colorado River study area were supersaturated with the mineral uraninite (UO 2 ), suggesting that UO 2 precipitation may be responsible for the low concentrations of U observed in these samples. Concentrations of strontium, which exists primarily in a cationic form (Sr 2+ ), were not significantly correlated with either groundwater age or pH. Strontium concentrations showed a strong positive correlation with total dissolved solids (TDS). Dissolved constituents, such as Sr, that interact with mineral surfaces through outer-sphere complexation become increasingly soluble with increasing TDS concentrations of groundwater. Boron concentrations also showed a significant positive correlation with TDS, indicating the B may interact to a large degree with mineral surfaces through outer-sphere complexation.

California↗
USGS2015

Initial characterization of the groundwater system near the Lower Colorado Water Supply Project, Imperial Valley, California

In 2009, the U.S. Geological Survey, in cooperation with the city of Needles, began a study of the hydrogeology along the All-American Canal, which conveys water from the Colorado River to the Imperial Valley. The focus of this study was to gain a better understanding of the effect of lining the All-American Canal, and other management actions, on future total dissolved solids concentrations in groundwater pumped by Lower Colorado Water Supply Project wells that is delivered to the All-American Canal. The study included the compilation and evaluation of previously published hydrogeologic and geochemical information, establishment of a groundwater-elevation and groundwater-quality monitoring network, results of monitoring groundwater elevations and groundwater quality from 2009 to 2011, site-specific hydrologic investigations of the Lower Colorado Water Supply Project area, examination of groundwater salinity by depth by using time-domain electromagnetic surveys, and monitoring of groundwater-storage change by using microgravity methods. Prior to the completion of the All-American Canal in 1940, groundwater in the study area flowed from east to west, and groundwater was recharged primarily by underflow from the Colorado River Valley. After construction of the All-American Canal, groundwater elevations were altered in the study area as seepage of Colorado River water from the All-American Canal and other canals became the dominant recharge source. By 2005, groundwater elevations had increased by as much as 50&ndash;70 feet along the All-American Canal. Superimposed on the east-to-west groundwater gradient was groundwater movement away from the All-American Canal to the north and, most likely, to the south into Mexico. After lining the All-American Canal, from 2007 to 2010, groundwater elevations declined as seepage from the All-American Canal decreased. Between 2005 (the last complete groundwater-elevation survey prior to lining the All-American Canal) and 2011, groundwater elevations declined 20&ndash;40 feet along the All-American Canal and as much as 40&ndash;45 feet in the vicinity of Lower Colorado Water Supply Project pumping wells. Water-quality and isotope data were used to differentiate historically recharged groundwater from groundwater more recently recharged by seepage of Colorado River surface water from the All-American Canal. Prior to the completion of the All-American Canal in 1940, groundwater in the southern part of the study area was primarily sodium-chloride/sulfate type water that had relatively low total dissolved solids concentrations (500&ndash;820 milligrams per liter). During 2007&ndash;11, groundwater in the southern part of the study area, near the All-American Canal, ranged from sodium-chloride type water to mixed-cation-sulfate type water that had total dissolved solids concentrations generally less than 879 milligrams per liter. The stable-isotopic signature of groundwater near the All-American Canal sampled in 2009&ndash;11 indicated inputs of Colorado River water that had been affected by evaporation, and radioactive isotopes indicated that a substantial fraction of water had been recharged recently, within the past 60 years. This contrasted with historically recharged groundwater near the All-American Canal, which had higher sodium and chloride concentrations, and lower calcium and sulfate concentrations, than recent recharge from the All-American Canal. Groundwater at a distance from the All-American Canal, in the East Mesa, Algodones Dunes, Pilot Knob Mesa, and Cargo Muchacho Mountains piedmont, was found to have higher total dissolved solids concentrations (generally greater than 1,000 milligrams per liter) than recently recharged groundwater near the All-American Canal. Time-domain electromagnetic data indicated that low-salinity groundwater was present down to about 377 feet below land surface near the All-American Canal; groundwater salinity at depth increased with distance north from the All-American Canal. Groundwater several miles or more from the canal also did not contain tritium and had a residence time on the order of thousands to tens of thousands of years. The groundwater in the piedmont of the Cargo Muchacho Mountains had a distinctly light stable-isotopic signature indicative of recharge by runoff from local precipitation, whereas the stable isotopic signature of groundwater in the East Mesa and the Algodones Dunes indicated a mixture of local precipitation and historic Colorado River recharge sources. During and after lining the All-American Canal (2007&ndash;11), groundwater elevations in the Lower Colorado Water Supply Project area declined, while total dissolved solids concentrations remained relatively constant. The total dissolved solids concentrations in well LCWSP-2 ranged from 650 to 800 milligrams per liter during this study. Depth-specific water-quality and isotope sampling at well LCWSP-2 indicated the groundwater pumped from the deeper part of the screened interval (240&ndash;280 feet below land surface) contained a greater proportion of historical groundwater than the groundwater pumped from the shallower part of the screened interval (350&ndash;385 feet below land surface). Age-tracer data at well LCWSP-2 indicated that all depths of the screened interval had received recent recharge from seepage of Colorado River water from the All-American Canal.

California↗
USGS2018

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

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

California↗
USGS2016

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967&ndash;2011, was divided into two sub-intervals because of varying water use, 1970&ndash;80 and 1981&ndash;2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981&ndash;2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967&ndash;2011; 20,724 acre-ft/yr during 1970&ndash;80, and 13,739 acre-ft/yr during 1981&ndash;2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967&ndash;2011; 19,963 acre-ft/yr during 1970&ndash;80; and 12,285 acre-ft/yr during 1981&ndash;2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967&ndash;2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980&ndash;2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980&ndash;2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7&ndash;279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4&ndash;279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981&ndash;2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994&ndash;2004) during the 1980&ndash;2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗
USGS2020

Updating data inputs, assessing trends, and evaluating a method to estimate probable high groundwater levels in selected areas of Massachusetts

A method to estimate the probable high groundwater level in Massachusetts, excluding Cape Cod and the islands, was developed in 1981. The method uses a groundwater measurement from a test site, groundwater measurements from an index well, and a distribution of high groundwater levels from wells in similar geologic and topographic settings. The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, conducted an update to the Frimpter method for estimating the probable high groundwater levels in Massachusetts. The study evaluated the potential changes to the method resulting from four decades of additional groundwater-level data and the expansion of the network of wells for monitoring groundwater levels. The differences and potential benefits of daily, as opposed to monthly, measurements in the application of the method were examined because of the increased availability of high-frequency (subdaily) groundwater-level data. The study also considered long-term trends in groundwater levels that may alter the accuracy of the method. Finally, the accuracy of the estimated high groundwater levels was evaluated, and improved implementation guidance was prepared. For this study, groundwater levels in 153 wells in Massachusetts and surrounding States with records with lengths of 16 to 78 years were analyzed. The highest recorded groundwater levels ranged from 1.2 feet (ft) above land surface (flooded conditions) to 45.8 ft below land surface, with a median of 4.6 ft below land surface. The maximum annual groundwater-level range was 1.4 to 17.9 ft, with a median of 5.5 ft. The within-month variation, maximum annual groundwater-level range, and highest recorded groundwater level were computed using daily mean groundwater-level values from 28 wells with continuous records. The use of daily data resulted in larger maximum annual groundwater-level ranges (0.02 to 2.94 ft larger, with a median of 0.58 ft larger) and shallower highest-recorded groundwater levels (0.0 to 1.60 ft shallower, with a median of 0.18 ft shallower) than computations based on monthly measurements in the same wells. Statistical tests showed moderate to strong evidence of trends in measurements of both high and low groundwater levels within most of the periods during which water levels were analyzed. High groundwater levels rose beneath the land surface at most sites during four of the six periods used for analysis (1966–2015, 1986–2015, 1991–2010, and 1981–2010). Low groundwater levels also increased at many sites during most of the periods evaluated, but this trend was less widespread than the similar trends in high groundwater levels, and the trend was to deeper low groundwater levels at more sites than the trend to deeper high groundwater levels. There was no clear trend in annual groundwater-level ranges at most sites during the six periods analyzed. In general, the Frimpter method predicted shallower (higher) high groundwater levels than were observed but correctly classified sites according to their suitabilities for unmounded septic systems. The mean error of the predictions (difference between the estimated and observed groundwater levels) ranged from −3.23 ft to −1.40 ft for various approaches to estimating the groundwater-level range and selecting an index well. The method correctly classified 83 to 86 percent of monitoring-well sites according to their suitability for an unmounded septic system for many approaches to estimating the annual groundwater-level range and selecting an index well. The approach selected for estimating the annual groundwater-level range and selecting an index well will depend upon the importance of an accurate estimate of the high groundwater level as compared to the importance of an estimated high groundwater level that is less likely to be exceeded.

Connecticut, Massachusetts, New Hampshire, Rhode I↗
USGS2021

Hydrogeologic framework and groundwater characterization in selected alluvial basins in the upper Rio Grande basin, Colorado, New Mexico, and Texas, United States, and Chihuahua, Mexico, 1980 to 2015

Increasing demand for the limited water resources of the United States continues to put pressure on resource management agencies to balance the competing needs of ecosystem health with municipal, agricultural, and other uses. To meet these needs, the U.S. Geological Survey conducted a multiyear study to evaluate water resources in the upper Rio Grande Basin in the southwestern United States. The upper Rio Grande Basin extends from south-central Colorado, through New Mexico, into west Texas near Fort Quitman, including parts of Chihuahua, Mexico. The upper Rio Grande Basin consists of a sequence of alluvial basins that formed in the Rio Grande rift approximately 30 million years ago. This report describes the hydrogeology of the upper Rio Grande Basin and how the groundwater resources in the basin have changed from 1980 to 2015. The hydrogeologic framework includes the horizontal delineation of the alluvial basins within the upper Rio Grande Basin from the headwaters in Colorado to Fort Quitman, Texas, including part of Mexico. Groundwater-level measurements from existing State and Federal data were used to construct groundwater-level altitude and groundwater-level change maps. Of the 2,699 wells with groundwater-level data used in this study, 1,055 wells had data for only a single 5-year period, 703 wells had data for 50 percent or more of the 35 years of the study, and only 57 wells have 5-year groundwater-level data for the entire study period. The median decline in water levels in the upper Rio Grande Basin was 0.13 foot (ft) per 5-year period, and declines were measured in 53 percent of the 703 wells that contained data for 50 percent or more of the study period. Rates of groundwater-level decline greater than 1 ft per 5-year period were measured in 17 percent of the wells, greater than 2 ft per 5-year period, in 3 percent of the wells, and greater than 3 ft per 5-year period, in 1 percent of the wells. Overall, groundwater levels rose in 6 percent of the 703 wells that contained data for 50 percent or more of the study period, and in 4 percent of the wells, groundwater levels rose by 1 ft or more per 5-year period. Groundwater-level changes in wells with consecutive 5-year measurement periods exhibited the most variability in the Española, Middle Rio Grande, and Mesilla/Conejos-Médanos alluvial basins. The largest declines in groundwater-level altitudes in individual wells were observed in the Española alluvial basin during 1995–2000, in the Palomas alluvial basin during 2010–2015, and in the Jornada del Muerto alluvial basin during 2005–10. The largest rises in groundwater-level altitudes in individual wells were observed in the Española alluvial basin during 2005–10, in the Middle Rio Grande alluvial basin during 1995–2000, and in the Mesilla/Conejos-Médanos alluvial basin during 1980–85. Changes in groundwater storage throughout the study period varied by alluvial basin, likely based largely on changes in groundwater withdrawals because of increased demands during drier periods and population growth. All alluvial basins except the Tularosa-Hueco alluvial basin were evaluated for changes in groundwater storage from 1980 to 2015. Extremely limited data availability in 2010–15 for the Tularosa-Hueco alluvial basin led to this 5-year period being dropped from the groundwater-level change map and storage analysis for this basin. In the San Luis Valley in southern Colorado, efforts to reverse groundwater depletion in the unconfined aquifer recovered approximately 250,000 acre-feet in storage between late 2013 and early 2018, following the implementation of a “pay-to-pump” groundwater program. However, severe drought that persists in the upper Rio Grande Basin, particularly in southern Colorado, has undone some of the conservation efforts. Within the Española alluvial basin, groundwater storage varied because municipal demand increased the demand on groundwater resources and conservation efforts were implemented. A groundwater-flow model evaluated for the Española alluvial basin indicated declines in groundwater storage from 1947 through 1982. Groundwater storage decreased in the Española alluvial basin in 1980–85, 1985–90, 1990–95, 1995–2000, and 2005–10 and increased in 2000–05 and 2010–15 leading to groundwater storage in 2015 about even with that in 1985. Based on gridded groundwater-level altitudes, groundwater storage decreased in the Middle Rio Grande Basin from 1980 to 2015, except for during the 1980–85, 2000–05, and 2010–15 periods with an overall cumulative storage decrease from 1980 to 2015. Groundwater-flow models evaluated for the Middle Rio Grande alluvial basin showed groundwater storage in the Middle Rio Grande alluvial basin has been reduced since the mid-1950s through the end of the study period except for a brief recovery (reduction in storage outflow) in the mid-1980s. Simulated groundwater storage has also decreased in parts of the Palomas and Mesilla/Conejos-Médanos alluvial basins, and the northern part of the Conejos-Médanos alluvial basin starting in 1995 (excluding 2005 and 2007) and in the Tularosa-Hueco alluvial basin from the early 1940s to the end of the study period. Groundwater storage increased in the Mesilla/Conejos-Médanos alluvial basin during 1980–85 and slightly during 1990–95 and then decreased in the other 5-year periods. Groundwater storage in the Tularosa-Hueco alluvial basin increased from 1985 to 1990, but otherwise decreased, leading to an overall net groundwater-level decline in this part of the basin from 1980 to 2010.

Colorado, New Mexico, Texas, Chihuahua↗
USGS2018

Groundwater discharge characteristics for selected streams within the Loup River Basin, Nebraska, 2014–16

Streams in the Loup River Basin are sensitive to groundwater withdrawals because of the close hydrologic connection between groundwater and surface water. Groundwater discharge is the primary component of streamflow in the Loup River Basin and constitutes more than 90 percent of streamflow in the central part of the Sand Hills. To improve the understanding of geologic controls and various climatic and land-use changes on groundwater discharge, the U.S. Geological Survey (USGS), in cooperation with the Upper Loup Natural Resources District (NRD), the Lower Loup NRD, and the Nebraska Environmental Trust, studied the spatial and temporal characteristics of groundwater discharge within the Loup River Basin. This report documents the methods of data collection and analysis, which include the collection of approximately 350 river miles of aerial thermal infrared imagery and continuous groundwater-level and temperature data from six streamflow-gaging stations within the Loup River Basin. The results from the stream reconnaissance and examination of aerial thermal infrared imagery demonstrated the influence of the surficial and subsurface geology on the spatial characteristics of groundwater discharge to streams in the Loup River Basin. At the headwaters of the South Loup River, streamflow is sustained and increased from focused groundwater discharge emanating from Quaternary deposits at many small (less than 0.1 cubic foot per second) focused points. The volume of water produced from this dense network of focused groundwater discharge points along the North Fork South Loup River is sufficient to provide approximately 40 percent of the flow measured at the South Loup River at Arnold, Nebraska streamflow-gaging station (USGS station 06781600) during the irrigation season. Approximately 5 miles downstream from the South Loup River at Arnold, Nebr., streamflow-gaging station, the river incises into Pliocene-age sand and gravel deposits, which provide additional groundwater discharge to the stream. The streamflow of the South Loup River increases by a factor of 5 across a 62-mile reach of the middle South Loup River. Increases in streamflow along the upper Dismal River result from a dense network of focused groundwater discharge points within semiconsolidated Pliocene-age deposits. Below the Dismal River near Thedford, Nebr., streamflow-gaging station (USGS station 06775900), the Dismal River incises into the Ogallala Formation over a short reach before flowing over coarser, more permeable Quaternary-age alluvial deposits. Diffuse groundwater discharge sustains and increases the streamflow of the lower Dismal River in this reach. Groundwater sapping was evident on some stream reaches and has increased the size and flow of focused groundwater discharge points. Previous researchers have documented streambed incision and groundwater sapping on the upper Dismal River that have created and enlarged focused groundwater discharge points capturing additional groundwater. Similar processes appear to have played a role in the formation of larger focused groundwater discharge points, which sustain the flow of the middle South Loup River. The constant flow of groundwater into the South Loup River has removed finer-grained Quaternary sediments and further exposed Pliocene-age gravel deposits. Headward erosion is evident where some of the large focused groundwater discharge points have incised their own draws and terminate in bowl-like depressions away from the stream. Within the Loup River NRDs, the percentage of groundwater-irrigated land in a stream basin is one factor that affects groundwater discharge to streams. A striking example was at the South Loup River at Saint Michael, Nebr., groundwater and streamflow-gaging station (USGS station 06784000) where the shallow groundwater levels declined below the level of the stream during the middle to late part of the growing season (July to September) when consumptive groundwater use was at its peak. The South Loup River Basin above the South Loup River at Saint Michael, Nebr., streamflow-gaging station has the highest percentage of groundwater-irrigated row crops of all the basins examined in this study. Continuous groundwater and surface-water levels measured at the North Loup River at the Taylor, Nebr., streamflow-gaging station (USGS station 06786000) indicate that the stream is receiving groundwater throughout the year; however, when consumptive groundwater use peaks during the middle to late part of the growing season (July to September), the difference in elevation between the groundwater level and the stream elevation decreases, which indicates a reduction in the amount of groundwater discharge received.

Nebraska↗
USGS2021

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗