Article révisé par les pairs
Résumé : 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.