Geology topics

Research about South Dakota

Source-linked reports with geographic coverage including South Dakota.

At least 163 records · Page 9Linked to original sources

USGS2011

Movements and survival of black-footed ferrets associated with an experimental translocation in South Dakota

Black-footed ferrets ( Mustela nigripes ) apparently were extirpated from all native habitats by 1987, and their repatriation requires a combination of captive breeding, reintroductions, and translocations among sites. Improvements in survival rates of released ferrets have resulted from experience in quasi-natural environments during their rearing. Reestablishment of a self-sustaining wild population by 1999 provided the 1st opportunity to initiate new populations by translocating wild-born individuals. Using radiotelemetry, we compared behaviors and survival of 18 translocated wild-born ferrets and 18 pen-experienced captive-born ferrets after their release into a prairie dog colony not occupied previously by ferrets. Translocated wild-born ferrets moved significantly less and had significantly higher short-term survival rates than their captive-born counterparts. Using mark–recapture methods, we also assessed potential impacts to the established donor population of removing 37% of its estimated annual production of kits. Annual survival rates for 30 ferret kits remaining at the donor subcomplex were higher than rates for 54 ferret kits at the control subcomplex (unmanipulated) for males (+82%) and females (+32%). Minimum survival of translocated kits did not differ significantly from survival of those at the control subcomplex. Direct translocation of young, wild-born ferrets from site to site appears to be an efficient method to establish new populations.

South Dakota↗

USGS2011

Interrelationships between fish tissue mercury concentrations and water quality for South Dakota natural lakes and impoundments

The purpose of this study was to determine whether water quality parameters commonly associated with primary productivity may be used to predict the susceptibility of a specific water body to exceed proposed fish consumption advisory limitation of 0.3 mg kg −1 . South Dakota currently has nine lakes and impoundments that exceed fish tissue mercury advisory limits of 1.0 mg kg −1 total mercury, far exceeding US Environmental Protection Agency and Food and Drug Administration 0.3 mg kg −1 consumption criteria. Previous studies suggest that increased aquatic productivity may mitigate the effects of biological production and subsequent uptake of methyl mercury through bio-dilution; however, it is uncertain whether these trends may exist within highly alkaline and highly productive aquatic conditions common to South Dakota lakes and impoundments. Water quality parameters and fish tissue mercury data for northern pike and walleye were collected and assessed using existing South Dakota Department of Environment and Natural Resources and Game Fish and Parks data. The data was initially screened using both parametric linear regression and non-parametric Mann–Whitney rank sum comparisons and further assessed using binary logistic regression and stepwise logistic regression methodology. Three separate phosphorus measurements (total, total dissolved, and Trophic State Index) and pH were determined to significantly correlate with increased mercury concentrations for the northern pike-in-impoundments model. However, phosphorus surprisingly was not a strong predictor for the remaining scenarios modeled. For the northern pike-in-natural lakes models, alkalinity was the most significant water quality parameter predicting increased mercury concentrations. Mercury concentrations for the walleye-in-natural lakes models were further influenced by pH and alkalinity. The water quality and fish tissue mercury interrelationships determined within this study suggest aquatic productivity, and consequential eutrophication processes appear to be reasonable indicators of fish tissue mercury susceptibility for aquatic conditions common to South Dakota and highlight the continuing need to minimize eutrophication through effective watershed management strategies.

South Dakota↗

USGS2011

Evaluation of a black-footed ferret resource utilization function model

Resource utilization function (RUF) models permit evaluation of potential habitat for endangered species; ideally such models should be evaluated before use in management decision‐making. We evaluated the predictive capabilities of a previously developed black‐footed ferret ( Mustela nigripes ) RUF. Using the population‐level RUF, generated from ferret observations at an adjacent yet distinct colony, we predicted the distribution of ferrets within a black‐tailed prairie dog ( Cynomys ludovicianus ) colony in the Conata Basin, South Dakota, USA. We evaluated model performance, using data collected during post‐breeding spotlight surveys (2007–2008) by assessing model agreement via weighted compositional analysis and count‐metrics. Compositional analysis of home range use and colony‐level availability, and core area use and home range availability, demonstrated ferret selection of the predicted Very high and High occurrence categories in 2007 and 2008. Simple count‐metrics corroborated these findings and suggested selection of the Very high category in 2007 and the Very high and High categories in 2008. Collectively, these results suggested that the RUF was useful in predicting occurrence and intensity of space use of ferrets at our study site, the 2 objectives of the RUF. Application of this validated RUF would increase the resolution of habitat evaluations, permitting prediction of the distribution of ferrets within distinct colonies. Additional model evaluation at other sites, on other black‐tailed prairie dog colonies of varying resource configuration and size, would increase understanding of influences upon model performance and the general utility of the RUF.

South Dakota↗

USGS2011

Multivariate analyses with end-member mixing to characterize groundwater flow: Wind Cave and associated aquifers

Principal component analysis (PCA) applied to hydrochemical data has been used with end-member mixing to characterize groundwater flow to a limited extent, but aspects of this approach are unresolved. Previous similar approaches typically have assumed that the extreme-value samples identified by PCA represent end members. The method presented herein is different from previous work in that (1) end members were not assumed to have been sampled but rather were estimated and constrained by prior knowledge; (2) end-member mixing was quantified in relation to hydrogeologic domains, which focuses model results on major hydrologic processes; (3) a method to select an appropriate number of end members using a series of cluster analyses is presented; and (4) conservative tracers were weighted preferentially in model calibration, which distributed model errors of optimized values, or residuals, more appropriately than would otherwise be the case. The latter item also provides an estimate of the relative influence of geochemical evolution along flow paths in comparison to mixing. This method was applied to groundwater in Wind Cave and the associated karst aquifer in the Black Hills of South Dakota, USA. The end-member mixing model was used to test a hypothesis that five different end-member waters are mixed in the groundwater system comprising five hydrogeologic domains. The model estimated that Wind Cave received most of its groundwater inflow from local surface recharge with an additional 33% from an upgradient aquifer. Artesian springs in the vicinity of Wind Cave primarily received water from regional groundwater flow.

South Dakota↗

USGS2011

Flood hydrology and dam-breach hydraulic analyses of four reservoirs in the Black Hills, South Dakota

Extensive information about the construction of dams or potential downstream hazards in the event of a dam breach is not available for many small reservoirs within the Black Hills National Forest. In 2009, the U.S. Forest Service identified the need for reconnaissance-level dam-breach assessments for four of these reservoirs within the Black Hills National Forest (Iron Creek, Horsethief, Lakota, and Mitchell Lakes) with the potential to flood downstream structures. Flood hydrology and dam-breach hydraulic analyses for the four selected reservoirs were conducted by the U.S. Geological Survey in cooperation with the U.S. Forest service to estimate the areal extent of downstream inundation. Three high-flow breach scenarios were considered for cases when the dam is in place (overtopped) and when a dam break (failure) occurs: the 100-year recurrence 24-hour precipitation, 500-year recurrence peak flow, and the probable maximum precipitation. Inundation maps were developed that show the estimated extent of downstream floodwaters from simulated scenarios. Simulation results were used to determine the hazard classification of a dam break (high, significant, or low), based primarily on the potential for loss of life or property damage resulting from downstream inundation because of the flood surge. The inflow design floods resulting from the two simulated storm events (100-year 24-hour and probable maximum precipitation) were determined using the U.S. Army Corps of Engineers Hydrologic Engineering Center Hydrologic Modeling System (HEC-HMS). The inflow design flood for the 500-year recurrence peak flow was determined by using regional regression equations developed for streamflow-gaging stations with similar watershed characteristics. The step-backwater hydraulic analysis model, Hydrologic Engineering Center's River Analysis System (HEC-RAS), was used to determine water-surface profiles of in-place and dam-break scenarios for the three inflow design floods that were simulated. Inundation maps for in-place and dam-break scenarios were developed for the area downstream from the dam to the mouth of each stream. Dam-break scenarios for three of the four reservoirs assessed in this study were rated as low hazards owing to absence of permanent structures downstream from the dams. Iron Creek Lake's downstream channel to its mouth does not include any permanent structures within the inundation flood plains. For the two reservoirs with the largest watershed areas, Lakota and Mitchell Lake, the additional floodwater surge resulting from a dam break would be minor relative to the magnitude of the large flood streamflow into the reservoirs, based on the similar areal extent of inundation for the in-place and dam-break scenarios as indicated by the developed maps. A dam-break scenario at Horsethief Lake is rated as a significant hazard because of potential lives-in-jeopardy in downstream dwellings and appreciable economic loss.

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↗

USGS2011

Nesting success and resource selection of greater sage grouse in South Dakota: Chapter 8

Declines of Greater Sage-Grouse ( Centrocercus urophasianus ) in South Dakota are a concern because further population declines may lead to isolation from populations in Wyoming and Montana. Furthermore, little information exists about reproductive ecology and resource selection of sage grouse on the eastern edge of their distribution. We investigated Greater Sage-Grouse nesting success and resource selection in South Dakota during 2006-2007. Radiomarked females were tracked to estimate nesting rates, nest success, and habitat resources selected for nesting. Nest initiation was 98.0%, with a maximum likelihood estimate of nest success of 45.6 ± 5.3%. Females selected nest sites that had greater sagebrush canopy cover and visual obstruction of the nest bowl compared to random sites. Nest survival models indicated that taller grass surrounding nests increased nest survival. Tall grass may supplement the low sagebrush cover in this area in providing suitable nest sites for Greater Sage-Grouse. Land managers on the eastern edge of Greater Sage-Grouse range could focus on increasing sagebrush density while maintaining tall grass by developing range management practices that accomplish this goal. To achieve nest survival rates similar to other populations, predictions from our models suggest 26 cm grass height would result in approximately 50% nest survival. Optimal conditions could be accomplished by adjusting livestock grazing systems and stocking rates.

South Dakota↗

USGS2011

Habitat selection by female swift foxes (Vulpes velox) during the pup-rearing season

The swift fox (Vulpes velox) was historically distributed in western South Dakota including the region surrounding Badlands National Park (BNP). The species declined during the mid-1800s, largely due to habitat loss and poisoning targeted at wolves (Canis lupis) and coyotes (C. latrans). Only a small population of swift foxes near Ardmore, South Dakota persisted. In 2003, an introduction program was initiated at BNP with swift foxes translocated from Colorado and Wyoming. We report on habitat use by female swift foxes during the pup-rearing season (May–July) in 2009. Analyses of location data from 13 radiomarked female foxes indicated disproportional use ( P < 0.001) of some habitats relative to their availability within swift fox home ranges. Swift foxes used grassland ( &#374 = 1.01), sparse vegetation ( &#374 = 1.43) and prairie dog towns ( &#374 = 1.18) in proportion to their availability, whereas they were less likely to use woodland ( &#374 = 0.00), shrubland ( &#374 = 0.14), pasture/agricultural-land ( &#374 = 0.25) and development ( &#374 = 0.16) relative to availability. Swift foxes typically are located in habitats that provide greater visibility, such as shortgrass prairie and areas with sparse vegetation; which allow detection of approaching coyotes (e.g., primary predator of swift foxes).

South Dakota↗

USGS2010

Transport and accumulation of radionuclides and stable elements in a Missouri River Reservoir

Several long sediment cores from the Cheyenne River Embayment of Lake Oahe, a 250-km-long Missouri River reservoir in South Dakota, have been analyzed for radionuclides and stable elements. The combination of fine-scale sampling and rapid sedimentation produces radionuclide distributions that can be used to estimate the detailed chronology of particle transport processes in the Oahe reservoir system. A self-consistent and quantitative treatment of the 137 Cs data suggests processes to which characteristic times may be associated. Times that characterize system-wide processes include (1) an integration time of several years reflecting retention of the sediment-bound tracer in regions within or external to the reservoir, (2) a relaxation time of approximately 15 years reflecting a decreasing rate of sediment accumulation ascribed to shoreline stabilization, (3) a time of a few months characterizing the breadth of riverine signatures in cores due to integration effects in the Cheyenne River system and deltaic deposits, and (4) times of a few years associated with propagation of riverine load signatures along the embayment. The distribution of total sedimentary arsenic confirms the validity of the variable sedimentation model. In 1977, a tailings retention facility was built at the Homestake Mine site, and the unrestricted input of As ceased. As a result of this remedial action, the concentration of sedimentary As decreased dramatically. In the upper section of the core, above the depth represented by the year 1976, the concentration of As decreases tenfold. In this same core the distribution of lithologically discriminating chemical elements, calcium and vanadium, relate to major flow events in the Cheyenne River basin. Because there is minimal diagenesis of chemical constituents in these rapidly accumulating sediments, stable element signatures, in addition to radiotracers, may be used to reconstruct hydrologic events in drainage basins that contribute sediment to lakes and reservoirs.

South Dakota↗

USGS2010

Influence of drought conditions on brown trout biomass and size structure in the Black Hills, South Dakota

We evaluated the influence of drought conditions on the biomass of brown trout Salmo trutta in Spearfish Creek, upper Rapid Creek, and lower Rapid Creek in the Black Hills of western South Dakota. Stream discharge, mean summer water temperature, the biomass of juvenile and adult brown trout, and brown trout size structure were compared between two time periods: early (2000–2002) and late drought (2005–2007). Mean summer water temperatures were similar between the early- and late-drought periods in Spearfish Creek (12.4°C versus 11.5°C), lower Rapid Creek (19.2°C versus 19.3°C), and upper Rapid Creek (9.8°C in both periods). In contrast, mean annual discharge differed significantly between the two time periods in Spearfish Creek (1.95 versus 1.50 m 3 /s), lower Rapid Creek (2.01 versus 0.94 m 3 /s), and upper Rapid Creek (1.41 versus 0.84 m 3 /s). The mean biomass of adult brown trout in all three stream sections was significantly higher in the early-drought than in the late-drought period (238 versus 69 kg/ha in Spearfish Creek, 272 versus 91 kg/ha in lower Rapid Creek, and 159 versus 32 kg/ha in upper Rapid Creek). The biomass of juvenile brown trout was similar (43 versus 23 kg/ha) in Spearfish Creek in the two periods, declined from 136 to 45 kg/ha in lower Rapid Creek, and increased from 14 to 73 kg/ha in upper Rapid Creek. Size structure did not differ between the early- and late-drought periods in lower Rapid and Spearfish creeks, but it did in upper Rapid Creek. In addition to drought conditions, factors such as angler harvest, fish movements, and the nuisance algal species Didymosphenia geminata are discussed as possible contributors to the observed changes in brown trout biomass and size structure in Black Hills streams.

South Dakota↗

USGS2010

Morning ambush attacks by black-footed ferrets on emerging prairie dogs

Black-footed ferrets ( Mustela nigripes ) often hunt at night, attacking normally diurnal prairie dogs ( Cynomys spp.) in underground burrow systems. While monitoring black-footed ferrets in South Dakota during morning daylight hours, we observed an adult female ferret ambush a black-tailed prairie dog ( C. ludovicianus ) emerging from a burrow. On a neighboring colony, we observed a second adult female ferret engaging in similar ambush behaviors on 12 occasions, although prey was not visible. We retrospectively assessed radio-telemetry data on white-tailed prairie dogs ( C. leucurus ) and a male and a female ferret to evaluate ferret activity in relation to timing of prairie dog emergence. Activity of radio-collared ferrets was high during the hourly period when prairie dogs first emerged and the following 2 hr, relative to later daylight hours. Such behavior is consistent with behaviors observed in South Dakota. Nighttime movements by ferrets might involve hunting but also reconnaissance of prey preparatory to morning ambush attacks.

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↗

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↗

USGS2010

Repeated use of an abandoned vehicle by nesting Turkey vultures ( Cathartes aura )

Turkey Vultures ( Cathartes aura ) lay their eggs on an existing substrate in the dark recesses of a variety of natural sites ( Kirk and Mossman 1998 ). Although an important requirement of Turkey Vulture nest-site selection is isolation from human disturbances ( Kirk and Mossman 1998 ), their nests have been reported in abandoned buildings since at least the early 1800s ( Nuttall 1832 ). Depopulation of rural areas in North America in recent decades has resulted in many abandoned buildings within the Turkey Vulture's breeding range ( Peck 2003 ). Increased use of abandoned buildings by nesting Turkey Vultures has been implicated in the species' recent northward range expansion ( Peck 2003 , Nelson et al. 2005 , Houston et al. 2007 ). Although abandoned or inoperative vehicles also are widespread in rural areas, we found no published literature documenting Turkey Vultures' use of these potential nest sites. Herein, we summarize the first documented incidence of a Turkey Vulture nesting in an abandoned vehicle.

South Dakota↗

USGS2009

Seasonal movements and home ranges of white-tailed deer in north-central South Dakota

Knowledge of movement patterns of white-tailed deer ( Odocoileus virginianus (Zimmermann, 1780)) inhabiting landscapes intensively modified by agricultural systems is important to the present and future understanding of deer ecology. Little information exists regarding daily and seasonal movements of white-tailed deer in north-central South Dakota. Therefore, our goal was to determine movement patterns and home-range use of female white-tailed deer in north-central South Dakota. From January 2005 to January 2007, 29 adult (>18 months) and 13 yearling (8–18 months) white-tailed deer were monitored for movement using radiotelemetry. We collected 2822 locations, calculated 76 home ranges, and documented 50 seasonal movements. Mean migration distance between summer and winter home ranges was 19.4 km (SE = 2.0 km). Mean 95% home-range size was 10.2 km 2 (SE = 1.2 km 2 , n = 27) during winter and 9.2 km 2 (SE = 1.0 km 2 , n = 49) during summer. Ambient temperature appeared to be a primary cause of seasonal migration. Additionally, movements exhibited by white-tailed deer in north-central South Dakota were influenced by a highly fragmented landscape dominated by row crops and pasture or grassland.

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↗

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↗

USGS2009

Role of invasive Melilotus officinalis in two native plant communities

This study examines the impact of the exotic nitrogen-fixing legume Melilotus officinalis (L.) Lam. on native and exotic species cover in two Great Plains ecosystems in Badlands National Park, South Dakota. Melilotus is still widely planted and its effects on native ecosystems are not well studied. Melilotus could have direct effects on native plants, such as through competition or facilitation. Alternatively, Melilotus may have indirect effects on natives, e.g., by favoring exotic species which in turn have a negative effect on native species. This study examined these interactions across a 4-year period in two contrasting vegetation types: Badlands sparse vegetation and western wheatgrass ( Pascopyrum smithii ) mixed-grass prairie. Structural equation models were used to analyze the pathways through which Melilotus , native species, and other exotic species interact over a series of 2-year time steps. Melilotus can affect native and exotic species both in the current year and in the years after its death (a lag effect). A lag effect is possible because the death of a Melilotus plant can leave an open, potentially nitrogen-enriched site on the landscape. The results showed that the relationship between Melilotus and native and exotic species varied depending on the habitat and the year. In Badlands sparse vegetation, there was a consistent, strong, and positive relationship between Melilotus cover and native and exotic species cover suggesting that Melilotus is acting as a nurse plant and facilitating the growth of other species. In contrast, in western wheatgrass prairie, Melilotus was acting as a weak competitor and had no consistent effect on other species. In both habitats, there was little evidence for a direct lag effect of Melilotus on other species. Together, these results suggest both facilitative and competitive roles for Melilotus , depending on the vegetation type it invades.

South Dakota↗