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

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