Difference between revisions of "GEOS-Chem Adjoint"

From Geos-chem
Jump to: navigation, search
Line 45: Line 45:
** Initial Conditions
** Initial Conditions
** Reaction Rate Parameters
** Reaction Rate Parameters
* [http://puck.che.caltech.edu/~daven/home/GC_adj/GC_adj_man.pdf, Code manual]
* [http://puck.che.caltech.edu/~daven/home/GC_adj/GC_adj_man.pdf Code manual]
* [http://puck.che.caltech.edu/~daven/home/GC_adj/Diagrams.pdf, Code diagrams]
* [http://puck.che.caltech.edu/~daven/home/GC_adj/Diagrams.pdf Code diagrams]
==== Primary Developer ====
==== Primary Developer ====

Revision as of 19:48, 16 October 2008


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.


  • Caltech: John Seinfeld
  • Carleton University: Amir Hakami
  • Columbia / NASA GISS / CU Boulder : Daven Henze
  • Dalhousie: Randall Martin, Chulkyu Lee
  • Harvard: Monika Kopacz, Lin Zhang, Daniel Jacob
  • JPL: Kevin Bowman, Qinbin Li, Xun Jiang, Changsub Shim, Yang Chen, Sunita Verma
  • U. of Toronto: Dylan Jones, Mark Parrington, Michael Seymour, Thomas Walker
  • 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


  • GEOS-3 and GEOS-4 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
      • now includes volcano SO2 (C. Lee)
      • now includes sensitivity w.r.t all emissions, though only some implemented for inverse modeling
    • Initial Conditions
    • Reaction Rate Parameters
  • Code manual
  • Code diagrams

Primary Developer

Daven Henze


Daven Henze (Columbia University), Chulkyu Lee (Dalhousie), Dylan Jones, Mark Parrington, Thomas Walker (U. Toronto), Changsub Shim (JPL)


  • Henze, D. K., A. Hakami and J. H. Seinfeld (2007), Development of the adjoint of GEOS-Chem, Atmos. Chem. Phys., 7, 2413-2433.
  • Henze, D. K., J. H. Seinfeld and D. T. Shindell (2008), Inverse modeling and mapping U.S. air quality influences of inorganic PM2.5 precursor emissions with the adjoint of GEOS-Chem, Atmos. Chem. Phys. Discuss., 8, 15033-15099.

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


  • 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


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


  • GEOS-4 met fields
  • 4x5 model resolution
  • Observational Operators
    • TES O3

Primary Developer

Kumaresh Singh