Skip to Main Content
Statistical modeling of ecosystem respiration using eddy covariance data: Maximum likelihood parameter estimation, and Monte Carlo simulation of model and parameter uncertainty, applied to three simple modelsAuthor(s): Andrew D. Richardson; David Y. Hollinger; David Y. Hollinger
Source: Agricultural and Forest Meteorology. 131: 191-208.
Publication Series: Miscellaneous Publication
PDF: Download Publication (1.77 MB)
DescriptionWhether the goal is to fill gaps in the flux record, or to extract physiological parameters from eddy covariance data, researchers are frequently interested in fitting simple models of ecosystem physiology to measured data. Presently, there is no consensus on the best models to use, or the ideal optimization criteria. We demonstrate that, given our estimates of the distribution of the stochastic uncertainty in nighttime flux measurements at the Howland (Maine, USA) AmeriFlux site, it is incorrect to fit ecosystem respiration models using ordinary least squares (OLS) optimization. Results indicate that the flux uncertainty follows a double-exponential (Laplace) distribution, and the standard deviation of the uncertainty (a@)) follows a strong seasonal pattern, increasing as an exponential function of temperature. These characteristics both violate OLS assumptions. We propose that to obtain maximum likelihood estimates of model parameters, fitting should be based on minimizing the weighted sum of the absolute deviations: C ]measured - modeledl/a(d). We examine in detail the effects of this fitting paradigm on the parameter estimates and model predictions for three simple but commonly used models of ecosystem respiration. The exponential Lloyd & Taylor model consistently provides the best fit to the measured data. Using the absolute deviation criterion reduces the estimated annud sum of respiration by about 10% (70-145 g C mP2 yP1) compared to OLS; this is comparable in magnitude but opposite in sign to the effect of filtering nighttime data using a range of plausible u* thresholds. The weighting scheme also influences the annual sum of respiration: specifying 4 6 ) as a function of air temperature consistently results in the smallest totals. However, annual sums are, in most cases, comparable (within uncertainty estimates) regardless of the model used. Monte Carlo simulations indicate that a 95% confidence interval for the annual sum of respiration is about f 20- 40 g C mP2 y-l, but varies somewhat depending on model, optimization criterion, and, most importantly, weighting scheme.
- Check the Northern Research Station web site to request a printed copy of this publication.
- Our on-line publications are scanned and captured using Adobe Acrobat.
- During the capture process some typographical errors may occur.
- Please contact Sharon Hobrla, firstname.lastname@example.org if you notice any errors which make this publication unusable.
- We recommend that you also print this page and attach it to the printout of the article, to retain the full citation information.
- This article was written and prepared by U.S. Government employees on official time, and is therefore in the public domain.
CitationRichardson, Andrew D.; Hollinger, David Y. 2005. Statistical modeling of ecosystem respiration using eddy covariance data: Maximum likelihood parameter estimation, and Monte Carlo simulation of model and parameter uncertainty, applied to three simple models. Agricultural and Forest Meteorology. 131: 191-208.
KeywordsAmeriFlux, Carbon flux, Eddy covariance, Flux uncertainty, Forest, Howland, Models, Maximum likelihood, Measurement error, Respiration
- Comparison of different objective functions for parameterization of simple respiration models
- Estimating parameters of a forest ecosystem C model with measurements of stocks and fluxes as joint constraints
- Comparing simple respiration models for eddy flux and dynamic chamber data
XML: View XML