Abstract:The simulations of hourly Radon 222 ( 222 Rn) concentrations are performed with the Unified EMEP model (Simpson et al., 2003) in order to validate different parameterization schemes for vertical mixing. In addition to the recently evaluated (Jeričević et al., 2010) operational EMEP vertical diffusion schemes K(z), the non-local O'Brien (1970) and local Blackadar (1979) schemes, as well as the nonlocal Grisogono scheme (e.g. Grisogono and Oerlemans, 2002), a new scheme which is local in stable boundary layer (SBL) and non-local in convective boundary layer (CBL) and based on total turbulent energy (TTE) closure (e.g., Mauritsen et. al, 2007)
INTRODUCTIONThe atmospheric boundary layer turbulence is the most important mechanism for the distribution of tracers. The parameterization of turbulent diffusion K(z) is an inevitable, traditional approach in description of turbulent processes and the estimation of turbulence effects in air quality models. Previous studies have already shown that the parameterizations of K(z) have significant impacts on simulated chemical concentrations (e.g. Oliviè et al., 2004). Various parameterizations, mainly first-order, non-local eddy diffusivity K schemes are proposed and widely used in practical applications (e.g. O'Brien, 1970;Holtslag and Moeng, 1991;Grisogono, 1995). A several modifications of the first-order schemes are proposed to overcome deficiencies to accurately simulate dispersion in different atmospheric stability conditions (e.g. Grisogono and Oerlemans, 2002;Mihailovic and Alapaty, 2007). In this paper, we evaluate the performance of the higher-order closure K scheme based on total turbulent energy (TTE) closure (e.g., Mauritsen et. al, 2007) in addition to recently evaluated schemes, the O'Brien, Blackadar and Grisogono schemes in the EMEP model (Jeričević et al. 2010). For the model evaluation available measurements 222 Rn from Europe are used. Radon is a radioactive gas which is found naturally in trace amounts in most rocks and soils. Since radionuclide 222 Rn has a half-life of 3.8 days and it is emitted primarily from the continents at a fairly constant emission rate between 0.8 and 1.3 atom cm -2 s -1 (Dentener et al., 1999) it is ideal to study the model sub-grid mixing schemes, numerical advection schemes or to compare different models. A considerable number of global and regional studies have been devoted to the simulation of