Processing arctic eddy-flux data using a simple carbon-exchange model embedded in the ensemble Kalman filter.

Continuous time-series estimates of net ecosystem carbon exchange (NEE) are routinely made using eddy covariance techniques. Identifying and compensating for errors in the NEE time series can be automated using a signal processing filter like the ensemble Kalman filter (EnKF). The EnKF compares each measurement in the time series to a model prediction and updates the NEE estimate by weighting the measurement and model prediction relative to a specified measurement error estimate and an estimate of the model-prediction error that is continuously updated based on model predictions of earlier measurements in the time series. Because of the covariance among model variables, the EnKF can also update estimates of variables for which there is no direct measurement. The resulting estimates evolve through time, enabling the EnKF to be used to estimate dynamic variables like changes in leaf phenology. The evolving estimates can also serve as a means to test the embedded model and reconcile persistent deviations between observations and model predictions. We embedded a simple arctic NEE model into the EnKF and filtered data from an eddy covariance tower located in tussock tundra on the northern foothills of the Brooks Range in northern Alaska, USA. The model predicts NEE based only on leaf area, irradiance, and temperature and has been well corroborated for all the major vegetation types in the Low Arctic using chamber-based data. This is the first application of the model to eddy covariance data. We modified the EnKF by adding an adaptive noise estimator that provides a feedback between persistent model data deviations and the noise added to the ensemble of Monte Carlo simulations in the EnKF. We also ran the EnKF with both a specified leaf-area trajectory and with the EnKF sequentially recalibrating leaf-area estimates to compensate for persistent model-data deviations. When used together, adaptive noise estimation and sequential recalibration substantially improved filter performance, but it did not improve performance when used individually. The EnKF estimates of leaf area followed the expected springtime canopy phenology. However, there were also diel fluctuations in the leaf-area estimates; these are a clear indication of a model deficiency possibly related to vapor pressure effects on canopy conductance.

[1]  A. Jazwinski.  Stochastic Processes and Filtering Theory , 1970 .

[2]  Peter C. Young,et al.  Systematic Identification of DO-BOD Model Structure , 1976 .

[3]  E. K. Webb,et al.  Correction of flux measurements for density effects due to heat and water vapour transfer , 1980 .

[4]  Robert Grover Brown,et al.  Introduction to random signal analysis and Kalman filtering , 1983 .

[5]  George M. Hornberger,et al.  Identification of photosynthesis-light models for aquatic systems II. Application to a macrophyte dominated stream , 1984 .

[6]  G. Hornberger,et al.  Identification of photosynthesis-light models for aquatic systems I. Theory and simulations , 1984 .

[7]  Peter C. Young,et al.  Recursive Estimation and Time-Series Analysis: An Introduction , 1984 .

[8]  I. E. Woodrow,et al.  A Model Predicting Stomatal Conductance and its Contribution to the Control of Photosynthesis under Different Environmental Conditions , 1987 .

[9]  D. Paslier,et al.  Net Exchange of CO2 in a Mid-Latitude Forest , 1993, Science.

[10]  E. Schulze,et al.  Leaf nitrogen, photosynthesis, conductance and transpiration : scaling from leaves to canopies , 1995 .

[11]  S. Wofsy,et al.  Modelling the soil-plant-atmosphere continuum in a Quercus-Acer stand at Harvard Forest : the regulation of stomatal conductance by light, nitrogen and soil/plant hydraulic properties , 1996 .

[12]  S. Wofsy,et al.  Physiological responses of a black spruce forest to weather , 1997 .

[13]  J. Tenhunen,et al.  Diurnal patterns of CO2 and H2O exchange of the Arctic sedges Eriophorum angustifolium and E. vaginatum (Cyperaceae). , 1998, American journal of botany.

[14]  B. Law,et al.  Carbon and water vapor exchange of an open-canopied ponderosa pine ecosystem , 1999 .

[15]  K. Hibbard,et al.  A Global Terrestrial Monitoring Network Integrating Tower Fluxes, Flask Sampling, Ecosystem Modeling and EOS Satellite Data , 1999 .

[16]  B Wang,et al.  Data assimilation and its applications. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[17]  E. Rastetter,et al.  PRIMARY PRODUCTION OF AN ARCTIC WATERSHED: AN UNCERTAINTY ANALYSIS , 2001 .

[18]  W. Oechel,et al.  Seasonality of ecosystem respiration and gross primary production as derived from FLUXNET measurements , 2001 .

[19]  J. Wilczak,et al.  Sonic Anemometer Tilt Correction Algorithms , 2001 .

[20]  S. Verma,et al.  Inversion of net ecosystem CO2 flux measurements for estimation of canopy PAR absorption , 2002 .

[21]  H. Schmid,et al.  An initial intercomparison of micrometeorological and ecological inventory estimates of carbon exchange in a mid‐latitude deciduous forest , 2002 .

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

[23]  W. Oechel,et al.  Phase and amplitude of ecosystem carbon release and uptake potentials as derived from FLUXNET measurements , 2002 .

[24]  P. Blanken,et al.  Airflows and turbulent flux measurements in mountainous terrain Part 1. Canopy and local effects , 2003 .

[25]  D. Baldocchi.  Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems: past, present and future , 2003 .

[26]  B. Law,et al.  Changes in carbon storage and fluxes in a chronosequence of ponderosa pine , 2003 .

[27]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[28]  G. Katul,et al.  Organised Motion and Radiative Perturbations in the Nocturnal Canopy Sublayer above an Even-Aged Pine Forest , 2004 .

[29]  M. Wijk,et al.  Identifying Differences in Carbon Exchange among Arctic Ecosystem Types , 2006, Ecosystems.

[30]  D. Hollinger,et al.  Uncertainty in eddy covariance measurements and its application to physiological models. , 2005, Tree physiology.

[31]  B. Law,et al.  An improved analysis of forest carbon dynamics using data assimilation , 2005 .

[32]  Niwot Ridge.  Model-data synthesis of diurnal and seasonal CO 2 fluxes at , 2006 .

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

[34]  T. Vesala,et al.  Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: algorithms and uncertainty estimation , 2006 .

[35]  J. Gove,et al.  Application of a dual unscented Kalman filter for simultaneous state and parameter estimation in problems of surface‐atmosphere exchange , 2006 .

[36]  R. Monson,et al.  Model‐data synthesis of diurnal and seasonal CO2 fluxes at Niwot Ridge, Colorado , 2006 .

[37]  E. Rastetter,et al.  Functional convergence in regulation of net CO2 flux in heterogeneous tundra landscapes in Alaska and Sweden , 2007 .

[38]  M. Bruce Beck,et al.  Understanding complex environmental systems : a dual approach , 2007 .

[39]  A. Grelle,et al.  Addressing the influence of instrument surface heat exchange on the measurements of CO2 flux from open‐path gas analyzers , 2008 .

[40]  B. Huntley,et al.  Net ecosystem exchange over heterogeneous Arctic tundra: Scaling between chamber and eddy covariance measurements , 2008 .

[41]  D. Hollinger,et al.  An improved state-parameter analysis of ecosystem models using data assimilation , 2008 .

[42]  L. White,et al.  Conditional inversion to estimate parameters from eddy-flux observations , 2009 .