Geology Reports⌕ Search

Geology topics

Daniel G. Driscoll

Publications and source records attributed to Daniel G. Driscoll.

At least 19 recordsLinked to original sources
USGS2021

Characterization of factors affecting groundwater levels in and near the former Lake Traverse Indian Reservation, South Dakota, water years 1956–2017

The U.S. Geological Survey (USGS), in cooperation with the Sisseton Wahpeton Oyate, completed a study to characterize water-level fluctuations in observation wells relative to driving factors that affect water levels in and near the historical 1867 boundary of the Lake Traverse Indian Reservation. The study investigated concerns regarding potential effects of groundwater withdrawals and climate conditions on groundwater levels within an area that includes the historical boundary of the reservation and a surrounding area that extends 10 miles in all directions within South Dakota. Characterization of water-level fluctuations in observation wells and relative driving factors was accomplished by statistical trend analysis. Monthly data from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) were aggregated to obtain annual and seasonal datasets for total precipitation, minimum air temperature ( T min ), and maximum air temperature ( T max ) for the study area and a surrounding buffer area. Trend tests for gridded data for total precipitation, T min , and T max were completed for annual and seasonal time series for water years 1956–2017, which is about 2 years before the earliest available water-level measurements. A 2-year offset was arbitrarily selected because scrutiny of water-level and precipitation data indicated that responses of groundwater levels for many of the observation wells lagged major changes in precipitation patterns by about 2 years. Statistically significant upward trends were detected for annual precipitation and annual T min for most of the study area and the surrounding buffer area. Statistically significant downward trends in T max were detected for only a few 2.5 arc-minute grid cells; however, the sparsity of the spatial coverage reduces confidence that these are true trends, in contrast to the near completeness of the spatial coverage in upward trends for T min . Spatial distributions of statistically significant trends in seasonal climate data were generally similar to the annual trends, but with substantial differences in the spatial density of the trends. Potential interactions among water levels in observation wells and streamflow were examined through correlation analyses of the annual median water level for each of 76 observation wells versus the annual mean streamflow for each of four area streamgages. Potential interactions among water levels in observation wells and lake levels were examined through correlation analyses involving 25 area lakes. Resulting correlation coefficients were used as part of an approach for selecting a lake to be plotted in conjunction with water-level and precipitation data for each observation well. Groundwater trends for 76 observation wells were analyzed for three separate water-level parameters (minimum, median, and maximum) because wells are measured sporadically, and data are biased towards more frequent measurements during periods of heaviest irrigation demand. Trends in the time series of annual precipitation (from PRISM) starting 2 years earlier than the associated water-level trend also were analyzed for the location of each individual observation well. Sen’s slope and Mann-Kendall p -values were computed for the three water-level parameters and for the annual precipitation time series. Graphs showing results of trend analyses for each observation well also showed changes with time in the sum of licensed groundwater withdrawals within six specified radii (0.5, 1.0, 2.0, 3.0, 4.0, and 5.0 miles) of each well as a qualitative indicator of proximal groundwater demand. Trends in groundwater levels in observation wells in the study area are predominantly upward, with 43 of 76 wells having significant upward trends for at least one of the three water-level parameters and only 8 wells having significant downward trends for at least one water-level parameter. The upward groundwater trends are driven by predominantly upward precipitation trends, with 43 wells (not all the same wells) also having significant upward trends and no wells having significant downward trends. Significant upward precipitation trends were detected for only two of the eight wells with significant downward groundwater trends. Groundwater levels in some observation wells likely are also substantially affected by interactions with surface water, especially with lakes. Water levels in many area lakes increased in response to wet conditions of the early 1990s and have maintained high water levels ever since. It is recognized that in many cases lakes that were selected for plotting with groundwater hydrographs likely are not hydraulically connected with a groundwater system or aquifer associated with an individual well; however, interactions also are plausible for numerous other lakes for which water-level records are not available.

South Dakota↗
USGS2021

Trends in groundwater levels in and near the Rosebud Indian Reservation, South Dakota, water years 1956–2017

The U.S. Geological Survey (USGS), in cooperation with the Rosebud Sioux Tribe, completed a study to characterize water-level fluctuations in observation wells to examine driving factors that affect water levels in and near the Rosebud Indian Reservation, which comprises all of Todd County. The study investigates concerns regarding potential effects of groundwater withdrawals and climate conditions on groundwater levels within an area that includes Todd County and a surrounding area that extends 10 miles north, east, and west of the county border. Characterization of water-level fluctuations in observation wells and relative driving factors was accomplished by statistical trend analysis. Two statistical methods were used for analysis of temporal trends for climatic and hydrologic data. To determine which trend analysis to use, applicable datasets were tested for statistically significant short-term persistence (STP). In the absence of significant STP, existence of statistical trends was determined using the standard Mann-Kendall test for probability values less than or equal to 0.10 (90-percent confidence level); however, a modified Mann-Kendall test was used for datasets where statistically significant STP was detected. Trend magnitudes were computed using the Sen’s slope estimator. Monthly data from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) were aggregated to obtain annual and seasonal datasets for total precipitation, minimum air temperature ( T min ), and maximum air temperature ( T max ) for the study area and a surrounding buffer area. Trend tests for total precipitation, T min , and T max were completed for annual and seasonal time series for water years 1956–2017, which is about 2 years before the earliest available water-level measurements. A 2-year offset was arbitrarily selected because scrutiny of water-level and precipitation data indicated that responses of groundwater levels for many of the observation wells lagged major changes in precipitation patterns by about 2 years. Statistically significant upward trends were detected for annual precipitation and annual T min for almost all of the study area and the surrounding buffer area. Statistically significant downward trends in T max were detected for a very small part of the study area; however, the sparse spatial coverage reduces confidence that these are true trends. Spatial distributions of statistically significant trends in seasonal climate data were generally similar to the annual trends, but with substantial differences in the spatial density of the trends. Groundwater trends for 58 observation wells were analyzed for three separate water-level parameters (minimum, median, and maximum) because wells are measured sporadically and data are biased towards more frequent measurements during periods of heaviest irrigation demand. Trends in the time series of annual precipitation (from PRISM) starting 2 years earlier than for the associated water-level trend also were analyzed for the location of each individual observation well. Sen’s slope and Mann-Kendall probability values (p-values) were computed for the three water-level parameters and for the annual precipitation time series. Graphs showing results of trend analyses for each observation well also showed changes over time in the sum of licensed groundwater withdrawals within six specified radii (0.5, 1, 2, 3, 4, and 5 miles) of each well as a qualitative indicator of proximal groundwater demand. Of all 58 observation wells considered, 28 wells had significant upward trends for at least one of the three water-level parameters, 11 wells had significant downward trends for at least one water-level parameter, and 19 wells did not have any significant trends. Significant upward trends in annual precipitation were detected for 48 of the 58 wells. Results of trend analyses likely show the effects of groundwater withdrawals on water levels in the Ogallala aquifer in areas of substantial demand. Precipitation trends are significantly upward for 43 of the 48 wells completed in the Ogallala aquifer that were analyzed. Of the 48 Ogallala aquifer wells, 24 had significant upward trends for at least one water-level parameter (17 with all 3); however, 10 wells had statistically significant downward trends for at least one water-level parameter (8 with all 3 parameters). All but one of the wells with significant downward trends are located in the south-central part of the study area where licensed irrigation withdrawals are concentrated.

South Dakota↗
USGS2019

Climate, streamflow, and lake-level trends in the Great Lakes Basin of the United States and Canada, water years 1960–2015

Water levels in the Great Lakes fluctuate substantially because of complex interactions among inputs (precipitation and streamflow), outputs (evaporation and outflow), and other factors. This report by the U.S. Geological Survey in cooperation with the Great Lakes Restoration Initiative was completed to describe trends in climate, streamflow, lake levels, and major water-budget components within the Great Lakes Basin for water years (WYs) 1960–2015 (study period). Resulting trends are applicable only to the study period and should not be considered indicative of longer-term trends. Analyses of climate trends used monthly data from the Parameter-elevation Regressions on Independent Slopes Model, which are available only for the United States. Trend tests were completed for annual and seasonal time series of monthly means for total precipitation, daily minimum air temperature ( T min ), and daily maximum air temperature ( T max ). Statistical significance for all time-trend tests (climate, streamflow, and lake levels) was determined using the Mann‑Kendall test for probability values less than or equal to 0.10. Trend analyses were completed without adjustments for serial correlation; however, a modified Mann-Kendall test was subsequently used to examine potential effects of short-term persistence in time-series data. Effects of short-term persistence were considered inconsequential for climate data and minor for streamflow data; however, the presence of short-term persistence in water-budget components had more substantial effects on trend analyses. Spatial distributions of trends in climatic data for WYs 1960–2015 for the U.S. part of the Great Lakes Basin (land only) indicate (1) generally ubiquitous upward trends in T min and (2) a sharp transition from neutral or downward trends in precipitation northwest of Lake Michigan to generally upward trends east of Lake Michigan. Trends in T max were not statistically significant. Analyses of annual climatic data aggregated for the U.S. land part of the Great Lakes Basin indicated statistically significant upward trends for precipitation and T min , and similar statistically significant trends existed for all the individual lake subbasins except Lake Superior. Of 103 U.S. Geological Survey streamgages analyzed for streamflow trends, 71 had significant annual trends (54 upward and 17 downward). Downward trends in annual streamflow are concentrated northwest of Lake Michigan (16 streamgages), and upward trends are concentrated east of Lake Michigan (53 streamgages). Of the 71 streamgages with significant annual trends, 70 had at least one season with a significant trend that matched the annual trend direction. Of 35 Environment and Climate Change Canada streamgages analyzed, 22 had significant upward trends in annual streamflow, and all but 1 of these 22 had at least one season with a significant upward trend. None of the Environment and Climate Change Canada streamgages had significant downward annual trends, and only one had a significant downward seasonal trend. Trends in lake levels and several major water-budget components affecting lake levels were analyzed for the study period. Significant downward trends in lake level and outflow for Lake Superior are driven primarily by low lake levels and outflows during WYs 1998–2014. A significant downward trend in runoff from the contributing drainage area also is indicated, which is consistent with numerous streamgages northwest of Lake Michigan with significant downward trends in annual streamflow. A significant upward trend in annual overlake evaporation also is indicated, which is consistent with the spatially distributed upward trends in annual T min . The sum of overlake precipitation and runoff from the contributing drainage area for each of the Great Lakes, less overlake evaporation, composes a variable called net basin supply (NBS). A significant downward trend in NBS is indicated for Lake Superior, which is consistent with significant trends for individual components of runoff (downward) and evaporation (upward) that contributed to a significant downward trend for lake outflow. Statistically significant upward trends in NBS for Lake Saint Clair and Lake Ontario offset the downward trend for Lake Superior and combine with nonsignificant upward trends in NBS for Lakes Michigan and Huron and Lake Erie to produce a neutral trend in NBS for the basin. A predictable pattern in monthly mean lake levels is noted for Lake Superior, with the minimum for each year usually during or near March and the maximum commonly during or near September or October. When an October lake level is in a period of substantial decline, potential for an ensuing short-term period of below-mean lake levels is enhanced. Downstream from Lake Superior, monthly lake levels have sawtooth patterns that somewhat resemble those for Lake Superior but with decreased predictability in timing. Similar to Lake Superior, Lakes Michigan and Huron, Lake Saint Clair, and Lake Erie all have a prolonged period of low lake levels around WYs 1998–2014; however, a significant downward trend is indicated only for Lakes Michigan and Huron. All these lakes also have a period of low lake levels before about WY 1968, when minimum lake levels were lower than during WYs 1998–2014. The significant downward trend of outflow from Lake Superior is carried downstream into Lakes Michigan and Huron; however, trends in outflow from the next three lakes downstream (Lakes Saint Clair, Erie, and Ontario) are offset by increased precipitation and runoff and are not significant.

Scientific Investigations Report↗
USGS2017

Application of paleoflood surveys for the southern Black Hills of South Dakota

Flood-frequency analyses for the Black Hills area have especially large uncertainties and are especially important for planning purposes because of a history of extremely large and damaging floods, such as the extreme floods of June 9–10, 1972. Geology, topography, and climatology are additional complicating factors for flood-frequency characterization for the area. Two previous paleoflood studies for the Black Hills area indicated good potential for improving flood-frequency analyses through implementation of paleoflood investigations. The objectives of this study (SD2010-04) for the southern Black Hills were to (1) develop long-term flood chronologies and associated peak-flow frequency analyses for selected stream reaches by applying paleoflood hydrology approaches; and (2) develop flood-frequency information regarding “high-elevation” stream reaches to help address questions regarding differential potential for generation of exceptionally strong rain-producing thunderstorms across elevation gradients in the area. Neither objective was accomplished because the study was terminated before planned completion. Substantial efforts could be applied for only 2 of the 12 research tasks prior to study termination. These were task 3 (preliminary reconnaissance) and task 4 (activities associated with Section 106 of the National Historic Preservation Act). Field reconnaissance conducted along 10 candidate streams indicated that conditions in the southern Black Hills appear quite favorable for conducting paleoflood investigations. All 10 candidate streams had moderate to good potential for favorable paleoflood evidence, and in general are well constrained in relatively narrow canyon reaches, which provides good sensitivity for changes in stage, relative to discharge. Task 4 (Section 106 activities) was needed because alcoves and rock shelters that are well suited for deposition and preservation of paleoflood evidence may have been used as shelters or cache locations by indigenous inhabitants and thus may be eligible for consideration as historic properties because of possible archaeological or cultural materials. The complexity of the Section 106 concerns and issues became progressively more apparent as the study evolved. The study eventually was terminated when it became apparent that the resources needed to address the Section 106 issues would overwhelm the resources available for study implementation. In the event of consideration of future re-implementation, approaches that might help expedite Section 106 issues could include (1) a partnership with another Federal agency that has substantial experience with the Section 106 process; (2) securing assistance from a consultant that could help with both the Section 106 process and the required archaeological component; and (3) partnering with a tribal college with archaeological or earth science/hydrology programs, which could help make this study become part of a learning exercise.

South Dakota↗
USGS2015

Late Holocene flood probabilities in the Black Hills, South Dakota with emphasis on the Medieval Climate Anomaly

A stratigraphic record of 35 large paleofloods and four large historical floods during the last 2000 years for four basins in the Black Hills of South Dakota reveals three long-term flooding episodes, identified using probability distributions, at A.D.: 120–395, 900–1290, and 1410 to present. During the Medieval Climate Anomaly (~ A.D. 900–1300) the four basins collectively experienced 13 large floods compared to nine large floods in the previous 800 years, including the largest floods of the last 2000 years for two of the four basins. This high concentration of extreme floods is likely caused by one or more of the following: 1) instability of air masses caused by stronger than normal westerlies; 2) larger or more frequent hurricanes in the Gulf of Mexico and Atlantic Ocean; and/or 3) reduced land covering vegetation or increased forest fires caused by persistent regional drought.

South Dakota↗
USGS2015

Hydrologic response for a high-elevation storm in the South Dakota Black Hills

A group of thunderstorms produced >4 in of rain during four periods of progressively more intense rainfall across a small part of a relatively high-elevation area of the northern Black Hills on 5 August 2014. The resulting hydrologic response was noteworthy in two very small headwater drainage basins, where the measured peak flows are by far the largest—relative to drainage area—ever documented for the high-elevation Limestone Plateau area. However, peak flows attenuated quickly in a downstream direction owing to the storms tracking perpendicular to the drainage direction, moderately dry antecedent conditions, and progressive widening of the valley bottoms.

South Dakota↗
USGS2014

Annual exceedance probabilities and trends for peak streamflows and annual runoff volumes for the Central United States during the 2011 floods

During 2011, excess precipitation resulted in widespread flooding in the Central United States with 33 fatalities and approximately $4.2 billion in damages reported in the Red River of the North, Souris, and Mississippi River Basins. At different times from late February 2011 through September 2011, various rivers in these basins had major flooding, with some locations having multiple rounds of flooding. This report provides broadscale characterizations of annual exceedance probabilities and trends for peak streamflows and annual runoff volumes for selected streamgages in the Central United States in areas affected by 2011 flooding. Annual exceedance probabilities (AEPs) were analyzed for 321 streamgages for annual peak streamflow and for 211 streamgages for annual runoff volume. Some of the most exceptional flooding was for the Souris River Basin, where of 11 streamgages considered for AEP analysis of peak streamflow, flood peaks in 2011 exceeded the next largest peak of record by at least double for 6 of the longest-term streamgages (75 to 108 years of peak-flow record). AEPs for these six streamgages were less than 1 percent. AEPs for 2011 runoff volumes were less than 1 percent for all seven Souris River streamgages considered for AEP analysis. Magnitudes of 2011 runoff volumes exceeded previous maxima by double or more for 5 of the 7 streamgages (record lengths 52 to 108 years). For the Red River of the North Basin, AEPs for 2011 runoff volumes were exceptional, with two streamgages having AEPs less than 0.2 percent, five streamgages in the range of 0.2 to 1 percent, and four streamgages in the range of 1 to 2 percent. Magnitudes of 2011 runoff volumes also were exceptional, with all 11 of the aforementioned streamgages eclipsing previous long-term (62 to 110 years) annual maxima by about one-third or more. AEPs for peak streamflows in the upper Mississippi River Basin were not exceptional, with no AEPs less than 1 percent. AEPs for annual runoff volumes indicated less frequent recurrence, with 11 streamgages having AEPs of less than 1 percent. The 2011 runoff volume for streamgage 05331000 (at Saint Paul, Minnesota) exceeded the previous record (112 years of record) by about 24 percent. An especially newsworthy feature was prolonged flooding along the main stem of the Missouri River downstream from Garrison Dam (located upstream from Bismarck, North Dakota) and extending downstream throughout the length of the Missouri River. The 2011 runoff volume for streamgage 06342500 (at Bismarck) exceeded the previous (1975) maximum by about 50 percent, with an associated AEP in the range of 0.2 to 1 percent. In the Ohio River Basin, peak-streamflow AEPs were less than 2 percent for only four streamgages. Runoff-volume AEPs were less than 2 percent for only three streamgages. Along the lower Mississippi River, the largest streamflow peak in 91 years was recorded for streamgage 07289000 (at Vicksburg, Mississippi), with an associated AEP of 0.8 percent. Trends in peak streamflow were analyzed for 98 streamgages, with 67 streamgages having upward trends, 31 with downward trends, and zero with no trend. Trends in annual runoff volume were analyzed for 182 streamgages, with 145 streamgages having upward trends, 36 with downward trends, and 1 with no trend. The trend analyses used descriptive methods that did not include measures of statistical significance. A dichotomous spatial distribution in trends was apparent for both peak streamflow and annual runoff volume, with a small number of streamgages in the northwestern part of the study area having downward trends and most streamgages in the eastern part of the study area having upward trends.

Professional Paper↗
USGS2012

Results of paleoflood investigations for Spring, Rapid, Boxedler, and Elk Creeks, Black Hills, western South Dakota

Flood-frequency analyses for the Black Hills area are especially important because of severe flooding of June 9–10, 1972, that was caused by a large mesoscale convective system and resulted in at least 238 deaths. This paper summarizes results of paleoflood investigations for six study reaches in the central Black Hills. Stratigraphic records and resulting long-term flood chronologies, locally extending more than 2,000 years, were combined with observed and historical flood information to derive flood-frequency estimates. Results indicate that floods as large as and even substantially larger than 1972 have affected most of the study reaches. Results of the paleoflood investigations provide better physically based information on low-probability floods than has been previously available, substantially improving estimates of the magnitude and frequency of large floods in the central Black Hills and reducing associated uncertainties. Collectively, the results provide insights regarding regional flood-generation processes and their spatial controls, enable approaches for extrapolation of results for hazard assessment beyond specific study reaches, and provide a millennial-scale perspective on the 1972 flooding.

South Dakota↗
USGS2011

Flood-frequency analyses from paleoflood investigations for Spring, Rapid, Boxelder, and Elk Creeks, Black Hills, western South Dakota

Flood-frequency analyses for the Black Hills area are important because of severe flooding of June 9-10, 1972, that was caused by a large mesoscale convective system and caused at least 238 deaths. Many 1972 peak flows are high outliers (by factors of 10 or more) in observed records that date to the early 1900s. An efficient means of reducing uncertainties for flood recurrence is to augment gaged records by using paleohydrologic techniques to determine ages and magnitudes of prior large floods (paleofloods). This report summarizes results of paleoflood investigations for Spring Creek, Rapid Creek (two reaches), Boxelder Creek (two subreaches), and Elk Creek. Stratigraphic records and resulting long-term flood chronologies, locally extending more than 2,000 years, were combined with observed and adjusted peak-flow values (gaged records) and historical flood information to derive flood-frequency estimates for the six study reaches. Results indicate that (1) floods as large as and even substantially larger than 1972 have affected most of the study reaches, and (2) incorporation of the paleohydrologic information substantially reduced uncertainties in estimating flood recurrence. Canyons within outcrops of Paleozoic rocks along the eastern flanks of the Black Hills provided excellent environments for (1) deposition and preservation of stratigraphic sequences of late-Holocene flood deposits, primarily in protected slack-water settings flanking the streams; and (2) hydraulic analyses for determination of associated flow magnitudes. The bedrock canyons ensure long-term stability of channel and valley geometry, thereby increasing confidence in hydraulic computations of ancient floods from modern channel geometry. Stratigraphic records of flood sequences, in combination with deposit dating by radiocarbon, optically stimulated luminescence, and cesium-137, provided paleoflood chronologies for 29 individual study sites. Flow magnitudes were estimated from elevations of flood deposits in conjunction with hydraulic calculations based on modern channel and valley geometry. Reach-scale paleoflood chronologies were interpreted for each study reach, which generally entailed correlation of flood evidence among multiple sites, chiefly based on relative position within stratigraphic sequences, unique textural characteristics, or results of age dating and flow estimation. The FLDFRQ3 and PeakfqSA analytical models (assuming log-Pearson Type III frequency distributions) were used for flood-frequency analyses for as many as four scenarios: (1) analysis of gaged records only; (2) gaged records with historical information; (3) all available data including gaged records, historical flows, paleofloods, and perception thresholds; and (4) the same as the third scenario, but ?top fitting? the distribution using only the largest 50 percent of gaged peak flows. The PeakfqSA model is most consistent with procedures adopted by most Federal agencies for flood-frequency analysis and thus was (1) used for comparisons among results for study reaches, and (2) considered by the authors as most appropriate for general applications of estimating low-probability flood recurrence. The detailed paleoflood investigations indicated that in the last 2,000 years all study reaches have had multiple large floods substantially larger than in gaged records. For Spring Creek, stratigraphic records preserved a chronology of at least five paleofloods in approximately (~) 1,000 years approaching or exceeding the 1972 flow of 21,800 cubic feet per second (ft 3 /s). The largest was ~700 years ago with a flow range of 29,300-58,600 ft 3 /s, which reflects the uncertainty regarding flood-magnitude estimates that was incorporated in the flood-frequency analyses. In the lower reach of Rapid Creek (downstream from Pactola Dam), two paleofloods in ~1,000 years exceeded the 1972 flow of 31,200 ft 3 /s. Those occurred ~440 and 1,000 years ago, with flows of 128,000-256,000 and 64,000-128,000 ft 3 /s, respectively. Five smaller paleofloods of 9,500-19,000 ft 3 /s occurred between ~200 and 400 years ago. In the upper reach of Rapid Creek (above Pactola Reservoir), the largest recorded floods are substantially smaller than for lower Rapid Creek and all other study reaches. Paleofloods of ~12,900 and 12,000 ft 3 /s occurred ~1,000 and 1,500 years ago. One additional paleoflood (~800 years ago) was similar in magnitude to the largest gaged flow of 2,460 ft 3 /s Boxelder Creek was treated as having two subreaches because of two tributaries that affect peak flows. During the last ~1,000 years, paleofloods of ~39,000-78,000 ft 3 /s and 40,000-80,000 ft 3 /s in the upstream subreach have exceeded the 1972 peak flow of 30,800 ft 3 /s. One other paleoflood was similar to the second largest gaged flow (16,400 ft 3 /s in 1907). For the downstream subreach, paleofloods of 61,300-123,000 ft 3 /s and 52,500-105,000 ft 3 /s in the last ~1,000 years have substantially exceeded the 1972 flood (50,500 ft 3 /s). Four additional paleofloods had flows between 14,200 and 33,800 ft 3 /s. The 1972 flow on Elk Creek (10,400 ft 3 /s) has been substantially exceeded at least five times in the last 1,900 years. The largest paleoflood (41,500-124,000 ft 3 /s) was ~900 years ago. Three other paleofloods between 37,500 and 120,000 ft 3 /s occurred between 1,100 and 1,800 years ago. A fifth paleoflood of 25,500-76,500 ft 3 /s was ~750 years ago. Considering analyses for all available data (PeakfqSA model) for all six study reaches, the 95-percent confidence intervals about the low-probability quantile estimates (100-, 200-, and 500-year recurrence intervals) were reduced by at least 78 percent relative to those for the gaged records only. In some cases, 95-percent uncertainty intervals were reduced by 99 percent or more. For all study reaches except the two Boxelder Creek subreaches, quantile estimates for these long-term analyses were larger than for the short-term analyses. The 1972 flow for the Spring Creek study reach (21,800 ft 3 /s) corresponds with a recurrence interval of ~400 years. Recurrence intervals are ~500 years for the 1972 flood magnitudes along the lower Rapid Creek reach and the upstream subreach of Boxelder Creek. For the downstream subreach of Boxelder Creek, the large 1972 flood magnitude (50,500 ft 3 /s) exceeds the 500-year quantile estimate by about 35 percent. The recurrence interval of ~100 years for 1972 flooding along the Elk Creek study reach is small relative to other study reaches along the eastern margin of the Black Hills. All of the paleofloods plot within the bounds of a national envelope curve, indicating that the national curve represents exceedingly rare floods for the Black Hills area. Elk Creek, lower Rapid Creek, and the downstream subreach of Boxelder Creek all have paleofloods that plot above a regional envelope curve; in the case of Elk Creek, by a factor of nearly two. The Black Hills paleofloods represent some of the largest known floods, relative to drainage area, for the United States. Many of the other largest known United States floods are in areas with physiographic and climatologic conditions broadly similar to the Black Hills-semiarid and rugged landscapes that intercept and focus heavy precipitation from convective storm systems. The 1972 precipitation and runoff patterns, previous analyses of peak-flow records, and the paleoflood investigations of this study support a hypothesis of distinct differences in flood generation within the central Black Hills study area. The eastern Black Hills are susceptible to intense orographic lifting associated with convective storm systems and also have high relief, thin soils, and narrow and steep canyons-factors favoring generation of exceptionally heavy rain-producing thunderstorms and promoting runoff and rapid concentration of flow into stream channels. In contrast, storm potential is smaller in and near the Limestone Plateau area, and storm runoff is further reduced by substantial infiltration into the limestone, gentle topography, and extensive floodplain storage. Results of the paleoflood investigations are directly applicable only to the specific study reaches and in the case of Rapid Creek, only to pre-regulation conditions. Thus, approaches for broader applications were developed from inferences of overall flood-generation processes, and appropriate domains for application of results were described. Example applications were provided by estimating flood quantiles for selected streamgages, which also allowed direct comparison with results of at-site flood-frequency analyses from a previous study. Several broad issues and uncertainties were examined, including potential biases associated with stratigraphic records that inherently are not always complete, uncertainties regarding statistical approaches, and the unknown applicability of paleoflood records to future watershed conditions. The results of the paleoflood investigations, however, provide much better physically based information on low-probability floods than has been available previously, substantially improving estimates of the magnitude and frequency of large floods in these basins and reducing associated uncertainty.

South Dakota↗
USGS2010

Thunderstorms and flooding of August 17, 2007, with a context provided by a history of other large storm and flood events in the Black Hills area of South Dakota

The Black Hills area of western South Dakota has a history of damaging flash floods that have resulted primarily from exceptionally strong rain-producing thunderstorms. The best known example is the catastrophic storm system of June 9-10, 1972, which caused severe flooding in several major drainages near Rapid City and resulted in 238 deaths. More recently, severe thunderstorms caused flash flooding near Piedmont and Hermosa on August 17, 2007. Obtaining a thorough understanding of peak-flow characteristics for low-probability floods will require a comprehensive long-term approach involving (1) documentation of scientific information for extreme events such as these; (2) long-term collection of systematic peak-flow records; and (3) regional assessments of a wide variety of peak-flow information. To that end, the U.S. Geological Survey cooperated with the South Dakota Department of Transportation and National Weather Service to produce this report, which provides documentation regarding the August 17, 2007, storm and associated flooding and provides a context through examination of other large storm and flood events in the Black Hills area. The area affected by the August 17, 2007, storms and associated flooding generally was within the area affected by the larger storm of June 9-10, 1972. The maximum observed 2007 precipitation totals of between 10.00 and 10.50 inches occurred within about 2-3 hours in a small area about 5 miles west of Hermosa. The maximum documented precipitation amount in 1972 was 15.0 inches, and precipitation totals of 10.0 inches or more were documented for 34 locations within an area of about 76 square miles. A peak flow of less than 1 cubic foot per second occurred upstream from the 2007 storm extent for streamflow-gaging station 06404000 (Battle Creek near Keystone); whereas, the 1972 peak flow of 26,200 cubic feet per second was large, relative to the drainage area of only 58.6 square miles. Farther downstream along Battle Creek, a 2007 flow of 26,000 cubic feet per second was generated entirely within an intervening drainage area of only 44.4 square miles. An especially large flow of 44,100 cubic feet per second was documented for this location in 1972. The 2007 peak flow of 18,600 cubic feet per second for Battle Creek at Hermosa (station 06406000) was only slightly smaller than the 1972 peak flow of 21,400 cubic feet per second. Peak-flow values from 2007 for three sites with small drainage areas (less than 1.0 square mile) plot close to a regional envelope curve, indicating exceptionally large flow values, relative to drainage area. Physiographic factors that affect flooding in the area were examined. The limestone headwater hydrogeologic setting (within and near the Limestone Plateau area on the western flank of the Black Hills) has distinctively suppressed peak-flow characteristics for small recurrence intervals. Uncertainty is large, however, regarding characteristics for large recurrence intervals (low-probability floods) because of a dearth of information regarding the potential for generation of exceptionally strong rain-producing thunderstorms. In contrast, the greatest potential for exceptionally damaging floods is around the flanks of the rest of the Black Hills area because of steep topography and limited potential for attenuation of flood peaks in narrow canyons. Climatological factors that affect area flooding also were examined. Area thunderstorms are largely terrain-driven, especially with respect to their requisite upward motion, which can be initiated by orographic lifting effects, thermally enhanced circulations, and obstacle effects. Several other meteorological processes are influential in the development of especially heavy precipitation for the area, including storm cell training, storm anchoring or regeneration, storm mergers, supercell development, and weak upper-level air flow. A composite of storm total precipitation amounts for 13 recent individual storm events indicates

South Dakota↗
USGS2009

Hydrologic characterization for Spring Creek and hydrologic budget and model scenarios for Sheridan Lake, South Dakota, 1962-2007

The U.S. Geological Survey cooperated with South Dakota Game, Fish and Parks to characterize hydrologic information relevant to management of water resources associated with Sheridan Lake, which is formed by a dam on Spring Creek. This effort consisted primarily of characterization of hydrologic data for a base period of 1962 through 2006, development of a hydrologic budget for Sheridan Lake for this timeframe, and development of an associated model for simulation of storage deficits and drawdown in Sheridan Lake for hypothetical release scenarios from the lake. Historically, the dam has been operated primarily as a 'pass-through' system, in which unregulated outflows pass over the spillway; however, the dam recently was retrofitted with an improved control valve system that would allow controlled releases of about 7 cubic feet per second (ft3/s) or less from a fixed depth of about 60 feet (ft). Development of a hydrologic budget for Sheridan Lake involved compilation, estimation, and characterization of data sets for streamflow, precipitation, and evaporation. The most critical data need was for extrapolation of available short-term streamflow records for Spring Creek to be used as the long-term inflow to Sheridan Lake. Available short-term records for water years (WY) 1991-2004 for a gaging station upstream from Sheridan Lake were extrapolated to WY 1962-2006 on the basis of correlations with streamflow records for a downstream station and for stations located along two adjacent streams. Comparisons of data for the two streamflow-gaging stations along Spring Creek indicated that tributary inflow is approximately proportional to the intervening drainage area, which was used as a means of estimating tributary inflow for the hydrologic budget. Analysis of evaporation data shows that sustained daily rates may exceed maximum monthly rates by a factor of about two. A long-term (1962-2006) hydrologic budget was developed for computation of reservoir outflow from Sheridan Lake for the historical pass-through operating system. Two inflow components (stream inflow and precipitation) and one outflow component (evaporation) were considered. The hydrologic budget uses monthly time steps within a computational year that includes two 6-month periods - May through October, for which evaporation is accounted for, and November through April, when evaporation is considered negligible. Results indicate that monthly evaporation rates can substantially exceed inflow during low-flow periods, and potential exists for outflows to begin approaching zero-flow conditions substantially prior to the onset of zero-inflow conditions, especially when daily inflow and evaporation are considered. Results also indicate that September may be the month for greatest potential benefit for enhancing fish habitat and other ecosystem values in downstream reaches of Spring Creek with managed releases of cool water. Computed monthly outflows from Sheridan Lake for September are less than 1.0 ft3/s for 8 of the 44 years (18 percent) and are less than 2.0 ft3/s for 14 of the 44 years (32 percent). Conversely, none of the computed outflows for May are less than 2.0 ft3/s. A short-term (July through September 2007) data set was used to calculate daily evaporation from Sheridan Lake and to evaluate the applicability of published pan coefficients. Computed values of pan coefficients of approximately 1.0 and 1.1 for two low-flow periods are larger than the mean annual pan coefficient of 0.74 for the area that is reported in the literature; however, the computed values are consistent with pan coefficients reported elsewhere for similar late summer and early fall periods. Thus, these results supported the use of variable monthly pan coefficients for the long-term hydrologic budget. A hydrologic model was developed using the primary components of the hydrologic budget and was used to simulate monthly storage deficits and drawdown for Sheridan Lake using hypothetical

South Dakota↗
USGS2008

Peak-flow frequency estimates based on data through water year 2001 for selected streamflow-gaging stations in South Dakota

Numerous users, including the South Dakota Department of Transportation, have continuing needs for peak-flow information for the design of highway infrastructure and many other purposes. This report documents results from a cooperative study between the South Dakota Department of Transportation and the U.S. Geological Survey to provide an update of peak-flow frequency estimates for South Dakota. Estimates of peak-flow magnitudes for 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals are reported for 272 streamflow-gaging stations, which include most gaging stations in South Dakota with 10 or more years of systematic peak-flow records through water year 2001. Recommended procedures described in Bulletin 17B were used as primary guidelines for developing peak-flow frequency estimates. The computer program PEAKFQ developed by the U.S. Geological Survey was used to run the frequency analyses. Flood frequencies for all stations were initially analyzed by using standard Bulletin 17B default procedures for fitting the log-Pearson III distribution. The resulting preliminary frequency curves were then plotted on a log-probability scale, and fits of the curves with systematic data were evaluated. In many cases, results of the default Bulletin 17B analyses were determined to be satisfactory. In other cases, however, the results could be improved by using various alternative procedures for frequency analysis. Alternative procedures for some stations included adjustments to skew coefficients or use of user-defined low-outlier criteria. Peak-flow records for many gaging stations are strongly influenced by low- or zero-flow values. This situation often results in a frequency curve that plots substantially above the systematic record data points at the upper end of the frequency curve. Adjustments to low-outlier criteria reduced the influence of very small peak flows and generally focused the analyses on the upper parts of the frequency curves (10- to 500-year recurrence intervals). The most common alternative procedures involved several different methods to extend systematic records, which was done primarily to address biases resulting from nonrepresentative climatic conditions during several specific periods of record and to reduce inconsistencies among multiple gaging stations along common stream channels with different periods of record. In some cases, records for proximal stations could be combined directly. In other cases, the two-station comparison procedure recommended in Bulletin 17B was used to adjust the mean and standard deviation of the logs of the systematic data for a target station on the basis of correlation with concurrent records from a nearby long-term index station. In some other cases, a 'mixed-station procedure' was used to adjust the log-distributional parameters for a target station, on the basis of correlation with one or more index stations, for the purpose of fitting the log-Pearson III distribution. Historical adjustment procedures were applied to peak-flow frequency analyses for 17 South Dakota gaging stations. A historical adjustment period extending back to 1881 (121 years) was used for 12 gaging stations in the James and Big Sioux River Basins, and various other adjustment periods were used for additional stations. Large peak flows that occurred in 1969 and 1997 accounted for 13 of the 17 historical adjustments. Other years for which historical peak flows were used include 1957, 1962, 1992, and 2001. A regional mixed-population analysis was developed to address complications associated with many high outliers for the Black Hills region. This analysis included definition of two populations of flood events. The population of flood events that composes the main body of peak flows for a given station is considered the 'ordinary-peaks population,' and the population of unusually large peak flows that plot substantially above the main body of peak flows on log-probability scale is co

South Dakota↗
USGS2004

Hydrologic Effects of the 1988 Galena Fire, Black Hills Area, South Dakota

The Galena Fire burned about 16,788 acres of primarily ponderosa pine forest during July 5-8, 1988, in the Black Hills area of South Dakota. The fire burned primarily within the Grace Coolidge Creek drainage basin and almost entirely within the boundaries of Custer State Park. A U.S. Geological Survey gaging station with streamflow records dating back to 1977 was located along Grace Coolidge Creek within the burned area. About one-half of the gaging station's 26.8-square-mile drainage area was burned. The drainage basin for Bear Gulch, which is tributary to Grace Coolidge Creek, was burned particularly severely, with complete deforestation occurring in nearly the entirety of the area upstream from a gaging station that was installed in 1989. A study to evaluate effects of the Galena Fire on streamflow, geomorphology, and water quality was initiated in 1988. The geomorphologic and water-quality components of the study were completed by 1990 and are summarized in this report. A data-collection network consisting of streamflow- and precipitation-gaging stations was operated through water year 1998 for evaluation of effects on streamflow characteristics, including both annual-yield and peak-flow characteristics, which are the main focus of this report. Moderately burned areas did not experience a substantial increase in the rate of surface erosion; however, severely burned areas underwent surficial erosion nearly twice that of the unburned areas. The sediment production rate of Bear Gulch estimated 8 to 14 months after the fire was 870 ft3/acre (44 tons/acre). Substantial degradation of stream channels within the severely burned headwater areas of Bear Gulch was documented. Farther downstream, channel aggradation resulted from deposition of sediments transported from the headwater areas. The most notable water-quality effect was on concentrations of suspended sediment, which were orders of magnitude higher for Bear Gulch than for the unburned control area. Effects on several other water-quality constituents, such as organic carbon and nitrogen and phosphorus nutrient constituents, probably were influenced by the large concentrations of suspended matter that were documented in initial post-fire, storm-flow events. The first post-fire stormflow produced the highest measured concentrations of specific conductance, nitrogen, phosphorus, organic carbon, calcium, magnesium, potassium, manganese, and sulfate in the burned areas. For most constituents sampled, differences in concentrations between burned and unburned areas were no longer discernible within about 1 year following the Galena Fire. The effects of the Galena Fire on annual-yield characteristics of Grace Coolidge Creek were evaluated primarily from comparisons with long-term streamflow records for Battle Creek, which is hydrogeologically similar and is located immediately to the north. Annual yield for Grace Coolidge Creek increased by about 20 percent as a result of the fire. This estimate was based on relations between annual yield for Grace Coolidge Creek and Battle Creek for pre- and post-burn periods. Many of the post-burn data points are well beyond the range of the pre-burn data, which is a source of uncertainty for this estimate. Substantial increases in peak-flow characteristics for severely burned drainages were visually apparent from numerous post-fire field observations. Various analyses of streamflow data indicated substantial increases in peak-flow response for burned drainage areas; however, quantification of effects was particularly difficult because peak-flow response diminished quickly and returned to a generally pre-burn condition by about 1991. Field observations of vegetation and analysis of remotely sensed data indicated that establishment of grasses and forbs occurred within a similar timeframe. Comparison of pre-fire peak flows to post-1991 peak flows indicates that these grasses and forbs were equally effective in suppressing peak flows

Water-Resources Investigations Report↗
USGS2003

Ground-Water Resources in the Black Hills Area, South Dakota

The availability of ground-water resources in the Black Hills area is influenced by many factors including location, local recharge and ground-water flow conditions, and structural features. Thus, the availability of ground water can be extremely variable throughout the Black Hills area, and even when water is available, it may not be suitable for various uses depending on the water quality. The major bedrock aquifers in the Black Hills area are the Deadwood, Madison, Minnelusa, Minnekahta, and Inyan Kara aquifers. Minor bedrock aquifers occur in other hydrogeologic units, including confining units, due to fracturing and interbedded permeable layers. Various information and maps are presented in this report that describe availability and quality of ground-water resources in the Black Hills area. However, there is no guarantee of obtaining usable water at any location due to the extreme potential variability in conditions that can affect the availability and quality of ground water in the area. Maps presented in this report include the distribution of hydrogeologic units; depth to the top of the five formations that contain major aquifers; thickness of the five formations that contain major aquifers; potentiometric maps for the five major aquifers; saturated thickness of the Madison and Minnelusa aquifers; water temperature in the Madison aquifer; specific conductance in the Madison, Minnelusa, and Inyan Kara aquifers; hardness in the Inyan Kara aquifer; sulfate concentrations in the Minnelusa aquifer; and radon concentrations in the Deadwood aquifer. Water quality of the major aquifers generally is very good in and near outcrop areas but deteriorates progressively with distance from the outcrops. In the Minnelusa aquifer, an abrupt increase in concentrations of dissolved sulfate occurs downgradient from outcrop areas, where a zone of active anhydrite dissolution occurs. Most limitations for the use of ground water are related to aesthetic qualities associated with hardness and high concentrations of chloride, sulfate, sodium, manganese, and iron. Very few health-related limitations exist for ground water; most limitations are for radionuclides, such as radon and uranium. In addition, high concentrations of arsenic have been measured in a few samples from the Minnelusa aquifer.

Water-Resources Investigations Report↗
USGS2002

Atlas of water resources in the Black Hills area, South Dakota

The Black Hills area is an important resource center that provides an economic base for western South Dakota through tourism, agriculture, the timber industry, and mineral resources. In addition, water originating from the area is used for municipal, industrial, agricultural, and recreational purposes throughout much of western South Dakota. The Black Hills area also is an important recharge area for aquifers in the northern Great Plains. Population growth, resource development, and periodic droughts have the potential to affect the quantity, quality, and availability of water within the Black Hills area. Growth has resulted in competing interests for available water supplies. The Black Hills Hydrology Study was initiated in 1990 to address these concerns. This long-term study is a cooperative effort between the U.S. Geological Survey (USGS), the South Dakota Department of Environment and Natural Resources, and the West Dakota Water Development District, which represents various local and county cooperators.

Hydrologic Atlas↗
USGS2001

Hydrologic budgets for the Madison and Minnelusa aquifers, Black Hills of South Dakota and Wyoming, water years 1987-96

The Madison and Minnelusa aquifers are two of the most important aquifers in the Black Hills area of South Dakota and Wyoming. Quantification and evaluation of various hydrologic budget components are important for managing and understanding these aquifers. Hydrologic budgets are developed for two scenarios, including an overall budget for the entire study area and more detailed budgets for subareas. Budgets generally are combined for the Madison and Minnelusa aquifers because most budget components cannot be quantified individually for the aquifers. An average hydrologic budget for the entire study area is computed for water years 1987-96, for which change in storage is approximately equal to zero. Annual estimates of budget components are included in detailed budgets for nine subareas, which consider periods of decreasing storage (1987-92) and increasing storage (1993-96). Inflow components include recharge, leakage from adjacent aquifers, and ground-water inflows across the study area boundary. Outflows include springflow (headwater and artesian), well withdrawals, leakage to adjacent aquifers, and ground-water outflow across the study area boundary. Leakage, ground-water inflows, and ground-water outflows are difficult to quantify and cannot be distinguished from one another. Thus, net ground-water flow, which includes these components, is calculated as a residual, using estimates for the other budget components. For the overall budget for water years 1987-96, net ground-water outflow from the study area is computed as 100 ft3/s (cubic feet per second). Estimates of average combined budget components for the Madison and Minnelusa aquifers are: 395 ft3/s for recharge, 78 ft3/s for headwater springflow, 189 ft3/s for artesian springflow, and 28 ft3/s for well withdrawals. Hydrologic budgets also are quantified for nine subareas for periods of decreasing storage (1987-92) and increasing storage (1993-96), with changes in storage assumed equal but opposite. Common subareas are identified for the Madison and Minnelusa aquifers, and previous components from the overall budget generally are distributed over the subareas. Estimates of net ground-water flow for the two aquifers are computed, with net ground-water outflow exceeding inflow for most subareas. Outflows range from 5.9 ft3/s in the area east of Rapid City to 48.6 ft3/s along the southwestern flanks of the Black Hills. Net groundwater inflow exceeds outflow for two subareas where the discharge of large artesian springs exceeds estimated recharge within the subareas. More detailed subarea budgets also are developed, which include estimates of flow components for the individual aquifers at specific flow zones. The net outflows and inflows from the preliminary subarea budgets are used to estimate transmissivity of flow across specific flow zones based on Darcy?s Law. For estimation purposes, it is assumed that transmissivities of the Madison and Minnelusa aquifers are equal in any particular flow zone. The resulting transmissivity estimates range from 90 ft2/d to about 7,400 ft2/d, which is similar to values reported by previous investigators. The highest transmissivity estimates are for areas in the northern and southwestern parts of the study area, and the lowest transmissivity estimates are along the eastern study area boundary. Evaluation of subarea budgets provides confidence in budget components developed for the overall budget, especially regarding precipitation recharge, which is particularly difficult to estimate. Recharge estimates are consistently compatible with other budget components, including artesian springflow, which is a dominant component in many subareas. Calculated storage changes for subareas also are consistent with other budget components, specifically artesian springflow and net ground-water flow, and also are consistent with water-level fluctuations for observation wells. Ground-water budgets and flowpaths are especially complex i

Water-Resources Investigations Report↗
USGS2001

Geochemistry of the Madison and Minnelusa aquifers in the Black Hills area, South Dakota

The Madison and Minnelusa aquifers are two of the most important aquifers in the Black Hills area because of utilization for water supplies and important influences on surface-water resources resulting from large springs and streamflow- loss zones. Examination of geochemical information provides a better understanding of the complex flow systems within these aquifers and interactions between the aquifers. Major-ion chemistry in both aquifers is dominated by calcium and bicarbonate near outcrop areas, with basinward evolution towards various other water types. The most notable differences in major-ion chemistry between the Madison and Minnelusa aquifers are in concentrations of sulfate within the Minnelusa aquifer. Sulfate concentrations increase dramatically near a transition zone where dissolution of anhydrite is actively occurring. Water chemistry for the Madison and Minnelusa aquifers is controlled by reactions among calcite, dolomite, and anhydrite. Saturation indices for gypsum, calcite, and dolomite for most samples in both the Madison and Minnelusa aquifers are indicative of the occurrence of dedolomitization. Because water in the Madison aquifer remains undersaturated with respect to gypsum, even at the highest sulfate concentrations, upward leakage into the overlying Minnelusa aquifer has potential to drive increased dissolution of anhydrite in the Minnelusa Formation. Isotopic information is used to evaluate ground-water flowpaths, ages, and mixing conditions for the Madison and Minnelusa aquifers. Distinctive patterns exist in the distribution of stable isotopes of oxygen and hydrogen in precipitation for the Black Hills area, with isotopically lighter precipitation generally occurring at higher elevations and latitudes. Distributions of 18O in ground water are consistent with spatial patterns in recharge areas, with isotopically lighter 18O values in the Madison aquifer resulting from generally higher elevation recharge sources, relative to the Minnelusa aquifer. Three conceptual models, which are simplifications of lumped-parameter models, are considered for evaluation of mixing conditions and general ground-water ages. For a simple slug-flow model, which assumes no mixing, measured tritium concentrations in ground water can be related through a first-order decay equation to estimated concentrations at the time of recharge. Two simplified mixing models that assume equal proportions of annual recharge over a range of years also are considered. An ?immediate-arrival? model is used to conceptually represent conditions in outcrop areas and a ?time-delay? model is used for locations removed from outcrops, where delay times for earliest arrival of ground water generally would be expected. Because of limitations associated with estimating tritium input and gross simplifying assumptions of equal annual recharge and thorough mixing conditions, the conceptual models are used only for general evaluation of mixing conditions and approximation of age ranges. Headwater springs, which are located in or near outcrop areas, have the highest tritium concentrations, which is consistent with the immediate-arrival mixing model. Tritium concentrations for many wells are very low, or nondetectable, indicating general applicability of the timedelay conceptual model for locations beyond outcrop areas, where artesian conditions generally occur. Concentrations for artesian springs generally are higher than for wells, which indicates generally shorter delay times resulting from preferential flowpaths that typically are associated with artesian springs. In the Rapid City area, a distinct division of isotopic values for the Madison aquifer corresponds with distinguishing 18O signatures for nearby streams, where large streamflow recharge occurs. Previous dye testing in this area documented rapid ground-water flow (timeframe of weeks) from a streamflow loss zone to sites located several miles away. These results are used to ill

Water-Resources Investigations Report↗