Biogeosciences www.biogeosciences.net 1726-4170 1726-4189 5 5 2008 10.5194/bg-5-1311-2008 http://www.biogeosciences.net/5/1311/2008/ http://www.biogeosciences.net/5/1311/2008/bg-5-1311-2008.html http://www.biogeosciences.net/5/1311/2008/bg-5-1311-2008.pdf 1311 1324 2008-09-17 Influences of observation errors in eddy flux data on inverse model parameter estimation G. Lasslop M. Reichstein J. Kattge D. Papale Max-Planck Institute for Biogeochemistry, Postfach 10 01 64, 07701 Jena, Germany Department of Forest Environment and Resources, DISAFRI, Univ. of Tuscia, Viterbo, Italy Eddy covariance data are increasingly used to estimate parameters of ecosystem models. For proper maximum likelihood parameter estimates the error structure in the observed data has to be fully characterized. In this study we propose a method to characterize the random error of the eddy covariance flux data, and analyse error distribution, standard deviation, cross- and autocorrelation of CO₂ and H₂O flux errors at four different European eddy covariance flux sites. Moreover, we examine how the treatment of those errors and additional systematic errors influence statistical estimates of parameters and their associated uncertainties with three models of increasing complexity – a hyperbolic light response curve, a light response curve coupled to water fluxes and the SVAT scheme BETHY. In agreement with previous studies we find that the error standard deviation scales with the flux magnitude. The previously found strongly leptokurtic error distribution is revealed to be largely due to a superposition of almost Gaussian distributions with standard deviations varying by flux magnitude. The crosscorrelations of CO₂ and H₂O fluxes were in all cases negligible (<i>R</i>² below 0.2), while the autocorrelation is usually below 0.6 at a lag of 0.5 h and decays rapidly at larger time lags. This implies that in these cases the weighted least squares criterion yields maximum likelihood estimates. To study the influence of the observation errors on model parameter estimates we used synthetic datasets, based on observations of two different sites. We first fitted the respective models to observations and then added the random error estimates described above and the systematic error, respectively, to the model output. This strategy enables us to compare the estimated parameters with true parameters. We illustrate that the correct implementation of the random error standard deviation scaling with flux magnitude significantly reduces the parameter uncertainty and often yields parameter retrievals that are closer to the true value, than by using ordinary least squares. The systematic error leads to systematically biased parameter estimates, but its impact varies by parameter. The parameter uncertainty slightly increases, but the true parameter is not within the uncertainty range of the estimate. This means that the uncertainty is underestimated with current approaches that neglect selective systematic errors in flux data. Hence, we conclude that potential systematic errors in flux data need to be addressed more thoroughly in data assimilation approaches since otherwise uncertainties will be vastly underestimated. Baldocchi, D D., Falge, E., Gu, L., Olson, R., Hollinger, D Y., Running, S W., Anthoni, P., Bernhofer, C., Davis, K J., Evans, R., Fuentes, J., Goldstein, A., Katul, G., Law, B E., Lee, X., Malhi, Y., Meyers, T P., Munger, J W., Oechel, W C., Paw~U, K T., Pilegaard, K., Schmid, H P., Valentini, R., Verma, S., Vesala, T., Wilson, K B., and Wofsy, S C.: FLUXNET: A New Tool to Study the Temporal and Spatial Variability of Ecosystem-Scale Carbon Dioxide, Water Vapor, and Energy Flux Densities, B. Am. Meteorol. Soc., 82, 2415–2434, 2001. Beer, C., Reichstein, M., Ciais, P., Farquhar, G D., and Papale, D.: Mean annual GPP of Europe derived from its water balance, Geophys. Res. Lett., 34, L05401, doi:10.1029/2006GL029006, 2007. Braswell, B., Sacks, W., Linder, E., and Schimel, D.: Estimating diurnal to annual ecosystem parameters by synthesis of a carbon flux model with eddy covariance net ecosystem exchange observations., Glob. Change Biol., 11, 335–355, 2005. Carvalhais, N., Reichstein, M., Seixas, J., Collatz, G.J., Santos Pereira, J., Berbigier, P. Carrara, A., Granier, A., Montagnani, L., Papale, D., Rambal, S., Sanz, M. J., Valentini, R.: Implications of Carbon Cycle Steady State Assumptions for Biogeochemical Modeling Performance and Inverse Parameter Retrieval, Global Biogeochem. Cy., 22, doi:10.1029/2007GB003033, 2008. Draper, N. and Smith, H.: Applied Regression Analysis, Wiley, New York, third edition, 736 pp., 1998. Evans, G T.: Defining misfit between biogeochemical models and data sets, J. Marine Syst., 40–41, 49–54, 2003. Farquhar, G D., von Caemmerer, S., and Berry, J A.: A biochemical model of photosynthesis in leaves of C3 species, Planta, 149, 78-90, 1980. Friend, A D., Arneth, A., Kiang, N Y., Lomas, M., Ogee, J., Rodenbeck, C., Running, S W., Santaren, J.-D., Sitch, S., Viovy, N., Woodwards, F I., and Zaehle, S.: FLUXNET and modelling the global carbon cycle, Glob. Change Biol., 13, 610–633, 2007. Hollinger, D. and Richardson, A.: Uncertainty in eddy covariance measurements and its application to physiological models, Tree Physiol., 25, 873–885, 2005. Hollinger, D Y., Aber, J., Dail, B., Davidson, E A., Goltz, S M., Hughes, H., Leclerc, M Y., Lee, J T., Richardson, A D., Rodrigues, C., Scott, N., Achuatavarier, D., and Walsh, J.: Spatial and temporal variability in forest-atmosphere CO₂ exchange, Glob. Change Biol., 10, 1689–1706, 2004. Knorr, W.: Annual and interannual CO₂ exchanges of the terrestrial biosphere: process-based simulations and uncertainties, Global Ecol. Biogeogr., 9, 225-252, 2000. Knorr, W. and Heimann, M.: Uncertainties in global terrestrial biosphere modelingUncertainties in global terrestrial biosphere modeling. Part I: a comprehensive sensitivity analysis with a new photosynthesis and energy balance scheme, Global Biogeochem. Cy., 15, 207-225, 2001. Knorr, W. and Kattge, J.: Inversion of terrestrial ecosystem model parameter values against eddy covariance measurements by Monte Carlo sampling, Glob. Change Biol., 11, 1333–1351, 2005. Moffat, A M., Papale, D., Reichstein, M., Hollinger, D Y., Richardson, A D., Barr, A G., Beckstein, C., Braswell, B H., Churkina, G., Desai, A R., Falge, E., Gove, J H., Heimann, M., Hui, D., Jarvis, A J., Kattge, J., Noormets, A., and Stauch, V J.: Comprehensive comparison of gap-filling techniques for eddy covariance net carbon fluxes, Agr. Forest Meteorol., 147, 209–232, 2007. Moncrieff, J., Malhi, Y., and Leuning, R.: The propagation of errors in long-term measurements of land-atmosphere fluxes of carbon and water, Glob. Change Biol., 2, 231–240, 1996. Monteith, J L.: Evaporation and environment, Symposium of the Society for Experimental Biology, 19, 205-234, 1965. Owen, K E., Tenhunen, J., Reichstein, M., Wang, Q., Falge, E., Geyer, R., Xiao, X., Stoy, P., Ammann, C., Arain, A., Aubinet, M., Aurela, M., Bernhofer, C., Chojnicki, B H., Granier, A., Gruenwald, T., Hadley, J., Heinesch, B., Hollinger, D., Knohl, A., Kutsch, W., Lohila, A., Meyers, T., Moors, E., Moureaux, C., Pilegaard, K., Saigusa, N., Verma, S., Vesala, T., and Vogel, C.: Linking flux network measurements to continental scale simulations: ecosystem carbon dioxide exchange capacity under non-water-stressed conditions, Glob. Change Biol., 13, 734–760, 2007. Papale, D., Reichstein, M., Aubinet, M., Canfora, E., Bernhofer, C., Kutsch, W., Longdoz, B., Rambal, S., Valentini, R., Vesala, T., and Yakir, D.: Towards a standardized processing of Net Ecosystem Exchange measured with eddy covariance technique: algorithms and uncertainty estimation, Biogeosciences, 3, 571–583, 2006. Raupach, M R., Rayner, P J., Barrett, D J., DeFries, R S., Heimann, M., Ojima, D S., Quegan, S., and Schmullius, C C.: Model-data synthesis in terrestrial carbon observation: methods, data requirements and data uncertainty specifications, Glob. Change Biol., 11, 378–397, 2005. Reichstein, M., Tenhunen, J., Roupsard, O., Ourcival, J M., Rambal, S., Miglietta, F., Peressotti, A., Pecchiari, M., Tirone, G., and Valentini, R.: Inverse modeling of seasonal drought effects on canopy CO₂/H₂O exchange in three Mediterranean ecosystems, J. Geophys. Res.-Atmos., 108, D23, 4726, doi:10.1029/2003JD003430, 2003. Reichstein, M., Falge, E., Baldocchi, D., Papale, D., Aubinet, M., Berbigier, P., Bernhofer, C., Buchmann, N., Gilmanov, T., Granier, A., Grunwald, T., Havrankova, K., Ilvesniemi, H., Janous, D., Knohl, A., Laurila, T., Lohila, A., Loustau, D., Matteucci, G., Meyers, T., Miglietta, F., Ourcival, J M., Pumpanen, J., Rambal, S., Rotenberg, E., Sanz, M., Tenhunen, J., Seufert, G., Vaccari, F., Vesala, T., Yakir, D., and Valentini, R.: On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm, Glob. Change Biol., 11, 1424–1439, 2005. Richardson, A D. and Hollinger, D Y.: 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, Agr. Forest Meteorol., 131, 191–208, 2005. Richardson, A D., Hollinger, D Y., Burba, G G., Davis, K J., Flanagan, L B., Katul, G G., Munger, J W., Ricciuto, D M., Stoy, P C., Suyker, A E., Verma, S B., and Wofsy, S C.: A multi-site analysis of random error in tower-based measurements of carbon and energy fluxes, Agr. Forest Meteorol., 136, 1–18, 2006. Richardson, A D., Mahecha, M D., Falge, E., Kattge, J., Moffat, A M., Papale, D., Reichstein, M., Stauch, V J., Braswell, B H., Churkina, G., Kruijt, B., and Hollinger, D Y.: Statistical properties of random CO₂ flux measurement uncertainty inferred from model residuals, Agr. Forest Meteorol., 148, 38–50, doi:10.1016/j.agrformet.2007.09001, 2007. Sacks, W J., Schimel, D S., Monson, R K., and Braswell, B H.: Model-data synthesis of diurnal and seasonal CO₂ fluxes at Niwot Ridge, Colorado, Glob. Change Biol., 12, 240–259, 2006. Tarantola, A.: Inverse Problem Theory: Methods for Data Fitting and Model Parameter Estimation, Elsevier, New York, 1st ed., 613 pp., 1987. Trudinger, C M., Raupach, M R., Rayner, P J., Kattge, J., Li, Q., Pak, B., Reichstein, M., Renzullo, L., Richardson, A D., Roxburgh, S H., Styles, J., Wang, Y P., Briggs, P., Barrett, D., and Nikolova, S.: OptIC project: An intercomparison of optimization techniques for parameter estimation in terrestrial biogeochemical models, J. Geophys. Res., 112, G2, G02027, doi:10.1029/2006/JG000367, 2007. Visual~Numerics, Inc.: PV-Wave 8.5 Reference Guide, Houston, http://www.vni.com/books/dod/pdf/wave85Docs/eReferenceGuide85.pdf, last access: 9 September 2008, 2005. Wang, Y P., Baldocchi, D., Leuning, R., Falge, E., and Vesala, T.: Estimating parameters in a land-surface model by applying nonlinear inversion to eddy covariance flux measurements from eight FLUXNET sites, Glob. Change Biol., 13, 652–670, 2007. Wilks, D S.: Statistical Methods in the Atmospheric Sciences, Academic Press, 467 pp., 1995. Wilson, K., Goldstein, A., Falge, E., Aubinet, M., Baldocchi, D., Berbigier, P., Bernhofer, C., Ceulemans, R., Dolman, H., Field, C., Grelle, A., Ibrom, A., Law, B., Kowalski, A., Meyers, T., Moncrieff, J., Monson, R., Oechel, W., Tenhunen, J., Valentini, R., and Verma, S.: Energy balance closure at FLUXNET sites, Agr. Forest Meteorol., 113, 21, 223–243, 2002.