Using atmospheric observations to evaluate the spatiotemporal variability of CO 2 fluxes simulated by terrestrial biospheric models

Terrestrial biospheric models (TBMs) are used to extrapolate local observations and process-level understanding of land-atmosphere carbon exchange to larger regions, and serve as predictive tools for examining carbon-climate interactions. Understanding the performance of TBMs is thus crucial to the carbon cycle and climate science communities. In this study, we present and assess an approach to evaluating the spatiotemporal patterns, rather than aggregated magnitudes, of net ecosystem exchange (NEE) simulated by TBMs using atmospheric CO 2 measurements. The approach is based on statistical model selection implemented within a high-resolution atmospheric inverse model. Using synthetic data experiments, we find that current atmospheric observations are sensitive to the underlying spatiotemporal flux variability at sub-biome scales for a large portion of North America, and that atmospheric observations can therefore be used to evaluate simulated spatiotemporal flux patterns as well as to differentiate between multiple competing TBMs. Experiments using real atmospheric observations and four prototypical TBMs further confirm the applicability of the method, and demonstrate that the performance of TBMs in simulating the spatiotemporal patterns of NEE varies substantially across seasons, with best performance during the growing season and more limited skill during transition seasons. This result is consistent with previous work showing that the ability of TBMs to model flux magnitudes is also seasonally-dependent. Overall, the proposed approach provides a new avenue for evaluating TBM performance based on sub-biome-scale flux patterns, presenting an opportunity for assessing and informing model development using atmospheric observations.

[1]  S. Wofsy,et al.  Assessment of ground-based atmospheric observations for verification of greenhouse gas emissions from an urban region , 2012, Proceedings of the National Academy of Sciences.

[2]  P. Kitanidis Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[3]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[4]  A. Scott Denning,et al.  Simulated and observed fluxes of sensible and latent heat and CO2 at the WLEF‐TV tower using SiB2.5 , 2003 .

[5]  Scott D. Miller,et al.  Seasonal drought stress in the Amazon: Reconciling models and observations , 2008 .

[6]  Atul K. Jain,et al.  A model-data comparison of gross primary productivity: Results from the North American Carbon Program site synthesis , 2012 .

[7]  Philippe Ciais,et al.  Evaluation of continental carbon cycle simulations with North American flux tower observations , 2013 .

[8]  T. A. Black,et al.  A model‐data intercomparison of CO2 exchange across North America: Results from the North American Carbon Program site synthesis , 2010 .

[9]  K. Davis,et al.  Large amplitude spatial and temporal gradients in atmospheric boundary layer CO2mole fractions detected with a tower-based network in the U.S. upper Midwest , 2012 .

[10]  J. Canadell,et al.  Soil organic carbon pools in the northern circumpolar permafrost region , 2009 .

[11]  Kevin Schaefer,et al.  Global monthly averaged CO2 fluxes recovered using a geostatistical inverse modeling approach: 2. Results including auxiliary environmental data , 2008 .

[12]  J. Randerson,et al.  An atmospheric perspective on North American carbon dioxide exchange: CarbonTracker , 2007, Proceedings of the National Academy of Sciences.

[13]  G. Powell,et al.  Terrestrial Ecoregions of the World: A New Map of Life on Earth , 2001 .

[14]  J. Randerson,et al.  Interannual variability in global biomass burning emissions from 1997 to 2004 , 2006 .

[15]  F. Woodward,et al.  Terrestrial Gross Carbon Dioxide Uptake: Global Distribution and Covariation with Climate , 2010, Science.

[16]  Samuel Hammer,et al.  Regional inversion of CO2 ecosystem fluxes from atmospheric measurements: reliability of the uncertainty estimates , 2013 .

[17]  Philippe Ciais,et al.  Terrestrial biosphere model performance for inter‐annual variability of land‐atmosphere CO2 exchange , 2012 .

[18]  Jeffrey L. Anderson,et al.  Toward reliable ensemble Kalman filter estimates of CO2 fluxes , 2012 .

[19]  Y. Xue,et al.  Terrestrial biosphere models need better representation of vegetation phenology: results from the North American Carbon Program Site Synthesis , 2012 .

[20]  John C. Lin,et al.  North American CO 2 exchange: inter-comparison of modeled estimates with results from a fine-scale atmospheric inversion , 2011 .

[21]  Atul K. Jain,et al.  North American Carbon Program (NACP) regional interim synthesis: Terrestrial biospheric model intercomparison , 2012 .

[22]  D. Huntzinger,et al.  A systematic approach for comparing modeled biospheric carbon fluxes across regional scales , 2011 .

[23]  Vineet Yadav,et al.  A geostatistical synthesis study of factors affecting gross primary productivity in various ecosystems of North America , 2010 .

[24]  Anna M. Michalak,et al.  Global monthly averaged CO2 fluxes recovered using a geostatistical inverse modeling approach: 1. Results using atmospheric measurements , 2008 .

[25]  Deborah N. Huntzinger,et al.  The utility of continuous atmospheric measurements for identifying biospheric CO2 flux variability , 2011 .

[26]  David P. Turner,et al.  Atmospheric inverse modeling to constrain regional-scale CO2 budgets at high spatial and temporal resolution , 2010 .

[27]  Warren B. Cohen,et al.  Multiple constraint analysis of regional land-surface carbon flux , 2011 .

[28]  Corinne Le Quéré,et al.  An international effort to quantify regional carbon fluxes , 2011 .

[29]  Martin Heimann,et al.  Comparing Lagrangian and Eulerian models for CO 2 transport – a step towards Bayesian inverse modeling using WRF/STILT-VPRM , 2012 .

[30]  John C. Lin,et al.  A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model , 2003 .

[31]  W. Oechel,et al.  An assessment of the carbon balance of Arctic tundra: comparisons among observations, process models, and atmospheric inversions , 2011 .

[32]  Vineet Yadav,et al.  Regional-scale geostatistical inverse modeling of North American CO 2 fluxes: a synthetic data study , 2009 .

[33]  Howard E. Epstein,et al.  High stocks of soil organic carbon in the North American Arctic region , 2008 .

[34]  Jennifer A Hoeting,et al.  Model selection for geostatistical models. , 2006, Ecological applications : a publication of the Ecological Society of America.

[35]  Yaxing Wei,et al.  Reconciling estimates of the contemporary North American carbon balance among terrestrial biosphere models, atmospheric inversions, and a new approach for estimating net ecosystem exchange from inventory‐based data , 2012 .

[36]  Vineet Yadav,et al.  A backward elimination discrete optimization algorithm for model selection in spatio-temporal regression models , 2013, Environ. Model. Softw..

[37]  B. Bolker,et al.  Tundra ecosystems observed to be CO2 sources due to differential amplification of the carbon cycle. , 2013, Ecology letters.

[38]  Ramakrishna R. Nemani,et al.  Simulating terrestrial carbon fluxes using the new biosphere model “biosphere model integrating eco‐physiological and mechanistic approaches using satellite data” (BEAMS) , 2005 .

[39]  A. Michalak,et al.  Detecting fossil fuel emissions patterns from subcontinental regions using North American in situ CO2 measurements , 2014, Geophysical research letters.

[40]  William C. Skamarock,et al.  A time-split nonhydrostatic atmospheric model for weather research and forecasting applications , 2008, J. Comput. Phys..

[41]  Vineet Yadav,et al.  Attributing the variability of eddy‐covariance CO2 flux measurements across temporal scales using geostatistical regression for a mixed northern hardwood forest , 2010 .

[42]  A. Mariotti,et al.  Terrestrial mechanisms of interannual CO2 variability , 2005 .

[43]  T. E. Osterkamp,et al.  The effect of permafrost thaw on old carbon release and net carbon exchange from tundra , 2009, Nature.

[44]  S. Houweling,et al.  Emissions of CH4 and N2O over the United States and Canada based on a receptor‐oriented modeling framework and COBRA‐NA atmospheric observations , 2008 .

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

[46]  P. Tans,et al.  A geostatistical approach to surface flux estimation of atmospheric trace gases , 2004 .

[47]  B. Poulter,et al.  Evaluating the potential of large-scale simulations to predict carbon fluxes of terrestrial ecosystems over a European Eddy Covariance network , 2013 .