Difference between revisions of "UKCA & UMUI Tutorial 7"

Revision as of 10:19, 24 June 2013

Adding Dry Deposition

UKCA uses two different dry-deposition schemes:

A simple 2D parameterisation described by Giannakopoulos (1999) [1], Ganzeveld and Lelieveld (1995)[2], and Sander and Crutzen 1996 (JGR 101 9121)[3].
A more detailed interactive parameterisation, based on the Wesely scheme (Wesely, 1989; Sanderson 2007)[4,5]

The default is to use the 2D scheme, although it is advisable to use the interactive scheme. Within the UKCA code, whether a species is dry deposited or not is controlled in the ukca_chem_scheme.F90 file. In the chch_defs_scheme array there are lines like

chch_t( 10,'HONO2     ',  1,'TR        ','          ',  1,  1,  0),  & !  10 DD: 7,WD: 4,       
chch_t( 11,'H2O2      ',  1,'TR        ','          ',  1,  1,  0),  & !  11 DD: 8,WD: 5,

Where the 1 in the 6th column turns on dry-deposition of that species.

References

Giannakopoulos, C., M. P. Chipperfield, K. S. Law, and J. A. Pyle (1999), Validation and intercomparison of wet and dry deposition schemes using 210Pb in a global three-dimensional off-line chemical transport model, J. Geophys. Res., 104(D19), 23761–23784, doi:10.1029/1999JD900392.
Ganzeveld, L., and J. Lelieveld (1995), Dry deposition parameterization in a chemistry general circulation model and its influence on the distribution of reactive trace gases, J. Geophys. Res., 100(D10), 20999–21012, doi:10.1029/95JD02266.
Sander, R., and P. J. Crutzen (1996), Model study indicating halogen activation and ozone destruction in polluted air masses transported to the sea, J. Geophys. Res., 101(D4), 9121–9138, doi:10.1029/95JD03793.
M.L. Wesely, Parameterization of surface resistances to gaseous dry deposition in regional-scale numerical models, Atmospheric Environment (1967), Volume 23, Issue 6, 1989, Pages 1293-1304, ISSN 0004-6981, http://dx.doi.org/10.1016/0004-6981(89)90153-4.
Sanderson, M. G., Collins, W. J., Hemming, D. L. and Betts, R. A. (2007), Stomatal conductance changes due to increasing carbon dioxide levels: Projected impact on surface ozone levels. Tellus B, 59: 404–411. doi: 10.1111/j.1600-0889.2007.00277.x

2D Dry Deposition Scheme

The deposition velocities for the 2D scheme are defined in the depvel_defs_scheme array, which is held in the ukca_chem_scheme.F90 module. This is a large array made up of size (6,5) blocks. These blocks mean

Summer (day) velocity over water	Summer (night) velocity over water	Summer (24h ave.) velocity over water	Winter (day) velocity over water	Winter (night) velocity over water	Winter (24h ave.) velocity over water
Summer (day) velocity over forest	Summer (night) velocity over forest	Summer (24h ave.) velocity over forest	Winter (day) velocity over forest	Winter (night) velocity over forest	Winter (24h ave.) velocity over forest
Summer (day) velocity over grass	Summer (night) velocity over grass	Summer (24h ave.) velocity over grass	Winter (day) velocity over grass	Winter (night) velocity over grass	Winter (24h ave.) velocity over grass
Summer (day) velocity over desert	Summer (night) velocity over desert	Summer (24h ave.) velocity over desert	Winter (day) velocity over desert	Winter (night) velocity over desert	Winter (24h ave.) velocity over desert
Summer (day) velocity over ice	Summer (night) velocity over ice	Summer (24h ave.) velocity over ice	Winter (day) velocity over ice	Winter (night) velocity over ice	Winter (24h ave.) velocity over ice

depvel_defs(6,5,jpdd)

see ASAD_MOD

6 = Number of time periods used in dry deposition i.e. summer(day,night,24h ave), winter(day,night,24h ave) - 24h aves not used

5 = Number of land use categories used in dry dep. (water,forest,grass,desert,ice) - desert not used

see ukca_ddeprt.F90 for allocation

Written by Luke Abraham 2013

@@ Line 8: / Line 8: @@
 * A more detailed interactive parameterisation, based on the Wesely scheme (Wesely, 1989; Sanderson 2007)[4,5]
-The default is to use the 2D scheme, although it is advisable to use the interactive scheme. Within the UKCA code, whether a species is dry deposited or not is controlled in the '''ukca_chem_<span style="color:blue">scheme</span>.F90''' file. In the '''chch_defs''' array there are lines like
+The default is to use the 2D scheme, although it is advisable to use the interactive scheme. Within the UKCA code, whether a species is dry deposited or not is controlled in the '''ukca_chem_<span style="color:blue">scheme</span>.F90''' file. In the '''chch_defs_<span style="color:blue">scheme</span>''' array there are lines like
  chch_t( 10,'HONO2     ',  1,'TR        ','          ',  <span style="color:red">'''1'''</span>,  1,  0),  & !  10 DD: 7,WD: 4,
  chch_t( 11,'H2O2      ',  1,'TR        ','          ',  <span style="color:red">'''1'''</span>,  1,  0),  & !  11 DD: 8,WD: 5,
+Where the <span style="color:red">'''1'''</span> in the 6th column turns on dry-deposition of that species.
 '''References'''
@@ Line 23: / Line 23: @@
 ==2D Dry Deposition Scheme==
+The deposition velocities for the 2D scheme are defined in the '''depvel_defs_<span style="color:blue">scheme</span>''' array, which is held in the ''ukca_chem_<span style="color:blue">scheme</span>.F90'' module. This is a large array made up of size (6,5) blocks. These blocks mean
+{| border="1"
+| Summer (day) velocity over water  || Summer (night) velocity over water || Summer (24h ave.) velocity over water || Winter (day) velocity over water || Winter (night) velocity over water || Winter (24h ave.) velocity over water
+|-
+| Summer (day) velocity over forest || Summer (night) velocity over forest || Summer (24h ave.) velocity over forest || Winter (day) velocity over forest || Winter (night) velocity over forest || Winter (24h ave.) velocity over forest
+|-
+| Summer (day) velocity over grass || Summer (night) velocity over grass || Summer (24h ave.) velocity over grass || Winter (day) velocity over grass || Winter (night) velocity over grass || Winter (24h ave.) velocity over grass
+|-
+| Summer (day) velocity over desert || Summer (night) velocity over desert || Summer (24h ave.) velocity over desert || Winter (day) velocity over desert || Winter (night) velocity over desert || Winter (24h ave.) velocity over desert
+|-
+| Summer (day) velocity over ice || Summer (night) velocity over ice || Summer (24h ave.) velocity over ice || Winter (day) velocity over ice || Winter (night) velocity over ice || Winter (24h ave.) velocity over ice
+|}
 depvel_defs(6,5,jpdd)