Articles | Volume 8, issue 11
Geosci. Model Dev., 8, 3497–3522, 2015
Geosci. Model Dev., 8, 3497–3522, 2015

Model description paper 03 Nov 2015

Model description paper | 03 Nov 2015

Construction of the SILAM Eulerian atmospheric dispersion model based on the advection algorithm of Michael Galperin

M. Sofiev1, J. Vira1, R. Kouznetsov1,2, M. Prank1, J. Soares1, and E. Genikhovich3 M. Sofiev et al.
  • 1Finnish Meteorological Institute, Helsinki, Finland
  • 2A.M. Obukhov Institute for atmospheric Physics, Moscow, Russia
  • 3Main Geophysical Observatory, St. Petersburg, Russia

Abstract. The paper presents the transport module of the System for Integrated modeLling of Atmospheric coMposition SILAM v.5 based on the advection algorithm of Michael Galperin. This advection routine, so far weakly presented in the international literature, is positively defined, stable at any Courant number, and efficient computationally. We present the rigorous description of its original version, along with several updates that improve its monotonicity and shape preservation, allowing for applications to long-living species in conditions of complex atmospheric flows. The scheme is connected with other parts of the model in a way that preserves the sub-grid mass distribution information that is a cornerstone of the advection algorithm. The other parts include the previously developed vertical diffusion algorithm combined with dry deposition, a meteorological pre-processor, and chemical transformation modules.

The quality of the advection routine is evaluated using a large set of tests. The original approach has been previously compared with several classic algorithms widely used in operational dispersion models. The basic tests were repeated for the updated scheme and extended with real-wind simulations and demanding global 2-D tests recently suggested in the literature, which allowed one to position the scheme with regard to sophisticated state-of-the-art approaches. The advection scheme performance was fully comparable with other algorithms, with a modest computational cost.

This work was the last project of Dr. Sci. Michael Galperin, who passed away on 18 March 2008.

Short summary
The paper presents a transport mechanism of SILAM CTM based on an algorithm of M. Galperin. We describe the original scheme and its updates needed for applications to long-living species, complex atmospheric flows, etc. The scheme is connected to vertical diffusion, chemical transformation and deposition algorithms. Quality of the advection routine is evaluated with a large set of tests, which showed performance fully comparable with state-of-the-art algorithms at much lower computational costs.