Evaluating dual-domain models for upscaling multicomponent reactive transport in mine waste rock.

Reactive transport models have proven abilities to simulate the quantity and quality of drainage from mine waste rock. Tracer experiments indicate the presence of fast and slow flow regimes in many heterogeneous waste-rock piles. Although multidomain models have been developed specifically for systems with such distinctive hydrodynamics, there have been limited applications of multidomain reactive transport models to simulate composite drainage chemistries from waste-rock piles to date. This work evaluated the ability of dual-domain multicomponent reactive transport models (DDMRTMs) to reproduce breakthrough curves of conservative (chloride) and reactive (molybdenum) solutes observed at a well-characterized experimental waste-rock pile at the Antamina Mine, Peru. We found that the DDMRTM simulations quantitatively matched eight-year-long records of conservative transport through the waste-rock pile when parameterized mainly with field-measured properties obtained from the site and limited calibration. The DDMRTM model also provided a reasonable match to field observations of the reactive solute. The limited calibrated parameters are physically realistic, corroborating the ability of these multidomain models to reproduce the complex reactive-transport processes governing polluted rock drainage from large-scale waste-rock piles.

Reactive transport modelling to investigate multi-scale waste rock weathering processes.

Prediction of drainage quantity and quality is critical to reduce the environmental risks associated with weathering mine waste rock. Reactive transport models can be effective tools to understand and disentangle the processes underlying waste-rock weathering and drainage, but their validity and applicability can be impaired by poor parametrization and the non-uniqueness conundrum. Here, a process-based multicomponent reactive transport model is presented to interpret and quantify the processes affecting drainage quantity and quality from 15 waste- rock experiments from the Antamina mine, Peru. The deployed uniform flow formulation and consistent set of geochemical rate equations could be calibrated almost exclusively with measured bulk waste-rock properties in experiments ranging from 2 kg to 6500 tons in size. The quantitative agreement between simulated dynamics and the observed drainage records, for systems with a variety of rock lithologies and over a wide range of pH, supports the proposed selection of processes. The controls of important physicochemical processes and feedbacks such as secondary mineral precipitation, surface passivation, oxygen limitations, were confirmed through sensitivity analyses. Our work shows that reactive transport models with a consistent formulation and evidence-based parametrization can be used to explain waste-rock drainage dynamics across laboratory to field scales.

Identification, spatial extent and distribution of fugitive gas migration on the well pad scale.

Global methane (CH4) emissions are becoming increasingly important due to the contribution of CH4 to global warming. Leaking oil and gas wells can lead to subsurface CH4 gas migration (GM), which can cause both aquifer contamination and atmospheric emissions. Despite the need to identify and quantify GM at oil and gas well pads, effective and reliable monitoring techniques are lacking. In this field study, we used CH4 and carbon dioxide (CO2) efflux measurements together with soil gas stable carbon isotopic signatures to identify the occurrence and to characterize the spatio-temporal migration of fugitive gas across 17 selected well pads in Northeastern British Columbia, Canada. At 13 of these sites, operators had previously reported the occurrence of GM; however, subsequent inspections based on visual, olfactory or auditory evidence only identified GM at two of these sites. Using the soil gas efflux method, evidence for GM was found at 15 of the 17 well pads with CH4 and CO2 effluxes ranging from 0.017 to 180µmolm-2s-1(0.024 to 250gCH4m-2d-1) and 0.50 to 32µmolm-2s-1 (1.9 to 122gCO2m-2d-1), respectively. Stable carbon isotopic composition was assessed at 10 of the 17 well pads with 9 well pads showing evidence of GM. The isotopic values indicated that CH4 in soil gas was from the same origin as CH4 in the surface casing vent flow gas. There was no correlation between CH4 effluxes and the distance from the well head; an equal portion of elevated effluxes were detected >10m from the well head as were detected <5m from the well head. In addition, CH4 effluxes varied temporally with values changing by up to an order of magnitude over 2h. Although the study was carried out in Northeastern British Columbia, the results are applicable on a global scale, suggesting that inspections mostly based on visual evidence (e.g. bubbling at the well head) are not reliable for the identification of GM and, that infrequent survey measurements at predefined locations close to the well head may overestimate, underestimate or even miss CH4 effluxes. Repetitive and relatively densely spaced gas efflux measurements using a dynamic closed chamber method proved to be a useful tool for detecting GM.

Multi-year CO2 efflux measurements for assessing natural source zone depletion at a large hydrocarbon-impacted site.

The changing landscape of fuel consumption related, in part, to increased engine efficiency and the inexpensive supply of natural gas, has led to the closure of multiple refineries. As the operational lifetime of many refineries exceeds 100 years, historical releases of oil and refined products is common. To evaluate remediation and rehabilitation options, there is a need to understand the rate and distribution of natural hydrocarbon mass losses across these large properties. Here, surficial CO2 flux measurements were used to evaluate naturally occurring hydrocarbon mass losses at a large-scale former refinery that has been closed since 1982. Natural source zone depletion (NSZD) rates over a five-year period (2012-2016) were derived from surficial CO2 efflux measurements on a high-resolution grid (N > 80). Results demonstrate substantial variations of mass loss rates across the site. Average site-wide mass loss rates ranged from 1.1-5.4 g TPH m-2 d-1 as C10H22 with a multi-year average of 4.0 g TPH m-2 d-1 as decane (C10H22), consistent with observations at other sites. Statistical analysis demonstrated that the same average mass loss rates would have been obtained with fewer measurement locations (N = 20-30). Comparing NSZD rates to site metadata show CO2 fluxes to be a reasonably good proxy for zones of subsurface hydrocarbon contamination - particularly with respect to vadose zone impacts. It is hypothesized that the observed decline of NSZD rates over the study period is related to rise of groundwater levels, leading to increased submergence of the smear zone. Overall, mass loss rates calculated from CO2 fluxes show NSZD can result in substantial contaminant removal, which may rival that obtained from engineered remediation, under some conditions.

Evaluation of single- and dual-porosity models for reproducing the release of external and internal tracers from heterogeneous waste-rock piles.

Accurate predictions of solute release from waste-rock piles (WRPs) are paramount for decision making in mining-related environmental processes. Tracers provide information that can be used to estimate effective transport parameters and understand mechanisms controlling the hydraulic and geochemical behavior of WRPs. It is shown that internal tracers (i.e. initially present) together with external (i.e. applied) tracers provide complementary and quantitative information to identify transport mechanisms. The analysis focuses on two experimental WRPs, Piles 4 and Pile 5 at the Antamina Mine site (Peru), where both an internal chloride tracer and externally applied bromide tracer were monitored in discharge over three years. The results suggest that external tracers provide insight into transport associated with relatively fast flow regions that are activated during higher-rate recharge events. In contrast, internal tracers provide insight into mechanisms controlling solutes release from lower-permeability zones within the piles. Rate-limited diffusive processes, which can be mimicked by nonlocal mass-transfer models, affect both internal and external tracers. The sensitivity of the mass-transfer parameters to heterogeneity is higher for external tracers than for internal tracers, as indicated by the different mean residence times characterizing the flow paths associated with each tracer. The joint use of internal and external tracers provides a more comprehensive understanding of the transport mechanisms in WRPs. In particular, the tracer tests support the notion that a multi-porosity conceptualization of WRPs is more adequate for capturing key mechanisms than a dual-porosity conceptualization.

Molybdenum and zinc stable isotope variation in mining waste rock drainage and waste rock at the Antamina mine, Peru.

The stable isotope composition of molybdenum (Mo) and zinc (Zn) in mine wastes at the Antamina Copper-Zn-Mo mine, Peru, was characterized to investigate whether isotopic variation of these elements indicated metal attenuation processes in mine drainage. Waste rock and ore minerals were analyzed to identify the isotopic composition of Mo and Zn sources, namely molybdenites (MoS2) and sphalerites (ZnS). Molybdenum and Zn stable isotope ratios are reported relative to the NIST-SRM-3134 and PCIGR-1 Zn standards, respectively. δ(98)Mo among molybdenites ranged from -0.6 to +0.6‰ (n=9) while sphalerites showed no δ(66)Zn variations (0.11±0.01‰, 2 SD, n=5). Mine drainage samples from field waste rock weathering experiments were also analyzed to examine the extent of isotopic variability in the dissolved phase. Variations spanned 2.2‰ in δ(98)Mo (-0.1 to +2.1‰) and 0.7‰ in δ(66)Zn (-0.4 to +0.3‰) in mine drainage over a wide pH range (pH2.2-8.6). Lighter δ(66)Zn signatures were observed in alkaline pH conditions, which was consistent with Zn adsorption and/or hydrozincite (Zn5(OH)6(CO3)2) formation. However, in acidic mine drainage Zn isotopic compositions reflected the value of sphalerites. In addition, molybdenum isotope compositions in mine drainage were shifted towards heavier values (0.89±1.25‰, 2 SD, n=16), with some overlap, in comparison to molybdenites and waste rock (0.13±0.82‰, 2 SD, n=9). The cause of heavy Mo isotopic signatures in mine drainage was more difficult to resolve due to isotopic heterogeneity among ore minerals and a variety of possible overlapping processes including dissolution, adsorption and secondary mineral precipitation. This study shows that variation in metal isotope ratios are promising indicators of metal attenuation. Future characterization of isotopic fractionation associated to key environmental reactions will improve the power of Mo and Zn isotope ratios to track the fate of these elements in mine drainage.

Biogeochemical processes controlling the mobility of major ions and trace metals in aquitard sediments beneath an oil sand tailing pond: laboratory studies and reactive transport modeling.

Increased production and expansion of the oil sand industry in Alberta are of great benefit to the economy, but they carry major environmental challenges. The volume of fluid fine tailings requiring storage is 840×10(6) m(3) and growing, making it imperative that we better understand the fate and transport of oil sand process-affected water (OSPW) seepage from these facilities. Accordingly, the current study seeks to characterize both a) the potential for major ion and trace element release, and b) the principal biogeochemical processes involved, as tailing pond OSPW infiltrates into, and interacts with, underlying glacial till sediments prior to reaching down gradient aquifers or surface waters. Objectives were addressed through a series of aqueous and solid phase experiments, including radial diffusion cells, an isotope analysis, X-ray diffraction, and sequential extractions. The diffusion cells were also simulated in a reactive transport framework to elucidate key reaction processes. The experiments indicate that the ingress and interaction of OSPW with the glacial till sediment-pore water system will result in: a mitigation of ingressing Na (retardation), displacement and then limited precipitation of exchangeable Ca and Mg (as carbonates), sulfate reduction and subsequent precipitation of the produced sulfides, as well as biodegradation of organic carbon. High concentrations of ingressing Cl (~375 mg L(-1)) and Na (~575 mg L(-1)) (even though the latter is delayed, or retarded) are expected to migrate through the till and into the underlying sand channel. Trace element mobility was influenced by ion exchange, oxidation-reduction, and mineral phase reactions including reductive dissolution of metal oxyhydroxides - in accordance with previous observations within sandy aquifer settings. Furthermore, although several trace elements showed the potential for release (Al, B, Ba, Cd, Mn, Pb, Si, Sr), large-scale mobilization is not supported. Thus, the present results suggest that in addition to the commonly cited naphthenic acids, remediation of OSPW-impacted groundwater will need to address high concentrations of major ions contributing to salinization.

Vadose zone attenuation of organic compounds at a crude oil spill site - interactions between biogeochemical reactions and multicomponent gas transport.

Contaminant attenuation processes in the vadose zone of a crude oil spill site near Bemidji, MN have been simulated with a reactive transport model that includes multicomponent gas transport, solute transport, and the most relevant biogeochemical reactions. Dissolution and volatilization of oil components, their aerobic and anaerobic degradation coupled with sequential electron acceptor consumption, ingress of atmospheric O(2), and the release of CH(4) and CO(2) from the smear zone generated by the floating oil were considered. The focus of the simulations was to assess the dynamics between biodegradation and gas transport processes in the vadose zone, to evaluate the rates and contributions of different electron accepting processes towards vadose zone natural attenuation, and to provide an estimate of the historical mass loss. Concentration distributions of reactive (O(2), CH(4), and CO(2)) and non-reactive (Ar and N(2)) gases served as key constraints for the model calibration. Simulation results confirm that as of 2007, the main degradation pathway can be attributed to methanogenic degradation of organic compounds in the smear zone and the vadose zone resulting in a contaminant plume dominated by high CH(4) concentrations. In accordance with field observations, zones of volatilization and CH(4) generation are correlated to slightly elevated total gas pressures and low partial pressures of N(2) and Ar, while zones of aerobic CH(4) oxidation are characterized by slightly reduced gas pressures and elevated concentrations of N(2) and Ar. Diffusion is the most significant transport mechanism for gases in the vadose zone; however, the simulations also indicate that, despite very small pressure gradients, advection contributes up to 15% towards the net flux of CH(4), and to a more limited extent to O(2) ingress. Model calibration strongly suggests that transfer of biogenically generated gases from the smear zone provides a major control on vadose zone gas distributions and vadose zone carbon balance. Overall, the model was successful in capturing the complex interactions between biogeochemical reactions and multicomponent gas transport processes. However, despite employing a process-based modeling approach, honoring observed parameter ranges, and generally obtaining good agreement between field observations and model simulations, accurate quantification of natural attenuation rates remains difficult. The modeling results are affected by uncertainties regarding gas phase saturations, tortuosities, and the magnitude of CH(4) and CO(2) flux from the smear zone. These findings highlight the need to better delineate gas fluxes at the model boundaries, which will help constrain contaminant degradation rates, and ultimately source zone longevity.

The importance of conceptual models in the reactive transport simulation of oxygen ingress in sparsely fractured crystalline rock.

Redox evolution in sparsely fractured crystalline rocks is a key, and largely unresolved, issue when assessing the geochemical suitability of deep geological repositories for nuclear waste. Redox zonation created by the influx of oxygenated waters has previously been simulated using reactive transport models that have incorporated a variety of processes, resulting in predictions for the depth of oxygen penetration that may vary greatly. An assessment and direct comparison of the various underlying conceptual models are therefore needed. In this work a reactive transport model that considers multiple processes in an integrated manner is used to investigate the ingress of oxygen for both single fracture and fracture zone scenarios. It is shown that the depth of dissolved oxygen migration is greatly influenced by the a priori assumptions that are made in the conceptual models. For example, the ability of oxygen to access and react with minerals in the rock matrix may be of paramount importance for single fracture conceptual models. For fracture zone systems, the abundance and reactivity of minerals within the fractures and thin matrix slabs between the fractures appear to provide key controls on O(2) attenuation. The findings point to the need for improved understanding of the coupling between the key transport-reaction feedbacks to determine which conceptual models are most suitable and to provide guidance for which parameters should be targeted in field and laboratory investigations.

A detailed field-based evaluation of naphthenic acid mobility in groundwater.

An anaerobic plume of process-affected groundwater was characterized in a shallow sand aquifer adjacent to an oil sands tailings impoundment. Based on biological oxygen demand measurements, the reductive capacity of the plume is considered minimal. Major dissolved components associated with the plume include HCO(3), Na, Cl, SO(4), and naphthenic acids (NAs). Quantitative and qualitative NA analyses were performed on groundwater samples to investigate NA fate and transport in the subsurface. Despite subsurface residence times exceeding 20 years, significant attenuation of NAs by biodegradation was not observed based on screening techniques developed at the time of the investigation. Relative to conservative tracers (i.e., Cl), overall NA attenuation in the subsurface is limited, which is consistent with batch sorption and microcosm studies performed by other authors. Insignificant biological oxygen demand and low concentrations of dissolved As (<10 microg L(-1)) in the plume suggest that the potential for secondary trace metal release, specifically As, via reductive dissolution reactions driven by ingress of process-affected water is minimal. It is also possible that readily leachable As is not present in significant quantities within the sediments of the study area. Thus, for similar plumes of process-affected groundwater in shallow sand aquifers which may occur as oil sands mining expands, a reasonable expectation is for NA persistence, but minimal trace metal mobilization.

Transport and reaction processes affecting the attenuation of landfill gas in cover soils.

Methane and trace organic gases produced in landfill waste are partly oxidized in the top 40 cm of landfill cover soils under aerobic conditions. The balance between the oxidation of landfill gases and the ingress of atmospheric oxygen into the soil cover determines the attenuation of emissions of methane, chlorofluorocarbons, and hydrochlorofluorocarbons to the atmosphere. This study was conducted to investigate the effect of oxidation reactions on the overall gas transport regime and to evaluate the contributions of various gas transport processes on methane attenuation in landfill cover soils. For this purpose, a reactive transport model that includes advection and the Dusty Gas Model for simulation of multicomponent gas diffusion was used. The simulations are constrained by data from a series of counter-gradient laboratory experiments. Diffusion typically accounts for over 99% of methane emission to the atmosphere. Oxygen supply into the soil column is driven exclusively by diffusion, whereas advection outward offsets part of the diffusive contribution. In the reaction zone, methane consumption reduces the pressure gradient, further decreasing the significance of advection near the top of the column. Simulations suggest that production of water or accumulation of exopolymeric substances due to microbially mediated methane oxidation can significantly reduce diffusive fluxes. Assuming a constant rate of methane production within a landfill, reduction of the diffusive transport properties, primarily due to exopolymeric substance production, may result in reduced methane attenuation due to limited O(2) -ingress.

Identification of key parameters controlling dissolved oxygen migration and attenuation in fractured crystalline rocks.

In the crystalline rocks of the Canadian Shield, geochemical conditions are currently reducing at depths of 500-1000 m. However, during future glacial periods, altered hydrologic conditions could potentially result in enhanced recharge of glacial melt water containing a relatively high concentration of dissolved oxygen (O2). It is therefore of interest to investigate the physical and geochemical processes, including naturally-occurring redox reactions, that may control O2 ingress. In this study, the reactive transport code MIN3P is used in combination with 2k factorial analyses to identify the most important parameters controlling oxygen migration and attenuation in fractured crystalline rocks. Scenarios considered are based on simplified conceptual models that include a single vertical fracture, or a fracture zone, contained within a rock matrix that extends from the ground surface to a depth of 500 m. Consistent with field observations, Fe(II)-bearing minerals are present in the fractures (i.e. chlorite) and the rock matrix (biotite and small quantities of pyrite). For the parameter ranges investigated, results indicate that for the single fracture case, the most influential factors controlling dissolved O2 ingress are flow velocity in the fracture, fracture aperture, and the biotite reaction rate in the rock matrix. The most important parameters for the fracture zone simulations are flow velocity in the individual fractures, pO2 in the recharge water, biotite reaction rate, and to a lesser degree the abundance and reactivity of chlorite in the fracture zone, and the fracture zone width. These parameters should therefore receive increased consideration during site characterization, and in the formulation of site-specific models intended to predict O2 behavior in crystalline rocks.

Integration of field measurements and reactive transport modelling to evaluate contaminant transport at a sulfide mine tailings impoundment.

Over a decade of field observations including geochemical, mineralogical and hydrological information are available on the generation of acid mine drainage from the Pistol Dam region of the P-area of Inco's tailings impoundment in Copper Cliff, Ontario. This work focuses on the integration and quantitative assessment of this data set using reactive transport modeling. The results of the reactive transport simulations are in general agreement with the field observations; however, exact agreement between the field and simulated results was not the objective of this study, and was not attained. Many factors contribute to the discrepancies between the field observations and simulation results including geochemical and hydrogeological complexities and necessary model simplifications. For example, fluctuating water levels observed at the site were averaged and described using a steady state flow system. In addition, the lack of representative thermodynamic and rate expression data contributed to the discrepancies between observations and simulation results, thus further research into the applicability of laboratory-derived thermodynamic and rate expression data to field conditions could minimize these discrepancies. Despite the discrepancies between the field observations and simulated results, integrating field observations with numerical modelling of the P-area tailings impoundment allowed for a more complete understanding of what affects the complex geochemical reactions.

Process-based reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage.

Reactive transport modeling of a permeable reactive barrier for the treatment of mine drainage was used to integrate a comprehensive data set including pore water chemistry and solid phase data from several sampling events over a >3-year time period. The simulations consider the reduction of sulfate by the organic carbon-based treatment material and the removal of sulfate and iron by precipitation of reduced mineral phases including iron monosulfides and siderite. Additional parameters constraining the model include dissolved H2S, alkalinity and pH, as well as a suite of solid phase S-fractions identified by extractions. Influences of spatial heterogeneity necessitated the use of a 2-dimensional modeling approach. Simulating observed seasonal fluctuations and long-term changes in barrier reactivity required the use of temperature dependent rate coefficients and a multimodal Monod-type rate expression accounting for the variable reactivity of different organic carbon fractions. Simulated dissolved concentrations of SO4, Fe, H2S, alkalinity and pH, as well as solid phase accumulations of reduced sulfur phases generally compare well to observed trends over 23 months. Spatial variations, seasonal fluctuations and the time-dependent decline in reactivity were also captured. The modeling results generally confirm, and further strengthen, the existing conceptual model for the site. Overall sulfate reduction and S-accumulation rates are constrained with confidence within a factor of 1.5.

Modelling the closure-related geochemical evolution of groundwater at a former uranium mine.

A newly developed reactive transport model was used to evaluate the potential effects of mine closure on the geochemical evolution in the aquifer downgradient from a mine site. The simulations were conducted for the Königstein uranium mine located in Saxony, Germany. During decades of operation, uranium at the former mine site had been extracted by in situ acid leaching of the ore underground, while the mine was maintained in a dewatered condition. One option for decommissioning is to allow the groundwater level to rise to its natural level, flooding the mine workings. As a result, pore water containing high concentrations of dissolved metals, radionuclides, and sulfate may be released. Additional contamination may arise due to the dissolution of minerals contained in the aquifer downgradient of the mine. On the other hand, dissolved metals may be attenuated by reactions within the aquifer. The geochemical processes and interactions involved are highly non-linear and their impact on the quality of the groundwater and surface water downstream of the mine is not always intuitive. The multicomponent reactive transport model MIN3P, which can describe mineral dissolution-precipitation reactions, aqueous complexation, and oxidation-reduction reactions, is shown to be a powerful tool for investigating these processes. The predictive capabilities of the model are, however, limited by the availability of key geochemical parameters such as the presence and quantities of primary and secondary mineral phases. Under these conditions, the model can provide valuable insight by means of sensitivity analyses.

Reactive transport modeling of processes controlling the distribution and natural attenuation of phenolic compounds in a deep sandstone aquifer.

Reactive solute transport modeling was utilized to evaluate the potential for natural attenuation of a contaminant plume containing phenolic compounds at a chemical producer in the West Midlands, UK. The reactive transport simulations consider microbially mediated biodegradation of the phenolic compounds (phenols, cresols, and xylenols) by multiple electron acceptors. Inorganic reactions including hydrolysis, aqueous complexation, dissolution of primary minerals, formation of secondary mineral phases, and ion exchange are considered. One-dimensional (1D) and three-dimensional (3D) simulations were conducted. Mass balance calculations indicate that biodegradation in the saturated zone has degraded approximately 1-5% of the organic contaminant plume over a time period of 47 years. Simulations indicate that denitrification is the most significant degradation process, accounting for approximately 50% of the organic contaminant removal, followed by sulfate reduction and fermentation reactions, each contributing 15-20%. Aerobic respiration accounts for less than 10% of the observed contaminant removal in the saturated zone. Although concentrations of Fe(III) and Mn(IV) mineral phases are high in the aquifer sediment, reductive dissolution is limited, producing only 5% of the observed mass loss. Mass balance calculations suggest that no more than 20-25% of the observed total inorganic carbon (TIC) was generated from biodegradation reactions in the saturated zone. Simulations indicate that aerobic biodegradation in the unsaturated zone, before the contaminant entered the aquifer, may have produced the majority of the TIC observed in the plume. Because long-term degradation is limited to processes within the saturated zone, use of observed TIC concentrations to predict the future natural attenuation may overestimate contaminant degradation by a factor of 4-5.