GEOS-Chem species: Wet deposition metadata

From Geos-chem
Revision as of 19:27, 12 May 2020 by Bmy (Talk | contribs) (Properties for gas-phase species)

Jump to: navigation, search

Previous | Next | Guide to species in GEOS-Chem

  1. Species in GEOS-Chem
  2. Henry's law parameters
  3. Dry deposition parameters
  4. Wet deposition parameters
  5. GEOS-Chem species database
  6. GEOS-Chem species units
  7. Adding passive species to GEOS-Chem
  8. Species indexing in GEOS-Chem


Overview

Each wet-depositing species in GEOS-Chem contains one or more of the following properties, depending on whether the species is a gas or aerosol:

Properties for aerosol species

Quantity Units Description
WD_Is_HNO3 1 Flag to determine if the species is HNO3.
  • Rainout/washout of HNO3 needs to be handled separately from other gas-phase species.
  • WD_Is_HNO3 will be automatically determined, depending on the value of Name.
WD_Is_SO2 1 Flag to determine if the species is SO2.
  • Rainout/washout of HNO3 separately from other gas-phase species.
  • WD_Is_SO2 will be automatically determined, depending on the value of Name.
WD_Is_H2SO4 1 Flag to determine if the species is H2SO4.
  • Rainout/washout of H2SO4 needs to be handled separately from other gas-phase species.
  • WD_Is_H2SO4 will be automatically determined, depending on the value of Name.
  • H2SO4 is a species in TOMAS aerosol microphysics only.
WD_CoarseAer 1 Set this flag to indicate that the species is a coarse aerosol (i.e. having a radius > 1 μm)
  • Washout of coarse aerosols is handled separately from the washout of fine aerosols.
WD_AerScavEff 1 Specifies the aerosol scavenging efficiency. This factor multiplies F, the fraction of aerosol species that is lost to convective updraft scavenging (cf Eq. 2 of Jacob et al 2000).
  • For most aerosol species, WD_AerScavEff = 1.0.
  • For secondary organic aerosol species, WD_AerScavEff = 0.8.
  • For hydrophobic species (e.g. MOPO), WD_AerScavEff = 0.0.

WD_AerScavEff is applied in routine COMPUTE_F of GEOS-Chem module wetscav_mod.F90.

WD_KcScaleFac 1 Specifies a temperature-dependent scale factor that is used to multiply K (aka Kc), the rate constant for conversion of cloud condensate to precipitation (cf Eq. 1 of Jacob et al, 2000).

WD_KcScaleFac is defined as a 3-element vector, where:

  • WD_KcScaleFac(1) multiplies K when T < 237 K.
  • WD_KcScaleFac(2) multiplies K when 237 K <= T < 258 K.
  • WD_KcScaleFac(3) multiplies K when T >= 258K.

This allows us to better simulate scavenging by snow, where, for most species, we need to halve the K rate when 237 K <= T < 258 K. This can be easily done by setting the value KcScaleFac(2) = 0.5.

WD_KcScaleFac is applied in routine F_AEROSOL of GEOS_Chem module GeosCore/wetscav_mod.F.

WD_RainoutEff 1 Specifies a temperature-dependent scale factor that is used to multiply Fi (aka RAINFRAC), the fraction of species scavenged by rainout (cf Eq. 1 of Jacob et al, 2000).

WD_RainoutEff is defined as a 3-element vector, where:

  • WD_RainoutEff(1) multiplies Fi when T < 237 K.
  • WD_RainoutEff(2) multiplies Fi when 237 K <= T < 258 K.
  • WD_RainoutEff(3) multiplies Fi when T >= 258K.

This allows us to better simulate scavenging by snow and impaction scavenging of BC. For most species, we need to be able to turn off rainout when 237 K <= T < 258 K. This can be easily done by setting RainoutEff(2) = 0.

WD_RainoutEff is applied in routine APPLY_RAINOUT_EFF of GEOS_Chem module GeosCore/wetscav_mod.F.

NOTE: For secondary organic aerosol species, the maximum value of WD_RainoutEff will be 0.8 instead of 1.0.

--Bob Yantosca (talk) 19:06, 12 May 2020 (UTC)

Properties for gas-phase species

Quantity Units Description
WD_LiqAndGas 1 True/False flag indicating whether the ice-to-gas ratio needs to be computed for by co-condensation

.

WD_ConvFacI2G 1 Specifies the conversion factor (i.e. the ratio of sticking coefficients on the ice surface) for computing the ice-to-gas-ratio by co-condensation (cf. Eq. 9 of Jacob et al, 2000).
  • WD_ConvFacI2G only needs to be specified for those species having WD_LiqAndGas = TRUE.
WD_RetFactor 1 Specifies the retention efficiency Ri of tracer in the liquid cloud condensate as it is converted to precipitation. (Ri < 1 accounts for volatization during riming.) The retention factor is defined in Eq. 1 of Jacob et al, 2000.

--Bob Yantosca (talk) 19:08, 12 May 2020 (UTC)

Table of wet deposition parameters

The table below shows the various wet deposition parameters for both gas-phase and aerosol species. The quantities in each column correspond to fields of the Species derived type, as described on our GEOS-Chem species database wiki page.

Co-condensation parameters

The species listed below have the ice-to-gas ratio computed by co-condensation (see Eq. 9 of Jacob et al, 2000).

Gas-phase Species Parameters
Species MW_g
[g mol-1]
EmMW_g
[g mol-1]
Molec
Ratio

[1]
WD_LiqAndGas WD_ConvFacI2G References
H2O2 34.2 - - YES 4.36564e-1 Jacob et al 2000
NH3 17.04 - - YES 6.17395e-1 Jacob et al 2000

All other parameters

The remaining wet deposition parameters are listed below:

Aerosol Species Parameters Gas-phase Species Parameters
Species MW_g
[g mol-1]
EmMW_g
[g mol-1]
Molec
Ratio

[1]
WD_Coarse
Aer

[?]
WD_Aer
ScavEff

[1]
WD_KcScaleFac
[1]
WD_RainoutEff
[1]
WD_Ret
Factor

[1]
References
T < 237 K 237 K <= T < 258 K T > 258 K T < 237 K 237 K <= T < 258 K T > 258 K
ACET 58.08 12 3 - - - - - - - - - -
ALD2 44.05 12 2 - - - - - - - - - -
ALK4 58.12 12 4 - - - - - - - - - -
ASOAN
ASOA1
ASOA2
ASOA3
150 150 - NO 0.8 1.0 0.5 1.0 0.8 0.0 0.8 - -
ASOG1
ASOG2
ASOG3
150 150 - - - - - - - - - 2.0e-2 -
Be7 201 201 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
Br2 160 160 - - - - - - - - - 0.0 -
BCPI 12.01 12 1 NO 1.0 1.0 1.0 0.5 1.0 0.0 1.0 - -
BCPO 12.01 12 1 NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
CH2O 30 30 - - - - - - - - - 2.0e-2 Jacob et al [2000]
DST1
DSTAL1
NITD1
SO4D1
29 29 - NO 1.0 1.0 1.0 1.0 1.0 1.0 1.0 - -
DST2
DSTAL2
NITD2
SO4D2
29 29 - YES 1.0 1.0 1.0 1.0 1.0 1.0 1.0 - -
DST3
DSTAL3
NITD3
SO4D3
29 29 - YES 1.0 1.0 1.0 1.0 1.0 1.0 1.0 - -
DST4
DSTAL4
NITD4
SO4D4
29 29 - YES 1.0 1.0 1.0 1.0 1.0 1.0 1.0 - -
GLYC 60 60 - - - - - - - - - 2.0e-2 -
GLYX 58 58 - - - - - - - - - - -
H2O2 34 34 - - - - - - - - - 5.00e-2 Jacob et al [2000]
HAC 74 74 - - - - - - - - - - -
HBr 81 81 - - - - - - - - - 1.0 -
HCl 36 36 - - - - - - - - - 1.0 -
Hg2 201 201 - - - - - - - - - 0.0 -
HgP 201 201 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
HNO3 63 63 - - - - - - - - - - -
HOBr 97 97 - - - - - - - - - 0.0 -
IEPOX 118 118 - - - - - - - - - 2.0e-2 -
ISOA1
ISOA2
ISOA3
150 150 - NO 0.8 1.0 0.5 1.0 0.8 0.0 0.8 - -
ISOG1
ISOG2
ISOG3
150 150 - - - - - - - - - 2.0e-2 -
ISOPN 147 147 - - - - - - - - - 2.0e-2 -
LIMO
MTPA
MTPO
136.23 136.23 - - - - - - - - - 2.0e-2 -
MAP 76 76 - - - - - - - - - 2.0e-2 -
MMN 149 149 - - - - - - - - - 2.0e-2 -
MOBA 114 114 - - - - - - - - - 2.0e-2 -
MOPI 12.01 12 1 NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
MOPO 12.01 12 1 NO 0.0 1.0 1.0 1.0 0.0 0.0 0.0 - -
MSA 96 96 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
NH3 17 17 - - - - - - - - - 5.0e-2 Jacob et al [2000]?
NH4 18 18 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
NIT 62 62 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
NITs 31.4 31.4 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
OCPI
POA1
12.01 12 1 NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
OCPO
POA2
12.01 12 1 NO 1.0 0.5 0.5 0.5 1.0 0.0 1.0 - -
OPOA1
OPOA2
12.01 12 1 NO 0.8 1.0 0.5 1.0 0.8 0.0 0.8 - -
OPOG1
OPOG2
12.01 12 1 - - - - - - - - 2.0e-2 -
Pb210 210 210 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
POG1
POG2
12.01 12 1 - - - - - - - - 2.0e-2 -
POPPBCPI
(PHE)
178.23 178.23 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
POPPBCPI
(PYR)
202.25 202.25 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
POPPBCPI
(BaP)
252.31 252.31 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
POPPBCPO
(PHE)
178.23 178.23 - NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
POPPBCPO
(PYR)
202.25 202.25 - NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
POPPBCPO
(BaP)
252.31 252.31 - NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
POPPOCPI
(PHE)
178.23 178.23 - NO 1.0 1.0 0.5 1.0 1.0 1.0 1.0 - -
POPPOCPI
(PHE)
202.25 202.25 - NO 1.0 1.0 0.5 1.0 1.0 1.0 1.0 - -
POPPOCPI (PHE) 252.31 252.31 - NO 1.0 1.0 0.5 1.0 1.0 1.0 1.0 - -
POPOCPO
(PHE)
178.23 178.23 - NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
POPOCPO
(PYR)
202.25 202.25 - NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
POPOCPO
(BaP)
252.31 252.31 - NO 1.0 1.0 1.0 0.5 1.0 1.0 0.0 - -
POPG (PHE) 178 178 - - - - - - - - - 0.0 -
POPG (PYR) 202 178 - - - - - - - - - 0.0 -
POPG (BaP) 252 178 - - - - - - - - - 0.0 -
PROPNN 119 119 - - - - - - - - - 2.0e-2 -
RIP 118 118 - - - - - - - - - 2.0e-2 -
SALA 31.4 31.4 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
SALC 31.4 31.4 - YES 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
SO2 64 64 - - - - - - - - - - -
SO4 96 96 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
SO4s 31.4 31.4 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
TSOA0
TSOA1
TSOA2
TSOA3
150 150 - NO 0.8 1.0 0.5 1.0 0.8 0.0 0.8 - -
TSOG0
TSOG1
TSOG2
TSOG3
150 150 - - - - - - - - - 2.0e-2 -
TOMAS microphysics species Aerosol Species Parameters Gas-phase Species Parameters
DUST1-
DUST40
100 100 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
ECIL1-
ECIL40
96 96 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
NK1-NK40 1 1 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
OCIL1-
OCIL40
12 12 1 NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
SF1-
SF40
96 96 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -
SS1-
SS40
58.5 58.5 - NO 1.0 1.0 0.5 1.0 1.0 0.0 1.0 - -

NOTES:

  1. We follow the algorithm of Qiaoqiao Wang et al (2014) for scavenging by snow and impaction scavenging for BC and IN, namely:
    • When 237 K <= T < 258 K:
      • Multiply Kc—the rate at which cloud condensate changes into precipitation— by 0.5, and
      • Allow black carbon, dust, HNO3, 210Pb, and 7Bethe following species, which are considered to be IN, to be rained out (aka "cold cloud scavenging").
      • Do not allow any other soluble aerosol species to rain out (i.e. set the rainout fraction to zero).

--Bob Y. (talk) 20:34, 23 October 2015 (UTC)


Previous | Next | Guide to species in GEOS-Chem