150 likes | 237 Views
Adapting GEOS_CHEM to Mars Photochemisty. Huiqun Wang 1 Kelly Chance 1 Daniel Jacob 2 Mark Richardson 3 Michael Mischna 4 1 Center for Astrophysics (CFA) 2 Harvard 3 CALTECH 4 JPL. Motivation. 1D photochemical model for Mars 3D Mars GCM
E N D
Adapting GEOS_CHEM to Mars Photochemisty Huiqun Wang1 Kelly Chance1 Daniel Jacob2 Mark Richardson3 Michael Mischna4 1 Center for Astrophysics (CFA) 2 Harvard 3 CALTECH 4 JPL
Motivation • 1D photochemical model for Mars • 3D Mars GCM • Observations of trace gases and aerosols in the Martian atmosphere
Mars atmospheric composition 1 prµm = 3.34e+18 cm-2; 1 µm-atm = 2.69e+15 cm-2 = 0.1DU
CO2 stability problem Observation: CO: 7.0e-4 O2: 1.3e-3 O3: < a few DU Prediction: CO: 7.72e-2 O2: 3.87e-2 O3: 126DU
HOx catalytic chemistry [McElroy & Donahue, 1972; Parkinson & Hunten, 1972]
HOx-NOx catalytic chemistry Heterogeneous chemistry
Mars GCMs CO2 IR heating Dust visible & IR heating Surface physics Subgrid-scale processes CO2 cycle Water cycle Dust cycle Atmospheric chemistry • NASA AMES C-grid, p coord, 0-80km • GFDL A-grid, σ-p coord, 20–40 layers, 0-80km • LMD/AOPP A-grid/spectral, σ coord, 32 layers, 0-120/240km • Planetary WRF C-grid, σ coord, 25 or 40 layers, 0-80km …
MEX simultaneous observations of water vapor and ozone [Perrier et al., 2006] SPICAM Water Vapor Column abundance [Fedorova et al., 2006]
3D Mars photochemical model • Lefevre et al., 2004. Three-dimensional modeling of ozone on Mars. J. Geophys. Res.,109, doi:10.1029/2004JE002268. • Hartogh et al., 2006. A new coupled 3D-model of the dynamics and chemistry of the Martian atmosphere. DPS 2006 abstract #60.05. • Moudden and McConnell, 2006. Three-dimensional on-line chemical modeling in a Mars general circulation model, submitted to Icarus.
Mars photochemistry • Kinetic reactions • Photolysis Mars meteorology GEOS_CHEM Mars GCM LIDORT Table for J-values
Interface with MGCM • Time Year, Day-of-Year, Planetary-second • Constants CMN_GCTM: eccentricity, g0, Rplanet, Rd, Cp, … • File format Input Netcdf: 40 days / GCM met file Output Bpch: restart file & diagnostics • Grid structure 64 longitudes, 36 latitudes, 25 levels
120nm – 800nm CO2, O2, O3, H2O, Rayleigh, dust Phot-freq Cross section Quantum yield Actinic flux LIDORT→Actinic Flux→Integration→Scaling Look-up table for J: SZA, Surface P, O3 column, H2O column, column dust, T profile type