Investigating spatial differentiation of model parameters in a carbon cycle data assimilation system

Better estimates of the net exchange of CO(2) between the atmosphere and the terrestrial biosphere are urgently needed to improve predictions of future CO(2) levels in the atmosphere. The carbon cycle data assimilation system (CCDAS) offers the capability of inversion, while it is at the same time based on a process model that can be used independent of observational data. CCDAS allows the assimilation of atmospheric CO(2) concentrations into the terrestrial biosphere model BETHY, constraining its process parameters via an adjoint approach. Here, we investigate the effect of spatial differentiation of a universal carbon balance parameter of BETHY on posterior net CO(2) fluxes and their uncertainties. The parameter, beta, determines the characteristics of the slowly decomposing soil carbon pool and represents processes that are difficult to model explicitly. Two cases are studied with an assimilation period of 1979 to 2003. In the base case, there is a separate beta for each plant functional type (PFT). In the regionalization case, beta is differentiated not only by PFT, but also according to each of 11 large continental regions as used by the TransCom project. We find that the choice of spatial differentiation has a profound impact not only on the posterior (optimized) fluxes and their uncertainties, but even more so on the spatial covariance of the uncertainties. Differences are most pronounced in tropical regions, where observations are sparse. While regionalization leads to an improved fit to the observations by about 20% compared to the base case, we notice large spatial variations in the posterior net CO(2) flux on a grid cell level. The results illustrate the need for universal process formulations in global-scale atmospheric CO(2) inversion studies, at least as long as the observational network is too sparse to resolve spatial fluctuations at the regional scale. (Less)

[1]  Gregg Marland,et al.  Global, Regional, and National Fossil-Fuel CO2 Emissions, 1751 - 2006 (published 2009) , 2009 .

[2]  Thomas Kaminski,et al.  Propagating uncertainty through prognostic carbon cycle data assimilation system simulations , 2007 .

[3]  Taro Takahashi,et al.  Net sea-air CO2 flux over the global oceans: An improved estimate based on the sea-air pCO2 difference , 1999 .

[4]  I. C. Prentice,et al.  A dynamic global vegetation model for studies of the coupled atmosphere‐biosphere system , 2005 .

[5]  I. C. Prentice,et al.  Evaluation of ecosystem dynamics, plant geography and terrestrial carbon cycling in the LPJ dynamic global vegetation model , 2003 .

[6]  R. Giering,et al.  Two decades of terrestrial carbon fluxes from a carbon cycle data assimilation system (CCDAS) , 2005 .

[7]  R. Houghton,et al.  Carbon Flux to the Atmosphere from Land-Use Changes 1850-2005 (NDP-050) , 2008 .

[8]  M. Heimann,et al.  Climate‐induced oceanic oxygen fluxes: Implications for the contemporary carbon budget , 2002 .

[9]  Robert J. Scholes,et al.  The Carbon Cycle and Atmospheric Carbon Dioxide , 2001 .

[10]  Thomas Kaminski,et al.  Recipes for adjoint code construction , 1998, TOMS.

[11]  A. Tarantola.  Inverse problem theory : methods for data fitting and model parameter estimation , 1987 .

[12]  Casper Labuschagne,et al.  Saturation of the Southern Ocean CO2 Sink Due to Recent Climate Change , 2007, Science.

[13]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[14]  K. Davis,et al.  A Bayesian calibration of a simple carbon cycle model: The role of observations in estimating and reducing uncertainty , 2008 .

[15]  Sander Houweling,et al.  CO 2 flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport , 2003 .

[16]  J. Berry,et al.  A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species , 1980, Planta.

[17]  Roger Fletcher,et al.  A Rapidly Convergent Descent Method for Minimization , 1963, Comput. J..

[18]  Wolfgang Knorr,et al.  Is the airborne fraction of anthropogenic CO2 emissions increasing? , 2009 .

[19]  W. Press.  Numerical recipes in Fortran 77 : the art of scientific computing : volume 1 of fortran numerical recipes , 1996 .

[20]  Andrew D. Friend,et al.  A process-based, terrestrial biosphere model of ecosystem dynamics (Hybrid v3.0) , 1997 .

[21]  Kevin R. Gurney,et al.  Interannual variations in continental‐scale net carbon exchange and sensitivity to observing networks estimated from atmospheric CO2 inversions for the period 1980 to 2005 , 2008 .

[22]  Thomas Kaminski,et al.  An Example of an Automatic Differentiation-Based Modelling System , 2003, ICCSA.

[23]  Corinne Le Quéré,et al.  Regional changes in carbon dioxide fluxes of land and oceans since 1980. , 2000, Science.

[24]  Wolfgang Knorr,et al.  Annual and interannual CO2 exchanges of the terrestrial biosphere: process-based simulations and uncertainties , 2000 .

[25]  Wolfgang Knorr,et al.  Uncertainties in global terrestrial biosphere modeling: 1. A comprehensive sensitivity analysis with a new photosynthesis and energy balance scheme , 2001 .

[26]  T Kaminski,et al.  Inverse modeling of atmospheric carbon dioxide fluxes. , 2001, Science.

[27]  Martin Heimann,et al.  The global atmospheric tracer model TM3 , 1995 .

[28]  Stephen Sitch,et al.  The Carbon Balance of the Terrestrial Biosphere: Ecosystem Models and Atmospheric Observations , 2000 .

[29]  M. Scholze.  Model studies on the response of the terrestrial carbon cycle to climate change and variability , 2003 .

[30]  H. Šantrůčková,et al.  The isotopic composition of soil organic carbon on a north–south transect in western Canada , 2002 .