Soil Organic Carbon (SOC) Equilibrium and Model Initialisation Methods: an Application to the Rothamsted Carbon (RothC) Model

Carbon (C) emissions from anthropogenic land use have accelerated climate change. To reduce C emissions, dynamic models can be used to assess the impact of human drivers on terrestrial C sequestration. Model accuracy requires
correct initialisation, since incorrect initialisation can influence the results obtained. Therefore, we sought to improve the initialisation of a process-based SOC model, RothC, which can estimate the effect of climate and land-use change
on SOC. The most common initialisation involves running the model until equilibrium (‘spin-up run’), when the SOC pools stabilise (method 1). However, this method does not always produce realistic results. At our experimental sites, the observed SOC was not at equilibrium after 10 years, suggesting that the commonly used spin-up initialisation method assuming equilibrium might be improved. In addition to method 1, we tested two alternative initialisations for RothC that involved adjusting the total or individual SOC pool equilibrium
values by regulating the C input during the entire spin-up initialisation period (method 2) and initialising each SOC pool with recently measured SOC values obtained by SOC fractionation (method 3). Analysis of the simulation accuracy for each model initialisation, quantified using the root mean square error (RMSE), indicated that a variant of method 2 that involved adjusting the equilibrium total SOC to observed values (method 2-T) generally showed less variation in the
individual SOC pools and total SOC. Furthermore, as total SOC is the sum of all SOC pools, and because total SOC data are more readily available than the individual SOC pool data, we conclude that method 2-T is best for initialising RothC.