Permafrost warming and potential soil carbon (SOC) release after thawing may amplify climate change, yet model estimates of present-day and future permafrost extent vary widely, partly due to ...uncertainties in simulated soil temperature. Here, we derive thermal diffusivity, a key parameter in the soil thermal regime, from depth-specific measurements of monthly soil temperature at about 200 sites in the high latitude regions. We find that, among the tested soil properties including SOC, soil texture, bulk density, and soil moisture, SOC is the dominant factor controlling the variability of diffusivity among sites. Analysis of the CMIP5 model outputs reveals that the parameterization of thermal diffusivity drives the differences in simulated present-day permafrost extent among these models. The strong SOC-thermics coupling is crucial for projecting future permafrost dynamics, since the response of soil temperature and permafrost area to a rising air temperature would be impacted by potential changes in SOC.
The high-latitude regions of the Northern Hemisphere are a nexus for the interaction between land surface physical properties and their exchange of carbon and energy with the atmosphere. At these ...latitudes, two carbon pools of planetary significance – those of the permanently frozen soils (permafrost), and of the great expanse of boreal forest – are vulnerable to destabilization in the face of currently observed climatic warming, the speed and intensity of which are expected to increase with time. Improved projections of future Arctic and boreal ecosystem transformation require improved land surface models that integrate processes specific to these cold biomes. To this end, this study lays out relevant new parameterizations in the ORCHIDEE-MICT land surface model. These describe the interactions between soil carbon, soil temperature and hydrology, and their resulting feedbacks on water and CO2 fluxes, in addition to a recently developed fire module. Outputs from ORCHIDEE-MICT, when forced by two climate input datasets, are extensively evaluated against (i) temperature gradients between the atmosphere and deep soils, (ii) the hydrological components comprising the water balance of the largest high-latitude basins, and (iii) CO2 flux and carbon stock observations. The model performance is good with respect to empirical data, despite a simulated excessive plant water stress and a positive land surface temperature bias. In addition, acute model sensitivity to the choice of input forcing data suggests that the calibration of model parameters is strongly forcing-dependent. Overall, we suggest that this new model design is at the forefront of current efforts to reliably estimate future perturbations to the high-latitude terrestrial environment.
Land surface models rarely incorporate the terrestrial phosphorus cycle and its interactions with the carbon cycle, despite the extensive scientific debate about the importance of nitrogen and ...phosphorus supply for future land carbon uptake. We describe a representation of the terrestrial phosphorus cycle for the ORCHIDEE land surface model, and evaluate it with data from nutrient manipulation experiments along a soil formation chronosequence in Hawaii. ORCHIDEE accounts for the influence of the nutritional state of vegetation on tissue nutrient concentrations, photosynthesis, plant growth, biomass allocation, biochemical (phosphatase-mediated) mineralization, and biological nitrogen fixation. Changes in the nutrient content (quality) of litter affect the carbon use efficiency of decomposition and in return the nutrient availability to vegetation. The model explicitly accounts for root zone depletion of phosphorus as a function of root phosphorus uptake and phosphorus transport from the soil to the root surface. The model captures the observed differences in the foliage stoichiometry of vegetation between an early (300-year) and a late (4.1 Myr) stage of soil development. The contrasting sensitivities of net primary productivity to the addition of either nitrogen, phosphorus, or both among sites are in general reproduced by the model. As observed, the model simulates a preferential stimulation of leaf level productivity when nitrogen stress is alleviated, while leaf level productivity and leaf area index are stimulated equally when phosphorus stress is alleviated. The nutrient use efficiencies in the model are lower than observed primarily due to biases in the nutrient content and turnover of woody biomass. We conclude that ORCHIDEE is able to reproduce the shift from nitrogen to phosphorus limited net primary productivity along the soil development chronosequence, as well as the contrasting responses of net primary productivity to nutrient addition.
Modeling of global soil organic carbon (SOC) is accompanied by large uncertainties. The heavy computational requirement limits our flexibility in disentangling uncertainty sources especially in high ...latitudes. We build a structured sensitivity analyzing framework through reorganizing the Organizing Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE)‐aMeliorated Interactions between Carbon and Temperature (MICT) model with vertically discretized SOC into one matrix equation, which brings flexibility in comprehensive sensitivity assessment. Through Sobol's method enabled by the matrix, we systematically rank 34 relevant parameters according to variance explained by each parameter and find a strong control of carbon input and turnover time on long‐term SOC storages. From further analyses for each soil layer and regional assessment, we find that the active layer depth plays a critical role in the vertical distribution of SOC and SOC equilibrium stocks in northern high latitudes (>50°N). However, the impact of active layer depth on SOC is highly interactive and nonlinear, varying across soil layers and grid cells. The stronger impact of active layer depth on SOC comes from regions with shallow active layer depth (e.g., the northernmost part of America, Asia, and some Greenland regions). The model is sensitive to the parameter that controls vertical mixing (cryoturbation rate) but only when the vertical carbon input from vegetation is limited since the effect of vertical mixing is relatively small. And the current model structure may still lack mechanisms that effectively bury nonrecalcitrant SOC. We envision a future with more comprehensive model intercomparisons and assessments with an ensemble of land carbon models adopting the matrix‐based sensitivity framework.
Key Points
One matrix equation reproduces spatial‐temporal dynamics of SOC from the vertically discretized ORCHIDEE‐MICT model
The matrix representation enables comprehensive sensitivity analyses (e.g., Sobol's method) for complex models
Active layer depth is critical for modeling SOC distribution in high latitudes
Current land surface models (LSMs) typically represent soils in a very simplistic way, assuming soil organic carbon (SOC) as a bulk, and thus impeding a correct representation of deep soil carbon ...dynamics. Moreover, LSMs generally neglect the production and export of dissolved organic carbon (DOC) from soils to rivers, leading to overestimations of the potential carbon sequestration on land. This common oversimplified processing of SOC in LSMs is partly responsible for the large uncertainty in the predictions of the soil carbon response to climate change. In this study, we present a new soil carbon module called ORCHIDEE-SOM, embedded within the land surface model ORCHIDEE, which is able to reproduce the DOC and SOC dynamics in a vertically discretized soil to 2 m. The model includes processes of biological production and consumption of SOC and DOC, DOC adsorption on and desorption from soil minerals, diffusion of SOC and DOC, and DOC transport with water through and out of the soils to rivers. We evaluated ORCHIDEE-SOM against observations of DOC concentrations and SOC stocks from four European sites with different vegetation covers: a coniferous forest, a deciduous forest, a grassland, and a cropland. The model was able to reproduce the SOC stocks along their vertical profiles at the four sites and the DOC concentrations within the range of measurements, with the exception of the DOC concentrations in the upper soil horizon at the coniferous forest. However, the model was not able to fully capture the temporal dynamics of DOC concentrations. Further model improvements should focus on a plant- and depth-dependent parameterization of the new input model parameters, such as the turnover times of DOC and the microbial carbon use efficiency. We suggest that this new soil module, when parameterized for global simulations, will improve the representation of the global carbon cycle in LSMs, thus helping to constrain the predictions of the future SOC response to global warming.
Bioenergy crop cultivation for lignocellulosic biomass is increasingly
important for future climate mitigation, and it is assumed on large scales in
integrated assessment models (IAMs) that develop ...future land use change
scenarios consistent with the dual constraint of sufficient food production
and deep decarbonization for low climate-warming targets. In most global
vegetation models, there is no specific representation of crops producing
lignocellulosic biomass, resulting in simulation biases of biomass yields and
other carbon outputs, and in turn of future bioenergy production. Here, we
introduced four new plant functional types (PFTs) to represent four major
lignocellulosic bioenergy crops, eucalypt, poplar and willow,
Miscanthus, and switchgrass, in the global process-based vegetation
model ORCHIDEE. New parameterizations of photosynthesis, carbon allocation,
and phenology are proposed based on a compilation of field measurements. A
specific harvest module is further added to the model to simulate the
rotation of bioenergy tree PFTs based on their age dynamics. The resulting
ORCHIDEE-MICT-BIOENERGY model is applied at 296 locations where field
measurements of harvested biomass are available for different bioenergy
crops. The new model can generally reproduce the global bioenergy crop yield
observations. Biases in the model results related to grid-based simulations
versus the point-scale measurements and the lack of fertilization and
fertilization management practices in the model are discussed. This study
sheds light on the importance of properly representing bioenergy crops for
simulating their yields. The parameterizations of bioenergy crops presented
here are generic enough to be applicable in other global vegetation models.
Current land surface models (LSMs) typically represent soils in a very simplistic way, assuming soil organic carbon (SOC) as a bulk, and thus impeding a correct representation of deep soil carbon ...dynamics. Moreover, LSMs generally neglect the production and export of dissolved organic carbon (DOC) from soils to rivers, leading to overestimations of the potential carbon sequestration on land. This common oversimplified processing of SOC in LSMs is partly responsible for the large uncertainty in the predictions of the soil carbon response to climate change. In this study, we present a new soil carbon module called ORCHIDEE-SOM, embedded within the land surface model ORCHIDEE, which is able to reproduce the DOC and SOC dynamics in a vertically discretized soil to 2 m. The model includes processes of biological production and consumption of SOC and DOC, DOC adsorption on and desorption from soil minerals, diffusion of SOC and DOC, and DOC transport with water through and out of the soils to rivers. We evaluated ORCHIDEE-SOM against observations of DOC concentrations and SOC stocks from four European sites with different vegetation covers: a coniferous forest, a deciduous forest, a grassland, and a cropland. The model was able to reproduce the SOC stocks along their vertical profiles at the four sites and the DOC concentrations within the range of measurements, with the exception of the DOC concentrations in the upper soil horizon at the coniferous forest. However, the model was not able to fully capture the temporal dynamics of DOC concentrations. Further model improvements should focus on a plant- and depth-dependent parameterization of the new input model parameters, such as the turnover times of DOC and the microbial carbon use efficiency. We suggest that this new soil module, when parameterized for global simulations, will improve the representation of the global carbon cycle in LSMs, thus helping to constrain the predictions of the future SOC response to global warming.