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.

Monitoring and describing the spatiotemporal variability in dust aerosols is crucial for understanding their multiple effects, related feedbacks, and impacts within the Earth system. This study describes the development of the ModIs Dust AeroSol (MIDAS) data set. MIDAS provides columnar daily dust optical depth (DOD) at 550 nm at a global scale and fine spatial resolution (0.1^∘ × 0.1^∘) over a 15-year period (2003–2017). This new data set combines quality filtered satellite aerosol optical depth (AOD) retrievals from MODIS-Aqua at swath level (Collection 6.1; Level 2), along with DOD-to-AOD ratios provided by the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) reanalysis to derive DOD on the MODIS native grid. The uncertainties of the MODIS AOD and MERRA-2 dust fraction, with respect to the AEronet RObotic NETwork (AERONET) and LIdar climatology of vertical Aerosol Structure for space-based lidar simulation (LIVAS), respectively, are taken into account for the estimation of the total DOD uncertainty. MERRA-2 dust fractions are in very good agreement with those of LIVAS across the dust belt in the tropical Atlantic Ocean and the Arabian Sea; the agreement degrades in North America and the Southern Hemisphere, where dust sources are smaller. MIDAS, MERRA-2, and LIVAS DODs strongly agree when it comes to annual and seasonal spatial patterns, with colocated global DOD averages of 0.033, 0.031, and 0.029, respectively; however, deviations in dust loading are evident and regionally dependent. Overall, MIDAS is well correlated with AERONET-derived DODs (R=0.89) and only shows a small positive bias (0.004 or 2.7 %). Among the major dust areas of the planet, the highest R values (>0.9) are found at sites of North Africa, the Middle East, and Asia. MIDAS expands, complements, and upgrades the existing observational capabilities of dust aerosols, and it is suitable for dust climatological studies, model evaluation, and data assimilation.

The regime of intense desert dust (DD) episodes over the broader Mediterranean Basin is studied for the period 2000–2007 at a complete spatial coverage. An objective and dynamic algorithm has been set up which uses daily measurements of various aerosol optical properties taken by different satellite databases, enabling the identification of DD episodes and their classification into strong and extreme ones. The algorithm's performance was tested against surface-based (in situ) particulate matter (PM) and (columnar) sun-photometric AERONET (AErosol RObotic NETwork) measurements from stations distributed across the Mediterranean. The comparisons have shown the reasonable ability of the algorithm to detect the DD episodes taking place within the study region. The largest disagreements with PM data were found in the western Mediterranean in summer, when African dust transport has a great vertical extent that cannot be satisfactorily captured by surface measurements. According to our results, DD episodes in the Mediterranean Basin are quite frequent (up to 11.4 episodes yr⁻¹), while there is a significant spatial and temporal variability in their frequency of occurrence and their intensity. Strong episodes occur more frequently in the western Mediterranean Basin, whilst extreme ones appear more frequently over central Mediterranean Sea areas. Apart from this longitudinal variation, there is a predominant latitudinal variability in both frequency and intensity, with decreasing values from south to north. A significant seasonal variation was also found for the frequency of DD episodes, with both strong and extreme episodes being more frequent during summer in the western Mediterranean Basin, but during spring in its central and eastern parts. In most cases (> 85%) the Mediterranean dust episodes last a bit longer than a day on average, although their duration can reach six days for strong episodes and four days for extreme episodes. A noticeable year-to-year variability was also found, especially for the frequency of the episodes.

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.

Cases of atmospheric circulation evolution favouring the occurrence of desert aerosol episodes (DAEs) over the broader Mediterranean region were investigated using an objective and dynamic algorithm, with daily satellite data for the period 2000–2013. After identifying strong and extreme DAEs, at a 1°×1° geographical‐cell level, 255 dust aerosol episode days (DAEDs) and 148 cases of consecutive DAEDs, namely desert aerosol episode cases (DAECs), are defined. For each DAEC, the evolution of the lower tropospheric circulation 1 and 2 days before, during the initiation and after the cessation of the DAEC, is considered. S‐mode factor analysis and k‐means cluster analysis are applied on mean sea‐level pressure and 700 hPa geopotential height fields obtained from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) Reanalysis Project, classifying the 148 cases of atmospheric‐circulation evolution into six homogeneous and discrete clusters. The mean intra‐annual variation of the DAECs reveals a primary maximum in May (18.9%), and their mean annual number is equal to 11.4 DAECs. On a seasonal basis, the highest percentage of the DAECs is found in spring (51.4%). Maximum duration of the DAECs is 7 days, with 58.8% lasting 1 day. Annually, the mean monthly number of the DAEs varies from 35.8 (September) to 58.0 (April). The western parts of the Mediterranean are affected by the DAEs when cyclonic conditions prevail in the western Mediterranean and northwestern Africa. In contrast, the central and eastern parts of the study region are affected by dust storms when a low‐pressure system in the central Mediterranean or central Europe and an anticyclone in the eastern Mediterranean prevail. As to the mean regional intensity (aerosol optical depth at 550 nm) the strong DAEs vary from 0.67 to 0.77, while the extremes vary from 1.14 to 2.06. Generally, strong DAEs are more frequent than extremes (in five out of six clusters).

Dust particles larger than 20 µm in diameter have been regularly observed to remain airborne during long-range transport. In this work, we modify the parameterization of the mineral dust cycle in the GOCART-AFWA dust scheme of WRFV4.2.1 to also include such coarse and giant particles, and we further discuss the underlying misrepresented physical mechanisms which hamper the model in reproducing adequately the transport of the coarse and giant mineral particles. The initial particle size distribution is constrained by observations over desert dust sources. Furthermore, the Stokes drag coefficient has been updated to account for realistic dust particle sizes (Re < 10⁵). The new code was applied to simulate dust transport over Cabo Verde in August 2015 (AER-D campaign). Model results are evaluated against airborne dust measurements and the CALIPSO-LIVAS pure dust product. The results show that the modeled lifetimes of the coarser particles are shorter than those observed. Several sensitivity runs are performed by reducing artificially the particles' settling velocities in order to compensate underrepresented mechanisms, such as the non-spherical aerodynamics, in the relevant parameterization schemes. Our simulations reveal that particles with diameters of 5.5–17 and 40–100 µm are better represented under the assumption of an 80 % reduction in the settling velocity (UR80), while particles with sizes ranging between 17 and 40 µm are better represented in a 60 % reduction in settling velocity (UR60) scenario. The overall statistical analysis indicates that the best agreement with airborne in situ measurements downwind (Cabo Verde) is achieved with a 40 % reduction in settling velocity (UR40). Moreover, the UR80 experiment improves the representation of the vertical structure of the dust layers as those are captured by the CALIPSO-LIVAS vertically resolved pure dust observations. The current study highlights the necessity of upgrading the existing model parameterization schemes of the dust life-cycle components towards improving the assessment of the dust-related impacts within the Earth–atmosphere system.