GEOS-Chem Adjoint Model

General

Historical Overview

Original work on the adjoint of GEOS-Chem began in 2003, focusing on the adjoint of the offline aerosol simulation. By 2005, the adjoint was expanded to include a tagged CO simulation and a full chemistry simulation as well as observational operators for MOPITT (CO) and IMPROVE network (aerosols). The adjoint code has been constructed in a hybrid fashion using a combination of automatic differentiation software (TAMC, KPP) and manual coding of both discrete and continuous adjoints. Current development aims to standardize the implementation of the adjoint model and make it more user friendly, with the eventual goal of making the adjoint publicly available as part of the standard GEOS-Chem code base.

Developers

• Caltech: John Seinfeld
• Carleton University: Amir Hakami
• Columbia / NASA GISS / CU Boulder : Daven Henze
• Dalhousie: Randall Martin, Chulkyu Lee
• Harvard: Monika Kopacz, Daniel Jacob
• JPL: Kevin Bowman, Qinbin Li, Xun Jiang, Changsub Shim, Yang Chen, Sunita Verma
• U. of Toronto: Dylan Jones, Mark Parrington, Michael Seymour
• Virginia Tech: Adrian Sandu, Kumaresh Singh, Paul Eller

Current Activities

• streamlining GEOS4 adjoint transport among the three code branches, some bug fixes to be implemented

Get Involved

Support and development of the GEOS-Chem adjoint has clearly become a multifaceted undertaking. In attempt to focus such developments, we've created a GC-adjoint mailing list. If you wish to become involved, you should add yourself to this list by sending an email to geos-chem-adjoint-requests@geos.as.harvard.edu with the word subscribe in the subject or body.

Kevin Bowman organizes bimonthly conference calls concerning current adjoint model activities.

Adjoint Model Branches

There are currently three main branches of the GEOS-Chem adjoint model.

• GEOS-Chem v6-02-05, full chemistry with online aerosols
• GEOS-Chem v6, tagged CO
• GEOS-Chem v7, full chemistry

GEOS-Chem v6-02-05, full chemistry with online aerosols

Features

• GEOS-3 met fields
• 4x5 model resolution
• Observational Operators
  • IMPROVE PM2.5 (NO3, SO4, OC, BC)
  • CASTNet (NH4+)
  • GOME / SCIAMACHY NO2 column
    • using KNMI retrieval (Henze)
    • using Dalhousie retrieval (Shim)
• Control parameters
  • Emissions (SOx, NOx, NH3, BC, OC)
    • now includes lightning NOx
    • now includes soil NOx
    • working on volcano SO₂ (C. Lee)
  • Initial Conditions
  • Reaction Rate Parameters
• Code manual
• Code diagrams

Primary Developer

Daven Henze

Publications

• Henze, D. K., A. Hakami and J. H. Seinfeld (2007), Development of the adjoint of GEOS-Chem, Atmos. Chem. Phys., 7, 2413-2433.

GEOS-Chem v6-02-05 (builds on dkh version), tagged CO

Features

• GEOS-3, GEOS-4 met fields
• 4x5, 2x2.5 model resolution
• Observational Operators
  • MOPITT CO column
  • SCIAMACHY CO colum
    • using Bremen retrieval
  • AIRS CO column
• Control parameters
  • Emissions (CO)
  • Initial Conditions
• computing Jacobians feature

Primary Developer

Monika Kopacz

Users

Monika Kopacz (Harvard), Guergana Guerova (U. Wollongong), Dylan Jones (U. Toronto), Lin Zhang (Harvard, Ox branch)

Publications and presentations

• Kopacz, M., D. J. Jacob, D. K. Henze, C. L. Heald, D. G. Streets, and Q. Zhang (2008), A comparison of analytical and adjoint Bayesian inversion methods for constraining Asian sources of CO using satellite (MOPITT) measurements of CO columns, J. Geophys. Res., doi:0.1029/2007JD009264.
• Adjoint inversion of CO sources using combined MOPITT, SCIAMACHY and AIRS CO columns, presented by Monika Kopacz at the COSPAR Scientific Assembly, Montreal, July 18, 2008. http://www-as.harvard.edu/chemistry/trop/presentations/powerpoints/mak2008/COSPAR_MKopacz_July2008.ppt

GEOS-Chem v7, full chemistry

NOTE: this version is currently under development

Features

• GEOS-3, GEOS-4, GEOS-5 met fields
• 4x5, 2x2.5 model resolution
• Observational Operators
  • TES O3

Primary Developer

Kumaresh Singh