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.     

Modeling preferential flow in reactive barriers: implications for performance and design

Reactive barriers are passive and in situ ground water treatment systems. Heterogeneities in hydraulic conductivity (K) within the aquifer-reactive barrier system will result in higher flux rates, and reduced residence times, through portions of the barrier. These spatial variations in residence time will affect the treatment capacity of the barrier. A numerical flow model was used to evaluate the effects of spatial variations in K on preferential flow through barriers. The simulations indicate that the impact of heterogeneities in K will be a function of their location and distribution; the more localized the high K zone, the greater the preferential flow. The geometry of the reactive barrier will also strongly influence flow distribution. Aquifer heterogeneities will produce greater preferential flow in thinner barriers compared to thicker barriers. If the barrier K is heterogeneous, greater preferential flow will occur in thicker barriers. The K of the barrier will affect the flow distribution; decreasing the K of the barrier can result in more even distribution of flow. Results indicate that less variable flow will be attained utilizing thicker, homogeneous barriers. The addition of homogeneous zones to thinner barriers will be effective at redistributing flow only if installed immediately adjacent to both the up- and downgradient faces of the barrier.     

安徽铜陵凤凰山夕卡岩型铜矿床中辉钼矿Re-Os同位素年龄及其地质意义

铜陵矿集区是我国长江中下游Cu-Au-Fe-Mo成矿带中最重要的有色金属基地之一,凤凰山矿床是铜陵矿集区的重要组成部分,为一个典型的夕卡岩型铜矿床。本文利用Re-Os同位素定年方法对凤凰山铜矿床进行了成矿时代测定,获得了辉钼矿的Re-Os同位素模式年龄范围为139.1±2.4～142.0±2.2Ma,等时线年龄为141.1±1.4Ma,与矿区内石英二长闪长岩和花岗闪长岩SHRIMP锆石U-Pb年龄(144.2±2.3Ma)相吻合,也与铜陵地区其他矿田的成矿时代基本一致,可能为岩石圈减薄事件的成矿响应。     

赣南营前岩体的年代学、地球化学、Sr-Nd-Hf同位素组成及其地质意义

赣南营前岩体由早期似斑状花岗闪长岩和晚期细粒黑云母花岗岩组成,在前者中普遍发育暗色基性包体。SHRIMP锆石U-Pb测年显示似斑状花岗闪长岩为172.2±3.1Ma,细粒黑云母花岗岩为168±3Ma,包体为167.6±3.3Ma。所有的岩石均为钙碱性系列,稀土元素显示LREE富集,HREE亏损,Eu负异常并不明显的特征,Nb-Ta负异常显示壳源物质的特征。晚期细粒黑云母花岗岩、早期似斑状花岗闪长岩、暗色包体的(87Sr/86Sr)i和εNd(t)值分别为0.70885～0.70965、0.71058～1.46393、0.70788～0.70923和-6.38～-5.31、-6.51～-5.17、-5.47～-4.31;锆石εHf(t)值分别为-12.4～-5.7、-9.2～-4.2、-9.4～+0.1。结合岩相学、地球化学研究,初步认为营前岩体的两期花岗质岩石源于地壳重熔,地幔底侵导致地壳熔融形成营前岩体,基性岩浆加入到未冷却的酸性岩浆中形成了目前所见到的暗色包体。结合区域上的研究资料,认为中侏罗世期间,沿"十杭"断裂带及与其相平行的一些深断裂发生了板块的撕裂,这些缝合带是岩石圈相对薄弱的地带,有利于幔源物质底侵上升以及壳幔相互作用。     

新疆阿尔泰克因布拉克铜锌矿床地质特征及成矿作用

克因布拉克中型铜锌矿床赋存于早二叠世花岗岩外接触带的上志留统-下泥盆统康布铁堡组黑云石英片岩、变质石英砂岩中。矿床的形成经历了夕卡岩期、热液期和表生期,铜锌矿主要形成于热液期。矿石中石英和方解石流体包裹体划分为H_2O-NaCl型和H_2O-CO_2(±CH_4/N_2)-NaCl型。成矿温度变化于146～448℃,在170℃、270℃和350℃出现峰值;流体盐度变化于0.2%～46.9%NaCl_(eq),峰值为1.5%NaCl_(eq)和5.5%NaCl_(eq);流体的密度0.55～1.19g/cm~3。硫化物的δ~(34)S集中变化于-8.4‰～1.9‰,峰值为0‰,表明硫来自岩浆。石英和方解石δD_(SMOW)介于-130‰～-79‰,δ~(18)O_(SMOW)值介于8.0‰～11.6‰,δ~(18)O_(H2O)值为-1.7‰～4.43‰,表明成矿流体主要是岩浆水,混合大气降水。方解石中δ~(13)C_(PDB)变化于-5.3‰～-1.1‰,暗示碳来自花岗质岩浆。成矿时代为早中二叠世,成矿作用与花岗质岩浆期后的热液活动有关。     

Numerical investigation of road salt impact on an urban wellfield

The impact of road salt on a wellfield in a complex glacial moraine aquifer system is studied by numerical simulation. The moraine underlies an extensive urban and industrial landscape, which draws its water supply from >20 wellfields, several of which are approaching or have exceeded the drinking water limit for chloride. The study investigates the mechanisms of road salt infiltration, storage, and transport in the subsurface and assesses the effectiveness of mitigation measures designed to reduce the impact. The three-dimensional transport model accounts for increases in salt loading, as well as growth of the urbanized area and road network over the past 50 years. The simulations, which focus on one impacted wellfield, show chloride plumes originating mainly at arterial roads and migrating through aquitard windows into the water supply aquifers. The results suggest that the aquifer system contains a large and heterogeneously distributed mass of chloride and that concentrations in the aquifer can be substantially higher than the concentrations in the well water. Future impact scenarios indicate that although the system responds rapidly to reductions in salt loading, the residual chloride mass may take decades to flush out, even if road salting were discontinued. The implications with respect to urban wellfields in typical snow-belt areas are discussed.     

Controls on permafrost thaw in a coupled groundwater-flow and heat-transport system: Iqaluit Airport,Nunavut, Canada

Numerical simulations of groundwater flow and heat transport are used to provide insight into the interaction between shallow groundwater flow and thermal dynamics related to permafrost thaw and thaw settlement at the Iqaluit Airport taxiway, Nunavut, Canada. A conceptual model is first developed for the site and a corresponding two-dimensional numerical model is calibrated to the observed ground temperatures. Future climate-warming impacts on the thermal regime and flow system are then simulated based on climate scenarios proposed by the Intergovernmental Panel on Climate Change (IPCC). Under climate warming, surface snow cover is identified as the leading factor affecting permafrost degradation, including its role in increasing the sensitivity of permafrost degradation to changes in various hydrogeological factors. In this case, advective heat transport plays a relatively minor, but non-negligible, role compared to conductive heat transport, due to the significant extent of low-permeability soil close to surface. Conductive heat transport, which is strongly affected by the surface snow layer, controls the release of unfrozen water and the depth of the active layer as well as the magnitude of thaw settlement and frost heave. Under the warmest climate-warming scenario with an average annual temperature increase of 3.23 °C for the period of 2011–2100, the simulations suggest that the maximum depth of the active layer will increase from 2 m in 2012 to 8.8 m in 2100 and, over the same time period, thaw settlement along the airport taxiway will increase from 0.11 m to at least 0.17 m.

     

Mass fluxes of xenobiotics below cities: challenges in urban hydrogeology

Urban areas are the focus of major ecological, social and economical activity. They are thus also prime locations of increasing conflict with regard to water use and water protection. As a direct and/or indirect consequence of urban land use and human activity, urban water systems are frequently polluted with organic contaminants including waste water-born xenobiotics such as pharmaceuticals, personal care products (collectively known as PPCPs) and endocrine-active substances. This study reviews new integrated methodologies including flux calculations as well as chemical investigations for determining the impact of human activities on urban water systems and on processes within the urban watershed. The use of indicator substances, representing different contaminant sources and pathways, integral pumping tests and mass balance approaches are suitable alternatives within these environments. The issues are explored using contaminant mass balance examples from Halle/Saale and Leipzig, Germany.     

Hydrogeology of a coal-seam gas exploration area, southeastern British Columbia, Canada: Part 1. Groundwater flow systems

Discovery of high contents of methane gas in coals of the Mist Mountain Formation in the Elk River valley, southeastern British Columbia, Canada, has led to increased exploration activity for coal-seam gas (CSG). CSG production requires groundwater abstraction to depressurize the coal beds and to facilitate methane flow to the production wells. Groundwater abstraction will have hydrodynamic effects on the flow system, and an understanding of the groundwater flow system is needed to evaluate these effects. The purpose of this paper is to describe the groundwater flow system in the area by means of a groundwater flow model and interpretation of hydrochemical and isotopic analyses of groundwater and surface water. Groundwater flow for the Weary Creek exploration area is modeled in two vertical sections. The model domains, based on classic upland–lowland conceptual flow models, are approximately 10,000 m long and 4,000 m deep. Each consists of a fixed water-table boundary and no-flow boundaries along the traces of major faults. Steady-state groundwater flow is calibrated to hydraulic-head, streamflow, and groundwater-recharge data. Simulated steady-state velocity fields define regional and local flow components consistent with the conceptual model. The results are consistent with regional trends in δ²H, δ¹⁸O, tritium, and TDS, which define two distinct groundwater groups (A and B) and a third of intermediate composition. An active, shallow, local flow component (group A) is recharged in beds cropping out along subdued ridges; this component discharges as seeps along lower and mid-slope positions in the southern part of the study area. The waters are tritiated, relatively enriched in δ²H and δ¹⁸O, and have low TDS. A deeper regional flow component (group B), which originates at a higher altitude and which discharges to the Elk River valley bottom, is characterized by non-tritiated groundwater with relatively depleted δ²H and δ¹⁸O, and higher TDS. Groundwater contributes less than 10% of the total direct flow to the Elk River, as indicated by flow measurements and by the absence of group A and group B characteristics in the river water. Thus it is hypothesized that groundwater extraction during CSG production will have little impact on the river. The groundwater flow model developed in this work is used in a companion paper to further test this hypothesis. Electronic Publication     

Well vulnerability: a quantitative approach for source water protection

The concept of vulnerability of drinking water sources is reviewed, and a quantitative approach for assessing well vulnerability for complex three-dimensional ground water systems is developed. The approach focuses on the relative expected impact of potential contaminant sources at unknown locations within a well capture zone, providing relative measures of intrinsic well vulnerability, including the expected times of arrival of a contaminant, the dispersion-related reduction in concentration, the time taken to breach a certain quality objective, and the corresponding exposure times. Thus, the result of the analysis includes the usual advective travel time information used in conventional wellhead protection analysis, plus a set of selected quantitative measures expressing the expected impact. The technique is based on adjoint theory and combines forward- and backward-in-time transport modeling using a standard numerical flow and transport code. The methodology is demonstrated using the case study of a complex glacial multiaquifer system in Ontario. The new approach will be useful in helping water managers develop more physically based and quantitative wellhead protection strategies.