A multidisciplinary-based conceptual model of a fractured sedimentary bedrock aquitard: improved prediction of aquitard integrity

A hydrogeologic conceptual model that improves understanding of variability in aquitard integrity is presented for a fractured sedimentary bedrock unit in the Cambrian-Ordovician aquifer system of midcontinent North America. The model is derived from multiple studies on the siliciclastic St. Lawrence Formation and adjacent strata across a range of scales and geologic conditions. These studies employed multidisciplinary techniques including borehole flowmeter logging, high-resolution depth-discrete multilevel well monitoring, fracture stratigraphy, fluorescent dye tracing, and three-dimensional (3D) distribution of anthropogenic tracers regionally. The paper documents a bulk aquitard that is highly anisotropic because of poor connectivity of vertical fractures across matrix with low permeability, but with ubiquitous bed parallel partings. The partings provide high bulk horizontal hydraulic conductivity, analogous to aquifers in the system, while multiple preferential termination horizons of vertical fractures serve as discrete low vertical hydraulic conductivity intervals inhibiting vertical flow. The aquitard has substantial variability in its ability to protect underlying groundwater from contamination. Across widespread areas where the aquitard is deeply buried by younger bedrock, preferential termination horizons provide for high aquitard integrity (i.e. protection). Protection is diminished close to incised valleys where stress release and weathering has enhanced secondary pore development, including better connection of fractures across these horizons. These conditions, along with higher hydraulic head gradients in the same areas and more complex 3D flow where the aquitard is variably incised, allow for more substantial transport to deeper aquifers. The conceptual model likely applies to other fractured sedimentary bedrock aquitards within and outside of this region.     

A comparison of pre‐ and post‐remediation water quality,Mineral Creek,Colorado

Pre‐ and post‐remediation data sets are used herein to assess the effectiveness of remedial measures implemented in the headwaters of the Mineral Creek watershed, where contamination from hard rock mining has led to elevated metal concentrations and acidic pH. Collection of pre‐ and post‐remediation data sets generally followed the synoptic mass balance approach, in which numerous stream and inflow locations are sampled for the constituents of interest and estimates of streamflow are determined by tracer dilution. The comparison of pre‐ and post‐remediation data sets is confounded by hydrologic effects and the effects of temporal variation. Hydrologic effects arise due to the relatively wet conditions that preceded the collection of pre‐remediation data, and the relatively dry conditions associated with the post‐remediation data set. This difference leads to a dilution effect in the upper part of the study reach, where pre‐remediation concentrations were diluted by rainfall, and a source area effect in the lower part of the study reach, where a smaller portion of the watershed may have been contributing constituent mass during the drier post‐remediation period. A second confounding factor, temporal variability, violates the steady‐state assumption that underlies the synoptic mass balance approach, leading to false identification of constituent sources and sinks. Despite these complications, remedial actions completed in the Mineral Creek headwaters appear to have led to improvements in stream water quality, as post‐remediation profiles of instream load are consistently lower than the pre‐remediation profiles over the entire study reach for six of the eight constituents considered (aluminium, arsenic, cadmium, copper, iron, and zinc). Concentrations of aluminium, cadmium, copper, lead, and zinc remain above chronic aquatic‐life standards, however, and additional remedial actions may be needed. Future implementations of the synoptic mass balance approach should be preceded by an assessment of temporal variability, and modifications to the synoptic sampling protocol should be made if necessary. Published in 2009 by John Wiley & Sons, Ltd.     

Classification of Thermal Patterns at Karst Springs and Cave Streams

Thermal patterns of karst springs and cave streams provide potentially useful information concerning aquifer geometry and recharge. Temperature monitoring at 25 springs and cave streams in southeastern Minnesota has shown four distinct thermal patterns. These patterns can be divided into two types: those produced by flow paths with ineffective heat exchange, such as conduits, and those produced by flow paths with effective heat exchange, such as small fractures and pore space. Thermally ineffective patterns result when water flows through the aquifer before it can equilibrate to the rock temperature. Thermally ineffective patterns can be either event‐scale, as produced by rainfall or snowmelt events, or seasonal scale, as produced by input from a perennial surface stream. Thermally effective patterns result when water equilibrates to rock temperature, and the patterns displayed depend on whether the aquifer temperature is changing over time. Shallow aquifers with seasonally varying temperatures display a phase‐shifted seasonal signal, whereas deeper aquifers with constant temperatures display a stable temperature pattern. An individual aquifer may display more than one of these patterns. Since karst aquifers typically contain both thermally effective and ineffective routes, we argue that the thermal response is strongly influenced by recharge mode.     

Determining long time‐scale hyporheic zone flow paths in Antarctic streams

In the McMurdo Dry Valleys of Antarctica, glaciers are the source of meltwater during the austral summer, and the streams and adjacent hyporheic zones constitute the entire physical watershed; there are no hillslope processes in these systems. Hyporheic zones can extend several metres from each side of the stream, and are up to 70 cm deep, corresponding to a lateral cross‐section as large as 12 m², and water resides in the subsurface year around. In this study, we differentiate between the near‐stream hyporheic zone, which can be characterized with stream tracer experiments, and the extended hyporheic zone, which has a longer time‐scale of exchange. We sampled stream water from Green Creek and from the adjacent saturated alluvium for stable isotopes of D and ¹⁸O to assess the significance and extent of stream‐water exchange between the streams and extended hyporheic zones over long time‐scales (days to weeks). Our results show that water residing in the extended hyporheic zone is much more isotopically enriched (up to 11‰ D and 2·2‰ ¹⁸O) than stream water. This result suggests a long residence time within the extended hyporheic zone, during which fractionation has occurred owing to summer evaporation and winter sublimation of hyporheic water. We found less enriched water in the extended hyporheic zone later in the flow season, suggesting that stream water may be exchanged into and out of this zone, on the time‐scale of weeks to months. The transient storage model OTIS was used to characterize the exchange of stream water with the extended hyporheic zone. Model results yield exchange rates (α) generally an order magnitude lower (10^?5 s^?1) than those determined using stream‐tracer techniques on the same stream. In light of previous studies in these streams, these results suggest that the hyporheic zones in Antarctic streams have near‐stream zones of rapid stream‐water exchange, where ‘fast’ biogeochemical reactions may influence water chemistry, and extended hyporheic zones, in which slower biogeochemical reaction rates may affect stream‐water chemistry at longer time‐scales. Copyright © 2003 John Wiley & Sons, Ltd.     

Application of iron and zinc isotopes to track the sources and mechanisms of metal loading in a mountain watershed

Here the hydrogeochemical constraints of a tracer dilution study are combined with Fe and Zn isotopic measurements to pinpoint metal loading sources and attenuation mechanisms in an alpine watershed impacted by acid mine drainage. In the tested mountain catchment, δ⁵⁶Fe and δ⁶⁶Zn isotopic signatures of filtered stream water samples varied by ∼3.5‰ and 0.4‰, respectively. The inherent differences in the aqueous geochemistry of Fe and Zn provided complimentary isotopic information. For example, variations in δ⁵⁶Fe were linked to redox and precipitation reactions occurring in the stream, while changes in δ⁶⁶Zn were indicative of conservative mixing of different Zn sources. Fen environments contributed distinctively light dissolved Fe (<−2.0‰) and isotopically heavy suspended Fe precipitates to the watershed, while Zn from the fen was isotopically heavy (>+0.4‰). Acidic drainage from mine wastes contributed heavier dissolved Fe (∼+0.5‰) and lighter Zn (∼+0.2‰) isotopes relative to the fen. Upwelling of Fe-rich groundwater near the mouth of the catchment was the major source of Fe (δ⁵⁶Fe ∼ 0‰) leaving the watershed in surface flow, while runoff from mining wastes was the major source of Zn. The results suggest that given a strong framework for interpretation, Fe and Zn isotopes are useful tools for identifying and tracking metal sources and attenuation mechanisms in mountain watersheds.