Abstract. Peatlands have often been neglected in Earth system models (ESMs).
Where they are included, they are usually represented via a separate, prescribed grid cell fraction that is given the physical characteristics of a peat (highly organic) soil. However, in reality soils vary on a spectrum between purely mineral soil (no organic material) and purely organic
soil, typically with an organic layer of variable thickness overlying mineral soil below. They are also dynamic, with organic layer thickness and its properties changing over time. Neither the spectrum
of soil types nor their dynamic nature can be captured by current ESMs. Here we present a new version of an ESM land surface scheme (Joint UK Land Environment Simulator, JULES) where soil organic matter accumulation – and thus peatland formation, degradation and stability – is integrated
in the vertically resolved soil carbon scheme. We also introduce the capacity to track soil carbon age as a function of depth in JULES and compare this to measured peat age–depth profiles. The new scheme is tested and evaluated at northern and temperate sites. This scheme simulates dynamic feedbacks between the soil organic material and its thermal and hydraulic characteristics. We show that draining the peatlands can lead to significant carbon loss, soil compaction and changes in peat properties. However, negative feedbacks can lead to the potential for peatlands to rewet themselves following drainage.
These ecohydrological feedbacks can also lead to peatlands maintaining themselves in climates where peat formation would not otherwise initiate in the model, i.e. displaying some degree of resilience. The new model produces similar results to the original model for mineral soils and realistic profiles of soil organic carbon for peatlands.
We evaluate the model against typical peat profiles based on 216 northern and temperate sites from a global dataset of peat cores.
The root-mean-squared error (RMSE) in the soil carbon profile is reduced by 35 %–80 % in the best-performing JULES-Peat simulations
compared with the standard JULES configuration. The RMSE in these JULES-Peat simulations is 7.7–16.7 kg C m−3 depending on climate zone, which is considerably smaller than the soil carbon itself (around 30–60 kg C m−3). The RMSE at mineral soil sites is also reduced
in JULES-Peat compared with the original JULES configuration (reduced by ∼ 30 %–50 %). Thus, JULES-Peat can be used as a complete scheme that simulates both organic and mineral soils. It does not require
any additional input data and introduces minimal additional variables to the model. This provides a new approach for improving the simulation of organic and peatland soils and
associated carbon-cycle feedbacks in ESMs.
.

Abstract. Most land surface models (LSMs), i.e. the land components of Earth system models
(ESMs), include representation of nitrogen (N) limitation on ecosystem
productivity. However, only a few of these models have incorporated phosphorus
(P) cycling. In tropical ecosystems, this is likely to be important as N
tends to be abundant, whereas the availability of rock-derived elements, such as
P, can be very low. Thus, without a representation of P cycling, tropical
forest response in areas such as Amazonia to rising atmospheric CO2
conditions remain highly uncertain. In this study, we introduced P dynamics
and its interactions with the N and carbon (C) cycles into the Joint UK Land
Environment Simulator (JULES). The new model (JULES-CNP) includes the
representation of P stocks in vegetation and soil pools, as well as key
processes controlling fluxes between these pools. We develop and evaluate
JULES-CNP using in situ data collected at a low-fertility site in the
central Amazon, with a soil P content representative of 60 % of soils
across the Amazon basin, to parameterize, calibrate, and evaluate JULES-CNP.
Novel soil and plant P pool observations are used for parameterization and
calibration, and the model is evaluated against C fluxes and stocks and
those soil P pools not used for parameterization or calibration. We then
evaluate the model at additional P-limited test sites across the Amazon and in
Panama and Hawaii, showing a significant improvement over the C- and CN-only
versions of the model. The model is then applied under elevated
CO2 (600 ppm) at our study site in the central Amazon to quantify the impact
of P limitation on CO2 fertilization. We compare our results against the
current state-of-the-art CNP models using the same methodology that was used
in the AmazonFACE model intercomparison study. The model is able to
reproduce the observed plant and soil P pools and fluxes used for evaluation
under ambient CO2. We estimate P to limit net primary productivity
(NPP) by 24 % under current CO2 and by 46 % under elevated
CO2. Under elevated CO2, biomass in simulations accounting for CNP
increase by 10 % relative to contemporary CO2 conditions, although it
is 5 % lower compared to CN- and C-only simulations. Our results
highlight the potential for high P limitation and therefore lower CO2 fertilization capacity in the Amazon rainforest with low-fertility soils.
.

The CO2 efflux from soil (soil respiration (SR)) is one of the largest fluxes in the global carbon (C) cycle and its response to climate change could strongly influence future atmospheric CO2 concentrations. Still, a large divergence of global SR estimates and its autotrophic (AR) and heterotrophic (HR) components exists among process based terrestrial ecosystem models. Therefore, alternatively derived global benchmark values are warranted for constraining the various ecosystem model output. In this study, we developed models based on the global soil respiration database (version 5.0), using the random forest (RF) method to generate the global benchmark distribution of total SR and its components. Benchmark values were then compared with the output of ten different global terrestrial ecosystem models. Our observationally derived global mean annual benchmark rates were 85.5 ± 40.4 (SD) Pg C yr-1 for SR, 50.3 ± 25.0 (SD) Pg C yr-1 for HR and 35.2 Pg C yr-1 for AR during 1982-2012, respectively. Evaluating against the observations, the RF models showed better performance in both of SR and HR simulations than all investigated terrestrial ecosystem models. Large divergences in simulating SR and its components were observed among the terrestrial ecosystem models. The estimated global SR and HR by the ecosystem models ranged from 61.4 to 91.7 Pg C yr-1 and 39.8 to 61.7 Pg C yr-1, respectively. The most discrepancy lays in the estimation of AR, the difference (12.0-42.3 Pg C yr-1) of estimates among the ecosystem models was up to 3.5 times. The contribution of AR to SR highly varied among the ecosystem models ranging from 18% to 48%, which differed with the estimate by RF (41%). This study generated global SR and its components (HR and AR) fluxes, which are useful benchmarks to constrain the performance of terrestrial ecosystem models.

The leaching of dissolved organic carbon (DOC) from soils to the river network is an overlooked component of the terrestrial soil C budget. Measurements of DOC concentrations in soil, runoff and drainage are scarce and their spatial distribution highly skewed towards industrialized countries. The contribution of terrestrial DOC leaching to the global-scale C balance of terrestrial ecosystems thus remains poorly constrained. Here, using a process based, integrative, modelling approach to upscale from existing observations, we estimate a global terrestrial DOC leaching flux of 0.28 ± 0.07 Gt C year-1 which is conservative, as it only includes the contribution of mineral soils. Our results suggest that globally about 15% of the terrestrial Net Ecosystem Productivity (NEP, calculated as the difference between Net Primary Production and soil respiration) is exported to aquatic systems as leached DOC. In the tropical rainforest, the leached fraction of terrestrial NEP even reaches 22%. Furthermore, we simulated spatial-temporal trends in DOC leaching from soil to the river networks from 1860 to 2010. We estimated a global increase in terrestrial DOC inputs to river network of 35 Tg C year-1 (14%) from 1860 to 2010. Despite their low global contribution to the DOC leaching flux, boreal regions have the highest relative increase (28%) while tropics have the lowest relative increase (9%) over the historical period (1860s compared to 2000s). The results from our observationally constrained model approach demonstrate that DOC leaching is a significant flux in the terrestrial C budget at regional and global scales.

Current global models of the carbon (C) cycle consider only vertical gas exchanges between terrestrial or oceanic reservoirs and the atmosphere, thus not considering the lateral transport of carbon from the continents to the oceans. Therefore, those models implicitly consider all of the C which is not respired to the atmosphere to be stored on land and hence overestimate the land C sink capability. A model that represents the whole continuum from atmosphere to land and into the ocean would provide a better understanding of the Earth's C cycle and hence more reliable historical or future projections. A first and critical step in that direction is to include processes representing the production and export of dissolved organic carbon in soils. Here we present an original representation of dissolved organic C (DOC) processes in the Joint UK Land Environment Simulator (JULES-DOCM) that integrates a representation of DOC production in terrestrial ecosystems based on the incomplete decomposition of organic matter, DOC decomposition within the soil column, and DOC export to the river network via leaching. The model performance is evaluated in five specific sites for which observations of soil DOC concentration are available. Results show that the model is able to reproduce the DOC concentration and controlling processes, including leaching to the riverine system, which is fundamental for integrating terrestrial and aquatic ecosystems. Future work should include the fate of exported DOC in the river system as well as DIC and POC export from soil.