Different proxies for the reactivity of aquatic sediments towards oxygen: A model assessment

Information on benthic carbon mineralization rates is often derived from the analysis of oxygen microprofiles in sediments. To enable a direct comparison of different sediment environments, it is often desirable to characterize sediments by a single proxy that expresses their “reactivity” towards oxygen. For this, there are three commonly used proxies: the oxygen penetration depth (OPD), the oxygen flux at the sediment-water interface (DOU), and the maximum volumetric oxygen consumption rate (R_max). The OPD can be directly determined from the oxygen depth profile, while the DOU is usually obtained by a linear fit to the oxygen gradient either in diffusive boundary layer. The oxygen consumption rate R_max requires the fitting of a reactive-transport model to the data profile. This article shows that the OPD alone is a suboptimal proxy, because it shows a strong dependence on the half-saturation constant K_s, and secondly, because it is sensitive to the particular re-oxidation conditions right above the oxic-anoxic interface. Similarly, the volumetric oxygen consumption rate R_max is rather strongly dependent on the kinetic model formulation employed. To show this we fitted three different (Bouldin, Blackman and Monod) kinetics to the same oxygen data profiles. When fitting these models, the R_max values obtained differed by 20% for exactly the same oxygen profile. Accordingly, if one reports R_max values, it is crucial to specify the kinetic model alongside. Overall, DOU emerges as sediment reactivity proxy which is the least model dependent.