The spatial distribution and fractional cover of plant functional types (PFTs) is a key uncertainty in land surface models (LSMs) that is closely linked to uncertainties in global carbon, hydrology ...and energy budgets. Land cover is considered to be an Essential Climate Variable because changes in it can result in local, regional or global scale impacts on climate. In LSMs, land cover (LC) class maps are converted to PFT fractional maps using a cross-walking (CW) table by prescribing the fraction of each PFT that occurs within each LC class. In this study we assess the largest plausible range of PFT uncertainty derived from remotely sensed LC maps produced under the European Space Agency Land Cover Climate Change Initiative on simulations of land surface fluxes using 3 leading LSMs. We evaluate the impact of uncertainty due to both LC classification algorithms, and CW procedure, on energy, moisture and carbon fluxes in LSMs. We investigate the maximum plausible range of uncertainty deriving from both LC and CW within the context of a potential biomass scale (bare ground-grass-shrub-tree), representing a gradient from low to high biomass PFTs. More specifically, plausible alternative land cover maps and associated PFT fractional distributions were produced to prioritise low or high biomass vegetation in the LC classification (uncertainty in LC), and subsequently in the assignment of PFT fractions for each LC class (uncertainty in CW), relative to a reference PFT distribution.
We examined the impact of PFT uncertainty on 3 key variables in the carbon, water and energy cycles (gross primary production (GPP), evapo-transpiration (ET), and albedo), for 3 LSMs (JSBACH, JULES and ORCHIDEE) at global scale. Results showed a greater uncertainty in PFT fraction due to CW as opposed to LC uncertainty, for all three variables. CW uncertainty in tree fraction was found to be particularly important in the northern boreal forests for simulated LSM albedo. Uncertainty in the balance between grass and bare soil fraction in arid parts of Africa, central Asia, and central Australia was also found to influence albedo and ET in all models. The spread due to PFT uncertainty for albedo was between 30 and 105% of inter-model uncertainty, for GPP between 20 and 90%, and for ET 0–30%. Each model had a different sensitivity to PFT uncertainty, for example, GPP in JSBACH was found to have a much higher sensitivity to PFT uncertainty in the tropics than JULES and ORCHIDEE, whereas the inverse was true for ET.
These results show that inter-model uncertainty for key variables in LSMs can be reduced by more accurate representation of PFT distributions. Future efforts in land cover mapping should therefore be focused on reducing CW uncertainty through better understanding of the fractional cover of PFTs within a land cover class. Efforts to reduce LC uncertainty should particularly be focused on more accurate mapping of grass and bare soil fractions in arid areas. In the context of Land Surface Models, these results demonstrate that prescribed vegetation distribution in models is a key source of uncertainty that is comparable to the spread between models for key model state variables.
•Plant functional types (PFT) mapped using land cover (LC) and cross-walking (CW).•African & South American savannahs are key to reducing carbon cycle uncertainty.•Uncertainty in bare soil & grass PFTs drives uncertainty in Europe & N America.•Northern high latitude tree cover uncertainty impacts surface albedo.•Clear need to better constrain CW fractions of PFTs using high resolution mapping.
Satellite remote sensing provides unmatched spatiotemporal information on vegetation gross primary productivity (GPP). Yet understanding of the relationship between GPP and remote sensing ...observations and how it changes with factors such as scale, biophysical constraint, and vegetation type remains limited. This knowledge gap is especially apparent for dryland ecosystems, which have characteristic high spatiotemporal variability and are under‐represented by long‐term field measurements. Here we utilize an eddy covariance (EC) data synthesis for southwestern North America in an assessment of how accurately satellite‐derived vegetation proxies capture seasonal to interannual GPP dynamics across dryland gradients. We evaluate the enhanced vegetation index, solar‐induced fluorescence (SIF), and the photochemical reflectivity index. We find evidence that SIF is more accurately capturing seasonal GPP dynamics particularly for evergreen‐dominated EC sites and more accurately estimating the full magnitude of interannual GPP dynamics for all dryland EC sites. These results suggest that incorporation of SIF could significantly improve satellite‐based GPP estimates.
Key Points
SIF captures seasonal GPP dynamics better than EVI and PRI, especially for evergreen forest sites
SIF is more sensitive than EVI and PRI to site‐level interannual GPP variability across all southwestern North America ecoregions
Incorporation of SIF could significantly improve satellite‐based GPP estimates for dryland ecosystems
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.
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
Several upcoming satellite missions have core science requirements to produce data for accurate forest aboveground biomass mapping. Largely because of these mission datasets, the number of available ...biomass products is expected to greatly increase over the coming decade. Despite the recognized importance of biomass mapping for a wide range of science, policy and management applications, there remains no community accepted standard for satellite-based biomass map validation. The Committee on Earth Observing Satellites (CEOS) is developing a protocol to fill this need in advance of the next generation of biomass-relevant satellites, and this paper presents a review of biomass validation practices from a CEOS perspective. We outline the wide range of anticipated user requirements for product accuracy assessment and provide recommendations for the validation of biomass products. These recommendations include the collection of new, high-quality in situ data and the use of airborne lidar biomass maps as tools toward transparent multi-resolution validation. Adoption of community-vetted validation standards and practices will facilitate the uptake of the next generation of biomass products.
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
Since 70 % of global forests are managed and forests impact the global carbon cycle and the energy exchange with the overlying atmosphere, forest management has the potential to mitigate climate ...change. Yet, none of the land-surface models used in Earth system models, and therefore none of today's predictions of future climate, accounts for the interactions between climate and forest management. We addressed this gap in modelling capability by developing and parametrising a version of the ORCHIDEE land-surface model to simulate the biogeochemical and biophysical effects of forest management. The most significant changes between the new branch called ORCHIDEE-CAN (SVN r2290) and the trunk version of ORCHIDEE (SVN r2243) are the allometric-based allocation of carbon to leaf, root, wood, fruit and reserve pools; the transmittance, absorbance and reflectance of radiation within the canopy; and the vertical discretisation of the energy budget calculations. In addition, conceptual changes were introduced towards a better process representation for the interaction of radiation with snow, the hydraulic architecture of plants, the representation of forest management and a numerical solution for the photosynthesis formalism of Farquhar, von Caemmerer and Berry. For consistency reasons, these changes were extensively linked throughout the code. Parametrisation was revisited after introducing 12 new parameter sets that represent specific tree species or genera rather than a group of often distantly related or even unrelated species, as is the case in widely used plant functional types. Performance of the new model was compared against the trunk and validated against independent spatially explicit data for basal area, tree height, canopy structure, gross primary production (GPP), albedo and evapotranspiration over Europe. For all tested variables, ORCHIDEE-CAN outperformed the trunk regarding its ability to reproduce large-scale spatial patterns as well as their inter-annual variability over Europe. Depending on the data stream, ORCHIDEE-CAN had a 67 to 92 % chance to reproduce the spatial and temporal variability of the validation data.
Space-borne retrievals of solar-induced chlorophyll fluorescence (SIF) over land surfaces have recently become a resource for studying and quantifying the broad scale dynamics of gross carbon uptake ...(gross primary productivity—GPP) across ecosystems. To prepare for the assimilation of SIF data in terrestrial biosphere models, we examine how differences between SIF products (due to differences in acquisition characteristics and processing chain) may affect the optimization of model parameters and the resultant GPP estimate. We compare recent daily mean SIF products (one from the Orbiting Carbon Observatory-2 OCO-2 and two from the Global Ozone Monitoring Experiment–2 GOME-2, GlobFluo GF and NASA-v28 N28, missions), averaged at 0.5° × 0.5° spatial resolution and 16-day temporal resolution, at the biome level. Phase differences between these products are relatively small. A first-order correction of the difference in spectral sampling between the two instruments shows that OCO-2 and N28 are consistent in terms of magnitude and amplitude, while GF is twice as large as the others. Using a bias-blind toy data assimilation framework, we analyze how biases between SIF products, and between model and products, can be partially alleviated by optimizing the slope and intercept parameters of a linear GPP-SIF operator. As observation biases can transfer to biases in other optimized process-based parameters and to modeled carbon fluxes— thereby resulting in unidentified inaccurate parameter values—we argue that potential SIF biases should be treated cautiously in real-world experiments in order to achieve realistic and reliable future simulations.
We investigate the benefits of assimilating in situ and satellite data of the fraction of photosynthetically active radiation (FAPAR) relative to eddy covariance flux measurements for the ...optimization of parameters of the ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystem) biosphere model. We focus on model parameters related to carbon fixation, respiration, and phenology. The study relies on two sites—Fontainebleau (deciduous broadleaf forest) and Puechabon (Mediterranean broadleaf evergreen forest)—where measurements of net carbon exchange (NEE) and latent heat (LE) fluxes are available at the same time as FAPAR products derived from ground measurements or derived from spaceborne observations at high (SPOT (Satellite Pour l′Observation de la Terre)) and medium (MERIS (MEdium Resolution Imaging Spectrometer)) spatial resolutions. We compare the different FAPAR products, analyze their consistency with the in situ fluxes, and then evaluate the potential benefits of jointly assimilating flux and FAPAR data. The assimilation of FAPAR data leads to a degradation of the model‐data agreement with respect to NEE at the two sites. It is caused by the change in leaf area required to fit the magnitude of the various FAPAR products. Assimilating daily NEE and LE fluxes, however, has a marginal impact on the simulated FAPAR. The results suggest that the main advantage of including FAPAR data is the ability to constrain the timing of leaf onset and senescence for deciduous ecosystems, which is best achieved by normalizing FAPAR time series. The joint assimilation of flux and FAPAR data leads to a model‐data improvement across all variables similar to when each data stream is used independently, corresponding, however, to different and likely improved parameter values.
Key Points
Compatibility of in situ NEE and LE data and various FAPAR products through model‐data fusion
FAPAR mainly constrains phenology; when assimilated alone, it may degrade the modeled fluxes
Combining the two data streams is preferable for improving a process‐based vegetation model
Land-use and land-cover change (LULCC) impacts local energy and
water balance and contributes on global scale to a net carbon
emission to the atmosphere. The newly released annual ESA CCI (climate ...change initiative) land
cover maps provide continuous land cover changes at 300 m
resolution from 1992 to 2015, and can be used in land surface models
(LSMs) to simulate LULCC effects on carbon stocks and on surface
energy budgets. Here we investigate the absolute areas and gross and
net changes in different plant functional types (PFTs) derived from
ESA CCI products. The results are compared with other
datasets. Global areas of forest, cropland and grassland PFTs from
ESA are 30.4, 19.3 and 35.7 million km2 in the year 2000. The
global forest area is lower than that from LUH2v2h (Hurtt et al.,
2011), Hansen et al. (2013) or Houghton and Nassikas (2017) while
cropland area is higher than LUH2v2h (Hurtt et al., 2011), in which
cropland area is from HYDE 3.2 (Klein Goldewijk et al., 2016). Gross
forest loss and gain during 1992–2015 are 1.5 and
0.9 million km2 respectively, resulting in a net forest
loss of 0.6 million km2, mainly occurring in South and
Central America. The magnitudes of gross changes in forest, cropland
and grassland PFTs in the ESA CCI are smaller than those in other
datasets. The magnitude of global net cropland gain for the whole
period is consistent with HYDE 3.2 (Klein Goldewijk et al., 2016),
but most of the increases happened before 2004 in ESA and after
2007 in HYDE 3.2. Brazil, Bolivia and Indonesia are the countries
with the largest net forest loss from 1992 to 2015, and the
decreased areas are generally consistent with those from Hansen
et al. (2013) based on Landsat 30 m resolution
images. Despite discrepancies compared to other datasets, and
uncertainties in converting into PFTs, the new ESA CCI products
provide the first detailed long-term time series of land-cover change and
can be implemented in LSMs to characterize recent carbon dynamics,
and in climate models to simulate land-cover change feedbacks on
climate. The annual ESA CCI land cover products can be downloaded
from http://maps.elie.ucl.ac.be/CCI/viewer/download.php (Land
Cover Maps – v2.0.7; see details in Sect. 5). The PFT map
translation protocol and an example in 2000 can be downloaded from
https://doi.org/10.5281/zenodo.834229. The annual ESA CCI PFT maps from
1992 to 2015 at 0.5∘×0.5∘ resolution can
also be downloaded from https://doi.org/10.5281/zenodo.1048163.