Abstract
The Organic Matter ENabled SEDiment model (OMEN-SED) is a one-dimensional, analytical reaction-transport model for early diagenesis in marine sediments. It explicitly resolves organic matter (OM) degradation and associated biogeochemical terminal electron acceptor, reduced species and nutrient dynamics in porous media under steady-state conditions. OMEN-SED has been specifically designed for coupling to global Earth system models and the analytical solution of the coupled set of mass conservation equations ensures the computational efficiency required for such a coupling. To find an analytical solution, OMEN-SED expresses all explicitly resolved biogeochemical processes as a function of OM degradation. The original version of OMEN-SED contains a relatively simple description of OM degradation based on two reactive OM classes, a so-called 2G model. However, such a simplified approach does not fully account for the widely observed continuous decrease in organic matter reactivity with burial depth/time. The reactive continuum model that accounts for the continuous distribution of organic compounds over the reactive spectrum represents an alternative and more realistic description but cannot be easily incorporated within the general OMEN-SED framework. Here, we extend the diagenetic framework of OMEN-SED with a multi-G approximation of the reactive continuum model (RCM) of organic matter degradation by using a finite but large number of OM fractions, each characterized by a distinct reactivity. The RCM and its multi-G approximation are fully constrained by only two free parameters, a and ν, that control the initial distribution of OM compounds over the reactivity spectrum. The new model is not only able to reproduce observed pore water profiles, sediment-water interface fluxes and redox zonation across a wide range of depositional environments but also provides a more realistic description of anaerobic degradation pathways. The added functionality extends the applicability of OMEN-SED to a broader range of environments and timescales, while requiring fewer parameters to simulate a wider spectrum of OM reactivities.
Original language | English |
---|---|
Pages (from-to) | 7155-7174 |
Number of pages | 20 |
Journal | Geoscientific Model Development |
Volume | 14 |
Issue number | 11 |
Early online date | 25 Nov 2021 |
DOIs | |
Publication status | Published - Nov 2021 |
Bibliographical note
Funding Information:2020 (grant no. C-CASCADES (643052)). Dominik Hülse was supported by a postdoctoral fellowship from the Simons Foundation (award ID 653829).
Funding Information:
Financial support. This research has been supported by Horizon
Funding Information:
Acknowledgements. We would like to thank Bernard Boudreau and one anonymous reviewer for their constructive critiques and insightful comments, which have improved the paper. Sandra Arndt and Philip Pika were supported by funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant (agreement no. 643052) (C-CASCADES).
Publisher Copyright:
© Copyright: