References

ACPD

Atmospheric Chemistry and Physics Discussions

ACPD

1680-7375

Copernicus GmbH

Göttingen, Germany

10.5194/acpd-11-3355-2011

A comparison of different inverse carbon flux estimation approaches for application on a regional domain

Tolk

L. F.

¹ Dolman

A. J.

¹ Meesters

A. G. C. A.

¹ Peters

Vrije Universiteit Amsterdam, Faculty of Earth and Life Science, Amsterdam, The Netherlands

Wageningen University, Dept. of Meteorology and Air Quality, Wageningen, The Netherlands

31 01 2011

11 1 3355 3398

This is an open-access article ditributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

This article is available from http://www.atmos-chem-phys-discuss.net/11/3355/2011/acpd-11-3355-2011.html

The full text article is available as a PDF file from http://www.atmos-chem-phys-discuss.net/11/3355/2011/acpd-11-3355-2011.pdf

We have implemented six different inverse carbon flux estimation methods in a regional carbon dioxide (CO2) flux modeling system for The Netherlands. The system consists of the Regional Atmospheric Mesoscale Modeling System (RAMS) coupled to a simple carbon flux scheme which is run in a coupled fashion on relatively high resolution (10 km). Using an Ensemble Kalman filter approach we try to estimate spatiotemporal carbon exchange patterns from atmospheric CO2 mole fractions over The Netherlands for a two week period in spring 2008. The focus of this work is the different strategies that can be employed to turn first-guess fluxes into optimal ones, which is known as a fundamental design choice that can affect the outcome of an inversion significantly. Different state-of-the-art approaches with respect to the estimation of net ecosystem exchange (NEE) are compared quantitatively: (1) where NEE is scaled by one linear multiplication factor per land-use type, (2) where the same is done for photosynthesis (GPP) and respiration (R) separately with varying assumptions for the correlation structure, (3) where we solve for those same multiplication factors but now for each grid box, and (4) where we optimize physical parameters of the underlying biosphere model for each land-use type. The pattern to be retrieved in this pseudo-data experiment is different in nearly all aspects from the first-guess fluxes, including the structure of the underlying flux model, reflecting the difference between the modeled fluxes and the fluxes in the real world. This makes our study a stringent test of the performance of these methods, which are currently widely used in carbon cycle inverse studies. Our results show that all methods struggle to retrieve the spatiotemporal NEE distribution, and none of them succeeds in finding accurate domain averaged NEE with correct spatial and temporal behavior. The main cause is the difference between the structures of the first-guess and true CO2 flux models used. Most methods display overconfidence in their estimate as a result. A commonly used daytime-only sampling scheme in the transport model leads to compensating biases in separate GPP and R scaling factors that are readily visible in the nighttime mixing ratio predictions of these systems. Overall, we recommend that the estimate of NEE scaling factors should not be used in this regional setup, while estimating bias factors for GPP and R for every grid box works relatively well. The biosphere parameter inversion is best at simultaneously producing space and time patterns of fluxes and CO2 mixing ratios, but non-linearity may significantly reduce the information content in the inversion if true parameter values are far from the prior estimate. Our results suggest that a carefully designed biosphere model parameter inversion or a pixel inversion of the respiration and GPP multiplication factors are the best tools to optimize spatiotemporal patterns of NEE.

References 1

Ahmadov,~R., Gerbig,~C., Kretschmer,~R., Körner,~S., Rödenbeck,~C., Bousquet,~P., and Ramonet,~M.: Comparing high resolution WRF-VPRM simulations and two global \chemCO_2 transport models with coastal tower measurements of \chemCO_2, Biogeosciences, 6, 807–817, doi:10.5194/bg-6-807-2009, 2009.

Arneth,~A., Lloyd,~J., Santruckova,~H., Bird,~M., Grigoryev,~S., Kalaschnikov,~Y N., Gleixner,~G., and Schulze,~E.-D.: Response of Central Siberian Scots Pine to soil water deficit and long-term trends in atmospheric \chemCO_2 concentration, Global Biogeochem. Cy., 16, 1005, \doi10.1029/2000GB001374,2002, 2002.

Baldocchi,~D., Falge,~E., Gu,~L., Olson,~R., Hollinger,~D., Running,~S., Anthoni,~P., Bernhofer,~C., Davis,~K., Evans,~R., Fuentes,~J., Goldstein,~A., Katul,~G., Law,~B., Lee,~X., Mahli,~Y., Meyers,~T., Munger,~W., Oechel,~W., Paw,~K T., Pilegaard,~K., Schmid,~H P., Valentini,~R., Verma,~S., Vesala,~T., Wilson,~K., and Wofsy,~S.: 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, \doi10.1175/1520-0477, 2001.

Bousquet,~P., Peylin,~P., Ciais,~P., Le Quéré,~C., Friedlingstein,~P., and Tans,~P P.: Regional changes in carbon dioxide fluxes of land and oceans since 1980, Science, 290, 1342–1346, 2000.

von Caemmerer,~S. and Farquhar,~G D.: Some relationships between the bio-chemistry of photosynthesis and the gas exchange of leaves, Planta, 153, 376–387, 1981.

%Carouge,~C., Rayner,~P J., Peylin,~P., Bousquet,~P., Chevallier,~F., and Ciais,~P.: What can we learn from European continuous atmospheric \chemCO_2 measurements to quantify regional fluxes – Part II: Sensitivity of flux accuracy to inverse setup, Atmospheric Chemistry and Physics, 10, 3119–3129, 2010. % ### SELF-REFERENCE ### Carouge, C., Rayner, P. J., Peylin, P., Bousquet, P., Chevallier, F., and Ciais, P.: What can we learn from European continuous atmospheric \chemCO_2 measurements to quantify regional fluxes – Part 2: Sensitivity of flux accuracy to inverse setup, Atmos. Chem. Phys., 10, 3119–3129, doi:10.5194/acp-10-3119-2010, 2010.

Carvalhais,~N., Reichstein,~M., Seixas,~J., Collatz,~G J., Pereira,~J S., Berbigier,~P., Carrara,~A., Granier,~A., Montagnani,~L., Papale,~D., and Rambal, S., Jos$\acute\rm e$ Sanz M., Valentini, R.: Implications of the carbon cycle steady state assumption for biogeochemical modeling performance and inverse parameter retrieval, Global Biogeochem. Cy., 22, GB2007, doi:10.1029/2007GB003033, 2008.

Ciais,~P., Rayner,~P., Chevallier,~F., Bousquet,~P., Logan,~M., Peylin,~P., and Ramonet,~M.: Atmospheric inversions for estimating \chemCO_2 fluxes: methods and perspectives, Climatic Change, 103, 69–92, \doi10.1007/s10584-010-9909-3, 2010.

Collatz,~G., Ball,~J., Grivet,~C., and Berry,~J.: Physiological and environmental regulation of stomatal conductance, photosynthesis and transpiration: a~model that includes a~laminar boundary layer, Agr. Forest Meteorol., 54, 107–136, 1991.

Farquhar,~G D., Caemmerer,~S V., and Berry,~J A.: A~biochemical model of photosynthetic \chemCO_2 assimilation in leaves of C 3 species, Planta, 149, 78–90, 1980.

%Gerbig,~C., Lin,~J C., Munger,~J W., and Wofsy,~S C.: What can tracer observations in the continental boundary layer tell us about surface-atmosphere fluxes?, Atmospheric Chemistry and Physics, 6, 539–554, 2006. % ### SELF-REFERENCE ### Gerbig, C., Lin, J. C., Munger, J. W., and Wofsy, S. C.: What can tracer observations in the continental boundary layer tell us about surface-atmosphere fluxes?, Atmos. Chem. Phys., 6, 539–554, doi:10.5194/acp-6-539-2006, 2006.

%Gerbig,~C., Körner,~S., and Lin,~J C.: Vertical mixing in atmospheric tracer transport models: error characterization and propagation, Atmospheric Chemistry and Physics, 8, 591–602, 2008. % ### SELF-REFERENCE ### Gerbig, C., Körner, S., and Lin, J. C.: Vertical mixing in atmospheric tracer transport models: error characterization and propagation, Atmos. Chem. Phys., 8, 591–602, doi:10.5194/acp-8-591-2008, 2008.

%Gourdji,~S M., Hirsch,~A I., Mueller,~K L., Yadav,~V., Andrews,~A E., and Michalak,~A M.: Regional-scale geostatistical inverse modeling of North American \chemCO_2 fluxes: a~synthetic data study, Atmos. Chem. Phys, 10, 6151–6167, 2010. % ### SELF-REFERENCE ### Gourdji, S. M., Hirsch, A. I., Mueller, K. L., Yadav, V., Andrews, A. E., and Michalak, A. M.: Regional-scale geostatistical inverse modeling of North American \chemCO_2 fluxes: a synthetic data study, Atmos. Chem. Phys., 10, 6151–6167, doi:10.5194/acp-10-6151-2010, 2010.

Groenendijk,~M., Dolman,~A J., van der Molen,~M K., Leuning,~R., Arneth,~A., Delpierre,~N., Gash,~J H C., Lindroth,~A., Richardson,~A D., Verbeeck,~H., and Wohlfahrt,~G.: Assessing parameter variability in a~photosynthesis model within and between plant functional types using global fluxnet eddy covariance data, Agr. Forest Meteorol. 151(15), 22–38, ISSN 0168-1923, doi:10.1016/j.agrformet.2010.08.013.

Gurney,~K R., Law,~R M., Denning,~A S., Rayner,~P J., Baker,~D., Bousquet,~P., Bruhwiler,~L., Chen,~Y H., Ciais,~P., Fan,~S M., Fung, I. Y., Gloor, M., Heimann, M., Higuchi, K. A. Z., John, J., Kowalczyk, E. V. A., Maki, T., Maksyutov, S., Peylin, P., Prather, M., Pak, B. C., Sarmiento, J., Taguchi, S., Takahashi, T., and Yuen, C.-W.: TransCom 3 \chemCO_2 inversion intercomparison: 1. annual mean control results and sensitivity to transport and prior flux information, Tellus B, 55, 555–579, 2003.

Kaminski,~T., Rayner,~P J., Heimann,~M., and Enting,~I.: On aggregation errors in atmospheric transport inversions, J. Geophys. Res., 106(D5), 4703–4715, 2001.

Law,~R M., Peters,~W., Rödenbeck,~C., Aulagnier,~C., Baker,~I., Bergmann, ~D J., Bousquet,~P., Brandt,~J., Bruhwiler,~L., Cameron-Smith,~P J., Christensen, J. H., Delage, F., Denning, A. S., Fan, S., Geels, C., Houweling, S., Imasu, R., Karstens, U., Kawa, S. R., Kleist, J., Krol, M. C., Lin, S. J., Lokupitiya, R., Maki, T., Maksyutov, S., Niwa, Y., Onishi, R., Parazoo, N., Patra, P. K., Pieterse, G., Rivier, L., Satoh, M., Serrar, S., Taguchi, S., Takigawa, M., Vautard, R., Vermeulen, A. T., and Zhu, Z.: TransCom model simulations of hourly atmospheric \chemCO_2: experimental overview and diurnal cycle results for 2002, Global Biogeochem. Cy., 22, GB3009, doi:10.1029/2007GB003050, 2008.

Lloyd,~J. and Taylor,~J A.: On the temperature dependence of soil respiration, Funct. Ecol., 8, 315–323, 1994.

Lokupitiya,~R S., Zupanski,~D., Denning,~A S., Kawa,~S R., Gurney,~K R., and Zupanski,~M.: Estimation of global \chemCO_2 fluxes at regional scale using the maximum likelihood ensemble filter, J. Geophys. Res., 113, D20110, doi:10.1029/2003JD004422, 2008.

Meesters,~A G C A., Tolk,~L F., and Dolman,~A J.: Mass conservation above slopes in the Regional Atmospheric Modelling System (RAMS), Environ. Fluid Mech., 8, 239–248, 2008.

Michalak,~A M., Bruhwiler,~L., and Tans,~P P.: A~geostatistical approach to surface flux estimation of atmospheric trace gases, J. Geophys. Res., 109, D14109, doi:10.1029/2007JD009679, 2004.

Mueller,~K L., Gourdji,~S M., Michalak,~A M.: Global monthly averaged \chemCO_2 fluxes recovered using a~geostatistical inverse modeling approach: 1. results using atmospheric measurements, J. Geophys. Res.-Atmos., 113(D21), D21114, doi:10.1029/2003JD004422, 2008.

Papale,~D. and Valentini,~R.: A~new assessment of European forests carbon exchanges by eddy fluxes and artificial neural network spatialization, Glob. Change Biol., 9, 525–535, 2003.

Patil,~D J., Hunt,~B R., Kalnay,~E., Yorke,~J A., and Ott,~E.: Local low dimensionality of atmospheric dynamics, Phys. Rev. Lett., 86, 5878–5881, 2001.

Peters,~W., Miller,~J B., Whitaker,~J., Denning,~A S., Hirsch,~A., Krol,~M C., Zupanski,~D., Bruhwiler,~L., and Tans,~P P.: An ensemble data assimilation system to estimate \chemCO_2 surface fluxes from atmospheric trace gas observations,~J. Geophys. Res., 110, D24304, doi:10.1029/2005JD006157, 2005.

Peters,~W., Jacobson,~A R., Sweeney,~C., Andrews,~A E., Conway,~T J., Masarie,~K., Miller,~J B., Bruhwiler,~L M P., Petron,~G., and Hirsch,~A I.: An atmospheric perspective on North American carbon dioxide exchange: CarbonTracker, P. Natl. Acad. Sci. USA, 104, 18925–18930, 2007.

Peters,~W., Krol,~M C., van der Werf,~G R., Houweling,~S., Jones,~C D., Hughes,~J., Scheafer,~K., Masarie,~K A., Jacobson,~A R., Miller,~J B., Cho,~C H., Ramonet, M., Schmidt, M., Ciattaglia, L., Apadula, F., Heltai, D., Meinhardt, F., Di Sarra, A. G., Piacentino, S., Sferlazzo, D., Aalto, T., Hatakka, J., Str�M, J., Haszpra, L., Meijer, H. A. J., Van Der Laan, S., Neubert, R. E. M., Jordan, A., Rod�, X., Morgu�, J. A., Vermeulen, A. T., Popa, E., Rozanski, K., Zimnoch, M., Manning, A. C., Leuenberger, M., Uglietti, C., Dolman, A. J., Ciais, P., Heimann, M., and Tans, P. P.: Seven years of recent European net terrestrial carbon dioxide exchange constrained by atmospheric observations, Glob. Change Biol., 16, 1317–1337, 2010.

%Peylin,~P., Rayner,~P., Bousquet,~P., Carouge,~C., Hourdin,~F., Heinrich,~P., and Ciais,~P.: Daily \chemCO_2 flux estimates over Europe from continuous atmospheric measurements: 1. inverse methodology, Atmos. Chem. Phys, 5, 3173–3186, 2005. % ### SELF-REFERENCE ### Peylin, P., Rayner, P. J., Bousquet, P., Carouge, C., Hourdin, F., Heinrich, P., Ciais, P., and AEROCARB contributors: Daily \chemCO_2 flux estimates over Europe from continuous atmospheric measurements: 1, inverse methodology, Atmos. Chem. Phys., 5, 3173–3186, doi:10.5194/acp-5-3173-2005, 2005.

Pielke,~R A., Cotton,~W R., Walko,~R L., Tremback,~C J., Lyons,~W A., Grasso,~L D., Nicholls,~M E., Moran,~M D., Wesley,~D A., Lee,~T J., and Copeland,~J H.: A~comprehensive meteorological modeling system RAMS, Meteorol. Atmos. Phys., 49, 69–91, 1992.

Pieterse,~G., Bleeker,~A., Vermeulen,~A T., Wu,~Y., and Erisman,~J W.: High resolution modelling of atmosphere canopy exchange of acidifying and eutrophying components and carbon dioxide for European forests, Tellus B, 59, 412–424, 2007.

Rayner,~P J., Scholze,~M., Knorr,~W., Kaminski,~T., Giering,~R., and Widmann,~H.: Two decades of terrestrial carbon fluxes from a~carbon cycle data assimilation system (CCDAS), Global Biogeochem. Cy., 19, GB2026, doi:10.1029/2004GB002254, 2005.

%Rödenbeck,~C., Houweling,~S., Gloor,~M., and Heimann,~M.: \chemCO_2 flux history 1982–2001 inferred from atmospheric data using a~global inversion of atmospheric transport, Atmospheric Chemistry and Physics, 3, 1919–1964, 2003. % ### SELF-REFERENCE ### Rödenbeck, C., Houweling, S., Gloor, M., and Heimann, M.: \chemCO_2 flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport, Atmos. Chem. Phys., 3, 1919–1964, doi:10.5194/acp-3-1919-2003, 2003.

Rüdiger,~C., Albergel,~C., Mahfouf,~J F., Calvet,~J C., and Walker,~J P.: Evaluation of the observation operator Jacobian for leaf area index data assimilation with an extended Kalman filter, J. Geophys. Res., 115, D09111, doi:10.1029/2009JD012912, 2010.

Scholze,~M., Kaminski,~T., Rayner,~P., Knorr,~W., and Giering,~R.: Propagating uncertainty through prognostic carbon cycle data assimilation system simulations,~J. Geophys. Res., 112, D17305, 2007.

Schuh,~A E., Denning,~A S., Uliasz,~M., and Corbin,~K D.: Seeing the forest through the trees: recovering large-scale carbon flux biases in the midst of small-scale variability, J. Geophys. Res., 114, G03007, doi:10.1029/2007JD008642, 2009.

Sellers,~P., Randall,~D., Collatz,~G., Berry,~J., Field,~C., Dazlich,~D., Zhang,~C., Collelo,~G., and Bounoua,~L.: A~revised land surface parameterization (SiB2) for atmospheric GCMs. Part I: model formulation, J. Climate, 9, 676–705, 1996.

Steeneveld,~G J., Mauritsen,~T., de Bruijn,~E I F., de Arellano,~J V G., Svensson,~G., and Holtslag,~A A M.: Evaluation of limited-area models for the representation of the diurnal cycle and contrasting nights in CASES-99, J. Appl. Meteorol. Clim., 47, 869–887, 2008.

%Tolk,~L F., Peters,~W., Meesters,~A G C.A.., Groenendijk,~M., Vermeulen,~A T., Steeneveld,~G J., and Dolman,~A J.: Modelling regional scale surface fluxes, meteorology and \chemCO_2 mixing ratios for the Cabauw tower in The Netherlands, Biogeosciences, 6, 2265–2280, 2009. % ### SELF-REFERENCE ### Tolk,~L F., Peters,~W., Meesters,~A G C A., Groenendijk,~M., Vermeulen,~A T., Steeneveld,~G J., and Dolman,~A J.: Modelling regional scale surface fluxes, meteorology and \chemCO_2 mixing ratios for the Cabauw tower in The Netherlands, Biogeosciences, 6, 2265–2280, doi:10.5194/bg-6-2265-2009, 2009.

Trudinger,~C M., Raupach,~M R., Rayner,~P J., Kattge,~J., Liu,~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, G02027, doi:10.1029/2008JG000842, 2007.

%Yadav,~V., Mueller,~K L., Dragoni,~D., and Michalak,~A M.: A~geostatistical synthesis study of factors affecting gross primary productivity in various ecosystems of North America, Biogeosciences, 7, 2655–2671, 2010. % ### SELF-REFERENCE ### Yadav,~V., Mueller,~K L., Dragoni,~D., and Michalak,~A M.: A geostatistical synthesis study of factors affecting gross primary productivity in various ecosystems of North America, Biogeosciences, 7, 2655–2671, doi:10.5194/bg-7-2655-2010, 2010.

Zupanski,~D., Denning,~A S., Uliasz,~M., Zupanski,~M., Schuh,~A E., Rayner,~P J., Peters,~W., and Corbin,~K D.: Carbon flux bias estimation employing maximum likelihood ensemble filter (MLEF), J. Geophys. Res., 112, D17107, doi:10.1029/2006JD008371, 2007.