Geology topics

Larry D. Putnam

Publications and source records attributed to Larry D. Putnam.

13 recordsLinked to original sources

USGS2015

Conceptual and numerical models of groundwater flow in the Ogallala and Arikaree aquifers, Pine Ridge Indian Reservation area, South Dakota, water years 1980-2009

The Ogallala and Arikaree aquifers are the largest sources of groundwater on the Pine Ridge Indian Reservation and are used extensively for irrigation and public and domestic water supplies. To assess the potential for decreased water levels and discharge to streams in the Pine Ridge Indian Reservation, conceptual and numerical models of groundwater flow in the Ogallala and Arikaree aquifers in southwestern South Dakota were developed by the U.S. Geological Survey in cooperation with the Oglala Sioux Tribe. The study area includes most of the Pine Ridge Reservation in Jackson and Shannon Counties and Indian trust lands in Bennett County in southwestern South Dakota. The High Plains aquifer, which includes the Ogallala and Arikaree aquifers, generally is less developed in South Dakota compared with other areas underlain by this aquifer; therefore, water levels in the High Plains aquifer in South Dakota generally fluctuated by less than 5 feet (ft) from 1980 to 1999. Despite minimal water-level changes in the High Plains aquifer in South Dakota, extensive withdrawals of groundwater for irrigation have caused water-level declines in many areas and increased concerns about the long-term sustainability of the aquifer; therefore, continued or increased withdrawals from the aquifer or prolonged drought may have the potential to affect water levels within the aquifer and discharge to important streams in the area. The Ogallala and Arikaree aquifers generally consist of poorly consolidated claystones, siltstones, sandstones, and shale deposited in fluvial and lacustrine environments. Saturated thicknesses ranged from 10 to 314 ft for the Ogllala aquifer and from 10 to 862 ft for the Arikaree aquifer. Previous hydraulic conductivity estimates ranged from less than 1 to 180 feet per day (ft/d) for the Ogallala aquifer and from less than 1 to 13 ft/d for the Arikaree aquifer. Recharge to the Ogallala and Arikaree aquifers is from precipitation on the outcrop areas, and discharge occurs through evapotranspiration, discharge to streams, and well withdrawals. Evapotranspiration generally occurs in topographically low areas along streams, and maximum evapotranspiration occurs when the water level is at the land surface. The generalized groundwater-flow direction is to the northeast with local flow towards streams. Precipitation for water years 1980–2009 ranged from about 11 to 39 inches per year (in/yr) and averaged about 19 in/yr. Estimated mean recharge for water years 1980–2009 was about 17.3 percent of precipitation for the Ogallala aquifer and 7.9 percent of precipitation for the Arikaree aquifer. The estimated mean maximum evapotranspiration for water years 1980–2009 was about 35 in/yr. Estimated mean base flow for gaged streams was about 0.06 cubic foot per second (ft 3 /s) per square mile of drainage area. Estimated mean total water use for water years 1980–2009 was 5.4 ft 3 /s from the Ogallala aquifer and 7.1 ft 3 /s from the Arikaree aquifer. A two-layer numerical groundwater-flow model was constructed using MODFLOW–NWT with a uniformly spaced grid consisting of 166 rows and 288 columns with cells 1,640 ft on a side. The numerical model of the Ogallala and Arikaree aquifers was used to simulate steady-state and transient conditions for water years 1980–2009. Model calibration was accomplished using the Parameter ESTimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. Aquifer boundaries were no-flow on the northern and western sides and constant-head on the southern and eastern sides. The mean arithmetic difference was 1.4 ft between the 731 simulated and observed hydraulic heads in the Ogallala aquifer and 9.8 ft between the 2,754 simulated and observed hydraulic heads in the Arikaree aquifer. Simulated mean discharge from the Ogallala and Arikaree aquifers to selected stream reaches was 92.1 ft 3 /s compared to estimated discharge of 88.7 ft 3 /s. Calibrated recharge for the transient simulation averaged 3.3 in/yr for the Ogallala aquifer and 1.1 in/yr for the Arikaree aquifer. The mean maximum potential evapotranspiration rate was 35.4 in/yr. Streambed conductance for perennial stream reaches averaged 530 feet squared per day. Horizontal hydraulic conductivity averaged 27 ft/d for the Ogallala aquifer and 1.0 ft/d for the Arikaree aquifer. The vertical hydraulic conductivity averaged 1.4 ft/d for the Ogallala aquifer and 0.004 ft/d for the Arikaree aquifer. Specific yield for the Ogallala aquifer was 0.15 (dimensionless) and averaged 0.02 for the Arikaree aquifer. Specific storage for the Arikaree aquifer was 1.7x10 -6 per foot. Simulated steady-state model inflow and outflow was 459 ft 3 /s. The percentages of inflows were 17 percent from constant-head boundaries, 9 percent from streams, and 74 percent from recharge. Percentages of outflow were 8 percent to constant-head boundaries, 1 percent to wells, 31 percent to streams, and 59 percent to evapotranspiration. Simulated net inflow from the Ogallala aquifer to the Arikaree aquifer ranged from about 22 ft 3 /s in dry years to about 37 ft 3 /s in wet years. Two hypothetical future stress scenarios were simulated using input from the 30-year calibrated simulation of water years 1980–2009. The first hypothetical scenario represented an increase in groundwater withdrawals from 50 hypothetical production wells completed in the Arikaree aquifer. At the end of the 30-year hypothetical increased pumping simulation, water levels declined as much as 66 ft in the Arikaree aquifer, decreased discharge to streams accounted for about 26 percent (2.6 ft 3 /s) of increased withdrawals, and decreased evapotranspiration accounted for about 53 (5.3 ft 3 /s) percent of increased withdrawals. The second hypothetical scenario represented a 30-year period of decreased recharge (drought) by decreasing recharge 0.2 inch (24 ft 3 /s) for each water year. At the end of the hypothetical drought simulation, water levels declined as much as 10.9 ft in the Arikaree aquifer, decreased discharge to streams accounted for about 23 percent (5.5 ft 3 /s) of decreased recharge, and decreased evapotranspiration accounted for about 72 percent (17.3 ft 3 /s) of decreased recharge. The numerical model is a tool that could be used to better understand the flow system of the Ogallala and Arikaree aquifers, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the Pine Ridge Reservation area in Bennett, Jackson, and Shannon Counties. The model also is useful to help assess the response of the aquifer to additional stress, including potential increased well withdrawals and potential drought conditions.

South Dakota↗

USGS2013

Conceptual and numerical models of groundwater flow in the Ogallala aquifer in Gregory and Tripp Counties, South Dakota, water years 1985--2009

The Ogallala aquifer is an important water resource for the Rosebud Sioux Tribe in Gregory and Tripp Counties in south-central South Dakota and is used for irrigation, public supply, domestic, and stock water supplies. To better understand groundwater flow in the Ogallala aquifer, conceptual and numerical models of groundwater flow were developed for the aquifer. A conceptual model of the Ogallala aquifer was used to analyze groundwater flow and develop a numerical model to simulate groundwater flow in the aquifer. The MODFLOW–NWT model was used to simulate transient groundwater conditions for water years 1985–2009. The model was calibrated using statistical parameter estimation techniques. Potential future scenarios were simulated using the input parameters from the calibrated model for simulations of potential future drought and future increased pumping. Transient simulations were completed with the numerical model. A 200-year transient initialization period was used to establish starting conditions for the subsequent 25-year simulation of water years 1985–2009. The 25-year simulation was discretized into three seasonal stress periods per year and used to simulate transient conditions. A single-layer model was used to simulate flow and mass balance in the Ogallala aquifer with a grid of 133 rows and 282 columns and a uniform spacing of 500 meters (1,640 feet). Regional inflow and outflow were simulated along the western and southern boundaries using specified-head cells. All other boundaries were simulated using no-flow cells. Recharge to the aquifer occurs through precipitation on the outcrop area. Model calibration was accomplished using the Parameter Estimation (PEST) program that adjusted individual model input parameters and assessed the difference between estimated and model-simulated values of hydraulic head and base flow. This program was designed to estimate parameter values that are statistically the most likely set of values to result in the smallest differences between simulated and observed values, within a given set of constraints. The potentiometric surface of the aquifer calculated during the 200-year initialization period established initial conditions for the transient simulation. Water levels for 38 observation wells were used to calibrate the 25-year simulation. Simulated hydraulic heads for the transient simulation were within plus or minus 20 feet of observed values for 95 percent of observation wells, and the mean absolute difference was 5.1 feet. Calibrated hydraulic conductivity ranged from 0.9 to 227 feet per day (ft/d). The annual recharge rates for the transient simulation (water years 1985–2009) ranged from 0.60 to 6.96 inches, with a mean of 3.68 inches for the Ogallala aquifer. This represents a mean recharge rate of 280.5 ft 3 /s for the model area. Discharge from the aquifer occurs through evapotranspiration, discharge to streams through river leakage and flow from springs and seeps, and well withdrawals. Water is withdrawn from wells for irrigation, public supply, domestic, and stock uses. Simulated mean discharge rates for water years 1985–2009 were about 185 cubic feet per second (ft 3 /s) for evapotranspiration, 66.7 ft 3 /s for discharge to streams, and 5.48 ft 3 /s for well withdrawals. Simulated annual evapotranspiration rates ranged from about 128 to 254 ft 3 /s, and outflow to streams ranged from 52.2 to 79.9 ft 3 /s. A sensitivity analysis was used to examine the response of the calibrated model to changes in model parameters for horizontal hydraulic conductivity, recharge, evapotranspiration, and spring and riverbed conductance. The model was most sensitive to recharge and maximum potential evapotranspiration and least sensitive to riverbed and spring conductances. Two potential future scenarios were simulated: a potential drought scenario and a potential increased pumping scenario. To simulate a potential drought scenario, a synthetic drought record was created, the mean of which was equal to 60 percent of the mean estimated recharge rate for the 25-year simulation period. Compared with the results of the calibrated model (non-drought simulation), the simulation representing a potential drought scenario resulted in water-level decreases of as much as 30 feet for the Ogallala aquifer. To simulate the effects of potential future increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 25-year simulation period. Compared with the results of the calibrated model, the simulation representing an increased pumping scenario resulted in water-level decreases of as much as 26 feet for the Ogallala aquifer. Groundwater budgets for the potential future scenario simulations were compared with the transient simulation representing water years 1985–2009. The simulation representing a potential drought scenario resulted in lower aquifer recharge from precipitation and decreased discharge from streams, springs, seeps, and evapotranspiration. The simulation representing a potential increased pumping scenario was similar to results from the transient simulation, with a slight increase in well withdrawals and a slight decrease in discharge from river leakage and evapotranspiration. This numerical model is suitable as a tool that could be used to better understand the flow system of the Ogallala aquifer, to approximate hydraulic heads in the aquifer, and to estimate discharge to rivers, springs, and seeps in the study area. The model also is useful to help assess the response of the aquifer to additional stresses, including potential drought conditions and increased well withdrawals.

South Dakota↗

USGS2012

Conceptual and numerical models of the glacial aquifer system north of Aberdeen, South Dakota

This U.S. Geological Survey report documents a conceptual and numerical model of the glacial aquifer system north of Aberdeen, South Dakota, that can be used to evaluate and manage the city of Aberdeen's water resources. The glacial aquifer system in the model area includes the Elm, Middle James, and Deep James aquifers, with intervening confining units composed of glacial till. The Elm aquifer ranged in thickness from less than 1 to about 95 feet (ft), with an average thickness of about 24 ft; the Middle James aquifer ranged in thickness from less than 1 to 91 ft, with an average thickness of 13 ft; and the Deep James aquifer ranged in thickness from less than 1 to 165 ft, with an average thickness of 23 ft. The confining units between the aquifers consisted of glacial till and ranged in thickness from 0 to 280 ft. The general direction of groundwater flow in the Elm aquifer in the model area was from northwest to southeast following the topography. Groundwater flow in the Middle James aquifer was to the southeast. Sparse data indicated a fairly flat potentiometric surface for the Deep James aquifer. Horizontal hydraulic conductivity for the Elm aquifer determined from aquifer tests ranged from 97 to 418 feet per day (ft/d), and a confined storage coefficient was determined to be 2.4x10 -5 . Estimates of the vertical hydraulic conductivity of the sediments separating the Elm River from the Elm aquifer, determined from the analysis of temperature gradients, ranged from 0.14 to 2.48 ft/d. Average annual precipitation in the model area was 19.6 inches per year (in/yr), and agriculture was the primary land use. Recharge to the Elm aquifer was by infiltration of precipitation through overlying outwash, lake sediments, and glacial till. The annual recharge for the model area, calculated by using a soil-water-balance method for water year (WY) 1975-2009, ranged from 0.028 inch in WY 1980 to 4.52 inches in WY 1986, with a mean of 1.56 inches. The annual potential evapotranspiration, calculated in soil-water-balance analysis, ranged from 21.8 inches in WY 1983 to 27.0 inches in WY 1985, with a mean of 24.6 inches. Water use from the glacial aquifer system primarily was from the Elm aquifer for irrigation, municipal, and suburban water supplies, and the annual rate ranged from 1.0 to 2.4 cubic feet per second (ft 3 /s). The MODFLOW-2005 numerical model represented the Elm aquifer, the Middle James aquifer, and the Deep James aquifer with model layers 1-3 respectively separated by confining layers 1-2 respectively. Groundwater flow was simulated with 75 stress periods beginning October 1, 1974, and ending September 30, 2009. Model grid spacing was 200 by 200 ft and boundaries were represented by specified-head boundaries and no-flow boundaries. The model used parameter estimation that focused on minimizing the difference between 954 observed and simulated hydraulic heads for 135 wells. Calibrated mean horizontal hydraulic conductivity values for model layers 1-3 were 94, 41, and 30 ft/d respectively. Vertical hydraulic conductivity values for confining layers 1 and 2 were 0.0002 and 0.0003 ft/d, respectively. Calibrated specific yield for model layer 1was 0.1 and specific storage ranged from 0.0003 to 0.0005 per foot. Calibrated mean recharge rates ranged from 2.5 in/yr where glacial till thickness was less than 10 ft to 0.8 in/yr where glacial till thickness was greater than 30 ft. Calibrated mean annual evapotranspiration rate was 8.8 in/yr. Simulated net streamflow gain from model layer 1 was 3.1 ft 3 /s.

South Dakota↗

USGS2010

Simulated groundwater flow in the Ogallala and Arikaree aquifers, Rosebud Indian Reservation area, South Dakota – Revisions with data through water year 2008 and simulations of potential future scenarios

The Ogallala and Arikaree aquifers are important water resources in the Rosebud Indian Reservation area and are used extensively for irrigation, municipal, and domestic water supplies. Drought or increased withdrawals from the Ogallala and Arikaree aquifers in the Rosebud Indian Reservation area have the potential to affect water levels in these aquifers. This report documents revisions and recalibration of a previously published three-dimensional, numerical groundwater-flow model for this area. Data for a 30-year period (water years 1979 through 2008) were used in steady-state and transient numerical simulations of groundwater flow. In the revised model, revisions include (1) extension of the transient calibration period by 10 years, (2) the use of inverse modeling for steady-state calibration, (3) model calibration to base flow for an additional four surface-water drainage basins, (4) improved estimation of transient aquifer recharge, (5) improved delineation of vegetation types, and (6) reduced cell size near large capacity water-supply wells. In addition, potential future scenarios were simulated to assess the potential effects of drought and increased groundwater withdrawals. The model comprised two layers: the upper layer represented the Ogallala aquifer and the lower layer represented the Arikaree aquifer. The model’s grid had 168 rows and 202 columns, most of which were 1,640 feet (500 meters) wide, with narrower rows and columns near large water-supply wells. Recharge to the Ogallala and Arikaree aquifers occurs from precipitation on the outcrop areas. The average recharge rates used for the steady-state simulation were 2.91 and 1.45 inches per year for the Ogallala aquifer and Arikaree aquifer, respectively, for a total rate of 255.4 cubic feet per second (ft 3 /s). Discharge from the aquifers occurs through evapotranspiration, discharge to streams as base flow and spring flow, and well withdrawals. Discharge rates for the steady-state simulation were 171.3 ft 3 /s for evapotranspiration, 74.4 ft 3 /s for net outflow to streams and springs, and 11.6 ft 3 /s for well withdrawals. Estimated horizontal hydraulic conductivity used for the numerical model ranged from 0.2 to 84.4 feet per day (ft/d) in the Ogallala aquifer and from 0.1 to 4.3 ft/d in the Arikaree aquifer. A uniform vertical hydraulic conductivity value of 4.2x10 -4 ft/d was estimated for the Ogallala aquifer. Vertical hydraulic conductivity was estimated for five zones in the Arikaree aquifer and ranged from 8.8x10 -5 to 3.7 ft/d. Average rates of recharge, maximum evapotranspiration, and well withdrawals were included in the steady-state simulation, whereas the time-varying rates were included in the transient simulation. Inverse modeling techniques were used for steady-state model calibration. These methods were designed to estimate parameter values that are, statistically, the most likely set of values to result in the smallest differences between simulated and observed hydraulic heads and base-flow discharges. For the steady-state simulation, the root mean square error for simulated hydraulic heads for all 383 wells was 27.3 feet. Simulated hydraulic heads were within ±50 feet of observed values for 93 percent of the wells. The potentiometric surfaces of the two aquifers calculated by the steady-state simulation established initial conditions for the transient simulation. For the transient simulation, the difference between the simulated and observed means for hydrographs was within ±40 feet for 98 percent of 44 observation wells. A sensitivity analysis was used to examine the response of the calibrated steady-state model to changes in model parameters including horizontal and vertical hydraulic conductivity, evapotranspiration, recharge, and riverbed conductance. The model was most sensitive to recharge and maximum evapotranspiration and least sensitive to riverbed and spring conductances. To simulate a potential future drought scenario, a synthetic recharge record was created, the mean of which was equal to 64 percent of the average estimated recharge rate for the 30-year calibration period. This synthetic recharge record was used to simulate the last 20 years of the calibration period under drought conditions. Compared with results of the calibrated model, decreases in hydraulic-head values for the drought scenario at the end of the simulation period were as much as 39 feet for the Ogallala aquifer. To simulate the effects of potential increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 30-year calibration period for the last 20 years of the calibration period. Compared with results of the calibrated model, decreases in hydraulic-head values for the scenario of increased pumping at the end of the simulation period were as much as 13 feet for the Ogallala aquifer. This numerical model is suitable as a tool to help understand the flow system, to help confirm that previous estimates of aquifer properties were reasonable, and to estimate aquifer properties in areas without data. The model also is useful to help assess the effects of drought and increases in pumping by simulations of these scenarios, the results of which are not precise but may be considered when making water management decisions.

South Dakota↗

USGS2009

Numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area, South Dakota

The city of Rapid City and other water users in the Rapid City area obtain water supplies from the Minnelusa and Madison aquifers, which are contained in the Minnelusa and Madison hydrogeologic units. A numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units in the Rapid City area was developed to synthesize estimates of water-budget components and hydraulic properties, and to provide a tool to analyze the effect of additional stress on water-level altitudes within the aquifers and on discharge to springs. This report, prepared in cooperation with the city of Rapid City, documents a numerical groundwater-flow model of the Minnelusa and Madison hydrogeologic units for the 1,000-square-mile study area that includes Rapid City and the surrounding area. Water-table conditions generally exist in outcrop areas of the Minnelusa and Madison hydrogeologic units, which form generally concentric rings that surround the Precambrian core of the uplifted Black Hills. Confined conditions exist east of the water-table areas in the study area. The Minnelusa hydrogeologic unit is 375 to 800 feet (ft) thick in the study area with the more permeable upper part containing predominantly sandstone and the less permeable lower part containing more shale and limestone than the upper part. Shale units in the lower part generally impede flow between the Minnelusa hydrogeologic unit and the underlying Madison hydrogeologic unit; however, fracturing and weathering may result in hydraulic connections in some areas. The Madison hydrogeologic unit is composed of limestone and dolomite that is about 250 to 610 ft thick in the study area, and the upper part contains substantial secondary permeability from solution openings and fractures. Recharge to the Minnelusa and Madison hydrogeologic units is from streamflow loss where streams cross the outcrop and from infiltration of precipitation on the outcrops (areal recharge). MODFLOW-2000, a finite-difference groundwater-flow model, was used to simulate flow in the Minnelusa and Madison hydrogeologic units with five layers. Layer 1 represented the fractured sandstone layers in the upper 250 ft of the Minnelusa hydrogeologic unit, and layer 2 represented the lower part of the Minnelusa hydrogeologic unit. Layer 3 represented the upper 150 ft of the Madison hydrogeologic unit, and layer 4 represented the less permeable lower part. Layer 5 represented an approximation of the underlying Deadwood aquifer to simulate upward flow to the Madison hydrogeologic unit. The finite-difference grid, oriented 23 degrees counterclockwise, included 221 rows and 169 columns with a square cell size of 492.1 ft in the detailed study area that surrounded Rapid City. The northern and southern boundaries for layers 1-4 were represented as no-flow boundaries, and the boundary on the east was represented with head-dependent flow cells. Streamflow recharge was represented with specified-flow cells, and areal recharge to layers 1-4 was represented with a specified-flux boundary. Calibration of the model was accomplished by two simulations: (1) steady-state simulation of average conditions for water years 1988-97 and (2) transient simulations of water years 1988-97 divided into twenty 6-month stress periods. Flow-system components represented in the model include recharge, discharge, and hydraulic properties. The steady-state streamflow recharge rate was 42.2 cubic feet per second (ft3/s), and transient streamflow recharge rates ranged from 14.1 to 102.2 ft3/s. The steady-state areal recharge rate was 20.9 ft3/s, and transient areal recharge rates ranged from 1.1 to 98.4 ft3/s. The upward flow rate from the Deadwood aquifer to the Madison hydrogeologic unit was 6.3 ft3/s. Discharge included springflow, water use, flow to overlying units, and regional outflow. The estimated steady-state springflow of 32.8 ft3/s from seven springs was similar to the simulated springflow of 31.6 ft3/s, which included 20.5 ft3

South Dakota↗

USGS2008

Algal and water-quality data for Rapid Creek and Canyon Lake near Rapid City, South Dakota, 2007

This report summarizes the results of algae and water-quality sampling on Rapid Creek and Canyon Lake during May and September 2007. The overall purpose of the study was to determine the algal community composition of Rapid Creek and Canyon Lake in relation to organisms that are known producers of unwanted tastes and odors in drinking-water supplies. Algal assemblage structure (phytoplankton and periphyton) was examined at 16 sites on Rapid Creek and Canyon Lake during May and September 2007, and actinomycetes bacteria were sampled at the Rapid City water treatment plant intake in May 2007, to determine if taste-and-odor producing organisms were present. During the May 2007 sampling, 3 Rapid Creek sites and 4 Canyon Lake sites were quantitatively sampled for phytoplankton in the water column, 7 Rapid Creek sites were quantitatively sampled for attached periphyton, and 4 lake and retention pond sites were qualitatively sampled for periphyton. Five Rapid Creek sites were sampled for geosmin and 2-methylisoborneol, two common taste-and-odor causing compounds known to affect water supplies. During the September 2007 sampling, 4 Rapid Creek sites were quantitatively sampled for attached periphyton, and 3 Canyon Lake sites were qualitatively sampled for periphyton. Water temperature, dissolved oxygen, pH, and specific conductance were measured during each sampling event. Methods of collection and sample analysis are presented for the various types of biological and chemical constituent samples. Diatoms comprised 91-100 percent of the total algal biovolume in periphyton samples collected during May and September. Cyanobacteria (also called blue-green algae) were detected in 7 of the 11 quantitative periphyton samples and ranged from 0.01 to 2.0 percent of the total biovolume. Cyanobacteria were present in 3 of the 7 phytoplankton samples collected in May, but the relative biovolumes were small (0.01-0.2 percent). Six of seven qualitative samples collected from Canyon Lake and retention ponds during May and September also contained cyanobacteria. Geosmin and 2-methylisoborneol concentrations were less than detection limits (0.005 ug/L) in all five of the Rapid Creek samples collected in May. Actinomycetes bacteria were present at the water treatment plant intake in May 2007, at a concentration of 6 colonies per milliliter. During this study, no taste-and-odor problems with the drinking water within the study area were reported. However, the presence of cyanobacterial taxa known to contain taste-and-odor producing strains (such as Leptolyngbya, Phormidium, and Anabaena) indicates the potential for taste-and-odor problems under certain physical and chemical conditions.

South Dakota↗

USGS2008

Water-quality effects and characterization of indicators of onsite wastewater disposal systems in the east-central Black Hills area, South Dakota, 2006-08

Onsite wastewater disposal systems (OWDS) are used extensively in the Black Hills of South Dakota where many of the watersheds and aquifers are characterized by fractured or solution-enhanced bedrock with thin soil cover. A study was conducted during 2006-08 to characterize water-quality effects and indicators of OWDS. Water samples were collected and analyzed for potential indicators of OWDS, including chloride, bromide, boron, nitrite plus nitrate (NO2+NO3), ammonia, major ions, nutrients, selected trace elements, isotopes of nitrate, microbiological indicators, and organic wastewater compounds (OWCs). The microbiological indicators were fecal coliforms, Escherichia coli (E. coli), enterococci, Clostridium perfringens (C. perfringens), and coliphages. Sixty ground-water sampling sites were located either downgradient from areas of dense OWDS or in background areas and included 25 monitoring wells, 34 private wells, and 1 spring. Nine surface-water sampling sites were located on selected streams and tributaries either downstream or upstream from residential development within the Precambrian setting. Sampling results were grouped by their hydrogeologic setting: alluvial, Spearfish, Minnekahta, and Precambrian. Mean downgradient dissolved NO2+NO3 concentrations in ground water for the alluvial, Spearfish, Minnekahta, and Precambrian settings were 0.734, 7.90, 8.62, and 2.25 milligrams per liter (mg/L), respectively. Mean downgradient dissolved chloride concentrations in ground water for these settings were 324, 89.6, 498, and 33.2 mg/L, respectively. Mean downgradient dissolved boron concentrations in ground water for these settings were 736, 53, 64, and 43 micrograms per liter (ug/L), respectively. Mean dissolved surface-water concentrations for NO2+NO3, chloride, and boron for downstream sites were 0.222 mg/L, 32.1 mg/L, and 28 ug/L, respectively. Mean values of delta-15N and delta-18O (isotope ratios of 14N to 15N and 18O to 16O relative to standard ratios) for nitrate in ground-water samples were 10.4 and -2.0 per mil (0/100), respectively, indicating a relatively small contribution from synthetic fertilizer and probably a substantial contribution from OWDS. The surface-water sample with the highest dissolved NO2+NO3 concentration of 1.6 mg/L had a delta-15N value of 12.36 0/100, which indicates warm-blooded animals (including humans) as the nitrate source. Fecal coliforms were detected in downgradient ground water most frequently in the Spearfish (19 percent) and Minnekahta (9.7 percent) settings. E. coli was detected most frequently in the Minnekahta (29 percent) and Spearfish (13 percent) settings. Enterococci were detected more frequently than other microbiological indicators in all four settings. Fecal coliforms and E. coli were detected in 73 percent and 95 percent of all surface-water samples, respectively. Enterococci, coliphages (somatic), and C. perfringens were detected in 50, 70, and 50 percent of surface-water samples, respectively. Of the 62 OWC analytes, 12 were detected only in environmental samples, 10 were detected in at least one environmental and one blank sample (not necessarily companion pairs), 2 were detected only in blank samples, and 38 were not detected in any blank, environmental, or replicate sample from either ground or surface water. Eleven different organic compounds were detected in ground-water samples at eight different sites. The most frequently occurring compound was DEET, which was found in 32 percent of the environmental samples, followed by tetrachloroethene, which was detected in 20 percent of the samples. For surface-water samples, 16 organic compounds were detected in 9 of the 10 total samples. The compound with the highest occurrence in surface-water samples was camphor, which was detected in 50 percent of samples. The alluvial setting was characterized by relatively low dissolved NO2+NO3 concentrations, detection of ammonia nitrogen, and relatively high concentr

South Dakota↗

USGS2007

Characterization of ground-water flow and water quality for the Madison and Minnelusa aquifers in northern Lawrence County, South Dakota

The Madison and Minnelusa aquifers are used extensively for water supplies for the city of Spearfish and other users in northern Lawrence County, South Dakota. Ground water in the Madison and Minnelusa aquifers in the study area generally flows north from outcrop areas where recharge from sinking streams and infiltration of precipitation occurs. Ground water that moves northward and eastward around the Black Hills enters the study area from the west and results in hydraulic heads that are several hundred feet higher on the western side of the study area than on the eastern side. The estimated average recharge rate of 38 cubic feet per second (ft3/s) on outcrops of the Madison Limestone and Minnelusa Formation is less than the total estimated average spring discharge rate of 51 ft3/s in the northwestern part of the study area. Sixteen pounds of fluorescein dye were injected into Spearfish Creek on March 25, 2003, when streamflow was 6.6 ft3/s. The dye was detected in water samples from four wells completed in the Madison aquifer ranging from 2.6 to 4.5 miles north of the injection site. First arrival times ranged from 5 to 169 days, and ground-water velocities ranged from about 0.1 to 0.5 mile per day. Sixty-four pounds of Rhodamine WT was injected into Spearfish Creek at the same location on May 9, 2003, when streamflow was 5.6 ft3/s. Rhodamine WT dye concentrations measured in samples from the same four wells were about an order of magnitude less than measured fluorescein concentrations. Oxygen- and deuterium-isotope values for samples from Cox Lake and McNenny Pond springs indicated a probable component of spring discharge that originates from outcrops of the Madison Limestone and Minnelusa Formation on the Limestone Plateau south of the study area. Oxygen- and deuterium-isotope values for samples from Mirror Lake spring indicated possible contributions from overlying aquifers and local recharge. Oxygen- and deuterium-isotope values for the combined springflow contributing to Crow Creek in the northwestern part of the study area indicated that the primary source of water is the Madison and Minnelusa aquifers. Oxygen- and deuterium-isotope values for Old Hatchery and Higgins Gulch springs, located north of Spearfish, indicated a source water originating from the outcrops of the Madison Limestone and Minnelusa Formation within the study area. Concentrations of three chlorofluorocarbons (CFC-11, CFC-12, and CFC-113) were used to characterize ground-water residence times in the study area. For the four wells where dye was detected, CFC-11 apparent ages ranged from 12 to 26 years, indicating that the wells contained months-old water mixed with years- to decades-old water. Logarithmic regression analysis of the CFC-11 apparent ages for water from 10 wells and distance to a possible conduit trending north through the area where dye was detected, yielded an r2 value of 0.71. Straight-line regression analysis of the CFC-11 apparent ages for the six wells closest to the possible conduit had an r2 value of 0.96. Two wells located relatively close to the outcrop areas had no or very low tritium values indicating relatively long residence times and diffuse ground-water flow. The tritium value of 7.2 TU in water from well COL where dye was detected, indicated that the water probably is a bimodal mixture, with a substantial portion that is older than 50 years. Water from well ELL, where dye was detected, had a tritium value of 19.7 TU and a CFC apparent age of 15 years, indicating that the sample from this well probably is a unimodal mixture with very little water older than 50 years. Comparison of the CFC apparent age for three spring sites (Cox Lake, 26 years; McNenny Pond, 26 years; Mirror Lake, 13 years) also indicated that Mirror Lake spring probably has a component of local recharge from formations that overlie the Minnelusa Formation. In the Madison aquifer, specific conductance ranges from 18 to 945 microsiemens per cen

South Dakota↗

USGS2007

Analysis of ground-water flow in the Madison aquifer using fluorescent dyes injected in Spring Creek and Rapid Creek near Rapid City, South Dakota, 2003-04

The Madison aquifer, which contains fractures and solution openings in the Madison Limestone, is used extensively for water supplies for the city of Rapid City and other suburban communities in the Rapid City, S. Dak., area. The 48 square-mile study area includes the west-central and southwest parts of Rapid City and the outcrops of the Madison Limestone extending from south of Spring Creek to north of Rapid Creek. Recharge to the Madison Limestone occurs when streams lose flow as they cross the outcrop. The maximum net loss rate for Spring and Rapid Creek loss zones are 21 and 10 cubic feet per second (ft3/s), respectively. During 2003 and 2004, fluorescent dyes were injected in the Spring and Rapid Creek loss zones to estimate approximate locations of preferential flow paths in the Madison aquifer and to measure the response and transit times at wells and springs. Four injections of about 2 kilograms of fluorescein dye were made in the Spring Creek loss zone during 2003 (sites S1, S2, and S3) and 2004 (site S4). Injection at site S1 was made in streamflow just upstream from the loss zone over a 12-hour period when streamflow was about equal to the maximum loss rate. Injections at sites S2, S3, and S4 were made in specific swallow holes located in the Spring Creek loss zone. Injection at site R1 in 2004 of 3.5 kilograms of Rhodamine WT dye was made in streamflow just upstream from the Rapid Creek loss zone over about a 28-hour period. Selected combinations of 27 wells, 6 springs, and 3 stream sites were monitored with discrete samples following the injections. For injections at sites S1-S3, when Spring Creek streamflow was greater than or equal to 20 ft3/s, fluorescein was detected in samples from five wells that were located as much as about 2 miles from the loss zone. Time to first arrival (injection at site S1) ranged from less than 1 to less than 10 days. The maximum fluorescein concentration (injection at site S1) of 120 micrograms per liter (ug/L) at well CO, which is located adjacent to the loss zone, was similar to the concentration in the stream. Fluorescein arrived at well NON (injection at site S1), which is located about 2 miles northeast of the loss zone, within about 1.6 days, and the maximum concentration was 44 ug/L. For injection at site S4, when streamflow was about 12 ft3/s, fluorescein was detected in samples from six wells and time to first arrival ranged from 0.2 to 16 days. Following injection at site S4 in 2004, the length of time that dye remained in the capture zone of well NON, which is located approximately 2 miles from the loss zone, was almost an order of magnitude greater than in 2003. For injection at site R1, Rhodamine WT was detected at well DRU and spring TI-SP with time to first arrival of about 0.5 and 1.1 days and maximum concentrations of 6.2 and 0.91 ug/L, respectively. Well DRU and spring TI-SP are located near the center of the Rapid Creek loss zone where the creek has a large meander. Measurable concentrations were observed for spring TI-SP as many as 109 days after the dye injection. The direction of a conduit flow path in the Spring Creek area was to the northeast with ground-water velocities that ranged from 770 to 6,500 feet per day. In the Rapid Creek loss zone, a conduit flow path east of the loss zone was not evident from the dye injection.

South Dakota↗

USGS2003

Simulated ground-water flow in the Ogallala and Arikaree aquifers, Rosebud Indian Reservation area, South Dakota

South Dakota↗

USGS2002

Flow-system analysis of the Madison and Minnelusa aquifers in the Rapid City area, South Dakota — Conceptual model

The conceptual model of the Madison and Minnelusa aquifers in the Rapid City area synthesizes the physical geography, hydraulic properties, and ground-water flow components of these important aquifers. The Madison hydrogeologic unit includes the karstic Madison aquifer, which is defined as the upper, more permeable 100 to 200 ft of the Madison Limestone, and the Madison confining unit, which consists of the lower, less permeable part of the Madison Limestone and the Englewood Formation. Overlying the Madison hydrogeologic unit is the Minnelusa hydrogeologic unit, which includes the Minnelusa aquifer in the upper, more permeable 200 to 300 ft and the Minnelusa confining unit in the lower, less permeable part. The Madison and Minnelusa hydrogeologic units outcrop in the study area on the eastern flank of the Black Hills where recharge occurs from streamflow losses and areal recharge. The conceptual model describes streamflow recharge, areal recharge, ground-water flow, storage in aquifers and confining units, unsaturated areas, leakage between aquifers, discharge from artesian springs, and regional outflow. Effective transmissivities estimated for the Madison aquifer range from 500 to 20,000 ft 2 /d and for the Minnelusa aquifer from 500 to 10,000 ft 2 /d. Localized anisotropic transmissivity in the Madison aquifer has tensor ratios as high as 45:1. Vertical hydraulic conductivities for the Minnelusa confining unit determined from aquifer tests range from 1.3x10 -3 to 3.0x10 -1 ft/d. The confined storage coefficient of the Madison and Minnelusa hydrogeologic units was estimated as 3x10 -4 ft/d. Specific yield was estimated as 0.09 for the Madison and Minnelusa aquifers and 0.03 for the Madison and Minnelusa confining units. Potentiometric surfaces for the Madison and Minnelusa aquifers have a general easterly gradient of about 70 ft/mi with local variations. Temporal change in hydraulic head in the Madison and Minnelusa aquifers ranged from about 5 to 95 ft in water years 1988-97. The unconfined areas were estimated at about 53 and 36 mi 2 for the Madison and Minnelusa hydrogeologic units, respectively, in contrast to an aquifer analysis area of 629 mi 2 . Dye-tracer tests, stable isotopes, and hydrogeologic features were analyzed conjunctively to estimate generalized ground-water flowpaths in the Madison aquifer and their influences on the Minnelusa aquifer. The western Rapid City area between Boxelder Creek and Spring Creek was characterized as having undergone extensive tectonic activity, greater brecciation in the Minnelusa Formation, large transmissivities, generally upward hydraulic gradients from the Madison aquifer to the Minnelusa aquifer, many karst springs, and converging flowpaths. Water-budget analysis included: (1) a dry-period budget for declining water levels; October 1, 1987, to March 31, 1993; (2) a wet-period budget for rising water levels, April 1, 1993, to September 30, 1997; and (3) a full 10-year period budget for water years 1988-97. By simultaneously balancing these water budgets, initial estimates of recharge, discharge, change in storage, and hydraulic properties were refined. Inflow rates for the 10-year budget included streamflow recharge of about 45 ft 3 /s or 61 percent of the total budget and areal recharge of 22 ft 3 /s or 30 percent. Streamflow recharge to the Madison hydrogeologic unit was about 86 percent of the total streamflow recharge. Outflow for the 10-year budget included springflow of 31 ft 3 /s or 42 percent of the total budget, water use of about 10 ft 3 /s or 14 percent, and regional outflow of 22 ft 3 /s or 30 percent. Ground-water storage increased 9 ft 3 /s during the 10-year period, and net ground-water movement from the Madison to Minnelusa hydrogeologic unit was about 8 ft 3 /s.

South Dakota↗

USGS2001

Hydrologic setting and geochemical characterization of free-phase hydrocarbons in the alluvial aquifer at Mandan, North Dakota, November 2000

Free-phase hydrocarbons are present in the alluvial aquifer at Mandan, North Dakota. A large contaminant body of the hydrocarbons [light nonaqueous phase liquid (LNAPL)] floats on the water table about 20 feet below land surface. The main LNAPL body is about 6 feet thick, and the areal extent is about 657,000 square feet. A study was conducted to describe the hydrologic setting and characterize the geochemical composition of the free-phase hydrocarbons in the alluvial aquifer. Most of the study area is underlain by alluvium of the Heart River Valley that ranges in thickness from about 25 to 109 feet. The alluvium can be divided into three stratigraphic units silty clay, silty sand, and sand and is underlain by shales and sandstones. Monitoring wells were installed prior to this study, to an average depth of about 29 feet. Regional ground-water flow in the Heart River aquifer generally may be from west-northwest to eastsoutheast and is influenced by hydraulic connections to the river. Hydraulic connections also are probable between the aquifer and the Missouri River. Ground-water flow across the north boundary of the aquifer is minimal because of adjacent shales and sandstones of relatively low permeability. Recharge occurs from infiltration of precipitation and is spatially variable depending on the thickness of overlying clays and silts. Although the general water-table gradient may be from west-northwest to east-southeast, the flow directions can vary depending on the river stage and recharge events. Any movement of the LNAPL is influenced by the gradients created by changes in water-level altitudes. LNAPL samples were collected from monitoring wells using dedicated bailers. The samples were transferred to glass containers, stored in the dark, and refrigerated before shipment for analysis by a variety of analytical techniques. For comparison purposes, reference-fuel samples provided by the refinery in Mandan also were analyzed. These reference-fuel samples included a current diesel fuel, a closely related but slightly broader refinery-cut fuel, a crude-oil composite, unleaded regular gasoline, and additives. Four principal analytical techniques were used for geochemical characterization: Purge-and-trap gas chromatography/mass spectrometry (volatile components); capillary gas chromatography/mass spectrometry (semivolatile components); isotope ratio mass spectrometry (carbon isotopes; whole oils); and liquid chromatography/mass spectrometry with electrospray ionization (additives and other organic components). Volatile analytes included solvents, disinfection byproducts, halogenated hydrocarbons, and alkylbenzenes, including benzene, toluene, ethylbenzene, and meta-, para-, and orf/zo-xylenes. Semivolatile analytes included rt-alkanes, isoprenoid alkanes, cycloalkanes, and polycyclic aromatic hydrocarbons and related compounds (naphthalenes, phenanthrenes, and dibenzothiophenes and their alkylated derivatives). Of the additives, only the diesel-fuel additive with the red dye marker was amenable to electrospray ionization. Results indicate the LNAPL consists of closely correlatable diesel fuel at various stages of degradation. All LNAPL samples contained the red dye marker for diesel fuel. None of the samples contained chlorinated solvents associated with industries such as drycleaning or automotive maintenance. Solvents such as acetone, dimethyl ether, and methylene chloride and the gasoline additives methyl-t-butyl ether (MTBE), ethyl-t-butyl ether (ETBE), and t-amyl-methyl ether (TAME) were not found. With one possible exception, no evidence of a different diesel or other hydrocarbon fuel contribution was identified. At one site near the north edge of the main LNAPL body, evidence exists for traces of possible gasoline components in addition to the diesel fuel. The geochemical analysis of the LNAPL and correlations with other fuel products and additives strongly suggest episodic releases of a single, local-source, diesel fuel into the aquifer over an extended period of time.

North Dakota↗

USGS2000

Sensitivity of ground water to contamination in Lawrence County, South Dakota

Ground-water supplies in Lawrence County, South Dakota, can be contaminated by agricultural, urban, suburban, commercial, and industrial land uses. To address this issue, the U.S. Geological Survey in cooperation with Lawrence County and the City of Spearfish mapped the sensitivity of ground water to contamination in Lawrence County. Sensitivity of ground water to contamination was determined by delineating hydrogeologic settings with common hydrogeologic characteristics as described in the DRASTIC method, developed by the U.S. Environmental Protection Agency and the National Water Well Association. Within the framework of 11 hydrogeologic settings, sensitivity to contamination was ranked for six intrinsic hydrogeologic characteristics: (1) aquifer media, (2) unsaturated media (3) hydraulic conductivity, (4) recharge rate, (5) depth to water, and (6) land-surface slope. The rating conventions of DRASTIC were modified to provide a relative ranking of hydrogeologic characteristics without assignment of a combined numerical score. Soil characteristics were not included as a map layer because detailed digital data were not available; however, the general distribution of two soil characteristics were shown. A total of 956 polygons were delineated and assigned a sensitivity-unit code that represented unique groups of sensitivity rank for the six intrinsic hydrogeologic characteristics. The polygons were created by overlaying and intersecting maps that describe the geology, precipitation, land-surface elevation, and depth to water using a geographic information system. Thirty drainage areas upstream from potential streamflow-loss zones were delineated to describe an additional mechanism of transport of potential contamination. The sensitivity of ground water to contamination was presented on a 1:100,000-scale map with code and label explanations. Limitations of the sensitivity map are described to facilitate appropriate use of the map as a screening tool to compare sensitivity to contamination.

Water-Resources Investigations Report↗