High-latitude ecosystems play an important role in the global carbon cycle and in regulating the climate system and are presently undergoing rapid environmental change. Accurate land cover data sets ...are required to both document these changes as well as to provide land-surface information for benchmarking and initializing Earth system models. Earth system models also require specific land cover classification systems based on plant functional types (PFTs), rather than species or ecosystems, and so post-processing of existing land cover data is often required. This study compares over Siberia, multiple land cover data sets against one another and with auxiliary data to identify key uncertainties that contribute to variability in PFT classifications that would introduce errors in Earth system modeling. Land cover classification systems from GLC 2000, GlobCover 2005 and 2009, and MODIS collections 5 and 5.1 are first aggregated to a common legend, and then compared to high-resolution land cover classification systems, vegetation continuous fields (MODIS VCFs) and satellite-derived tree heights (to discriminate against sparse, shrub, and forest vegetation). The GlobCover data set, with a lower threshold for tree cover and taller tree heights and a better spatial resolution, tends to have better distributions of tree cover compared to high-resolution data. It has therefore been chosen to build new PFT maps for the ORCHIDEE land surface model at 1 km scale. Compared to the original PFT data set, the new PFT maps based on GlobCover 2005 and an updated cross-walking approach mainly differ in the characterization of forests and degree of tree cover. The partition of grasslands and bare soils now appears more realistic compared with ground truth data. This new vegetation map provides a framework for further development of new PFTs in the ORCHIDEE model like shrubs, lichens and mosses, to represent the water and carbon cycles in northern latitudes better. Updated land cover data sets are critical for improving and maintaining the relevance of Earth system models for assessing climate and human impacts on biogeochemistry and biophysics. The new PFT map at 5 km scale is available for download from the PANGAEA website at 10.1594/PANGAEA.810709.
The carbon cycle strongly interacts with the nitrogen cycle. Several observations show that the effects of global change on primary production and carbon storage in plant biomass and soils are ...partially controlled by N availability. Nevertheless, only a small number of terrestrial biosphere models represent explicitly the nitrogen cycle, despite its importance on the carbon cycle and on climate. These models are difficult to evaluate at large spatiotemporal scales because of the scarcity of data at the global scale over a long time period. In this study, we benchmark the capacity of the O–CN global terrestrial biosphere model to reproduce temporal changes in leaf area index (LAI) at the global scale observed by NOAA_AVHRR satellites over the period 1982–2002. Using a satellite LAI product based on the normalized difference vegetation index of global inventory monitoring and modelling studies dataset, we estimate the long-term trend of LAI and we compare it with the results from the terrestrial biosphere models, either with (O–CN) or without (O–C) a dynamic nitrogen cycle coupled to the carbon–water-energy cycles. In boreal and temperate regions, including a dynamic N cycle (O–CN) improved the fit between observed and modeled temporal changes in LAI. In contrast, in the tropics, simulated LAI from the model without the dynamic N cycle (O–C) better matched observed changes in LAI over time. Despite differential regional trends, the satellite estimate suggests an increase in the global average LAI during 1982–2002 by 0.0020 m
2
m
−2
y
−1
. Both versions of the model substantially overestimated the rate of change in LAI over time (0.0065 m
2
m
−2
y
−1
for O–C and 0.0057 m
2
m
−2
y
−1
for O–CN), suggesting that some additional limitation mechanisms are missing in the model. We also estimated the relative importance of climate, CO
2
and N deposition as potential drivers of the temporal changes in LAI. We found that recent climate change better explained temporal changes in LAI when the dynamic N cycle was included in the model (higher ranked fit for O–CN vs. O–C). Using the O–C configuration to estimate the direct effect of climate on LAI, we quantified the importance of climate-N cycle feedbacks in explaining the LAI response. We found that the warming-induced release of N from soil organic matter decomposition explains 17.5 % of the global trend in LAI over time, however, reaching up to 40.9 % explained variance in the boreal zone, which is a more important contribution than increasing anthropogenic nitrogen deposition. Our analysis supports a strong connection between warming, N cycling, and vegetation productivity. These findings underscore the importance of including N cycling in global-scale models of vegetation response to environmental change.
Winter CO2 fluxes represent an important component of the annual carbon budget in northern ecosystems. Understanding winter respiration processes and their responses to climate change is also central ...to our ability to assess terrestrial carbon cycle and climate feedbacks in the future. However, the factors influencing the spatial and temporal patterns of winter ecosystem respiration (Reco) of northern ecosystems are poorly understood. For this reason, we analyzed eddy covariance flux data from 57 ecosystem sites ranging from ~35° N to ~70° N. Deciduous forests were characterized by the highest winter Reco rates (0.90 0.39 g C m−2 d−1), when winter is defined as the period during which daily air temperature remains below 0 °C. By contrast, arctic wetlands had the lowest winter Reco rates (0.02 0.02 g C m−2 d−1). Mixed forests, evergreen needle-leaved forests, grasslands, croplands and boreal wetlands were characterized by intermediate winter Reco rates (g C m−2 d−1) of 0.70(±0.33), 0.60(±0.38), 0.62(±0.43), 0.49(±0.22) and 0.27(±0.08), respectively. Our cross site analysis showed that winter air (Tair) and soil (Tsoil) temperature played a dominating role in determining the spatial patterns of winter Reco in both forest and managed ecosystems (grasslands and croplands). Besides temperature, the seasonal amplitude of the leaf area index (LAI), inferred from satellite observation, or growing season gross primary productivity, which we use here as a proxy for the amount of recent carbon available for Reco in the subsequent winter, played a marginal role in winter CO2 emissions from forest ecosystems. We found that winter Reco sensitivity to temperature variation across space (QS) was higher than the one over time (interannual, QT). This can be expected because QS not only accounts for climate gradients across sites but also for (positively correlated) the spatial variability of substrate quantity. Thus, if the models estimate future warming impacts on Reco based on QS rather than QT, this could overestimate the impact of temperature changes.
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.
Carbonyl sulfide (COS), a trace gas showing striking similarity to CO2 in terms of biochemical diffusion pathway into leaves, has been recognized as a promising indicator of the plant gross primary ...production (GPP), the amount of carbon dioxide that is absorbed through photosynthesis by terrestrial ecosystems. However, large uncertainties about the other components of its atmospheric budget prevent us from directly relating the atmospheric COS measurements to GPP. The largest uncertainty comes from the closure of its atmospheric budget, with a source component missing. Here, we explore the benefit of assimilating both COS and CO2 measurements into the LMDz atmospheric transport model to obtain consistent information on GPP, plant respiration and COS budget. To this end, we develop an analytical inverse system that optimizes biospheric fluxes for the 15 plant functional types (PFTs) defined in the ORCHIDEE global land surface model. Plant uptake of COS is parameterized as a linear function of GPP and of the leaf relative uptake (LRU), which is the ratio of COS to CO2 deposition velocities in plants. A possible scenario for the period 2008–2019 leads to a global biospheric sink of 800 GgS yr−1, with higher absorption in the high latitudes and higher oceanic emissions between 400 and 600 GgS yr−1 most of which is located in the tropics. As for the CO2 budget, the inverse system increases GPP in the high latitudes by a few GtC yr−1 without modifying the respiration compared to the ORCHIDEE fluxes used as a prior. In contrast, in the tropics the system tends to weaken both respiration and GPP. The optimized components of the COS and CO2 budgets have been evaluated against independent measurements over North America, the Pacific Ocean, at three sites in Japan and at one site in France. Overall, the posterior COS concentrations are in better agreement with the COS retrievals at 250 hPa from the MIPAS satellite and with airborne measurements made over North America and the Pacific Ocean. The system seems to have rightly corrected the underestimated GPP over the high latitudes. However, the change in seasonality of GPP in the tropics disagrees with solar-induced fluorescence (SIF) data. The decline in biospheric sink in the Amazon driven by the inversion also disagrees with MIPAS COS retrievals at 250 hPa, highlighting the lack of observational constraints in this region. Moreover, the comparison with the surface measurements in Japan and France suggests misplaced sources in the prior anthropogenic inventory, emphasizing the need for an improved inventory to better partition oceanic and continental sources in Asia and Europe.
The first satellite-based global retrievals of terrestrial sun-induced
chlorophyll fluorescence (SIF) were achieved in 2011. Since then, a number of
global SIF datasets with different spectral, ...spatial, and temporal sampling characteristics have become available to the scientific community. These
datasets have been useful to monitor the dynamics and productivity of a range
of vegetated areas worldwide, but the coarse spatiotemporal sampling and low
signal-to-noise ratio of the data hamper their application over small or
fragmented ecosystems. The recent advent of the Copernicus Sentinel-5P TROPOMI
mission and the high quality of its data products promise to alleviate this
situation, as TROPOMI provides daily global measurements at a much denser
spatial and temporal sampling than earlier satellite instruments. In this
work, we present a global SIF dataset produced from TROPOMI measurements
within the TROPOSIF project funded by the European Space Agency. The current
version of the TROPOSIF dataset covers the time period between May 2018 and
April 2021. Baseline SIF retrievals are derived from the 743–758 nm
window. A secondary SIF dataset derived from an extended fitting window
(735–758 nm window) is included. This provides an enhanced signal-to-noise
ratio at the expense of a higher sensitivity to atmospheric effects. Spectral
reflectance spectra at seven 3 nm windows devoid of atmospheric absorption
within the 665–785 nm range are also included in the TROPOSIF dataset as an
important ancillary variable to be used in combination with SIF. The
methodology to derive SIF and ancillary data as well as results from an
initial data quality assessment are presented in this work. The TROPOSIF
dataset is available through the following digital object identifier (DOI):
https://doi.org/10.5270/esa-s5p_innovation-sif-20180501_20210320-v2.1-202104
(Guanter et al., 2021).
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.
Greenland ice sheet mass loss continues to accelerate as global temperatures increase. The surface albedo of the ice sheet determines the amount of absorbed solar energy, which is a key factor in ...driving surface snow and ice melting. Satellite-retrieved snow albedo allows us to compare and optimise modelled albedo over the entirety of the ice sheet. We optimise the parameters of the albedo scheme in the ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) land surface model for 3 random years taken over the 2000–2017 period and validate over the remaining years. In particular, we want to improve the albedo at the edges of the ice sheet, since they correspond to ablation areas and show the greatest variations in runoff and surface mass balance. By giving a larger weight to points at the ice sheet's edge, we improve the model–data fit by reducing the root-mean-square deviation by over 25 % for the whole ice sheet for the summer months. This improvement is consistent for all years, even those not used in the calibration step. We also show the optimisation successfully improves the model–data fit at 87.5 % of in situ sites from the PROMICE (Programme for Monitoring of the Greenland Ice Sheet) network. We conclude by showing which additional model outputs are impacted by changes to the albedo parameters, encouraging future work using multiple data streams when optimising these parameters.
Land surface modellers need measurable proxies to
constrain the quantity of carbon dioxide (CO2) assimilated by
continental plants through photosynthesis, known as gross primary production
(GPP). ...Carbonyl sulfide (COS), which is taken up by leaves through their
stomates and then hydrolysed by photosynthetic enzymes, is a candidate GPP
proxy. A former study with the ORCHIDEE land surface model used a fixed
ratio of COS uptake to CO2 uptake normalised to respective ambient
concentrations for each vegetation type (leaf relative uptake, LRU) to
compute vegetation COS fluxes from GPP. The LRU approach is known to have
limited accuracy since the LRU ratio changes with variables such as
photosynthetically active radiation (PAR): while CO2 uptake slows under
low light, COS uptake is not light limited. However, the LRU approach has
been popular for COS–GPP proxy studies because of its ease of application
and apparent low contribution to uncertainty for regional-scale
applications. In this study we refined the COS–GPP relationship and
implemented in ORCHIDEE a mechanistic model that describes COS uptake by
continental vegetation. We compared the simulated COS fluxes against
measured hourly COS fluxes at two sites and studied the model behaviour and
links with environmental drivers. We performed simulations at a global scale,
and we estimated the global COS uptake by vegetation to be −756 Gg S yr−1,
in the middle range of former studies (−490 to −1335 Gg S yr−1). Based
on monthly mean fluxes simulated by the mechanistic approach in ORCHIDEE, we
derived new LRU values for the different vegetation types, ranging between
0.92 and 1.72, close to recently published averages for observed values of
1.21 for C4 and 1.68 for C3 plants. We transported the COS using the monthly
vegetation COS fluxes derived from both the mechanistic and the LRU
approaches, and we evaluated the simulated COS concentrations at NOAA sites.
Although the mechanistic approach was more appropriate when comparing to
high-temporal-resolution COS flux measurements, both approaches gave similar
results when transporting with monthly COS fluxes and evaluating COS
concentrations at stations. In our study, uncertainties between these two
approaches are of secondary importance compared to the uncertainties in the
COS global budget, which are currently a limiting factor to the potential of
COS concentrations to constrain GPP simulated by land surface models on the
global scale.