Isotope Ratio – Discharge Relationships of Solutes Derived From Weathering Reactions

To date, the vast majority of studies seeking to link discharge to solute concentrations have been based on representations of fluid age distributions in watersheds that are time-invariant. As increasingly detailed spatial and temporal datasets become available for weathering-derived riverine solute concentrations, the capacity to link this mass flux to transient routing of reactive fluids through Critical Zone environments is vital to quantitative interpretation. Relationships between fluid age distributions and the stable isotope ratios of these geogenic solutes are even less developed, yet these signatures are vital to parsing the suite of water-rock-life interactions that create concentration-discharge relationships. Here we offer the first merging of a hydrological model featuring time-variant fluid age distributions with a geochemical model for isotopically fractionating weathering reactions. Using SiO2(aq) and the corresponding silicon isotope ratio δ30Si as an example, we show that the stable isotope signatures of riverine solutes produced by weathering reactions reflect a component of the fluid age distribution that is unique to the corresponding solute concentrations. This distinct sensitivity is the result of a stronger link between isotope ratios and the age distribution parameters describing a given watershed. This novel modeling framework is used to provide a quantitative basis for the interpretation of SiO2(aq) and δ30Si in six low-order streams spread across a diversity of climates, geologies, and ecosystems. To our knowledge, this is the first forward and process-based model to describe the isotopic signatures of solutes derived from weathering reactions in watersheds subject to time-varying discharge

An emulation-based approach for interrogating reactive transport models

We present a new approach to understand the interactions among different chemical and biological processes modelled in environmental reactive transport models (RTMs) and explore how the parameterisation of these processes influences the results of multi-component RTMs. We utilize a previously published RTM consisting of 20 primary species, 20 secondary complexes, 17 mineral reactions and 2 biologically-mediated reactions which describes bio-stimulation using sediment from a contaminated aquifer. We choose a subset of the input parameters to vary over a range of values. The result is the construction of a new dataset that describes the model behaviour over a range of environmental conditions. Using this dataset to train a statistical model creates an emulator of the underlying RTM. This is a condensed representation of the original RTM that facilitates rapid exploration of a broad range of environmental conditions and sensitivities. As an illustration of this approach, we use the emulator to explore how varying the boundary conditions in the RTM describing the aquifer impacts the rates and volumes of mineral precipitation. A key result of this work is the recognition of an unanticipated dependency of pyrite precipitation on pCO2 in the injection fluid due to the stoichiometry of the microbially-mediated sulphate reduction reaction. This complex relationship was made apparent by the emulator, while the underlying RTM was not specifically constructed to create such a feedback. We argue that this emulation approach to sensitivity analysis for RTMs may be useful in discovering such new coupled sensitives in geochemical systems and for designing experiments to optimise environmental remediation. Finally, we demonstrate that this approach can maximise specific mineral precipitation or dissolution reactions by using the emulator to find local maxima, which can be widely applied in environmental systems.</p

Modelling the Effects of Non-Steady State Transport Dynamics on the Sulfur and Oxygen Isotope Composition of Sulfate in Sedimentary Pore Fluids

We present the results of an isotope-enabled reactive transport model of a sediment column undergoing active microbial sulfate reduction to explore the response of the sulfur and oxygen isotopic composition of sulfate under perturbations to steady state. In particular, we test how perturbations to steady state influence the cross plot of δ34S and δ18O for sulfate. The slope of the apparent linear phase (SALP) in the cross plot of δ34S and δ18O for sulfate has been used to infer the mechanism, or metabolic rate, of microbial metabolism, making it important that we understand how transient changes might influence this slope. Tested perturbations include changes in boundary conditions and changes in the rate of microbial sulfate reduction in the sediment. Our results suggest that perturbations to steady state influence the pore fluid concentration of sulfate and the δ34S and δ18O of sulfate but have a minimal effect on SALP. Furthermore, we demonstrate that a constant advective flux in the sediment column has no measurable effect on SALP. We conclude that changes in the SALP after a perturbation are not analytically resolvable after the first 5% of the total equilibration time. This suggests that in sedimentary environments the SALP can be interpreted in terms of microbial metabolism and not in terms of environmental parameters

The Carbon-Sulfur Link in the Remineralization of Organic Carbon in Surface Sediments

Here we present the carbon isotopic composition of dissolved inorganic carbon (DIC) and the sulfur isotopic composition of sulfate, along with changes in sulfate concentrations, of the pore fluid collected from a series of sediment cores located along a depth transect on the Iberian Margin. We use these data to explore the coupling of microbial sulfate reduction (MSR) to organic carbon oxidation in the uppermost (up to nine meters) sediment. We argue that the combined use of the carbon and sulfur isotopic composition, of DIC and sulfate respectively, in sedimentary pore fluids, viewed through a δ13CDIC vs. δ34SSO4 cross plot, reveals significant insight into the nature of carbon-sulfur coupling in marine sedimentary pore fluids on continental margins. Our data show systemic changes in the carbon and sulfur isotopic composition of DIC and sulfate (respectively) where, at all sites, the carbon isotopic composition of the DIC decreases before the sulfur isotopic composition of sulfate increases. We compare our results to global data and show that this behavior persists over a range of sediment types, locations and water depths. We use a reactive-transport model to show how changes in the amount of DIC in seawater, the carbon isotopic composition of organic matter, the amount of organic carbon oxidation by early diagenetic reactions, and the presence and source of methane influence the carbon and sulfur isotopic composition of sedimentary pore fluids and the shape of the δ13CDIC vs. δ34SSO4 cross plot. The δ13C of the DIC released during sulfate reduction and sulfate-driven anaerobic oxidation of methane is a major control on the minimum δ13CDIC value in the δ13CDIC vs. δ34SSO4 cross plot, with the δ13C of the organic carbon being important during both MSR and combined sulfate reduction, sulfate-driven AOM and methanogenesis

Isotopic insights into microbial sulfur cycling in oil reservoirs

Microbial sulfate reduction in oil reservoirs (biosouring) is often associated with secondary oil production where seawater containing high sulfate concentrations (~28 mM) is injected into a reservoir to maintain pressure and displace oil. The sulfide generated from biosouring can cause corrosion of infrastructure, health exposure risks, and higher production costs. Isotope monitoring is a promising approach for understanding microbial sulfur cycling in reservoirs, enabling early detection of biosouring, and understanding the impact of souring. Microbial sulfate reduction is known to result in large shifts in the sulfur and oxygen isotope compositions of the residual sulfate, which can be distinguished from other processes that may be occurring in oil reservoirs, such as precipitation of sulfate and sulfide minerals. Key to the success of this method is using the appropriate isotopic fractionation factors for the conditions and processes being monitored. For a set of batch incubation experiments using a mixed microbial culture with crude oil as the electron donor, we measured a sulfur fractionation factor for sulfate reduction of −30‰. We have incorporated this result into a simplified 1D reservoir reactive transport model to highlight how isotopes can help discriminate between biotic and abiotic processes affecting sulfate and sulfide concentrations. Modeling results suggest that monitoring sulfate isotopes can provide an early indication of souring for reservoirs with reactive iron minerals that can remove the produced sulfide, especially when sulfate reduction occurs in the mixing zone between formation waters (FW) containing elevated concentrations of volatile fatty acids (VFAs) and injection water (IW) containing elevated sulfate. In addition, we examine the role of reservoir thermal, geochemical, hydrological, operational and microbiological conditions in determining microbial souring dynamics and hence the anticipated isotopic signatures

A reactive transport model for geochemical mitigation of CO2 leaking into a confined aquifer

Long-term storage of anthropogenic CO2 in the subsurface generally assumes that caprock formations will serve as physical barriers to upward migration of CO2. However, as a precaution and to provide assurances to regulators and the public, monitoring is used detect any unexpected leakage from the storage reservoir. If a leak is found, the ability to rapidly deploy mitigation measures is needed. Here we use the TOUGHREACT code to develop a series of two-dimensional reactive transport simulations of the hydrogeochemical characteristics of a newly formed CO2 leak into an overlying aquifer. Using this model, we consider: (1) geochemical shifts in formation water indicative of a leak; (2) hydrodynamics of pumping wells in the vicinity of a leak; and (3) delivery of a sealant to a leak through an adjacent well bore.Our results demonstrate that characteristic shifts in pH and dissolved inorganic carbon can be detected in the aquifer prior to the breakthrough of supercritical CO2, and could offer a potential means of identifying small and newly formed leaks. Pumping water into the aquifer in the vicinity of the leak provides a hydrodynamic control that can temporarily mitigate the flux rate of CO2 and facilitate delivery of a sealant to the location of the caprock defect. Injection of a fluid-phase sealant through the pumping well is demonstrated by introduction of a silica-bearing alkaline flood, resulting in precipitation of amorphous silica in areas of neutral to acidic pH. Results show that a decrease in permeability of several orders of magnitude can be achieved with a high molar volume sealant, such that CO2 flux rate is decreased. However, individual simulation results are highly contingent upon both the properties of the sealant, the porosity-permeability relationship employed in the model, and the relative flux rates of CO2 and alkaline flood introduced into the aquifer. These conclusions highlight the need for both experimental data and controlled field tests to constrain modelling predictions

Isotopic communication across the water - rock interface: Preservation in solids and signatures in fluids

International audienceWater - rock interactions are characteristically associated with isotopic fractionation between phases, yet their remains significant uncertainty in the degree to which these isotope ratios record the conditions of mineral formation or equilibration, and thus their fidelity as proxy records. Here, we present a novel generic modeling framework allowing sub-grid scale tracking of the isotopic signatures recorded in newly formed minerals over the timescale of precipitation and continued isotopic exchange between phases. We demonstrate that the surface area of the solid phase strongly influences the extent to which isotopes reflect the conditions of formation through time

Isotopic communication across the water – rock interface: Preservation in solids and signatures in fluids

Water – rock interactions are characteristically associated with isotopic fractionation between phases, yet their remains significant uncertainty in the degree to which these isotope ratios record the conditions of mineral formation or equilibration, and thus their fidelity as proxy records. Here, we present a novel generic modeling framework allowing sub-grid scale tracking of the isotopic signatures recorded in newly formed minerals over the timescale of precipitation and continued isotopic exchange between phases. We demonstrate that the surface area of the solid phase strongly influences the extent to which isotopes reflect the conditions of formation through time

Multi-Scale Microfluidics for Transport in Shale Fabric

We develop a microfluidic experimental platform to study solute transport in multi-scale fracture networks with a disparity of spatial scales ranging between two and five orders of magnitude. Using the experimental scaling relationship observed in Marcellus shales between fracture aperture and frequency, the microfluidic design of the fracture network spans all length scales from the micron (1 &mu;) to the dm (10 dm). This intentional `tyranny of scales&rsquo; in the design, a determining feature of shale fabric, introduces unique complexities during microchip fabrication, microfluidic flow-through experiments, imaging, data acquisition and interpretation. Here, we establish best practices to achieve a reliable experimental protocol, critical for reproducible studies involving multi-scale physical micromodels spanning from the Darcy- to the pore-scale (dm to &mu;m). With this protocol, two fracture networks are created: a macrofracture network with fracture apertures between 5 and 500 &mu;m and a microfracture network with fracture apertures between 1 and 500 &mu;m. The latter includes the addition of 1 &mu;m &lsquo;microfractures&rsquo;, at a bearing of 55&deg;, to the backbone of the former. Comparative analysis of the breakthrough curves measured at corresponding locations along primary, secondary and tertiary fractures in both models allows one to assess the scale and the conditions at which microfractures may impact passive transport