USGS Science⌕ Search

Geology topics

R.L. Runkel

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

At least 19 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

Naturally acidic surface and ground waters draining porphyry-related mineralized areas of the Southern Rocky Mountains, Colorado and New Mexico

Acidic, metal-rich waters produced by the oxidative weathering and resulting leaching of major and trace elements from pyritic rocks can adversely affect water quality in receiving streams and riparian ecosystems. Five study areas in the southern Rocky Mountains with naturally acidic waters associated with porphyry mineralization were studied to document variations in water chemistry and processes that control the chemical variations. Study areas include the Upper Animas River watershed, East Alpine Gulch, Mount Emmons, and Handcart Gulch in Colorado and the Red River in New Mexico. Although host-rock lithologies in all these areas range from Precambrian gneisses to Cretaceous sedimentary units to Tertiary volcanic complexes, the mineralization is Tertiary in age and associated with intermediate to felsic composition, porphyritic plutons. Pyrite is ubiquitous, ranging from ???1 to >5 vol.%. Springs and headwater streams have pH values as low as 2.6, SO4 up to 3700 mg/L and high dissolved metal concentrations (for example: Fe up to 400 mg/L; Cu up to 3.5 mg/L; and Zn up to 14.4 mg/L). Intensity of hydrothermal alteration and presence of sulfides are the primary controls of water chemistry of these naturally acidic waters. Subbasins underlain by intensely hydrothermally altered lithologies are poorly vegetated and quite susceptible to storm-induced surface runoff. Within the Red River study area, results from a storm runoff study documented downstream changes in river chemistry: pH decreased from 7.80 to 4.83, alkalinity decreased from 49.4 to <1 mg/L, SO4 increased from 162 to 314 mg/L, dissolved Fe increased from to 0.011 to 0.596 mg/L, and dissolved Zn increased from 0.056 to 0.607 mg/L. Compared to mine drainage in the same study areas, the chemistry of naturally acidic waters tends to overlap but not reach the extreme concentrations of metals and acidity as some mine waters. The chemistry of waters draining these mineralized but unmined areas can be used to estimate premining conditions at sites with similar geologic and hydrologic conditions. For example, the US Geological Survey was asked to estimate premining ground-water chemistry at the Questa Mo mine, and the proximal analog approach was used because a mineralized but unmined area was located adjacent to the mine property. By comparing and contrasting water chemistry from different porphyry mineralized areas, this study not only documents the range in concentrations of constituents of interest but also provides insight into the primary controls of water chemistry.

Applied Geochemistry↗
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↗
USGS2009

Copper isotope fractionation in acid mine drainage

We measured the Cu isotopic composition of primary minerals and stream water affected by acid mine drainage in a mineralized watershed (Colorado, USA). The δ 65 Cu values (based on 65 Cu/ 63 Cu) of enargite (δ 65 Cu = −0.01 ± 0.10‰; 2 σ ) and chalcopyrite (δ 65 Cu = 0.16 ± 0.10‰) are within the range of reported values for terrestrial primary Cu sulfides (−1‰ < δ 65 Cu < 1‰). These mineral samples show lower δ 65 Cu values than stream waters (1.38‰ ⩽ δ 65 Cu ⩽ 1.69‰). The average isotopic fractionation (Δ aq-min = δ 65 Cu aq − δ 65 Cu min , where the latter is measured on mineral samples from the field system), equals 1.43 ± 0.14‰ and 1.60 ± 0.14‰ for chalcopyrite and enargite, respectively. To interpret this field survey, we leached chalcopyrite and enargite in batch experiments and found that, as in the field, the leachate is enriched in 65 Cu relative to chalcopyrite (1.37 ± 0.14‰) and enargite (0.98 ± 0.14‰) when microorganisms are absent. Leaching of minerals in the presence of Acidithiobacillus ferrooxidans results in smaller average fractionation in the opposite direction for chalcopyrite ( &#x394; aq-min o = - 0.57 &#xB1; 0.14 &#x2030; ">‰ Δaq-mino=-0.57±0.14‰ , where min o refers to the starting mineral) and no apparent fractionation for enargite ( &#x394; aq-min o = 0.14 &#xB1; 0.14 &#x2030; ">‰ Δaq-mino=0.14±0.14‰ ). Abiotic fractionation is attributed to preferential oxidation of 65 Cu + at the interface of the isotopically homogeneous mineral and the surface oxidized layer, followed by solubilization. When microorganisms are present, the abiotic fractionation is most likely not seen due to preferential association of 65 Cu aq with A. ferrooxidans cells and related precipitates. In the biotic experiments, Cu was observed under TEM to occur in precipitates around bacteria and in intracellular polyphosphate granules. Thus, the values of δ 65 Cu in the field and laboratory systems are presumably determined by the balance of Cu released abiotically and Cu that interacts with cells and related precipitates. Such isotopic signatures resulting from Cu sulfide dissolution should be useful for acid mine drainage remediation and ore prospecting purposes.

Geochimica et Cosmochimica Acta↗
USGS2008

Eulerian-Lagrangian numerical scheme for simulating advection, dispersion, and transient storage in streams and a comparison of numerical methods

Past applications of one-dimensional advection, dispersion, and transient storage zone models have almost exclusively relied on a central differencing, Eulerian numerical approximation to the nonconservative form of the fundamental equation. However, there are scenarios where this approach generates unacceptable error. A new numerical scheme for this type of modeling is presented here that is based on tracking Lagrangian control volumes across a fixed (Eulerian) grid. Numerical tests are used to provide a direct comparison of the new scheme versus nonconservative Eulerian numerical methods, in terms of both accuracy and mass conservation. Key characteristics of systems for which the Lagrangian scheme performs better than the Eulerian scheme include: nonuniform flow fields, steep gradient plume fronts, and pulse and steady point source loadings in advection-dominated systems. A new analytical derivation is presented that provides insight into the loss of mass conservation in the nonconservative Eulerian scheme. This derivation shows that loss of mass conservation in the vicinity of spatial flow changes is directly proportional to the lateral inflow rate and the change in stream concentration due to the inflow. While the nonconservative Eulerian scheme has clearly worked well for past published applications, it is important for users to be aware of the scheme’s limitations.

Journal of Environmental Engineering↗
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↗
USGS2005

Sensitivity analysis of conservative and reactive stream transient storage models applied to field data from multiple-reach experiments

The transient storage model (TSM) has been widely used in studies of stream solute transport and fate, with an increasing emphasis on reactive solute transport. In this study we perform sensitivity analyses of a conservative TSM and two different reactive solute transport models (RSTM), one that includes first-order decay in the stream and the storage zone, and a second that considers sorption of a reactive solute on streambed sediments. Two previously analyzed data sets are examined with a focus on the reliability of these RSTMs in characterizing stream and storage zone solute reactions. Sensitivities of simulations to parameters within and among reaches, parameter coefficients of variation, and correlation coefficients are computed and analyzed. Our results indicate that (1) simulated values have the greatest sensitivity to parameters within the same reach, (2) simulated values are also sensitive to parameters in reaches immediately upstream and downstream (inter-reach sensitivity), (3) simulated values have decreasing sensitivity to parameters in reaches farther downstream, and (4) in-stream reactive solute data provide adequate data to resolve effective storage zone reaction parameters, given the model formulations. Simulations of reactive solutes are shown to be equally sensitive to transport parameters and effective reaction parameters of the model, evidence of the control of physical transport on reactive solute dynamics. Similar to conservative transport analysis, reactive solute simulations appear to be most sensitive to data collected during the rising and falling limb of the concentration breakthrough curve.

Advances in Water Resources↗
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

Inorganic N and P dynamics of Antarctic glacial meltwater streams as controlled by hyporheic exchange and benthic autotrophic communities

The McMurdo Dry Valleys of South Victoria Land, Antarctica, contain numerous glacial meltwater streams that drain into lakes on the valley floors. Many of the streams have abundant perennial mats of filamentous cyanobacteria. The algal mats grow during streamflow in the austral summer and are in a dormant freeze-dried state during the rest of the year. NO 3 and soluble reactive P (SRP) concentrations were lower in streams with abundant algal mats than in streams with sparse algal mats. NO 3 and SRP concentrations were higher in the hyporheic zone of a stream with abundant algal mats than in the stream itself. An experimental injection of LiCl, NaNO 3 , and K 3 PO 4 was conducted in Green Creek, which has abundant algal mats. Substantial hyporheic exchange occurred. The NO 3 and PO 4 concentrations at 50 m below the injection were 55 μM and 18 μM, respectively, during the experiment. NO 3 and PO 4 concentrations were below the detection limit of 1 to 2 μM at a site 497 m below the injection during the Cl tracer arrival, indicating a high capacity for nutrient uptake by algal communities. NO 2 and NH 4 were present at sites 226 and 327 m below the injection, indicating that, in addition to denitrification and algal uptake, dissimilatory NO 3 reduction to NO 2 and NH 4 may be a NO 3 sink during transport. Transport modelling with nutrient uptake represented as a 1 st -order process yielded reach-scale parameters of 4.3 × 10 −5 to 3.9 × 10 −4 /s and 1.4 × 10 −4 to 3.8 × 10 −4 /s for uptake of NO 3 and PO 4 , respectively. The best match with the observed data was a model in which PO 4 uptake occurred only in the main channel and NO 3 uptake occurred in the main channel and in the hyporheic zone. Hyporheic NO 3 uptake was 7 to 16% of the total uptake for the different stream reaches. These results demonstrate that nutrient flux to the lakes is controlled by hyporheic exchange and nutrient uptake by algal mats in dry valley streams. Streams without algal mats contribute more nutrients to the lakes than streams with algal mats.

Journal of the North American Benthological Societ↗
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

Evaluating remedial alternatives for an acid mine drainage stream: Application of a reactive transport model

A reactive transport model based on one-dimensional transport and equilibrium chemistry is applied to synoptic data from an acid mine drainage stream. Model inputs include streamflow estimates based on tracer dilution, inflow chemistry based on synoptic sampling, and equilibrium constants describing acid/base, complexation, precipitation/dissolution, and sorption reactions. The dominant features of observed spatial profiles in pH and metal concentration are reproduced along the 3.5-km study reach by simulating the precipitation of Fe(III) and Al solid phases and the sorption of Cu, As, and Pb onto freshly precipitated iron(III) oxides. Given this quantitative description of existing conditions, additional simulations are conducted to estimate the streamwater quality that could result from two hypothetical remediation plans. Both remediation plans involve the addition of CaCO 3 to raise the pH of a small, acidic inflow from ∼2.4 to ∼7.0. This pH increase results in a reduced metal load that is routed downstream by the reactive transport model, thereby providing an estimate of post-remediation water quality. The first remediation plan assumes a closed system wherein inflow Fe(II) is not oxidized by the treatment system; under the second remediation plan, an open system is assumed, and Fe(II) is oxidized within the treatment system. Both plans increase instream pH and substantially reduce total and dissolved concentrations of Al, As, Cu, and Fe(II+III) at the terminus of the study reach. Dissolved Pb concentrations are reduced by ∼18% under the first remediation plan due to sorption onto iron(III) oxides within the treatment system and stream channel. In contrast, iron(III) oxides are limiting under the second remediation plan, and removal of dissolved Pb occurs primarily within the treatment system. This limitation results in an increase in dissolved Pb concentrations over existing conditions as additional downstream sources of Pb are not attenuated by sorption.

Environmental Science & Technology↗
USGS2002

Assessment of metal loads in watersheds affected by acid mine drainage by using tracer injection and synoptic sampling: Cement Creek, Colorado, USA

Watersheds in mineralized zones may contain many mines, each of which can contribute to acidity and the metal load of a stream. In this study the authors delineate hydrogeologic characteristics determining the transport of metals from the watershed to the stream in the watershed of Cement Creek, Colorado. Combining the injection of a chemical tracer, to determine a discharge, with synoptic sampling, to obtain chemistry of major ions and metals, spatially detailed load profiles are quantified. Using the discharge and load profiles, the authors (1) identified sampled inflow sources which emanate from undisturbed as well as previously mined areas; (2) demonstrate, based on simple hydrologic balance, that unsampled, likely dispersed subsurface, inflows are significant; and (3) estimate attenuation. For example, along the 12-km study reach, 108 kg per day of Zn were added to Cement Creek. Almost half of this load came from 10 well-defined areas that included both mined and non-mined parts of the watershed. However, the combined effect of many smaller inflows also contributed a substantial load that could limit the effectiveness of remediation. Of the total Zn load, 58.3 kg/day came from stream segments with no visible inflow, indicating the importance of contributions from dispersed subsurface inflow. The subsurface inflow mostly occurred in areas with substantial fracturing of the bedrock or in areas downstream from tributaries with large alluvial fans. Despite a pH generally less than 4.5, there was 58.4 kg/day of Zn attenuation that occurred in mixing zones downstream from inflows with high pH. Mixing zones can have local areas of pH that are high enough for sorption and precipitation reactions to have an effect. Principal component analysis classified inflows into 7 groups with distinct chemical signatures that represent water-rock interaction with different mineral-alteration suites in the watershed. The present approach provides a detailed snapshot of metal load for the watershed to support remediation decisions, and quantifies processes affecting metal transport.

Colorado↗