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

Kyle W. Davis; William G. Eldridge; Joshua F. Valder; Kristen J. Valseth

doi:10.3133/sir20195117

Home/Archive/USGS

USGS · sir20195117

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

Kyle W. Davis·William G. Eldridge·Joshua F. Valder·Kristen J. Valseth

View original source ↗Read source document ↗Export metadata ↓

Abstract

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.

Explore related subjects

Scientific Investigations Report Report

Source-reported bounding extent: 43.29919735147067° to 43.757208878849376° latitude; -97.06146240234375° to -96.42425537109375° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Kyle W. Davis, William G. Eldridge, Joshua F. Valder, Kristen J. Valseth. 2019-11-27. Groundwater-flow model and analysis of groundwater and surface-water interactions for the Big Sioux aquifer, Sioux Falls, South Dakota. https://doi.org/10.3133/sir20195117

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2026

Monitoring recreation on federally managed lands and waters—Aspects of visitor use

Federally managed public lands and waters receive around 1 billion recreational visits each year. Data on these visitors can aid in guiding policy decisions, managing resources effectively, and communicating the economic contributions of lands and waters. This report explores the methods used by agencies to collect data on aspects of recreational visitor use to Federal lands and waters (apart from visitation numbers, which are the focus of a companion publication). Aspects of recreational visitor use include visitor demographics, recreational activity participation, visitor satisfaction, visitor attitudes and experiences, trip characteristics, and economic contributions. We review practices used to understand aspects of visitor use across seven Federal agencies, revealing similarities such as the use of visitor intercept surveys and coverage of similar topic areas, and differences in how survey programs are operationalized and how specific questions on visitor surveys are worded. We also evaluate emerging technologies, such as geolocated social media and mobile device location data, for their potential to aid in understanding aspects of visitor use. This report concludes with potential opportunities to enhance data collection and coordination, ensuring cost-effective data collection and informed decision making.

Scientific Investigations Report↗

USGS2026

Evaluation of best management practices at an edge-of-field site in the Eagle Creek watershed, Ohio, 2012–20

In 2010, the U.S. Geological Survey worked in partnership with the Great Lakes Restoration Initiative and the Natural Resources Conservation Service to identify farm fields in priority watersheds—watersheds critically important to Great Lake health. Intensive best management practices were implemented in five States to test the efficacy of best management practices, also referred to as agricultural conservation practices, to reduce sediment and nutrient runoff on agricultural fields. These fields were chosen as representatives of the priority watersheds because their farming practices and geographical conditions were common among farms in those watersheds. A farm site located near Findlay, Ohio, in the Eagle Creek watershed used two conservation practices, cover crop (Natural Resources Conservation Service practice standard 340) and variable rate technology (VRT; Natural Resources Conservation Service practice standard 590), which represent common practices in the Maumee priority watershed. This study monitored surface runoff and subsurface tile runoff at the Eagle Creek watershed site. The effects of cover crop and VRT nutrient application on sediment and nutrient runoff were assessed from October 2012 to September 2020. Cover crops were applied in fall of 2016, 2017, and 2018, coinciding with VRT applications. The following parameters were analyzed as part of the study: runoff discharge; peak discharge; and concentrations and loads of total phosphorus, particulate phosphorus, dissolved reactive phosphorus, and suspended sediment. High peak runoff events disproportionally affected nutrient concentrations and loads in runoff. A threshold of the 85th percentile of the peak discharge was chosen to assess the effect of cover crop and VRT nutrient application on nutrient mitigation at each gage. Runoff events below the 85th percentile were considered small, and runoff events over the 85th percentile were considered large. Concentrations from collected water-quality samples were used for the analysis; estimated concentrations were excluded from the statistical evaluation. All surface runoff parameters during small runoff events, except surface runoff volume, were significantly lower in the period with cover crops and VRT implementation. Contrarily, subsurface runoff parameters during small events were not statistically different between the periods with and without cover crops and VRT implementation. Surface and subsurface parameters during large events were also not statistically different between the periods with and without cover crops and VRT implementation. These results indicate cover crops and VRT may improve water quality during small runoff events, but additional best management practices that mitigate large runoff events may lead to greater water-quality improvement given the contribution of large events to overall losses.

Ohio↗

USGS2026

Continuous monitoring and temporal variability of fluorescence of dissolved organic matter in karst groundwater of the Edwards aquifer, south-central Texas, 2019–24

Widespread urbanization on the Edwards aquifer recharge zone in south-central Texas has prompted the development of new approaches and tools for evaluating the current and future status of water quality in the San Antonio and Barton Springs segments of the Edwards aquifer. The U.S. Geological Survey, in cooperation with the San Antonio Water System and the City of Austin, applied continuous monitoring of fluorescence of dissolved organic matter (fDOM) to characterize the sources and transport of dissolved organic matter and associated constituents in the Edwards aquifer. Continuous fDOM time-series data were adjusted using empirically derived temperature, turbidity, and inner-filter-effect correction factors prior to analysis and interpretation. Fully corrected fDOM time-series data from the study sites were evaluated for short- and long-term temporal variability in fDOM in the context of the hydrologic and climatic conditions that occurred during the 5-year study period to understand vulnerability of the aquifer to potential contaminants in recharge. Study results indicate that fDOM is an effective proxy for dissolved organic carbon that, in turn, is indicative of the influx of recent surface water and associated contaminants, and thus, aquifer vulnerability to those contaminants. Although other continuously monitored water-quality parameters provide insights into changes in water quality in response to varying hydrologic conditions, the quantitative relation between fDOM and dissolved organic carbon directly indicates the timing and magnitude of pulses of organic constituents derived from the land surface, whereas other water-quality parameters do not.

Texas↗