Terrestrial ecosystems play a vital role in regulating the accumulation of carbon (C) in the atmosphere. Understanding the factors controlling land C uptake is critical for reducing uncertainties in ...projections of future climate. The relative importance of changing climate, rising atmospheric CO
, and other factors, however, remains unclear despite decades of research. Here, we use an ensemble of land models to show that models disagree on the primary driver of cumulative C uptake for 85% of vegetated land area. Disagreement is largest in model sensitivity to rising atmospheric CO
which shows almost twice the variability in cumulative land uptake since 1901 (1 s.d. of 212.8 PgC vs. 138.5 PgC, respectively). We find that variability in CO
and temperature sensitivity is attributable, in part, to their compensatory effects on C uptake, whereby comparable estimates of C uptake can arise by invoking different sensitivities to key environmental conditions. Conversely, divergent estimates of C uptake can occur despite being based on the same environmental sensitivities. Together, these findings imply an important limitation to the predictability of C cycling and climate under unprecedented environmental conditions. We suggest that the carbon modeling community prioritize a probabilistic multi-model approach to generate more robust C cycle projections.
Land evapotranspiration (ET) estimates are available from several global data sets. Here, monthly global land ET synthesis products, merged from these individual data sets over the time periods ...1989–1995 (7 yr) and 1989–2005 (17 yr), are presented. The merged synthesis products over the shorter period are based on a total of 40 distinct data sets while those over the longer period are based on a total of 14 data sets. In the individual data sets, ET is derived from satellite and/or in situ observations (diagnostic data sets) or calculated via land-surface models (LSMs) driven with observations-based forcing or output from atmospheric reanalyses. Statistics for four merged synthesis products are provided, one including all data sets and three including only data sets from one category each (diagnostic, LSMs, and reanalyses). The multi-annual variations of ET in the merged synthesis products display realistic responses. They are also consistent with previous findings of a global increase in ET between 1989 and 1997 (0.13 mm yr−2 in our merged product) followed by a significant decrease in this trend (−0.18 mm yr−2), although these trends are relatively small compared to the uncertainty of absolute ET values. The global mean ET from the merged synthesis products (based on all data sets) is 493 mm yr−1 (1.35 mm d−1) for both the 1989–1995 and 1989–2005 products, which is relatively low compared to previously published estimates. We estimate global runoff (precipitation minus ET) to 263 mm yr−1 (34 406 km3 yr−1) for a total land area of 130 922 000 km2. Precipitation, being an important driving factor and input to most simulated ET data sets, presents uncertainties between single data sets as large as those in the ET estimates. In order to reduce uncertainties in current ET products, improving the accuracy of the input variables, especially precipitation, as well as the parameterizations of ET, are crucial.
Quantification of global land evapotranspiration (ET) has long been associated with large uncertainties due to the lack of reference observations. Several recently developed products now provide the ...capacity to estimate ET at global scales. These products, partly based on observational data, include satellite based products, land surface model (LSM) simulations, atmospheric reanalysis output, estimates based on empirical upscaling of eddycovariance flux measurements, and atmospheric water balance datasets. The LandFlux-EVAL project aims to evaluate and compare these newly developed datasets. Additionally, an evaluation of IPCC AR4 global climate model (GCM) simulations is presented, providing an assessment of their capacity to reproduce flux behavior relative to the observations based products. Though differently constrained with observations, the analyzed reference datasets display similar large-scale ET patterns. ET from the IPCC AR4 simulations was significantly smaller than that from the other products for India (up to 1 mm/d) and parts of eastern South America, and larger in the western USA, Australia and China. The inter-product variance is lower across the IPCC AR4 simulations than across the reference datasets in several regions, which indicates that uncertainties may be underestimated in the IPCC AR4 models due to shared biases of these simulations.
Correct representation of seasonal leaf dynamics is crucial for terrestrial biosphere models (TBMs), but many such models cannot accurately reproduce observations of leaf onset and senescence. Here ...we optimised the phenology-related parameters of the ORCHIDEE TBM using satellite-derived Normalized Difference Vegetation Index data (MODIS NDVI v5) that are linearly related to the model fAPAR. We found the misfit between the observations and the model decreased after optimisation for all boreal and temperate deciduous plant functional types, primarily due to an earlier onset of leaf senescence. The model bias was only partially reduced for tropical deciduous trees and no improvement was seen for natural C4 grasses. Spatial validation demonstrated the generality of the posterior parameters for use in global simulations, with an increase in global median correlation of 0.56 to 0.67. The simulated global mean annual gross primary productivity (GPP) decreased by ~ 10 PgC yr−1 over the 1990–2010 period due to the substantially shortened growing season length (GSL – by up to 30 days in the Northern Hemisphere), thus reducing the positive bias and improving the seasonal dynamics of ORCHIDEE compared to independent data-based estimates. Finally, the optimisations led to changes in the strength and location of the trends in the simulated vegetation productivity as represented by the GSL and mean annual fraction of absorbed photosynthetically active radiation (fAPAR), suggesting care should be taken when using un-calibrated models in attribution studies. We suggest that the framework presented here can be applied for improving the phenology of all global TBMs.
Assimilation of in situ and satellite data in mechanistic terrestrial ecosystem models helps to constrain critical model parameters and reduce uncertainties in the simulated energy, water and carbon ...fluxes. So far the assimilation of eddy covariance measurements from flux-tower sites has been conducted mostly for individual sites ("single-site" optimization). Here we develop a variational data assimilation system to optimize 21 parameters of the ORCHIDEE biogeochemical model, using net CO2 flux (NEE) and latent heat flux (LE) measurements from 12 temperate deciduous broadleaf forest sites. We assess the potential of the model to simulate, with a single set of inverted parameters, the carbon and water fluxes at these 12 sites. We compare the fluxes obtained from this "multi-site" (MS) optimization to those of the prior model, and of the "single-site" (SS) optimizations. The model-data fit analysis shows that the MS approach decreases the daily root-mean-square difference (RMS) to observed data by 22%, which is close to the SS optimizations (25% on average). We also show that the MS approach distinctively improves the simulation of the ecosystem respiration (Reco), and to a lesser extent the gross primary productivity (GPP), although we only assimilated net CO2 flux. A process-oriented parameter analysis indicates that the MS inversion system finds a unique combination of parameters which is not the simple average of the different SS sets of parameters. Finally, in an attempt to validate the optimized model against independent data, we observe that global-scale simulations with MS optimized parameters show an enhanced phase agreement between modeled leaf area index (LAI) and satellite-based observations of normalized difference vegetation index (NDVI).
Remote sensing instruments are key players to map land surface temperature (LST) at large temporal and spatial scales. In this paper, we present how we combine passive microwave and thermal infrared ...data to estimate LST during summer snow-free periods over northern high latitudes. The methodology is based on the SSM/I–SSMIS 37GHz measurements at both vertical and horizontal polarizations on a 25km×25km grid size. LST is retrieved from brightness temperatures introducing an empirical linear relationship between emissivities at both polarizations as described in Royer and Poirier (2010). This relationship is calibrated at pixel scale, using cloud-free independent LST data from MODIS instruments. The SSM/I–SSMIS and MODIS data are synchronized by fitting a diurnal cycle model built on skin temperature reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). The resulting temperature dataset is provided at 25km scale and at an hourly time step during the ten-year analysis period (2000–2011). This new product was locally evaluated at five experimental sites of the EU-PAGE21 project against air temperature measurements and meteorological model reanalysis, and compared to the MODIS LST product at both local and circumpolar scale. The results giving a mean RMSE of the order of 2.2K demonstrate the usefulness of the microwave product, which is unaffected by clouds as opposed to thermal infrared products and offers a better resolution compared to model reanalysis. The dataset can be downloaded from the PANGAEA website: http://doi.pangaea.de/10.1594/PANGAEA.833409.
•Summer Arctic land surface temperature was estimated from SSM/I–SSMIS passive microwave sensors.•The product is available during snow-free periods for a ten-year period (2000–2011).•Method is based on modeling 37GHz emissivity polarization ratio at pixel scale.•Method calibration is based on the concurrent MODIS LST data.•Continuous time series at 25km resolution and hourly time step were produced.
Fire is an important global ecological process that influences the distribution of biomes, with consequences for carbon, water, and energy budgets. Therefore it is impossible to appropriately model ...the history and future of the terrestrial ecosystems and the climate system without including fire. This study incorporates the process-based prognostic fire module SPITFIRE into the global vegetation model ORCHIDEE, which was then used to simulate burned area over the 20th century. Special attention was paid to the evaluation of other fire regime indicators such as seasonality, fire size and fire length, next to burned area. For 2001–2006, the simulated global spatial extent of fire agrees well with that given by satellite-derived burned area data sets (L3JRC, GLOBCARBON, GFED3.1), and 76–92% of the global burned area is simulated as collocated between the model and observation, depending on which data set is used for comparison. The simulated global mean annual burned area is 346 Mha yr−1, which falls within the range of 287–384 Mha yr−1 as given by the three observation data sets; and is close to the 344 Mha yr−1 by the GFED3.1 data when crop fires are excluded. The simulated long-term trend and variation of burned area agree best with the observation data in regions where fire is mainly driven by climate variation, such as boreal Russia (1930–2009), along with Canada and US Alaska (1950–2009). At the global scale, the simulated decadal fire variation over the 20th century is only in moderate agreement with the historical reconstruction, possibly because of the uncertainties of past estimates, and because land-use change fires and fire suppression are not explicitly included in the model. Over the globe, the size of large fires (the 95th quantile fire size) is underestimated by the model for the regions of high fire frequency, compared with fire patch data as reconstructed from MODIS 500 m burned area data. Two case studies of fire size distribution in Canada and US Alaska, and southern Africa indicate that both number and size of large fires are underestimated, which could be related with short fire patch length and low daily fire size. Future efforts should be directed towards building consistent spatial observation data sets for key parameters of the model in order to constrain the model error at each key step of the fire modelling.
Over the last few years, solar‐induced chlorophyll fluorescence (SIF) observations from space have emerged as a promising resource for evaluating the spatio‐temporal distribution of gross primary ...productivity (GPP) simulated by global terrestrial biosphere models. SIF can be used to improve GPP simulations by optimizing critical model parameters through statistical Bayesian data assimilation techniques. A prerequisite is the availability of a functional link between GPP and SIF in terrestrial biosphere models. Here we present the development of a mechanistic SIF observation operator in the ORCHIDEE (Organizing Carbon and Hydrology In Dynamic Ecosystems) terrestrial biosphere model. It simulates the regulation of photosystem II fluorescence quantum yield at the leaf level thanks to a novel parameterization of non‐photochemical quenching as a function of temperature, photosynthetically active radiation, and normalized quantum yield of photochemistry. It emulates the radiative transfer of chlorophyll fluorescence to the top of the canopy using a parametric simplification of the SCOPE (Soil Canopy Observation Photosynthesis Energy) model. We assimilate two years of monthly OCO‐2 (Orbiting Carbon Observatory‐2) SIF product at 0.5° (2015–2016) to optimize ORCHIDEE photosynthesis and phenological parameters over an ensemble of grid points for all plant functional types. The impact on the simulated GPP is considerable with a large decrease of the global scale budget by 28 GtC/year over the period 1990–2009. The optimized GPP budget (134/136 GtC/year over 1990–2009/2001–2009) remarkably agrees with independent GPP estimates, FLUXSAT (137 GtC/year over 2001–2009) in particular and FLUXCOM (121 GtC/year over 1990–2009). Our results also suggest a biome dependency of the SIF‐GPP relationship that needs to be improved for some plant functional types.
Key Points
We developed a process‐based SIF observation operator in a terrestrial biosphere model
We assimilated monthly OCO‐2 SIF products to optimize model photosynthesis and phenology‐related parameters
The optimized GPP is considerably reduced with spatio‐temporal patterns in closer agreement with independent products
Predicting terrestrial carbon, C, budgets and carbon‐climate feedbacks strongly relies on our ability to accurately model interactions between vegetation, C and water cycles, and the atmosphere. ...However, C fluxes simulated by global, process‐based terrestrial biosphere models (TBMs) remain subject to large uncertainties, partly due to unknown or poorly calibrated parameters. This is because TBMs have not routinely been confronted against C cycle related datasets within a statistical data assimilation (DA) system. In this review, we present 15 years' development of a C cycle DA system for optimizing C cycle parameters of the ORCHIDEE TBM. We analyze the impact of assimilating multiple different C cycle related datasets on regional to global‐scale gross and net CO2 fluxes. We find that assimilating atmospheric CO2 data is crucial for improving (increasing) ORCHIDEE predictions of the terrestrial land C sink. The improvement is predominantly due to the global‐scale constraint these data provide for optimizing initial soil C stocks, which are likely in error due to inaccurate assumptions about steady state spin‐up and incomplete knowledge of land use change histories. When comparing the data‐constrained ORCHIDEE land C sink estimates to the CAMS atmospheric inversion, we show that while the two approaches agree on the global C sink magnitude, they continue to differ in how the global C sink is partitioned between the northern hemisphere and tropics. We also discuss technical challenges faced in our C cycle DA studies, in particular the difficulty in characterizing the error covariance matrix due to unknown observation biases and/or model‐data inconsistencies. We offer our perspectives on how to tackle these challenges that we hope can serve as a roadmap for other TBM groups wishing to develop C cycle DA systems.
Key Points
Considerable progress has been made in constraining ORCHIDEE terrestrial biosphere model regional to global scale CO2 fluxes within a data assimilation system
Results highlight the importance of optimizing initial C stocks ‐ in addition to C cycle related parameters ‐ using global‐scale datasets
Challenges remain in utilizing the wide variety of available data, particularly when characterizing the observation error covariance matrix
The sensitivity of global carbon and water cycling to climate variability is coupled directly to land cover and the distribution of vegetation. To investigate biogeochemistry-climate interactions, ...earth system models require a representation of vegetation distributions that are either prescribed from remote sensing data or simulated via biogeography models. However, the abstraction of earth system state variables in models means that data products derived from remote sensing need to be post-processed for model-data assimilation. Dynamic global vegetation models (DGVM) rely on the concept of plant functional types (PFT) to group shared traits of thousands of plant species into usually only 10–20 classes. Available databases of observed PFT distributions must be relevant to existing satellite sensors and their derived products, and to the present day distribution of managed lands. Here, we develop four PFT datasets based on land-cover information from three satellite sensors (EOS-MODIS 1 km and 0.5 km, SPOT4-VEGETATION 1 km, and ENVISAT-MERIS 0.3 km spatial resolution) that are merged with spatially-consistent Köppen-Geiger climate zones. Using a beta (ß) diversity metric to assess reclassification similarity, we find that the greatest uncertainty in PFT classifications occur most frequently between cropland and grassland categories, and in dryland systems between shrubland, grassland and forest categories because of differences in the minimum threshold required for forest cover. The biogeography-biogeochemistry DGVM, LPJmL, is used in diagnostic mode with the four PFT datasets prescribed to quantify the effect of land-cover uncertainty on climatic sensitivity of gross primary productivity (GPP) and transpiration fluxes. Our results show that land-cover uncertainty has large effects in arid regions, contributing up to 30% (20%) uncertainty in the sensitivity of GPP (transpiration) to precipitation. The availability of PFT datasets that are consistent with current satellite products and adapted for earth system models is an important component for reducing the uncertainty of terrestrial biogeochemistry to climate variability.