Geology ReportsSearch

Geology topics

Joshua F. Valder

Publications and source records attributed to Joshua F. Valder.

At least 19 recordsLinked to original sources

Tapwater exposures, effects potential, and residential risk management in Northern Plains Nations

In the United States (US), private-supply tapwater (TW) is rarely monitored. This data gap undermines individual/community risk-management decision-making, leading to an increased probability of unrecognized contaminant exposures in rural and remote locations that rely on private wells. We assessed point-of-use (POU) TW in three northern plains Tribal Nations, where ongoing TW arsenic (As) interventions include expansion of small community water systems and POU adsorptive-media treatment for Strong Heart Water Study participants. Samples from 34 private-well and 22 public-supply sites were analyzed for 476 organics, 34 inorganics, and 3 in vitro bioactivities. 63 organics and 30 inorganics were detected. Arsenic, uranium (U), and lead (Pb) were detected in 54%, 43%, and 20% of samples, respectively. Concentrations equivalent to public-supply maximum contaminant level(s) (MCL) were exceeded only in untreated private-well samples (As 47%, U 3%). Precautionary health-based screening levels were exceeded frequently, due to inorganics in private supplies and chlorine-based disinfection byproducts in public supplies. The results indicate that simultaneous exposures to co-occurring TW contaminants are common, warranting consideration of expanded source, point-of-entry, or POU treatment(s). This study illustrates the importance of increased monitoring of private-well TW, employing a broad, environmentally informative analytical scope, to reduce the risks of unrecognized contaminant exposures.

North Dakota, South Dakota

Assessment of streamflow trends in the eastern Dakotas, water years 1960–2019

Hydrologic extremes, whether periods of drought or flooding, are occurring more frequently with greater severity and can have substantial economic impacts. Along with flooding, the timing and volume of streamflow also is changing across the United States. The focus of this report is to characterize a unique trend in mean annual streamflow occurring in eastern North and South Dakota, hereafter referred to as the eastern Dakotas, that is not being observed anywhere else in the conterminous United States. Streamflow records for 1,853 U.S. Geological Survey streamgages obtained from the U.S. Geological Survey National Water Information System database with a continuous record of mean annual streamflow during water years 1960–2019 were included in this study. Using a Kendall tau statistical test ( p -value less than or equal to 0.10), 573 streamgages had a statistically significant upward trend in mean annual streamflow and are primarily located in the Midwest and northeastern United States. Of the streamgages, 182 had a statistically significant downward trend and are located primarily in the western and southeastern States. Several sites had increases in streamflow between 100 and 500 percent. Most of the streamgages with the highest increases in mean annual streamflow are along the same rivers in the eastern Dakotas, regardless of basin size. A comparison of mean annual streamflow of the last decade (2010–19) to the first decade (1960–69) of the study period shows that the largest increases in annual streamflow volumes in the United States also are in the eastern Dakotas. Among all 1,853 streamgages in the United States, the Sheyenne River near Warwick, North Dakota (U.S. Geological Survey station 05056000), has the greatest percent change, with an increase of 486 percent. Several factors may be contributing to increasing trends in streamflow in the eastern Dakotas and may include, in part, precipitation changes owing to climatic variation within the region, geologic makeup of the subsurface, and land-use changes. A better understanding of these research areas will help producers, resource managers, and infrastructure engineers to make more informed environmental and economic decisions.

North Dakota, South Dakota

Virtual training prepared for the former Afghanistan Ministry of Energy and Water—Streamgaging, fluvial sediment sampling, bathymetry, and streamflow and sediment modeling

The U.S. Geological Survey (USGS) created a virtual training series for the Afghanistan Ministry of Energy and Water (MEW), now known as the National Water Affairs Regulation Authority (NWARA), to provide critical hydrological training as an alternative to an in-person training. The USGS was scheduled to provide in-person surface-water training for NWARA during 2020; however, travel was halted because of the Coronavirus disease 2019 (COVID–19) pandemic. The virtual training consisted of prerecorded and live presentations that were scheduled during 4 weeks in August 2021. However, the training was halted after the second week due to the collapse of the Afghan Government. Fortunately, the prerecorded presentations and training materials were delivered before the trainings were halted, so they can be viewed or shared by the participants in the future. A benefit to having produced prerecorded trainings is that USGS can leverage or adapt the trainings for nongovernmental organizations (NGOs) involved in humanitarian water relief efforts in Afghanistan or can be used for other international training efforts.

Fact Sheet

Estimates of water use associated with continuous oil and gas development in the Permian Basin, Texas and New Mexico, 2010–19, with comparisons to the Williston Basin, North Dakota and Montana

The Permian Basin, in west Texas and southeastern New Mexico is one of the largest conventional oil and gas reservoirs in the United States and is becoming one of the world’s largest continuous oil and gas (COG) reservoirs. Advances in technology have enabled oil and gas to be extracted from reservoirs that historically were developed using conventional, or vertical, well drilling techniques. Conventional oil and gas reservoirs have discrete deposits that are well defined and are typically trapped by an overlying geologic formation or caprock, whereas COG reservoirs contain deposits that are distributed evenly throughout the geologic formation, typically have much lower permeability (the capacity of a porous rock to transmit a fluid) than the con­ventional deposits, and require specialized horizontal extraction techniques. The methods to extract the oil from the two different reservoirs require differing amounts of water, and the horizontal extraction methods typically require substantially more water. In 2015, the U.S. Geological Survey started a topical study to quantify water used during COG development. The Permian Basin, which contains both types of reservoirs (continuous and conventional), was the second basin in the United States in the U.S. Geological Survey’s topical study to quantify water used during COG development.

Montana, New Mexico, North Dakota, Texas

Conceptual and numerical groundwater flow model of the Cedar River alluvial aquifer system with simulation of drought stress on groundwater availability near Cedar Rapids, Iowa, for 2011 through 2013

Between July 2011 and February 2013, the City of Cedar Rapids observed water level declines in their horizontal collector wells approaching 11 meters. As a result, pumping from these production wells had to be halted, and questions were raised about the reliability of the alluvial aquifer under future drought conditions. The U.S. Geological Survey, in cooperation with the City of Cedar Rapids, completed a study to better understand the effects of drought stress on the Cedar River alluvial aquifer using a numerical groundwater flow model. Previously published groundwater flow models were combined with newly collected airborne, waterborne, down-hole, and land-based geophysical survey data and provided a detailed three-dimensional lithologic model of the Cedar River alluvial aquifer and surrounding area. An improved conceptual model for the groundwater flow system and a lithologic model were used to build and inform a numerical groundwater flow model capable of simulating water levels observed in the City of Cedar Rapids horizontal collector wells during the 2012 drought. Model performance was assessed primarily on the ability of the model to simulate water table elevation at six monitoring wells. Statistical tests were used to assess the numerical model during the calibration period, and results varied from satisfactory to unsatisfactory, likely because of stage changes in the Cedar River and production well withdrawal rates near monitoring wells. Simulated water levels during the 2012 drought indicated a depression near the horizontal collector wells, although simulated elevations at these locations and at monitoring wells were generally overestimated compared to measured values. The numerical groundwater flow model was modified to account for a decrease in seepage rate caused by low flow in the Cedar River and increased production. With seepage rate modification, model results improved; the simulated water table elevations were like those observed in horizontal collector and monitoring wells. Results demonstrated the ability of the model to simulate water levels observed in the horizontal collector wells during the 2012 drought when accounting for a decrease in infiltration from the Cedar River.

Iowa

Estimates of water use associated with continuous oil and gas development in the Permian Basin, Texas and New Mexico, 2010–19

In 2015, the U.S. Geological Survey started a topical study to quantify water use in areas of continuous oil and gas (COG) development. The first phase of the study was completed in 2019 and analyzed the Williston Basin. The second phase of the study analyzed the Permian Basin using the same techniques and approaches used for the Williston Basin analysis. The Permian Basin was selected for the second phase of water-use analysis for the following reasons: (1) the basin has the largest undiscovered technically recoverable oil and gas resource in the United States, (2) the basin has a continuous resource in tight shale that primarily produces oil, and (3) the basin is within the contiguous United States. This study used data from 60 counties in Texas and New Mexico with spatial coverage based on the Permian Basin extent defined by the U.S. Energy Information Administration, a representation of the geologically defined Permian Basin. Data from several sources were used in the analysis of direct, indirect, and ancillary water use associated with COG development in the Permian Basin and are available in an associated data release. Hydraulic fracturing water-use data were used to determine the start of the recent (before 2019) COG development boom in oil production in the Permian Basin in the same way that the data were used for the Williston Basin study. Water-use data were aggregated by county and year, which were the sampling units used in the analysis. The water-use analysis of the Permian Basin contained three elements: (1) estimates of water use, in million gallons, by county and year; (2) coefficients of water use from regression models, in million gallons per developed oil and gas well; and (3) performance (based on goodness-of-fit metrics) of the regression models in estimating the observed water use. Coefficients from the linear and quantile regression models of direct, indirect, and ancillary water use in the Permian Basin were produced as aggregate values for the counties and years. The mean estimate of direct water use had a 95-percent confidence interval of 4.13–5.45 million gallons (Mgal) per developed oil and gas well. The coefficient from the linear regression model of indirect water use was 0.111 Mgal per well, with a 95-percent confidence interval of 0.104–0.117 Mgal per well. The mean estimate of ancillary water use in the Permian Basin was 1.09 Mgal per well, with a 95-percent confidence interval of 1.05–1.13 Mgal per well. Model performance was evaluated with goodness-of-fit metrics including coefficient of determination ( R 2 ), root mean square error, and the ratio of root mean square error to standard deviation of observations computed from leave-one-out cross validation of the linear and quantile regression models of direct, indirect, and ancillary water use. Model performance for direct water use was acceptable, with an R 2 value of 0.91. The model performance of indirect water use was acceptable, with an R 2 value of 0.89. Values of R 2 for the ancillary water-use categories were at least 0.89. Annual mean estimates for hydraulic fracturing, cementing, drilling, indirect, and ancillary water use per well for the years 2010–17 were comparable between the Permian and Williston Basins. Hydraulic fracturing water use increased similarly from 2010 to 2015 in the Permian Basin and the Williston Basin, increasing from 0.6 Mgal per well in 2010 to 5.4 Mgal per well in 2015 in the Permian Basin and from 1.4 Mgal per well in 2010 to 4.7 Mgal per well in 2015 in the Williston Basin. By design, the Permian water-use assessment is a simplification of a complex and continually developing system and therefore has uncertainty and limitations in the interpretation of results. Despite the modeling limitations, the results summarized in the report, when compared to other studies, compare well with water-use estimations. The favorable comparison highlights the transferability of the water-use methodology to other areas of COG development.

New Mexico, Texas

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

Estimates of water use associated with continuous oil and gas development in the Williston Basin, North Dakota and Montana, 2007–17

This study of water use associated with development of continuous oil and gas resources in the Williston Basin is intended to provide a preliminary model-based analysis of water use in major regions of production of continuous oil and gas resources in the United States. Direct, indirect, and ancillary water use associated with development of continuous oil and gas resources in the Williston Basin was estimated in North Dakota and Montana from 2007 to 2017. Water-use data were aggregated by county and year, which were the sampling units used in this analysis. Linear and quantile regression models of water use in relation to the number of oil and gas wells developed were fit for the direct, indirect, and ancillary water-use categories for each State. A 95-percent confidence interval for each parameter estimate from the linear regression models was computed as a measure of uncertainty. Additional information on uncertainty can be gained from modeling other distribution parameters, so quantile regression models of the 5th, 50th, and 95th percentiles also were fit. To assess uncertainty in the estimates from the regression models of direct, indirect, and ancillary water use, leave-one-out cross-validation was used. Model performance was evaluated with three goodness-of-fit metrics used to compare the estimates and observations of water use. Mean annual direct and indirect water use for development of continuous oil and gas resources in North Dakota was estimated at 4,512 million gallons (Mgal) per year (Mgal/yr), with a 95-percent confidence interval of 4,021–5,152 Mgal/yr, and in Montana was estimated at 196 Mgal/yr, with a 95-percent confidence interval of 189–203 Mgal/yr. Ancillary water use (for domestic and public supply) had an estimated annual mean of 2,753 Mgal/yr in North Dakota and 396 Mgal/yr in Montana. The coefficient from the linear regression model of direct water use was 3.86 Mgal per well and hydraulic fracturing water use was 3.70 Mgal per well for North Dakota. The mean estimate of direct water use had a 95-percent confidence interval of 3.48–4.23 Mgal per well. For North Dakota, the coefficient from the linear regression model of indirect water use was 0.453 Mgal per well, with a 95-percent confidence interval of 0.415–0.492 Mgal per well. Direct and indirect water use had a mean estimate of about 4.31 Mgal per well in North Dakota. The mean estimate of ancillary water use (for domestic and public supply) in North Dakota was 2.03 Mgal per well, with a 95-percent confidence interval of 1.76–2.31 Mgal per well. For Montana, the linear regression model of hydraulic fracturing water use had a mean estimate of 2.04 Mgal per well. The 95-percent confidence interval for the mean estimate was 1.80–2.28 Mgal per well. Direct and indirect water use in Montana had a mean estimate of 2.49 Mgal per well. The mean estimate of ancillary water use (for domestic and public supply) in Montana was 2.43 Mgal per well, with a 95-percent confidence interval of 1.76–3.11 Mgal per well.

Montana, North Dakota, South Dakota

Groundwater characterization of the Madison aquifer near Jewel Cave National Monument, South Dakota

Jewel Cave National Monument in the Black Hills of southwestern South Dakota has more than 200 miles of mapped cave passages and several subterranean lakes that have been discovered since 2015. Jewel Cave is one of the world’s longest known caves and its natural beauty and unique natural cave features led U.S. President Theodore Roosevelt to designate the cave as a national monument in 1908. Jewel Cave was naturally formed in the regionally extensive Madison Limestone, which is characterized as a carbonate karst environment (containing caves and sinkholes) with extensive subterranean cave networks and losing streams at the land surface. Preserving and protecting the cave is an important element of the National Park Service mission, and understanding the hydrogeologic connection between the surface and the subsurface is essential for ensuring the preservation and protection of the cave for future generations. A component in preserving and protecting the park includes the improved understanding of groundwater flow and vulnerability of the subsurface, which allows scientists, park managers, the visiting public, and the surrounding communities to better manage, protect, and preserve the site and its unique natural features.

South Dakota

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

Analytical framework to estimate water use associated with continuous oil and gas development

An analytical framework was designed to estimate water use associated with continuous oil and gas (COG) development in support of the U.S. Geological Survey Water Availability and Use Science Program. This framework was developed to better understand the relation between the production of COG resources for energy and the amount of water needed to sustain this type of energy development in the United States. The total mean undiscovered, technically recoverable volume of COG has increased, highlighting the continued need to develop approaches to better characterize water use associated with COG development. The analytical framework can be used to estimate water use associated with COG development for three water-use components—direct, indirect, and ancillary water use—that are related to the life cycle of COG development. Direct water use is defined as water used in a wellbore to complete a well, including the water used for drilling, cementing, stimulating, and maintaining the well during production. Indirect water use is the water used at or near the well site, including water used for dust abatement, for cleaning equipment, and for crew and staff use. Ancillary water use is all other water used during the life cycle of COG development that is not categorized as direct or indirect, such as additional local or regional water use resulting from a change (for example, population) related to COG development. The analytical framework includes the data inputs, the processes involved in estimating the water-use coefficients and analyzing their uncertainties, and the outputs. The analytical framework was developed as an R script, which contains the statistical models used to estimate water-use components. The availability of data across COG reservoirs in the United States is variable and presents challenges for estimating water use for extracting COG from their reservoirs; thus, the R script can be modified for the types of data available within a COG reservoir, the extent and resolution of data available for each water-use component, and the desired output of the water-use assessment. The script was written so that the units of the data in the script were standardized. Water-use estimates were simulated for the mean and 10th, 50th, and 90th percentiles of the data distributions. Uncertainties were quantified with confidence intervals for the estimated coefficients. Uncertainty for estimated or simulated data can be calculated with the R script by providing a range of representative values that are within the appropriate confidence intervals of the mean of the data.

Lower 48 states

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

Delineation of selected lithologic units using airborne electromagnetic data near Cedar Rapids, Iowa

The U.S. Geological Survey, in cooperation with the City of Cedar Rapids, began a study in 2013 to better understand the effects of drought stress on the Cedar River alluvial aquifer. After an evaluation of the existing groundwater-flow models for the alluvial aquifer, a plan was begun to construct an updated groundwater-flow model capable of evaluating the effect of prolonged drought and increased demand. As part of the effort to update the existing groundwater-flow model, data were collected during an airborne electromagnetic (AEM) survey in May 2017. The study area for the AEM survey encompasses about 53 square kilometers of the Cedar River Basin, west of Cedar Rapids, Iowa, and includes a 19-kilometer reach of the Cedar River. The AEM survey of the Cedar River alluvial aquifer and adjacent areas was completed to characterize the subsurface geology of the area to refine a lithologic framework. The collected AEM data were postprocessed by numerical inversion using the program EM1DFM to produce subsurface apparent resistivity cross sections. Changes observed in resistivity profile values with depth were used to infer lithologic changes and delineate three of the four lithologic units designated in the lithologic framework for this area: alluvial deposits, glacial till, and bedrock; hereafter referred to as the “lithologic framework.” The fourth unit, composed of surficial eolian sediments, was not delineated in these profiles because these units are thin and discontinuous and are not reliably distinguishable from flood plain alluvial deposits. For the purposes of delineating lithologic units using the AEM data, bedrock was assumed to be the lowest unit in a profile, glacial till was deposited on a bedrock surface, and alluvium was deposited on erosional till or bedrock surfaces. A three-dimensional fence diagram was created as part of the lithologic framework to further define the extent and thickness of the lithologic units near the Cedar River alluvial aquifer. The fence diagram shows a three-dimensional perspective of unit thickness, extent, and orientation of the delineated lithologic framework. A lithologic framework, by design, is intended to represent a simplification of a more complex natural system through data interpolation between known points, which usually are lithologic logs. The resistivity profiles produced from the AEM survey allow for continuous mapping and accurate interpolation of lithology between lithologic logs; however, the apparent resistivity value may reflect several characteristics of subsurface materials including variations in lithology, porosity, water quality, grain sorting, and degree of saturation. In this study, the only variables considered were those related to changes in the subsurface material.

Iowa

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

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

As the Nation’s principal earth-science information agency, the U.S. Geological Survey (USGS) is depended upon to collect accurate data and produce factual and impartial interpretive reports. Methods for data collection and analysis that were developed by the USGS have become standard techniques used by numerous Federal, State, and local agencies and by private enterprises. The USGS has implemented a program designed to ensure that all scientific work done by or for USGS Water Science Centers is done in accordance with a quality-assurance plan. The implementation of a groundwater quality-assurance plan will enhance groundwater data collected by the USGS. This report is a quality-assurance plan for groundwater activities conducted by the USGS Dakota Water Science Center and is meant to complement qualityassurance plans for surface-water and water-quality activities and similar plans for the Dakota Water Science Center and general project activities throughout the USGS.

Open-File Report

Conceptual model to assess water use associated with the life cycle of unconventional oil and gas development

As the demand for energy increases in the United States, so does the demand for water used to produce many forms of that energy. Technological advances, limited access to conventional oil and gas accumulations, and the rise of oil and gas prices resulted in increased development of unconventional oil and gas (UOG) accumulations. Unconventional oil and gas is developed using a method that combines directional drilling and hydraulic fracturing techniques, allowing for greater oil and gas production from previously unrecoverable reservoirs. Quantification of the water resources required for UOG development and production is difficult because of disparate data sources, variable reporting requirements across boundaries (local, State, and national), and incomplete or proprietary datasets. A topical study was started in 2015 under the U.S. Geological Survey’s Water Availability and Use Science Program, as part of the directive in the Secure Water Act for the U.S. Geological Survey to conduct a National Water Census, to better understand the relation between production of UOG resources for energy and the amount of water needed to produce and sustain this type of energy development in the United States. The Water Availability and Use Science Program goal for this topical study is to develop and apply a statistical model to better estimate the water use associated with UOG development, regardless of the location and target geologic formation. As a first step, a conceptual model has been developed to characterize the life cycle of water use in areas of UOG development. Categories of water use and the way water-use data are collected might change over time; therefore, a generic approach was used in developing the conceptual model to allow for greater flexibility in adapting to future changes or newly available data. UOG development can be summarized into four stages: predrilling construction, drilling, hydraulic fracturing, and ongoing production. The water used in UOG production can be categorized further as direct, indirect, or ancillary water use. Direct water use is defined as the water used for drilling and hydraulic fracturing a well and for maintaining the well during ongoing production. Indirect water use is defined as the water used at or near a well pad. The water used for dust abatement also is considered an indirect use but may be applied away from the well pad. Ancillary water use is defined as the additional local or regional water use resulting from a change (for example, population) directly related to UOG development throughout the life cycle that is not used directly in the well or indirectly for any other purpose at the well pad. The conceptual model presented in this report consists of five elements: (1) input data, (2) processes, (3) decisions, (4) output data, and (5) outcomes. The input data requirements for estimating water use associated with UOG development are somewhat onerous, and obtaining suitable datasets can be challenging because local, State, and Federal agencies do not collect data similarly. The quality of a water-use assessment that uses the conceptual model presented in this report is dependent on the quality and quantity of data that are available for a UOG play. The conceptual model can be used for an assessment with sparse data; however, having sparse data likely will result in greater uncertainty in the water-use estimates. The conceptual model presented in this report is designed to be robust to characterize and simulate the data processing to estimate water use associated with UOG development. Although the results of an analysis that includes missing data have greater uncertainty, the analysis still can be insightful because it can establish a baseline estimate of UOG water use that may be refined further as more data become available. Analysis of models that include missing data also could aid in identifying the data most needed for future water-use estimates. Characterizing individual model limitations is important because the conceptual model can be used in future water-use studies to facilitate data compiling, data processing, estimating, and assessing UOG activities regardless of location.

Scientific Investigations Report

Hydrogeologic framework and groundwater conditions of the Ararat Basin in Armenia

Armenia is a landlocked country located in the mountainous Caucasus region between Asia and Europe. It shares borders with the countries of Georgia on the north, Azerbaijan on the east, Iran on the south, and Turkey and Azerbaijan on the west. The Ararat Basin is a transboundary basin in Armenia and Turkey. The Ararat Basin (or Ararat Valley) is an intermountain depression that contains the Aras River and its tributaries, which also form the border between Armenia and Turkey and divide the basin into northern and southern regions. The Ararat Basin also contains Armenia’s largest agricultural and fish farming zone that is supplied by high-quality water from wells completed in the artesian aquifers that underlie the basin. Groundwater constitutes about 40 percent of all water use, and groundwater provides 96 percent of the water used for drinking purposes in Armenia. Since 2000, groundwater withdrawals and consumption in the Ararat Basin of Armenia have increased because of the growth of aquaculture and other uses. Increased groundwater withdrawals caused decreased springflow, reduced well discharges, falling water levels, and a reduction of the number of flowing artesian wells in the southern part of Ararat Basin in Armenia. In 2016, the U.S. Geological Survey and the U.S. Agency for International Development (USAID) began a cooperative study in Armenia to share science and field techniques to increase the country’s capabilities for groundwater study and modeling. The purpose of this report is to describe the hydrogeologic framework and groundwater conditions of the Ararat Basin in Armenia based on data collected in 2016 and previous hydrogeologic studies. The study area includes the Ararat Basin in Armenia. This report was completed through a partnership with USAID/Armenia in the implementation of its Science, Technology, Innovation, and Partnerships effort through the Advanced Science and Partnerships for Integrated Resource Development program and associated partners, including the Government of Armenia, Armenia’s Hydrogeological Monitoring Center, and the USAID Global Development Lab and its GeoCenter. The hydrogeologic framework of the Ararat Basin includes several basin-fill stratigraphic units consisting of interbedded dense clays, gravels, sands, volcanic basalts, and andesite deposits. Previously published cross sections and well lithologic logs were used to map nine general hydrogeologic units. Hydrogeologic units were mapped based on lithology and water-bearing potential. Water-level data measured in the water-bearing hydrogeologic units 2, 4, 6, and 8 in 2016 were used to create potentiometric surface maps. In hydrogeologic unit 2, the estimated direction of groundwater flow is from the west to north in the western part of the basin (away from the Aras River) and from north to south (toward the Aras River) in the eastern part of the basin. In hydrogeologic unit 4, the direction of groundwater flow is generally from west to east and north to south (toward the Aras River) except in the western part of the basin where groundwater flow is toward the north or northwest. Hydrogeologic unit 6 has the same general pattern of groundwater flow as unit 4. Hydrogeologic unit 8 is the deepest of the water-bearing units and is confined in the basin. Groundwater flow generally is from the south to north (away from the Aras River) in the western part of the basin and from west to east and north to south (toward the Aras River) elsewhere in the basin. In addition to water levels, personnel from Armenia’s Hydrogeological Monitoring Center also measured specific conductance at 540 wells and temperature at 2,470 wells in the Ararat Basin using U.S. Geological Survey protocols in 2016. The minimum specific conductance was 377 microsiemens per centimeter (μS/cm), the maximum value was 4,000 μS/cm, and the mean was 998 μS/cm. The maximum water temperature was 24.2 degrees Celsius. An analysis between water temperature and well depth indicated no relation; however, spatially, most wells with cooler water temperatures were within the 2016 pressure boundary or in the western part of the basin. Wells with generally warmer water temperatures were in the eastern part of the basin. Samples were collected from four groundwater sites and one surface-water site by the U.S. Geological Survey in 2016. The stable-isotope values were similar for all five sites, indicating similar recharge sources for the sampled wells. The Hrazdan River sample was consistent with the groundwater samples, indicating the river could serve as a source of recharge to the Ararat artesian aquifer.

Ararat Basin

Construction of a groundwater-flow model for the Big Sioux Aquifer using airborne electromagnetic methods, Sioux Falls, South Dakota

The city of Sioux Falls is the fastest growing community in South Dakota. In response to this continued growth and planning for future development, Sioux Falls requires a sustainable supply of municipal water. Planning and managing sustainable groundwater supplies requires a thorough understanding of local groundwater resources. The Big Sioux aquifer consists of glacial outwash sands and gravels and is hydraulically connected to the Big Sioux River, which provided about 90 percent of the city’s source-water production in 2015. Managing sustainable groundwater supplies also requires an understanding of groundwater availability. An effective mechanism to inform water management decisions is the development and utilization of a groundwater-flow model. A groundwater-flow model provides a quantitative framework for synthesizing field information and conceptualizing hydrogeologic processes. These groundwater-flow models can support decision making processes by mapping and characterizing the aquifer. Accordingly, the city of Sioux Falls partnered with the U.S. Geological Survey to construct a groundwater-flow model. Model inputs will include data from advanced geophysical techniques, specifically airborne electromagnetic methods.

South Dakota