Search

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

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 redu

The role of soil in current climate models is reviewed and discussed, with a focus on developments over the last two decades. Soil modeling may be divided into three major parts: simulation of soil hydrological dynamics, soil biogeochemistry and the soil thermal environment. Each of these three major parts is summarized with a brief description of current best practice and developments. Specific issues and modifications relevant to four extreme environments are highlighted: drylands, tropical moist and wet forests, cold regions, and peatlands and wetlands. Finally, current advances in the areas of hyperresolution and coupled model environments are discussed, which we see as the two leading edges of current soil model development. This is an update of Smith, P. (2005). Climate models, role of soil. In Daniel Hillel (ed.), Encyclopedia of soils in the environment (pp 262-268). Amsterdam: Academic Press. ISBN 9780123485304.

Microtopography can be a key driver of heterogeneity in the ground thermal and hydrological regime of permafrost landscapes. In turn, this heterogeneity can influence plant communities, methane fluxes, and the initiation of abrupt thaw processes. Here we have implemented a two-tile representation of microtopography in JULES (the Joint UK Land Environment Simulator), where tiles are representative of repeating patterns of elevation difference. Tiles are coupled by lateral flows of water, heat, and redistribution of snow, and a surface water store is added to represent ponding. Simulations are performed of two Siberian polygon sites, (Samoylov and Kytalyk) and two Scandinavian palsa sites (Stordalen and Iškoras). The model represents the observed differences between greater snow depth in hollows vs. raised areas well. The model also improves soil moisture for hollows vs. the non-tiled configuration (“standard JULES”) though the raised tile remains drier than observed. The modelled differences in snow depths and soil moisture between tiles result in the lower tile soil temperatures being warmer for palsa sites, as in reality. However, when comparing the soil temperatures for July at 20 cm depth, the difference in temperature between tiles, or “temperature splitting”, is smaller than observed (3.2 vs. 5.5 ∘C). Polygons display small (0.2 ∘C) to zero temperature splitting, in agreement with observations. Consequently, methane fluxes are near identical (+0 % to 9 %) to those for standard JULES for polygons, although they can be greater than standard JULES for palsa sites (+10 % to 49 %). Through a sensitivity analysis we quantify the relative importance of model processes with respect to soil moisture and temperatures, identifying which parameters result in the greatest uncertainty in modelled temperature. Varying the palsa elevation between 0.5 and 3 m has little effect on modelled soil temperatures, showing that using only two tiles can still be a valid representation o

Microtopography can be a key driver of heterogeneity in the ground thermal and hydrological regime of permafrost landscapes. In turn, this heterogeneity can influence plant communities, methane fluxes, and the initiation of abrupt thaw processes. Here we have implemented a two-tile representation of microtopography in JULES (the Joint UK Land Environment Simulator), where tiles are representative of repeating patterns of elevation difference. Tiles are coupled by lateral flows of water, heat, and redistribution of snow, and a surface water store is added to represent ponding. Simulations are performed of two Siberian polygon sites, (Samoylov and Kytalyk) and two Scandinavian palsa sites (Stordalen and Iškoras). The model represents the observed differences between greater snow depth in hollows vs. raised areas well. The model also improves soil moisture for hollows vs. the non-tiled configuration (“standard JULES”) though the raised tile remains drier than observed. The modelled differences in snow depths and soil moisture between tiles result in the lower tile soil temperatures being warmer for palsa sites, as in reality. However, when comparing the soil temperatures for July at 20 cm depth, the difference in temperature between tiles, or “temperature splitting”, is smaller than observed (3.2 vs. 5.5 ∘C). Polygons display small (0.2 ∘C) to zero temperature splitting, in agreement with observations. Consequently, methane fluxes are near identical (+0 % to 9 %) to those for standard JULES for polygons, although they can be greater than standard JULES for palsa sites (+10 % to 49 %). Through a sensitivity analysis we quantify the relative importance of model processes with respect to soil moisture and temperatures, identifying which parameters result in the greatest uncertainty in modelled temperature. Varying the palsa elevation between 0.5 and 3 m has little effect on modelled soil temperatures, showing that using only two tiles can still be a valid representation o

Scenarios avoiding global warming greater than 1.5 or 2 ∘C, as stipulated in the Paris Agreement, may require the combined mitigation of anthropogenic greenhouse gas emissions alongside enhancing negative emissions through approaches such as afforestation–reforestation (AR) and biomass energy with carbon capture and storage (BECCS). We use the JULES land surface model coupled to an inverted form of the IMOGEN climate emulator to investigate mitigation scenarios that achieve the 1.5 or 2 ∘C warming targets of the Paris Agreement. Specifically, within this IMOGEN-JULES framework, we focus on and characterise the global and regional effectiveness of land-based (BECCS and/or AR) and anthropogenic methane (CH4) emission mitigation, separately and in combination, on the anthropogenic fossil fuel carbon dioxide (CO2) emission budgets (AFFEBs) to 2100. We use consistent data and socio-economic assumptions from the IMAGE integrated assessment model for the second Shared Socioeconomic Pathway (SSP2). The analysis includes the effects of the methane and carbon–climate feedbacks from wetlands and permafrost thaw, which we have shown previously to be significant constraints on the AFFEBs. Globally, mitigation of anthropogenic CH4 emissions has large impacts on the anthropogenic fossil fuel emission budgets, potentially offsetting (i.e. allowing extra) carbon dioxide emissions of 188–212 Gt C. This is because of (a) the reduction in the direct and indirect radiative forcing of methane in response to the lower emissions and hence atmospheric concentration of methane and (b) carbon-cycle changes leading to increased uptake by the land and ocean by CO2-based fertilisation. Methane mitigation is beneficial everywhere, particularly for the major CH4-emitting regions of India, the USA, and China. Land-based mitigation has the potential to offset 51–100 Gt C globally, the large range reflecting assumptions and uncertainties associated with BECCS. The ranges for CH4 reduction and BECCS i

Microtopography can be a key driver of heterogeneity in the ground thermal and hydrological regime of permafrost landscapes. In turn, this heterogeneity can influence plant communities, methane fluxes and the initiation of abrupt thaw processes. Here we have implemented a two-tile representation of microtopography in JULES (the Joint UK Land Environment Simulator), where tiles are representative of repeating patterns of elevation difference. We evaluate the model against available spatially resolved observations at four sites, gauge the importance of explicitly representing microtopography for modelling methane emissions and quantify the relative importance of model processes and the model’s sensitivity its parameters. Tiles are coupled by lateral flows of water, heat and redistribution of snow. A surface water store is added to represent ponding. The model is parametrised using characteristic dimensions of landscape features at sites. Simulations are performed of two Siberian polygon sites, Samoylov and Kytalyk, and two Scandinavian palsa sites, Stordalen and Iškoras. The model represents the observed differences between greater snow depth in hollows vs raised areas well. The model also improves soil moisture for hollows vs the non-tiled configuration (‘standard JULES’) though the raised tile remains drier than observed. For the two palsa sites, it is found that drainage needs to be impeded from the lower tile, representing the non-permafrost mire, to achieve the observed soil saturation. This demonstrates the need for the landscape-scale drainage to be correctly modelled. Causes of moisture heterogeneity between tiles are decreased runoff from the low tile, differences in snowmelt, and high to low-tile water flow. Unsaturated flows between tiles are negligible, suggesting the adequacy of simpler water-table based models of lateral flow in wetland environments. The modelled differences in snow depths and soil moistures betwee

Climate change is destabilizing permafrost landscapes, affecting infrastructure, ecosystems and human livelihoods. The rate of permafrost thaw is controlled by surface and subsurface properties and processes, all of which are potentially linked with each other. Yet, no standardized protocol exists for measuring permafrost thaw and related processes and properties in a linked manner. The permafrost thaw action group of the Terrestrial Multidisciplinary distributed Observatories for the Study of the Arctic Connections (T-MOSAiC) project has developed a protocol, for use by non-specialist scientists and technicians, citizen scientists and indigenous groups, to collect standardized metadata and data on permafrost thaw. The protocol introduced here addresses the need to jointly measure permafrost thaw and the associated surface and subsurface environmental conditions. The parameters measured along transects are: snow depth, thaw depth, vegetation height, soil texture, and water level. The metadata collection includes data on timing of data collection, geographical coordinates, land surface characteristics (vegetation, ground surface, water conditions), as well as photographs. Our hope is that this openly available dataset will also be highly valuable for validation and parameterization of numerical and conceptual models, thus to the broad community represented by the T-MOSAIC project.

Methane emissions from natural wetlands tend to increase with temperature and therefore may lead to a positive feedback under future climate change. However, their temperature response includes confounding factors and appears to differ on different time scales. Observed methane emissions depend strongly on temperature on a seasonal basis, but if the annual mean emissions are compared between sites, there is only a small temperature effect. We hypothesize that microbial dynamics are a major driver of the seasonal cycle and that they can explain this apparent discrepancy. We introduce a relatively simple model of methanogenic growth and dormancy into a wetland methane scheme that is used in an Earth system model. We show that this addition is sufficient to reproduce the observed seasonal dynamics of methane emissions in fully saturated wetland sites, at the same time as reproducing the annual mean emissions. We find that a more complex scheme used in recent Earth system models does not add predictive power. The sites used span a range of climatic conditions, with the majority in high latitudes. The difference in apparent temperature sensitivity seasonally versus spatially cannot be recreated by the non-microbial schemes tested. We therefore conclude that microbial dynamics are a strong candidate to be driving the seasonal cycle of wetland methane emissions. We quantify longer-term temperature sensitivity using this scheme and show that it gives approximately a 12% increase in emissions per degree of warming globally. This is in addition to any hydrological changes, which could also impact future methane emissions.

Peatlands are wetland ecosystems with great significance as natural habitats and as major global carbon stores. They have been subject to widespread exploitation and degradation with resulting losses in characteristic biota and ecosystem functions such as climate regulation. More recently, large-scale programmes have been established to restore peatland ecosystems and the various services they provide to society. Despite significant progress in peatland science and restoration practice, we lack a process-based understanding of how soil microbiota influence peatland functioning and mediate the resilience and recovery of ecosystem services, to perturbations associated with land use and climate change. We argue that there is a need to: in the short-term, characterise peatland microbial communities across a range of spatial and temporal scales and develop an improved understanding of the links between peatland habitat, ecological functions and microbial processes; in the medium term, define what a successfully restored ‘target’ peatland microbiome looks like for key carbon cycle related ecosystem services and develop microbial-based monitoring tools for assessing restoration needs; and in the longer term, to use this knowledge to influence restoration practices and assess progress on the trajectory towards ‘intact’ peatland status. Rapid advances in genetic characterisation of the structure and functions of microbial communities offer the potential for transformative progress in these areas, but the scale and speed of methodological and conceptual advances in studying ecosystem functions is a challenge for peatland scientists. Advances in this area require multidisciplinary collaborations between peatland scientists, data scientists and microbiologists and ultimately, collaboration with the modelling community. Developing a process-based understanding of the resilience and recovery of peatlands to perturbations, such as climate extremes, fires, and drainage, will be key

Carbon cycle feedbacks represent large uncertainties in climate change projections, and the response of soil carbon to climate change contributes the greatest uncertainty to this. Future changes in soil carbon depend on changes in litter and root inputs from plants and especially on reductions in the turnover time of soil carbon (tau(s)) with warming. An approximation to the latter term for the top one metre of soil (Delta C-s,C-tau) can be diagnosed from projections made with the CMIP6 and CMIP5 Earth System Models (ESMs), and is found to span a large range even at 2 degrees C of global warming (-196 +/- 117 PgC). Here, we present a constraint on Delta C-s,C-tau, which makes use of current heterotrophic respiration and the spatial variability of tau(s) inferred from observations. This spatial emergent constraint allows us to halve the uncertainty in Delta C-s,C-tau at 2 degrees C to -232 +/- 52 PgC.