1 Introduction
Earth system models (ESMs) are essential in our need to understand and examine the global carbon (C) cycle, investigate the influence of environmental and human activities on it, and simulate possible future changes. One part of these complex models is the land-surface components that, for example, describe the C stored in soils as well as the various interactions in soils and between vegetation and the atmospheric boundary layer . Total C residing in soils has been estimated to be roughly 3000 Pg , but this estimate has large uncertainties as exact quantities are difficult to measure globally. Nevertheless, soils are the second-largest global carbon pool, after oceans. The changing climate conditions are expected to influence the processes controlling C allocation into soils and the permanence of these storage pools . A simple approach to track changes in these processes and pools is to examine the ratio of carbon isotopes in, for example, soil litter input and respiration. In order to fully utilise this approach, the underlying model must be capable of reliably simulating carbon-13 isotope (C)-related processes in the soil.
The soil C pool can be divided into different fractions based on their chemical composition, physical characteristics, or assumed turnover or residence times . Soil processes in general are complex as biological, chemical and physical drivers act simultaneously. For modelling purposes, the fate of C can be used to test alternative model formulations and to give valuable additional information of the optimal parameter values as the C signatures are sensitive indicators of changes in processes. Soil organic matter (SOM) consists of molecules with different carbon isotopes. In theory, molecules with lighter C atoms have lower activation (kinetic) energy requirements than those with C. This leads to easier decomposition of C-bearing compounds and enrichment of C in residual organic molecules . Estimating C in different fractions of SOM or varying residence times and adding CC reaction kinetics into the models would allow verification of the model functioning and improve model predictions. C-CO measurements associated with gas flux measurements provide a promising way to link soil models to ecosystem models and allow further integration to ESMs, where C isotopes are used to detect large-scale C-cycling patterns . Furthermore, the use of C as a natural tracer enables more detailed examination of underlying C-cycling processes .
Additionally, soil carbon sequestration has been a serious topic of interest for several decades as a promising method to mitigate the rising atmospheric carbon dioxide (CO) concentrations. These types of methods aim to increase the soil carbon sink by, for example, different soil tilling , crop rotation or fertilisation practises . The fundamental problem related to C sequestration is how to demonstrate that the proposed management practice and land use change increase soil C stock size, as well as under what conditions and for how long the C will remain in the soil. The quantification of small changes in soil C stocks is challenging due to large natural variability in soils and the large standing C stock. It has been estimated that the detection of 0.1 kg C m change ( %) in an agricultural field in Finland where C stock size ranges between 8.4 and 9.8 kg m in the top 30 cm requires hundreds of soil samples to be analysed . For this reason, an efficient verification system based on a combination of measurements (C stock size, CO exchange, remote sensing, etc.) and modelling is required, and a new global vision of MRV (monitoring, reporting, verification) platform was proposed by . An overlooked approach to improve MRV is to examine the C isotope compositions in the soil and in heterotrophic respiration. The reason behind the lack of these examinations is simple – such approaches require a model that can reliably represent the soil organic carbon (SOC) dynamics for different carbon isotopes while still retaining relatively straightforward structure. The latter is especially important when we take into account the lack of good-quality calibration and validation data.
In this paper we introduce a simple C isotopic circulation into the recently re-calibrated SOC model Yasso . In our approach, the decomposition of C-specific soil organic matter (C-SOM) is linearly dependent on the default Yasso model parameters, the carbon isotope fraction CC and a new scaling factor , which represents the change to the decomposition rate between the carbon isotopes. The underlying hypothesis behind this design is that since C has a larger atomic weight, it is therefore not as reactive as C, but environmental factors should still affect the decomposition of SOM, containing either isotope, similarly. We calibrate the new C-related decomposition parameters () and assess the model functionality both on a short- and long-term (100-year simulation) basis.
Our aim is to improve the Yasso20 model parameterisation to include CC reaction kinetics in the model by using empirically measured SOM and C data. The purpose of this paper is to provide a proof of concept, and we hypothesise that measuring C in soil organic matter fractions (1) detects differences in the pool C content supporting the C-fractionation and enrichment theory and (2) allows model development for future improvements in SOM decomposition predictions as C can be used as a tracer to detect changes in the underlying C decomposition processes.
2 Materials and methods
2.1 Measurements
The SOC measurements were derived from experiments described in , and , where different types of plant litter were left to decompose inside litterbags in a natural environment at Lakkasuo, a raised bog complex located at the Hyytiälä forestry field station of the University of Helsinki ( N, E; 150 m a.s.l.). We utilised data detailing the conditions for pine-branch- and pine-needle-specific litterbags. In addition to determining the initial states for both litter types, 14 litterbags describe the soil conditions for pine branches and 7 for pine needles at later stages of decomposition during the 4-year-long experiment.
The litter was characterised by dividing it into carbon fractions by sequential extractions and hydrolysis according to
1
The isotopic composition of carbon was measured on a NC2500 elemental analyser coupled to a Thermo Scientific Delta V Plus isotope ratio mass spectrometer at the Laboratory of Chronology, Finnish Museum of Natural History. The raw isotope data were normalised with a multi-point calibration using certified isotopic reference materials (USGS-40, USGS-41, IAEA-CH3 and IAEA-CH7). The mean measured (pre-normalisation) C values for calibration references were for USGS-40, for USGS-41, for IAEA-CH3 and for IAEA-CH7, with an of between measured and expected values. Replicate analyses of quality control reference materials analysed alongside the unknowns indicate a internal precision of . For the purpose of model calibration, all samples were scaled to represent the same amount of original matter (we use 1000 mass units of original matter – the given C and C values are in relation to this value).
In addition to litterbag measurements, we have utilised peat profile measurements, detailed in , to examine long-term carbon decomposition with the Yasso model. The corresponding profile is given in Table , and the age of the three topmost layers is based on radiocarbon signatures , which could not be used to accurately estimate the other two layers. We have naively set their age to 20 years each. The fractionated C content of these layers is taken directly from Table 2 in , and their relation to the Yasso fractions will be given in the results.
Table 1Peat profile used to examine long-term carbon decomposition.
Depth | Layer age |
---|---|
0–20 cm | 0–30 years |
20–26 cm | 30–40 years |
26–32 cm | 40–50 years |
32–38 cm | 50–70 years |
38–44 cm | 70–90 years |
The meteorological variables required to run the Yasso model were extracted from a nearby weather station measurements , located at Hyytiälä ( N, E; 180 m a.s.l.). We gathered monthly temperature and annual precipitation from the beginning of year 2005 to the end of 2008. Additionally we calculated averaged monthly temperature and averaged annual precipitation from years 2000–2014 to be used in simulating the long-term carbon decomposition.
2.2 Yasso modelThe soil carbon model Yasso was originally built to simulate the soil carbon stock and changes in this stock for forest soils . Previous versions of the Yasso model have also been incorporated into land-surface models; see e.g. for Yasso implementation to JSBACH – the land-surface component of the Max Planck Institutes Earth System Model (MPI-ESM). We generate the soil carbon pools utilising a recently re-calibrated version of the model called Yasso20 . This calibration utilised multiple global datasets, and we present relevant model parameter values in Appendix Table . We introduce a new C-specific decomposition into this model, which utilises new parameters (, explained below). When all , both C and C use the same standard decomposition parameters; hence, we call this model the default Yasso model and the new version, with optimised parameters, as Yasso-C13 or optimised Yasso , when we need to make a distinction. Otherwise, both models are simply referred to as the Yasso model.
Yasso decomposes litter into different pools that represent acid, water and ethanol (A, W and E) soluble matter and a lignin-like pool (N), with all the pools having different decomposition rates. The decomposed C is released back to the atmosphere as heterotrophic respiration, shifted between the AWEN pools or transferred to inactive humus (H) pool. The model is driven with monthly temperature and annual precipitation. The SOC pool decomposition in the Yasso model can be represented by the following equation.
2
The state vector (), representing the C content in AWEN pools at time , is calculated by operating the state transition matrix () on the state vector of the previous time step () and adding litter input (), which in our simulations is set to zero (as we assume no litter is added into the litterbags after the beginning). The model initial state (in our simulations) is set to match the first measurements. The matrix is a product of fluxes () between the pools and decomposition (). Pool-specific decomposition rates , presented in Eq. (), are dependent on certain parameters (), base decomposition rates (), and functions depicting the effect of woody litter diameter , precipitation () and temperature (). See for details. 3
We introduce C-SOM decomposition into the Yasso model by adding separate C-specific storages for each AWEN pool and including an additional C-specific SOM decomposition step. The input data (essentially carbon content) is first separated into C-specific content and the rest of the carbon. The Yasso model is first run normally as in Eq. (), which is followed by C decomposition using a modified version of the pool-specific base decomposition rate: 4
The new coefficients replace in Eq. (), which essentially introduces a dependency for the mass ratio of the carbon isotopes (CC) as well as a free parameter to the Yasso model for each C-AWEN pool separately. We did not include a parameter for the humus pool (H) as we did not have measurements to calibrate the related parameter. Additionally, we do not include the humus pool in the bulk carbon examinations for the same reason – bulk carbon refers to the total sum of carbon in AWEN pools.
2.3 Model calibrationWe calibrated the four parameters related to the decomposition of each AWEN pool C-SOM. The objective function () of the calibration is the cumulative squared error of the observed and modelled C values:
5 Here the summation is taken over all AWEN pools and available litterbag measurements (with measurements indicating zero concentration for total carbon content removed from the calculations). The unnormalised (pointwise) parameter likelihood is calculated as .
Since we had only four parameters to calibrate, we produced a parameter grid by fixing an increment and giving each parameter a initial value of zero. Then we added multiple new values symmetrically around the initial value with the given increment. We ran the model with every member of the parameter grid to get an estimate of the overall shape of the parameter likelihood. This process was repeated several times with grid re-centred to the point with maximum likelihood and with readjusted (smaller) increment. The results were also verified with independent runs using different initial values. All experiments were run on a eight-core laptop utilising RStudio version 1.4.1103. We used the R interface of Yasso (see code and data availability) in addition to R.utils version 2.10.1 (no other libraries were needed).
Figure 1
Presented are the pointwise parameter likelihoods (when setting the other parameters at their optimal values), divided by the maximum likelihood value. The vertical lines indicate the parameter value with the highest likelihood.
[Figure omitted. See PDF]
3 ResultsThe Yasso model calibration resulted in strictly unimodal parameter probability distributions (Fig. ). This was not unexpected as each calibrated parameter could only directly affect a single AWEN pool. The optimised parameter values are , , and (we also note that the precision of the calibration was left at the third decimal as likelihood values started to plateau). When we examine the parameter combinations with highest likelihoods (top 250 values), the strongest correlations () are present between and , , and , as well as and .
Figure 2
Shown are time series produced by the default and optimised Yasso model versions for the different AWEN pools together with assimilated observations. The model results for C and C concentrations are on top of each other.
[Figure omitted. See PDF]
The default and optimised parameter values were used to generate SOM decomposition and related C, C and C time series from the given initial states (Fig. ). The differences between the simulated C concentrations are too small to be evident (C concentrations are identical), but we get a clear signal from the C values. The default model depicts changing C, which is driven by differences in the initial isotopic fractions between the pools and the rate of the flow of carbon between the pools.
Both model versions tend to underestimate the speed of SOM decomposition (the C and C concentrations) at Lakkasuo (Fig. ) for the A and N pools and overestimate for the W pool. The bulk carbon content behaves similarly to A and N pools as these contain the most carbon. The default Yasso model is reducing the relative C content (reducing the C values) for the A and W pools and deviating from the observations, whereas the optimised model version seems to be increasing the relative C content and following the observations more closely. There is no apparent difference for the E pool, but the calibration has lowered the rate of C enrichment for the N pool. We have calculated and gathered the mean and standard deviation of the corresponding pointwise C model bias values (model observations) for the individual AWEN pools to Table .
Table 2Calculated default and optimised model bias (model observations) mean and standard deviation for the different AWEN pools.
A-C | W-C | E-C | N-C | Bulk-C | ||||||
---|---|---|---|---|---|---|---|---|---|---|
Default | 1.05 | 0.74 | 0.36 | 0.47 | 0.49 | 0.44 | ||||
Optimised | 0.76 | 0.63 | 0.37 | 0.4 | 0.44 |
Figure 3
Time series of simulated C values of the different AWEN pools for the default and optimised Yasso model versions. Scatterplot shows the peat column C values at different depths, positioned at the midpoint of assumed age, along with corresponding trend lines.
[Figure omitted. See PDF]
The Lakkasuo initial states and generated average year (averaged monthly temperature and annual precipitation) from years 2000–2014 were used to simulate a 100-year-long carbon decomposition (Fig. ). This simulation can be compared to Lakkasuo peat column C values at different depths given in Table and in . The holocellulose values are comparable to the A pool, Klason to the N pool, nonpolar extracts to the E pool and polar extracts to the W pool. A noteworthy detail is that in the short term (Fig. ) the default model increased the relative C content (C values) of E and N more than the optimised version, but on longer timescale this situation is reversed (Fig. ). The difference in C values for bulk C, between the default and optimised model versions at the end of the 100-year-long simulation, translates to approximately 0.3 ‰ difference in C mass percentage.
4 DiscussionWe have introduced simple modifications to the Yasso model in order to account for C-SOM decomposition. Incorporation of C on SOM decomposition models is a necessary step towards integration of Earth system and dynamic land ecosystem models. The C values of different organic compounds or chemical fractions of mixed organic material can be used as natural tracers which provide a unique tool to investigate and uncover complex decomposition processes in the soil. Examples of such approaches include, for example, examination of tree carbon use by chasing C-CO pulse in the respiration of leaves, whole crowns, roots, and soil ; an analysis of how stable isotope fractionation might be used to identify microbial function without incubation in soil samples ; and assessment of carbon uptake and respiration fluxes via the use of C depleted CO .
In the current study, we introduced new parameters to account for C-SOM decomposition in the Yasso model. The calibration of these parameters only depends on the C values, i.e. the relative carbon isotope concentrations, and revealed unimodal distributions for all four AWEN pool related parameters. Considering the acquired optima and taking into account that generally the ratio , then the new C-SOM decomposition utilises values that differ at maximum 3 ‰ (for ) from the default decomposition parameter values. Therefore, it is not surprising that both default and optimised model versions generate nearly identical SOM decomposition both on a short- (Fig. ) and long-term basis. Moreover, the 0.3 ‰ relative difference in bulk C between the model versions at the end of the 100-year simulation is too small to be used directly to improve bulk C estimates. The benefits come in the various ways C can be used as a natural tracer.
The acquired optima for , and are all negative, which is consistent with the initial hypothesis and C-fractionation and enrichment theory. In the model, this translates to reduced C-SOM decomposition rate (Fig. ). Likewise, the positive value for implies increased C-SOM decomposition when compared to the default model. However, the reduction in C values, when compared to the default model version, is only true on shorter timescales (Fig. ) as each pool has a trend to increase relative C content during the 100-year-long simulation (Fig. ). This is due to the reduced C-SOM decomposition in other pools – as there is more C present in these pools, there is more C available to be transported into the N pool, which compensates for the increased decomposition. The positive goes against the initial hypothesis but is a direct result of model structure. It could be worthwhile to investigate how modifications to the model, such as adding multiple soil layers with differences in the effects of temperature, precipitation or Q10 (see e.g. Fig. 1 in ), would improve the model capabilities.
The straightforward changes to the Yasso model have improved the model capabilities in reproducing observed C values in short (Fig. ) and longer timescales (Fig. ). Results from the 100-year-long simulation seem to corroborate the initial hypothesis for A, W and N pools that the relative C content in soil (larger C values) increases with time. The optimised model even yields a positive trend for E pool C, whereas the default model tends to converge the C values of all pools to roughly . The optimised model behaviour follows the trend of these measurements, and the results are highly encouraging, even though the model is driven with a single averaged year representing the meteorological conditions from the beginning of the 21st century.
Estimation and modelling of soil organic matter decomposition, but also C sequestration, are current scientific challenges. We have demonstrated how C can be implemented into a soil carbon model, so that carbon isotope signals could then be used to analyse carbon cycles in more detail and to improve model capabilities, accuracy and predictability. The required model modifications were straightforward and resulted in drastic improvement of modelled C values of SOM extracts. Although we emphasise the preliminary nature of our results due to limited calibration dataset, we foresee the model to act as a truly important tool to understand the role of isotopic fingerprints within soil carbon decomposition. The experiments demonstrated here should be viewed as a proof of concept, but further research is needed to verify the model capabilities on other sites, ecosystems and larger areas. Future research is expected to also include different management practices and croplands. Since the production of AWEN extractions with C measurements is labour-intensive, future research will likely rely on inverse calibration on larger C datasets that are rather straightforward to produce, or meta-analysis using literature-based values could be also used for further evaluation across varying scales (local, regional, global).
5 Conclusions
We have demonstrated how to incorporate C-SOM decomposition into the Yasso model and calibrate it. The model modifications were simple and straightforward and resulted in significantly improved simulated C values. The results support the initial hypothesis of C-fractionation and enrichment theory. The capability of a model to simulate soil C content and isotope-specific SOM decomposition improves the applicability of the Yasso-C13 model to scale processes from the ecosystem level to the regional and global scale using C as a tracer. Conceptually the presented work is on solid ground, but the lack of suitable calibration and validation data urges further studies with new, precise experimental C data suitable for Yasso-C13 model calibration and validation.
Appendix A Yasso model parameters
Table A1Utilised Yasso model parameter values.
Description | Parameter | Value |
---|---|---|
Base decomposition rate for A pool | 0.51 | |
Base decomposition rate for W pool | 5.19 | |
Base decomposition rate for E pool | 0.13 | |
Base decomposition rate for N pool | 0.1 | |
Mass transfer fraction from W to A | 0.5 | |
Mass transfer fraction from E to A | 0 | |
Mass transfer fraction from N to A | 1.0 | |
Mass transfer fraction from A to W | 1.0 | |
Mass transfer fraction from E to W | 0.99 | |
Mass transfer fraction from N to W | 0 | |
Mass transfer fraction from A to E | 0 | |
Mass transfer fraction from W to E | 0 | |
Mass transfer fraction from N to E | 0 | |
Mass transfer fraction from A to N | 0 | |
Mass transfer fraction from W to N | 0.163 | |
Mass transfer fraction from E to N | 0 | |
First-order temperature impact parameter for AWE | 0.158 | |
Second-order temperature impact parameter for AWE | ||
First-order temperature impact parameter for N | 0.17 | |
Second-order temperature impact parameter for N | ||
First-order temperature impact parameter for H | 0.067 | |
Second-order temperature impact parameter for H | 0 | |
Precipitation impact parameter for AWE | ||
Precipitation impact parameter for N |
Code and data availability
The data required to reproduce the results are available at Zenodo portal . The Yasso model source code and R interface are also available at Zenodo or as the “C13” branch at
Author contributions
JM was responsible for the model modifications, simulations and preparing the manuscript, with JH, JL and TV contributing to the experiment design. LA was responsible for the C analysis and HF and KK for the AWEN extractions. MO provided premises and equipment for the isotopic measurements.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
Special thanks go to Anneli Rautiainen for the work on AWEN extractions and Hanna Turunen and Igor Shevchuk for stable isotopic sample preparation and measurement.
Financial support
This research has been supported by the Strategic Research Council (grant no. 327214; MULTA), the Academy of Finland (grant no. 337552; ACCC), and the Suomalainen Tiedeakatemia (grant no. 297350; ACCUSOIL).
Review statement
This paper was edited by Ben Bond-Lamberty and reviewed by three anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2022. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Soils account for the largest share of carbon found in terrestrial ecosystems, and their status is of considerable interest for the global carbon cycle budget and atmospheric carbon concentration. The decomposition of soil organic matter depends on environmental conditions and human activities, which raises the question of how permanent are these carbon storages under changing climate. One way to get insight into carbon decomposition processes is to analyse different carbon isotope concentrations in soil organic matter. In this paper we introduce a carbon-13-isotope-specific soil organic matter decomposition add-on into the Yasso soil carbon model and assess its functionality. The new
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details




1 Finnish Meteorological Institute, P.O. Box 503, 00101 Helsinki, Finland
2 Finnish Museum of Natural History (LUOMUS), University of Helsinki, P.O. Box 64, 00014 Helsinki, Finland
3 Natural Resources Institute Finland, P.O. Box 18, 01301, Vantaa, Finland
4 Department of Microbiology and Institute for Atmospheric and Earth System Research (INAR), Faculty of Agriculture and Forestry, University of Helsinki, P.O. Box 56, 00014 Helsinki, Finland
5 Department of Forest Sciences, Faculty of Agriculture and Forestry, University of Helsinki, P.O. Box 56, 00014 Helsinki, Finland
6 Natural Resources Institute Finland, P.O. Box 2, 00791 Helsinki, Finland