Evaluation of dual carbon isotope constraints (&#948;13C and &#916;14C) on the parameterisation of a mechanistic, depth explicit soil organic carbon model
Authors: Van de Broek Marijn , Govers Gerard , Schrumpf Marion , Six Johan 
Published: March 2022

&lt;p&gt;Over the past years, many mechanistic soil organic carbon models have been developed. While these models offer a clear improvement in process representation compared to classic SOC models (e.g. of microbial mineralization and mineral protection of OC), the increase in model complexity, and thus the number of uncertain model parameters, is generally not supported by sufficient available data. Thus, model parameters can generally not be sufficiently constrained through a calibration with typically used data, like soil carbon contents, because multiple, similarly good solutions are possible (referred to as the equifinality problem).&lt;/p&gt;&lt;p&gt;A common approach to reduce equifinality in SOC models is to use &lt;sup&gt;14&lt;/sup&gt;C isotopes of SOC as an additional constraint during model calibration. While this approach has shown to improve parameter identifiability (i.e., the identification of unique parameters that lead to optimal results), data on the &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C of SOC is often not available due to high measurement costs. Therefore, we explored the potential of constraining a SOC model with the more widely available ratio of stable carbon isotopes of SOC (&amp;#948;&lt;sup&gt;13&lt;/sup&gt;C). While simulations of &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C generally better constrain the turnover time of slowly cycling SOC pools, simulations of &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C allow to constrain additional processes, such as mixing of OC from aboveground and belowground sources and recycling of SOC through microbes. To do so, we developed a novel mechanistic, microbially driven model that simulates depth profiles of SOC (SOILcarb). In addition to total OC, the model simulates the &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C and &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C value of SOC, by incorporating multiple processes affecting these isotopic ratios in soils.&lt;/p&gt;&lt;p&gt;Our results show that SOILcarb was able to accurately simulate depth profiles of total SOC, &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C and &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C of a forest soil. To optimally explore the parameter space, we used a differential evolution calibration algorithm and extracted all parameter sets that led to a reliable simulation of depth profiles given a fixed error margin. The results showed that simulations calibrated only on total OC data did not results in a good fit of simulated depth profiles of either &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C or &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C. In contrast, simulations using either &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C or &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C as an additional model constraint led to accurate simulations of depth profiles of total OC, &amp;#948;&lt;sup&gt;13&lt;/sup&gt;C and &amp;#916;&lt;sup&gt;14&lt;/sup&gt;C by reducing the range in, and absolute value of, parameter values related to mainly vertical transport and protection rates of OC. Notably, parameters related to microbial OC uptake rates and microbial turnover were not better constrained by either isotopic ratio. Our results show that additional constraints on parameter values, in addition to total SOC, are necessary to increase confidence in model parameters of mechanistic SOC models, while more work is needed to better constrain microbial processes in these models.&lt;/p&gt;