Geology topics

William G. Eldridge

Publications and source records attributed to William G. Eldridge.

13 recordsLinked to original sources

USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗

USGS2025

Hydrologic budgets and water availability of six bedrock aquifers in the Black Hills area, South Dakota and Wyoming, 1931–2022

Population growth and recurring droughts in the Black Hills region raised interest in water resources and future availability. The Black Hills hydrology study (BHHS) was initiated in the early 1990s to address questions regarding water resources. Since completion of the BHHS in the early 2000s, the population of the Black Hills region increased by about 39 percent, which has renewed interest in water demand and availability in the Black Hills. The U.S. Geological Survey, in cooperation with the Western Dakota Regional Water System, completed a study to update hydrologic budgets from the BHHS for six of the most used aquifers in the Black Hills. Water availability was determined by comparing results from hydrologic budgets to modern well withdrawals (2003–22) and water rights information. Key updates to the BHHS budgets included adding available data from 1999 to 2022 and determining hydrologic budgets for six aquifers in nine smaller areas (called “subareas”). Inflows for the hydrologic budget included recharge from precipitation and streamflow losses to aquifers. Total mean annual recharge for the six aquifers in the study area was estimated at 278,900 acre-feet, with 205,100 acre-feet from precipitation recharge and 73,800 acre-feet from streamflow recharge. Mean annual precipitation recharge for the Madison and Minnelusa aquifers together accounted for 76 percent of the total mean annual precipitation recharge, with the Madison aquifer contributing 57,000 acre-feet and the Minnelusa aquifer contributing 98,100 acre-feet. Outflow components estimated for the hydrologic budget include artesian springflow and well withdrawals. Total mean annual artesian springflow in the study area was estimated as 166,100 acre-feet for the combined Madison and Minnelusa aquifers. Mean total annual well withdrawals for 2003–22 in the study area were about 50,000 acre-feet. No increased well withdrawal patterns corresponding to population increases were observed between 2003 and 2022. Water availability was determined by comparing total annual appropriations and mean and maximum annual well withdrawals for 2003–22 to mean annual recharge for 1931–2022 for each aquifer in subareas 1–9. Modern well withdrawals (mean and maximum for 2003–22) exceeded mean annual recharge for only the Deadwood and Inyan Kara aquifers in subareas 9 and 4, respectively. Additionally, total annual appropriations did not exceed mean annual recharge in most subareas, except most notably in subarea 4 (Rapid City area) where appropriations exceeded recharge for the Madison, Minnelusa, and Inyan Kara aquifers. Total annual appropriations also exceeded mean annual recharge for the Inyan Kara aquifer in subareas 3 and 5. In addition to recharge, water availability includes the water stored in pore spaces of aquifer materials. Estimates of total volume of recoverable water in storage were updated as part of this study to include the portion of aquifers in Wyoming, which were omitted during the BHHS. In total, the estimated total amount of recoverable water in storage in the study area was 356.9 million acre-feet for six major aquifers in the Black Hills area of South Dakota and Wyoming.

South Dakota, Wyoming↗

USGS2023

Maps of elevation of top of Pierre Shale and surficial deposit thickness with hydraulic properties from borehole geophysics and aquifers tests within and near Ellsworth Air Force Base, South Dakota, 2020–21

The U.S. Geological Survey, in cooperation with the U.S. Air Force Civil Engineer Center, collected borehole geophysical data and completed simple aquifer tests to estimate the thickness and hydraulic properties of surficial deposits. The purpose of data collection was to create generalized contour maps of Pierre Shale elevation and surficial deposit thickness within and near Ellsworth Air Force Base (study area). Natural gamma and electromagnetic induction data were collected to refine or determine surficial deposit thickness at selected wells. Additionally, data from previous geophysical studies and driller logs were compiled and combined with results from natural gamma and electromagnetic induction data to provide a more spatially complete image of the subsurface. Borehole nuclear magnetic resonance (bNMR) data were collected to estimate hydraulic conductivity and water content of surficial deposits overlying Pierre Shale. Simple aquifer tests using water slugs (slug tests) were completed to estimate hydraulic conductivity of surficial deposits, and results were compared to hydraulic conductivity estimates from bNMR data. All data used to construct maps and estimate hydraulic properties are provided in an accompanying U.S. Geological Survey data release (available at https://doi.org/10.5066/P9FLR79F ). Generalized contour maps were constructed using results from 26 borehole geophysical logs, 35 geophysical transects from previous studies, and 304 wells with driller logs. Pierre Shale elevation generally followed land-surface topography, sloping from high elevations in the north to lower elevations in the south. Topographic highs of Pierre Shale, where present, could act as groundwater divides that potentially affect groundwater flow direction. Surficial deposit thickness varied spatially and ranged from 0 to 86 feet. Surficial deposits generally were thickest in higher elevation areas near ephemeral streams in the northern part of the study area. Hydraulic conductivity estimated from bNMR results using two analytical methods ranged from 0.1 to 2,314 feet per day, whereas hydraulic conductivity estimated from slug tests ranged from 0.001 to 193 feet per day. Hydraulic conductivity estimates from slug tests were plotted with surficial deposit thickness contours instead of bNMR estimates because bNMR estimates were determined to overestimate hydraulic conductivity. Hydraulic conductivity values generally were greater in the southwestern part of study area than the northeastern part.

South Dakota↗

USGS2021

Borehole analysis, single-well aquifer testing, and water quality for the Burnpit well, Mount Rushmore National Memorial, South Dakota

Mount Rushmore National Memorial (hereafter referred to as “the memorial”), in western South Dakota, is maintained by the National Park Service (NPS) and includes 1,278 acres of land in the east-central part of the Black Hills. An ongoing challenge for NPS managers at the memorial is providing water from sustainable and reliable sources for operations, staff, and the increasing number of visitors. In 2020, the U.S. Geological Survey (USGS) and NPS completed a hydrological study of the Burnpit well (well 5), a 580-foot-deep open hole groundwater well completed in metamorphic (crystalline) rock at the memorial. The purpose of this study was to estimate the geological and hydraulic properties of the aquifer supplying the well and to determine the water quality of the groundwater from the well. The study provides NPS staff and managers background information for assessing future uses for the well. Methods for data collection and analysis for the study included borehole and video camera analysis in 2020, aquifer testing by the NPS in 2009 and the USGS in 2020, and water-quality sampling in 2020. Borehole camera video generally matched the lithology recorded in the well log. Fractures recorded in the well log and observed with the borehole camera, including more than 20 less prominent fractures and rough sidewall areas, indicated a fractured aquifer. The fractures are the primary conduits for groundwater flow through the rock and into the well. Transmissivity was estimated for the upper and lower water-level drawdown zones at the Burnpit well with data from the NPS and USGS using the Theis and Cooper-Jacob methods. Transmissivity for the NPS test using the Theis method was 9.0 and 11 feet squared per day (ft 2 /d) for the upper and lower drawdown zones, respectively. Using the Cooper-Jacob method, the transmissivity was 22 and 14 ft 2 /d for the upper and lower drawdown zones of the aquifer, respectively. Transmissivity estimates from data from the USGS test were similar. The Theis method, applied to the upper and lower drawdown zones of the aquifer, produced transmissivity estimates of 7.7 and 10 ft 2 /d, and the Cooper-Jacob method produced estimates of 9.7 and 12 ft 2 /d, respectively. Storativity (specific yield) estimated using the Theis method for the NPS aquifer-test data was 0.85 and 0.92 for the upper and lower drawdown zones of the aquifer, respectively. The Cooper-Jacob method applied to the NPS aquifer-test data produced storativity estimates of 0.11 and 0.50 for the upper and lower drawdown zones, respectively. The Theis method applied to the USGS aquifer-test data estimated storativity values of 0.77 and 1.0 for the upper and lower drawdown zones, respectively. The Cooper-Jacob method estimated storativity of 0.50 and 0.60 for the upper and lower drawdown zones of the USGS aquifer test, respectively. The estimated storativity values from the NPS and USGS aquifer tests for the upper and lower drawdown zones were higher than expected for limestones and schists. The hypothetical equilibrium drawdown for the Burnpit well was estimated after the NPS test in 2009 at no more, and possibly less, than 35 gallons per minute. The NPS noted that the sustainable yield likely was overestimated because the water level did not stabilize during the NPS aquifer test. The specific capacity for the NPS aquifer test in 2009 was 0.16 gallon per minute per foot ([gal/min]/ft) of drawdown at 3 hours, and the specific capacity for the USGS aquifer test in 2020 was 0.13 (gal/min)/ft of drawdown at 3 hours. The rate of water-level recovery after pumping ceased was 0.017 and 0.013 (gal/min)/ft for the NPS and USGS aquifer tests, respectively. The water-level recovery rate was nearly an order of magnitude less than the specific capacity estimated during pumping, indicating that water levels in the Burnpit well may not recover quickly enough during pumping to provide for a continuous source of water. Water-quality samples were collected at the Burnpit well on June 24 and July 23, 2020, and analyzed for field-measured properties, major ions, metals, nutrients, and perchlorate. Iron, zinc, and lithium concentrations for unfiltered samples in the well were at least three times greater than the mean filtered sample concentrations reported for crystalline aquifers in the Black Hills. Manganese concentrations were less than the mean concentration for crystalline aquifers but exceeded the U.S. Environmental Protection Agency (EPA) secondary drinking-water standards. The iron concentration from the June 24 sample was about 11 times greater than the EPA secondary drinking-water standards and mean concentrations from crystalline aquifers in the Black Hills. Arsenic concentrations in Burnpit well samples collected in 2020 were greater than the EPA primary drinking-water standard and the mean concentration for crystalline aquifers in the Black Hills. Arsenic occurs naturally in the rock of crystalline aquifers, and concentrations from samples in the Black Hills commonly exceed the EPA primary drinking-water standard of 10 micrograms per liter. High concentrations of arsenic, iron, and manganese metals in the Burnpit well make groundwater from the well in its natural state unusable as a drinking-water source, and water treatment would be necessary to reduce the trace element concentrations to less than the EPA primary and secondary drinking-water standards. However, if the memorial has immediate nonpotable water requirements, such as for construction and fire suppression, groundwater from the Burnpit well could provide water without causing additional stress to current (2021) drinking-water sources.

South Dakota↗

USGS2021

Spring types and contributing aquifers from water-chemistry and multivariate statistical analyses for seeps and springs in Theodore Roosevelt National Park, North Dakota, 2018

Water resources in Theodore Roosevelt National Park, North Dakota, support wildlife, visitors, and staff, and play a vital role in supporting the native ecology of the park. The U.S. Geological Survey, in cooperation with the National Park Service, completed field work in 2018 for a study to address concerns about water availability and possible sources of groundwater contamination for seeps and springs in Theodore Roosevelt National Park. The objective of the study was to improve hydrologic knowledge and determine the water composition of 11 seeps and springs in the park by collecting water-chemistry data at springs, streams, wells, and rain collectors. Water samples were collected at 26 sites at springs, streams, wells, and rain collectors in the North and South Units of Theodore Roosevelt National Park. Samples in the North Unit were collected at 5 springs, 1 stream, 2 wells, and 1 rain collector. Samples in the South Unit were collected at 6 springs, 2 streams, 8 wells, and 1 rain collector. Samples from springs, streams, and wells were collected in May, July, and September 2018. Samples from rain collectors were collected when enough daily precipitation accumulated in the collectors. Sampled precipitation events during the study period were in May, June, July, August, and September 2018. Physical properties of sampled water—temperature, pH, and specific conductance—were measured in the field. Water samples were analyzed for stable isotopes of oxygen and hydrogen and for chloride concentration. Recharge rates for aquifers supplying springs were determined using precipitation volume and chloride concentrations for a 12-day period before the sample-collection date. Multivariate statistical analysis methods used on water-chemistry data included principal component analysis, cluster analysis, and end-member mixing analysis. Water composition was used to determine the spring type and contributing aquifers for 11 springs in the North and South Units of Theodore Roosevelt National Park from analyses of water-chemistry data between May and September 2018. In the North Unit, Achenbach Spring was classified as a filtration spring with water from an unconfined part of the upper Fort Union aquifer and infiltration of precipitation. Hagen Spring, Mandal Spring, and Stevens Spring were classified as contact springs supplied by semiconfined parts of the upper Fort Union aquifer. Overlook Spring at one time may have been a natural spring or seep but now is a developed spring that behaves like a flowing artesian well completed in a confined part of the upper Fort Union aquifer. In the South Unit, six springs were classified into two spring types: filtration and contact springs. Boicourt Spring and Sheep Butte Spring were classified as filtration springs that have water supplied by unconfined parts of the upper Fort Union aquifer and infiltrated precipitation. Big Plateau Spring, Lone Tree Spring, Sheep Pasture Spring, and Southeast Corner Spring were classified as contact springs that receive waters from a semiconfined part of the upper Fort Union aquifer.

North Dakota↗

USGS2020

Hydrogeology and groundwater flow in alluvial deposits, north Summerset, South Dakota

The city of Summerset is a growing community in west South Dakota. The Sun Valley Estates subdivision in the north part of the city was developed on unconsolidated deposits surrounded by steep terrain. During years with greater than normal precipitation, particularly in 2019, groundwater levels increased in the unconsolidated deposits and caused damage to stormwater systems, sewer infrastructure, and houses with basements. The U.S. Geological Survey, in cooperation with the City of Summerset, completed a study of the hydrogeology and groundwater flow in the alluvial aquifer part of the unconsolidated deposits in north Summerset to understand the groundwater system in the area and to provide hydrogeologic information in support of future development planning. The study area included most of the Sun Valley Estates subdivision in the north part of the city of Summerset in the east Black Hills of west South Dakota. About 0.7 square mile of water-bearing alluvial deposits is included in the study area. Precipitation in the study area from 2017 to 2019 was compared to the monthly normal values at a nearby climate site. The largest departure from normal was in May 2019 with precipitation exceeding the monthly normal by about 5 inches (in.). All months in 2019, except March, exceeded the monthly normal precipitation. Cumulative departure from normal precipitation in 2019 increased from about 4 in. greater than normal in January to about 18 in. greater than normal in December. The geologic setting of the study area is characterized by the surrounding Black Hills. Unconsolidated Quaternary-age deposits overlie consolidated to partially consolidated Mesozoic-age and Paleozoic-age shales, sandstones, and limestones. Surficial deposits of alluvium and other unconsolidated deposits are the primary surficial geologic units in the study area and form the components of the alluvium hydrogeologic unit of the study area. Results from previous studies of alluvium along nearby Rapid Creek estimated hydraulic conductivity to range from 89 to 2,292 feet per day (ft/d), transmissivity to range from 1,001 to 32,083 feet squared per day, and storage coefficients to range from 0.0002 to 0.16. Hydraulic conductivity and transmissivity generally decreased downstream along Rapid Creek (west to east). Slug tests were completed August 16, 2019, at two observation wells completed in the alluvial aquifer in the Sun Valley Estates subdivision to determine hydraulic conductivity. Hydraulic conductivity estimated from AQTESOLV curve-fitting analysis using the Bouwer-Rice method for all slug-in and slug-out trials from two observation wells in the study ranged from 0.20 to 0.26 ft/d for well 441318103220001 (SunValley1 well) and from 0.54 to 14 ft/d for well 441319103215701 (SunValley2 well). The mean, median, and standard deviation of all trials at both wells were 4.3 ft/d, 0.8 ft/d, and 5.6 ft/d, respectively. The extent of the alluvial aquifer was determined by geologic maps and lithologic logs. Alluvial deposits in the study area extend to about 1 mile in the north–south direction and about 1.5 miles in the southeast–northwest direction. The direction of groundwater flow was estimated using water-level records and topographic maps. The resulting potentiometric map indicated that groundwater in the alluvial aquifer under the Sun Valley Estates subdivision originates from higher elevations of the west part of the area of interest and from streams in the southeast part. Recharge and evapotranspiration estimates were results from a Soil-Water Balance model that calculated a matrix of recharge for 2019 with values ranging from 0 to 11.4 in. and an annual mean value of 5.1 in. across the study area. Soil-Water Balance-estimated potential evapotranspiration for 2019 ranged from 28.90 to 28.75 in. and the estimated annual mean was 28.86 in. across the study area. Estimated groundwater budget components for the alluvial aquifer in the area of interest included inflows and outflows. Total estimated groundwater budget components for inflows for 2019 were about 66 percent from recharge, 33 percent from streamflow, and 1 percent from inflow from adjacent aquifers. Total estimated outflows were about 99-percent evapotranspiration and less than 1-percent outflow to adjacent aquifers.

South Dakota↗

USGS2020

Stratigraphic units of shallow unconsolidated deposits in Deadwood, South Dakota, delineated by real-time kinematic surveys

The City of Deadwood, South Dakota, has been working on a new archeological investigation in preparation for economic growth and expansion within the city limits, through the Deadwood Historic Preservation Office. During the excavation process, buried artifacts and historical features from the late 1800s have been uncovered. The stratigraphy of shallow unconsolidated deposits in the city of Deadwood, S. Dak., was surveyed on January 29, 2020, using real-time kinematic survey methods and described to identify variations in geologic material, thickness, and depth from the land surface in support of archeological studies by the city. The findings of the study will provide city managers and the public with reliable and impartial information for their use by advancing field or analytical methodology and understanding of hydrologic processes in the study area. The primary excavation site was surveyed, and stratigraphic units were delineated from changes in material properties or depositional environment. The primary excavation site consisted of nine stratigraphic units; however, some units were not consistent along the length of the excavation and pinched out along the cross section. Survey data points also were collected for artifacts and other sites of interest. The shallow surficial geology in the study area was affected by human construction, fires, and flooding.

South Dakota↗

USGS2020

Inventory and analysis of groundwater resources: Theodore Roosevelt National Park, North Dakota

Industrial and commercial developments in western North Dakota potentially could affect the sources of water that contribute to wells, spring flow, and seeps within Theodore Roosevelt National Park. Without basic water resources data, accurately predicting the effects of water withdrawals and water quality concerns related to industrial and commercial developments near the park would be challenging. Water resources in the park include surface water and groundwater. The Little Missouri River and its tributaries cross all three Theodore Roosevelt National Park units and are the primary surface-water features in the park. Groundwater resources include well discharges, springs, and seeps. The geology and hydrogeology of Theodore Roosevelt National Park are defined by the surrounding Williston Basin. Four aquifers are sources of groundwater to the park: unconsolidated aquifers including alluvial systems, the upper Fort Union aquifer, the lower Fort Union aquifer, and the Fox Hills-lower Hell Creek aquifer. Data used for wells, springs, seeps, and water quality in this report were compiled from the U.S. Geological Survey National Water Information System or from the North Dakota State Water Commission. An inventory of 16 wells was completed for sites within the boundaries of the park. In addition to well data, an inventory of 11 springs and seeps was completed. The groundwater-quality analysis had two objectives: (1) to characterize the groundwater chemistry in aquifers underlying the park and (2) to spatially map selected physical properties and chemical constituents of interest. Groundwater-quality data from the North Dakota State Water Commission were summarized, mapped, and used to characterize groundwater for each aquifer in the study area. Spatial concentration distribution maps were constructed for selected physical properties and chemical constituents using summary statistics and exceedances. Piper diagrams were used to classify and characterize groundwater for each aquifer. Future research to help fill data gaps in water resources information for Theodore Roosevelt National Park, including recommendations from previous studies, consists of the following: (1) evaluating the variability in discharge from springs and seeps in comparison to changes in precipitation or other recharge sources, (2) evaluating flow control measures for flowing artesian wells, (3) completing a water rights review, and (4) performing routine water-quality monitoring for wells and springs.

North Dakota↗

USGS2019

Groundwater-flow model and analysis of groundwater and surface-water interactions for the Big Sioux aquifer, Sioux Falls, South Dakota

The city of Sioux Falls, in southeastern South Dakota, is the largest city in South Dakota. The U.S. Geological Survey (USGS), in cooperation with the city of Sioux Falls, completed a groundwater-flow model to use for improving the understanding of groundwater-flow processes, estimating hydrogeologic properties, and analyzing groundwater and surface-water interactions for the Big Sioux aquifer in the model area. The model area includes the Big Sioux aquifer and the underlying hydrogeologic units from Dell Rapids, South Dakota, to the confluence of the Big Sioux River and the outlet of the Sioux Falls Diversion Channel in eastern Sioux Falls, S. Dak. The Big Sioux aquifer is the primary aquifer in the model area and the focus of the groundwater-flow model. The Big Sioux River is the largest stream in the model area and is in hydraulic connection with the Big Sioux aquifer. A conceptual model for the area was constructed and includes a characterization of the hydrogeologic framework, analysis and construction of potentiometric surfaces, and summary of estimated water budget components in the model area. The primary hydrogeologic units in the model area consist of (1) the Big Sioux aquifer, (2) a glacial till confining unit, and (3) bedrock aquifers (Split Rock Creek and Sioux Quartzite aquifers). Sources of groundwater recharge included infiltration of precipitation, stream seepage, and groundwater exchanges among the hydraulically connected Big Sioux aquifer, glacial till confining unit, and bedrock aquifers. Groundwater losses included evapotranspiration, groundwater discharge to streams, and groundwater withdrawal to supply water-use needs. A numerical groundwater-flow model (numerical model) was constructed and was used to simulate all aspects of the conceptual model for predevelopment (steady-state) and time-varying (transient) monthly conditions for 1950–2017. The numerical model was constructed using the USGS modular hydrologic simulation program, MODFLOW–6, and was calibrated using the Parameter ESTimation software, PEST++. The transient numerical model was calibrated for steady-state and transient monthly conditions for 1950–2017. Calibration targets were observations of hydraulic head, changes in hydraulic head, monthly mean streamflow (as a rate), and cumulative monthly stream discharge (as a volume). Parameters adjusted during model calibration were horizontal and vertical hydraulic conductivity, specific storage, specific yield, recharge and evapotranspiration multipliers, and streambed hydraulic conductivity. Horizontal and vertical hydraulic conductivity were estimated at pilot points distributed within the model area; specific storage and specific yield were assigned to uniform values in each layer in the model area; recharge and evapotranspiration multipliers were assigned uniformly for every stress period in the numerical model; and streambed hydraulic conductivity values were assigned uniformly between stream confluences. The final calibrated parameter values of horizontal and vertical hydraulic conductivity, specific yield, specific storage, streambed hydraulic conductivity, recharge, and evapotranspiration were considered reasonable for the hydrogeologic materials and conditions in the model area for 1950–2017. Overall, simulated hydraulic head altitudes had a linear regression coefficient of determination (R 2 ) of 0.48. Hydraulic head altitude residuals for the glacial till confining unit and bedrock aquifers were typically greater in magnitude when compared to residuals in the Big Sioux aquifer, but simulated hydraulic head altitudes in the Big Sioux aquifer compared favorably with mean observed hydraulic head altitudes and had a linear regression R 2 of 0.93. Simulated streamflow hydrographs matched the general trends of observed increases and decreases in streamflow for USGS streamgages 06482000 (Big Sioux River at Sioux Falls, S. Dak.) and 06482020 (Big Sioux River at North Cliff Avenue at Sioux Falls, S. Dak.), but larger streamflows were overestimated at the first streamgage and underestimated at the second streamgage. The numerical model reasonably estimated cumulative monthly stream discharge for the first 10–15 years of available streamflow records at both USGS streamgages. After the first 10–15 years of available streamflow record, cumulative monthly stream discharge was closely estimated for USGS streamgage 06482000 and underestimated at USGS streamgage 06482020. Composite sensitivities without regularization were calculated by PEST++ for the calibrated numerical model parameters and were averaged by parameter group. The parameter group with the highest mean composite sensitivity was the recharge multiplier parameter group. Model simplifications, assumptions, and limitations were necessary for construction of the conceptual and numerical models and for calibration efficiency. Spatial simplification of hydraulic properties could cause the numerical model to misrepresent reactions to changes in localized stresses, such as additional demands for groundwater withdrawal. The numerical model was temporally discretized into monthly periods and required scaling daily rates into representative monthly rates for model input and calibration targets. Based on the comparison between the observed and simulated groundwater levels, monthly mean streamflow and cumulative monthly stream discharge, and general groundwater distribution and flow, the numerical model favorably simulated the flow in the Big Sioux aquifer. Eventual capture was calculated in the model area using a steady-state numerical groundwater-flow model. The eventual capture map shows areas of higher streamflow capture adjacent to the Big Sioux River north of the city of Sioux Falls and along the lower part of the Sioux Falls Diversion Channel, and areas of lower streamflow capture along aquifer boundaries and near the southern Sioux Quartzite barrier. The timing of capture was determined using a transient numerical groundwater-flow model to determine the likely captured water sources for 30 years of groundwater withdrawal at three hypothetical wells using three continuous withdrawal rates (112.5, 450.0, and 900.0 gallons per minute). Supply for all three hypothetical wells became capture-dominated after only a short period of continuous withdrawal. Capture stabilized after about 10–15 years for well A, and after 20–25 years for well B, and after about 10–15 years for well C. The groundwater-flow model is a suitable tool to use for improving the understanding of groundwater-flow processes, estimating hydrogeologic properties, and analyzing groundwater and surface-water interactions for the Big Sioux aquifer near Sioux Falls, S. Dak. The numerical model can be used to simulate hydrologic scenarios, advance understanding of groundwater budgets, compute system response to stress, and determine likely sources of water supplied to wells.

South Dakota↗

USGS2019

Generalized potentiometric-surface map and groundwater flow directions in the Madison aquifer near Jewel Cave National Monument, South Dakota

A generalized potentiometric-surface map of the Madison aquifer near Jewel Cave National Monument was constructed using water levels measured from calendar years 1988 to 2019 in 24 groundwater wells and 4 subterranean cave lakes interpreted to be in hydraulic connection with the aquifer. The map indicated that groundwater near Jewel Cave National Monument originates from recharge sources to the Madison aquifer in the higher elevations in the north-central area of the map, flows west to south-southwest through the Jewel Cave network, then southeast. Hydrographs were constructed using water levels from four observation wells and one subterranean lake (Hourglass Lake) in the Jewel Cave network to evaluate historical and current groundwater recharge to the Madison aquifer in the study area. Hydrographs from 1992 through 2018 indicated water levels were lowest from the early to mid-1990s, increased through the late 1990s, peaked in the early 2000s, decreased until 2010, and then increased to the highest levels during 2016–18. A visual comparison of the Hourglass Lake hydrograph with cumulative precipitation, and quantitative (statistical) comparison of lake-water levels with cumulative precipitation, indicated that lake-water levels increased as cumulative precipitation increased, most likely due to some degree of precipitation recharge to hydraulically connected Madison Limestone outcrops. Comparing the potentiometric-surface map constructed for this study, with a map by Strobel and others (2000) of the same region and aquifer, indicated similarity and, therefore, provided some validation of map construction. The potentiometric-surface map constructed for this study could be used by park managers and others as a tool to evaluate the hydrogeologic characteristics of the Madison aquifer in the study area.

South Dakota↗

USGS2019

Hydraulic conductivity estimates from slug tests in the Big Sioux aquifer near Sioux Falls, South Dakota

Hydraulic conductivity estimates were made for 15 observation wells using slug-out (rising-head) tests in the Big Sioux aquifer near Sioux Falls, South Dakota, as part of a cooperative study with the City of Sioux Falls to characterize the hydrogeology and the extent of the Big Sioux aquifer north of the city. Well and aquifer data were collected from field measurements and drillers’ logs. Multiple slug tests were completed at each observation well with a transducer to record the change in water level and a U.S. Geological Survey standard mechanical slug to displace the well’s water column. In total, 110 slug-out test trials were completed among the 15 observation wells. Hydraulic conductivity was estimated by curve fitting with AQTESOLV Pro version 4.50.002. Hydraulic conductivity estimates ranged from 64 to 379 feet per day (ft/d). The mean, standard deviation, and median hydraulic conductivity for the 110 slug-out test trials were 171 ft/d, 73 ft/d, and 157 ft/d, respectively. The mean hydraulic conductivity calculated for each well ranged from 88 to 270 ft/d, the standard deviation ranged from 7 to 66 ft/d, and the median hydraulic conductivity ranged from 86 to 256 ft/d.

South Dakota↗

USGS2019

Multivariate analysis of hydrochemical data for Jewel Cave, Wind Cave, and surrounding areas

Jewel Cave National Monument and Wind Cave National Park in South Dakota contain two of the six longest caves worldwide. These caves contain subterranean lakes that are important points of intersection between the water table of the Madison aquifer and the caves. During 2015 to 2017, several subterranean lakes were discovered in Jewel Cave, which previously was thought to be above the regional water table. The objectives of this study were to better understand the hydrology of the recently discovered lakes in Jewel Cave and to evaluate their relation or possible connection to similar subterranean lakes in Wind Cave. Both objectives align with National Park Service resource management purposes. Multivariate analysis, consisting of principal component analysis (PCA), cluster analysis, and end member mixing, was applied to hydrochemical data for 70 sites within and surrounding Jewel Cave and Wind Cave. Hydrochemical data consisted of the contents of major ions (Ca, Mg, Na, HCO3, Cl, Si, SO4), arsenic (As), strontium (Sr), uranium (U), stable isotopes of oxygen and hydrogen (δ18O, δ2H), radiogenic isotope ratios of strontium and uranium (87Sr/86Sr and 234U/238U), pH, and electrical conductivity (EC) in water samples. Five hydrogeologic domains were identified on the basis of PCA and cluster analysis for the area encompassing Jewel Cave and Wind Cave in the southern Black Hills. The Artesian 1 and Artesian 2 domains represent primarily artesian springs and surrounding groundwater, the East and West domains represent areas where Madison and Minnelusa aquifer rocks are at or near the land surface, and the Precambrian domain represents the Precambrian aquifer. Multivariate analysis indicates that the Jewel Cave area is part of the West domain and that Wind Cave is part of the East domain. End member mixing was applied to estimate that groundwater in the Jewel Cave area primarily was derived from the West domain and secondarily from the Precambrian domain. Jewel Cave and Wind Cave contain lakes that are well connected to regional groundwater flow in the Madison aquifer.

South Dakota↗

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↗