scholarly journals A new step-wise Carbon Cycle Data Assimilation System using multiple data streams to constrain the simulated land surface carbon cycle

2016 ◽  
Author(s):  
P. Peylin ◽  
C. Bacour ◽  
N. MacBean ◽  
S. Leonard ◽  
P. J. Rayner ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Carbon Cycle ◽  
Data Streams ◽  
Land Surface ◽  
Atmospheric Co2 ◽  
Model Parameters ◽  
Modis Ndvi ◽  
Optimized Model ◽  
Data Assimilation System ◽  
Assimilation System

Abstract. Large uncertainties in Land surface models (LSMs) simulations still arise from inaccurate forcing, incorrect model parameter values and incomplete representation of biogeochemical processes. The recent increase in the number and type of carbon cycle related observations, including both in situ and remote sensing measurements, has opened a new road to optimize model parameters via robust statistical model-data integration techniques, in order to reduce the simulated carbon fluxes and stocks uncertainties. In this study we present a Carbon Cycle Data Assimilation System (CCDAS) that assimilates three major data streams, namely MODIS-NDVI observations of vegetation activity, net ecosystem exchange (NEE) and latent heat (LE) flux measurements at more than 70 sites (FLUXNET), and atmospheric CO2 concentrations at 53 surface stations, in order to optimize the main parameters of the ORCHIDEE LSM (around 180 parameters in total). The system relies on a step-wise approach that assimilates each data stream in turn, propagating the information gained on the parameters from one step to the next. Overall, the ORCHIDEE model is able to achieve a consistent fit to all three data streams, which suggests that current LSMs have reached the level of development to assimilate these observations. The assimilation of MODIS-NDVI (step 1) reduced the growing season length in ORCHIDEE for temperate and boreal ecosystems, thus decreasing the global mean annual gross primary production (GPP). Using FLUXNET data (step 2) led to large improvements in the seasonal cycle of the NEE and LE fluxes for all ecosystems (i.e., increased amplitude for temperate ecosystems). The assimilation of atmospheric CO2, using the atmospheric transport model LMDz (step 3), provides an overall constraint (i.e., constraint on large scale net CO2 fluxes), resulting in an improvement of the fit to the observed atmospheric CO2 growth rate. Thus the optimized model predicts a land C sink of around 2.2 PgC.yr−1 (for the 2000–2009 period), which is more compatible with current estimates from the Global Carbon Project (GCP) than the prior value. The consistency of the step-wise approach is evaluated with backcompatibility checks. The final optimized model (after step 3) does not significantly degrade the fit to MODIS-NDVI and FLUXNET data that were assimilated in the first two steps, suggesting that a stepwise approach can be used instead of the more “challenging” implementation of a simultaneous optimization in which all data streams are assimilated together. Most parameters, including the scalar of the initial soil carbon pool size, changed during the optimization with a large error reduction. This work opens new perspectives for better predictions of the land carbon budgets.

scholarly journals A new stepwise carbon cycle data assimilation system using multiple data streams to constrain the simulated land surface carbon cycle

2016 ◽  
Vol 9 (9) ◽  
pp. 3321-3346 ◽  
Author(s):  
Philippe Peylin ◽  
Cédric Bacour ◽  
Natasha MacBean ◽  
Sébastien Leonard ◽  
Peter Rayner ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Carbon Cycle ◽  
Data Streams ◽  
Land Surface ◽  
Atmospheric Co2 ◽  
Stepwise Approach ◽  
Modis Ndvi ◽  
Optimized Model ◽  
Data Assimilation System ◽  
Assimilation System

Abstract. Large uncertainties in land surface models (LSMs) simulations still arise from inaccurate forcing, poor description of land surface heterogeneity (soil and vegetation properties), incorrect model parameter values and incomplete representation of biogeochemical processes. The recent increase in the number and type of carbon cycle-related observations, including both in situ and remote sensing measurements, has opened a new road to optimize model parameters via robust statistical model–data integration techniques, in order to reduce the uncertainties of simulated carbon fluxes and stocks. In this study we present a carbon cycle data assimilation system that assimilates three major data streams, namely the Moderate Resolution Imaging Spectroradiometer (MODIS)-Normalized Difference Vegetation Index (NDVI) observations of vegetation activity, net ecosystem exchange (NEE) and latent heat (LE) flux measurements at more than 70 sites (FLUXNET), as well as atmospheric CO2 concentrations at 53 surface stations, in order to optimize the main parameters (around 180 parameters in total) of the Organizing Carbon and Hydrology in Dynamics Ecosystems (ORCHIDEE) LSM (version 1.9.5 used for the Coupled Model Intercomparison Project Phase 5 (CMIP5) simulations). The system relies on a stepwise approach that assimilates each data stream in turn, propagating the information gained on the parameters from one step to the next. Overall, the ORCHIDEE model is able to achieve a consistent fit to all three data streams, which suggests that current LSMs have reached the level of development to assimilate these observations. The assimilation of MODIS-NDVI (step 1) reduced the growing season length in ORCHIDEE for temperate and boreal ecosystems, thus decreasing the global mean annual gross primary production (GPP). Using FLUXNET data (step 2) led to large improvements in the seasonal cycle of the NEE and LE fluxes for all ecosystems (i.e., increased amplitude for temperate ecosystems). The assimilation of atmospheric CO2, using the general circulation model (GCM) of the Laboratoire de Météorologie Dynamique (LMDz; step 3), provides an overall constraint (i.e., constraint on large-scale net CO2 fluxes), resulting in an improvement of the fit to the observed atmospheric CO2 growth rate. Thus, the optimized model predicts a land C (carbon) sink of around 2.2 PgC yr−1 (for the 2000–2009 period), which is more compatible with current estimates from the Global Carbon Project (GCP) than the prior value. The consistency of the stepwise approach is evaluated with back-compatibility checks. The final optimized model (after step 3) does not significantly degrade the fit to MODIS-NDVI and FLUXNET data that were assimilated in the first two steps, suggesting that a stepwise approach can be used instead of the more “challenging” implementation of a simultaneous optimization in which all data streams are assimilated together. Most parameters, including the scalar of the initial soil carbon pool size, changed during the optimization with a large error reduction. This work opens new perspectives for better predictions of the land carbon budgets.

scholarly journals Constraining a land surface model with multiple observations by application of the MPI-Carbon Cycle Data Assimilation System

2016 ◽  
Author(s):  
G. J. Schürmann ◽  
T. Kaminski ◽  
C. Köstler ◽  
N. Carvalhais ◽  
M. Voßbeck ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Carbon Cycle ◽  
Land Surface ◽  
Data Stream ◽  
Atmospheric Co2 ◽  
Surface Model ◽  
Terrestrial Biosphere ◽  
Multiple Observations ◽  
Data Assimilation System ◽  
Assimilation System

Abstract. We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the land surface scheme of the MPI-Earth System Model v1 (JSBACH). The simulated terrestrial biosphere processes (phenology and carbon balance) were constrained by observations of the fraction of photosynthetically active radiation (TIP-FAPAR product) and by observations of atmospheric CO2 at a global set of monitoring stations for the years 2005–2009. The system successfully, and computationally efficiently, improved average foliar area and northern extra-tropical seasonality of foliar area when constrained by TIP-FAPAR. Global net and gross carbon fluxes were improved when constrained by atmospheric CO2, although the system tended to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, therefore overall increasing the appropriateness of the resultant parameter values and biosphere dynamics. Our study thus highlights the role of the TIP-FAPAR product in stabilising the underdetermined atmospheric inversion problem and demonstrates the value of multiple-data stream assimilation for the simulation of terrestrial biosphere dynamics. The constraint on regional gross and net CO2 flux patterns is limited through the parametrisation of the biosphere model. We expect improvement on that aspect through a refined initialisation strategy and inclusion of further biosphere observations as constraints.

scholarly journals Constraining a land-surface model with multiple observations by application of the MPI-Carbon Cycle Data Assimilation System V1.0

2016 ◽  
Vol 9 (9) ◽  
pp. 2999-3026 ◽  
Author(s):  
Gregor J. Schürmann ◽  
Thomas Kaminski ◽  
Christoph Köstler ◽  
Nuno Carvalhais ◽  
Michael Voßbeck ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Carbon Cycle ◽  
Land Surface ◽  
Data Stream ◽  
Atmospheric Co2 ◽  
Surface Model ◽  
Max Planck ◽  
Multiple Observations ◽  
Data Assimilation System ◽  
Assimilation System

Abstract. We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the JSBACH land-surface scheme, which is part of the MPI-Earth System Model v1. The simulated phenology and net land carbon balance were constrained by globally distributed observations of the fraction of absorbed photosynthetically active radiation (FAPAR, using the TIP-FAPAR product) and atmospheric CO2 at a global set of monitoring stations for the years 2005 to 2009. When constrained by FAPAR observations alone, the system successfully, and computationally efficiently, improved simulated growing-season average FAPAR, as well as its seasonality in the northern extra-tropics. When constrained by atmospheric CO2 observations alone, global net and gross carbon fluxes were improved, despite a tendency of the system to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, thereby increasing the overall appropriateness of the simulated biosphere dynamics and underlying parameter values. Our study thus demonstrates the value of multiple-data-stream assimilation for the simulation of terrestrial biosphere dynamics. It further highlights the potential role of remote sensing data, here the TIP-FAPAR product, in stabilising the strongly underdetermined atmospheric inversion problem posed by atmospheric transport and CO2 observations alone. Notwithstanding these advances, the constraint of the observations on regional gross and net CO2 flux patterns on the MPI-CCDAS is limited through the coarse-scale parametrisation of the biosphere model. We expect improvement through a refined initialisation strategy and inclusion of further biosphere observations as constraints.

scholarly journals Three decades of simulated global terrestrial carbon fluxes from a data assimilation system confronted with different periods of observations

2019 ◽  
Vol 16 (15) ◽  
pp. 3009-3032 ◽  
Author(s):  
Karel Castro-Morales ◽  
Gregor Schürmann ◽  
Christoph Köstler ◽  
Christian Rödenbeck ◽  
Martin Heimann ◽  
...  
Keyword(s):  
Data Assimilation ◽  
Carbon Cycle ◽  
Atmospheric Co2 ◽  
Term Trend ◽  
Atmospheric Co2 Concentrations ◽  
Long Term Trend ◽  
Global Carbon ◽  
Data Assimilation System ◽  
Assimilation System

Abstract. During the last decade, carbon cycle data assimilation systems (CCDAS) have focused on improving the simulation of seasonal and mean global carbon fluxes over a few years by simultaneous assimilation of multiple data streams. However, the ability of a CCDAS to predict longer-term trends and variability of the global carbon cycle and the constraint provided by the observations have not yet been assessed. Here, we evaluate two near-decade-long assimilation experiments of the Max Planck Institute – Carbon Cycle Data Assimilation System (MPI-CCDAS v1) using spaceborne estimates of the fraction of absorbed photosynthetic active radiation (FAPAR) and atmospheric CO2 concentrations from the global network of flask measurement sites from either 1982 to 1990 or 1990 to 2000. We contrast these simulations with independent observations from the period 1982–2010, as well as a third MPI-CCDAS assimilation run using data from the full 1982–2010 period, and an atmospheric inversion covering the same data and time. With 30 years of data, MPI-CCDAS is capable of representing land uptake to a sufficient degree to make it compatible with the atmospheric CO2 record. The long-term trend and seasonal amplitude of atmospheric CO2 concentrations at station level over the period 1982 to 2010 is considerably improved after assimilating only the first decade (1982–1990) of observations. After 15–19 years of prognostic simulation, the simulated CO2 mixing ratio in 2007–2010 diverges by only 2±1.3 ppm from the observations, the atmospheric inversion, and the MPI-CCDAS assimilation run using observations from the full period. The long-term trend, phenological seasonality, and interannual variability (IAV) of FAPAR in the Northern Hemisphere over the last 1 to 2 decades after the assimilation were also improved. Despite imperfections in the representation of the IAV in atmospheric CO2, model–data fusion for a decade of data can already contribute to the prognostic capacity of land carbon cycle models at relevant timescales.

Sign in / Sign up

Export Citation Format

Share Document