Estimation of biogeochemical rates from concentration profiles: A novel inverse method

A numerical procedure is proposed to estimate net microbiological and chemical consumption and production rates from measured solute profiles in limnic and marine sediments, as well as pelagic aquatic systems. The estimation process builds on the one-dimensional steady state diagenetic transport-reaction equation for dissolved compounds (e.g. electron acceptors, nutrients, or metabolic products). This equation is inverted by a technique called Tikhonov regularization which is a common and robust technique for solving ill-conditioned inverse problems. w procedure (REC model, Rate Estimation from Concentrations) is tested with both some artificial test cases and real measurements of sulfate and methane profiles in marine sediments. Furthermore, the technique is compared with results obtained from the models of Berg et al. (1998) and Wang et al. (2008). C model results demonstrate that the proposed numerical procedure is well suited to estimate robust and reasonable rate profiles for relative concentration measurement uncertainties not larger than ±5%. At larger uncertainties, the model provides an estimation of a nearly constant mean rate profile. w approach has a similar performance as the PROFILE model of Berg et al. (1998) and the model of Wang et al. (2008). However, the REC model has the advantage of providing smooth rate profiles and estimation of the position of discrete zones of enhanced biogeochemical activities, more accurately.