Diavik Waste Rock Project: Evolution of Mineral Weathering, Element Release, and Acid Generation and Neutralization during a Five-Year Humidity Cell Experiment

Diavik Waste Rock Project: Simulation of the geochemical evolution of a large test pile using a scaled temperature and sulfide-content dependent reactive transport model

The Diavik Waste Rock Project (DWRP) project included four principal components focused on the development of techniques for assessing the environmental impacts of waste rock at mine sites. These components were small-volume laboratory experiments, intermediate- and large-volume field experiments, and assessment of the operational-scale waste-rock stockpiles, which facilitated characterization of waste-rock weathering at different scales. The heavily instrumented large-scale field experiments (test piles) were constructed to replicate, as closely as practicable, the temperature, water flow, and gas transport regimes of a waste-rock pile that is exposed to annual freezing and thawing cycles and to facilitate characterization of the long-term weathering of a low-sulfide waste rock. An integrated conceptual model of sulfide-bearing waste-rock weathering, developed at the small scale, was applied to assess the capacity of the conceptual model to capture the geochemical evolution of the waste rock at the large field-scale test-pile experiment. The integrated conceptual model was implemented using reactive transport code MIN3P, taking into account scale-dependent mechanisms. The test-pile mineralogy was similar to the small-scale laboratory experiments and included low-sulfide waste rock with an S content of 0.053 wt% (primarily pyrrhotite). The flow regime of the test pile was simulated using parameters measured as part of other DWRP investigations, including temporally variable infiltration estimates that represented the measured precipitation events at the site. The temporally and spatially variable temperature of the test pile was interpolated from values measured using instrumentation installed at the beginning of the experiment and was included in the simulation to refine the temperature dependence of the geochemical reactions. To allow continuous, multi-year simulation, freezing was also simulated to represent the conditions experienced at the test-pile experiment. Normalized root mean square error analysis of the large-scale field experiment simulation results indicated most parameters compare well to measured daily mass flux (i.e., the fraction of the range of annual values encompassed in the residual was less than 0.5 for SO₄, Fe, Ni, Si, Ca, K, Mg, Na, and pH and 1.0 or less for all parameters except Cu). The method of using an integrated conceptual model developed from the results of humidity cell experiments to implement a mechanistic approach for assessing the primary geochemical processes of waste-rock weathering on a large scale was shown to provide reasonable results; however, the results are specific to the study site and the approach requires application to various sites under different geological and climatological conditions to facilitate further refinement.

Travel time-based modelling of nitrate reduction in a fractured limestone aquifer by pyrite and iron carbonates under pore size limitation

We investigate denitrification in a ferric iron-containing fractured micritic limestone aquifer (Triassic Upper Muschelkalk) in south-west Germany by numerical simulations. Low porosity values (average value of 1%), partly small pore sizes of the rock matrix (~ 0.1 μm), and thus potential absence of microbial activity in the rock matrix suggest that denitrification is taking place solely in the fracture. A key question is whether the nitrate reduction derived from groundwater observations at 25 locations in the study area can be explained by a model that restricts microbial denitrification to the fractures. A travel time-based reactive transport model is developed to efficiently simulate long-term nitrate reduction on the catchment scale. The model employs a 2-D numerical reaction model describing the fracture-rock matrix system and parametric travel time distributions. The role of (i) biotic and abiotic iron oxidation, (ii) the type and amount of iron bearing minerals, and (iii) mass transfer between matrix and fracture are investigated. The simulations show that pyrite and siderite (used as surrogate for iron carbonates) together as a source of electron donors provide enough reduction potential to decrease the nitrate concentrations as observed in the field. This confirms the hypothesis that diffusion-controlled mass transfer of electron donors from the matrix to the fracture is sufficient to establish considerable denitrification in the fracture. Uncertainty in modelled concentrations is demonstrated as a result of both the geochemical aquifer properties and the unknown shape of travel time distributions.

Diavik Waste Rock Project: Geostatistical Analysis of Sulfur, Carbon, and Hydraulic Conductivity Distribution in a Large-Scale Experimental Waste Rock Pile

One of the large-scale field waste rock experiments (test piles) conducted as part of the Diavik Waste Rock Project was deconstructed, providing a spatially located set of geochemical, mineralogical, and particle-size distribution samples. Geostatistical analyses were conducted for sulfur and carbon content and saturated hydraulic conductivity, which affect the geochemical evolution of waste rock, to investigate the spatial dependence of these parameters. Analyses included population statistics, experimental semi-variogram estimation, and theoretical semi-variogram fitting. Population statistics were calculated for additional data sets from samples collected during the construction of the test piles. The population statistical analyses indicated that log-normal distribution provided the best fit for all investigated data sets. Experimental semi-variograms were estimated for the spatially located data set (test pile deconstruction) using the classical estimator, and theoretical semi-variograms were fitted. This investigation showed that the spatial distribution of sulfur, carbon, and hydraulic conductivity within the core of the test-pile experiments can be approximated using a log-normal distribution with a mean and standard deviation calculated using the samples collected during construction of the piles, and that little to no spatial relationship was present for these parameters at the scale of sampling. That the saturated hydraulic conductivity of the matrix material can be represented by the same statistical distribution throughout the test pile is significant because water flow, as well as mineral surface area and reactivity are dominantly controlled within the matrix portion of the test pile. Reactive transport simulations are included to demonstrate the influence of the matrix material on effluent geochemistry.

Experimental Investigation on the Energy Properties and Failure Process of Thermal Shock Treated Sandstone Subjected to Coupled Dynamic and Static Loads

Thermal shock (TS) is known as the process where fractures are generated when rocks go through sudden temperature changes. In the field of deep rock engineering, the rock mass can be subjected to the TS process in various circumstances. To study the influence of TS on the mechanical behaviors of rock, sandstone specimens are heated at different high temperatures and three cooling methods (stove cooling, air cooling, and freezer cooling) are adopted to provide different cooling rates. The coupled dynamic and static loading tests are performed on the heated sandstone through a modified split Hopkinson pressure bar (SHPB) system. The influence of heating level and cooling rate on the dynamic compressive strength, energy dissipations, and fracturing characteristics is investigated based on the experimental data. The development of the microcracks of the sandstone specimens after the experiment is analyzed utilizing a scanning electron microscope (SEM). The extent of the development of the microcracks serves to explain the variation pattern of the mechanical responses and energy dissipations of the specimens obtained from the loading test. The findings of this study are valuable for practices in rock engineering involving high temperature and fast cooling.

Dynamic Properties of Thermal Shock Treated Sandstone Subjected to Coupled Dynamic and Static Loads

In deep rock engineering, the rock mass can be subjected to thermal stress caused by sudden changes in temperature, which is referred to as thermal shock (TS). To study the effect of TS on heated sandstone, three cooling methods are used to provide different cooling rates. Then the coupled dynamic and static loading tests are carried out on the heated sandstone by means of a modified split Hopkinson pressure bar (SHPB) system. The test results show that as the heating level increases, the dry density, P-wave velocity, and the dynamic combined strength of the heated sandstone decrease, while specimen porosity increases. Particularly, a sharp change in the physical properties of sandstone can be observed at 650 °C, which is believed to be caused by the α-β transition of quartz at 573 °C. At each heating level of the test, the damage caused by the higher cooling rate to the heated sandstone is more than that caused by the lower cooling rate. The different failure modes of sandstone with increasing temperature are analyzed. The mechanism of TS acting on heated sandstone is discussed, and two typical fracture patterns reflecting the action of TS are identified through SEM.

Copper release from waste rocks in an abandoned mine (NE, Brazil) and its impacts on ecosystem environmental quality

This study aimed to estimate the impact of an abandoned copper (Cu) mine on ecosystem environmental quality, using integrated ecological and biogeochemical analyses. Through a controlled experiment, the amount of Cu released by waste rocks, Cu adsorbed in soils, Cu geochemical behaviour and its leached amount were measured. Furthermore, to investigate the impacts of mine drainage on the adjacent ecosystem, samples of sediments, water and aquatic macroinvertebrates were analysed. We found that waste rocks still have high Cu concentration even after 30 years under weathering, ranging from 7782 to 8717 mg kg^-1, associated mainly with carbonates, amorphous oxides and sulphides. It was estimated that 7.2 tonnes of Cu were released by waste rocks into the environment over last 30 years. The concentration of Cu observed in Ubari stream water was (<dl to 90 μg L^-1), in sediments (28.0-1185 mg kg^-1) and in macroinvertebrates (1.3-28.9 mg kg^-1 d/w). The ecological indexes showed that near mine discharge a significance decrease in the density of aquatic macroinvertebrates and a significance increase of Cu in biological tissues occurs, causing disturbances in biodiversity. The results showed that, even after long periods, the waste rocks from abandoned mines still contain high levels of metal, that are gradually released into the environment through weathering and erosion, representing a potential source of environmental pollution and a clear threat to the environmental quality of adjacent ecosystems.

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.

Mine Waste Rock: Insights for Sustainable Hydrogeochemical Management

Mismanagement of mine waste rock can mobilize acidity, metal (loid)s, and other contaminants, and thereby negatively affect downstream environments. Hence, strategic long-term planning is required to prevent and mitigate deleterious environmental impacts. Technical frameworks to support waste-rock management have existed for decades and typically combine static and kinetic testing, field-scale experiments, and sometimes reactive-transport models. Yet, the design and implementation of robust long-term solutions remains challenging to date, due to site-specificity in the generated waste rock and local weathering conditions, physicochemical heterogeneity in large-scale systems, and the intricate coupling between chemical kinetics and mass- and heat-transfer processes. This work reviews recent advances in our understanding of the hydrogeochemical behavior of mine waste rock, including improved laboratory testing procedures, innovative analytical techniques, multi-scale field investigations, and reactive-transport modeling. Remaining knowledge-gaps pertaining to the processes involved in mine waste weathering and their parameterization are identified. Practical and sustainable waste-rock management decisions can to a large extent be informed by evidence-based simplification of complex waste-rock systems and through targeted quantification of a limited number of physicochemical parameters. Future research on the key (bio)geochemical processes and transport dynamics in waste-rock piles is essential to further optimize management and minimize potential negative environmental impacts.

Long-term monitoring of waste-rock weathering at the Antamina mine, Peru

The weathering of mine waste rock can cause release of metal-laden and acidic drainage that requires long-term and costly environmental management. To identify and quantify the geochemical processes and physical transport mechanisms controlling drainage quality, we monitored the weathering of five large-scale (20,000 t) instrumented waste-rock piles of variable and mixed-composition at the Antamina mine, Peru, in a decade-long monitoring program. Fine-grained, sulfidic waste rock with low-carbonate content exhibited high sulfide oxidation rates (>1 g S kg^-1 waste rock yr^-1) and within 7 years produced acidic (pH < 3) drainage with high Cu and Zn concentrations in the g L^-1 range. In contrast, drainage from coarse, carbonate-rich waste rock remained neutral for >10 years and had significantly lower metal loads. Efficient metal retention (>99%) caused by sorption and secondary mineral formation of e.g., gypsum, Fe-(oxy)hydroxides, and Cu/Zn-hydroxysulfates enforced strong (temporary) controls on drainage quality. Furthermore, reactive waste-rock fractions, as small as 10% of total mass, dominated the overall drainage chemistry from the waste-rock piles through internal mixing. This study demonstrates that a reliable prediction of the timing and quality of waste-rock drainage on practice-relevant spatiotemporal scales requires a quantitative understanding of the prevailing in-situ porewater conditions, secondary mineralogy, and spatial distribution of reactive waste-rock fractions in composite piles.

Forecasting Geoenvironmental Risks: Integrated Applications of Mineralogical and Chemical Data

Management of solid mine wastes requires detailed material characterisation at the start of a project to minimize opportunities for the generation of acid and metalliferous drainage (AMD). Mine planning must focus on obtaining a thorough understanding of the environmental properties of the future waste rock materials. Using drill core obtained from a porphyry Cu project in Northern Europe, this study demonstrates the integrated application of mineralogical and geochemical data to enable the construction of enviro-geometallurgical models. Geoenvironmental core logging, static chemical testing, bulk- and hyperspectral mineralogical techniques, and calculated mineralogy from assay techniques were used to critically evaluate the potential for AMD formation. These techniques provide value-adding opportunities to existing datasets and provide robust cross-validation methods for each technique. A new geoenvironmental logging code and a new geoenvironmental index using hyperspectral mineralogical data (Hy-GI) were developed and embedded into the geochemistry-mineralogy-texture-geometallurgy (GMTG) approach for waste characterisation. This approach is recommended for new mining projects (i.e., early life-of-mine stages) to ensure accurate geoenvironmental forecasting, therefore facilitating the development of an effective waste management plan that minimizes geoenvironmental risks posed by the mined materials.

Novel Microbial Assemblages Dominate Weathered Sulfide-Bearing Rock from Copper-Nickel Deposits in the Duluth Complex, Minnesota, USA

The Duluth Complex in northeastern Minnesota hosts economically significant deposits of copper, nickel, and platinum group elements (PGEs). The primary sulfide mineralogy of these deposits includes the minerals pyrrhotite, chalcopyrite, pentlandite, and cubanite, and weathering experiments show that most sulfide-bearing rock from the Duluth Complex generates moderately acidic leachate (pH 4 to 6). Microorganisms are important catalysts for metal sulfide oxidation and could influence the quality of water from mines in the Duluth Complex. Nevertheless, compared with that of extremely acidic environments, much less is known about the microbial ecology of moderately acidic sulfide-bearing mine waste, and so existing information may have little relevance to those microorganisms catalyzing oxidation reactions in the Duluth Complex. Here, we characterized the microbial communities in decade-long weathering experiments (kinetic tests) conducted on crushed rock and tailings from the Duluth Complex. Analyses of 16S rRNA genes and transcripts showed that differences among microbial communities correspond to pH, rock type, and experimental treatment. Moreover, microbial communities from the weathered Duluth Complex rock were dominated by taxa that are not typically associated with acidic mine waste. The most abundant operational taxonomic units (OTUs) were from the genera Meiothermus and Sulfuriferula, as well as from diverse clades of uncultivated Chloroflexi, Acidobacteria, and Betaproteobacteria Specific taxa, including putative sulfur-oxidizing Sulfuriferula spp., appeared to be primarily associated with Duluth Complex rock, but not pyrite-bearing rocks subjected to the same experimental treatment. We discuss the implications of these results for the microbial ecology of moderately acidic mine waste with low sulfide content, as well as for kinetic testing of mine waste.IMPORTANCE Economic sulfide mineral deposits in the Duluth Complex may represent the largest undeveloped source of copper and nickel on Earth. Microorganisms are important catalysts for sulfide mineral oxidation, and research on extreme acidophiles has improved our ability to manage and remediate mine wastes. We found that the microbial assemblages associated with weathered rock from the Duluth Complex are dominated by organisms not widely associated with mine waste or mining-impacted environments, and we describe geochemical and experimental influences on community composition. This report will be a useful foundation for understanding the microbial biogeochemistry of moderately acidic mine waste from these and similar deposits.

Geochemical and mineralogical characterization of sulfur and iron in coal waste rock, Elk Valley, British Columbia, Canada

Exposure of coal waste rock to atmospheric oxygen can result in the oxidation of sulfide minerals and the release of sulfate (SO₄^2-) and associated trace elements (e.g., Se, As, Cd, and Zn) to groundwaters and surface waters. Similarly, reduced iron minerals such as siderite, ankerite, and the sulfide, pyrite, present in the waste rock can also undergo oxidation, resulting in the formation of iron oxyhydroxides that can adsorb trace elements released from the oxidation of the sulfide minerals. Characterization and quantification of the distribution of sulfide and iron minerals, their oxidation products, as well as leaching rates are critical to assessing present-day and future impacts of SO₄^2- and associated trace elements on receiving waters. Synchrotron-based X-ray absorption near edge spectroscopic analysis of coal waste rock samples from the Elk Valley, British Columbia showed Fe present as pyrite (mean 6.0%), siderite (mean 44.3%), goethite (mean 35.4%), and lepidocrocite (mean 14.3%) with S present as sulfide (mean 26.9%), organic S (mean 58.7%), and SO₄^2- (mean 14.4%). Squeezed porewater samples from dump solids yielded mean concentrations of 0.28mg/L Fe and 1246mg/L SO₄^2-. Geochemical modeling showed the porewaters in the dumps to be supersaturated with respect to Fe oxyhydroxides and undersaturated with respect to gypsum, consistent with solids analyses. Coupling Fe and S mineralogical data with long-term water quality and quantity measurements from the base of one dump suggest about 10% of the sulfides (which represent 2% of total S) in the dump were oxidized over the past 30years. The S from these oxidized sulfides was released to the receiving surface water as SO₄^2- and the majority of the Fe precipitated as secondary Fe oxyhydroxides (only 3.0×10^-5% of the Fe was released to the receiving waters over the past 30years). Although the data suggest that the leaching of SO₄^2- from the waste rock dump could continue for about 300years, assuming no change in the rate of oxidation of sulfides, SO₄^2- is currently not a concern in receiving surface waters as the concentration levels are below regulatory limits.