Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

Follow this topic with RSS

At least 73 records · Page 4Linked to original sources

USGS2013

Simulation of groundwater flow, effects of artificial recharge, and storage volume changes in the Equus Beds aquifer near the city of Wichita, Kansas well field, 1935–2008

The Equus Beds aquifer is a primary water-supply source for Wichita, Kansas and the surrounding area because of shallow depth to water, large saturated thickness, and generally good water quality. Substantial water-level declines in the Equus Beds aquifer have resulted from pumping groundwater for agricultural and municipal needs, as well as periodic drought conditions. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project to store and later recover groundwater, and to form a hydraulic barrier to the known chloride-brine plume near Burrton, Kansas. In October 2009, the U.S. Geological Survey, in cooperation with the city of Wichita, began a study to determine groundwater flow in the area of the Wichita well field, and chloride transport from the Arkansas River and Burrton oilfield to the Wichita well field. Groundwater flow was simulated for the Equus Beds aquifer using the three-dimensional finite-difference groundwater-flow model MODFLOW-2000. The model simulates steady-state and transient conditions. The groundwater-flow model was calibrated by adjusting model input data and model geometry until model results matched field observations within an acceptable level of accuracy. The root mean square (RMS) error for water-level observations for the steady-state calibration simulation is 9.82 feet. The ratio of the RMS error to the total head loss in the model area is 0.049 and the mean error for water-level observations is 3.86 feet. The difference between flow into the model and flow out of the model across all model boundaries is -0.08 percent of total flow for the steady-state calibration. The RMS error for water-level observations for the transient calibration simulation is 2.48 feet, the ratio of the RMS error to the total head loss in the model area is 0.0124, and the mean error for water-level observations is 0.03 feet. The RMS error calculated for observed and simulated base flow gains or losses for the Arkansas River for the transient simulation is 7,916,564 cubic feet per day (91.6 cubic feet per second) and the RMS error divided by (/) the total range in streamflow (7,916,564/37,461,669 cubic feet per day) is 22 percent. The RMS error calculated for observed and simulated streamflow gains or losses for the Little Arkansas River for the transient simulation is 5,610,089 cubic feet per day(64.9 cubic feet per second) and the RMS error divided by the total range in streamflow (5,612,918/41,791,091 cubic feet per day) is 13 percent. The mean error between observed and simulated base flow gains or losses was 29,999 cubic feet per day (0.34 cubic feet per second) for the Arkansas River and -1,369,250 cubic feet per day (-15.8 cubic feet per second) for the Little Arkansas River. Cumulative streamflow gain and loss observations are similar to the cumulative simulated equivalents. Average percent mass balance difference for individual stress periods ranged from -0.46 to 0.51 percent. The cumulative mass balance for the transient calibration was 0.01 percent. Composite scaled sensitivities indicate the simulations are most sensitive to parameters with a large areal distribution. For the steady-state calibration, these parameters include recharge, hydraulic conductivity, and vertical conductance. For the transient simulation, these parameters include evapotranspiration, recharge, and hydraulic conductivity. The ability of the calibrated model to account for the additional groundwater recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project was assessed by using the U.S. Geological Survey subregional water budget program ZONEBUDGET and comparing those results to metered recharge for 2007 and 2008 and previous estimates of artificial recharge. The change in storage between simulations is the volume of water that estimates the recharge credit for the aquifer storage and recovery system. The estimated increase in storage of 1,607 acre-ft in the basin storage area compared to metered recharge of 1,796 acre-ft indicates some loss of metered recharge. Increased storage outside of the basin storage area of 183 acre-ft accounts for all but 6 acre-ft or 0.33 percent of the total. Previously estimated recharge credits for 2007 and 2008 are 1,018 and 600 acre-ft, respectively, and a total estimated recharge credit of 1,618 acre-ft. Storage changes calculated for this study are 4.42 percent less for 2007 and 5.67 percent more for 2008 than previous estimates. Total storage change for 2007 and 2008 is 0.68 percent less than previous estimates. The small difference between the increase in storage from artificial recharge estimated with the groundwater-flow model and metered recharge indicates the groundwater model correctly accounts for the additional water recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project. Small percent differences between inflows and outflows for all stress periods and all index cells in the basin storage area, improved calibration compared to the previous model, and a reasonable match between simulated and measured long-term base flow indicates the groundwater model accurately simulates groundwater flow in the study area. The change in groundwater level through recent years compared to the August 1940 groundwater level map has been documented and used to assess the change of storage volume of the Equus Beds aquifer in and near the Wichita well field for three different areas. Two methods were used to estimate changes in storage from simulation results using simulated change in groundwater levels in layer 1 between stress periods, and using ZONEBUDGET to calculate the change in storage in the same way the effects of artificial recharge were estimated within the basin storage area. The three methods indicate similar trends although the magnitude of storage changes differ. Information about the change in storage in response to hydrologic stresses is important for managing groundwater resources in the study area. The comparison between the three methods indicates similar storage change trends are estimated and each could be used to determine relative increases or decreases in storage. Use of groundwater level changes that do not include storage changes that occur in confined or semi-confined parts of the aquifer will slightly underestimate storage changes; however, use of specific yield and groundwater level changes to estimate storage change in confined or semi-confined parts of the aquifer will overestimate storage changes. Using only changes in shallow groundwater levels would provide more accurate storage change estimates for the measured groundwater levels method. The value used for specific yield is also an important consideration when estimating storage. For the Equus Beds aquifer the reported specific yield ranges between 0.08 and 0.35 and the storage coefficient (for confined conditions) ranges between 0.0004 and 0.16. Considering the importance of the value of specific yield and storage coefficient to estimates of storage change over time, and the wide range and substantial overlap for the reported values for specific yield and storage coefficient in the study area, further information on the distribution of specific yield and storage coefficient within the Equus Beds aquifer in the study area would greatly enhance the accuracy of estimated storage changes using both simulated groundwater level, simulated groundwater budget, or measured groundwater level methods.

Kansas↗

USGS2020

Water Quality of groundwater used for public supply in principal aquifers of the western United States

Groundwater provides nearly half of the Nation’s drinking water. As the Nation’s population grows, the importance of (and need for) high-quality drinking-water supplies increases. As part of a national-scale effort to assess groundwater quality in principal aquifers (PAs) that supply most of the groundwater used for public supply, the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Project staff sampled six principal aquifers in the western United States between 2013 and 2017: (1) the Basin and Range carbonate-rock aquifers, (2) Basin and Range basin-fill aquifers, (3) Rio Grande aquifer system, (4) High Plains aquifer, (5) Colorado Plateaus aquifers, and (6) Columbia Plateau basaltic-rock aquifers. These six PAs supply a large part of the Nation’s drinking water and cover a large geographic extent of the western conterminous United States. Groundwater samples were analyzed for a large suite of water-quality constituents including major ions, nutrients, trace elements, volatile organic compounds (VOCs), pesticide compounds, radioactive constituents, age tracers, and, in selected PAs, perchlorate. Two types of assessments were made: (1) a status assessment that describes the quality of the groundwater resource at time of collection and (2) an understanding assessment that evaluates relations between groundwater quality and potential explanatory factors that represent characteristics of the aquifer system. The assessments characterize untreated groundwater quality, which might be different than the quality of drinking water delivered to consumers. The assessments are based on water-quality data collected from 352 wells and 6 springs using an equal-area grid sampling design. This sampling approach allows for the estimation of the proportion of high, moderate, or low concentrations relative to federal water-quality benchmarks of selected constituents in the area of each PA. Results were compared to established benchmarks for drinking-water quality to provide context for evaluating the quality of untreated groundwater: Federal regulatory benchmarks for protecting human health, non-regulatory human-health benchmarks, and non-regulatory benchmarks for nuisance chemicals. Not all constituents that were analyzed have benchmarks and thus were not considered for assessments. Concentrations are characterized as high if they are greater than their benchmark. Concentrations are considered moderate if they are greater than one-half their benchmark (for inorganic constituents), or greater than one-tenth their benchmark (for organic constituents). Concentrations are considered low if they are less than moderate or the constituent was not detected. Status assessment results indicated that inorganic constituents more commonly occurred at high and moderate concentrations in the six PAs than organic constituents, and organic constituents predominately occurred at low concentrations. Inorganic constituents that exceeded health-based benchmarks (high concentrations) were present in all six PAs; aquifer-scale proportion were 30 percent in the Rio Grande aquifer system, 22 percent in the Basin and Range basin-fill aquifers, 20 percent in the Basin and Range carbonate-rock aquifers, 19 percent in the High Plains aquifer, 16 percent in the Colorado Plateaus aquifers, and 8 percent in the Columbia Plateau basaltic-rock aquifers. Arsenic, fluoride, manganese, and total dissolved solids were the constituents most commonly present at high concentrations. Organic constituents with human-health benchmarks (pesticide compounds and VOCs) did not occur at high concentrations and moderate concentrations were infrequent; aquifer-scale proportions ranged from 0 to 5 percent. Detections of organic compounds at low concentrations, however, occurred in all six PAs, with detection frequencies ranging from 10 to 26 percent for pesticide compounds and from 10 to 46 percent for VOCs. Specific organic constituents with detection frequencies greater than 10 percent were four herbicides (atrazine, didealkylatrazine, bromoform, and propazine), one insecticide (propoxur), and two VOCs (the trihalomethanes chloroform and bromodichloromethane). Where collected—in the Rio Grande aquifer system and High Plains aquifer—perchlorate did not occur at high concentrations; moderate aquifer-scale proportions were 3 and 11 percent, respectively. The understanding assessment included statistical tests to evaluate relations between constituent concentrations and potential explanatory factors to identify natural and human factors that affect groundwater quality. Potential explanatory factors included depth to bottom of well perforation, groundwater age category, land use, aquifer lithology, hydrologic conditions, and geochemical conditions. Higher concentrations of trace elements, radioactive constituents, and constituents with non-health-based benchmarks generally were associated with unconsolidated sand and gravel aquifer lithologies, premodern groundwater age, greater aridity, and more alkaline pH. Organic constituents with detection frequencies greater than 10 percent generally were associated with urban land use, shallower well depths, and higher total dissolved solids concentrations. The results for the six western PAs provide important insights into the quality of groundwater that is used for drinking water in the western United States, as well as natural and human factors that affect groundwater quality in this region.

Arizona, California, Colorado, Idaho, Kansas, Mont↗

USGS2020

An update of hydrologic conditions and distribution of selected constituents in water, Eastern Snake River Plain aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2016–18

Since 1952, wastewater discharged to infiltration ponds (also called percolation ponds) and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain (ESRP) aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, maintains groundwater-monitoring networks at the INL to determine hydrologic trends and to delineate the movement of radiochemical and chemical wastes in the aquifer and in perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from the ESRP aquifer and perched groundwater wells in the USGS groundwater monitoring networks during 2016–18. From March–May 2015 to March–May 2018, water levels in wells completed in the ESRP aquifer declined in the northern part of the INL and increased in the southwestern part. Water-level decreases ranged from 0.5 to 3.0 feet (ft) in the northern part of the INL and increases ranged from 0.5 to 3.0 ft in the southwestern part. Detectable concentrations of radiochemical constituents in water samples from wells in the ESRP aquifer at the INL generally decreased or remained constant during 2016–18. Decreases in concentrations were attributed to radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In 2018, concentrations of tritium in water samples collected from 46 of 111 aquifer wells were greater than the reporting level of three times the sample standard deviation and ranged from 260±50 to 5,100±190 picocuries per liter (pCi/L). Tritium concentrations in water from 10 wells completed in deep perched groundwater above the ESRP aquifer near the Advanced Test Reactor (ATR) Complex generally were greater than or equal to the reporting level during at least one sampling event during 2016–18, and concentrations ranged from 150 ±50 to 12,900 ±200 pCi/L. Concentrations of strontium-90 in water from 17 of 60 ESRP aquifer wells sampled during April or October 2018 exceeded the reporting level, ranging from 2.2±0.7 to 363±19 pCi/L. Strontium-90 was not detected in the ESRP aquifer beneath the ATR Complex. During at least one sampling event during 2016–18, concentrations of strontium-90 in water from eight wells completed in deep perched groundwater above the ESRP aquifer at the ATR Complex equaled or exceeded the reporting levels, and concentrations ranged from 0.57±0.17 to 34.3±1.2 pCi/L. During 2016–18, concentrations of cesium-137 were less than the reporting level in all but one ESRP aquifer well, and concentrations of plutonium-238, -239, and -240 (undivided), and americium-241 were less than the reporting level in water samples from all ESRP aquifer wells. In April 2009, the dissolved chromium concentration in water from one ESRP aquifer well, USGS 65, south of the ATR Complex equaled the maximum contaminant level (MCL) of 100 micrograms per liter (μg/L). In April 2018, the concentration of chromium in water from that well had decreased to 76.0 μg/L, less than the MCL. Concentrations in water samples from 62 other ESRP aquifer wells sampled ranged from less than 0.6 to 21.6 μg/L. During 2016–18, dissolved chromium was detected in water from all wells completed in deep perched groundwater above the ESRP aquifer at the ATR Complex, and concentrations ranged from 4.2 to 98.8 μg/L. In 2018, concentrations of sodium in water from most ESRP aquifer wells in the southern part of the INL were greater than the western tributary background concentration of 8.3 milligrams per liter (mg/L). After the new percolation ponds were put into service in 2002 southwest of the Idaho Nuclear Technology and Engineering Center (INTEC), concentrations of sodium in water samples from the Rifle Range well increased steadily until 2008, when concentrations generally began decreasing. The increases and decreases were attributed to disposal variability in the new percolation ponds. During 2016–18, dissolved sodium concentrations in water from 18 wells completed in deep perched groundwater above the ESRP aquifer at the ATR Complex ranged from 6.37 to 143 mg/L. In 2018, concentrations of chloride in most water samples from ESRP aquifer wells south of the INTEC and at the Central Facilities Area exceeded the background concentrations. Chloride concentrations in water from wells south of the INTEC generally have decreased since 2002 when chloride disposal to the old percolation ponds was discontinued. After the new percolation ponds southwest of the INTEC were put into service in 2002, concentrations of chloride in water samples from one well rose steadily until 2008 then began decreasing. During 2016–18, dissolved chloride concentrations in deep perched groundwater above the ESRP aquifer from 18 wells at the ATR Complex ranged from 3.89 to 176 mg/L. In 2018, sulfate concentrations in water samples from ESRP aquifer wells in the south-central part of the INL exceeded the background concentration of sulfate and ranged from 22 to 151 mg/L. The greater-than-background concentrations in water from these wells probably resulted from sulfate disposal at the ATR Complex infiltration ponds or the old INTEC percolation ponds. In 2018, sulfate concentrations in water samples from wells near the Radioactive Waste Management Complex (RWMC) mostly were greater than background concentrations and could have resulted from well construction techniques and (or) waste disposal at the RWMC or the ATR complex. The maximum dissolved sulfate concentration in shallow perched groundwater above the ESRP aquifer near the ATR Complex was 215 mg/L in well CWP 3 in April 2016. During 2018, dissolved sulfate concentrations in water from wells completed in deep perched groundwater above the ESRP aquifer near the cold-waste ponds at the ATR Complex ranged from 65.8 to 171 mg/L. In 2018, concentrations of nitrate in water from most ESRP aquifer wells at and near the INTEC exceeded the western tributary background concentration of 0.655 mg/L. Concentrations of nitrate in wells southwest of the INTEC and farther away from the influence of disposal areas and the Big Lost River show a general decrease in nitrate concentration through time. Two wells south of the INTEC show increasing trends that could be the result of wastewater beneath the INTEC tank farm being mobilized to the aquifer. During 2016–18, water samples from several ESRP aquifer wells were collected and analyzed for volatile organic compounds (VOCs). Sixteen VOCs were detected. At least 1 and as many as 7 VOCs were detected in water samples from 15 wells. The primary VOCs detected include carbon tetrachloride, trichloromethane, tetrachloroethene, 1,1,1-trichloroethane, and trichloroethene. In 2016–18, concentrations for all VOCs were less than their respective MCLs for drinking water, except carbon tetrachloride in water from two wells and trichloroethene in one well. During 2016–18, variability and bias were evaluated from 37 replicate and 15 blank quality-assurance samples. Results from replicate analyses were investigated to evaluate sample variability. Constituents with acceptable reproducibility were major ions, trace elements, nutrients, and VOCs. All radiochemical constituents had acceptable reproducibility except for gross alpha- and beta-particle radioactivity. The gross alpha- and beta-particle radioactivity samples that did not meet reproducibility criteria had low concentrations. Bias from sample contamination was evaluated from equipment, field, and source-solution blanks. Cadmium had a concentration slightly greater than its reporting level in a source-solution blank, and chloride and ammonia had concentrations that were slightly greater than their respective reporting levels in field and equipment blanks. Subtracting concentrations of chloride and ammonia in field blanks from the concurrently collected equipment blank indicates that adjusted concentrations for chloride and ammonia in the equipment blanks were less than their respective reporting levels. Therefore, no sample bias was observed for any of the sample periods.

Idaho↗

USGS2020

Groundwater availability of the Northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

Executive Summary The Northern High Plains aquifer underlies about 93,000 square miles of Colorado, Kansas, Nebraska, South Dakota, and Wyoming and is the largest subregion of the nationally important High Plains aquifer. Irrigation, primarily using groundwater, has supported agricultural production since before 1940, resulting in nearly $50 billion in sales in 2012. In 2010, the High Plains aquifer had the largest groundwater withdrawals of any major aquifer system in the United States. Nearly one-half of those withdrawals were from the Northern High Plains aquifer, which has little hydrologic interaction with parts of the aquifer farther south. Land-surface elevation ranges from more than 7,400 feet (ft) near the western edge to less than 1,100 ft near the eastern edge. Major stream primarily flow west to east and include the Big Blue River, Elkhorn River, Loup River, Niobrara River, Republican River and Platte River with its two forks—the North Platte River and South Platte River. Population in the Northern High Plain aquifer area is sparse with only 2 cities having a population greater than 30,000. Droughts across much of the area from 2001 to 2007, combined with recent (2004–18) legislation, have heightened concerns regarding future groundwater availability and highlighted the need for science-based water-resource management. Groundwater models with the capability to provide forecasts of groundwater availability and related stream base flows from the Northern High Plains aquifer were published recently (2016) and were used to analyze groundwater availability. Stream base flows are generally the dominant component of total streamflow in the Northern High Plains aquifer, and total streamflows or shortages thereof define conjunctive management triggers, at least in Nebraska. Groundwater availability was evaluated through comparison of aquifer-scale water budgets compared for periods before and after major groundwater development and across selected future forecasts. Groundwater-level declines and the forecast amount of groundwater in storage in the aquifer also were examined. Major Findings Aquifer losses to irrigation withdrawals increased greatly from 1940 to 2009 and were the largest average 2000–9 outflow (49 percent of total). Basin to basin groundwater flows were not a large part of basin water budgets. Development of irrigated land and associated withdrawals were not uniform across the Northern High Plains aquifer, and different parts of the Northern High Plains aquifer responded differently to agricultural development. For the Northern High Plains aquifer, areas with high recharge and low evapotranspiration had the most streamflow, and most streams only remove water from the aquifer. Results of a baseline future forecast indicated that groundwater levels declined overall, indicating an overdraft of the aquifer when climate was about average and agricultural development was held at the same state as 2009. Results of two human stresses future forecasts indicated that increases of 13 percent or 23 percent in agricultural development, mostly near areas of previous development, caused increases in groundwater pumping of 8 percent or 11 percent, and resulted in continued groundwater-level declines, at rates 0.3 or 0.5 million acre-feet per year larger than the baseline forecast. Results of environmental stresses forecasts (generated from two downscalings of global climate model outputs) compared with the baseline forecast indicated that even though annual precipitation was nearly the same, differences in temperature and a redistribution of precipitation from the spring to the growing season (from about May 1 through September 30), created a large (12–15 percent) decrease in recharge to the aquifer. For the two environmental stresses forecasts, temperature and precipitation were distributed about the same among basins of the Northern High Plains aquifer, but the amounts were different.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗

USGS2013

An update of hydrologic conditions and distribution of selected constituents in water, eastern Snake River Plain aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2009–11

Since 1952, wastewater discharged to infiltration ponds (also called percolation ponds) and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain (ESRP) aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, maintains groundwater monitoring networks at the INL to determine hydrologic trends, and to delineate the movement of radiochemical and chemical wastes in the aquifer and in perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from aquifer, multilevel monitoring system (MLMS), and perched groundwater wells in the USGS groundwater monitoring networks during 2009–11. Water in the ESRP aquifer primarily moves through fractures and interflow zones in basalt, generally flows southwestward, and eventually discharges at springs along the Snake River. The aquifer primarily is recharged from infiltration of irrigation water, infiltration of streamflow, groundwater inflow from adjoining mountain drainage basins, and infiltration of precipitation. From March–May 2009 to March–May 2011, water levels in wells generally declined in the northern part of the INL. Water levels generally rose in the central and eastern parts of the INL. Detectable concentrations of radiochemical constituents in water samples from aquifer wells or MLMS equipped wells in the ESRP aquifer at the INL generally decreased or remained constant during 2009–11. Decreases in concentrations were attributed to radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In 2011, concentrations of tritium in groundwater from 50 of 127 aquifer wells were greater than or equal to the reporting level and ranged from 200±60 to 7,000±260 picocuries per liter. Tritium concentrations from one or more discrete zones from four wells equipped with MLMS were greater than or equal to reporting levels in water samples collected at various depths. Tritium concentrations in water from wells completed in shallow perched groundwater at the Advanced Test Reactor Complex (ATR Complex) were less than the reporting levels. Tritium concentrations in deep perched groundwater at the ATR Complex equaled or exceeded the reporting level in 12 wells during at least one sampling event during 2009–11 at the ATR Complex. Concentrations of strontium-90 in water from 20 of 76 aquifer wells sampled during April or October 2011 exceeded the reporting level. Strontium-90 was not detected within the ESRP aquifer beneath the ATR Complex. During at least one sampling event during 2009–11, concentrations of strontium-90 in water from 10 wells completed in deep perched groundwater at the ATR Complex equaled or exceeded the reporting levels. During 2009–11, concentrations of plutonium-238, and plutonium-239, -240 (undivided), and americium-241 were less than the reporting level in water samples from all aquifer wells and in all wells equipped with MLMS. Concentrations of cesium-137 were equal to or slightly above the reporting level in 8 aquifer wells and from 2 wells equipped with MLMS. The concentration of chromium in water from one well south of the ATR Complex was 97 micrograms per liter (μg/L) in April 2011, just less than the maximum contaminant level (MCL) of 100 μg/L. Concentrations of chromium in water samples from 69 other wells sampled ranged from 0.8 μg/L to 25 μg/L. During 2009–11, dissolved chromium was detected in water from 15 wells completed in perched groundwater at the ATR Complex. In 2011, concentrations of sodium in water from most wells in the southern part of the INL were greater than the background concentration of 10 milligrams per liter (mg/L); the highest concentrations were at or near the Idaho Nuclear Engineering and Technology Center (INTEC). After the newpercolation ponds were put into service in 2002 southwest of the INTEC, concentrations of sodium in water samples from the Rifle Range well rose steadily until 2008, when the concentrations generally began decreasing. The increases and decreases were attributed to disposal variability in the new percolation ponds. Concentrations of sodium in most wells equipped with MLMS generally were consistent with depth. During 2011, dissolved sodium concentrations in water from 17 wells completed in deep perched groundwater at the ATR Complex ranged from 6 to 146 mg/L. In 2011, concentrations of chloride in most water samples from aquifer wells south of the INTEC and at the Central Facilities Area exceeded the background concentrations of 15 mg/L, but were less than the secondary MCL of 250 mg/L. Chloride concentrations in water from wells south of the INTEC have generally increased because of increased chloride disposal to the old percolation ponds since 1984 when discharge of wastewater to the INTEC disposal well was discontinued. After the new percolation ponds were put into service in 2002 southwest of the INTEC, concentrations of chloride in water samples from one well rose steadily until 2008 then began decreasing. Chloride concentrations in water from all but one well completed in the ESRP aquifer at or near the ATR Complex were less than background and ranged between 10 and 14 mg/L during 2011, similar to concentrations detected during the 2006–08 reporting period. During 2011, chloride concentrations in water from two aquifer wells at the Radioactive Waste Management Complex (RWMC) were slightly greater than concentrations detected during the 2006–08 reporting period. The vertical distribution of chloride concentrations in wells equipped with MLMS were generally consistent within zones during 2009–11 and ranged from about 8 to 20 mg/L. During April 2011, dissolved chloride concentrations in shallow perched groundwater at the ATR Complex ranged from 7 to 13 mg/L in water from three wells. Dissolved chloride concentrations in deep perched groundwater at the ATR Complex during 2011 ranged from 4 to 54 mg/L. In 2011, sulfate concentrations in water samples from 11 aquifer wells in the south-central part of the INL equaled or exceeded the background concentration of sulfate and ranged from 40 to 167 mg/L. The greater-than-background concentrations in water from these wells probably resulted from sulfate disposal at the ATR Complex infiltration ponds or the old INTEC percolation ponds. In 2011, sulfate concentrations in water samples from two wells near the RWMC were greater than background levels and could have resulted from well construction techniques and (or) waste disposal at the RWMC. The vertical distribution of sulfate concentrations in three wells near the southern boundary of the INL was generally consistent with depth, and ranged between 19 and 25 mg/L. The maximum dissolved sulfate concentration in shallow perched groundwater near the ATR Complex was 400 mg/L in well CWP 1 in April 2011. During 2009–11, the maximum concentration of dissolved sulfate in deep perched groundwater at the ATR Complex was 1,550 mg/L in a well located west of the chemical-waste pond. In 2011, concentrations of nitrate in water from most wells at and near the INTEC exceeded the regional background concentrations of 1 mg/L and ranged from 1.6 to 5.95 mg/L. Concentrations of nitrate in wells south of INTEC and farther away from the influence of disposal areas and the Big Lost River show a general decrease in nitrate concentrations through time. During 2009–11, water samples from 30 wells were collected and analyzed for volatile organic compounds (VOCs). Six VOCs were detected. At least one and up to five VOCs were detected in water samples from 10 wells. The primary VOCs detected include carbon tetrachloride, chloroform, tetrachloroethylene, 1,1,1-trichloroethane, and trichloroethylene. In 2011, concentrations for all VOCs were less than their respective MCL for drinking water, except carbon tetrachloride in water from two wells. During 2009–11, variability and bias were evaluated from 56 replicate and 16 blank quality-assurance samples. Results from replicate analyses were investigated to evaluate sample variability. Constituents with acceptable reproducibility were stable isotope ratios, major ions, nutrients, and VOCs. All radiochemical constituents and trace metals had acceptable reproducibility except for gross beta-particle radioactivity, aluminum, antimony, and cobalt. Bias from sample contamination was evaluated from equipment, field, container, and source-solution blanks. No detectable constituent concentrations were reported for equipment blanks of the thief samplers and sampling pipes or for the source-solution and field blanks. Equipment blanks of bailers had detectable concentrations of strontium-90, sodium, chloride, and sulfate, and the container blank had a detectable concentration of dichloromethane.

Idaho↗

USGS2014

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

The Edwards-Trinity aquifer is a vital groundwater resource for agricultural, industrial, and public supply uses in the Pecos County region of western Texas. The U.S. Geological Survey completed a comprehensive, integrated analysis of available hydrogeologic data to develop a numerical groundwater-flow model of the Edwards-Trinity and related aquifers in the study area in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. The active model area covers about 3,400 square miles of the Pecos County region of Texas west of the Pecos River, and its boundaries were defined to include the saturated areas of the Edwards-Trinity aquifer. The model is a five-layer representation of the Pecos Valley, Edwards-Trinity, Dockum, and Rustler aquifers. The Pecos Valley aquifer is referred to as the alluvial layer, and the Edwards-Trinity aquifer is divided into layers representing the Edwards part of the Edwards-Trinity aquifer and the Trinity part of the Edwards-Trinity aquifer, respectively. The calibration period of the simulation extends from 1940 to 2010. Simulated hydraulic heads generally were in good agreement with observed values; 1,684 out of 2,860 (59 percent) of the simulated values were within 25 feet of the observed value. The average root mean square error value of hydraulic head for the Edwards-Trinity aquifer was 34.2 feet, which was approximately 4 percent of the average total observed change in groundwater-level altitude (groundwater level). Simulated spring flow representing Comanche Springs exhibits a pattern similar to observed spring flow. Independent geochemical modeling corroborates results of simulated groundwater flow that indicates groundwater in the Edwards-Trinity aquifer in the Leon-Belding and Fort Stockton areas is a mixture of recharge from the Barilla and Davis Mountains and groundwater that has upwelled from the Rustler aquifer. The model was used to simulate groundwater-level altitudes resulting from prolonged pumping to evaluate sustainability of current and projected water-use demands. Each of three scenarios utilized a continuation of the calibrated model. Scenario 1 extended recent (2008) irrigation and nonirrigation pumping values for a 30-year period from 2010 to 2040. Projected groundwater-level changes in and around the Fort Stockton area under scenario 1 change little from current conditions, indicating that the groundwater system is near equilibrium with respect to recent (2008) pumping stress. Projected groundwater-level declines in the eastern part of the model area ranging from 5.0 to 15.0 feet are likely the result of nonequilibrium conditions associated with recent increases in pumping after a prolonged water-level recovery period of little or no pumping. Projected groundwater-level declines (from 15.0 to 31.0 feet) occurred in localized areas by the end of scenario 1 in the Leon-Belding area. Scenario 2 evaluated the effects of extended recent (2008) pumping rates as assigned in scenario 1 with year-round maximum permitted pumping rates in the Belding area. Results of scenario 2 are similar in water-level decline and extent as those of scenario 1. The extent of the projected groundwater-level decline in the range from 5.0 to 15.0 feet in the Leon-Belding irrigation area expanded slightly (about a 2-percent increase) from that of scenario 1. Maximum projected groundwater-level declines in the Leon-Belding irrigation area were approximately 31.3 feet in small isolated areas. Scenario 3 evaluated the effects of periodic increases in pumping rates over the 30-year extended period. Results of scenario 3 are similar to those of scenario 2 in terms of the areas of groundwater-level decline; however, the maximum projected groundwater-level decline increased to approximately 34.5 feet in the Leon-Belding area, and the extent of the decline was larger in area (about a 17-percent increase) than that of scenario 2. Additionally, the area of projected groundwater-level declines in the eastern part of the model area increased from that of scenario 2—two individual areas of decline coalesced into one larger area. The localized nature of the projected groundwater-level declines is a reflection of the high degree of fractured control on storage and hydraulic conductivity in the Edwards-Trinity aquifer. Additionally, the finding that simulated spring flow is highly dependent on the transient nature of hydraulic heads in the underlying aquifer indicates the importance of adequately understanding and characterizing the entire groundwater system.

Texas↗

USGS2012

A conceptual model of the hydrogeologic framework, geochemistry, and groundwater-flow system of the Edwards-Trinity and related aquifers in the Pecos County region, Texas

A conceptual model of the hydrogeologic framework, geochemistry, and groundwater-flow system of the Edwards-Trinity and related aquifers, which include the Pecos Valley, Igneous, Dockum, Rustler, and Capitan Reef aquifers, was developed as the second phase of a groundwater availability study in the Pecos County region in west Texas. The first phase of the study was to collect and compile groundwater, surface-water, water-quality, geophysical, and geologic data in the area. The third phase of the study involves a numerical groundwater-flow model of the Edwards-Trinity aquifer in order to simulate groundwater conditions based on various groundwater-withdrawal scenarios. Resource managers plan to use the results of the study to establish management strategies for the groundwater system. The hydrogeologic framework is composed of the hydrostratigraphy, structural features, and hydraulic properties of the groundwater system. Well and geophysical logs were interpreted to define the top and base surfaces of the Edwards-Trinity aquifer units. Elevations of the top and base of the Edwards-Trinity aquifer generally decrease from the southwestern part of the study area to the northeast. The thicknesses of the Edwards-Trinity aquifer units were calculated using the interpolated top and base surfaces of the hydrostratigraphic units. Some of the thinnest sections of the aquifer were in the eastern part of the study area and some of the thickest sections were in the Pecos, Monument Draw, and Belding-Coyanosa trough areas. Normal-fault zones, which formed as growth and collapse features as sediments were deposited along the margins of more resistant rocks and as overlying sediments collapsed into the voids created by the dissolution of Permian-age evaporite deposits, were delineated based on the interpretation of hydrostratigraphic cross sections. The lowest aquifer transmissivity values were measured in the eastern part of the study area; the highest transmissivity values were measured in a faulted area of the Monument Draw trough. Hydraulic conductivity values generally exhibited the same trends as the transmissivity values. Groundwater-quality data and groundwater-level data were used in context with the hydrogeologic framework to assess the chemical characteristics of water from different sources, regional groundwater-flow paths, recharge sources, the mixing of water from different sources, and discharge in the study area. Groundwater-level altitudes generally decrease from southwest to northeast and regional groundwater flow is from areas of recharge south and west to the north and northeast. Four principal sources of recharge to the Edwards-Trinity aquifer were identified: (1) regional flow that originated as recharge northwest of the study area, (2) runoff from the Barilla, Davis, and Glass Mountains, (3) return flow from irrigation, and (4) upwelling from deeper aquifers. Results indicated Edwards-Trinity aquifer water in the study area was dominated by mineralized, regional groundwater flow that most likely recharged during the cooler, wetter climates of the Pleistocene with variable contributions of recent, local recharge. Groundwater generally flows into the down-dip extent of the Edwards-Trinity aquifer where it discharges into overlying or underlying aquifer units, discharges from springs, discharges to the Pecos River, follows a regional flow path east out of the study area, or is withdrawn by groundwater wells. Structural features such as mountains, troughs, and faults play a substantial role in the distribution of recharge, local and regional groundwater flow, spring discharge, and aquifer interaction.

Texas↗

USGS2018

Simulation of groundwater flow and analysis of projected water use for the Rush Springs aquifer, western Oklahoma

The U.S. Geological Survey, in cooperation with the Bureau of Reclamation and the Oklahoma Water Resources Board, (1) quantified the groundwater resources of the Rush Springs aquifer in western Oklahoma by developing a numerical groundwater-flow model, (2) evaluated the effects of estimated equal-proportionate-share (EPS) pumping rates on aquifer storage and streamflow for time periods of 20, 40, and 50 years into the future, (3) assessed the uncertainty in the EPS scenario results, and (4) evaluated the effects of (a) projected groundwater-use rates extended 50 years into the future and (b) sustained hypothetical drought conditions over a 10-year period on stream base flow and groundwater in storage. The Rush Springs aquifer is an important source of water for municipal and irrigation use by many communities and agricultural users in the study area. The study area is composed of about 4,970 square miles (3,181,003 acres) of Rush Springs aquifer bedrock deposits located in 14 counties. The study area also includes the alluvium and terrace deposits of the Canadian and Washita Rivers, as well as alluvium along the Little Washita River, Deer Creek, and a number of smaller tributaries of the Washita River that overlie the bedrock. A numerical groundwater-flow model of the Rush Springs aquifer was constructed by using MODFLOW with the Newton solver. Groundwater flow was simulated for January 1979–December 2015 by using monthly stress periods, and an initial steady-state stress period was configured to represent mean annual inflows and outflows. The model was calibrated to groundwater-level observations at selected wells, monthly base flow at nine streamgages, stream seepage as estimated for the conceptual water budget, and Fort Cobb Reservoir stage. The EPS scenarios for the Rush Springs 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.82, 0.49, and 0.43 acre-foot per acre per year, respectively. Given the 2,954,545-acre aquifer area used for the EPS scenarios, the 20-year rate corresponds to an annual yield of about 2,422,727 acre-feet per year. Groundwater storage at the end of the 20-year EPS scenario was about 13,321,000 acre-feet, or about 31,516,437 acre-feet (70 percent) less than the starting EPS scenario storage. This decrease in storage was equivalent to a mean groundwater-level decline of about 152 feet. Water availability under the EPS pumping rate was primarily from the western area of the model. Saturation was sustained though the entire EPS scenario where the aquifer was sufficiently thick or a shallow hydraulic gradient was present. Fort Cobb Reservoir stage was below the dead-pool stage after about 5 years of 20-year EPS pumping. An uncertainty analysis was conducted to assess the uncertainty in the EPS scenario results. An ensemble of 400 random sets of possible parameter values was performed for the uncertainty analysis by using a multivariate normal distribution centered on the calibrated parameter values. The parameter bounds for the uncertainty analysis were determined by using the posterior covariance matrix, which allows for the incorporation of knowledge gained during the calibration process as well as observation uncertainty and the correlation between estimated parameters. The uncertainty results indicate a 95-percent confidence interval for the 20-year EPS pumping rate between 0.73 and 0.95 acre-foot per acre per year. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Rush Springs aquifer. The effects of well withdrawals were evaluated by comparing changes in groundwater storage between four 50-year scenarios using (1) no groundwater use, (2) mean groundwater use for the study period (1979–2015), (3) increasing groundwater use, and (4) groundwater use at the 2015 rate. The increasing-use scenario assumed a 38-percent increase in pumping over 50 years on the basis of 2010–60 demand projections for western Oklahoma. Simulated groundwater storage changes ranged between an increase of 6.3 percent for the scenario with no groundwater use, and 0.9 percent for the scenario with 2015 groundwater-use rates. For the Fort Cobb Reservoir surface watershed, simulated groundwater storage changes ranged between an increase of 23.6 percent for the scenario with no groundwater use and a decrease of 4.0 percent for the increasing groundwater-use scenario. Groundwater-level changes were generally greater in areas with a large concentration of groundwater wells and groundwater use such as the Fort Cobb Reservoir surface watershed. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on the Rush Springs aquifer groundwater storage and Fort Cobb Reservoir stage and storage. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model. To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–1992), and upstream inflows to the Canadian and Washita Rivers and associated tributaries were reduced by 37 percent. Groundwater storage at the end of the hypothetical drought period in December 1992 was about 42,983,000 acre-feet, or about 3,525,000 acre-feet (7.6 percent) less than the groundwater storage of the calibrated numerical model. This change in groundwater storage is equivalent to a mean groundwater-level decline of 15.8 feet. Simulated mean base-flow declines at the Canadian and Washita River streamgages were between 39 and 59 percent during the drought period. The minimum stage in Fort Cobb Reservoir at the end of the hypothetical drought period was 1,311 feet, indicating a storage capacity of only 10 percent of active conservation pool storage. The Fort Cobb Reservoir storage declines mostly resulted from reduced base flows in Cobb, Lake, and Willow Creeks upstream from the reservoir.

Oklahoma↗

USGS2021

Status and understanding of groundwater quality in the northern Sierra Nevada foothills domestic-supply aquifer study units, 2015–17—California GAMA Priority Basin Project

Groundwater quality in the northern Sierra Nevada foothills region of California was investigated as part of California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring Assessment Priority Basin Project (GAMA-PBP). The region was divided into two study units: the Yuba-Bear watersheds (YBW) study unit and the American-Cosumnes-Mokelumne watersheds (ACMW) study unit. The GAMA-PBP made a spatially unbiased assessment of aquifer systems used for domestic drinking-water supply in the study region, which are predominantly composed of fractured, hard-rock aquifers of varying lithology. These assessments characterized the quality of raw groundwater to evaluate ambient conditions in the domestic-supply aquifer and not the quality of treated drinking water. The study included three components: (1) a status assessment , which characterized the quality of groundwater resources used for domestic drinking-water supply in the YBW and ACMW study units; (2) an understanding assessment , which evaluated natural and anthropogenic explanatory factors that could potentially affect groundwater quality in the study region; and (3) a comparative assessment between the groundwater resources used for domestic and public drinking-water supply in the study region. The status assessment was based on data collected by the GAMA-PBP from 74 sites in the YBW study unit during 2015–16 and 67 sites in the ACMW study unit from 2016 to 2017. To contextualize water-quality results, concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. The status assessment used a grid-based method to estimate proportions of groundwater resources with concentrations approaching or exceeding benchmark thresholds. This method provides spatially unbiased results and allows inter-comparability with similar groundwater-quality assessments. Inorganic constituents with health-based benchmarks were present at high relative concentration (RC), meaning they exceeded the benchmark threshold, in 5.4 and 10 percent of domestic-supply aquifer systems in the YBW and ACMW study units, respectively. Inorganic constituents with aesthetic-based benchmarks were detected at high-RCs in 20 and 28 percent of the YBW and ACMW study units, respectively. The inorganic constituents present at high RC were arsenic, barium, boron, molybdenum, strontium, nitrate, adjusted gross-alpha particle activity, chloride, total dissolved solids, specific conductance, iron, manganese, and hardness. Groundwater samples were tested for presence or absence of three microbial indicators (total coliform, Escherichia coli , and Enterococci ). At least one microbial indicator was present in 26 and 28 percent of the YBW and ACMW study units, respectively. At least one organic constituent was detected in 30 and 42 percent of the YBW and ACMW study units, respectively. Organic constituents were not present at high RC, but tetrachloroethene (PCE), trichloroethene (TCE), and toluene were detected in the YBW study unit at moderate RC (between the benchmark concentration and one-tenth of the benchmark concentration). Methyl tert -butyl ether (MTBE) and chloroform were present at low RC (less than one-tenth of the benchmark concentration) in the YBW and ACMW study units with detection frequencies greater than 10 percent. Perchlorate, a constituent of special interest, was detected in 31 and 41 percent of the YBW and ACMW study units, respectively, at either low or moderate RCs. Relations among select water-quality constituents and potential explanatory factors were evaluated using statistical and graphical approaches. Nitrate, microbial indicators, and perchlorate were all correlated to elevation-dependent variables relating to climate, land use, and recharge condition. Isotopic and dissolved noble-gas tracers indicated these water-quality constituents are associated with recharge conditions associated with irrigation during the summer dry-season, which is common in areas of rural-residential or agricultural land uses. Higher concentrations of iron and manganese were primarily associated with anoxic groundwater in aquifers of metasedimentary lithology. Increased hardness was primarily associated with anoxic groundwater in aquifers of mafic-ultramafic or metavolcanics lithologies at lower elevations in the study region in the Melones fault zone. Chloroform and MTBE were associated with shallow groundwater (wells depths less than 130 m) under oxic and anoxic redox conditions, respectively. The comparative assessment evaluated differences between the aquifer systems used for domestic- and public-supply in study region based on (1) well-construction characteristics, and (2) water quality. Analysis of over 60,000 well-completion reports in the study region showed that although domestic-supply wells span the deepest depth zones in regional aquifers, median depths for public-supply wells were significantly greater than those of domestic-supply wells in both study units. Water-quality data from more than 300 public-supply wells in the study region were assessed using a spatially weighted method for calculation aquifer-scale proportions and compared with the domestic-supply assessment results. Detections of inorganic constituents at high RC and detection frequencies for organic constituents were generally similar between the domestic- and public-supply aquifer systems in both study units, with a few notable exceptions in the ACMW study unit: nitrate was greater for the public- compared to domestic-supply aquifer system and both manganese, hardness, and MTBE were greater in the domestic- compared to public-supply aquifer system. These differences are likely related to contrasting land uses, aquifer lithologies, landscape positions, and depths characterizing domestic- and public-supply wells in the ACMW study unit. Overall, fewer samples from domestic-supply wells in the northern Sierra Nevada foothills exceeded health-based benchmarks compared to aesthetic-based benchmarks for groundwater quality. Exceedences of health-based benchmarks were primarily caused by nitrate and coliform bacteria, which were associated with recharge from diverted surface water used primarily for irrigation. Exceedences of aesthetic-based benchmarks were primarily caused by iron, managanese, and hardness, which were associated with geologic factors. Regional irrigation practices and aquifer lithology can affect groundwater quality in fractured-rock aquifers in the northern Sierra Nevada foothills used for domestic drinking-water supply.

California↗

USGS2018

Revised groundwater-flow model of the glacial aquifer system north of Aberdeen, South Dakota, through water year 2015

The city of Aberdeen, in northeastern South Dakota, requires an expanded and sustainable supply of water to meet current and future demands. Conceptual and numerical models of the glacial aquifer system in the area north of Aberdeen were developed by the U.S. Geological Survey in cooperation with the City of Aberdeen in 2012. The U.S. Geological Survey, in cooperation with the City of Aberdeen, completed a study to revise the original numerical groundwater-flow model using data through water year (WY) 2015 to aid the City of Aberdeen in their development of plans and strategies for a sustainable water supply and to increase understanding of the glacial aquifer system and groundwater-flow system near Aberdeen. The original model was revised to improve the fit between model-simulated values and observed (measured or estimated) data, provide greater insight into surface-water interactions, and improve the usefulness of the model for water-supply planning. The revised groundwater-flow model (hereafter referred to as the “revised model”) presented in this report supersedes the original model. The purpose of this report is to describe a revised groundwater-flow model including data collection, model calibration, and model results for the glacial aquifer system including the Elm, Middle James, and Deep James aquifers north of Aberdeen, South Dakota, using updated hydrologic data through WY 2015. The original numerical model was revised in several ways. The model was modified by adding four new layers, which included a surficial layer, two intervening confining layers, and a shale bedrock layer. The revised model provides an improved understanding of the groundwater-flow system in comparison to the original model. The principal aquifers of the model area include portions of the Elm, Middle James, and Deep James aquifers. The lithologic information used to define and describe the aquifers in the model area was unaltered; however, aquifer properties and boundary conditions were reviewed and updated using geological information reported by the South Dakota Department of Environmental and Natural Resources and information obtained from geophysical investigations for this study. The horizontal extent of the Elm, Middle James, and Deep James aquifers was unaltered from the original model. The thickness of the Deep James aquifer was modified based on interpretations from the geophysical investigations. In general, groundwater in the Elm aquifer flowed from northwest to southeast and locally towards rivers and streams. Similarly, in the Middle James and Deep James aquifers, groundwater also typically flowed southeast. The revisions made to the original model include use of the following MODFLOW stress packages: Recharge, Evapotranspiration, Time-Variant Specified Head, Wells, Drains, and Stream Flow Routing, all of which were updated from the original model except for the Stream Flow Routing Package, which replaced the River Package used in the original model. Model calibration is the process of estimating model parameters to minimize the differences, or residuals, between observed data and simulated values; therefore, Parameter ESTimation (PEST) software was used to optimize model input parameters by matching model-simulated values to observed data. Calibration parameters included horizontal hydraulic conductivity, vertical hydraulic conductivity, specific yield, specific storage, and vertical streambed conductance for stream and drain cells. Multipliers were used to calibrate the recharge and evapotranspiration stresses. Evapotranspiration extinction depth also was adjusted during model calibration. Comparisons to the original model are described to highlight the changes made in the revised model. In general, the revised model adequately simulates the natural system and compares favorably with observed hydrologic data. Simulated water levels were evaluated by comparing them to single water-level observations at selected well locations. The selected wells were the same wells used in the original model. The coefficient of determination value between simulated and observed water levels for the revised model was 0.89 and included simulated and observed values from October 1, 1974 (WY 1975), through September 30, 2015 (WY 2015). The coefficient of determination value for the original model was 0.94 and included simulated and observed values from October 1, 1974, through September 30, 2009. The difference may indicate that the original model could have been overfit to hydraulic head observations because base flow was not simulated. The additional data used in the revised model included some climatically wetter, more extreme periods, such as 2011, in which annual precipitation was 30.9 inches. Average annual precipitation for the original model timeframe, which included data from WYs 1975–2009, was 20.26 inches. Additional precipitation data for WYs 2010–15, included in the revised model timeframe, resulted in an average annual precipitation for WYs 1975–2015 in the model area of 20.6 inches. The larger variability in climate data coupled with the additional water-level data could explain the lower coefficient of determination for water levels in the revised model. The revised model was used to calculate various groundwater-budget components for steady-state and transient conditions for WYs 1975–2015. The time-variant specified-head cells in the revised model had the largest change when compared to the original steady-state model for inflows and outflows. Comparing the transient budget components between the original and the revised models indicated that inflow from recharge and time-variant specified-head cells had the greatest effect on groundwater inflows, and outflow from storage had the greatest effect on groundwater outflows. The simulated potentiometric contours from the revised model were compared with (1) the observed (interpreted) potentiometric surface (layer 2) and the hydraulic head values (layers 4 and 6) and (2) the simulated contours from the original model. The simulated hydraulic gradients and general direction of groundwater flow in the Elm aquifer in the revised model generally matched the observed potentiometric contours, the simulated potentiometric contours from the original model, and general flow directions interpreted to be perpendicular to the contours. Minor discrepancies between simulated potentiometric contours from the revised model and the observed potentiometric contours may be due to the lack of observed data in the model area. The revised model was designed to reduce the limitations of the original model. The revisions were validated by comparing the results of the original model with the revised model. A primary benefit of the revised model is the inclusion of the surficial deposits and the confining units as explicit layers in the model. The addition of the surficial layer was beneficial for three primary reasons: (1) more accurate representation of recharge from precipitation, (2) more accurate representation of groundwater evapotranspiration, and (3) more accurate representation of groundwater and surface-water interactions. The groundwater model is a numeric approximation of a complex physical hydrologic system, and the revised model data were interpolated in regions with sparse data. Additionally, model discretization included averaged and interpolated values for water use, withdrawal rates, and hydraulic conductivity. The revised model provides a useful estimate for hydraulic gradients, groundwater-flow directions, and aquifer response to groundwater withdrawals.

South Dakota↗

USGS2009

Simulation of Groundwater Flow in the Coastal Plain Aquifer System of Virginia

The groundwater model documented in this report simulates the transient evolution of water levels in the aquifers and confining units of the Virginia Coastal Plain and adjacent portions of Maryland and North Carolina since 1890. Groundwater withdrawals have lowered water levels in Virginia Coastal Plain aquifers and have resulted in drawdown in the Potomac aquifer exceeding 200 feet in some areas. The discovery of the Chesapeake Bay impact crater and a revised conceptualization of the Potomac aquifer are two major changes to the hydrogeologic framework that have been incorporated into the groundwater model. The spatial scale of the model was selected on the basis of the primary function of the model of assessing the regional water-level responses of the confined aquifers beneath the Coastal Plain. The local horizontal groundwater flow through the surficial aquifer is not intended to be accurately simulated. Representation of recharge, evapotranspiration, and interaction with surface-water features, such as major rivers, lakes, the Chesapeake Bay, and the Atlantic Ocean, enable simulation of shallow flow-system details that influence locations of recharge to and discharge from the deeper confined flow system. The increased density of groundwater associated with the transition from fresh to salty groundwater near the Atlantic Ocean affects regional groundwater flow and was simulated with the Variable Density Flow Process of SEAWAT (a U.S. Geological Survey program for simulation of three-dimensional variable-density groundwater flow and transport). The groundwater density distribution was generated by a separate 108,000-year simulation of Pleistocene freshwater flushing around the Chesapeake Bay impact crater during transient sea-level changes. Specified-flux boundaries simulate increasing groundwater underflow out of the model domain into Maryland and minor underflow from the Piedmont Province into the model domain. Reported withdrawals accounted for approximately 75 percent of the total groundwater withdrawn from Coastal Plain aquifers during the year 2000. Unreported self-supplied withdrawals were simulated in the groundwater model by specifying their probable locations, magnitudes, and aquifer assignments on the basis of a separate study of domestic-well characteristics in Virginia. The groundwater flow model was calibrated to 7,183 historic water-level observations from 497 observation wells with the parameter-estimation codes UCODE-2005 and PEST. Most water-level observations were from the Potomac aquifer system, which permitted a more complex spatial distribution of simulated hydraulic conductivity within the Potomac aquifer than was possible for other aquifers. Zone, function, and pilot-point approaches were used to distribute assigned hydraulic properties within the aquifer system. The good fit (root mean square error = 3.6 feet) of simulated to observed water levels and reasonableness of the estimated parameter values indicate the model is a good representation of the physical groundwater flow system. The magnitudes and temporal and spatial distributions of residuals indicate no appreciable model bias. The model is intended to be useful for predicting changes in regional groundwater levels in the confined aquifer system in response to future pumping. Because the transient release of water stored in low-permeability confining units is simulated, drawdowns resulting from simulated pumping stresses may change substantially through time before reaching steady state. Consequently, transient simulations of water levels at different future times will be more accurate than a steady-state simulation for evaluating probable future aquifer-system responses to proposed pumping.

Scientific Investigations Report↗

USGS2017

Estimation of the groundwater resources of the bedrock aquifers at the Kettle Moraine Springs State Fish Hatchery, Sheboygan County, Wisconsin

Groundwater resources information was needed to understand regional aquifer systems and water available to wells and springs for rearing important Lake Michigan fish species at the Kettle Moraine Springs State Fish Hatchery in Sheboygan County, Wisconsin. As a basis for estimating the groundwater resources available, an existing groundwater-flow model was refined, and new groundwater-flow models were developed for the Kettle Moraine Springs State Fish Hatchery area using the U.S. Geological Survey (USGS) finite-difference code MODFLOW. This report describes the origin and construction of these groundwater-flow models and their use in testing conceptual models and simulating the hydrogeologic system. The study area is in the Eastern Ridges and Lowlands geographical province of Wisconsin, and the hatchery property is situated on the southeastern edge of the Kettle Moraine, a north-south trending topographic high of glacial origin. The bedrock units underlying the study area consist of Cambrian, Ordovician, and Silurian units of carbonate and siliciclastic lithology. In the Sheboygan County area, the sedimentary bedrock sequence reaches a thickness of as much as about 1,600 feet (ft). Two aquifer systems are present at the Kettle Moraine Springs State Fish Hatchery. A shallow system is made up of Silurian bedrock, consisting chiefly of dolomite, overlain by unconsolidated Quaternary-age glacial deposits. The glacial deposits of this aquifer system are the typical source of water to local springs, including the springs that have historically supplied the hatchery. The shallow aquifer system, therefore, consists of the unconsolidated glacial aquifer and the underlying bedrock Silurian aquifer. Most residential wells in the area draw from the Silurian aquifer. A deeper confined aquifer system is made up of Cambrian- and Ordovician-age bedrock units including sandstone formations. Because of its depth, very few wells are completed in the Cambrian-Ordovician aquifer system (COAS) near the Kettle Moraine Springs State Fish Hatchery. Three groundwater-flow models were used to estimate the water resources available to the hatchery from bedrock aquifers under selected scenarios of well placement and seasonal water requirements and subject to constraints on the effects of pumping on neighboring wells, local springs, and creeks. Model input data (recharge, water withdrawal, and boundary conditions) for these models were compiled from a number of data and information sources. The first model, named the “KMS model,” (KMS stands for Kettle Moraine Springs) is an inset model derived from a published USGS regional Lake Michigan Basin model and was constructed to simulate groundwater pumping from the semiconfined Silurian aquifer. The second model, named the “Pumping Test model,” was constructed to evaluate an aquifer pumping test conducted in the COAS as part of this project. The Pumping Test model was also used to simulate the local effects of 20 years of groundwater pumping from this deep bedrock aquifer for future hatchery operations. The third model, named the “LMB modified model,” is a version of the published Lake Michigan Basin (LMB) model that was modified with aquifer parameters refined in an area around the hatchery (approximately a 5-mile radius circle, corresponding to the area stressed by the aquifer pumping test). This LMB modified model was applied to evaluate regional effects of pumping from the confined COAS. The available Silurian aquifer groundwater resource was estimated using the KMS model with three scenarios—named “AllConstraints,” “Constraints2,” and “Constraints3”—that specified local water-level and flow constraints such as drawdown at nearby household wells, water levels inside pumping well boreholes, and flow in local streams and springs. Each scenario utilized the MODFLOW Groundwater Management Process (GWM) to select three locations from six candidate locations that provided the greatest combined flow while satisfying the constraints. The three constraint scenarios provided estimates of 430 gallons per minute (gal/min), 480 gal/min, and 520 gal/min pumping from three wells—AllConstraints, Constraints2, and Constraints3, respectively. The same three wells were selected for the scenarios that estimated 480 gal/min and 520 gal/min; the scenario that estimated 430 gal/min shared two of these same wells, but the third selected well was different. The available COAS groundwater resource was estimated by two scenarios with each conducted over a period of 20 years with the Pumping Test model and the LMB modified model . The Pumping Test model was used to simulate local effects of pumping, and the LMB modified model was used to simulate regional effects of pumping. The scenarios simulate a range of total and seasonal pumping rates potentially linked to site activities. Scenario 1 simulates two wells completed in the Cambrian-Ordovician aquifer system, each pumping for 8 months at 300 gal/min, followed by pumping for 4 months at 600 gal/min. The average yearly pumping rate of Scenario 1 is 800 gal/min. Scenario 2 simulates three wells completed in the Cambrian-Ordovician aquifer system pumping for 8 months at 200 gal/min, followed by pumping for 4 months at 500 gal/min. The average yearly pumping rate of Scenario 2 is 900 gal/min. The Pumping Test model simulations confirmed that drawdown in the boreholes of the pumping wells at the selected 2-well or 3-well rates will meet the desired condition that the pumping water level remains at least 100 ft above the highest Cambrian-Ordovician unit open to the well. The LMB modified model was used to evaluate the regional drawdown of the pumping from the confined COAS under the same 2-well and 3-well scenarios. At the nearest known existing COAS well, Campbellsport production well #4, the simulated drawdown for Scenario 1 after 20 years of cyclical pumping with two pumping wells averaging a total of 800 gal/min is 16.9 ft, whereas the simulated drawdown for Scenario 2 after 20 years of pumping with three pumping wells averaging a total of 900 gal/min is 19.0 ft. The total deep aquifer thickness at the Campbellsport location is on the order of 620 ft, meaning that the simulated drawdown for either scenario is about 3 percent of the confined aquifer thickness. The models developed as part of this project are archived in the project data release. The archive includes the model input and output files as well as MODFLOW source code and executables. ( Haserodt and others, 2017 ).

Wisconsin↗

USGS2010

Conceptual understanding and groundwater quality of selected basin-fill aquifers in the Southwestern United States

The National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey has been conducting a regional analysis of water quality in the principal aquifer systems in the southwestern United States (hereinafter, “Southwest”) since 2005. Part of the NAWQA Program, the objective of the Southwest Principal Aquifers (SWPA) study is to develop a better understanding of water quality in basin-fill aquifers in the region by synthesizing information from case studies of 15 basins into a common set of important natural and human-related factors found to affect groundwater quality. The synthesis consists of three major components: 1. Summary of current knowledge about the groundwater systems, and the status of, changes in, and influential factors affecting quality of groundwater in basin-fill aquifers in 15 basins previously studied by NAWQA (this report). 2. Development of a conceptual model of the primary natural and human-related factors commonly affecting groundwater quality, thereby building a regional understanding of the susceptibility and vulnerability of basin-fill aquifers to contaminants. 3. Development of statistical models that relate the concentration or occurrence of specific chemical constituents in groundwater to natural and human-related factors linked to the susceptibility and vulnerability of basin-fill aquifers to contamination. Basin-fill aquifers occur in about 200,000 mi 2 of the 410,000 mi 2 SWPA study area and are the primary source of groundwater supply for cities and agricultural communities. Four of the principal aquifers or aquifer systems of the United States are included in the basin-fill aquifers of the study area: (1) the Basin and Range basin-fill aquifers in California, Nevada, Utah, and Arizona; (2) the Rio Grande aquifer system in New Mexico and Colorado; (3) the California Coastal Basin aquifers; and (4) the Central Valley aquifer system in California. Because of the generally limited availability of surface-water supplies in the arid to semiarid climate, cultural and economic activities in the Southwest are particularly dependent on supplies of good-quality groundwater. Irrigation and public-supply withdrawals from basin-fill aquifers in the study area account for about one quarter of the total withdrawals from all aquifers in the United States. Many factors influence the quality of groundwater in the 15 case-study basins, but some common factors emerge from the basin summaries presented in this report. These factors include the chemical composition of the recharge water, consolidated rock geology and composition of aquifer materials derived from consolidated rock, and land and water use. The major water-quality issues in many of the developed case-study basins are increased concentrations of dissolved solids, nitrate, and VOCs in groundwater as a result of human activities. The information presented and the citations listed in this report serve as a resource for those interested in the groundwater-flow systems in the NAWQA case-study basins. The summaries of water-development history, hydrogeology, conceptual understanding of the groundwater system under both predevelopment and modern conditions, and effects of natural and human-related factors on groundwater quality presented in the sections on each basin also serve as a foundation for the synthesis and modeling phases of the SWPA regional study.

Arizona, California, Colorado, Nevada, New Mexico,↗

USGS2015

The quality of our Nation's waters: Groundwater quality in the Columbia Plateau and Snake River Plain basin-fill and basaltic-rock aquifers and the Hawaiian volcanic-rock aquifers, Washington, Idaho, and Hawaii, 1993-2005

The Columbia Plateau, Snake River Plain, and Hawaii are large volcanic areas in the western United States and mid-Pacific ocean that contain extensive regional aquifers of a hard, gray, volcanic rock called basalt. Residents of the Columbia Plateau, the Snake River Plain, and the island of Oahu depend on groundwater as their primary source of drinking water. Although the depth to the water table can be several hundred feet, the groundwater is highly vulnerable to contamination because the permeable sediments and rocks allow contaminants to move readily down to the water table. Intense agricultural and urban activities occur above the drinking-water supply and are increasing in some areas. Contaminants, such as nitrate, pesticides, and volatile organic compounds, associated with agricultural and urban activities, have adversely affected groundwater quality.

Hawaii, Idaho, Washington↗

USGS2015

Groundwater-level and storage-volume changes in the Equus Beds aquifer near Wichita, Kansas, predevelopment through January 2015

Development of the Wichita well field began in the 1940s in the Equus Beds aquifer to provide the city of Wichita, Kansas, a new water-supply source. After development of the Wichita well field began, groundwater levels began to decline. Extensive development of irrigation wells that began in the 1970s also contributed to substantial groundwater-level declines. Groundwater-level declines likely enhance movement of brine from past oil and gas production near Burrton, Kansas, and natural saline water from the Arkansas River into the Wichita well field. Groundwater levels reached a historical minimum in 1993 because of drought conditions, irrigation, and the city of Wichita’s withdrawals from the aquifer. In 1993, the city of Wichita adopted the Integrated Local Water Supply Program to ensure that Wichita’s water needs would be met through the year 2050 and beyond as part of its efforts to manage the part of the Equus Beds aquifer Wichita uses. A key component of the Integrated Local Water Supply Program was the Equus Beds Aquifer Storage and Recovery project. The Aquifer Storage and Recovery project’s goal is to store and eventually recover groundwater and help protect the Equus Beds aquifer from oil-field brine water near Burrton, Kansas, and saline water from the Arkansas River. Since 1940, the U.S. Geological Survey has monitored groundwater levels and storage-volume changes in the Equus Beds aquifer to provide data to the city of Wichita in order to better manage its water supply. Groundwater mostly flowed from west to east in the shallow and deep parts of the Equus Beds aquifer in January 2015. A large area of declines greater than 10 feet in the shallow part of the Equus Beds aquifer from predevelopment (before substantial pumpage began in the area in September 1940) to January 2015 covered most of the central part of the study area, where the city of Wichita well field is located, and extended beyond it. Groundwater-level rises of greater than 10 feet from 1993 (the historical minimum groundwater levels) to January 2015 covered most of the central part of the study area in the shallow and deep parts of the Equus Beds aquifer; rises of greater than 20 feet mostly were within the north-central part of the study area. The 1993 to January 2015 recovery of storage volume previously lost from predevelopment to 1993 was about 46 percent (55,200 acre-feet) for the central part of the study area and the percentage recovery was larger than the 31 percent (59,800 acre-feet) recovery for the entire study area. Groundwater-level rises and the larger percentage recovery of storage volume in the central part of the study area was most likely a result of the city of Wichita adopting the Integrated Local Water Supply Program strategy which reduced Wichita’s pumpage from the Equus Beds aquifer in 2014 to the smallest amount since 1940. January 2015 storage volumes were about 96 percent (3,057,000 acre-feet) and 94 percent (960,000 acre-feet) of total aquifer storage for the study area and the central part of the study area, respectively. Groundwater levels from January 2014 to January 2015 in the central part of the study area rose about 3 feet in some places, probably because Wichita reduced its withdrawals from the aquifer in 2014 by more than 50 percent. Groundwater levels probably recovered less than anticipated because of decreased recharge and net groundwater flow and increased agricultural pumpage. A volumetric water budget for the central part of the study area between 2013 and 2014 showed that the substantial decrease in total pumping (10,412 acre-feet) did not result in an increase in storage volume because it was more than offset by decreased recharge (6,502 acre-feet; artificial and from precipitation) and an even greater decrease in net groundwater flow (11,710 acre-feet).

Kansas↗

USGS2012

Groundwater flow and water budget in the surficial and Floridan aquifer systems in east-central Florida

A numerical transient model of the surficial and Floridan aquifer systems in east-central Florida was developed to (1) increase the understanding of water exchanges between the surficial and the Floridan aquifer systems, (2) assess the recharge rates to the surficial aquifer system from infiltration through the unsaturated zone and (3) obtain a simulation tool that could be used by water-resource managers to assess the impact of changes in groundwater withdrawals on spring flows and on the potentiometric surfaces of the hydrogeologic units composing the Floridan aquifer system. The hydrogeology of east-central Florida was evaluated and used to develop and calibrate the groundwater flow model, which simulates the regional fresh groundwater flow system. The U.S. Geological Survey three-dimensional groundwater flow model, MODFLOW-2005, was used to simulate transient groundwater flow in the surficial, intermediate, and Floridan aquifer systems from 1995 to 2006. The east-central Florida transient model encompasses an actively simulated area of about 9,000 square miles. Although the model includes surficial processes-rainfall, irrigation, evapotranspiration, runoff, infiltration, lake water levels, and stream water levels and flows-its primary purpose is to characterize and refine the understanding of groundwater flow in the Floridan aquifer system. Model-independent estimates of the partitioning of rainfall into evapotranspiration, streamflow, and aquifer recharge are provided from a water-budget analysis of the surficial aquifer system. The interaction of the groundwater flow system with the surface environment was simulated using the Green-Ampt infiltration method and the MODFLOW-2005 Unsaturated-Zone Flow, Lake, and Streamflow-Routing Packages. The model is intended to simulate the part of the groundwater system that contains freshwater. The bottom and lateral boundaries of the model were established at the estimated depths where the chloride concentration is 5,000 milligrams per liter in the Floridan aquifer system. Potential flow across the interface represented by this chloride concentration is simulated by the General Head Boundary Package. During 1995 through 2006, there were no major groundwater withdrawals near the freshwater and saline-water interface, making the general head boundary a suitable feature to estimate flow through the interface. The east-central Florida transient model was calibrated using the inverse parameter estimation code, PEST. Steady-state models for 1999 and 2003 were developed to estimate hydraulic conductivity (K) using average annual heads and spring flows as observations. The spatial variation of K was represented using zones of constant values in some layers, and pilot points in other layers. Estimated K values were within one order of magnitude of aquifer performance test data. A simulation of the final two years (2005-2006) of the 12-year model, with the K estimates from the steady-state calibration, was used to guide the estimation of specific yield and specific storage values. The final model yielded head and spring-flow residuals that met the calibration criteria for the 12-year transient simulation. The overall mean residual for heads, defining residual as simulated minus measured value, was -0.04 foot. The overall root-mean square residual for heads was less than 3.6 feet for each year in the 1995 to 2006 simulation period. The overall mean residual for spring flows was -0.3 cubic foot per second. The spatial distribution of head residuals was generally random, with some minor indications of bias. Simulated average evapotranspiration (ET) over the 1995 to 2006 period was 34.5 inches per year, compared to the calculated average ET rate of 36.6 inches per year from the model-independent water-budget analysis. Simulated average net recharge to the surficial aquifer system was 3.6 inches per year, compared with the calculated average of 3.2 inches per year from the model-independent waterbudget analysis. Groundwater withdrawals from the Floridan aquifer system averaged about 800 million gallons per day, which is equivalent to about 2 inches per year over the model area and slightly more than half of the simulated average net recharge to the surficial aquifer system over the same period. Annual net simulated recharge rates to the surficial aquifer system were less than the total groundwater withdrawals from the Floridan aquifer system only during the below-average rainfall years of 2000 and 2006.

Florida↗

USGS2012

Groundwater flow and water budget in the surficial and Floridan aquifer systems in east-central Florida

A numerical transient model of the surficial and Floridan aquifer systems in east-central Florida was developed to (1) increase the understanding of water exchanges between the surficial and the Floridan aquifer systems, (2) assess the recharge rates to the surficial aquifer system from infiltration through the unsaturated zone and (3) obtain a simulation tool that could be used by water-resource managers to assess the impact of changes in groundwater withdrawals on spring flows and on the potentiometric surfaces of the hydrogeologic units composing the Floridan aquifer system. The hydrogeology of east-central Florida was evaluated and used to develop and calibrate the groundwater flow model, which simulates the regional fresh groundwater flow system. The U.S. Geological Survey three-dimensional groundwater flow model, MODFLOW-2005, was used to simulate transient groundwater flow in the surficial, intermediate, and Floridan aquifer systems from 1995 to 2006. The East-Central Florida Transient model encompasses an actively simulated area of about 9,000 square miles. Although the model includes surficial processes-rainfall, irrigation, evapotranspiration (ET), runoff, infiltration, lake water levels, and stream water levels and flows-its primary purpose is to characterize and refine the understanding of groundwater flow in the Floridan aquifer system. Model-independent estimates of the partitioning of rainfall into ET, streamflow, and aquifer recharge are provided from a water-budget analysis of the surficial aquifer system. The interaction of the groundwater flow system with the surface environment was simulated using the Green-Ampt infiltration method and the MODFLOW-2005 Unsaturated-Zone Flow, Lake, and Streamflow-Routing Packages. The model is intended to simulate the part of the groundwater system that contains freshwater. The bottom and lateral boundaries of the model were established at the estimated depths where the chloride concentration is 5,000 milligrams per liter in the Floridan aquifer system. Potential flow across the interface represented by this chloride concentration is simulated by the General Head Boundary Package. During 1995 through 2006, there were no major groundwater withdrawals near the freshwater and saline-water interface, making the general head boundary a suitable feature to estimate flow through the interface. The east-central Florida transient model was calibrated using the inverse parameter estimation code, PEST. Steady-state models for 1999 and 2003 were developed to estimate hydraulic conductivity (K) using average annual heads and spring flows as observations. The spatial variation of K was represented using zones of constant values in some layers, and pilot points in other layers. Estimated K values were within one order of magnitude of aquifer performance test data. A simulation of the final two years (2005-2006) of the 12-year model, with the K estimates from the steady-state calibration, was used to guide the estimation of specific yield and specific storage values. The final model yielded head and spring-flow residuals that met the calibration criteria for the 12-year transient simulation. The overall mean residual for heads, defining residual as simulated minus measured value, was -0.04 foot. The overall root-mean square residual for heads was less than 3.6 feet for each year in the 1995 to 2006 simulation period. The overall mean residual for spring flows was -0.3 cubic foot per second. The spatial distribution of head residuals was generally random, with some minor indications of bias. Simulated average ET over the 1995 to 2006 period was 34.47 inches per year, compared to the calculated average ET rate of 36.39 inches per year from the model-independent water-budget analysis. Simulated average net recharge to the surficial aquifer system was 3.58 inches per year, compared with the calculated average of 3.39 inches per year from the model-independent water-budget analysis. Groundwater withdrawals from the Floridan aquifer system averaged about 920 million gallons per day, which is equivalent to about 2 inches per year over the model area and slightly more than half of the simulated average net recharge to the surficial aquifer system over the same period. Annual net simulated recharge rates to the surficial aquifer system were less than the total groundwater withdrawals from the Floridan aquifer system only during the below-average rainfall years of 2000 and 2006.

Florida↗

USGS2020

Hydrogeology and groundwater geochemistry of till confining units and confined aquifers in glacial deposits near Litchfield, Cromwell, Akeley, and Olivia, Minnesota, 2014–18

Confined (or buried) aquifers of glacial origin overlain by till confining units provide drinking water to hundreds of thousands of Minnesota residents. The sustainability of these groundwater resources is not well understood because hydraulic properties of till that control vertical groundwater fluxes (leakage) to underlying aquifers are largely unknown. The U.S. Geological Survey, Iowa State University, Minnesota Geological Survey, and Minnesota Department of Health investigated hydraulic properties and groundwater flow through till confining units using field studies and heuristic MODFLOW simulations. Till confining units in the following late-Wisconsinan stratigraphic units (with locations in parentheses) were characterized: Des Moines lobe till of the New Ulm Formation (Litchfield, Minnesota), Superior lobe till of the Cromwell and Aitkin Formations (Cromwell, Minn.), and Wadena lobe till of the Hewitt Formation (hydrogeology field camp [HFC] near Akeley, Minn.). Pre-Illinoian till of the Good Thunder formation (Olivia, Minn.) was also characterized. Hydraulic and geochemical field data were collected from sediment cores and a series of five piezometer nests. Each nest consisted of five to eight piezometers screened at short vertical intervals in hydrostratigraphic units including (if present) surficial aquifers, till confining units, confined/buried aquifers, and underlying bedrock. Till hydraulic conductivity was estimated from slug tests (horizontal [ K h ]) and constant-rate aquifer tests in the confined aquifer (vertical [ K v ]). Travel times through the till were evaluated with Darcy’s law and stable isotope concentrations. A series of heuristic MODFLOW simulations were used to evaluate groundwater fluxes through till across the range of till hydraulic properties and pumping rates observed at the field sites. The field data demonstrated variability in hydraulic properties between and within till stratigraphic units horizontally and vertically. The variability in hydraulic properties within and between sites resulted in substantial differences in groundwater flux through till. A conceptual understanding that emerges from the vertical till profiles is that they are not homogeneous hydrostratigraphic units with uniform properties; rather, each vertical sequence is a heterogeneous mixture of glacial sediment with differing abilities to transmit water. Till thicknesses varied from 60 to 166 feet, and till textures ranged from a sandy loam (Hewitt Formation, HFC site) to a silt loam/clay loam (Good Thunder formation, Olivia site). Till Kh varied by one to three orders of magnitude within each piezometer nest. Four piezometer nests had downward hydraulic gradients ranging from 0.04 to 0.56, and one nest had a slight upward hydraulic gradient of 0.02. The Cromwell, HFC, and Litchfield 1 sites were examples of “leaky” tills with high Kv (0.001 to 1.1 feet per day [ft/d]) and geometric mean Kh (0.03 to 0.07 ft/d) and extensive vertical hydraulic connectivity between the confined aquifer and the overlying till. Estimated groundwater travel times through these sites ranged from 1 to 81 years, and two of these sites had tritium throughout their till profiles. The tills at the other two sites, Olivia and Litchfield 2, were effective confining units that had low Kv (0.001 to 0.0005 ft/d) and geometric mean Kh (0.0002 to 0.004 ft/d). The till piezometers at these sites had no drawdown response to short-term (up to 10 hours for Olivia and up to 5 days for Litchfield) high-capacity pumping from the confined aquifer. Estimated groundwater travel times through the tills at these sites ranged from 165 to nearly 1,800 years, and tritium was only detected in the upper one-third of these till profiles. Across all sites, the till vertical anisotropy (ratio of K h to K v ) ranged by four orders of magnitude from 0.05 at the Cromwell nest to 70 at the Litchfield 1 nest. Stable isotopes of oxygen and hydrogen indicate that groundwater throughout all five till profiles is younger than the last glacial advance into Minnesota at about 11,000 years ago. The heuristic modeling demonstrated that, for understanding sustainability of groundwater pumping from confined aquifers, knowledge of till hydraulic properties is just as important as knowledge of aquifer hydraulic properties. Substantial differences in groundwater fluxes into and through till were observed across hydrogeologic settings representative of the field sites. Over long periods of time (hundreds of years), pumping-induced hydraulic gradients are established in confined aquifer systems and, even in low hydraulic conductivity tills, these pumping-induced hydraulic gradients increase leakage into and through till compared to ambient conditions. In conclusion, groundwater flowing vertically downward through till confining units (leakage) replenishes water pumped from confined aquifers. Till hydraulic properties, such as those presented in this report, provide important information that can be used to quantify leakage rates through till. Till hydraulic properties are variable over short distances and profoundly affect leakage rates, demonstrating the importance of site-specific till hydraulic data for evaluating the sustainability of groundwater withdrawals from confined aquifers.

Minnesota↗