Attributing the variability of eddy‐covariance CO2 flux measurements across temporal scales using geostatistical regression for a mixed northern hardwood forest

[1] The relationships between terrestrial carbon dioxide flux and its primary environmental drivers are uncertain because the processes controlling CO2 cycling, especially at ecosystem scales, are not well understood. This uncertainty is compounded by the fact that the importance of controlling processes, and therefore environmental drivers, may differ across temporal scales. This paper presents and applies a geostatistical regression (GR) approach that can be used with eddy‐covariance data to investigate the relationships between carbon flux and environmental variables at multiple time scales, ranging from monthly to daily. The approach uses an adaptation of the Bayes Information Criterion to identify an optimal set of environmental variables that are able to explain the observed variability in carbon flux. In addition, GR quantifies the temporal correlation in the portion of the flux signal that cannot be explained by the selected variables and directly accounts for this correlation in the analysis. This GR approach was applied to data from the University of Michigan Biological Station (UMBS) AmeriFlux site to (i) identify the dominant explanatory variables for Net Ecosystem Exchange (NEE), Gross Ecosystem Exchange (GEE), and heterotrophic and autotrophic respiration (Rh+a )a t different temporal scales, (ii) evaluate whether environmental variables can be used to isolate the GEE and Rh+a signals from the NEE measurements, and (iii) determine the impact of temporal scale on the inferred relationships between environmental variables and CO2 flux. The results confirm the strong correlation between respiration and temperature and the influence of solar radiation on carbon uptake during the growing season. In addition, results highlight the influence of variables such as precipitation, vapor pressure deficit, and the fraction of photosynthetically active radiation (fPAR) in carbon cycling at UMBS. Many relationships between flux and auxiliary variables are found to be scale‐dependent. Site‐specific and remote‐sensing leaf area index and fPAR data are not found to be interchangeable at finer temporal scales. Results also show that a linear GR model is able to capture what may initially appear to be nonlinear relationships between flux and environmental variables, because this apparent nonlinearity is found to be explained by the covariability among key auxiliary variables. Finally, results indicate that GR can be used to identify variables that partially isolate GEE and Rh+a from the NEE signal at finer temporal scales at UMBS.

[1]  W. Oechel,et al.  On the use of MODIS EVI to assess gross primary productivity of North American ecosystems , 2006 .

[2]  H. Schmid,et al.  Controls on Annual Forest Carbon Storage: Lessons from the Past and Predictions for the Future , 2008 .

[3]  Timothy J. Robinson,et al.  Linear Models With R , 2005, Technometrics.

[4]  P. Ciais,et al.  Temperature sensitivity of decomposition in relation to soil organic matter pools: critique and outlook , 2005 .

[5]  G. Vining,et al.  Data Analysis: A Model-Comparison Approach , 1989 .

[6]  F. F. Blackman.  Optima and Limiting Factors , 1905 .

[7]  Clare E. Kazanski,et al.  Coarse woody debris and the carbon balance of a north temperate forest , 2007 .

[8]  J. Roden.  Modeling the light interception and carbon gain of individual fluttering aspen (Populus tremuloides Michx) leaves , 2003, Trees.

[9]  John Pastor,et al.  Factors Controlling Nitrogen Cycling and Nitrogen Saturation in Northern Temperate Forest Ecosystems. , 1991, Ecological applications : a publication of the Ecological Society of America.

[10]  S. Wofsy,et al.  Factors controlling CO2 exchange on timescales from hourly to decadal at Harvard Forest , 2007 .

[11]  Markus Reichstein,et al.  Modeling temporal and large‐scale spatial variability of soil respiration from soil water availability, temperature and vegetation productivity indices , 2003 .

[12]  Hans Peter Schmid,et al.  Ecosystem‐atmosphere exchange of carbon dioxide over a mixed hardwood forest in northern lower Michigan , 2003 .

[13]  Eric A. Davidson,et al.  Minimizing artifacts and biases in chamber-based measurements of soil respiration , 2002 .

[14]  E. Dorrepaal.  Are plant growth‐form‐based classifications useful in predicting northern ecosystem carbon cycling feedbacks to climate change? , 2007 .

[15]  D. Hodáňová.  An introduction to environmental biophysics , 1979, Biologia Plantarum.

[16]  W. Oechel,et al.  Environmental controls over carbon dioxide and water vapor exchange of terrestrial vegetation , 2002 .

[17]  K. Hibbard,et al.  An analysis of soil respiration across northern hemisphere temperate ecosystems , 2005 .

[18]  Hans Peter Schmid,et al.  Local-scale heterogeneity of photosynthetically active radiation (PAR), absorbed PAR and net radiation as a function of topography, sky conditions and leaf area index , 2006 .

[19]  A. Grelle,et al.  Measurement of net ecosystem exchange, productivity and respiration in three spruce forests in Sweden shows unexpectedly large soil carbon losses , 2008 .

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

[21]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[22]  H. Schmid,et al.  Respiratory carbon losses and the carbon-use efficiency of a northern hardwood forest, 1999-2003. , 2005, The New phytologist.

[23]  William W. Hargrove,et al.  New analysis reveals representativeness of the AmeriFlux network , 2003 .

[24]  A variance-ratio test for supporting a variable mean in kriging , 1997 .

[25]  Hsin-Cheng Huang,et al.  Optimal Geostatistical Model Selection , 2007 .

[26]  Maosheng Zhao,et al.  A Continuous Satellite-Derived Measure of Global Terrestrial Primary Production , 2004 .

[27]  S. Running,et al.  Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data , 2002 .

[28]  P. Kitanidis,et al.  Geostatistical interpolation of chemical concentration , 1996 .

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

[30]  Ü. Rannik,et al.  Productivity overshadows temperature in determining soil and ecosystem respiration across European forests , 2001 .

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

[32]  Noel A. C. Cressie,et al.  Statistics for Spatial Data: Cressie/Statistics , 1993 .

[33]  Yiqi Luo,et al.  Partitioning interannual variability in net ecosystem exchange between climatic variability and functional change. , 2003, Tree physiology.

[34]  W. Oechel,et al.  FLUXNET: A New Tool to Study the Temporal and Spatial Variability of Ecosystem-Scale Carbon Dioxide, Water Vapor, and Energy Flux Densities , 2001 .

[35]  K. Davis,et al.  A multi-site analysis of random error in tower-based measurements of carbon and energy fluxes , 2006 .

[36]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[37]  Michael G. Ryan,et al.  Seasonal and annual respiration of a ponderosa pine ecosystem , 1999 .

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

[39]  D. Beerling,et al.  A new technique for estimating rates of carboxylation and electron transport in leaves of C3 plants for use in dynamic global vegetation models , 1995 .

[40]  E. Ward.  A review and comparison of four commonly used Bayesian and maximum likelihood model selection tools , 2008 .

[41]  Markus Reichstein,et al.  Biosphere-atmosphere exchange of CO2 in relation to climate: a cross-biome analysis across multiple time scales , 2009 .

[42]  John L. Monteith,et al.  The Photosynthesis and Transpiration of Crops , 1966, Experimental Agriculture.

[43]  Charles T. Garten,et al.  Separating root and soil microbial contributions to soil respiration: A review of methods and observations , 2000 .

[44]  Hans Peter Schmid,et al.  Biometric and eddy-covariance based estimates of annual carbon storage in five eastern North American deciduous forests , 2002 .

[45]  E. Davidson,et al.  On the variability of respiration in terrestrial ecosystems: moving beyond Q10 , 2006 .

[46]  Allison L. Dunn,et al.  A satellite‐based biosphere parameterization for net ecosystem CO2 exchange: Vegetation Photosynthesis and Respiration Model (VPRM) , 2008 .

[47]  Martha C. Anderson,et al.  An analytical model for estimating canopy transpiration and carbon assimilation fluxes based on canopy light-use efficiency , 2000 .

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

[49]  M. Bhaskara Rao,et al.  Model Selection and Inference , 2000, Technometrics.

[50]  Patricia L. Wiberg,et al.  Relative importance of local and regional controls on coupled water, carbon, and energy fluxes , 2001 .

[51]  David R. Anderson,et al.  Model Selection and Inference: A Practical Information-Theoretic Approach , 2001 .