Online coupled mesoscale meteorology atmospheric chemistry models have undergone a rapid evolution in recent years. Although mainly developed by the air quality modelling community, these models are also of interest for numerical weather prediction and regional climate modelling as they can consider not only the effects of meteorology on air quality, but also the potentially important effects of atmospheric composition on weather. Two ways of online coupling can be distinguished: online integrated and online access coupling. Online integrated models simulate meteorology and chemistry over the same grid in one model using one main time step for integration. Online access models use independent meteorology and chemistry modules that might even have different grids, but exchange meteorology and chemistry data on a regular and frequent basis. This article offers a comprehensive review of the current research status of online coupled meteorology and atmospheric chemistry modelling within Europe. Eighteen regional online coupled models developed or being used in Europe are described and compared. Topics discussed include a survey of processes relevant to the interactions between atmospheric physics, dynamics and composition; a brief overview of existing online mesoscale models and European model developments; an analysis on how feedback processes are treated in these models; numerical issues associated with coupled models; and several case studies and model performance evaluation methods. Finally, this article highlights selected scientific issues and emerging challenges that require proper consideration to improve the reliability and usability of these models for the three scientific communities: air quality, numerical meteorology modelling (including weather prediction) and climate modelling. This review will be of particular interest to model developers and users in all three fields as it presents a synthesis of scientific progress and provides recommendations for future research directions and priorities in the development, application and evaluation of online coupled models.

During the months of March and April 2020 we witnessed the largest-scale experiment in history in terms of air quality in cities. Any prediction of this experiment's results may be obvious to science, as it was totally expected, the air quality has improved substantially. Simply stated, it comes as no surprise. The lockdown has made it possible to quantify the limit of decrease in pollution in light of this drastic reduction in traffic, in Madrid and Barcelona showed a significant decrease of the order of 75%. In the case of Spain's two largest cities, the reductions of NO2 concentrations were 62% and 50%, respectively. Hourly measurements were obtained from 24 and 9 air quality stations from the monitoring networks during the month of March 2020. These results allow us to see the limits that can be achieved by implementing low emission zones (LEZ), as well as the amount of contamination that must be eliminated, which in the cases of Madrid and Barcelona, represent 55%. This value defines the levels of effort and scope of actions to be taken in order to ensure that both cities achieve a clean and healthy atmosphere in terms of NO2.

We demonstrate improvements in CALIPSO (Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations) dust extinction retrievals over northern Africa and Europe when corrections are applied regarding the Saharan dust lidar ratio assumption, the separation of the dust portion in detected dust mixtures, and the averaging scheme introduced in the Level 3 CALIPSO product. First, a universal, spatially constant lidar ratio of 58 sr instead of 40 sr is applied to individual Level 2 dust-related backscatter products. The resulting aerosol optical depths show an improvement compared with synchronous and collocated AERONET (Aerosol Robotic Network) measurements. An absolute bias of the order of −0.03 has been found, improving on the statistically significant biases of the order of −0.10 reported in the literature for the original CALIPSO product. When compared with the MODIS (Moderate-Resolution Imaging Spectroradiometer) collocated aerosol optical depth (AOD) product, the CALIPSO negative bias is even less for the lidar ratio of 58 sr. After introducing the new lidar ratio for the domain studied, we examine potential improvements to the climatological CALIPSO Level 3 extinction product: (1) by introducing a new methodology for the calculation of pure dust extinction from dust mixtures and (2) by applying an averaging scheme that includes zero extinction values for the nondust aerosol types detected. The scheme is applied at a horizontal spatial resolution of 1° × 1° for ease of comparison with the instantaneous and collocated dust extinction profiles simulated by the BSC-DREAM8b dust model. Comparisons show that the extinction profiles retrieved with the proposed methodology reproduce the well-known model biases per subregion examined. The very good agreement of the proposed CALIPSO extinction product with respect to AERONET, MODIS and the BSC-DREAM8b dust model makes this dataset an ideal candidate for the provision of an accurate and robust multiyear dust climatology over northern Africa and Europe.

The EURODELTA III exercise has facilitated a comprehensive intercomparison and evaluation of chemistry transport model performances. Participating models performed calculations for four 1-month periods in different seasons in the years 2006 to 2009, allowing the influence of different meteorological conditions on model performances to be evaluated. The exercise was performed with strict requirements for the input data, with few exceptions. As a consequence, most of differences in the outputs will be attributed to the differences in model formulations of chemical and physical processes. The models were evaluated mainly for background rural stations in Europe. The performance was assessed in terms of bias, root mean square error and correlation with respect to the concentrations of air pollutants (NO₂, O₃, SO₂, PM₁₀ and PM_2.5), as well as key meteorological variables. Though most of meteorological parameters were prescribed, some variables like the planetary boundary layer (PBL) height and the vertical diffusion coefficient were derived in the model preprocessors and can partly explain the spread in model results. In general, the daytime PBL height is underestimated by all models. The largest variability of predicted PBL is observed over the ocean and seas. For ozone, this study shows the importance of proper boundary conditions for accurate model calculations and then on the regime of the gas and particle chemistry. The models show similar and quite good performance for nitrogen dioxide, whereas they struggle to accurately reproduce measured sulfur dioxide concentrations (for which the agreement with observations is the poorest). In general, the models provide a close-to-observations map of particulate matter (PM_2.5 and PM₁₀) concentrations over Europe rather with correlations in the range 0.4–0.7 and a systematic underestimation reaching −10 µg m⁻³ for PM₁₀. The highest concentrations are much more underestimated, particularly in wintertime. Further evaluation of the mean diurnal cycles of PM reveals a general model tendency to overestimate the effect of the PBL height rise on PM levels in the morning, while the intensity of afternoon chemistry leads formation of secondary species to be underestimated. This results in larger modelled PM diurnal variations than the observations for all seasons. The models tend to be too sensitive to the daily variation of the PBL. All in all, in most cases model performances are more influenced by the model setup than the season. The good representation of temporal evolution of wind speed is the most responsible for models' skillfulness in reproducing the daily variability of pollutant concentrations (e.g. the development of peak episodes), while the reconstruction of the PBL diurnal cycle seems to play a larger role in driving the corresponding pollutant diurnal cycle and hence determines the presence of systematic positive and negative biases detectable on daily basis.

The main aim of the present study is to describe the vertical structure of the intense Mediterranean dust outbreaks, based on the use of satellite and surface-based retrievals/measurements. Strong and extreme desert dust (DD) episodes are identified at 1° × 1° spatial resolution, over the period March 2000–February 2013, through the implementation of an updated objective and dynamic algorithm. According to the algorithm, strong DD episodes occurring at a specific place correspond to cases in which the daily aerosol optical depth at 550 nm (AOD_550 nm) exceeds or equals the long-term mean AOD_550 nm (Mean) plus two standard deviations (SD), which is also smaller than Mean+4 × SD. Extreme DD episodes correspond to cases in which the daily AOD_550 nm value equals or exceeds Mean+4 × SD. For the identification of DD episodes, additional optical properties (Ångström exponent, fine fraction, effective radius and aerosol index) derived by the MODIS-Terra & Aqua (also AOD retrievals), OMI-Aura and EP-TOMS databases are used as inputs. According to the algorithm using MODIS-Terra data, over the period March 2000–February 2013, strong DD episodes occur more frequently (up to 9.9 episodes year⁻¹) over the western Mediterranean, while the corresponding frequencies for the extreme ones are smaller (up to 3.3 episodes year⁻¹, central Mediterranean Sea). In contrast to their frequency, dust episodes are more intense (AODs up to 4.1), over the central and eastern Mediterranean Sea, off the northern African coasts. Slightly lower frequencies and higher intensities are found when the satellite algorithm operates based on MODIS-Aqua retrievals, for the period 2003–2012. The consistency of the algorithm is successfully tested through the application of an alternative methodology for the determination of DD episodes, which produced similar features of the episodes' frequency and intensity, with just slightly higher frequencies and lower intensities. The performance of the satellite algorithm is assessed against surface-based daily data from 109 sun-photometric (AERONET) and 22 PM₁₀ stations. The agreement between AERONET and MODIS AOD is satisfactory (R = 0.505 − 0.750) and improves considerably when MODIS level 3 retrievals with higher sub-grid spatial representativeness and homogeneity are considered. Through the comparison against PM₁₀ concentrations, it is found that the presence of dust is justified in all ground stations with success scores ranging from 68 to 97 %. However, poor agreement is evident between satellite and ground PM₁₀ observations in the western parts of the Mediterranean, which is attributed to the desert dust outbreaks' vertical extension and the high altitude of dust presence. The CALIOP vertical profiles of pure and polluted dust observations and the associated total backscatter coefficient at 532 nm (β_532 nm), indicate that dust particles are mainly detected between 0.5 and 6 km, though they can reach 8 km between the parallels 32 and 38° N in warm seasons. An increased number of CALIOP dust records at higher altitudes is observed with increased latitude, northwards to 40° N, revealing an ascending mode of the dust transport. However, the overall intensity of DD episodes is maximum (up to 0.006 km⁻¹ sr⁻¹) below 2 km and at the southern parts of the study region (30–34° N). Additionally, the average thickness of dust layers gradually decreases from 4 to 2 km, moving from south to north. In spring, dust layers of moderate-to-high β_532 nm values ( ∼ 0.004 km⁻¹ sr⁻¹) are detected over the Mediterranean (35–42° N), extending from 2 to 4 km. Over the western Mediterranean, dust layers are observed between 2 and 6 km, while their base height is decreased down to 0.5 km for increasing longitudes underlying the role of topography and thermal convection. The vertical profiles of CALIOP β_532 nm confirm the multilayered structure of the Mediterranean desert dust outbreaks on both annual and seasonal bases, with several dust layers of variable geometrical characteristics and intensities. A detailed analysis of the vertical structure of specific DD episodes using CALIOP profiles reveals that the consideration of the dust vertical structure is necessary when attempting comparisons between columnar MODIS AOD retrievals and ground PM₁₀ concentrations.

One of the major sources of uncertainty in model estimates of the global sea-salt aerosol distribution is the emission parameterization. We evaluate a new sea-salt aerosol life cycle module coupled to the online multiscale chemical transport model NMMB/BSC-CTM. We compare 5 yr global simulations using five state-of-the-art sea-salt open-ocean emission schemes with monthly averaged coarse aerosol optical depth (AOD) from selected AERONET sun photometers, surface concentration measurements from the University of Miami's Ocean Aerosol Network, and measurements from two NOAA/PMEL cruises (AEROINDOEX and ACE1). Model results are highly sensitive to the introduction of sea-surface-temperature (SST)-dependent emissions and to the accounting of spume particles production. Emission ranges from 3888 Tg yr⁻¹ to 8114 Tg yr⁻¹, lifetime varies between 7.3 h and 11.3 h, and the average column mass load is between 5.0 Tg and 7.2 Tg. Coarse AOD is reproduced with an overall correlation of around 0.5 and with normalized biases ranging from +8.8% to +38.8%. Surface concentration is simulated with normalized biases ranging from −9.5% to +28% and the overall correlation is around 0.5. Our results indicate that SST-dependent emission schemes improve the overall model performance in reproducing surface concentrations. On the other hand, they lead to an overestimation of the coarse AOD at tropical latitudes, although it may be affected by uncertainties in the comparison due to the use of all-sky model AOD, the treatment of water uptake, deposition and optical properties in the model and/or an inaccurate size distribution at emission.