USGS Science⌕ Search

SEARCH · USGS Science

Results for “groundwater aquifer”

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 541 records · Page 30Linked to original sources

Hydrogeology and simulated effects of ground-water development on an unconfined aquifer in the Closed Basin Division, San Luis Valley, Colorado

Wells completed in an unconfined aquifer in the Closed Basin Division of the San Luis Valley Project, Colorado, are expected to provide about 101,800 acre-ft of groundwater/year to the Rio Grande when this project is completed. Lowering of groundwater levels in the unconfined aquifer is expected to decrease the quantity of groundwater that is lost by evapotranspiration. The aquifer system, which consists of an unconfined aquifer that is 50 to 130 ft thick, overlies a thick, leaky confined aquifer. Groundwater moves from the edge of the valley toward a topographic low near the center of the Closed Basin Division, where it is lost by evapotranspiration. A two-dimensional groundwater flow model was used to evaluate the effects of projected withdrawal of about 141 cu ft/sec by 168 wells throughout a 20-year period. The simulated pumpage resulted in a projected drawdown greater than 0.1 ft in the water-levels of the unconfined aquifer over an area of about 370 sq mi. Maximum simulated drawdown was 25 ft. Simulations indicate that about 66 % of the water to be withdrawn from the unconfined aquifer would be derived from decreases of evapotranspiration, 26% from induced leakage from an underlying confined aquifer, and 8% from storage of the unconfined aquifer. Model simulations were based only on withdrawals from wells completed in the unconfined aquifer. Pumpage from the confined aquifer was not simulated. Upward leakage from the confined aquifer predicted by the model, results from the simulated declines of the potentiometric surface in the unconfined aquifer. (USGS)

Water-Resources Investigations Report↗

Groundwater chemistry and water-level elevations in bedrock aquifers of the Piceance and Yellow Creek watersheds, Rio Blanco County, Colorado, 2013–16

The Piceance and Yellow Creek watersheds in Rio Blanco County, Colorado, are known to contain important energy resources (oil shale and natural gas) and mineral resources (nahcolite). The primary sources of fresh groundwater in the Piceance and Yellow Creek watersheds are bedrock aquifers in the Uinta and Green River Formations. The aquifers are divided into an upper and lower aquifer separated by a regionally extensive semiconfining layer. These aquifers provide water to streams and springs in the watersheds and are an important resource to people living and working in the area. Development of energy and mineral resources has the potential to affect the quality of groundwater in several ways. The Bureau of Land Management and the U.S. Geological Survey began groundwater monitoring in 2010 to characterize the groundwater quality and water-level elevations of shallow bedrock aquifers in the Piceance and Yellow Creek watersheds. The purpose of this report is to present ground-water chemistry and water-level elevations collected during 2013–16. Comparisons are made to data that were collected from the bedrock aquifers from 2010 to 2012 to identify the potential for changes in water quality and water-level elevations. Appreciable changes in water-level elevations and hydraulic gradient were observed in early April 2015 in two wells completed in the upper and lower aquifers. The hydraulic gradient between the two wells was consistently downward from the upper aquifer to the lower aquifer during 2010–15; however, in early April 2015, the gradient changed from downward to upward between the two aquifers. Overall, water-level elevations declined by about 14 and 11 feet in the upper and lower aquifers, respectively, from 2013 to 2016. Previously published data estimated groundwater ages at 1,200 years old in the upper aquifer and 9,600 years old in the lower aquifer. These groundwater ages indicate that ground-water was recharged over thousands of years. With such long periods of time for aquifer recharge, declines in water-level elevation over short time steps (a few months) have important implications for sustainable management of this resource. Solution mining activities or drilling for oil and natural gas in the area could be related to the changes observed in water-level elevations in these wells; however, further investigation would be needed to evaluate causation. Changes in major-ion chemistry were evaluated in the bedrock aquifer using time series plots of select major-ion data from 2010 to 2016. Major-ion chemistry was variable for a single well from 2010 to 2016 where alkalinity and sulfate were the most variable constituents. One possible explanation for the observed changes in major-ion chemistry may be that the sample depth for that well no longer represents the most appreciable flow in the borehole. On a larger scale, potential changes in flow within the borehole may indicate changes in the regional flow system. Methane and volatile organic compound concentrations were evaluated using a similar approach to that of major ions and had similar findings. Methane concentrations in wells sampled from 2010 to 2016 were generally constant. The only exception was observed at a single well where the range of methane concentrations was from 57.4 (2010) to 4.02 milligrams per liter (2013). This is the same well where changes in water-level elevation, hydraulic gradient, and major-ion chemistry were observed, providing multiple lines of evidence to indicate change in the bedrock aquifers. Sampling of a well located in an area with little energy development but where faults or fractures could provide a path for the migration of fluids indicate mixing of groundwater between the upper and lower aquifers.

Colorado↗

Hydrogeologic framework of the uppermost principal aquifer systems in the Williston and Powder River structural basins, United States and Canada

The glacial, lower Tertiary, and Upper Cretaceous aquifer systems in the Williston and Powder River structural basins within the United States and Canada are the uppermost principal aquifer systems and most accessible sources of groundwater for these energy-producing basins. The glacial aquifer system covers the northeastern part of the Williston structural basin. The lower Tertiary and Upper Cretaceous aquifer systems are present in about 91,300 square miles (mi 2 ) of the Williston structural basin and about 25,500 mi 2 of the Powder River structural basin. Directly under these aquifer systems are 800 to more than 3,000 feet (ft) of relatively impermeable marine shale that serves as a basal confining unit. The aquifer systems in the Williston structural basin have a shallow (less than 2,900 ft deep), wide, and generally symmetrical bowl shape. The aquifer systems in the Powder River structural basin have a very deep (as much as 8,500 ft deep), narrow, and asymmetrical shape. The Williston structural basin has been an important oil and natural gas producing region since the 1950s, and production has increased substantially since the mid-2000s due to improved drilling and hydraulic fracturing methods from deep formations, such as the Bakken and Three Forks Formations. These improved methods require considerable volumes of freshwater mostly from shallow aquifers or surface water. Coal, lignite, and coal-bed natural gas are additional sources of energy in both basins that can affect the quality and quantity of shallow aquifers through strip mining and groundwater depletion. In 2011, the U.S. Geological Survey initiated a regional study of the glacial, lower Tertiary, and Upper Cretaceous aquifer systems in the Williston and Powder River structural basins with the goal to quantify groundwater availability. This report, together with a companion report of the conceptual flow model, provides an improved understanding of the groundwater flow systems and a basis for a numerical, regional groundwater-flow model. This study combines the lithostratigraphic units of the glacial, lower Tertiary, and Upper Cretaceous aquifer systems in the United States and Canada into 7 regional hydrogeologic units—glacial deposits, 4 bedrock aquifers, and 2 bedrock confining units—using general hydraulic properties. The glacial deposits are composed of till and glacial outwash sands and gravels with areas of cobbles and boulders. The four bedrock aquifers are the upper Fort Union, lower Fort Union, lower Hell Creek, and Fox Hills aquifers and are contained primarily in sandstone layers. The two confining units are the middle Fort Union hydrogeologic unit (shale) and upper Hell Creek hydrogeologic unit (contains less sandstone than the underlying lower Hell Creek aquifer). Water from hydrogeologic units in these three aquifer systems is relatively fresh and potable, whereas withdrawals seldom occur from units below the basal confining unit because of great depths (greater than 800 ft) and poor water quality. Analysis of about 300 electric (resistivity) and lithologic logs in the Williston structural basin and numerous existing publications for the Powder River structural basin were used to develop a three-dimensional hydrogeologic framework for both basins. Interpolated thicknesses of the glacial deposits, the lower Tertiary aquifer system, and the Upper Cretaceous aquifer system in the Williston structural basin are less than about 750; 2,250; and 1,050 ft, respectively. Interpolated thicknesses of the lower Tertiary aquifer system and the Upper Cretaceous aquifer system in the Powder River structural basin are less than about 7,180 and 5,070 ft, respectively. Interpolated horizontal hydraulic conductivity values for the Williston structural basin were as much as 25 feet per day (ft/d) in the glacial deposits and had smaller ranges in the lower Tertiary aquifer system (0.01–9.8 ft/d) and in the Upper Cretaceous aquifer system (0.06–5.5 ft/d). In the Powder River structural basin, the lower Tertiary aquifer system had a greater range of interpolated horizontal hydraulic conductivity values (0.10–11 ft/d) than the Upper Cretaceous aquifer system (0.02–5.7 ft/d). Transmissivity is greatest in the gravel zones of the glacial deposits (2,120 feet squared per day) and generally decreases with depth into the bedrock units. Regionally, water in the lower Tertiary and Upper Cretaceous aquifer systems flows in a northerly or northeasterly direction from the Powder River structural basin to the Williston structural basin. Groundwater flow in the Williston structural basin generally is easterly or northeasterly. Flow in the uppermost hydrogeologic units generally is more local and controlled by topography where unglaciated in the Williston structural basin than is flow in the glaciated part and in underlying aquifers. Groundwater flow in the Powder River structural basin generally is northerly with local variations greatest in the uppermost aquifers. Groundwater is confined, and flow is regional in the underlying aquifers.

Powder River Basin;Williston Basin↗

Groundwater-quality data in the Cascade Range and Modoc Plateau study unit, 2010-Results from the California GAMA Program

Groundwater quality in the 39,000-square-kilometer Cascade Range and Modoc Plateau (CAMP) study unit was investigated by the U.S. Geological Survey (USGS) from July through October 2010, as part of the California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring and Assessment (GAMA) Program’s Priority Basin Project (PBP). The GAMA PBP was developed in response to the California Groundwater Quality Monitoring Act of 2001 and is being conducted in collaboration with the SWRCB and Lawrence Livermore National Laboratory (LLNL). The CAMP study unit is the thirty-second study unit to be sampled as part of the GAMA PBP. The GAMA CAMP study was designed to provide a spatially unbiased assessment of untreated-groundwater quality in the primary aquifer system and to facilitate statistically consistent comparisons of untreated-groundwater quality throughout California. The primary aquifer system is defined as that part of the aquifer corresponding to the open or screened intervals of wells listed in the California Department of Public Health (CDPH) database for the CAMP study unit. The quality of groundwater in shallow or deep water-bearing zones may differ from the quality of groundwater in the primary aquifer system; shallow groundwater may be more vulnerable to surficial contamination. In the CAMP study unit, groundwater samples were collected from 90 wells and springs in 6 study areas (Sacramento Valley Eastside, Honey Lake Valley, Cascade Range and Modoc Plateau Low Use Basins, Shasta Valley and Mount Shasta Volcanic Area, Quaternary Volcanic Areas, and Tertiary Volcanic Areas) in Butte, Lassen, Modoc, Plumas, Shasta, Siskiyou, and Tehama Counties. Wells and springs were selected by using a spatially distributed, randomized grid-based method to provide statistical representation of the study unit (grid wells). Groundwater samples were analyzed for field water-quality indicators, organic constituents, perchlorate, inorganic constituents, radioactive constituents, and microbial indicators. Naturally occurring isotopes and dissolved noble gases also were measured to provide a dataset that will be used to help interpret the sources and ages of the sampled groundwater in subsequent reports. In total, 221 constituents were investigated for this study. Three types of quality-control samples (blanks, replicates, and matrix spikes) were collected at approximately 10 percent of the wells in the CAMP study unit, and the results for these samples were used to evaluate the quality of the data for the groundwater samples. Blanks rarely contained detectable concentrations of any constituent, suggesting that contamination from sample collection procedures was not a significant source of bias in the data for the groundwater samples. Replicate samples generally were within the limits of acceptable analytical reproducibility. Matrix-spike recoveries were within the acceptable range (70 to 130 percent) for approximately 90 percent of the compounds. This study did not attempt to evaluate the quality of water delivered to consumers; after withdrawal from the ground, untreated groundwater typically is treated, disinfected, and (or) blended with other waters to maintain water quality. Regulatory benchmarks apply to water that is served to the consumer, not to untreated groundwater. However, to provide some context for the results, concentrations of constituents measured in the untreated groundwater were compared with regulatory and non-regulatory health-based benchmarks established by the U.S. Environmental Protection Agency (USEPA) and CDPH, and to non-regulatory benchmarks established for aesthetic concerns by CDPH. Comparisons between data collected for this study and benchmarks for drinking water are for illustrative purposes only and are not indicative of compliance or non-compliance with those benchmarks. All organic constituents and most inorganic constituents that were detected in groundwater samples from the 90 grid wells in the CAMP study unit were detected at concentrations less than drinking-water benchmarks. Of the 148 organic constituents analyzed, 27 were detected in groundwater samples; concentrations of all detected constituents were less than regulatory and nonregulatory health-based benchmarks, and all were less than 1/10 of benchmark levels. One or more organic constituents were detected in 52 percent of the grid wells in the CAMP study unit: VOCs were detected in 30 percent, and pesticides and pesticide degradates were detected in 31 percent. Trace elements, major ions, nutrients, and radioactive constituents were sampled for at 90 grid wells in the CAMP study unit, and most detected concentrations were less than health-based benchmarks. Exceptions include three detections of arsenic greater than the USEPA maximum contaminant level (MCL-US) of 10 micrograms per liter (µg/L), two detections of boron greater than the CDPH notification level (NL-CA) of 1,000 µg/L, two detections of molybdenum greater than the USEPA lifetime health advisory level (HAL-US) of 40 µg/L, two detections of vanadium greater than the CDPH notification level (NL-CA) of 50 µg/L, one detection of nitrate, as nitrogen, greater than the MCL-US of 10 milligrams per liter (mg/L), two detections of uranium greater than the MCL-US of 30 µg/L and the MCL-CA of 20 picocuries per liter (pCi/L), one detection of radon-222 greater than the proposed MCL-US of 4,000 pCi/L, and two detections of gross alpha particle activity greater than the MCL-US of 15 pCi/L. Results for inorganic constituents with non-regulatory benchmarks set for aesthetic concerns showed that iron concentrations greater than the CDPH secondary maximum contaminant level (SMCL-CA) of 300 µg/L were detected in four grid wells. Manganese concentrations greater than the SMCL-CA of 50 µg/L were detected in nine grid wells. Chloride and TDS were detected at concentrations greater than the upper SMCL-CA benchmarks of 500 mg/L and 1,000 mg/L, respectively, in one grid well. Microbial indicators (total coliform and Escherichia coli [E. coli]) were detected in 11 percent of the 83 grid wells sampled for these analyses in the CAMP study unit. The presence of total coliform was detected in nine grid wells, and the presence of E. coli was detected in one of these same grid wells.

California↗

Projected ground-water development, ground-water levels, and stream-aquifer leakage in the South Fork Solomon River Valley between Webster Reservoir and Waconda Lake, north-central Kansas, 1979-2020

A two-dimensional finite difference computer model was used to project changes in the potentiometric surface, saturated thickness, and stream aquifer leakage in an alluvial aquifer resulting from four instances of projected groundwater development. The alluvial aquifer occurs in the South Fork Solomon River valley between Webster Reservoir and Waconda Lake in north-central Kansas. In the first two projections, pumpage for irrigation was held constant at 1978 rates throughout the projection period (1979-2020). In the second two projections, the 1978 pumpage was progressively increased each yr through 2020. In the second and fourth projections, surface water diversions in the Osborne Irrigation Canal were decreased by 50 %. For the third and fourth projections, each grid-block in the modeled area was classified initially as one of six types according to whether it represented irrigable or nonirrigable land, to its saturated thickness, to its location inside or outside the canal-river area, and to its pumping rate. The projected base-flow rates (leakage from the aquifer to the river) were lower during the irrigation season (June, July, and August) than during the other months of the yr because of the decline in hydraulic head produced by groundwater pumpage. Stream depletion, calculated as a decrease below the average (1970-78) estimated winter base-flow rate of 16.5 cu ft/sec, varied inversely with base flow. For the first two projections, a constant annual cycle of well pumpage and recharge was used throughout the projection period. Aquifer leakage to the river was nearly constant by the mid-to-late 1990's, implying that flow conditions had attained a stabilized annual cycle. The third and fourth projections never attained an annual stabilized cycle because the irrigation pumpage rate was increased each year. By the early 1980's, the hydraulic head had fallen below river stage, reversing the hydraulic gradient at the stream-aquifer interface and resulting in net leakage from the river to the aquifer during the summer months. By the early 1990 's, the projected potentiometric surface of the aquifer was lower than the river stage even during the winter and spring months. (Author 's abstract)

Water-Resources Investigations Report↗

Status and understanding of groundwater quality in the Madera, Chowchilla Study Unit, 2008: California GAMA Priority Basin Project

Groundwater quality in the approximately 860-square-mile Madera and Chowchilla Subbasins (Madera-Chowchilla study unit) of the San Joaquin Valley Basin was investigated as part of the Priority Basin Project of the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is located in California's Central Valley region in parts of Madera, Merced, and Fresno Counties. The GAMA Priority Basin Project is being conducted by the California State Water Resources Control Board in collaboration with the U.S. Geological Survey (USGS) and the Lawrence Livermore National Laboratory. The Project was designed to provide statistically robust assessments of untreated groundwater quality within the primary aquifer systems in California. The primary aquifer system within each study unit is defined by the depth of the perforated or open intervals of the wells listed in the California Department of Public Health (CDPH) database of wells used for municipal and community drinking-water supply. The quality of groundwater in shallower or deeper water-bearing zones may differ from that in the primary aquifer system; shallower groundwater may be more vulnerable to contamination from the surface. The assessments for the Madera-Chowchilla study unit were based on water-quality and ancillary data collected by the USGS from 35 wells during April-May 2008 and water-quality data reported in the CDPH database. Two types of assessments were made: (1) status , assessment of the current quality of the groundwater resource, and (2) understanding , identification of natural factors and human activities affecting groundwater quality. The primary aquifer system is represented by the grid wells, of which 90 percent (%) had depths that ranged from about 200 to 800 feet (ft) below land surface and had depths to the top of perforations that ranged from about 140 to 400 ft below land surface. Relative-concentrations (sample concentrations divided by benchmark concentrations) were used for evaluating groundwater quality for those constituents that have Federal or California regulatory or non-regulatory benchmarks for drinking-water quality. A relative-concentration (RC) greater than 1.0 indicates a concentration above a benchmark. RCs for organic constituents (volatile organic compounds and pesticides) and special-interest constituents (perchlorate) were classified as "high" (RC is greater than 1.0), "moderate" (RC is less than or equal to 1.0 and greater than 0.1), or "low" (RC is less than or equal to 0.1). For inorganic constituents (major and minor ions, trace elements, nutrients, and radioactive constituents), the boundary between low and moderate RCs was set at 0.5. The assessments characterize untreated groundwater quality, not the quality of treated drinking water delivered to consumers by water purveyors; drinking-water benchmarks, and thus relative-concentrations, are used to provide context for the concentrations of constituents measured in groundwater. Aquifer-scale proportion was used in the status assessment as the primary metric for evaluating regional-scale groundwater quality. High aquifer-scale proportion is defined as the percentage of the area of the primary aquifer system with RCs greater than 1.0 for a particular constituent or class of constituents; moderate and low aquifer-scale proportions are defined as the percentages of the area of the primary aquifer system with moderate and low RCs, respectively. Percentages are based on an areal, rather than a volumetric basis. Two statistical approaches--grid-based, which used one value per grid cell, and spatially weighted, which used multiple values per grid cell--were used to calculate aquifer-scale proportions for individual constituents and classes of constituents. The spatially weighted estimates of high aquifer-scale proportions were within the 90% confidence intervals of the grid-based estimates for all constituents except iron. The status assessment showed that inorganic constituents had greater high and moderate aquifer-scale proportions in the Madera-Chowchilla study unit than did organic constituents. RCs for inorganic constituents with health-based benchmarks were high in 37% of the primary aquifer system, moderate in 30%, and low in 33%. The inorganic constituents contributing most to the high aquifer-scale proportion were arsenic (13%), uranium (17%), gross alpha particle activity (20%), nitrate (6.7%), and vanadium (3.3%). RCs for inorganic constituents with non-health-based benchmarks were high in 6.7% of the primary aquifer system, and the constituent contributing most to the high aquifer-scale proportion was total dissolved solids (TDS). RCs for organic constituents with health-based benchmarks were high in 10% of the primary aquifer system, moderate in 3.3%, and low in 40%; organic constituents were not detected in 47% of the primary aquifer system. The fumigant 1,2-dibromo-3-chloropropane (DBCP) was the only organic constituent detected at high RCs. Seven organic constituents were detected in 10% or more of the primary aquifer system: DBCP; the fumigant additive 1,2,3-trichloropropane; the herbicides simazine, atrazine, and diuron; the trihalomethane chloroform; and the solvent tetrachloroethene (PCE). RCs for the special-interest constituent perchlorate were moderate in 20% of the primary aquifer system. The second component of this study, the understanding assessment , identified the natural and human factors that may affect groundwater quality by evaluating statistical correlations between water-quality constituents and potential explanatory factors, such as land use, position relative to important geologic features, groundwater age, well depth, and geochemical conditions in the aquifer. Results of the statistical evaluations were used to explain the distribution of constituents in the study unit. Depth to the top of perforations in the well and groundwater age were the most important explanatory factors for many constituents. High and moderate RCs of nitrate, uranium, and TDS and the presence of herbicides, trihalomethanes, and solvents were all associated with depths to the top of perforations less than 235 ft and modern- and mixed-age groundwater. Positive correlations between uranium, bicarbonate, TDS, and the proportion of calcium and magnesium in the total cations suggest that downward movement of recharge from irrigation water contributed to the elevated concentrations of these constituents in the primary aquifer system. High and moderate RCs of arsenic were associated with depths to the top of perforations greater than 235 ft, mixed- and pre-modern-age groundwater, and location in sediments from the Chowchilla River alluvial fan, suggesting that increased residence time and appropriate aquifer materials were needed for arsenic to accumulate in the groundwater. High and moderate RCs of fumigants were associated with depths to the top of perforations of less than 235 ft and location south of the city of Madera; low RCs of fumigants were detected in wells dispersed across the study unit with a range of depths to top of perforations.

California↗

Simulation of groundwater withdrawal scenarios for the Redwall-Muav and Coconino Aquifer Systems of northern and central Arizona

The Northern Arizona Regional Groundwater Flow Model was used to estimate the hydrologic changes, including water-level change and groundwater discharge to streams and springs, that may result from future changes in groundwater withdrawals in and near the Coconino Plateau Water Advisory Council study area, Coconino and Navajo Counties, Arizona. Three future groundwater withdrawal scenarios for tribal and nontribal uses were developed by the Coconino Plateau Water Advisory Council and were simulated for the period representing the years from 2006 through 2105. Scenario 1 assumes no major changes in groundwater use except for increased demand based on population projections. Scenario 2 assumes that a pipeline will provide a source of surface water from Lake Powell to areas near Cameron and Moenkopi that would replace local groundwater withdrawals. Scenario 3 assumes that the pipeline extends to the Flagstaff and Williams areas, and would replace groundwater demands for water in the area. The Coconino Plateau Water Advisory Council withdrawal scenarios primarily influence water levels and groundwater discharge in the Coconino Plateau basin, near the western margin of the Little Colorado River Plateau basin, and the Verde Valley subbasin. Simulated effects of the withdrawal scenarios are superimposed on effects of previous variations in groundwater withdrawals and artificial and incidental recharge. Pre-scenario variations include changes in water-levels in wells; groundwater storage; discharge to streams and springs; and evapotranspiration by plants that use groundwater. Future variations in groundwater discharge and water-levels in wells will continue to occur as a result of both the past and any future changes. Water-level variations resulting from post-2005 stresses, including groundwater withdrawals and incidental and artificial recharge, in the area of the withdrawal scenarios are primarily localized and superimposed on the regional changes caused by variations in stresses that occurred since the beginning of the initial stresses in the early 1900s through 2005. Withdrawal scenario 1 produced a broad region on the Coconino Plateau where water-levels declined 3–5 feet by 2105, and local areas with water-level declines of 100 feet or more where groundwater withdrawals are concentrated, near the City of Flagstaff Woody Mountain and Lake Mary well fields, and the towns of Tusayan, Williams, and Moenkopi. Water-level rises of 100 feet or more were simulated at areas of incidental recharge near wastewater treatment facilities near Flagstaff, Tusayan, Grand Canyon South Rim, Williams, and Munds Park. Simulated water-level change from 2006 through 2105 for scenarios 2 and 3 is mostly different from water-level change simulated for scenario 1 at the local level. For scenarios 2 and 3, water levels near Cameron in 2105 where 1–3 feet higher than simulated for scenario 1. Water levels at Moenkopi are more than 100 feet higher due to the elimination of a proposed withdrawal well that was simulated in scenario 1. Scenario 3 eliminates more groundwater withdrawals in the Flagstaff and Williams areas, simulates 1–3 feet less water-level decline than scenario 1 across much of the Coconino Plateau, and water levels that are as much as 50 feet higher than simulated by scenario 1 near withdrawal wells in the Williams and Flagstaff areas. Scenario 1 simulated the most change in groundwater discharge for the Little Colorado River below Cameron and for Oak Creek above Page Springs where declines in discharge of about 1.3 and 0.9 cubic feet per second (ft 3 /s), respectively, were simulated. Other simulated changes in discharge through 2105 in scenario 1 are losses of less than 0.4 ft 3 /s at the Upper Verde River, losses of less than 0.3 ft 3 /s at Havasu Creek and at Colorado River below Havasu Creek, losses of less than 0.1 ft 3 /s at Clear Creek, and increases in flow at the south rim springs and Chevelon Creek of less than 0.1 and 0.3 ft 3 /s, respectively. Simulated changes in discharge for scenarios 2 and 3 are less than for scenario 1 because of lower rates of groundwater withdrawal. Scenario 3 resulted in greater groundwater discharge than scenarios 1 and 2 at all major groundwater discharge features from 2006 through 2105 except for Clear and Chevelon Creeks, where the same groundwater discharge was simulated by each of the three scenarios. Changes in groundwater discharge are expected to occur after 2105 to all major surface features that discharge from the Redwall-Muav and Coconino aquifers because change in aquifer storage was occurring at the end of the simulation in 2105. The accuracy of simulated changes resulting from the Coconino Plateau Water Advisory Council groundwater withdrawal scenarios is dependent on the persistence of several hydrologic assumptions that are inherent in the Northern Arizona Regional Groundwater Flow Model including, but not limited to, the reasonably accurate simulation of (1) transmissivity distributions, (2) distributions of vertical hydraulic properties, (3) distributions of spatial rates of withdrawal and incidental recharge, (4) aquifer extents, and (5) hydrologic barriers and conduits.

Arizona↗

Recent advances in simulation-optimization groundwater management modeling

In recent years, the aquifer simulation model has been combined with techniques of optimization to address important groundwater management problems. The combined simulation and optimization model accounts for the complex behavior of the groundwater system and identifies the best management strategy under consideration of the management objectives and constraints. Simulation-optimization groundwater management models have been developed for a variety of applications, such as restoration of contaminated groundwater, control of aquifer hydraulics, allocation of ground and surface water resources, and evaluation of groundwater policies [see reviews by Gorelick , 1983, 1990; Yeh , 1992; Ahlfeld and Heidari , 1994; Bredehoeft et al , 1994]. Gorelick [1983] divides groundwater management models into three categories: groundwater hydraulic management, groundwater quality management, and groundwater policy evaluation and allocation. This review focuses in detail on the recent advances made in groundwater hydraulic and groundwater quality management modeling, which encompasses most of the aquifer management work done in the past four years. For a review of recent research in groundwater policy evaluation and allocation, see Bredehoeft et al . [1994].

Reviews of Geophysics↗

Strontium isotope geochemistry of groundwater in the central part of the Dakota (Great Plains) aquifer, USA

The Dakota aquifer of the central and eastern Great Plains of the United States is an important source of water for municipal supplies, irrigation and industrial use. Although the regional flow system can be characterized generally as east to northeasterly from the Rocky Mountains towards the Missouri River, locally the flow systems are hydrologically complex. This study uses Sr isotopic data from groundwater and leached aquifer samples to document the complex subsystems within the Dakota aquifer in Nebraska and Kansas. The interaction of groundwater with the geologic material through which it flows has created spatial patterns in the isotopic measurements that are related to: long-term water-rock interaction, during which varying degrees of isotopic equilibrium between water and rock has been achieved; and the alteration of NaCl fluids by water-rock interaction. Specifically, Sr isotopic data distinguish brines from Kansas and western Nebraska from those in eastern Nebraska: the former are interpreted to reflect interaction with Permian rocks, whereas the latter record interaction with Pennsylvanian rocks. The Sr isotopic composition of groundwater from other parts of Nebraska and Kansas are a function of the dynamic interaction between groundwater and unlithified sediments (e.g., glacial till and loess), followed by interaction with oxidized and unoxidized sediments within the Dakota Formation. This study illustrates the power of combining Sr chemistry with more conventional geochemical data to obtain a more complete understanding of groundwater flow systems within regional aquifer systems where extensive monitoring networks do not exist. ?? 2003 Elsevier Ltd. All rights reserved.

Applied Geochemistry↗

Groundwater quality in the Western San Joaquin Valley study unit, 2010: California GAMA Priority Basin Project

Water quality in groundwater resources used for public drinking-water supply in the Western San Joaquin Valley (WSJV) was investigated by the USGS in cooperation with the California State Water Resources Control Board (SWRCB) as part of its Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. The WSJV includes two study areas: the Delta–Mendota and Westside subbasins of the San Joaquin Valley groundwater basin. Study objectives for the WSJV study unit included two assessment types: (1) a status assessment yielding quantitative estimates of the current (2010) status of groundwater quality in the groundwater resources used for public drinking water, and (2) an evaluation of natural and anthropogenic factors that could be affecting the groundwater quality. The assessments characterized the quality of untreated groundwater, not the quality of treated drinking water delivered to consumers by water distributors. The status assessment was based on data collected from 43 wells sampled by the U.S. Geological Survey for the GAMA Priority Basin Project (USGS-GAMA) in 2010 and data compiled in the SWRCB Division of Drinking Water (SWRCB-DDW) database for 74 additional public-supply wells sampled for regulatory compliance purposes between 2007 and 2010. To provide context, concentrations of constituents measured in groundwater were compared to U.S. Environmental Protection Agency (EPA) and SWRCB-DDW regulatory and non-regulatory benchmarks for drinking-water quality. The status assessment used a spatially weighted, grid-based method to estimate the proportion of the groundwater resources used for public drinking water that has concentrations for particular constituents or class of constituents approaching or above benchmark concentrations. This method provides statistically unbiased results at the study-area scale within the WSJV study unit, and permits comparison of the two study areas to other areas assessed by the GAMA Priority Basin Project statewide. Groundwater resources used for public drinking water in the WSJV study unit are among the most saline and most affected by high concentrations of inorganic constituents of all groundwater resources used for public drinking water that have been assessed by the GAMA Priority Basin Project statewide. Among the 82 GAMA Priority Basin Project study areas statewide, the Delta–Mendota study area ranked above the 90th percentile for aquifer-scale proportions of groundwater resources having concentrations of total dissolved solids (TDS), sulfate, chloride, manganese, boron, chromium(VI), selenium, and strontium above benchmarks, and the Westside study area ranked above the 90th percentile for TDS, sulfate, manganese, and boron. In the WSJV study unit as a whole, one or more inorganic constituents with regulatory or non-regulatory, health-based benchmarks were present at concentrations above benchmarks in about 53 percent of the groundwater resources used for public drinking water, and one or more organic constituents with regulatory health-based benchmarks were detected at concentrations above benchmarks in about 3 percent of the resource. Individual constituents present at concentrations greater than health-based benchmarks in greater than 2 percent of groundwater resources used for public drinking water included: boron (51 percent, SWRCB-DDW notification level), chromium(VI) (25 percent, SWRCB-DDW maximum contaminant level (MCL)), arsenic (10 percent, EPA MCL), strontium (5.1 percent, EPA Lifetime health advisory level (HAL)), nitrate (3.9 percent, EPA MCL), molybdenum (3.8 percent, EPA HAL), selenium (2.6 percent, EPA MCL), and benzene (2.6 percent, SWRCB-DDW MCL). In addition, 50 percent of the resource had TDS concentrations greater than non-regulatory, aesthetic-based SWRCB-DDW upper secondary maximum contaminant level (SMCL), and 44 percent had manganese concentrations greater than the SWRCB-DDW SMCL. Natural and anthropogenic factors that could affect the groundwater quality were evaluated by using results from statistical testing of associations between constituent concentrations and values of potential explanatory factors, inferences from geochemical and age-dating tracer results, and by considering the water-quality results in the context of the hydrogeologic setting of the WSJV study unit. Natural factors, particularly the lithologies of the source areas for groundwater recharge and of the aquifers, were the dominant factors affecting groundwater quality in most of the WSJV study unit. However, where groundwater resources used for public supply included groundwater recharged in the modern era, mobilization of constituents by recharge of water used for irrigation also affected groundwater quality. Public-supply wells in the Westside study area had a median depth of 305 m and primarily tapped groundwater recharged hundreds to thousands of years ago, whereas public-supply wells in the Delta–Mendota study area had a median depth of 85 m and primarily tapped either groundwater recharged within the last 60 years or groundwater consisting of mixtures of this modern recharge and older recharge. Public-supply wells in the WSJV study unit are screened in the Tulare Formation and zones above and below the Corcoran Clay Member are used. The Tulare Formation primarily consists of alluvial sediments derived from the Coast Ranges to the west, except along the valley trough at the eastern margin of the WSJV study unit where the Tulare Formation consists of fluvial sands derived from the Sierra Nevada to the east. Groundwater from wells screened in the Sierra Nevada sands had manganese-reducing or manganese- and iron-reducing oxidation-reduction (redox) conditions. These redox conditions commonly were associated with elevated arsenic or molybdenum concentrations, and the dominance of arsenic(III) in the dissolved arsenic supports reductive dissolution of iron and manganese oxyhydroxides as the mechanism. In addition, groundwater from many wells screened in Sierra Nevada sands contained low concentrations of nitrite or ammonium, indicating reduction of nitrate by denitrification or dissimilatory processes, respectively. Geology of the Coast Ranges westward of the study unit strongly affects groundwater quality in the WSJV. Elevated concentrations of TDS, sulfate, boron, selenium and strontium in groundwater were primarily associated with aquifer sediments and recharge derived from areas of the Coast Ranges dominated by Cretaceous-to-Miocene age, organic-rich, reduced marine shales, known as the source of selenium in WSJV soils, surface water, and groundwater. Low sulfur-isotopic values (δ34S) of dissolved sulfate indicate that the sulfate was largely derived from oxidation of biogenic pyrite from the shales, and correlations with trace element concentrations, geologic setting, and groundwater geochemical modeling indicated that distributions of sulfate, strontium, and selenium in groundwater were controlled by dissolution of secondary sulfate minerals in soils and sediments. Elevated concentrations of chromium(VI) were primarily associated with aquifer sediments and recharge derived from areas of the Coast Ranges dominated by the Franciscan Complex and ultramafic rocks. The Franciscan Complex also has boron-rich, sodium-chloride dominated hydrothermal fluids that contribute to elevated concentrations of boron and TDS. Groundwater from wells screened in Coast Ranges alluvium was primarily oxic and relatively alkaline (median pH value of 7.55) in the Delta–Mendota study area, and primarily nitrate-reducing or suboxic and alkaline (median pH value of 8.4) in the Westside study area. Many groundwater samples from those wells have elevated concentrations of arsenic(V), molybdenum, selenium, or chromium(VI), consistent with desorption of metal oxyanions from mineral surfaces under those geochemical conditions. High concentrations of benzene were associated with deep wells located in the vicinity of petroleum deposits at the southern end of the Westside study area. Groundwater from these wells had premodern age and anoxic geochemical conditions, and the ratios among concentrations of hydrocarbon constituents were different from ratios found in fuels and combustion products, which is consistent with a geogenic source for the benzene rather than contamination from anthropogenic sources. Water stable-isotope compositions, groundwater recharge temperatures, and groundwater ages were used to infer four types of groundwater: (1) groundwater derived from natural recharge of water from major rivers draining the Sierra Nevada; (2) groundwater primarily derived from natural recharge of water from Coast Ranges runoff; (3) groundwater derived from recharge of pumped groundwater applied to the land surface for irrigation; and (4) groundwater derived from recharge during a period of much cooler paleoclimate. Water previously used for irrigation was found both above and below the Corcoran Clay, supporting earlier inferences that this clay member is no longer a robust confining unit. Recharge of water used for irrigation has direct and indirect effects on groundwater quality. Elevated nitrate concentrations and detections of herbicides and fumigants in the Delta–Mendota study area generally were associated with greater agricultural land use near the well and with water recharged during the last 60 years. However, the extent of the groundwater resource affected by agricultural sources of nitrate was limited by groundwater redox conditions sufficient to reduce nitrate. The detection frequency of perchlorate in Delta–Mendota groundwater was greater than expected for natural conditions. Perchlorate, nitrate, selenium, and strontium concentrations were correlated with one another and were greater in groundwater inferred to be recharge of previously pumped groundwater used for irrigation. The source of the perchlorate, selenium, and strontium appears to be salts deposited in the soils and sediments of the arid WSJV that are dissolved and flushed into groundwater by the increased amount of recharge caused by irrigation. In the Delta–Mendota study area, the groundwater with elevated concentrations of selenium was found deeper in the aquifer system than it was reported by a previous study 25 years earlier, suggesting that this transient front of groundwater with elevated concentrations of constituents derived from dissolution of soil salts by irrigation recharge is moving down through the aquifer system and is now reaching the depth zone used for public drinking water supply.

California↗

Groundwater availability of the Columbia Plateau Regional Aquifer System, Washington, Oregon, and Idaho

The Columbia Plateau Regional Aquifer System (CPRAS) covers about 44,000 square miles of southeastern Washington, northeastern Oregon, and western Idaho. The area supports a $6-billion per year agricultural industry, leading the Nation in production of apples, hops, and eight other commodities. Groundwater pumpage and surface-water diversions supply water to croplands that account for about 5 percent of the Nation’s irrigated lands. Groundwater also is the primary source of drinking water for the more than 1.3 million people in the study area. Increasing competitive demands for water for municipal, fisheries/ecosystems, agricultural, domestic, hydropower, and recreational uses must be met by additional groundwater withdrawals and (or) by changes in the way water resources are allocated and used throughout the hydrologic system. As of 2014, most surface-water resources in the study area were either over allocated or fully appropriated, especially during the dry summer season. In response to continued competition for water, numerous water-management activities and concerns have gained prominence: water conservation, conjunctive use, artificial recharge, hydrologic implications of land-use change, pumpage effects on streamflow, and effects of climate variability and change. An integrated understanding of the hydrologic system is important in order to implement effective water-resource management strategies that address these concerns. To provide information to stakeholders involved in water-management activities, the U.S. Geological Survey (USGS) Groundwater Resources Program assessed the groundwater availability as part of a national study of regional systems (U.S. Geological Survey, 2008). The CPRAS assessment includes: The present status of groundwater resources, How these resources have changed over time, and Development and application of tools to estimate system responses to stresses from future uses and climate variability and change. This effort builds on previous investigations, especially the USGS Columbia Plateau Regional Aquifer-System Analysis study (CP-RASA). A major product of this new assessment is a numerical groundwater-flow model of the system. The model was used to estimate water-budget components of the hydrogeologic units composing the groundwater system, and to evaluate groundwater availability under existing land- and water-use conditions and a possible future climate scenario representing an increase in pumpage demand due to a warming climate. Information from this study also allowed for analysis of: The CPRAS for predevelopment times (pre-1920), Variations from 1920 through 2007, Conditions during 1985–2007 (referred to as “existing conditions”), and Changes in the system from predevelopment times. The model also is a useful tool for investigating water supply, water demand, management strategies, groundwater-surface water exchanges, and potential effects of changing climate on the hydrologic system.

Idaho, Oregon, Washington↗

Quality-assurance plan for groundwater activities, U.S. Geological Survey, Washington Water Science Center

This report documents the standard procedures, policies, and field methods used by the U.S. Geological Survey’s (USGS) Washington Water Science Center staff for activities related to the collection, processing, analysis, storage, and publication of groundwater data. This groundwater quality-assurance plan changes through time to accommodate new methods and requirements developed by the Washington Water Science Center and the USGS Office of Groundwater. The plan is based largely on requirements and guidelines provided by the USGS Office of Groundwater, or the USGS Water Mission Area. Regular updates to this plan represent an integral part of the quality-assurance process. Because numerous policy memoranda have been issued by the Office of Groundwater since the previous groundwater quality assurance plan was written, this report is a substantial revision of the previous report, supplants it, and contains significant additional policies not covered in the previous report. This updated plan includes information related to the organization and responsibilities of USGS Washington Water Science Center staff, training, safety, project proposal development, project review procedures, data collection activities, data processing activities, report review procedures, and archiving of field data and interpretative information pertaining to groundwater flow models, borehole aquifer tests, and aquifer tests. Important updates from the previous groundwater quality assurance plan include: (1) procedures for documenting and archiving of groundwater flow models; (2) revisions to procedures and policies for the creation of sites in the Groundwater Site Inventory database; (3) adoption of new water-level forms to be used within the USGS Washington Water Science Center; (4) procedures for future creation of borehole geophysics, surface geophysics, and aquifer-test archives; and (5) use of the USGS Multi Optional Network Key Entry System software for entry of routine water-level data collected as part of long-term water-level monitoring networks.

Open-File Report↗

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer, a major aquifer in the Pecos County region of western Texas, is a vital groundwater resource for agricultural, industrial, and public supply uses. Resource managers would like to better understand the future availability of water in the Edwards-Trinity aquifer in the Pecos County region and the effects of the possible increase or temporal redistribution of groundwater withdrawals. To that end, the U.S. Geological Survey (USGS), in cooperation with the Middle Pecos Groundwater Conservation District, Pecos County, City of Fort Stockton, Brewster County, and Pecos County Water Control and Improvement District No. 1, completed a comprehensive, integrated analysis of available hydrogeologic data to develop a groundwater-flow model of the Edwards-Trinity and related aquifers in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. Following calibration, the model was used to evaluate the sustainability of recent (2008) and projected water-use demands on groundwater resources in the study area.

Texas↗

Improving groundwater model calibration with repeat microgravity measurements

Groundwater-flow models depend on hydraulic head and flux observations for evaluation and calibration. A different type of observation—change in storage measured using repeat microgravity—can also be used for parameter estimation by simulating the expected change in gravity from a groundwater model and including the observation misfit in the objective function. The method is demonstrated using new software linked to MODFLOW input and output files and field data from the vicinity of the All American Canal in southeast California, USA. Over a 10-year period following lining of the previously highly permeable canal with concrete, gravity decreased by over 100 μGal (equivalent to about 2.5 m of free-standing water) at some locations as seepage decreased and the remnant groundwater mound dissipated into the aquifer or was removed by groundwater pumping. Simulated gravity from a MODFLOW model closely matched observations, and repeat microgravity data proved useful for constraining both hydraulic conductivity and specific yield estimates. Specific yield estimated using the infinite-horizontal slab approximation agreed well with model-derived values, and the departure from the linear, flat-water-table approximation was small, less than 2%, despite relatively large and dynamic water-table slope. First-order second-moment parameter uncertainty analysis shows reduction in uncertainty for all hydraulic conductivity and specific yield parameter estimates with the addition of repeat microgravity data, as compared to drawdown data alone.

California↗

Simulation of groundwater flow in the Brunswick Area, Georgia, for 2004 and 2015, and selected groundwater-management scenarios

The Upper Floridan aquifer (UFA) is the principal water source for industrial and public supply in Glynn County, Georgia. Wells in active pumping centers that tap the UFA for industries near the city of Brunswick have created an upward hydraulic-head gradient in the Floridan aquifer system, which has allowed high chloride (saline) groundwater from the Fernandina permeable zone of the Lower Floridan aquifer (LFA) to migrate upward into freshwater zones. Chloride concentrations of more than 250 milligrams per liter—the State and Federal secondary drinking-water standard—have been measured in a 2-square-mile area near downtown Brunswick. An existing regional U.S. Geological Survey modular finite-difference groundwater-flow model (MODFLOW-2000) was modified using greater horizontal and vertical resolution to enable more detailed simulation of the effects of pumping in the vicinity of chloride contamination. Modifications to the regional model consisted of (1) limiting grid size to a maximum of 500 feet (ft) per side in the vicinity of the chloride plume; (2) representing the upper and lower Brunswick aquifers with distinct model layers; (3) similarly, representing upper and lower water-bearing zones of the UFA with distinct model layers in Glynn and Camden Counties, Ga.; and (4) establishing new hydraulic-property geographic zones in the UFA within Glynn County. The revised groundwater-flow model was calibrated to steady-state conditions that were assumed to exist during 2000 and 2004. The calibration and framework of the revised groundwater-flow model were documented in a separate report. For the current study, steady-state conditions were calibrated using October 2015 pumping rates in the Brunswick/Glynn County area as a 2015 Base Case. The 2015 Base Case simulation was used as the basis to evaluate seven groundwater-management scenarios in the Brunswick/Glynn County area. Seven groundwater management-scenarios were developed on the basis of short- and long-term groundwater-use projections for the UFA in the Brunswick/Glynn County area. Scenarios A and B simulated additional pumping in the upper water-bearing zone (UWBZ) of the UFA at existing public-supply wells located near a chloride plume and planned public-supply wells to be constructed north of downtown Brunswick. Scenario C simulated a shutdown at Brunswick Cellulose Inc. and Pinova Inc. and the resulting deactivation of nine production wells, with a combined total pumping of 31.3 million gallons per day (Mgal/d) for the 2015 Base Case simulation. Scenario D (three scenarios) simulated 12.5, 25, and 50 percent (designated Scenarios D1, D2, and D3) of the total pumping of 31.3 Mgal/d at Brunswick Cellulose and Pinova. The objective of Scenario D was to determine pumping rates that may reverse groundwater-flow directions toward the Brunswick Cellulose well field and potentially allow groundwater with higher chloride concentration to migrate toward nearby public-supply wells. Scenario E simulated an additional pumping of 5 Mgal/d from the UWBZ of the UFA at a recently constructed production well within the Brunswick Cellulose well field. Backward particle-tracking (MODPATH) analysis in public-supply wells located just outside the chloride plume to the north shows that predominant groundwater-flow directions are from the northeast toward the Brunswick Cellulose well field. The analysis covered 20- and 50-year periods for the 2015 Base Case and Scenario C simulations with 100 percent of backtracked particles remaining in the UWBZ and lower water-bearing zone of the UFA. Groundwater-flow directions are characterized by some vertical movement and dominant horizontal movement away from the chloride plume in the northern Brunswick area. For the 2015 Base Case simulation, the mean rate of particle movement ranged from 268 to 413 feet per year. For the Scenario C simulation, the mean rate of particle movement ranged from 89 to 182 feet per year with 50 percent of particles migrating from the chloride plume area. The rate of particle movement is influenced most by the horizontal hydraulic-head gradient in the UWBZ of the UFA. The revised groundwater-flow model is subject to the limitations documented in the original model. In addition, the values used for the specified-head boundaries in the Floridan aquifer system for the 2004 calibrated model were based on the sparse data available and were not changed for the 2015 update to the model. These model boundaries control 80 percent of the inflows and about 60 percent of the outflows. Composite-scaled sensitivities of the model parameters indicate the revised model is most sensitive to pumping rates, followed by the horizontal hydraulic conductivity in the UFA for zones along coastal Georgia.

Georgia↗

A comparison of recharge rates in aquifers of the United States based on groundwater-age data

An overview is presented of existing groundwater-age data and their implications for assessing rates and timescales of recharge in selected unconfined aquifer systems of the United States. Apparent age distributions in aquifers determined from chlorofluorocarbon, sulfur hexafluoride, tritium/helium-3, and radiocarbon measurements from 565 wells in 45 networks were used to calculate groundwater recharge rates. Timescales of recharge were defined by 1,873 distributed tritium measurements and 102 radiocarbon measurements from 27 well networks. Recharge rates ranged from < 10 to 1,200 mm/yr in selected aquifers on the basis of measured vertical age distributions and assuming exponential age gradients. On a regional basis, recharge rates based on tracers of young groundwater exhibited a significant inverse correlation with mean annual air temperature and a significant positive correlation with mean annual precipitation. Comparison of recharge derived from groundwater ages with recharge derived from stream base-flow evaluation showed similar overall patterns but substantial local differences. Results from this compilation demonstrate that age-based recharge estimates can provide useful insights into spatial and temporal variability in recharge at a national scale and factors controlling that variability. Local age-based recharge estimates provide empirical data and process information that are needed for testing and improving more spatially complete model-based methods.

Hydrogeology Journal↗

A parametric numerical analysis of factors controlling ground ruptures caused by groundwater pumping

A modeling analysis is used to investigate the relative susceptibility of various hydrogeologic configurations to aseismic rupture generation due to deformation of aquifer systems accompanying groundwater pumping. An advanced numerical model (GEPS3D) is used to simulate rupture generation and propagation for three typical processes: (i) reactivation of a preexisting fault, (ii) differential compaction due to variations in thickness of aquifer/aquitard layers constituting the aquifer system, and (iii) tensile fracturing above a bedrock ridge that forms the base of the aquifer system. A sensitivity analysis is developed to address the relative importance of various factors, including aquifer depletion, aquifer thickness, the possible uneven distribution and depth below land surface of the aquifer/aquitard layers susceptible to aquifer-system compaction, and the height of bedrock ridges beneath the aquifer system which contributes to thinning of the aquifer system. The rupture evolution is classified in two occurrences. In one, the rupture develops at the top of the aquifer or at land surface and does not propagate. In the other, the developed rupture propagates from the aquifer top toward the land surface and/or from the land surface downward. The aquifer depth is the most important factor controlling rupture evolution. Specifically, the probability of a significant rupture propagation is higher when the aquifer top is near land surface. The numerical results are processed by a statistical regression analysis to provide a general methodology for a preliminary evaluation of possible ruptures development in exploited aquifer systems susceptible to aquifer-system compaction and accompanying land subsidence. A comparison with a few representative case studies in Arizona, USA, China, and Mexico supports the study outcomes.

Arizona↗

Evaluation of groundwater quality and selected hydrologic conditions in the South Coast aquifer, Santa Isabel area, Puerto Rico, 2008–09

The source of drinking water in the Santa Isabel and Coamo areas of Puerto Rico (Molina and Gómez-Gómez, 2008) is the South Coast aquifer (hereafter referred to as the aquifer), which supplies about 30,700 cubic meters per day (m³/d) to Puerto Rico Aqueduct and Sewer Authority (PRASA) public-supply wells. In addition, approximately 45 wells provide an estimated 33,700 m³/d of groundwater to irrigate crops in the area. In 1967, baseline nitrate concentrations in groundwater throughout most of the aquifer were generally less than 6 milligrams per liter (mg/L) as nitrogen in collected water samples (U.S. Geological Survey, 2012). In 2007, elevated nitrate concentrations were detected in the aquifer, near Santa Isabel and the foothills north of the coastal plain at Santa Isabel as part of a regional groundwater-quality assessment conducted by the U.S. Geological Survey (USGS) during 2007 (Rodríguez and Gómez-Gómez, 2008). The increase in nitrate concentrations has been of concern to local government agencies because of its potential effect on public supply. To address public-supply concerns, the USGS, in cooperation with the Puerto Rico Department of Natural and Environmental Resources (PRDNER), evaluated groundwater quality in the aquifer near the Santa Isabel area between January 2008 and May 2009. The objectives of the study were to (1) define the groundwater-quality conditions of the aquifer, with emphasis on the distribution of nitrate concentrations; (2) identify potential sources leading to elevated nitrate concentrations; (3) estimate the nitrate loads from major sources identified; and (4) estimate the groundwater withdrawals by principal-use categories in the area. Results of this study will be used by Commonwealth of Puerto Rico and Federal agencies in developing strategies that can result in containment of high nitrate groundwater to minimize degradation of fresh groundwater in the aquifer.

Santa Isabel↗