Geology Reports⌕ Search

Geology topics

R.L. Runkel

Publications and source records attributed to R.L. Runkel.

16 recordsLinked to original sources
USGS2010

An approach to quantify sources, seasonal change, and biogeochemical processes affecting metal loading in streams: Facilitating decisions for remediation of mine drainage

Historical mining has left complex problems in catchments throughout the world. Land managers are faced with making cost-effective plans to remediate mine influences. Remediation plans are facilitated by spatial mass-loading profiles that indicate the locations of metal mass-loading, seasonal changes, and the extent of biogeochemical processes. Field-scale experiments during both low- and high-flow conditions and time-series data over diel cycles illustrate how this can be accomplished. A low-flow experiment provided spatially detailed loading profiles to indicate where loading occurred. For example, SO 4 2 - was principally derived from sources upstream from the study reach, but three principal locations also were important for SO 4 2 - loading within the reach. During high-flow conditions, Lagrangian sampling provided data to interpret seasonal changes and indicated locations where snowmelt runoff flushed metals to the stream. Comparison of metal concentrations between the low- and high-flow experiments indicated substantial increases in metal loading at high flow, but little change in metal concentrations, showing that toxicity at the most downstream sampling site was not substantially greater during snowmelt runoff. During high-flow conditions, a detailed temporal sampling at fixed sites indicated that Zn concentration more than doubled during the diel cycle. Monitoring programs must account for diel variation to provide meaningful results. Mass-loading studies during different flow conditions and detailed time-series over diel cycles provide useful scientific support for stream management decisions.

Applied Geochemistry↗
USGS2009

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.

Colorado↗
USGS2009

Spatially detailed quantification of metal loading for decision making: Metal mass loading to American fork and Mary Ellen Gulch, Utah

Effective remediation requires an understanding of the relative contributions of metals from all sources in a catchment, and that understanding must be based on a spatially detailed quantification of metal loading. A traditional approach to quantifying metal loading has been to measure discharge and chemistry at a catchment outlet. This approach can quantify annual loading and the temporal changes in load, but does not provide the needed spatial detail to evaluate specific sources, which is needed to support remediation decisions. A catchment or mass-loading approach provides spatial detail by combining tracer-injection and synoptic-sampling methods to quantify loading. Examples of studies in American Fork, Utah, and its tributary Mary Ellen Gulch illustrate this different approach. The mass-loading study in American Fork treated Mary Ellen Gulch as a single inflow. From that point of view, Mary Ellen Gulch was one of the greatest sources of Fe, Mn, Zn, and colloidal Pb loads to American Fork. But when Mary Ellen Gulch was evaluated in a separate catchment study, the detailed locations of metal loading were identified, and the extent of metal attenuation upstream from the mouth of Mary Ellen Gulch was quantified. The net, instantaneous load measured at the mouth of Mary Ellen Gulch for remediation planning would greatly underestimate the contributions of principal sources within the catchment. Extending the detailed sampling downstream from Mary Ellen Gulch indicated the possibility of diffuse groundwater inflow from Mary Ellen Gulch to American Fork. Comparing loads for Mary Ellen Gulch in the two studies indicates that metal loads could be substantially underestimated for planning purposes without the detailed catchment approach for the low-flow conditions in these studies. A mass-loading approach provides both the needed quantification of metal loading and the spatial detail to guide remediation decisions that would be the most effective in the catchments.

Mine Water and the Environment↗
USGS2009

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, δ 56 Fe and δ 66 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 δ 56 Fe were linked to redox and precipitation reactions occurring in the stream, while changes in δ 66 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 ( δ 56 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.

Applied Geochemistry↗
USGS2006

Multiple injected and natural conservative tracers quantify mixing in a stream confluence affected by acid mine drainage near Silverton, Colorado

The acidic discharge from Cement Creek, containing elevated concentrations of dissolved metals and sulphate, mixed with the circumneutral-pH Animas River over a several hundred metre reach (mixing zone) near Silverton, CO, during this study. Differences in concentrations of Ca, Mg, Si, Sr, and SO42- between the creek and the river were sufficiently large for these analytes to be used as natural tracers in the mixing zone. In addition, a sodium chloride (NaCl) tracer was injected into Cement Creek, which provided a Cl- 'reference' tracer in the mixing zone. Conservative transport of the dissolved metals and sulphate through the mixing zone was verified by mass balances and by linear mixing plots relative to the injected reference tracer. At each of seven sites in the mixing zone, five samples were collected at evenly spaced increments of the observed across-channel gradients, as determined by specific conductance. This created sets of samples that adequately covered the ranges of mixtures (mixing ratios, in terms of the fraction of Animas River water, %AR). Concentrations measured in each mixing zone sample and in the upstream Animas River and Cement Creek were used to compute %AR for the reference and natural tracers. Values of %AR from natural tracers generally showed good agreement with values from the reference tracer, but variability in discharge and end-member concentrations and analytical errors contributed to unexpected outlier values for both injected and natural tracers. The median value (MV) %AR (calculated from all of the tracers) reduced scatter in the mixing plots for the dissolved metals, indicating that the MV estimate reduced the effects of various potential errors that could affect any tracer.

Hydrological Processes↗
USGS2006

Evaluation of metal loading to streams near Creede, Colorado, August and September 2000

Decisions about remediation of mine drainage on the watershed scale require an understanding of metal contributions from all sources to be able to choose the best sites for remediation. A hydrologic framework to study metal loading in the Willow Creek watershed, a tributary to the Rio Grande River, was established by conducting a series of tracer-injection studies. Each study used the tracer-dilution method in conjunction with synoptic sampling to determine the spatial distribution of discharge and concentration. Discharge and concentration data were then used to develop mass-loading curves for the metals of interest. The discharge and load profiles (1) identify the principal sources of load to the streams; (2) demonstrate the scale of unsampled, dispersed subsurface inflows; and (3) estimate the amount of natural attenuation. The greatest source of metal loads was from the Nelson Tunnel on West Willow Creek, which contributed 158 kilograms per day of zinc to the stream. Additional loading from other dispersed, subsurface inflows along West Willow Creek added substantial loads, but these were small in comparison to the loads from the Nelson Tunnel. No significant contributions of metal load from potential sources occurred along East Willow Creek. The lack of measurable loading may be a result of previous remedial actions along that stream. The lower Willow Creek section had relatively small contributions of load compared to what had been contributed upstream. This watershed approach provides a detailed snapshot of metal load for the watershed to support remediation decisions and quantifies processes that affect metal transport.

Colorado↗
USGS2004

Denitrification and hydrologic transient storage in a glacial meltwater stream, McMurdo Dry Valleys, Antarctica

In extreme environments, retention of nutrients within stream ecosystems contributes to the persistence of aquatic biota and continuity of ecosystem function. In the McMurdo Dry Valleys, Antarctica, many glacial meltwater streams flow for only 5–12 weeks a year and yet support extensive benthic microbial communities. We investigated NO 3 − uptake and denitrification in Green Creek by analyzing small‐scale microbial mat dynamics in mesocosms and reach‐scale nutrient cycling in two whole‐stream NO 3 − enrichment experiments. Nitrate uptake results indicated that microbial mats were nitrogen (N)‐limited, with NO 3 − uptake rates as high as 16 nmol N cm −2 h −1 . Denitrification potentials associated with microbial mats were also as high as 16 nmol N cm −2 h −1 . During two whole‐stream NO 3 − −enrichment experiments, a simultaneous pulse of NO 2 − was observed in the stream water. The one‐dimensional solute transport model with inflow and storage was modified to simulate two storage zones: one to account for short time scale hydrologic exchange of stream water into and out of the benthic microbial mat, the other to account for longer time scale hydrologic exchange with the hyporheic zone. Simulations indicate that injected NO 3 − was removed both in the microbial mat and in the hyporheic zone and that as much as 20% of the NO 3 − that entered the microbial mat and hyporheic zone was transformed to NO 2 − by dissimilatory reduction. Because of the rapid hydrologic exchange in microbial mats, it is likely that denitrification is limited either by biotic assimilation, reductase limitation, or transport limitation (reduced NO 2 − is transported away from reducing microbes).

Limnology and Oceanography↗
USGS2004

Fate of volatile organic compounds in constructed wastewater treatment wetlands

The fate of volatile organic compounds was evaluated in a wastewater-dependent constructed wetland near Phoenix, AZ, using field measurements and solute transport modeling. Numerically based volatilization rates were determined using inverse modeling techniques and hydraulic parameters established by sodium bromide tracer experiments. Theoretical volatilization rates were calculated from the two-film method incorporating physicochemical properties and environmental conditions. Additional analyses were conducted using graphically determined volatilization rates based on field measurements. Transport (with first-order removal) simulations were performed using a range of volatilization rates and were evaluated with respect to field concentrations. The inverse and two-film reactive transport simulations demonstrated excellent agreement with measured concentrations for 1,4-dichlorobenzene, tetrachloroethene, dichloromethane, and trichloromethane and fair agreement for dibromochloromethane, bromodichloromethane, and toluene. Wetland removal efficiencies from inlet to outlet ranged from 63% to 87% for target compounds.

Environmental Science & Technology↗
USGS2004

Evaluating remedial alternatives for the Alamosa River and Wightman Fork, near Summitville Mine, Colorado: Application of a reactive transport model to low- and high-flow simulations

Reactive-transport processes in Wightman Fork and the Alamosa River downstream from the Summitville Mine, south-central Colorado, were simulated at low and high flow using the OTEQ reactive-transport model. The simulations were calibrated using data from synoptic studies conducted during October 1998 and June 1999. Discharge over the 30-km reach from just below the mine site to the Alamosa River above Terrace Reservoir ranged from 0.077 to 1.3 m3/s at low flow and from 1.17 to 17.0 m3/s at high flow. Travel time was about 28 hours at low flow and about 8.5 hours at high flow; pH ranged from 4.6 to 5.7 at low flow and from 3.7 to 6.7 at high flow. Simulations revealed that pH, Fe, Al, and Cu were non-conservative. Simulations included Fe(II) oxidation, constrained using measured values of Fe(II) and Fe(total). Precipitation of hydrous Fe oxides and hydrous Al oxides and hydroxysulfates match observed conditions more closely in simulations that included Fe(II) oxidation and Fe(III) precipitation than in simulations without Fe(II) oxidation or Fe(III) precipitation. Simulation results indicate that sorption is controlling Cu concentrations in the Alamosa River. The calibrated models were used to evaluate nine remediation alternatives.

Book chapter↗
USGS2003

Automated calibration of a stream solute transport model: Implications for interpretation of biogeochemical parameters

The hydrologic processes of advection, dispersion, and transient storage are the primary physical mechanisms affecting solute transport in streams. The estimation of parameters for a conservative solute transport model is an essential step to characterize transient storage and other physical features that cannot be directly measured, and often is a preliminary step in the study of reactive solutes. Our study used inverse modeling to estimate parameters of the transient storage model OTIS (One dimensional Transport with Inflow and Storage). Observations from a tracer injection experiment performed on Uvas Creek, California, USA, are used to illustrate the application of automated solute transport model calibration to conservative and nonconservative stream solute transport. A computer code for universal inverse modeling (UCODE) is used for the calibrations. Results of this procedure are compared with a previous study that used a trial-and-error parameter estimation approach. The results demonstrated 1) importance of the proper estimation of discharge and lateral inflow within the stream system; 2) that although the fit of the observations is not much better when transient storage is invoked, a more randomly distributed set of residuals resulted (suggesting non-systematic error), indicating that transient storage is occurring; 3) that inclusion of transient storage for a reactive solute (Sr2+) provided a better fit to the observations, highlighting the importance of robust model parameterization; and 4) that applying an automated calibration inverse modeling estimation approach resulted in a comprehensive understanding of the model results and the limitation of input data.

Journal of the North American Benthological Societ↗
USGS2002

A new metric for determining the importance of transient storage

A review of various metrics used to characterize transient storage indicates that none of the existing measures successfully integrate the interaction between advective velocity and the transient storage parameters (storage zone area, storage zone exchange coefficient). Further, 2 existing metrics are related to mean travel time, a quantity that is independent of the storage zone exchange coefficient, ??. This interaction and the effect of ?? on travel time are important considerations when determining the mass of solute entering the storage zone within a given reach. A new metric based on median reach travel time is therefore proposed. Median reach travel time due to advection-dispersion and transient storage, and median reach travel time due solely to advection-dispersion are computed based on numerical simulations. These 2 travel times are used to determine Fmed, the fraction of median travel time due to transient storage. Application of the new metric to 53 existing parameter sets indicates that transient storage accounts for 0.12% to 68.0% of total reach travel time. Rankings of storage zone importance based on the new metric are substantially different from rankings based on storage zone residence time, storage exchange flux, and the hydrological retention factor. These differences result from the ability of the new metric to characterize the interaction between advecrive velocity and transient storage, and the resultant effects on reach travel time and mass transport.

Journal of the North American Benthological Societ↗
USGS2001

pH dependence of iron photoreduction in a rocky mountain stream affected by acid mine drainage

The redox speciation of dissolved iron and the transport of iron in acidic, metal‐enriched streams is controlled by precipitation and dissolution of iron hydroxides, by photoreduction of dissolved ferric iron and hydrous iron oxides, and by oxidation of the resulting dissolved ferrous iron. We examined the pH dependence of these processes in an acidic mine‐drainage stream, St Kevin Gulch, Colorado, by experimentally increasing the pH of the stream from about 4·0 to 6·5 and following the downstream changes in iron species. We used a solute transport model with variable flow to evaluate biogeochemical processes controlling downstream transport. We found that at pH 6·4 there was a rapid and large initial loss of ferrous iron concurrent with the precipitation of aluminium hydroxide. Below this reach, ferrous iron was conservative during the morning but there was a net downstream loss of ferrous iron around noon and in the afternoon. Calculation of net oxidation rates shows that the noontime loss rate was generally much faster than rates for the ferrous iron oxidation at pH 6 predicted by Singer and Stumm (1970. Science 167 : 1121). The maintenance of ferrous iron concentrations in the morning is explained by the photoreduction of photoreactive ferric species, which are then depleted by noon.

Hydrological Processes↗
USGS1998

Analysis of transient storage subject to unsteady flow: Diel flow variation in an Antarctic stream

Transport of dissolved material in streams and small rivers may be characterized using tracer-dilution methods and solute transport models. Recent studies have quantified stream/substream interactions using models of transient storage. These studies are based on tracer-dilution data obtained during periods of steady flow. We present a modeling framework for the analysis of transient storage in stream systems with unsteady flows. The framework couples a kinematic wave routing model with a solute transport model that includes transient storage. The routing model provides time-varying flows and cross-sectional areas that are used as input to the solute transport model. The modeling framework was used to quantify stream/substream interaction in Huey Creek, an Antarctic stream fed exclusively by glacial meltwater. Analysis of tracer-dilution data indicates that there was substantial interaction between the flowing surface water and the hyporheic (substream) zone. The ratio of storage zone area to stream cross-sectional area (A<sub>s</sub>/A) was >1 in all stream reaches, indicating that the substream area contributing to hyporheic exchange was large relative to stream cross-sectional area. The rate of exchange, as governed by the transient storage exchange coefficient (α), was rapid because of a high stream gradient and porous alluvial materials. Estimates of α generally exceed those determined for other small streams. The high degree of hyporheic exchange supports the hypothesis that weathering reactions within the hyporheos account for observed increases in solute concentration with stream length, as noted in other studies of Antarctic streams.

Journal of the North American Benthological Societ↗
USGS1996

Reactive solute transport in an acidic stream: Experimental pH increase and simulation of controls on pH, aluminum, and iron

Solute transport simulations quantitatively constrained hydrologic and geochemical hypotheses about field observations of a pH modification in an acid mine drainage stream. Carbonate chemistry, the formation of solid phases, and buffering interactions with the stream bed were important factors in explaining the behavior of pH, aluminum, and iron. The precipitation of microcrystalline gibbsite accounted for the behavior of aluminum; precipitation of Fe(OH)3 explained the general pattern of iron solubility. The dynamic experiment revealed limitations on assumptions that reactions were controlled only by equilibrium chemistry. Temporal variation in relative rates of photoreduction and oxidation influenced iron behavior. Kinetic limitations on ferrous iron oxidation and hydrous oxide precipitation and the effects of these limitations on field filtration were evident. Kinetic restraints also characterized interaction between the water column and the stream bed, including sorption and desorption of protons from iron oxides at the sediment-water interface and post-injection dissolution of the precipitated aluminum solid phase.

Environmental Science & Technology↗
USGS1995

Simulation models for conservative and nonconservative solute transport in streams

Solute transport in streams is governed by a suite of hydrologic and chemical processes. Interactions between hydrologic processes and chemical reactions may be quantified through a combination of field-scale experimentation and simulation modeling. Two mathematical models that simulate conservative and nonconservative solute transport in streams are presented. A model for conservative solutes that considers One Dimensional Transport with Inflow and Storage (OTIS) may be used in conjunction with tracer-dilution methods to quantify hydrologic transport processes (advection, dispersion, lateral inflow and transient storage). For nonconservative solutes, a model known as OTEQ may be used to quantify chemical processes within the context of hydrologic transport. OTEQ combines the transport mechanisms in OTIS with a chemical equilibrium sub-model that considers complexation, precipitation/dissolution and sorption. OTEQ has been used to quantify processes affecting trace metals in two streams in the Rocky Mountains of Colorado, USA.

Effects of scale on interpretation and management ↗