Atmospheric Chemistry and Physics Discussions
www.atmos-chem-phys-discuss.net
1680-7367
1680-7375
8
3
2008
10.5194/acpd-8-12023-2008
http://www.atmos-chem-phys-discuss.net/8/12023/2008/
http://www.atmos-chem-phys-discuss.net/8/12023/2008/acpd-8-12023-2008.html
http://www.atmos-chem-phys-discuss.net/8/12023/2008/acpd-8-12023-2008.pdf
12023
12052
2008-06-18
Technical Note: Four-dimensional variational data assimilation for inverse modelling of atmospheric methane emissions: method and comparison with synthesis inversion
J. F. Meirink
meirink@knmi.nl
P. Bergamaschi
M. C. Krol
Institute for Marine and Atmospheric research Utrecht (IMAU), University of Utrecht, Utrecht, The Netherlands
Royal Netherlands Meteorological Institute (KNMI), De Bilt, The Netherlands
European Commission Joint Research Centre, Institute for Environment and Sustainability (EC JRC IES), Ispra (VA), Italy
Wageningen University and Research Centre (WUR), Wageningen, The Netherlands
Netherlands Institute for Space Research (SRON), Utrecht, The Netherlands
A four-dimensional variational (4D-var) data assimilation system for inverse modelling
of atmospheric methane emissions is presented. The system is based on the TM5
atmospheric transport model. It can be used for assimilating large volumes of measurements, in particular satellite observations
and quasi-continuous in-situ observations, and at the same time it enables the optimization of a large number of model parameters,
specifically grid-scale emission rates. Furthermore, the variational method allows to estimate uncertainties in
posterior emissions. Here, the system is applied to optimize monthly methane emissions
over a 1-year time window on the basis of surface observations from the
NOAA-ESRL network. The results are rigorously compared with an analogous inversion by Bergamaschi et al. (2007), which was based
on the traditional synthesis approach. The posterior emissions
as well as their uncertainties obtained in both inversions show
a high degree of consistency. At the same time we illustrate the advantage of 4D-Var in reducing aggregation errors by
optimizing emissions at the grid scale of the transport model.
The full potential of the assimilation system is exploited in
Meirink et al. (2008), who use satellite observations of column-averaged methane
mixing ratios to optimize emissions at high spatial resolution,
taking advantage of the zooming capability of the TM5 model.
Baker, D F., Doney, S C., and Schimel, D S.: Variational data assimilation for atmospheric CO<sub>2</sub>, Tellus, 58B, 359–365, doi:10.1111/j.1600-0889.2006.00218.x, 2006.
Bergamaschi, P., Krol, M., Dentener, F., Vermeulen, A., Meinhardt, F., Graul, R., Ramonet, M., Peters, W., and Dlugokencky, E J.: Inverse modelling of national and European CH<sub>4</sub> emissions using the atmospheric zoom model TM5, Atmos. Chem. Phys., 5, 2431–2460, 2005.
Bergamaschi, P., Meirink, J F., Müller, J F., Körner, S., Heimann, M., Dlugokencky, E J., Kaminski, U., Marcazzan, G., Vecchi, R., Meinhardt, F., Ramonet, M., Sartorius, H., and Zahorowski, W.: Model inter-comparison on transport and chemistry – report on model inter-comparison performed within European Commission FP5 project EVERGREEN (Global satellite observation of greenhouse gas emissions), Sci. and Tech. Res. Ser. EUR 22241 EN, 53 pp., Eur. Comm., DG Joint Res. Cent., Inst. for Environ.\ and Sustainability, Ispra, Italy, 2006.
Bergamaschi, P., Frankenberg, C., Meirink, J F., Krol, M., Dentener, F., Wagner, T., Platt, U., Kaplan, J O., Körner, S., Heimann, M., Dlugokencky, E J., and Goede, A.: Satellite chartography of atmospheric methane from SCIAMACHY on board ENVISAT: 2. Evaluation based on inverse model simulations, J. Geophys. Res., 112, D02304, doi:10.1029/2006JD007268, 2007.
Bousquet, P., Ciais, P., Miller, J. B., Dlugokencky, E. J., Hauglustaine, D. A., Prigent, C., van der Werf, G. R., Peylin, P., Brunke, E. G., Carouge, C., Langenfelds, R. L., Lathière, J., Papa, F., Ramonet, M., Schmidt, M., Steele, L. P., Tyler, S. C., and White, J.: Contribution of anthropogenic and natural sources to atmospheric methane variability, Nature, 443, 439–443, doi:10.1038/nature05132, 2006.
Buchwitz, M., de~Beek, R., Burrows, J. P., Bovensmann, H., Warneke, T., Notholt, J., Meirink, J. F., Goede, A. P. H., Bergamaschi, P., Körner, S., Heimann, M., and Schulz, A.: Atmospheric methane and carbon dioxide from SCIAMACHY satellite data: initial comparison with chemistry and transport models, Atmos. Chem. Phys., 5, 941–962, 2005.
Buchwitz, M., de~Beek, R., No\"el, S., Burrows, J. P., Bovensmann, H., Schneising, O., Khlystova, I., Bruns, M., Bremer, H., Bergamaschi, P., Körner, S., and Heimann, M.: Atmospheric carbon gases retrieved from SCIAMACHY by WFM-DOAS: version 0.5 CO and CH<sub>4</sub> and impact of calibration improvements on CO<sub>2</sub> retrieval, Atmos. Chem. Phys., 6, 2727–2751, 2006.
Chen, Y.-H. and Prinn, R. G.: Estimation of atmospheric methane emissions between 1996 and 2001 using a three-dimensional global chemical transport model, J. Geophys. Res., 111, D10307, doi:10.1029/2005JD006058, 2006.
Chevallier, F., Fisher, M., Peylin, P., Serrar, S., Bousquet, P., Bréon, F.-M., Chédin, A., and Ciais, P.: Inferring CO<sub>2</sub> sources and sinks from satellite observations: Method and application to TOVS data, J. Geophys. Res., 110, D24309, doi:10.1029/2005JD006390, 2005.
Chevallier, F., Bréon, F.-M., and Rayner, P. J.: Contribution of the Orbiting Carbon Observatory to the estimation of CO<sub>2</sub> sources and sinks: Theoretical study in a variational data assimilation framework, J. Geophys. Res., 112, D09307, doi:10.1029/2006JD007375, 2007.
Courtier, P., Thépaut, J.-N., and Hollingsworth, A.: A strategy for operational implementation of 4D-Var, using an incremental approach, Q. J. Roy. Meteor. Soc., 120, 1367–1387, 1994.
Dlugokencky, E. J., Steele, L. P., Lang, P. M., and Masarie, K. A.: The growth rate and distribution of atmospheric methane, J. Geophys. Res., 99, 17 021–17 044, doi:10.1029/94JD01245, 1994.
Dlugokencky, E J., Houweling, S., Bruhwiler, L., Masarie, K A., Lang, P M., Miller, J B., and Tans, P P.: Atmospheric methane levels off: Temporary pause or a new steady-state?, Geophys. Res. Let., 30, 1992, doi:10.1029/2003GL018126, 2003.
Elbern, H., Strunk, A., Schmidt, H., and Talagrand, O.: Emission rate and chemical state estimation by 4-dimensional variational inversion, Atmos. Chem. Phys., 7, 3749–3769, 2007.
Engelen, R J., Denning, A S., Gurney, K R., and TransCom3 modelers: On error estimation in atmospheric CO<sub>2</sub> inversions, J. Geophys. Res., 107, 4635, doi:10.1029/2002JD002195, 2002.
Enting, I G.: Inverse problems in atmospheric constituent transport, Cambridge University Press, New York, 2002.
Evensen, G.: Sequential Data Assimilation with a Nonlinear Quasi-Geostrophic Model Using Monte-Carlo Methods to Forecast Error Statistics, J. Geophys. Res., 99, 10 143–10 162, 1994.
Fisher, M. and Courtier, P.: Estimating the covariance matrices of analysis and forecast error in variational data assimilation, Technical Memorandum 220, ECMWF, Reading, UK, 1995.
Frankenberg, C., Meirink, J F., van Weele, M., Platt, U., and Wagner, T.: Assessing methane emissions from global space-borne observations, Science, 308, 1010–1014, 2005.
Frankenberg, C., Meirink, J F., Bergamaschi, P., Goede, A. P H., Heimann, M., Körner, S., Platt, U., van Weele, M., and Wagner, T.: Satellite chartography of atmospheric methane from SCIAMACHY on board ENVISAT: Analysis of the years 2003 and 2004, J. Geophys. Res., 111, D07303, doi:10.1029/2005JD006235, 2006.
Giering, R. and Kaminski, T.: Recipes for adjoint code construction, ACM Trans. On Math. Software, 24, 437–474, 1998.
Hein, R., Crutzen, P J., and Heimann, M.: An inverse modelling approach to investigate the global atmospheric methane cycle, Global Biogeochem. Cy., 11, 43–76, 1997.
Houweling, S., Kaminski, T., Dentener, F J., Lelieveld, J., and Heimann, M.: Inverse modelling of methane sources and sinks using the adjoint of a global transport model, J. Geophys. Res., 104, 26 137–26 160, 1999.
Kaminski, T., Heimann, M., and Giering, R.: A coarse grid three-dimensional global inverse model of the atmospheric transport, 2, Inversion of the transport of CO<sub>2</sub> in the 1980s, J. Geophys. Res., 104, 18 555–18 581, 1999.
Kaminski, T., Rayner, P J., Heimann, M., and Enting, I G.: On aggregation errors in atmospheric transport inversions, J. Geophys. Res., 106, 4703–4716, doi:10.1029/2000JD900581, 2001.
Krol, M., Houweling, S., Bregman, B., van den Broek, M., Segers, A., van Velthoven, P., Peters, W., Dentener, F., and Bergamaschi, P.: The two-way nested global chemistry – transport zoom model TM5: algorithm and applications, Atmos. Chem. Phys., 5, 417–432, 2005.
Krol, M. C., Meirink, J. F., Bergamaschi, P., Mak, J. E., Lowe, D., Jöckel, P., Houweling, S., and Röckmann, T.: What do $^14$CO measurements tell us about OH?, Atmos. Chem. Phys. Discuss., 7, 10 405–10 438, 2007.
Lanczos, C.: An iteration method for the solution of the eigenvalue problem of linear differential and integral operators, J. Res. Nat. Bur. Standards, 45, 255–282, 1950.
Meirink, J. F., Eskes, H. J., and Goede, A. P. H.: Sensitivity analysis of methane emissions derived from SCIAMACHY observations through inverse modelling, Atmos. Chem. Phys., 6, 1275–1292, 2006.
Meirink, J F., Bergamaschi, P., Frankenberg, C., d'Amelio, M. T S., Dlugokencky, E J., Gatti, L V., Houweling, S., Miller, J B., Röckmann, T., Villani, M G., and Krol, M.: Four-dimensional variational data assimilation for inverse modelling of atmospheric methane emissions: Analysis of SCIAMACHY observations, J. Geophys. Res., accepted, 2008.
Mikaloff Fletcher, S. E M., Tans, P. P., Bruhwiler, L. M., Miller, J. B., and Heimann, M.: CH<sub>4</sub> sources estimated from atmospheric observations of CH<sub>4</sub> and its $^13$C/$^12$C isotopic ratios: 1 Inverse modeling of source processes, Global Biogeochem. Cy., 18, GB4004, doi:10.1029/2004GB002223, 2004.
Parrish, D. F. and Derber, J. C.: The National Meteorological Center's spectral statistical-interpolation analysis system, Mon. Weather Rev., 120, 1747–1763, 1992.
Peters, W., Miller, J., Whitaker, J., Denning, S., Hirsch, A., Krol, M., Zupanski, D., Bruhwiler, L., and Tans, P.: An ensemble data assimilation system to estimate CO<sub>2</sub> surface fluxes from atmospheric trace gas observations, J. Geophys. Res., 110, D24304, doi:10.1029/2005JD006157., 2005.
Russell, G. A. and Lerner, J. A.: A new finite-differencing scheme for the tracer transport equation, J. Appl. Meteorol., 20, 1482–1498, 1981.
Stavrakou, T. and Müller, J.-F.: Grid-based versus big region approach for inverting CO emissions using Measurement of Pollution in the Troposphere (MOPITT) data, J. Geophys. Res., 111, D15304, doi:10.1029/2005JD006896, 2006.
Talagrand, O. and Courtier, P.: Variational assimilation of meteorological observations with the adjoint vorticity equation. I: Theory, Q. J. Roy. Meteor. Soc., 113, 1311–1328, 1987.
Tarantola, A.: Inverse problem theory and methods for model parameter estimation, Soc. for Ind. and Appl. Math., Philadelphia PA, USA, 2005.
Yumimoto, K., Uno, I., Sugimoto, N., Shimizu, A., and Satake, S.: Adjoint inverse modeling of dust emission and transport over East Asia, Geophys. Res. Let., 34, L08806, doi:10.1029/2006GL028551, 2007.