geoschem / CAM

Community Atmosphere Model
Other
0 stars 0 forks source link

Fix dry and wet deposition for GEOS-Chem non-gas species. #24

Closed jimmielin closed 1 year ago

jimmielin commented 1 year ago

This update restores dry and wet deposition for non-gas (aerosol) species in GEOS-Chem chemical mechanism.

This bug was likely introduced due to: 1) update in the species list leading to a new way to generate the species lists which inadvertently removed the non-gas species from <gas_wetdep_list> and <drydep_list>. (more on this later), and 2) update in the MCT -> NUOPC transition to remove the Fritz et al. approach of performing deposition for aerosols in a manner similar to HNO3, which was originally coded in seq_drydep_mod.F90 (in mct, now removed) and mo_neu_wetdep.F90 (changes removed for some reason)

This update restores: 1) the aerosol species into the dry and wet deposition lists. Note: for some reason, the aerosols currently go into the gas-phase lists. This was how it was always done; the repercussions are unknown but the aer_ lists are all MAM4-species, and the deposition is handled within aero_model. I suspect they won't work if they're put into the aerosol list, so we maintain the old approach. 2) <dep_data_file> which uses the appropriate parameters and coefficients for aerosol species in dry and wet deposition, and uses parameters from HNO3 if otherwise not available. Complimentary repository and commit to the updated <dep_data_file> generation is here: https://github.com/CESM-GC/drydep_henrys_law_nc/commit/a9f85e4abb4b274a533a768ebbe068487576c921

The results of this update have been verified to lead the model to much better simulate nitrate aerosol, and NOx levels compared to the offline GEOS-Chem simulation and better reproduce NOx levels compared to observations.

Global Mass: Time slice 2        LM_left: 1:18 (~966.03 to 212.62 hPa)       LM_right: 1:18 (~966.04 to 212.62 hPa)
Reference model: mdl_ref_gc_benchmark/mdl_ref_gc_benchmark_13.4.0rc.4-trop-201904.mat
Total Mass [Gg]        cesm2.3gc14.1.1.pnit.2016MAM         cesm2.3gc14.1.1.pnit.no_dep.2016MAM
Species                Left_Model               Right_Model             % diff (L-R)          Reference_Model
--------------------------------------------------------------------------------------------------------------
N2O                    2080863.322839       2080837.466562        +0.001            2020800.125000
N2O5                   41.383915        67.660355        -38.836 !!!        27.337704
NIT                    406.795531       1985.468125          -79.511 !!!        409.717499
NO                     64.227130        91.946825        -30.148 !!!        49.911560
NO2                    323.534085       433.519836       -25.370 !!!        303.669586
NO3                    11.976234        16.162134        -25.899 !!!        13.659160
NOx (dummy)            9.172012         12.486169        -26.543 !!!        -1.000000
NOy (dummy)            68.457629        102.500587       -33.212 !!!        -1.000000
O3                     368486.586416        417136.442371        -11.663 !!         297461.093750
OH                     0.220583         0.264629         -16.644 !!         0.187330
PAN                    2258.947760          2746.872443          -17.763 !!         1453.664185
SO2                    450.174673       420.840783        +6.970 !      324.507172
SO4                    1826.706777          1833.872111           -0.391            1379.099609