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

Author(s):  
G. J. Schürmann ◽  
T. Kaminski ◽  
C. Köstler ◽  
N. Carvalhais ◽  
M. Voßbeck ◽  
...  

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.

2016 ◽  
Vol 9 (9) ◽  
pp. 2999-3026 ◽  
Author(s):  
Gregor J. Schürmann ◽  
Thomas Kaminski ◽  
Christoph Köstler ◽  
Nuno Carvalhais ◽  
Michael Voßbeck ◽  
...  

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.


2016 ◽  
Vol 9 (9) ◽  
pp. 3321-3346 ◽  
Author(s):  
Philippe Peylin ◽  
Cédric Bacour ◽  
Natasha MacBean ◽  
Sébastien Leonard ◽  
Peter Rayner ◽  
...  

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.


2016 ◽  
Author(s):  
P. Peylin ◽  
C. Bacour ◽  
N. MacBean ◽  
S. Leonard ◽  
P. J. Rayner ◽  
...  

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.


2017 ◽  
Vol 53 (11) ◽  
pp. 8941-8965 ◽  
Author(s):  
Sujay V. Kumar ◽  
Shugong Wang ◽  
David M. Mocko ◽  
Christa D. Peters-Lidard ◽  
Youlong Xia

2008 ◽  
Vol 9 (6) ◽  
pp. 1249-1266 ◽  
Author(s):  
Jonathan L. Case ◽  
William L. Crosson ◽  
Sujay V. Kumar ◽  
William M. Lapenta ◽  
Christa D. Peters-Lidard

Abstract This manuscript presents an assessment of daily regional simulations of the Weather Research and Forecasting (WRF) numerical weather prediction (NWP) model initialized with high-resolution land surface data from the NASA Land Information System (LIS) software versus a control WRF configuration that uses land surface data from the National Centers for Environmental Prediction (NCEP) Eta Model. The goal of this study is to investigate the potential benefits of using the LIS software to improve land surface initialization for regional NWP. Fifty-eight individual nested simulations were integrated for 24 h for both the control and experimental (LISWRF) configurations during May 2004 over Florida and the surrounding areas: 29 initialized at 0000 UTC and 29 initialized at 1200 UTC. The land surface initial conditions for the LISWRF runs came from an offline integration of the Noah land surface model (LSM) within LIS for two years prior to the beginning of the month-long study on an identical grid domain to the subsequent WRF simulations. Atmospheric variables used to force the offline Noah LSM integration were provided by the North American Land Data Assimilation System and Global Data Assimilation System gridded analyses. The LISWRF soil states were generally cooler and drier than the NCEP Eta Model soil states during May 2004. Comparisons between the control and LISWRF runs for one event suggested that the LIS land surface initial conditions led to an improvement in the timing and evolution of a sea-breeze circulation over portions of northwestern Florida. Surface verification statistics for the entire month indicated that the LISWRF runs produced a more enhanced and accurate diurnal range in 2-m temperatures compared to the control as a result of the overall drier initial soil states, which resulted from a reduction in the nocturnal warm bias in conjunction with a reduction in the daytime cold bias. Daytime LISWRF 2-m dewpoints were correspondingly drier than the control dewpoints, again a manifestation of the drier initial soil states in LISWRF. The positive results of the LISWRF experiments help to illustrate the importance of initializing regional NWP models with high-quality land surface data generated at the same grid resolution.


Sign in / Sign up

Export Citation Format

Share Document