Can Karst Conduit Models Be Calibrated? A Dual Approach Using Dye Tracing and Temperature

Modeling flow and transport using both temperature and dye tracing provides constraints that can improve understanding of karst networks. A laminar flow and transport model using the finite element subsurface flow model simulated the conduit connection between a sinking stream and spring in central Pennsylvania to evaluate how conduit morphology might affect dye transport. Single and overly tortuous conduit models resulted in high concentrations as dye flowed back into the conduit from the matrix after dye injections ceased. A forked conduit model diverted flow from the main conduit, reducing falling limb dye concentration. Latin hypercube sampling was performed to evaluate the sensitivity of 52 parameter combinations (conduit hydraulic conductivity, conduit cross-sectional area, matrix transmissivity, matrix porosity, and dispersivity) for four conduit geometry scenarios. Sensitivity of arrival time for 50% of the dye indicated no parameter combinations which simulate falling limb dye concentrations for tortuous geometries, confirming the importance of the forked geometry regardless of other parameters. Temperature data from high-resolution loggers were then incorporated into the forked conduit model to reproduce seasonal spring temperature using variable sink inflow. Unlike the dye trace models, the thermal models were sensitive to other model parameters, such as conduit cross-sectional area and matrix transmissivity. These results showed this dual approach (dye and temperature) to karst network modeling is useful for (1) exploring the role of conduit and matrix interaction for contaminant storage, (2) constraining karst conduit geometries, which are often poorly understood, and (3) quantifying the effect of seasonal trends on karst aquifers.     

Observing lingering hyporheic storage using electrical resistivity: variations around stream restoration structures,Crabby Creek,PA

Time‐lapse geophysical surveys can map lingering hyporheic storage by detecting changes in response to saline tracer. Tracer tests were conducted in Crabby Creek, an urban stream outside Philadelphia, to examine the influence of stream restoration structures and variable sediment thickness. We compared electrical resistivity surveys with extensive well sampling (57 wells) in two 13.5‐m‐long reaches, each with a step drop created by a J‐hook. The two step drops varied in tracer behaviour, based on both the well data and the geophysical data. The well data showed more variation in arrival time where the streambed sediment was thick and was more uniform where sediment was thin. The resistivity in the reach with thin sediment showed lingering tracer in the hyporheic zone both upstream and downstream from the J‐hook. In the second reach where the sediment was thicker, the lingering tracer in the hyporheic zone was more extensive downstream from the J‐hook. The contrasting results between the two reaches from both methods suggested that sediments influenced hyporheic exchange more than the step at this location. Resistivity inversion differed from well data in both reaches in that it showed evidence for tracer after well samples had returned to background, mapping lingering tracer either upstream or downstream of a step. We conclude that resistivity surveys may become an important tool for hyporheic zone characterization because they provide information on the extent of slow moving fluids in the hyporheic zone, which have the potential to enhance chemical reactions. Copyright © 2012 John Wiley & Sons, Ltd.     

Characterizing snowmelt anomalies in hydrochemographs of a karst spring,Cumberland Valley,Pennsylvania (USA): evidence for multiple recharge pathways

Conductivity, water level, air temperature, and depth of snowpack were monitored during a 26-day melt period of 88-cm-deep snowpack at a karst spring to characterize internal runoff and diffuse infiltration. Chloride from road salt provided a tracer and the snowpack a recharge source during the melt period. The melt period was divided into phases based on air temperature and chemograph pattern. For the first and third phases, mean air temperatures were below freezing, but above freezing during the second and fourth phases. During the first phase when the temperature peaked above freezing, conductivity typically spiked 10–50 μS/cm, suggesting input of road salt from conduits. When the snowpack continuously melted, conductivity and water-level trends were upward with smaller daily spikes in conductivity indicating infiltration from the dilute snowpack. This pattern suggests that road salt input continued when snowmelt recharged through the epikarst, but at lower concentrations than the conduit input. Refreezing of the snowpack and shallow subsurface for a brief period interrupted the recharge, and there was no longer a sawtooth pattern of conductivity. It is apparent that frozen conditions did not cease recharge because a dual recharge process was evident. While dual recharge from internal runoff and diffuse infiltration occurred, the portions varied because of changing melting rates. Observed patterns indicated internal runoff dominated during frozen periods because recharge water moved as overland flow across a frozen surface to focused pathways. Diffuse infiltration became available during warmer periods because subsurface thawing allowed the snowmelt to penetrate the epikarst. Results of snowmelt monitoring in spring discharge indicated the transient nature of karst recharge.     

Modeling alternative paths of chemical evolution of Na-HCO3-type groundwater near Oak Ridge, Tennessee, USA

 This paper demonstrates that both cation exchange, a commonly invoked mechanism, and silicate hydrolysis, which is less commonly considered, can produce Na-HCO₃-type water in sedimentary rocks. Evolution of Na-HCO₃ groundwater beneath the Oak Ridge Reservation, Tennessee, USA, was studied by comparing observed end-member groundwater composition from multiport samplers to compositions generated by reaction-path geochemical models. Observed groundwater compositions could be reproduced by either the silicate-hydrolysis model or the cation-exchange model. Secondary minerals precipitated in the silicate-hydrolysis model are similar to those present along fractures in the shale and carbonate host rocks, and observed molar Sr²⁺/Ca²⁺ ratios more closely resemble evolution from shale weathering. Both mechanisms should be considered to understand the evolution of Na-HCO₃ groundwater. Received, April 1998 · Revised, January 1999 · Accepted, March 1999     

Geophysical and Hydrologic Studies of Lake Seepage Variability

Variations in lake seepage were studied along a 130 m shoreline of Mirror Lake NH. Seepage was downward from the lake to groundwater; rates measured from 28 seepage meters varied from 0 to ?282 cm/d. Causes of this variation were investigated using electrical resistivity surveys and lakebed sediment characterization. Two‐dimensional (2D) resistivity surveys showed a transition in lakebed sediments from outwash to till that correlated with high‐ and low‐seepage zones, respectively. However, the 2D survey was not able to predict smaller scale variations within these facies. In the outwash, fast seepage was associated with permeability variations in a thin (2 cm) layer of sediments at the top of the lakebed. In the till, where seepage was slower than that in the outwash, a three‐dimensional resistivity survey mapped a point of high seepage associated with heterogeneity (lower resistivity and likely higher permeability). Points of focused flow across the sediment–water interface are difficult to detect and can transmit a large percentage of total exchange. Using a series of electrical resistivity geophysical methods in combination with hydrologic data to locate heterogeneities that affect seepage rates can help guide seepage meter placement. Improving our understanding of the causes and types of heterogeneity in lake seepage will provide better data for lake budgets and prediction of mass transfer of solutes or contaminants between lakes and groundwater.