Parameterization of gaseous dry deposition in atmospheric chemistry models: Sensitivity to aerodynamic resistance formulations under statically stable conditions

Publication date: December 2016
Source:Atmospheric Environment, Volume 147
Author(s): Kenjiro Toyota, Ashu P. Dastoor, Andrei Ryzhkov
Turbulence controls the vertical transfer of momentum, heat and trace constituents in the atmospheric boundary layer. In the lowest 10% of this layer lies the surface boundary layer (SBL) where the vertical fluxes of transferred quantities have been successfully parameterized using the Monin-Obukhov similarity theory in weather forecast, climate and atmospheric chemistry models. However, there is a large degree of empiricism in the stability-correction parameterizations used to formulate eddy diffusivity and aerodynamic resistance particularly under strongly stable ambient conditions. Although the influence of uncertainties in stability-correction parameterizations on eddy diffusivity is actively studied in boundary-layer meteorological modeling, its impact on dry deposition in atmospheric chemistry modeling is not well characterized. In this study, we address this gap by providing the mathematical basis for the relationship between the formulations of vertical surface flux used in meteorological and atmospheric chemistry modeling communities, and by examining the sensitivity of the modeled dry deposition velocities in statically stable SBL to the choice of stability-correction parameterizations used in three operational and/or research environmental models (GEM/GEM-MACH, ECMWF IFS and CMAQ-MM5). Aerodynamic resistances (ra) calculated by the three sets of parameterizations are notably different from each other and are also different from those calculated by a "z-less" scaling formulation under strongly stable conditions (the bulk Richardson number > 0.2). Furthermore, we show that many atmospheric chemistry models calculate ra using formulations which are inconsistent with the derivation of micro-meteorological parameters. Finally, practical implications of the differences in stability-correction algorithms are discussed for the computations of dry deposition velocities of SO2, O3 and reactive bromine compounds for specific cases of stable SBL.

Graphical abstract

http://ift.tt/2dQDGEG