Geochemical transport modelling of drainage from experimental mine tailings cells covered by capillary barriers

The use of covers with capillary barrier effects (CCBEs) for reducing acid mine drainage (AMD) from sulphidic mine tailings is simulated using the MIN3P finite volume model for coupled groundwater flow, O₂ diffusion and multi-component reactive transport. The model is applied to simulate five pilot-scale in situ test cells containing reactive tailings from the Manitou mine site, Val d’Or, Que., Canada. Four of the cells were constructed with CCBEs over the tailings, while the fifth tailings cell was left uncovered. Observed and simulated discharge from the base of each cell showed that the capillary barrier covers significantly reduced sulphide oxidation and AMD. Compared to acidic discharge from the uncovered cell, discharge from the four CCBE-covered cells had neutral pH levels and 1–7 orders of magnitude lower concentrations of SO₄, Fe, Zn, Cu and Al. The simulations showed that the moisture retaining layer of the CCBEs reduced AMD by inhibiting O₂ diffusion into the underlying reactive wastes. Provided the moisture-retention layer of the CCBE remains close to saturation, its thickness had a relatively minor effect. Under such near-saturated conditions, O₂ availability is limited by its diffusion rate through the bulk porous medium and not by the diffusion rate through the oxidized grain shells. The model is providing important new insights for comparing design alternatives for reducing or controlling AMD.     

The presence of anomalous trace element levels in present day Jamaican soils and the geochemistry of Late-Miocene or Pliocene phosphorites

High levels of Cd and Zn in Jamaican soils observed in geochemical surveys are related to the presence of phosphorites of possible Late-Miocene or Pliocene age. The trace element and REE geochemistry of the phosphorites, together with SEM studies, indicate a guano origin for the phosphorites. No specific host minerals for Cd could be identified in the fossiliferous phosphorite which is characterized by uniquely high levels of Cd, Zn, Ag, Be, U and Y. However, in the soil Cd is present in lithiophorite and a complex history of pedological development is preserved in the aluminous–goethite present in the soil. The unique guano signature is preserved in the soil despite the fact that guanos themselves have either not been observed or have been destroyed by continuing karst and soil development. The phosphorite geochemical signature can be traced in the data of a 1988 island-wide soil geochemical survey, identifying areas where the Palaeo-environment that supported bird ‘rookeries’ existed in the Late-Miocene or Pliocene.     

Geochemical characterization of surface waters and groundwater resources in the Managua area (Nicaragua,Central America)

This paper reports new geochemical data on dissolved major and minor constituents in surface waters and ground waters collected in the Managua region (Nicaragua), and provides a preliminary characterization of the hydrogeochemical processes governing the natural water evolution in this area. The peculiar geological features of the study site, an active tectonic region (Nicaragua Depression) characterized by active volcanism and thermalism, combined with significant anthropogenic pressure, contribute to a complex evolution of water chemistry, which results from the simultaneous action of several geochemical processes such as evaporation, rock leaching, mixing with saline brines of natural or anthropogenic origin. The effect of active thermalism on both surface waters (e.g., saline volcanic lakes) and groundwaters, as a result of mixing with variable proportions of hyper-saline geothermal Na–Cl brines (e.g., Momotombo geothermal plant), accounts for the high salinities and high concentrations of many environmentally-relevant trace elements (As, B, Fe and Mn) in the waters. At the same time the active extensional tectonics of the Managua area favour the interaction with acidic, reduced thermal fluids, followed by extensive leaching of the host rock and the groundwater release of toxic metals (e.g., Ni, Cu). The significant pollution in the area, deriving principally from urban and industrial waste-water, probably also contributes to the aquatic cycling of many trace elements, which attain concentrations above the WHO recommended limits for the elements Ni (∼40 μg/l) and Cu (∼10 μg/l) limiting the potential utilisation of Lake Xolotlan for nearby Managua.     

Thermal history of the unsaturated zone at Yucca Mountain,Nevada, USA

Secondary calcite, silica and minor amounts of fluorite deposited in fractures and cavities record the chemistry, temperatures, and timing of past fluid movement in the unsaturated zone at Yucca Mountain, Nevada, the proposed site of a high-level radioactive waste repository. The distribution and geochemistry of these deposits are consistent with low-temperature precipitation from meteoric waters that infiltrated at the surface and percolated down through the unsaturated zone. However, the discovery of fluid inclusions in calcite with homogenization temperatures (T_h) up to ∼80 °C was construed by some scientists as strong evidence for hydrothermal deposition. This paper reports the results of investigations to test the hypothesis of hydrothermal deposition and to determine the temperature and timing of secondary mineral deposition. Mineral precipitation temperatures in the unsaturated zone are estimated from calcite- and fluorite-hosted fluid inclusions and calcite δ¹⁸O values, and depositional timing is constrained by the ²⁰⁷Pb/²³⁵U ages of chalcedony or opal in the deposits. Fluid inclusion T_h from 50 samples of calcite and four samples of fluorite range from ∼35 to ∼90 °C. Calcite δ¹⁸O values range from ∼0 to ∼22‰ (SMOW) but most fall between 12 and 20‰. The highest T_h and the lowest δ¹⁸O values are found in the older calcite. Calcite T_h and δ¹⁸O values indicate that most calcite precipitated from water with δ¹⁸O values between −13 and −7‰, similar to modern meteoric waters.     

The arsenic source term for an in-pit uranium mine tailings facility and its long-term impact on the regional groundwater

Detailed field sampling and analyses and laboratory-based diffusion-cell experiments were used in conjunction with 3-D reactive transport modeling (MODFLOW and MT3D99) to quantify the fate and long-term (10 ka) transport of As in the Rabbit Lake In-pit Tailings Management Facility (RLITMF), northern Saskatchewan, Canada. The RLITMF (300 m × 425 m × 90 m thick) was engineered to ensure solute transport within the RLITMF is dominated by diffusion. Concentrations of As in the tailings pore fluids ranged from 0.24 to 140 mg/L (n = 43). Arsenic speciation analyses indicate 90% of this arsenic exists as As⁵⁺. This observation is supported by pH–Eh measurements of pore fluids (n = 135). Geochemical analyses yielded a strong inverse correlation between the Fe/As molar ratio in the tailings solids and the corresponding concentration of dissolved As, which is attributed to the adsorption of As to secondary 2-line ferrihydrite present in the tailings. Diffusion-cell testing yielded values for the effective diffusion coefficient, sorption coefficient, and effective porosity of As in the tailings of 4.5 × 10⁻¹⁰ m²/s, 2–4 cm³/g and 0.36, respectively. Reactive transport simulations using the field and laboratory data show adsorption of As to the tailings and diffusive transport of dissolved As in the tailings should reduce the source term concentration of As to between 40% and 70% of the initial concentrations over the 10 ka simulation period. Based on these simulations, the As concentrations in the regional groundwater, 50 m down gradient of the tailings facility, should be maintained at background concentrations of 0.001 mg/L over the 10 ka period. These findings suggest the engineered in-pit disposal of U mine tailings can provide long-term protection for the local groundwater regime from As contamination.     

The hydrochemistry of a semi-arid pan basin case study: Sua Pan,Makgadikgadi, Botswana

This study presents results on the fluid and salt chemistry for the Makgadikgadi, a substantial continental basin in the semi-arid Kalahari. The aims of the study are to improve understanding of the hydrology of such a system and to identify the sources of the solutes and the controls on their cycling within pans. Sampling took place against the backdrop of unusually severe flooding as well as significant anthropogenic extraction of subsurface brines. This paper examines in particular the relationship between the chemistry of soil leachates, fresh stream water, salty lake water, surface salts and subsurface brines at Sua Pan, Botswana with the aim of improving the understanding of the system’s hydrology. Occasionally during the short wet season (December–March) surface water enters the saline environment and precipitates mostly calcite and halite, as well as dolomite and traces of other salts associated with the desiccation of the lake. The hypersaline subsurface brine (up to TDS 190,000 mg/L) is homogenous with minor variations due to pumping by BotAsh mine (Botswana Ash (Pty) Ltd.), which extracts 2400 m³ of brine/h from a depth of 38 m. Notable is the decrease in TDS as the pumping rate increases which may be indicative of subsurface recharge by less saline water. Isotope chemistry for Sr (⁸⁷Sr/⁸⁶Sr average 0.722087) and S (δ³⁴S average 34.35) suggests subsurface brines have been subject to a lithological contribution of undetermined origin. Recharge of the subsurface brine from surface water including the Nata River appears to be negligible.     

Geochemistry and mineralogy of arsenic in (natural) anaerobic groundwaters

Here new data from field bioremediation experiments and geochemical modeling are reported to illustrate the principal geochemical behavior of As in anaerobic groundwaters. In the field bioremediation experiments, groundwater in Holocene alluvial aquifers in Bangladesh was amended with labile water-soluble organic C (molasses) and MgSO₄ to stimulate metabolism of indigenous SO₄-reducing bacteria (SRB). In the USA, the groundwater was contaminated by Zn, Cd and SO₄, and contained <10 μg/L As under oxidized conditions, and a mixture of sucrose and methanol were injected to stimulate SRB metabolism. In Bangladesh, groundwater was under moderately reducing conditions and contained ∼10 mg/L Fe and ∼100 μg/L As. In the USA experiment, groundwater rapidly became anaerobic, and dissolved Fe and As increased dramatically (As > 1000 μg/L) under geochemical conditions consistent with bacterial Fe-reducing conditions. With time, groundwater became more reducing and biogenic SO₄ reduction began, and Cd and Zn were virtually completely removed due to precipitation of sphalerite (ZnS) and other metal sulfide mineral(s). Following precipitation of chalcophile elements Zn and Cd, the concentrations of Fe and As both began to decrease in groundwater, presumably due to formation of As-bearing FeS/FeS₂. By the end of the six-month experiment, dissolved As had returned to below background levels. In the initial Bangladesh experiment, As decreased to virtually zero once biogenic SO₄ reduction commenced but increased to pre-experiment level once SO₄ reduction ended. In the ongoing experiment, both SO₄ and Fe(II) were amended to groundwater to evaluate if FeS/FeS₂ formation causes longer-lived As removal. Because As-bearing pyrite is the common product of SRB metabolism in Holocene alluvial aquifers in both the USA and Southeast Asia, it was endeavored to derive thermodynamic data for arsenian pyrite to better predict geochemical processes in naturally reducing groundwaters. Including the new data for arsenian pyrite into Geochemist’s Workbench, its stability field completely dominates in reducing Eh–pH space and “displaces” other As-sulfides (orpiment, realgar) that have been implied to be important in previous modeling exercises and reported in rare field conditions.     

National-scale estimation of potentially harmful element ambient background concentrations in topsoil using parent material classified soil:stream–sediment relationships

Regulatory authorities require estimates of ambient background concentrations (ABCs) of potentially harmful elements (PHEs) in topsoil; such data are currently not available in many countries. High resolution soil geochemical data exist for only part of England and Wales, whilst stream sediment data cover the entire landscape. A novel methodology is presented for estimating soil equivalent ABCs for PHEs from high-resolution (HR) stream sediment geochemical data grouped by common parent materials (PM), using arsenic (As) as an example. Geometric mean (GM) values for local PM groups are used to investigate different approaches for transforming sediment to soil equivalent concentrations. Holdout validation is used to assess: (i) the optimum number of samples for calculating local GM values, and (ii) the optimum scale at which to group data when using linear regression analysis to estimate GM soil ABCs from local sediment geochemical values. Holdout validation showed that the smallest differences were generally observed when five observations were used to calculate the GM and that these should be grouped over the smallest possible area in order to encompass soils over PMs with elevated GM As concentrations. Geometric mean ABCs are estimated and mapped for As in mineral soil across all of England and Wales within delineations of PM polygons. Errors for the estimation of soil equivalent GM As ABCs based on sediment data for an independent validation set were of a similar magnitude to those from holdout validation applied to the original data suggesting the approach is robust. The estimates of soil equivalent ABCs suggest that As exceeds the regulatory threshold used in risk assessments for residential land use (20 mg kg⁻¹) across 16% of the landscape of England and Wales. The applicability of the method for cognate landscapes, and potential refinements is discussed.