Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Hydrogeology Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9Linked to original sources
USGS2019

Hydrogeology of an alpine talus aquifer: Cordillera Blanca, Peru

The dramatic loss of glacial mass in low latitudes is causing shifts in downstream water availability and use during the driest months of the year. The world’s largest concentration of tropical glaciers lies in the Cordillera Blanca range of Peru, where glacial runoff is declining and regional stresses are emerging over water resources. Throughout the Cordillera Blanca, groundwater inputs from alpine meadow–talus complexes, locally known as pampas, supply proglacial streams with up to 80% of their flow during the region’s dry season. Structural knowledge of the pampa aquifers is needed to estimate their drainable groundwater storage capacity and residence time, to elucidate the role and importance of alpine groundwater storage in the regional water budget of the Cordillera Blanca. To understand the structure of these proglacial aquifers, multiple near-surface geophysical methods were implemented in a proglacial valley near dense networks of spring-fed tributaries. Geophysical results and borehole logs suggest groundwater is stored in a confined aquifer composed of buried talus deposits overlain by lacustrine clay, while deeper portions of the unit, 10–15 m in depth, are relatively clay-free and more hydraulically conductive. Based on these findings and assumptions of aquifer porosity, the pampas of the Callejon de Huaylas may store from 0.006 to 0.02 km 3 of groundwater. Furthermore, these findings suggest that the talus aquifers of the Cordillera Blanca were formed in proglacial lakes, followed by infilling with fine lacustrine sediments that confine lower units and allow for groundwater discharge to springs via macropores and preferential flow.

Hydrogeology Journal↗
USGS2020

Conjoint use of hydraulic head and groundwater age data to detect hydrogeologic barriers

Hydraulic head and groundwater age data are effective in building understanding of groundwater systems. Yet their joint role in detecting and characterising low-permeability geological structures, i.e. hydrogeologic barriers such as faults and dykes, has not been widely studied. Here, numerical flow and transport models, using MODFLOW-NWT and MT3D-USGS, were developed with different hydrogeologic barrier configurations in a hypothetical aquifer. Computed hydraulic head and groundwater age distributions were compared to those without a barrier. The conjoint use of these datasets helps in detecting vertically-oriented barriers. Two forms of recharge were compared: (1) applied across the entire aquifer surface (uniform), and (2) applied to the upstream part of the aquifer (upgradient). The hydraulic head distribution is significantly impacted by a barrier that penetrates the aquifer’s full vertical thickness. This barrier also perturbs the groundwater age distribution when upgradient recharge prevails; however, with uniform recharge, groundwater age is not successful in detecting the barrier. When a barrier is buried, such as by younger sediment, hydraulic head data also do not clearly identify the barrier. Groundwater age data could, on the other hand, prove to be useful if sampled at depth-specific intervals. These results are important for the detection and characterisation of hydrogeologic barriers, which may play a significant role in the compartmentalisation of groundwater flow, spring dynamics, and drawdown and recovery associated with groundwater extraction.

Hydrogeology Journal↗
USGS2020

High-resolution mapping of the freshwater-brine interface using deterministic and Bayesian inversion of airborne electromagnetic data at Paradox Valley, USA

Salt loads in the Colorado River Basin are a primary water quality concern. Natural groundwater brine discharge to the Dolores River where it passes through the collapsed salt anticline of the Paradox Valley in western Colorado is a significant source of salt to the Colorado River. An airborne electromagnetic survey of Paradox Valley has provided insights into the 3D distribution of brine in the surficial aquifer. A combination of stochastic and deterministic resistivity inversions were used to interpret the top of the freshwater-brine interface and to qualitatively describe the vertical salinity gradients across the interface. Low-resistivity regions indicative of brine occur near the land surface where brine discharges to the Dolores River and increase in depth several kilometers up-gradient along the axis of the valley. The most conductive parts of the brine plume are found in the areas below and adjacent to the river, suggesting that the brine becomes shallower and more concentrated as it reaches its natural discharge location. A significant freshwater lens overlying the brine west of the Dolores River is spatially correlated to the intermittent West Paradox Creek and agricultural irrigation. Below this lens, the transition from freshwater to brine appears to occur abruptly over a few meters and correlates to available well information. However, away from these regions and particularly with distance from the river, the freshwater-brine interface appears to be more diffuse.

Colorado↗
USGS2020

Seasonal manganese transport in the hyporheic zone of a snowmelt-dominated river (East River, Colorado)

Manganese (Mn) plays a critical role in river-water quality because Mn-oxides serve as sorption sites for contaminant metals. The aim of this study is to understand the seasonal cycling of Mn in an alpine streambed that experiences large spring snowmelt events and the potential responses to changes in snowmelt timing and magnitude. To address this goal, annual variations in river-water/groundwater interaction and Mn (aq) transport were measured and modeled in the bed of East River, Colorado, USA. In observations and numerical models, oxygenated river water containing dissolved organic carbon (DOC) mixes with groundwater rich in Mn (aq) in the streambed. The mixing depth increases during spring snowmelt when river discharge increases, leading to a greater DOC supply to the hyporheic zone and net respiration of Mn-oxides, despite an enhanced supply of oxygen. As groundwater upwelling resumes during the subsequent baseflow period, Mn (aq) -rich groundwater mixes with oxygenated river water, resulting in net accumulation of Mn-oxides until the bed freezes in winter. To explore potential responses of Mn transport to different climate-induced hydrological regimes, three hydrograph scenarios were numerically modeled (historic, low-snow, and storm) for the Rocky Mountain region. In a warming climate, Mn (aq) export to the river decreases, and Mn (aq) oxidation is favored in the upper streambed sediments over more of the year. One important implication is that the streambed may have an increased sorption capacity for metals over more of the year, leading to potential changes in river-water quality.

Colorado↗
USGS2020

The influence of snow cover, air temperature, and groundwater flow on the active-layer thermal regime of Arctic hillslopes drained by water tracks

Permafrost in Arctic watersheds limits soil biological activity to a thin, seasonally thawed active layer that contributes water to streams. In many hillslopes, relatively wet drainage features called water tracks have distinct freeze-thaw patterns that affect groundwater flow and storage, and thus the export of heat and solutes to Arctic streams. This study uses groundwater flow and energy transport models to examine potential controls on the timing and duration of freeze–thaw conditions and the magnitude of temperature fluctuations within water tracks and their adjacent hillslopes. The simulated length of the active-layer thaw season varies by 1 month over the range of snow-cover and mean annual air-temperature scenarios simulated. The timing and duration of freezing is particularly sensitive to depth and duration of snow cover. Thus, the deeper snowpack covers that can accumulate in water tracks contribute to their more persistent thaw conditions and their ability to conduct groundwater downslope. A three-dimensional simulation shows that during the summer thaw season, the water track captures groundwater laterally from half way across the hillslope. The models presented here elucidate key mechanisms driving small-scale variation in the active-layer thermal regime of tundra hillslopes, which may be responsible for changes in drainage-network geometry and Arctic biogeochemical fluxes under a warming climate.

Alaska↗
USGS2020

Geochemical characterization of groundwater evolution south of Grand Canyon, Arizona (USA)

Better characterization of the geochemical evolution of groundwater south of Grand Canyon, Arizona (USA), is needed to understand natural conditions and assess potential effects from breccia-pipe uranium mining in the region. Geochemical signatures of groundwater at 28 sampling locations were evaluated; baseline concentrations for select trace elements (As, B, Ba, Cr, Li, Mo, Rb, Se, Sr, Th, Tl, U, V) were established, and anomalous chemistry characteristics were identified. Concentrations at some groundwater sites exceeded the USEPA drinking water standard for As of 10 μg/L (Red Canyon, Miners, JT, Havasu, and Warm Springs) and U of 30 μg/L (Salt Creek Spring). Four springs from the study area (Blue, Havasu, Fern, and Warm Springs) had unique chemistry, which may indicate a deep flow path or potential contribution of fluids from lower in the crust. Other springs in the study area were distinguished by major anion water type: sulfate, bicarbonate, and a mixture of the two. Water type distinctions were somewhat spatially segregated, with sulfate type present on the western side of the study area, bicarbonate type on the eastern side, and a mixture of the two interspersed between the endmember sites. Sulfate-type water from this study area had low strontium isotopic ratio (87Sr/86Sr) values. The location of spring discharge within single drainages of the Grand Canyon may influence chemistry, as groundwater discharging from bedrock was altered after flowing through alluvial material. Geochemical analysis of groundwater in Grand Canyon indicates the importance of continued monitoring and better understanding of short-term chemical fluctuations.

Arizona↗
USGS2020

Rethinking groundwater flow on the South Rim of the Grand Canyon, USA: Characterizing recharge sources and flow paths with environmental tracers

In the arid landscape south of the Grand Canyon, natural springs and seeps are a critical resource for endemic species and Native American tribes. Groundwater is potentially threatened by expanding populations, visitations, and mineral extraction activities. Environmental tracers including noble gases, stable isotopes of hydrogen and oxygen in water, tritium, and carbon-14 were used to characterize recharge sources and flow paths in South Rim aquifers. Results confirm the regional Redwall-Muav Aquifer as the primary groundwater source to springs. However, a second local recharge source is required to explain the detection of tritium. Two probable sources are identified as: low-elevation infiltration of surface run-off with warm noble gas recharge temperatures, high excess air, and relatively low fractions of winter recharge, and high-elevation plateau recharge with cool recharge temperatures, low excess air, and fraction of winter recharge of ~ 1. Previous investigators have linked spring occurrence with regional faults and fractures. We show such features are also the likely control chemical mixing between the regional and local groundwater sources, the transport of deeply sourced and local recharge fluids, groundwater age, and thus the relative vulnerability of groundwater to depletion and contamination. The new conceptual model of groundwater sources and flow paths suggest many South Rim springs may respond on the order of 10s to 100s of years to groundwater depletion and contamination, even though the majority of groundwater flow is along longer flow paths with longer lag times. The magnitude of response to short term changes in the flow system remains unclear.

Arizona↗
USGS2020

Critical evaluation of stable isotope mixing end-members for estimating groundwater recharge sources: Case study from the South Rim of the Grand Canyon, Arizona, USA

Springs and groundwater seeps along the South Rim of the Grand Canyon serve an important function for the region’s ecosystems, residents (both human and wild animal), and economy. However, these springs and seeps are potentially vulnerable to contamination, increased groundwater extraction, or reduced recharge due to climate change. Protection of South Rim groundwater resources requires improved understanding of the regional groundwater system. In this study, statistical methods are used to investigate δ2H and δ18O in precipitation, surface water, and groundwater. A mixing model for δ18O is developed using statistically distinct seasonal end-members represented by modeled winter (Nov-Apr.) precipitation and summer (May-Oct.) surface water run-off. The calculated fraction of winter recharge (Fwin) indicates that South Rim groundwater is primarily sourced from snow-melt and winter rains with an average Fwin of 0.97 ± 0.09. Groundwater sourced from the highest elevations of the study area are more depleted than the winter end-member suggesting values of Fwin are overestimated or a meaningful portion of recharge occurs at lower elevations. Lower elevation recharge from the Coconino Plateau is supported by consistent spatial trends in δ2H and δ18O with respect to longitude, Fwin values less than 0.9 for 9 of the 50 samples, and age tracer data indicating young groundwater discharging from springs which is distinct from old groundwater observed in the regional flow system. These results suggest a new conceptual model is needed to account for recharge sources from low elevation and summer precipitation. Results imply resource managers need to reconsider current land-use and water management practices on the South Rim to protect future water quantity and quality.

Arizona↗
USGS2020

Evaluating the effects of downscaled climate projections on groundwater storage and simulated base-flow contribution to the North Fork Red River and Lake Altus, southwest Oklahoma (USA)

Potential effects of projected climate variability on base flow and groundwater storage in the North Fork Red River aquifer, Oklahoma (USA), were estimated using downscaled climate model data coupled with a numerical groundwater-flow model. The North Fork Red River aquifer discharges groundwater to the North Fork Red River, which provides inflow to Lake Altus. To approximate future conditions, Coupled Model Intercomparison Project Phase 5 climate data were downscaled to the watershed and a time-series of scaling factors were developed and interpolated for three climate scenarios (central tendency, warmer and drier, and less warm and wetter) representing future climate conditions for the period 2045–2074. These scaling factors were then applied to a soil-water-balance model to produce groundwater recharge and evapotranspiration estimates. A MODFLOW groundwater-flow model of the North Fork Red River aquifer used the scaled recharge and evapotranspiration data to estimate changes in base flow and water-surface elevation of Lake Altus. Compared to a baseline scenario, the mean percent change in annual base flow during 2045–2074 was −10.8 and −15.9% for the central tendency and warmer/drier scenarios, respectively; the mean percent change in annual base flow for the less-warm/wetter scenario was +15.7%. The mean annual percent change in groundwater storage for the central tendency, warmer/drier, and less-warm/wetter climate scenarios and the baseline are −2.7, −3.2, and +3.0%, respectively. The range of outcomes from the climate scenarios may be influenced by variability in the downscaled climate data for precipitation more than for temperature.

Oklahoma↗
USGS2021

Review: “Jacob’s Zoo”— How using Jacob’s method for aquifer testing leads to more intuitive understanding of aquifer characteristics

The interpretation of aquifer responses to pumping tests is an important tool for assessing aquifer geometry and properties, which are critical in the assessment of water resources or in environmental remediation. However, the responses of aquifers, measured by time-drawdown relationships in monitoring wells, are nonunique solutions that are affected by many factors. Jacob’s Zoo is a collection of graphical interpretations that allows students and practitioners to develop an intuitive feel for how natural hydrogeological systems work, and develop a set of skills that provide a better understanding of aquifer properties far beyond interpretation of pumping tests. Jacob’s Zoo, based on the work of Jacob (1950), fosters a deeper understanding, although few practitioners realize the full utility of the method. Jacob CE (1950) Flow of groundwater, In: Rouse H (ed) Engineering Hydraulics, Wiley, New York. P 321–386.

Hydrogeology Journal↗
USGS2021

Stratigraphic and structural controls on groundwater salinity variations in the Poso Creek Oil Field, Kern County, California, USA

Groundwater total dissolved solids (TDS) distribution was mapped with a three-dimensional (3D) model, and it was found that TDS variability is largely controlled by stratigraphy and geologic structure. General TDS patterns in the San Joaquin Valley of California (USA) are attributed to predominantly connate water composition and large-scale recharge from the adjacent Sierra Nevada. However, in smaller areas, stratigraphy and faulting play an important role in controlling TDS. Here, the relationship of stratigraphy and structure to TDS concentration was examined at Poso Creek Oil Field, Kern County, California. The TDS model was constructed using produced water TDS samples and borehole geophysics. The model was used to predict TDS concentration at discrete locations in 3D space and used a Gaussian process to interpolate TDS over a volume. In the overlying aquifer, TDS is typically <1,000 mg/L and increases with depth to ~1,200–3,500 mg/L in the hydrocarbon zone below the Macoma claystone—a regionally extensive, fine-grained unit—and reaches ~7,000 mg/L in isolated places. The Macoma claystone creates a vertical TDS gradient in the west where it is thickest, but control decreases to the east where it pinches out and allows freshwater recharge. Previously mapped normal faults were found to exhibit inconsistent control on TDS. In one case, high-density faulting appears to prevent recharge from flushing higher-TDS connate water. Elsewhere, the high-throw segments of a normal fault exhibit variable behavior, in places blocking lower-TDS recharge and in other cases allowing flushing. Importantly, faults apparently have differential control on oil and groundwater.

California↗
USGS2022

Identification of fresh submarine groundwater off the coast of San Diego, USA, using electromagnetic methods

Climate change has a pronounced effect on water resources in many semiarid climates, causing populated areas such as San Diego County (USA), to become more vulnerable to water shortages in the coming decades. To prepare for decreased water supply, San Diego County is adopting policies to decrease water use and to develop additional local sources of water. One new local source of freshwater is produced by a desalination facility that purifies brackish groundwater from the coastal San Diego Formation. This formation has been studied extensively onshore, but little is known about the geology or groundwater quality offshore in the adjacent continental shelf. Because most groundwater systems are interconnected and complex, further analysis is needed to identify offshore geology, possible sequestration of freshwater in the shelf, and potential pathways for saltwater intrusion. This comprehensive understanding is important because seawater intrusion may limit use of the San Diego Formation and longevity of desalination facilities. Controlled-source electromagnetic methods are uniquely suited to detecting offshore groundwater as they are sensitive to changes in pore fluids such as the transition from fresh to brackish groundwater. This paper describes results from surface-towed electromagnetic surveys that mapped the pore-fluid salinity and possible fluid pathways in the continental shelf off the coast of San Diego. The results indicate a considerable volume of fresh-to-brackish groundwater sequestered in the shelf, both in continuous lenses and isolated pockets, that appear influenced by fault systems and shallow stratigraphy.

Baja California, California↗
USGS2022

Effects of structure and volcanic stratigraphy on groundwater and surface water flow: Hat Creek basin, California, USA

Hydrogeologic systems in the southern Cascade Range in California (USA) develop in volcanic rocks where morphology, stratigraphy, extensional structures, and attendant basin geometry play a central role in groundwater flow paths, groundwater/surface-water interactions, and spring discharge locations. High-volume springs (greater than 3 m 3 /s) flow from basin-filling (<800 ka) volcanic rocks in the Hat Creek and Fall River tributaries and contribute approximately half of the average annual flow of the Pit River, the largest tributary to Shasta Lake. A hydrogeologic conceptual framework is constructed for the Hat Creek basin combining new geologic mapping, water-well lithologic logs, a database of active faults, LiDAR mapping of faults and volcanic landforms, streamflow measurements and airborne thermal infrared remote sensing of stream temperature. These data are used to integrate the geologic structure and the volcanic and volcaniclastic stratigraphy to create a three-dimensional interpretation of the hydrogeology in the basin. Two large streamflow gains from focused groundwater discharge near Big Spring and north of Sugarloaf Peak result from geologic barriers that restrict lateral groundwater flow and force water into Hat Creek. The inferred groundwater-flow barriers divide the aquifer system into at least three leaky compartments. The two downstream compartments lose streamflow in the upstream reaches (immediately downstream of the groundwater-flow barriers) and gain in downstream reaches with the greatest inflows immediately upstream of the barriers.

California↗
USGS2023

A model of transmissivity and hydraulic conductivity from electrical resistivity distribution derived from airborne electromagnetic surveys of the Mississippi River Valley Alluvial Aquifer, Midwest USA

Groundwater-flow models require the spatial distribution of the hydraulic conductivity parameter. One approach to defining this spatial distribution in groundwater-flow model grids is to map the electrical resistivity distribution by airborne electromagnetic (AEM) survey and establish a petrophysical relation between mean resistivity calculated as a nonlinear function of the resistivity layering and thicknesses of the layers and aquifer transmissivity compiled from historical aquifer tests completed within the AEM survey area. The petrophysical relation is used to transform AEM resistivity to transmissivity and to hydraulic conductivity over areas where the saturated thickness of the aquifer is known. The US Geological Survey applied this approach to a gain better understanding of the aquifer properties of the Mississippi River Valley alluvial aquifer. Alluvial-aquifer transmissivity data, compiled from 160 historical aquifer tests in the Mississippi Alluvial Plain (MAP), were correlated to mean resistivity calculated from 16,816 line-kilometers (km) of inverted resistivity soundings produced from a frequency-domain AEM survey of 95,000 km 2 of the MAP. Correlated data were used to define petrophysical relations between transmissivity and mean resistivity by omitting from the correlations the aquifer-test and AEM sounding data that were separated by distances greater than 1 km and manually calibrating the relation coefficients to slug-test data. The petrophysical relation yielding the minimum residual error between simulated and slug-test data was applied to 2,364 line-km of AEM soundings in the 1,000-km 2 Shellmound (Mississippi) study area to calculate hydraulic property distributions of the alluvial aquifer for use in future groundwater-flow models.

Mississippi River Valley Alluvial Aquifer↗
USGS2023

Quantifying stream-loss recovery in a spring using dual-tracer injections in the Snake Creek drainage, Great Basin National Park, Nevada, USA

Simultaneous short-pulse injections of two tracers (sodium bromide [Br – ] and fluorescein dye) were made in a losing reach of Snake Creek in Great Basin National Park, Nevada, USA, to evaluate the quantity of stream loss through permeable carbonates that resurfaces at a spring approximately 10 km down drainage. A revised hydrogeologic cross section for a possible flow path of the infiltrated Snake Creek water is presented, and the results may inform water management in the region. First arrival and peak concentration of the two tracers occurred at 9.5 and 12.7 days after injection, respectively. Fracture transport simulations indicate that Br – preferentially diffuses into immobile regions of the aquifer, and this diffusive flux is likely responsible for the major differences in the breakthrough curves. When considering the diffusive tracer flux, total apparent Br – and fluorescein dye recoveries were 16.9–22.1% and 21.7–24.3%, respectively. These findings imply that consideration of diffusive flux and long-term monitoring in fracture-dominated flow may support accurate quantification of tracer recovery. In addition, the apparent power law slopes of the breakthrough tails for both tracers were steeper at early times than have been attributed to heterogeneous advection or channeling in meter-scale tests, but the late-time Br – power law slope becomes less steep than has been attributed to diffusive exchange. These deviations may reflect fracture transport patterns that occur at larger scales.

Nevada↗
USGS2023

Detection and quantification of preferential flow using artificial rainfall with multiple experimental approaches

Preferential flow in the unsaturated zone strongly influences important hydrologic processes, such as infiltration, contaminant transport, and aquifer recharge. Because it entails various combinations of physical processes arising from the interactions of water, air, and solid particles in a porous medium, preferential flow is highly complex. Major research is needed to improve the ability to understand, quantify, model, and predict preferential flow. Toward a solution, a combination of diverse experimental measurements at multiple scales, from laboratory scale to mesoscale, has been implemented to detect and quantify preferential paths in carbonate and karstic unsaturated zones. This involves integration of information from (1) core samples, by means of mercury intrusion porosimeter, evaporation, quasi-steady centrifuge and dewpoint potentiometer laboratory methods, to investigate the effect of pore-size distribution on hydraulic characteristics and the potential activation of preferential flow, (2) field plot experiments with artificial sprinkling, to visualize preferential pathways related to secondary porosity, through use of geophysical measurements, and (3) mesoscale evaluation of field data through episodic master recession modeling of episodic recharge. This study demonstrates that preferential flow processes operate from core scale to two different field scales and impact on the qualitative and quantitative groundwater status, by entailing fast flow with subsequent effects on recharge rate and contaminant mobilizing. The presented results represent a rare example of preferential flow detection and numerical modeling by reducing underestimation of the recharge and contamination risks.

Bari↗
USGS2023

Arizona Groundwater Explorer: Interactive maps for evaluating the historical and current groundwater conditions in wells in Arizona, USA

Groundwater is an important water source in Arizona, accounting for about 41% of water use in this mostly arid-to-semiarid state in the southwestern United States, and the availability of groundwater resources in the state is a concern. To provide accessible information from depth-to-groundwater data, a series of web-based interactive maps were developed, called the Arizona Groundwater Explorer (AGEx). Scripts were written to harmonize and synthesize groundwater datasets from the two largest publicly available sources, subset these data to address different groundwater availability questions, and display the results in online, interactive maps. The combined dataset contained 1,820,122 depth-to-groundwater measurements from 1891 through 2022 from 41,918 wells in Arizona. Data views are provided for 20 topics, including recent (2020 or later) depth to groundwater (4,569 wells), historical (pre-1950) depth to groundwater (4,287 wells), wells with long-term (≥50 years) records (1,183 wells), wells with recent groundwater level decline (277 wells), wells with recent groundwater level rise (120 wells), and linear trends in groundwater levels over ten 10-year periods (number of wells ranging from 341 in 1978–1987 to 1,208 in 2003–2012), among others. With ongoing drought in the region resulting in declining surface-water supplies in Arizona, groundwater may play an even larger role in satisfying water needs in the state. The AGEx series of maps provides a nonspecialist audience with an improved understanding of historical, current, and changes in groundwater levels in Arizona.

Arizona↗
USGS2025

Controls on water quality below a reclaimed surface coal mine, southeastern Montana

Coal mining and reclamation can have a profound influence on hydrogeologic systems, with clear consequences for groundwater quality, yet their long-term influence on downgradient water quality over time following reclamation is less well documented. Geochemical trends were evaluated in water quality downgradient of a fully reclaimed landscape at the former Big Sky Mine in the Rosebud Creek watershed (southeastern Montana, USA), over a 3-year period (2020–2022), including bond release in 2022. Within 6 km downgradient from the reclaimed area, sulfate concentrations decreased from approximately 3500 to 1800 mg l −1 within the Miller Coulee alluvial aquifer. Major ions, δ 34 S SO4 values, and residence time tracers suggest that the observed decreases in sulfate concentration result from a combination of dilution by mixed-age inflows and incomplete transit of the high salinity plume from the mine boundary. Both bedrock and alluvial aquifers of the Rosebud Creek corridor contained contributions of millennia-old regional groundwater, which may serve to mitigate mine-derived high salinity waters. Rosebud Creek, which traverses the outflow zone of Miller Coulee in the study area, exhibited high sulfate concentrations during low flows and consistent downgradient increases in sulfate concentration. The possibility of plume dynamics in Miller Coulee suggests that the greatest water quality impacts may not yet have reached Rosebud Creek.

Montana↗