Quantification of the dust optical depth across spatiotemporal scales with the MIDAS global dataset (2003–2017)

Quantifying the dust optical depth (DOD) and its uncertainty across spatiotemporal scales is key to understanding and constraining the dust cycle and its interactions with the Earth System. This study quantifies the DOD along with its monthly and year-to-year variability between 2003 and 2017 at global and regional levels based on the MIDAS (ModIs Dust AeroSol) dataset, which combines Moderate Resolution Imaging Spectroradiometer (MODIS)-Aqua retrievals and Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), reanalysis products. We also describe the annual and seasonal geographical distributions of DOD across the main dust source regions and transport pathways. MIDAS provides columnar mid-visible (550 nm) DOD at fine spatial resolution (0.1× 0.1), expanding the current observational capabilities for monitoring the highly variable spatiotemporal features of the dust burden. We obtain a global DOD of 0.032±0.003 – approximately a quarter (23.4%± 2.4%) of the global aerosol optical depth (AOD) – with about 1 order of magnitude more DOD in the Northern Hemisphere (0.056± 0.004; 31.8%± 2.7%) than in the Southern Hemisphere (0.008±0.001; 8.2%±1.1%) and about 3.5 times more DOD over land (0.070±0.005) than over ocean (0.019± 0.002). The Northern Hemisphere monthly DOD is highly correlated with the corresponding monthly AOD (R2 = 0.94) and contributes 20 % to 48 % of it, both indicating a dominant dust contribution. In contrast, the contribution of dust to the monthly AOD does not exceed 17 % in the Southern Hemisphere, although the uncertainty in this region is larger. Among the major dust sources of the planet, the maximum DODs (∼ 1.2) are recorded in the Bodélé Depression of the northern Lake Chad Basin, whereas moderate-to-high intensities are encountered in the Western Sahara (boreal summer), along the eastern parts of the Middle East (boreal summer) and in the Taklamakan Desert (spring). Over oceans, major long-range dust transport is observed primarily along the tropical Atlantic (intensified during boreal summer) and secondarily in the North Pacific (intensified during boreal spring). Our calculated global and regional averages and associated uncertainties are consistent with some but not all recent observation-based studies. Our work provides a simple yet flexible method to estimate consisPublished by Copernicus Publications on behalf of the European Geosciences Union. 3554 A. Gkikas et al.: Quantification of the dust optical depth across spatiotemporal scales with MIDAS tent uncertainties across spatiotemporal scales, which will enhance the use of the MIDAS dataset in a variety of future studies.


Introduction
Mineral dust particles are emitted throughout the year across the arid and semi-arid regions of the planet, when winds exceed a threshold velocity mainly determined by soil texture, soil moisture and surface roughness. While dust aerosols have mainly a natural origin, the contribution of anthropogenic land use is estimated to be between 10 % and 25 % (Tegen et al., 2004;Stanelle et al., 2014;Ginoux et al., 2012). Dust is mobilized by microscale-to synoptic-scale phenomena, from dust devils developed under strong surface heating (Koch and Renno, 2005) to haboobs formed by intense cold-pool downdrafts related to deep moist convection (Knippertz et al., 2007) to synoptic patterns associated with intensified pressure gradients (Klose et al., 2010) and low-level jets (LLJs; Fiedler et al., 2013). Meteorology also plays a key role in the dust transport over maritime areas taking place mainly across the tropical Atlantic Ocean (Prospero and Mayol-Bracero, 2013;Yu et al., 2015), the northern Pacific Ocean (Husar et al., 2001), the Mediterranean (Flaounas et al., 2015;, the Arabian Sea (Ramaswamy et al., 2017) and the southern Atlantic Ocean (Gassó and Stein, 2007). Dust perturbs the radiation budget through direct (Sokolik and Toon, 1996), semi-direct (Huang et al., 2006) and indirect (Haywood and Bucher, 2000) processes, leading to impacts upon weather (Pérez et al., 2006;Gkikas et al., 2018Gkikas et al., , 2019 and climate (Lambert et al., 2013;Nabat et al., 2015). Upon deposition, nutrient-rich dust particles can increase the productivity of oceanic waters (Jickells et al., 2005) and terrestrial ecosystems (Okin et al., 2004) and perturb the carbon cycle (Jickells et al., 2014). Dust has been associated with epidemics of meningococcal meningitis in the African Sahel (Pérez García-Pando et al., 2014a, b) and with air quality degradation in urban areas (Kanakidou et al., 2011) causing respiratory (Kanatani et al., 2010) and cardiovascular (Du et al., 2016) disease when the population is exposed to high dust concentrations (Querol et al., 2019). Other socio-economic sectors can be regionally affected by dust storms (Middleton, 2017), including transportation (Weinzierl et al., 2012), agriculture (Stefanski and Sivakumar, 2009) and solar energy production (Kosmopoulos et al., 2018).
Satellite measurements and numerical simulations have repeatedly shown the remarkable contrast in dust load between the two hemispheres. The substantially higher dust load in the Northern Hemisphere is associated with the wider deserts extending across the so-called dust belt (Prospero et al., 2002;Ginoux et al., 2012), in contrast to the smaller sources in Australia, South Africa and South America. At a global scale, most of the entrained dust loads in the atmosphere originate from tropical and sub-tropical arid regions; yet, it is estimated that up to 5 % of the global dust budget consists of particles emitted from high-latitude sources (Bullard and Austin, 2011;Bullard et al., 2016). Given the key role of dust aerosols in the Earth system, it is imperative to monitor and understand the global dust cycle along with its multi-scale spatiotemporal variability over long time periods and fine spatial resolution. This task can be fulfilled to a certain degree using contemporary satellite instruments providing accurate retrievals and global coverage over extended time periods. With this approach, one of the key challenges is to discriminate dust from other aerosols. Several studies have combined AOD and aerosol index (AI) (e.g., Middleton and Goudie, 2001;Prospero et al., 2002) or AOD, single scattering albedo (SSA), and Ångström exponent (AE) (Ginoux et al., 2012) to identify the most active dust sources worldwide. Other studies have focused on the dust load and its variability in specific regions such as the Atlantic Ocean and the Arabian Sea (Peyridieu et al., 2013), the Sistan Basin (Rashki et al., 2015), the Mediterranean (Gkikas et al., 2016), Europe and North Africa (Marinou et al., 2017), and east Asia , among others. Liu et al. (2008) described the three-dimensional structure of dust aerosols at a global scale based on Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) vertically resolved retrievals acquired during the first operational year of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite mission. A more advanced approach has been introduced by Amiridis et al. (2013) and Marinou et al. (2017), who applied a more realistic lidar ratio for the Saharan dust and a series of quality filters on the CALIOP vertical profiles, in order to provide information about the vertical structure of dust layers at a global scale and coarse resolution in the LIVAS dataset (Amiridis et al., 2015). Ridley et al. (2016) quantified the global average DOD and its uncertainty for the period 2004-2008 based on AOD retrievals from passive spaceborne sensors (Moderate Resolution Imaging Spectroradiometer, MODIS; Multi-angle Imaging SpectroRadiometer, MISR), groundbased (Aerosol Robotic Network, AERONET) and shipborne (Maritime Aerosol Network, MAN) measurements from sunphotometers, and numerical simulations. Voss and Evan (2020) provided a long-term DOD climatology over the tropics and mid-latitudes at a coarse spatial resolution (1 • × 1 • ) based on MODIS and Advanced Very-High-Resolution Radiometer (AVHRR) observations, where DOD was estimated based on AOD, SSA and AE over land following Ginoux et al. (2012) and AOD, fine and coarse AOD (AERONET) and MERRA-2 winds over ocean. Based on vertically resolved CALIOP retrievals and columnar MODIS optical properties, Song et al. (2021) provided a long-term 4D global dust optical depth dataset, excluding the polar regions, over the period 2007-2019. In their approach, they took advantage of spaceborne observations that can be used for the discrimination/identification of dust aerosols characterized by their aspherical shape, coarse size and absorption.
Our study provides a global and regional quantification and description of the DOD based on the new ModIs Dust AeroSol (MIDAS) dataset . The powerful and innovative elements of the MIDAS DOD dataset are the (i) daily availability and fine spatial resolution (0.1 • × 0.1 • ), (ii) full global coverage including the sources and downwind areas (both over land and sea), (iii) 15-year temporal range (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017) using the most updated MODIS data collection, and (iv) grid-cell-level uncertainty quantification. In this contribution, we first describe the annual and seasonal geographical distribution of DOD across the main dust source regions and transport pathways (Sect. 4.1). We then quantify the average DOD and its monthly and year-toyear variability at global, hemispherical and regional levels, along with its fractional contribution to the AOD (Sect. 4.2). We summarize the main findings in Sect. 5.

ModIs Dust AeroSol (MIDAS) dataset
Our study is based on the MIDAS global fine-resolution dataset described in detail in Gkikas et al. (2021). We analyze the DOD at 550 nm and at 0.1 • × 0.1 • spatial resolution between 2003 and 2017. The MIDAS DOD results from the combination of the quality-filtered MODIS aerosol optical depth (AOD, Collection 6.1, Level 2; Levy et al., 2013) and the MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, version 2; Gelaro et al., 2017) fraction of AOD that is due to dust (MERRA-2 dust fraction, MDF). In Gkikas et al. (2021), the MDF was evaluated against the dust fraction obtained from quality-assured dust and non-dust CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization; Winker et al., 2009) profiles, available from the LIVAS database (Amiridis et al., 2015;Marinou et al., 2017;Proestakis et al., 2018). The MDF compares well with the LIVAS dust fraction over the dust-abundant areas extending across the NH dust belt, with maximum underestimations of 10 % in Asian deserts. The agreement is more limited in North America and the Southern Hemisphere (Figs. 1 and 2 in Gkikas et al., 2021). Overall, the MIDAS DOD is well correlated with AERONET dust-dominant retrievals (R = 0.89 at a global scale), and the absolute biases are mainly below 0.12 at stations near sources (Figs. 3 and 4 in Gkikas et al., 2021). The MIDAS DOD dataset was further verified against the LIVAS DOD and compared with MERRA-2 DODs (Fig. 5 in Gkikas et al., 2021). Among the three datasets, there is good agreement on the monthly vari-ability of the global and hemispherical DODs as well as on their long-term averages ( Fig. 6 and Table 1 in Gkikas et al., 2021). Moreover, the annual and seasonal DOD patterns are broadly similar in the three datasets throughout the period 2007-2015. Nevertheless, regional differences are found due to the different techniques (passive and active remote sensing, numerical simulations) applied for the DOD derivation.

Spatiotemporal averaging and propagation of grid-cell-level uncertainties
In Sect. 4.2 we provide DOD estimates that are averaged in space (regionally and globally) and in time (over months, seasons and years) along with their respective uncertainties. Averaging is performed according to the upper branch of Fig. 5 in Levy et al. (2009); i.e., spatial averaging is performed after grid cell temporal averaging for any of the timescales considered. The uncertainties of the DOD averages at the different spatiotemporal scales are based on the propagation of the daily grid cell uncertainties provided within the MIDAS dataset and presented in Gkikas et al. (2021). In short, the daily grid cell uncertainties combine the uncertainties of the MODIS AOD and the MERRA-2 MDF with respect to AERONET and LIVAS, respectively. The former is based on linear equations expressing the uncertainty with respect to AERONET AOD over ocean (Levy et al., 2013) and land (Levy et al., 2010;Sayer et al., 2013) with updated coefficients for C061 data depending on vegetated and arid surface types (see Eqs. 4 to 7 in Gkikas et al., 2021). The latter is based on a quartic (fourth degree) polynomial equation expressing the uncertainty with respect to the LIVAS dust fraction (see Eq. 8 in Gkikas et al., 2021). In order to estimate the uncertainties of the spatiotemporal averages we first assume that each of the daily grid cell uncertainties are composed of (1) a fraction that is completely random in time and space, (2) a fraction that is systematic (correlated) in time and random in space, and (3) a fraction that is systematic (correlated) in space and random in time. Our framework also assumes that the fraction of the daily grid cell uncertainty that is correlated both in space and time, for instance an instrument bias, is very small and therefore neglected. Under this framework, the propagation of uncertainty fraction (1) is negligible across the spatiotemporal scales considered, the propagation of uncertainty fraction (2) depends upon the size of the domain considered but is negligible at a global scale and across most of the regional domains considered in this study, and the propagation of fraction (3) accounts for most of the total average uncertainty. Since we cannot know fractions (1), (2), and (3) and (1) and (2) are negligible or small, we assume that (3) represents 100 % of the uncertainty -i.e., the grid cell uncertainty is systematic (correlated) in space and random in time -to provide an upper limit on the uncertainty. In addition, we also take into account the sampling uncertainty when temporally  averaging over each grid cell using the standard error; i.e., we take the standard deviation divided by the square root of the number of measurements.
In practice, when averaging the daily values for every grid cell i over months, seasons or years, the uncertainty σ i is obtained by adding in quadrature the daily uncertainties σ 2 N i and dividing by the number of available daily measurements N i : In addition, we add in quadrature σ i and the standard error SE i to obtain the total uncertainty of the temporal average σ i for every grid cell: where SD i is the standard deviation of the daily values in grid cell i. The standard error measures how far the sample could be from the true population mean. Finally, when spatially averaging globally or regionally, under the assumption that the errors are correlated across space, the overall uncertainty is calculated by averaging σ i across the N j grid cells in spatial domain j weighted by the grid cell area fraction with respect to the total area (i.e., grid cell/total area = w i ) with available retrievals:

Results
Our analysis is divided in two main parts. In the first one (Sect. 4.1) we assess the annual and seasonal climatological DOD maps for nine distinct regions. In the second one (Sect. 4.2), emphasis is given on the quantification of DOD averages along with their monthly and interannual variability of the fractional contribution to the AOD, from a global to hemispherical level as well as for specific regional domains.

North Africa, tropical Atlantic Ocean and Mediterranean
According to the long-term average map ( Fig. 1), the maximum DODs (up to 1.2) are recorded in the Bodélé depression, which is considered the most active individual dust source of the planet (Washington et al., 2003;Koren et al., 2006;Ginoux et al., 2012). Over the area, the prevailing strong winds are intensified further between the Tibesti mountains and the Ennedi plateau (Washington et al., 2009), forming a low-level jet (Washington and Todd, 2005). This dominant wind pattern, affected by the local topography (Washington et al., 2009), acts as the driving force mobilizing mineral particles from arid and erodible soils of the region (Tegen et al., 2006). Under these favorable conditions, dust aerosols are easily uplifted and accumulated in the atmosphere, thus causing the very high DODs (>0.5) observed in the broader area (Chad, Niger). Throughout the year, the high DOD levels are quite persistent, exhibiting, however, a seasonal variation with more intense loads recorded during DJF (Fig. S1i) and MAM (Fig. S1ii) following the annual cycle of source activation (Washington et al., 2009). The second hotspot in North Africa is situated between the northern parts of Nigeria and the southern parts of Niger with annual DODs reaching up to 0.7 ( Fig. 1) while on a seasonal basis varying from 0.4 (SON; Fig. S1iv) to 0.8 (JJA; Fig. S1iii). MIDAS DODs match well with those presented by Rajot et al. (2008), who relied on ground-based sun-photometric measurements of AOD obtained at the Banizoumbou AERONET site. Very high DODs are also evident along the coasts of the Gulf of Guinea, which may be unrealistic considering that dust aerosols are mainly transported there and are mixed with anthropogenic and biomass burning (Knippertz et al., 2015). Along this area of high DODs, MERRA-2 also overestimates the dust fraction compared to LIVAS , thus resulting in higher intensities according to the applied methodology (Sect. 2). Moreover, the temporal availability of DODs in the region is very limited (<10 %; Fig. 8c in Gkikas et al., 2021); the DOD uncertainty is large; and AOD outliers, either realistic or cloud contaminated, can yield exceptional high DODs in this complex environment where aerosol and clouds are spatially correlated (Andrew Sayer, personal communication, 2020). This abrupt reduction of DOD levels, from inland to the nearby maritime environment, reveals an artifact of the MIDAS dataset mainly introduced by the raw MODIS AOD retrievals, which are obtained by retrieval algorithms built on different assumptions/considerations depending on the underlying surface type. Across the Sahara Desert, there is a distinct longitudinal contrast with more intense dust loads in western North Africa than in eastern North Africa (Fig. 1). In the former sector, the DODs range mainly from 0.3 to 0.6, while over the eastern parts of the Sahara the corresponding limits are bounded between 0.1 and 0.3 without revealing significant intra-annual variation. During MAM (Fig. S1ii), along the southern Sahel, the activation of dust sources results in DODs which locally can exceed 0.8, while during boreal summer ( Fig. S1iii) a vast area of the Western Sahara is under the impact of heavy dust loadings (DOD >0.5). According to Ginoux et al. (2012), in the former region, dust is mainly produced by agricultural activities (cultivation, overgrazing) disturbing soils in which alluvial sediments have been accumulated. Northwards, dust has natural origin, and the accumulation of mineral particles is favored by the development of the Saharan heat low (SHL) affecting also the prevailing airflow (harmattan winds) as well as the West African monsoon (WAM) . Under these meteorological conditions, several dynamic processes, from microscale to mesoscale, are taking place, triggering dust emission (Knippertz and Todd, 2012) from highly active sources (Schepanski et al., 2007).
Under the impact of the trade winds, Saharan dust can travel across the tropical Atlantic Ocean, reaching the Caribbean Sea, the southern United States and northeastern South America (Prospero, 1999;Prospero et al., 2014). The signal of this long-range transport is evident on the annual climatological pattern ( Fig. 1) with DODs up to 0.6 (off the Western Saharan coasts) fading down to 0.1 at the maximum distance. Within the course of the year, the Saharan dust plume varies in terms of intensity, range and latitudinal position, as it is depicted in Fig. S1. During boreal summer (Fig. S1iii), the corridor of the transatlantic dust transport is bounded between 10 and 20 • N latitudes, whereas both the intensity (DODs up to 0.6) and the range are maximized. During boreal winter (Fig. S1i), the dust zone migrates southwards (between the Equator and 10 • N), while maximum (up to 0.6) and considerable (0.1-0.2) DODs are observed over the Gulf of Guinea and mid-Atlantic (45 • W), respectively. Between the transition seasons ( Fig. S1ii, iv), dust loads are stronger in MAM (∼ 0.45), mainly residing within 5 and 20 • N latitudes, in contrast to SON (∼ 0.3) when they are shifted northwards (10 and 25 • N). According to the existing literature, several factors modulate the westwards propagation of dust plumes, originating in the Western Sahara and the Bodélé Depression, over the tropical Atlantic. For instance, the south-north displacement of the Saharan plumes is driven by the location of the Intertropical Convergence Zone (ITCZ) and the disturbances of the African easterly jet (Knippertz and Todd, 2012;Doherty et al., 2012). Teleconnection patterns, such as the El Niño-Southern Oscillation (ENSO; Prospero and Lamb, 2003), the North Atlantic Oscillation (NAO; Ginoux et al., 2004) and the North African Dipole Index (NAFDI; Rodríguez et al., 2015), have been also studied in order to interpret the decadal variations of dust concentrations over the Atlantic. Likewise, the vegetation coverage across the Sahel and the wind speeds, determined by the prevailing atmospheric circulation, over the Sahara play a key role on the amount of the emitted dust particles.
Due to the vicinity of the largest deserts of the planet, the Mediterranean is affected by dust outbreaks throughout the year (Gkikas et al., , 2016Marinou et al., 2017). Mineral particles originating primarily from north African and secondarily from Middle Eastern deserts are transported towards the Mediterranean mainly under the prevalence of cyclonic systems . The intensity of dust loads decreases for increasing latitudes, forming a distinct south-north gradient with DODs up to 0.20 between the gulfs of Gabes (Tunisia) and Sidra (Libya), according to the annual pattern ( Fig. 1). Among seasons (Fig. S1), DODs vary on the locations where the maximum levels are recorded as well as on their magnitude, attributed to the position of the prevailing synoptic systems . The central and eastern Mediterranean sectors are affected by dust loads mainly in spring (DODs up to 0.3; Fig. S1ii) and winter (DODs up to 0.12; Fig. S1i). In summer (Fig. S1iii), dust activity is more pronounced in the western parts with optical depths up to 0.18 (Alboran Sea), while thanks to the fine-resolution product, hotspots of similar DODs can be identified in the southern parts (Andalucia) of Spain. In SON (Fig. S1iv), dust loads are found in the central Mediterranean with DODs lower than 0.12 off the Tunisian and Libyan coasts.

Middle East
In the Middle East, there is a zone of moderate-to-high DODs (locally up to 0.8) extending from the historical region of Mesopotamia to the southern parts of Saudi Arabia, where one of the largest sand deserts of the world (Rub' al Khali) (Hamidi et al., 2013) is situated (Fig. 2). Based on Ginoux et al. (2012), the origin of mineral particles between Tigris and Euphrates as well as across the Rub' al Khali Desert is mainly natural, while in the intermediate part (Ad-Dahna Desert) dust accumulation is attributed to the mixing of anthropogenic and hydrological sources. Slightly higher maximum DODs (up to 0.7; Fig. 2) are recorded in Oman and particularly between Dhofar and Al Wusta, in contrast to previous studies (Pease et al., 1998) which have identified the Sharqiya Sands (formerly known as the Wahiba Sands) area as a major dust source or the coastal areas of Yemen (Ginoux et al., 2012). On a seasonal basis, the intensity of mineral loads exhibits a strong variability with minimum DODs (up to 0.4) during DJF (Fig. S2i) and SON (Fig. S2iv) and maximum (up to 1) during the dry period of the year (Fig. S2ii, iii), being in agreement with the results presented in Yu et al. (2013). More specifically, across the Arabian Peninsula, the increase in DOD levels is getting evident in boreal spring, and it is further intensified during summer months. Dust storms emanating in Iraq and the eastern parts of Saudi Arabia favor dust transport towards the Persian Gulf (Gianakopoulou and Toumi, 2012) and account for the considerably high DOD levels (>0.6) found there. Due to convergence of the northern-northwesterly Shamal winds (Yu et al., 2016) and the airflow from the subtropical anticyclone, in JJA, mineral particles are traveling at even longer distances towards the northern Arabian Sea (Ramaswamy et al., 2017), as indicated by the intense dust loads (DODs up to 0.5; Fig. S2iii) contributing about half of the AOD (Jin et al., 2018). Likewise, during boreal summer, short-range dust transport takes place off the coasts of Oman and Yemen (Gulf of Aden). Among seas in the vicinity of the Arabian Peninsula, the most intense dust loads are observed in the Red Sea, forming a clear latitudinal gradient on annual (Fig. 2) and summer (Fig. S2iii) geographical DOD patterns, as it has been noted also in Brindley et al. (2015) and Banks et al. (2017). Due to its location, the southern sector of the Red Sea receives dust aerosols either originating from the Republic of the Sudan or from the Arabian Peninsula, depending on the zonal airflow (Banks et al., 2017). Dusty air masses traveling westwards are uplifted when they are crossing the mountain range in the southwestern Arabian Peninsula, and for this reason dust loads over the southern basin are suspended above 2 km (Banks et al., 2017). On the contrary, low-elevated dust layers are recorded when winds blow from west, triggering dust emission from the Tokar Gap (Sudanese coasts), and subsequently dust outflows into the southern Red Sea (Banks et al., 2017).

Central and southwest Asia
Northwards and eastwards of the Caspian Sea, various deserts are situated in the central segments of the Asian continent. Most of Turkmenistan is occupied by the Karakum Desert, while the Kyzylkum Desert is located in Uzbekistan. Other arid regions stretch between the Caspian and Aral seas (Ustyurt Plateau), in the eastern and southern flanks of the Aral Sea and in the lowlands of western Kazakhstan and southeastern Russia (Ryn Desert) (Elguindi et al., 2016). Based on our seasonal spatial patterns (Fig. S3), the major dust activity is recorded in the Ustyurt Plateau (Li and Sokolik, 2018) and in the large lagoon embayment of Garabogazköl (Shen et al., 2016), a gulf of Turkmenistan dried into a salt-covered playa (Gill, 1996), with minimum (in DJF and SON) and maximum (in MAM and JJA) DODs equal to ∼ 0.2 and ∼ 0.4, respectively. In the rest of areas, the corresponding upper limits can reach up to 0.8-0.9, during boreal summer, in localized spots (Chimboy Lake, Sarygamysh Lake) across the Karakum and Kyzylkum deserts. In the same season, moderate dust loadings (DOD up to 0.25) are encountered in the southern Caspian Sea (Elguindi et al., 2016) as the result of transported mineral particles mainly coming from the sandy deserts of Turkmenistan (Xi and Sokolik, 2015), under the impact of east-ern/southeastern winds (Shen et al., 2016). Since the 1960s, the anthropogenic intervention (agricultural activities, overirrigation) caused the retreat of the Aral Sea and the formation of the Aralkum Desert (Saiko and Zonn, 2000;Micklin, 2007), from which large amounts of aeolian dust are emitted and travel distances of hundreds of kilometers (Indoitu et al., 2015). According to the annual climatological map (Fig. 3), extremely high DODs (>1) are found in the southeastern parts of the Aralkum Desert (Fig. 3) which are also persistent among the seasons (Fig. S3). Nevertheless, these are not trustworthy as it has been thoroughly discussed in Gkikas et al. (2021) In the Sistan Basin, extending between Iran, Pakistan and Afghanistan, the long-term average JJA DODs can reach up to 1.1 (Fig. S3iii) in the Margo Desert (Afghanistan), due to the frequent occurrence of dust storms (Middleton, 1986) triggered by the northerly Levar winds blowing from June to September (Alizadeh- Choobari et al., 2014). These maximum DOD levels are substantially higher than the annual mean (0.8; Fig. 3) as well as against the corresponding averages for the other seasons. Thanks to the high-resolution MIDAS DOD, we identify the borders of other active arid regions, surrounded by mountain ranges, such as the Rīgestān (Afghanistan), the Balochistan (Pakistan), the Dasht-e Kavir (Iran), the Dasht-e Lut (Iran) and the Jazmurian drainage basin (Iran). In the aforementioned topographic lows, the magnitude of the dust loads is significantly lower than those observed in the Margo Desert and can be as large as 0.6 (Balochistan) during hot-dry months (Fig. S3iii). The presence of absorbing mineral particles, over the area and in the northernmost part of the Arabian Sea, is also confirmed by the high AI values, especially in June-July, discussed by Rashki et al. (2015), who relied on long-term records obtained by the Ozone Monitoring Instrument (OMI) and Total Ozone Mapping Spectrometer (TOMS) spaceborne sensors.

Indian subcontinent
In the Indian subcontinent, the maximum annual DODs (∼ 0.5; Fig. 4) are observed along the Indus River basin, in the western side of the Thar Desert, whereas a branch of gradually decreasing DODs, along the Indo-Gangetic Plain towards eastwards directions, is also evident. Ginoux et al. (2012) stated that much of dust activity in the Indus River basin is attributed to the suspension of soil particles originating primarily from agricultural land use and to a lesser extent from the desiccation of ephemeral water bodies. The strong presence of absorbing coarse particles over the area is further supported by the coexistence of considerably high aerosol index (AI) values (Alam et al., 2011). As indicated by the seasonal patterns (Fig. S4), the processes regulating the suspended dust loads are highly variable during the year, causing a remarkable temporal variability of DOD, which is low (<0.3) in DJF and SON, moderate in MAM (<0.5), and maximum in JJA (<0.8). Similar seasonal vari-ability is evident in the Thar Desert, in agreement with the findings of Proestakis et al. (2018) and Dey and Di Girolamo (2010), who used vertically resolved (CALIOP) and multi-angle (MISR) satellite retrievals, respectively. Nevertheless, our climatological DODs are higher with respect to the CALIOP corresponding values and the MISR nonspherical AODs, particularly when dust activity over the area is pronounced. During the pre-monsoon season, westerly to northwesterly winds are blowing over the Thar Desert, mobilizing dust particles which subsequently are advected towards the Indo-Gangetic basin (Dey et al., 2004;Srivastava et al., 2011). According to our results, between the Haryana state and the eastern parts of the plain, DODs fade down from ∼ 0.6-0.7 to ∼ 0.1-0.2, forming a NW-SE gradient (Fig.  S4ii, iii). Such high DODs are attributed to the eastwards propagation of intense dust storms having a strong signature on the optical, microphysical and radiative properties derived by AERONET stations operating in the region Eck et al., 2010).

East Asia and North Pacific Ocean
Northwards of the Tibetan Plateau, the Tarim Basin (northwest China) is situated, in which one of the largest natural dust sources, the Taklamakan Desert, resides. This elevated desert area (average elevation 1.1 km) is surrounded by the Pamir Plateau (average elevation 5.5 km) to its west, by the Kunlun Shan range (average elevation 5.5 km) in its southern flanks and by the Tian Shan range (average elevation 4.8 km) along its northern boundaries, while the ground elevation is low only in its eastern margin (Ge et al., 2014). DODs are maximized in spring (Fig. S5ii), yielding values up to 1 along the foothills of the Tian Shan and Kunlun Shan ranges, attributed to the role of the topography on winds strengthening (Ge et al., 2014). Similar values are recorded in JJA (Fig. S5iii), but the geographical distribution reveals that the highest DODs are less widespread in contrast to spring. Throughout the year, the weaker dust loads are recorded during winter and autumn. Our results are consistent with relevant studies that rely on active and passive satellite retrievals either of pure dust load  or AOD Sogacheva et al., 2018).
A common feature in the seasonal DOD patterns is the reduction of dust loads' intensity towards the interior parts of the Taklamakan Desert, as it has been also documented by Ge et al. (2014), who utilized MISR retrievals. The highresolution of the MIDAS DOD dataset provides in detail the spatial information of these geographical patterns. During spring, similar high DODs to those found over the Taklamakan Desert are recorded in the Qaidam Basin (northeast side of the Tibetan Plateau), surrounded by the Altun, Kunlun and Qilian mountain ranges, attributed to strong downslope winds causing the erosion of soil particles (Rohrmann et al., 2013) and their entrainment into the atmosphere. The intensity of dust loads over the Gobi Desert (north China-south Mongolia) hardly exceeds 0.3 on an annual basis (Fig. 5), while it can reach up to 0.4 during spring (Fig. S5iii). The remarkable deviations in dust abundance between Taklamakan and Gobi during springtime are interpreted by variations in soil characteristics. More specifically, Taklamakan is composed mainly by fine sand particles in contrast to the rocky soils of the Gobi Desert (Sun et al., 2013). Due to these differences in soil textures, dust particles from the former desert region can be emitted even with low wind speeds, while they are uplifted at higher elevations in the troposphere, as it has been shown with MISR stereo observations (Yu et al., 2019) and CALIOP lidar profiles . The injection of Taklamakan dust particles at higher altitudes increases their residence time inducing also their entrainment into the upper-level westerly airflow, around at 4 a.m.s.l., both contributing to the higher potential for long-range transport (Yu et al., 2019), in contrast to Gobi dust, towards continental East Asia and the northern Pacific Ocean. Under the impact of cold fronts propagating eastwards (Eguchi et al., 2009) in spring, air masses carrying mineral particles during the first 2 d of dust transport affect a wide area of China (Yu et al., 2019), from near sources to its eastern parts through the Hexi Corridor and the Loess Plateau (DODs ranging from 0.2 to 0.4; Fig. S5iii). Subsequently, the Asian dust plumes are suspended over the Yellow Sea, the Korean Peninsula and further eastwards, in a latitudinal band bounded between the parallels 30 • N and 45 • , reaching the west coasts of the United States (Yu et al., 2008). Across this belt, where the trans-Pacific dust transport is taking place, the springtime DODs decrease smoothly from 0.15 to 0.05 (Fig. S5ii). In summer (Fig. S5iii), DODs up to 0.05 are observed between 40 and 60 • N, indicating a northwards displacement of the Asian dust layers (mainly originating from the Gobi Desert) due to the weakening and northwards shift of the polar jet streams (Yu et al., 2019).

North America
Across North America, major dust activity is detected in the southwestern United States and in northwest Mexico, with annual and seasonal DODs hardly exceeding 0.15, as illustrated in Figs. 6 and S6, respectively. These weak dust load intensities are mainly recorded in the Sonoran and the Mojave deserts, while lower values are found in the Chihuahuan Desert in which isolated spots (e.g., White Sands desert) become visible thanks to the high-resolution of the MIDAS DOD dataset. Low-to-moderate DODs are evident in the Great Plains with local maxima (exceeding 0.2 in spring; Fig. S6ii) in the Great Salt Lake Desert and in the surrounding area as well as in the Baja Californian Desert (Mexico; DODs up to 0.14), residing in the western side of the Gulf of California. Our annual spatial distribution of DOD (Fig. 6) is highly consistent with those of frequency of observation (FoO) of DOD (Ginoux et al., 2012;Baddock et al., 2016) and AI given by Prospero et al. (2002).
Moreover, the increase in dust loads' concentration in MAM (Fig. S6ii) has also been documented by Hand et al. (2017) and Tong et al. (2017), both relying on aerosol observations acquired at numerous stations of the Interagency Monitoring of Protected Visual Environments (IMPROVE) network. During springtime, dust emission over the broader area is associated with the transmit of Pacific cold fronts inducing dust-entraining winds as the result of pressure gradient enhancement (Rivera Rivera et al., 2009). The geomorphological soil characteristics are determinant for dust emission with the most prominent natural sources being ephemeral and dry lakes , while anthropogenic dust aerosols are mainly emitted in the Great Plains and in the eastern side of the Gulf of California (Ginoux et al., 2012).

Australia
Earlier studies based on unconstrained numerical simulations (Tanaka and Chiba, 2006;Wagener et al., 2008) have shown that among the desert areas of the Southern Hemisphere, the largest contribution of dust particles arises from Australia. However, a more recent assessment (Kok et al., 2021b) in which dust models have been constrained by observations revealed that the emitted dust amounts from South America are slightly higher than those of Australia. Due to the fairly bright landmasses and the predominance of weak aerosol loadings, there is minimal contrast between surface and atmosphere, leading to systematic algorithm uncertainties, which can explain the slightly lower land DODs compared to those recorded in the surrounding oceanic regions (Figs. 7 and S7). Nevertheless, in the sources as well as in areas affected by dust plumes, the atmospheric signal becomes evident. In particular, the highest dust emissions are encountered in the Lake Eyre Basin (LEB; Prospero et al., 2002) composed by ephemeral lakes, alluvial channels, gibber (stone-covered plains), aeolian sand deposits and bedrock (Bullard et al., 2008). Based on the annual climatological pattern (Fig. 7), DODs can locally exceed 0.2 (in the southern parts) but in general vary between 0.06 and 0.12. From a seasonal perspective (Fig. S7), the highest DODs (mainly up to 0.18 in the Warburton River estuary; few exceedances above 0.4 are found in local spots) are recorded during austral summer (DJF; Fig. S7i) and spring (SON; Fig. S7iv). Similar seasonal variation in ground-based sun-photometric observations at nearby sites (Birdsville, Tinga Tingana), with slightly lower AODs, has been reported by Mitchell et al. (2017). Southwards of the LEB, three spots of notable DODs (up to 0.2 in SON; Fig. S7iv) are identified in the lakes Gairdner, Torrens and Frome, while northeastwards (Lake Yamma Yamma) and northwards (Simpson Desert) from the basin the suspended dust loads exhibit optical depths as large as 0.12 during the driest months of the year. Similar maximum DODs are recorded in the Northern Territory and in the western side of the Great Dividing Range (Queensland), and in contrast to Ginoux et al. (2012) these levels appear in DJF instead of SON. In the southwestern coastal parts of the Australian landmass as well as in Riverina (southeast), during austral spring (Fig. S7iv) very low DODs are evidently associated with anthropogenic dust originating from agricultural activities (Ginoux et al., 2012). Finally, during the same season, weak signals (DODs up to 0.05) of dust transport are revealed over the Tasman and Timor seas, attributed to the eastward movement of cyclonic frontal systems causing the entrainment of mineral particles in air masses that can travel long distances (Knight et al., 1995;Alizadeh Choobari et al., 2012).

South Africa
Dust activity in South Africa is mainly related with shortrange and short-lived plumes (Vickery et al., 2013) that are suspended at low tropospheric altitudes (below 600 hPa) due to the predominance of anticyclonic circulations inhibiting the vertical extension of dust layers (Piketh et al., 1999). Mineral aerosol loadings are mainly originating from the ephemeral lake basins of the Etosha Pan (Namibia) and Makgadikgadi Pan (Botswana) and the Namib Desert (Bryant et al., 2007;Vickery et al., 2013). In the aforementioned source areas, the maximum annual (Fig. 8) and seasonal (Fig. S8) DODs are equal to 0.1 and 0.16, respectively. Throughout the year, the increase in DODs in Etosha and Makgadikgadi pans is evident primarily in DJF (Fig. S8i) and secondarily in SON (Fig. S8iv). Our results are consistent with those provided by Ginoux et al. (2012) and Bryant et al. (2007) for the former region (including also the Kalahari Desert in which very weak dust loads are recorded), in contrast to the latter one and opposite to the findings of Vickery et al. (2013) for both sources. In these arid areas dust emission is linked with lakes' inundation, characterized by strong intra-annual variability, playing an important role when different time periods are considered. However, the moderate performance of the MERRA-2 dust portion with respect to LIVAS in South Africa as well as in most desert areas of the Southern Hemisphere must be also taken into account . Along the Namibian coastline, the deviations of DOD between the high-and low-dust seasons are small, indicating that dust activity remains relatively constant within the course of the year (Ginoux et al., 2012). Soil particles from salt pans and dry river beds of the Namib Desert are emitted from aeolian processes related to berg winds (katabatic winds) blowing in the escarpment, from the Central Plateau down to the coasts (Eckardt and Kuring, 2005). Dust outflow towards the southern Atlantic Ocean, with a SE-NW orientation; it is shown between 18 and 9 • S during austral winter (DODs up to 0.08; Fig. S8iii), becoming more evident in SON (Fig. S8iv), in agreement with the geographical distributions provided by Voss and Evan (2020). Such transport is favored by the propagation of barotropic low-level easterly waves formed between continental highpressure systems and the semi-permanent South Atlantic an-  ticyclone (Tyson et al., 1996). Finally, weak signals of DODs are recorded in the croplands north of Cape Town, with annual and DJF DODs not exceeding 0.1.

South America
In South America, the most intense dust loads are encountered in the Patagonia Desert, where the most active dust sources are situated in the river basins of the Rio Negro and Chubut provinces and in its southern end. Among these areas, higher DODs (up to 0.16 in DJF; Fig. S9i) are found along the Rio Negro, attributed to anthropogenic dust originating from overgrazing, irrigation and oil prospecting (Mc-Connell, et al., 2007;Mazzonia and Vazquez, 2009). In southern latitudes, mineral particles originate from glacier washout plains. Under favorable meteorological conditions, aeolian dust from Patagonia travels either towards the southern Atlantic Ocean, contributing to iron concentrations and marine biological productivity in the surface waters (Johnson et al., 2011), or towards the Antarctica peninsula (Gassó et al., 2010), as it has been found in ice core samples (Basile et al., 1997). Both transport pathways are not visible in our climatological patterns (Figs. 9 and S9) since dust outbreaks are not so strong (Foth et al., 2019), while the extended cloud coverage over the region results in large observational gaps of the spaceborne retrievals (Gassó and Torres, 2019). Along the western side of Andes, dust emission arises from natural sources located in the Sechura (Peru), Nazca (Peru) and Atacama (Chile) deserts (Ginoux et al., 2012). In the aforementioned regions, the annual DODs (Fig. 9) can reach up to 0.1, 0.08 and 0.06, respectively, while the intra-annual variability is characterized as weak (Fig. S9). During MAM (Fig. S9ii), DODs up to 0.16 appear in Guyana, Suriname and French Guiana as well as over their offshore areas, while similar intensities are evident in the northern parts of the Amazon rainforest (around the Equator and bounded between 65 and 60 • W). The presence of coarse mineral particles (Moran-Zuloaga et al., 2018) over these distant areas from deserts is attributed to the long-range dust transport from North Africa across the Atlantic Ocean (Yu et al., 2015), under the impact of the trade winds taking place northwards of the convective precipitation zone formed around the ITCZ. Finally, the latitudinal zone of weak DODs in the western parts of Brazil, fading down abruptly eastwards of ∼ 58 • W, indicates an artifact of the MIDAS product that becomes more evident in SON (Fig. S9iv). This peculiar pattern is induced by the MERRA-2 dust fraction (results not shown here), which is used for the derivation of MIDAS DOD from the MODIS AOD. An additional deficiency is the relatively large DODs over an area where biomass burning particles, emitted at enormous amounts by extended wildfires, clearly dominate over other aerosol species. Under these conditions, the nondust AODs are very high as well as their relevant uncertainties (Eqs. 5-7 in Gkikas et al. 2021) while the reliability of the MERRA-2 dust fraction decreases there (see Fig. 2 in Gkikas et al., 2021).

DOD averages and variability at global, hemispherical and regional scales
In this section, we discuss the average AOD and DOD along with their monthly and interannual variability at global, hemispherical and regional scales. The left column of Fig. 10 shows the interannual time series of AOD (black curve) and DOD (red curve) averaged over the whole globe (upper panel; GLB), the Northern Hemisphere (middle panel; NHE) and the Southern Hemisphere (bottom panel; SHE). The right column of Fig. 10 depicts the monthly seasonal cycle of AOD and DOD along with the DOD-to-AOD ratio (blue curve), while the shaded areas correspond to the total uncertainty (see Sect. 3.2 in Gkikas et al., 2021, and Sect. 3 in the current study). The significant role of dust particles in the global aerosol budget becomes evident by visually inspecting the AOD and DOD interannual time series (Fig. 10a.i). The monthly contribution of suspended dust to the total AOD varies from 14 % to 39 %, with minimum values mainly in DJF and maximum values in MAM or JJA, depending on the year. Monthly DODs range from 0.016 ± 0.013 (December 2005)  (Table 1). The global DOD mean, computed here from the fine-resolution data, is almost identical with those obtained by the coarse-spatialresolution MERRA-2 and MIDAS DODs and slightly higher than those calculated based on LIVAS-CALIOP (0.029) (see Table 1 in Gkikas et al., 2021; it is noted that the three datasets were collocated based on the spatial resolution and the temporal availability of the LIVAS dataset). Likewise, our global average and uncertainty computed over the period 2004-2008 (0.033 ± 0.004) is close to the one obtained in Ridley et al. (2016) (0.030 ± 0.005), despite the different methods applied for the derivation of DOD and its uncertainty. Our global DOD long-term average is very close to the CALIOP-derived value (0.029) and about half of the MODIS-derived one (0.067) reported by Song et al. (2021). Table 1. Annual and seasonal DOD averages, representative for the period 2003-2017, along with the associated uncertainty. The first three rows refer to the whole globe (GLB), the global land (GLB-land) and global ocean (GLB-ocean). In the fourth and fifth line the results for Northern Hemisphere (NHE) and Southern Hemisphere (SHE) DODs are given, whereas in the remaining 17 entries the corresponding results for selected subregions (denoted with colored rectangles in Fig. 7 in Gkikas et al., 2021)  Our continental (0.070 ± 0.005) and oceanic (0.019 ± 0.002) mean DODs (see Table 1) are substantially lower than those obtained in Voss and Evan (2020) (land: 0.1; ocean: 0.03). This difference may be attributed to the different averaging approaches, which can have an important impact on the calculations as it is shown in Levy et al. (2009) (see their Fig. 5). Based on our method, we are giving the same weight at each grid cell (regardless of the amount of available data in that grid cell throughout the study period) when we are calculating the domain (from regional to global) average. Therefore, we are avoiding an overestimation of the spatial average since MIDAS data availability is larger over/nearby deserts (see Fig. 8c in Gkikas et al., 2021) where the higher DODs are observed. To be more specific, when we are calculating the global long-term DOD average based on the second branch (i.e., straight, the standard approach for the calculation of the average value by considering all the available values in space and time) in Levy et al. (2009), we obtain a climatological value equal to 0.047. Such different approaches for the calculation of the long-term DOD averages might interpret the deviations found between this study and Song et al. (2021). Finally, the computed global mean MIDAS DOD is somewhat higher than those simulated by most AeroCom Phase I models (Huneeus et al., 2011), being about 40 % higher than the median (0.023); nevertheless, it must be taken into account that most models consider the diurnal variation of DOD in contrast to the single measurements taken during MODIS overpass.
As expected, the interannual GLB DOD time series is driven by the variability in the NHE DOD (Fig. 10a.ii) since the most widespread and intense dust sources are located in the Northern Hemisphere. This is justified by their high temporal co-variation, while a positive NHE-GLB offset is constantly observed, being lower during boreal winter and autumn (up to 0.035) and maximum during the high dust seasons (0.058). The fraction of monthly NHE AOD attributed to dust particles ranges from 20 % to 48 %, and the R 2 value between monthly AOD and DOD is equal to 0.94, both indicating a dominant dust contribution. Over the study period (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017), the NHE DOD yields a climatological mean equal to 0.056 ± 0.004 (Table 1) ranging from 0.024 ± 0.015 (December 2005) to 0.121±0.050 (March 2012). In contrast, marine and biomass burning aerosols, rather than dust, regulate AOD in the Southern Hemisphere (Fig. 10a.iii). SHE DODs are estimated to be low (0.008±0.001), with the maximum value (0.016 ± 0.016) recorded in February 2016. The contribution of dust aerosols to the total aerosol load does not exceed 17 % throughout the study period (Fig. 10a.iii), and on average it is equal to 8.2 % ± 1.1 %, which is in very good agreement with the findings by Kok et al. (2021b).
A better view of the seasonal cycles of AOD, DOD and the DOD-to-AOD ratio can be obtained by investigating their climatological patterns, representative for the period of interest (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017). On a global scale (Fig. 10b.i), DODs peak between March and June (∼ 0.045) and then decline until November (0.018) before rising during boreal winter. Despite the monthly shifts between maximum AOD and DOD averages, the seasonal cycles of the total aerosol and dust burdens are similar to a large extent, whereas the contribution of mineral particles to the total extinction ranges from 16 % (November) to 33 % (March-June). The MIDAS global DOD-to-AOD ratio (∼ 23 %) is close to the values reported by Gelaro et al. (2017) and Kinne et al. (2006), ∼ 22 % and ∼ 26 %, respectively, but higher than most of the modelderived estimations (12 %-28 %) from the AeroCom Phase III (Gliß et al., 2021). These discrepancies, excluding the aerosol parametrizations, may be partly due to the different sampling between single-overpass satellite observations and reanalyses (Gelaro et al., 2017) or models (Kinne et al., 2006) where the diurnal aerosol variability Yu et al., 2021) is included. In the NHE (Fig. 10b.ii), the mean seasonal trend of DODs remains relatively unchanged when compared with GLB; however, the hemispheric means (0.030-0.088) and the dust fraction (24 %-41 %) are higher. On the contrary, the weak signal of aeolian dust in SHE (Fig. 10b.iii) interprets the very low DODs (0.005-0.011) and their minor impact (6 %-12 %) upon AOD magnitude.
The analysis presented above has also been conducted for each one of the 17 subregions illustrated in Fig. 7 in Gkikas et al. (2021), and the main findings are summarized in this paragraph. Among the regional domains, a persistency of high DODs (>0.3), both at interannual and seasonal scales, is found only in BOD, which yields a long-term average value equal to 0.533 ± 0.009, being almost double that of WSA (0.302±0.006) and TAK (0.246±0.020) as illustrated in Table 1. However, when focus is given to individual months, the maximum DODs over the study period ( Fig. 11a.vi) and on their climatological levels are recorded in the Taklamakan Desert and can be as high as 0.868 (April 2007) and 0.600 (April), respectively. Comparable or even higher DODs than those computed in BOD are also evident for specific months in THA (Fig. 11a.vii), GOG (Fig. 11a.xii) and SSA (Fig. 11a.xv) as well as on the monthly time series (THA; Fig. 11b.vii). Mineral particles' contribution to the total AOD (i.e., blue curves in the seasonal cycle plots) is at least 50 % over dust sources or dust-abundant areas in North Africa, the Middle East and Asia, and it is constantly higher than 70 %, reaching up to 95 %, in BOD ( Fig. 11b.i), WSA (Fig. 11b.viii) and TAK (Fig. 11b.vi). Over downwind regions, such as EAS (Fig. 11b.ix), GOG (Fig. 11b.xii), MED (Fig. 11b.xiii) and SSA (Fig. 11b.xv), the dust contribution can prevail over the non-dust portion (GOG, MED, SSA) while in EAS it does not exceed 30 %, due to the predominance of anthropogenic aerosols. In the oceanic areas of the tropical Atlantic and North Pacific, where large-scale dust transport is taking place, AOD and DOD co-vary, indicating that the dust activity regulates the temporal variations of aerosols' load, except during summer months in WNP ( Fig. 11a.xvi, b.xvi). Regarding the seasonal cycle of DOD, the maximum values are recorded either during boreal spring (GOB, CAS, NME, SUS, TAK, EAS, ENP, GOG, MED, WNP and SSA) or during boreal summer (THA, WSA, ETA, SME and WTA) or are similar between the two high-dust seasons (BOD).
A final intercomparison of the MIDAS DODs against those derived by Ridley et al. (2016) and Adebiyi et al. (2020), on a seasonal basis over the period 2004-2008, has been performed for 15 regions defined in Kok et al. (2021a) (see their Fig. 2b and Table 2). The obtained results are illustrated in Fig. 12. For the Southern Hemisphere regions (Figs. 12xiii,xiv,xv) as well as for North America (Fig. 12xii), MIDAS DODs are compared versus those from Adebiyi et al. (2020), while for the remaining 11 domains (Fig. 12i-xi) the results from Ridley et al. (2016) have been utilized. As an overview, it is noted that the seasonal cycle among the three databases is commonly reproduced, with a few exceptions (Mali-Niger, Kyzylkum, southern Africa), whereas the DOD uncertainties (represented by the error bars) are comparable. Regarding the magnitudes, MIDAS DODs are mainly somewhat lower than those of Ridley et al. (2016) across the dust belt in contrast to the outflow region of the mid-Atlantic (Fig. 12i). The obtained differences are mainly attributed to the consideration of different models for accounting for the non-dust portion, the different treatment of AODs (bias correction vs. quality filtering), the different versions of MODIS retrievals (C006 vs C061), the consideration of multi-satellite observations instead of relying only on MODIS-Aqua retrievals and to the different spatial scales (coarse vs. fine). In relative terms, the largest deviations are found in the desert areas of the Southern Hemisphere where models struggle to represent adequately the dust sources and the emitted amounts of mineral particles, thus affecting the dust fraction ratio provided by MERRA-2.

Summary and conclusions
The current study presents a scientific exploitation of the MI-DAS dataset , which provides columnar mid-visible (550 nm) dust optical depth (DOD) at fine spatial resolution (0.1 • × 0.1 • ) and over a 15-year period (2003-2017). Taking advantage of the global coverage of the MIDAS DOD product, we analyzed the contribution of dust aerosols to AOD at various spatial and temporal scales. More specifically, we focused on nine regions that account for the majority of the global dust budget, encompassing sources and downwind areas with the main dust transport pathways. Such regions comprise the deserts extending across the dust belt, North America, Australia, South Africa, and South America as well as maritime areas (tropical Atlantic Ocean, Mediterranean, North Pacific Ocean) receiving constantly large amounts of mineral particles from the nearby   Table 2). The error bars represent the estimated uncertainties. From (i) to (xi), the blue bars correspond to the Ridley et al. (2016) results, whereas for the remaining regions MIDAS DODs are compared against the corresponding levels obtained by Adebiyi et al. (2020). deserts. As a further step, the interannual and intra-annual time series of DODs along with their contribution to the total aerosol load (AOD) were investigated at global, hemispherical and regional level.
According to our findings, the global long-term DOD average over the study period (2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016)(2017) is equal to 0.032 ± 0.003, yielding a strong contrast between the contributions from the Northern Hemisphere (0.056±0.004) and Southern Hemisphere (0.008 ± 0.001). Our global estimations are almost identical with those given by Ridley et al. (2016) and the CALIOP-derived estimate of Song et al. (2021), in contrast to the MODIS-based average reported in the latter study. Nevertheless, when the global averages are calculated separately over land (0.070 ± 0.005) and ocean (0.019 ± 0.002), our results differ substantially from those found in Voss and Evan (2020), who reported continental and maritime DODs equal to 0.100 and 0.030, respectively. Such large deviations are attributed to the different applied methodologies and averaging procedures followed. Moreover, we find very good agreement, in terms of DOD magnitude and uncertainty, in the MIDAS seasonal DODs versus those of Ridley et al. (2016) and Adebiyi et al. (2020) for 15 regions defined in Kok et al. (2021a). Considering that the long-term DOD averages can be utilized for constraining global dust in climate models or can be used in several other applications, a detailed analysis is required for enlightening the factors resulting in disagreements among studies. Likewise, our computed global DOD average resides around the middle of the AeroCom (Huneeus et al., 2011) limits, being higher than the median (0.023) and mean (0.028). However, in the modelbased calculations the diurnal variability is taken into account in contrast to the satellite-based estimations relying on single-overpass measurements per day.
Regarding the dust contribution to the total aerosol optical depth, the DOD-to-AOD ratio from 32 % in the Northern Hemisphere drops down to 8 % in the Southern Hemisphere, while at a global scale it is about one-quarter (23 %). The contradiction found between the two hemispheres, both for DOD and dust fraction, is interpreted by the most pronounced dust activity recorded in the Bodélé Depression of the northern Lake Chad Basin (DODs up to ∼ 1.2), across the Sahel (DODs up to 0.8), in western parts of the Sahara Desert (DODs up to 0.6), in the eastern parts of the Arabian Peninsula (DODs up to ∼ 1), along the Indus River basin (DODs up to 0.8) and in the Taklamakan Desert (DODs up to ∼ 1). On the contrary, the weaker emission mechanisms triggering dust mobilization over the spatially limited sources of Patagonia, South Africa and interior arid areas of Australia do not favor the accumulation of mineral particles at large amounts (DODs up to 0.4 at local hotspots), even during high-dust seasons. Except for the Bodélé Depression, where the seasonal variability of the intense dust loads is relatively weak, in the other dust sources of the Northern Hemisphere, DODs exhibit a strong seasonal cycle with maximum levels either during boreal spring or summer and minimum in boreal winter.
Over oceans, the main pathways of long-range dust transport are observed along the tropical Atlantic and the northern Pacific, revealing a remarkable variation, within the course of the year, in terms of intensity, latitudinal position and range. Saharan dust plumes, reaching the Caribbean Sea in summer under the impact of the trade winds, are more abundant with respect to Asian dust, arriving at the western coasts of the United States in spring under the impact of midlatitude cyclones. Due to the convergence of the Shamal winds, blowing over the Arabian Peninsula, and the wind flow from the subtropical anticyclone, dust aerosols originating in the Middle East can reach the western Indian coasts in summer, crossing the Arabian Sea. Dust loads in the southern parts of the Red Sea are maximized during boreal summer when Saharan or Middle East dust is transported, depending on the zonal airflow. The intensity of dust burden in the Mediterranean forms a south-north gradient, whereas a seasonal longitudinal shift of the maximum DODs, off the northern African coasts, is evidently attributed to the prevailing synoptic circulation.
Despite the strong capabilities of the MIDAS dataset, we have also identified some limitations, thoroughly discussed here and in Gkikas et al. (2021), attributed either to inherent weaknesses of the raw MODIS AOD retrievals or to deficiencies of MDF, resulting in not-too-realistic patterns in specific regions (e.g., South America, Gulf of Guinea, Aral Sea) of the planet. Thanks to this detailed analysis, potential users are aware of any issue that may rise when utilizing the MI-DAS DOD product.
Concerning the DOD uncertainties presented here, in the MIDAS dataset, MODIS AOD retrievals, obtained based on different assumptions in the respective algorithms, and MERRA-2 products are mixed. Therefore, the AOD and MDF errors, with combined DOD uncertainty and spatially and temporally averaged, are more likely heterogeneous and quite difficult to quantify. Actually, the evaluation of spaceborne retrievals and numerical outputs can be much more complex, and further work is definitely needed towards optimizing the confidence margins of total (speciated) optical depth levels. Accurately quantifying satellite-based aerosol uncertainties is still an open issue, and it is among our priorities to minimize the impacts of the aforementioned drawbacks and misrepresentations in the future versions of the MIDAS dataset.
As already mentioned, a variety of research studies can rely on the MIDAS dataset. MIDAS has been already used for the investigation of DOD trends (Logothetis et al., 2021), whereas in a follow-up study the mechanisms contributing to the temporal variations of dust burden will be investigated. Likewise, the MIDAS DOD product has been utilized in radiative transfer studies (Fountoulakis et al., 2021;Masoom et al., 2021) focusing on the impacts on solar energy production. Moreover, taking advantage of the fine spatial resolution of the MIDAS dataset and of its extended temporal avail-

3572
A. Gkikas et al.: Quantification of the dust optical depth across spatiotemporal scales with MIDAS ability, the dataset can be used for the identification of dust sources worldwide, similarly to the analysis done in Ginoux et al. (2012). Finally, we have provided a simple yet flexible method (independent from other datasets) to calculate consistent uncertainties across spatiotemporal scales, which will ease the use of the MIDAS dataset in data assimilation applications.
Data availability. The MIDAS dataset has been developed in the framework of the DUST-GLASS project (grant no. 749461; European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie Actions) and is available at https://doi.org/10.5281/zenodo.4244106 (Gkikas et al., 2020).
Author contributions. AG was responsible for the design of the study and the whole analysis, with support from SK and CPGP. EP, VA, EDT, EM, NH and JFK provided feedback on the performed analysis. All authors contributed to editing the paper. AG prepared the replies to the three anonymous reviewers, including feedback from SK, CPGP, JFK and EDT.
Competing interests. At least one of the (co-)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

Disclaimer.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement.
This article is part of the special issue "Dust aerosol measurements, modeling and multidisciplinary effects (AMT/ACP inter-journal SI)". It is not associated with a conference.  Entrepreneurship, andInnovation (NSRF 2014-2020) and cofinanced by Greece and the European Union (under the European Regional Development Fund). NOA members acknowledge support from the Stavros Niarchos Foundation (SNF). The authors acknowledge support by the COST Action InDust (grant no. CA16202), supported by COST (European Cooperation in Science and Technology). The authors would like to thank Andrew Mark Sayer for his valuable and constructive comments. The authors would like also to thank Thanasis Georgiou for developing the ftp server on which the MIDAS dataset is stored.

Acknowledgements. Antonis
Financial support. This research has been supported by the H2020 Marie Skłodowska-Curie Actions (grant no. 749461).
Review statement. This paper was edited by Armin Sorooshian and reviewed by three anonymous referees.