Modeling of dispersion of air pollutants in the atmosphere is one of the most important and challenging scientific problems. There are several natural and anthropogenic events where passive or chemically active compounds are emitted into the atmosphere. The effect of these chemical species can have serious impacts on our environment and human health. Modeling the dispersion of air pollutants can predict this effect. Therefore, development of various model strategies is a key element for the governmental and scientific communities. We provide here a brief review on the mathematical modeling of the dispersion of air pollutants in the atmosphere. We discuss the advantages and drawbacks of several model tools and strategies, namely Gaussian, Lagrangian, Eulerian and CFD models. We especially focus on several recent advances in this multidisciplinary research field, like parallel computing using graphical processing units, or adaptive mesh refinement.

During the GRAMINAE Integrated Experiment between 20 May and 15 June 2000, the ozone flux was measured by the eddy covariance method above intensively managed grassland in Braunschweig, northern Germany. Three different phases of vegetation were covered during the measuring campaign: tall grass canopy before cut (29 May 2000), short grass after cut, and re-growing vegetation after fertilization (5 June 2000). Results show that beside weather conditions, the agricultural activities significantly influenced the O₃ fluxes. After the cut the daytime average of the deposition velocity (v_d) decreased from 0.44 cm s⁻¹ to 0.26 cm s⁻¹ and increased again to 0.32 cm s⁻¹ during the third period. Detailed model calculations were carried out to estimate deposition velocity and ozone flux. The model captures the general diurnal patter of deposition, with v_d daytime values of 0.52, 0.24, and 0.35 cm s⁻¹ in the first, second and third period, respectively. Thus the model predicts a stronger response to the cut than the measurements, which is nevertheless smaller than expected on the basis of change in leaf area. The results show that both cut and fertilization have complex impacts on fluxes. Reduction of vegetation by cutting decreased the stomatal flux initially greatly, but the stomatal flux recovered to 80% of its original value within a week. At the same time, the non-stomatal flux appears to have increased directly after the cut, which the model partially explains by an increase in the deposition to the soil. A missing sink after the cut may be the chemical interaction with biogenic volatile organic compounds released after the cut and exposed senescent plant parts, or the increase in soil NO emissions after fertilization. Increased canopy temperatures may also have promoted ozone destruction on leaf surfaces. These results demonstrate the importance of canopy structure and non-stomatal pathways on O₃ fluxes.

Long-range simulation of the dispersion of air pollutants in the atmosphere is one of the most challenging tasks in geosciences. Application of precise and fast numerical models in risk management and decision support can save human lives and can diminish consequences of an accidental release. Disaster at Fukushima Daiichi nuclear power plant has been the most serious event in the nuclear technology and industry in the recent years. We present and discuss the results of the numerical simulations on dispersion of Fukushima-derived particulate ¹³¹I and ¹³⁷Cs using a global scale Lagrangian particle model. We compare concentrations and arrival times, using two emission scenarios, with the measured data obtained from 182 monitoring stations located all over the Northern Hemisphere. We also investigate the homogenization of isotopes in the atmosphere. Peak concentrations were predicted with typical accuracy of one order of magnitude showing a general underestimation in the case of ¹³¹I but not for ¹³⁷Cs. Tropical and Arctic plumes, as well as the early detections in American and European midlatitudes were generally well predicted, however, the later regional-scale mixing could not be captured by the model. Our investigation highlights the importance of the parameterization of free atmospheric turbulence.