Permeability and Groundwater Flow Dynamics in Deep‐Reaching Orogenic Faults Estimated From Regional‐Scale Hydraulic Simulations

Numerical modeling is used to understand the regional scale flow dynamics of the fault-hosted orogenic geothermal system at the Grimsel Mountain Pass in the Swiss Alps. The model is calibrated against observations from thermal springs discharging in a tunnel some 250 m underneath Grimsel Pass to derive estimates for the bulk permeability of the fault. Simulations confirm that without the fault as a hydraulic conductor the thermal springs would not exist. Regional topography alone drives meteoric water in a single pass through the fault plane where it penetrates to depths exceeding 10 km and acquires temperatures in excess of 250°C. Thermal constraints from the thermal springs at Grimsel Pass suggest bulk fault permeabilities in the range of 2e−15 m2–4.8e−15 m2. Reported residence times of >30,000 and 7 years for the deep geothermal and shallow groundwater components in the thermal spring water, respectively, suggest fault permeabilities of around 2.5e−15 m2. We show that the long residence time of the deep geothermal water is likely a consequence of low recharge rates during the last glaciation event in the Swiss Alps, which started some 30,000 years ago. Deep groundwater discharging at Grimsel Pass today thus infiltrated the Grimsel fault prior to the last glaciation event. The range of permeabilities estimated from observational constraints is fully consistent with a subcritical single-pass flow system in the fault plane

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal-hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Previous studies have shown that fault-bounded circulation of meteoric water is driven to depths of ∼10 km by the high alpine topography. Simulations suggest that the current single-pass flow is typical of interglacial periods, during which (a) meteoric recharge into the fault is high (above tens of centimeters per year), (b) conditions are at or somewhat below the critical Rayleigh number, and (c) hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to <1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system, with single-pass flow confined to shallow depths while non-Rayleigh convection occurred deeper in the fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits flux through the fault

Reaction Mechanism and Water/Rock Ratios Involved in Epidosite Alteration of the Oceanic Crust

Epidosites are a prominent type of subseafloor hydrothermal alteration of basalts in ophiolites and greenstone belts, showing an end-member mineral assemblage of epidote + quartz + titanite + Fe-oxide. Epidosites are known to form within crustal-scale upflow zones and their fluids have been proposed as deep equivalents of black-smoker seafloor vent fluids. Proposals of the mass of fluid per mass of rock (W/R ratio) needed to form epidosites are contradictory, varying from 20 (Sr isotopes) to > 1,000 (Mg mobility). To test these proposals we have conducted a petrographic, geochemical and reactive-transport numerical simulation study of the chemical reaction that generates km3-size epidosite zones within the lavas and sheeted dike complex of the Samail ophiolite, Oman. At 250–400°C the modeled epidosite-forming fluid has near-neutral pH (∼ 5.2), high fO2, low sulfur and very low Fe (10−6 mol/kg) contents. These features argue against a genetic link with black-smoker fluids. Chemical buffering by the epidosite fluid enriches the precursor spilites in Ca and depletes them in Na and Mg. Completion of the spilite-to-epidosite reaction requires enormous W/R ratios of 700–∼40,000, depending on initial Mg content and temperature. Collectively, the variably altered rocks in the Samail epidosite zones record flow of ∼1015 kg of fluid through each km3 of precursor spilite rock. This fluid imposed on the epidosite an Sr-isotope signature inherited from the previous rock-buffered chemical evolution of the fluid through the oceanic crust, thereby explaining the apparently contradictory low W/R ratios based on Sr isotopes

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions

3D modelling of long-term sulfide corrosion of copper canisters in a spent nuclear fuel repository

Copper canisters are a central component in the safety of the Finnish spent fuel repository concept (KBS-3), where the main corrodent potentially affecting the canister integrity is sulfide. In this study, a 3D numerical model is developed to assess the evolution of sulfide fluxes and the spatially resolved canister corrosion depths for the Finnish spent nuclear fuel repository concept. The backfilled tunnel and the disposal hole are implemented using repository geometries, with sulfide being produced at their interface with the rock (excavation damaged zone) by sulfate reducing bacteria (SRB). Recent experimental findings regarding the microbial sulfate reduction process as well as the scavenging of sulfide via iron (oxy)hydroxides are incorporated in the reactive transport model. Long-term simulations are performed, predicting a heterogeneous corrosion of the canister with a max. corrosion depth of 1.3 mm at the bottom corner after one million years. The evolution of sulfide fluxes shows two main phases, depending on the source of sulfate: first sulfate is supplied by the dissolution of gypsum from the bentonite barriers, followed by a steady, low-level supply from the groundwater. Sensitivity cases demonstrate that both the organic carbon and Fe(III) oxide contents in the bentonite are critical to the corrosion evolution, by being the main electron donor for SRB activities and the major sulfide scavenger in the bentonite, respectively. The backfilled tunnel contributes little to the flux of corrosive sulfide to the canister due to the attenuation by Fe(III)-oxides/hydroxides but induces a notable flux of sulfate into the disposal hole

Simulating Donnan equilibria based on the Nernst-Planck equation

Understanding ion transport through clays and clay membranes is important for many geochemical and environmental applications. Ion transport is affected by electrostatic forces exerted by charged clay surfaces. Anions are partly excluded from pore water near these surfaces, whereas cations are enriched. Such effects can be modeled by the Donnan approach. Here we introduce a new, comparatively simple way to represent Donnan equilibria in transport simulations. We include charged surfaces as immobile ions in the balance equation and calculate coupled transport of all components, including the immobile charges, with the Nernst-Planck equation. This results in an additional diffusion potential that influences ion transport, leading to Donnan ion distributions while maintaining local charge balance. The validity of our new approach was demonstrated by comparing Nernst-Planck simulations using the reactive transport code Flotran with analytical solutions available for simple Donnan systems. Attention has to be paid to the numerical evaluation of the electrochemical migration term in the Nernst-Planck equation to obtain correct results for asymmetric electrolytes. Sensitivity simulations demonstrate the influence of various Donnan model parameters on simulated anion accessible porosities. It is furthermore shown that the salt diffusion coefficient in a Donnan pore depends on local concentrations, in contrast to the aqueous salt diffusion coefficient. Our approach can be easily implemented into other transport codes. It is versatile and facilitates, for instance, assessing the implications of different activity models for the Donnan porosity

Insights into the evolution of an oceanic hydrothermal system and a method for constraining estimates of the vigor of hydrothermal convection

The permeability of the oceanic crust is the primary hydrologic parameter that controls the geometry, the vigor, and the duration of hydrothermal fluid flow at mid-ocean ridges. The coupling of fluid flow, temperature, and chemistry and the effect of the permeability on these coupled processes are assessed to determine whether geochemical data can be used to constrain estimates of basement permeabilities and the vigor of convection. The coupling of flow, temperature, and chemistry is investigated for an open, sedimented ridge setting from the onset of fluid motion after an initial state of a conductive temperature distribution and fluid stagnation, to a near steady-state convective system. Fluid residence times, physical water/rock ratios, temperature conditions, recharge and discharge rates, flow geometries, and the degree of fluid mixing are calculated for the evolving hydrothermal system and the influence of these parameters on mineral alteration and aqueous phase concentrations are discussed. The chemical evolution of the system suggests that despite differences in alteration patterns and the intensity of alteration for different basement permeabilities, the nature of the alteration reactions is unlikely to be a useful parameter for constraining the vigor of convection. Estimates of fluid velocity can be obtained using a formulation that relates the mass transfer between the solid and the fluid phase along the fluid's flowpath to the average velocity along the flowpath. Different regions in the hydrothermal system and a range of chemical species are examined to assess their usefulness in constraining estimates of average flow velocities. The results of these calculations suggest that the mass transfer of aqueous silica may be useful for estimating fluid flow velocities in hydrothermal systems, in particular in those regions of the system at or near quartz equilibrium so that aqueous silica concentration is buffered by quartz and correlated with the temperature distribution.Science, Faculty ofEarth, Ocean and Atmospheric Sciences, Department ofGraduat

Quantification of 3‐D Thermal Anomalies From Surface Observations of an Orogenic Geothermal System (Grimsel Pass, Swiss Alps)

Geothermal systems in amagmatic orogens involve topography‐driven infiltration of meteoric water up to 10 km deep into regional‐scale faults and exfiltration of the heated water in surface springs. The thermal anomalies along the upflow zones have not been quantified, yet they are key to estimating the geothermal exploitation potential of such systems. Here we quantify the three‐dimensional heat anomaly below the orogenic geothermal system at Grimsel Pass, Swiss Alps, where warm springs emanate from an exhumed, fossil hydrothermal zone. We use discharge rates and temperatures of the springs, temperature measurements along a shallow tunnel, and the formation temperature and depth of the fossil system to constrain coupled thermal–hydraulic numerical simulations of the upflow zone. The simulations reveal that upflow rates act as a first‐order control on the temperature distribution and that the site is underlain by an ellipsoidal thermal plume enclosing 102–103 PJ of anomalous heat per km depth. When the fossil system was active (3.3 Ma), the thermal plume was double its present size, corresponding to a theoretical petrothermal power output of 30–220 MW, with the 120 °C threshold for geothermal electricity production situated at less than 2‐km depth. We conclude that mountainous orogenic belts without igneous activity and even with only low background geothermal gradients typical of waning orogens are surprisingly promising plays for petrothermal power production. Our study implies exploration should focus on major valley floors because there the hydraulic head gradients and thus upflow rates and heat anomalies reach maximum values

Reactive mass transport modelling of a three-dimensional vertical fault zone with a finger-like convective flow regime

For a three-dimensional vertically-oriented fault zone, we consider the coupled effects of fluid flow, heat transfer and reactive mass transport, to investigate the patterns of fluid flow, temperature distribution, mineral alteration and chemically induced porosity changes. We show, analytically and numerically, that finger-like convection patterns can arise in a vertically-oriented fault zone. The onset and patterns of convective fluid flow are controlled by the Rayleigh number which is a function of the thermal properties of the fluid and the rock, the vertical temperature gradient, and the height and the permeability of the fault zone. Vigorous fluid flow causes low temperature gradients over a large region of the fault zone. In such a case, flow across lithological interfaces becomes the most important mechanism for the formation of sharp chemical reaction fronts. The degree of rock buffering, the extent and intensity of alteration, the alteration mineralogy and in some cases the formation of ore deposits are controlled by the magnitude of the flow velocity across these compositional interfaces in the rock. This indicates that alteration patterns along compositional boundaries in the rock may provide some insights into the convection pattern. The advective mass and heat exchanges between the fault zone and the wallrock depend on the permeability contrast between the fault zone and the wallrock. A high permeability contrast promotes focussed convective flow within the fault zone and diffusive exchange of heat and chemical reactants between the fault zone and the wallrock. However, a more gradual permeability change may lead to a regional-scale convective flow system where the flow pattern in the fault affects large-scale fluid flow, mass transport and chemical alteration in the wallrock