Search

In this study, we developed a novel algorithm based on the collocated Moderate Resolution Imaging Spectroradiometer (MODIS) thermal infrared (TIR) observations and dust vertical profiles from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) to simultaneously retrieve dust aerosol optical depth at 10 µm (DAOD10μm) and the coarse-mode dust effective diameter (Deff) over global oceans. The accuracy of the Deff retrieval is assessed by comparing retrieved Deff with the in-situ measured dust particle size distributions (PSDs) from the AER-D, SAMUM-2 and SALTRACE field campaigns through case studies. The new DAOD10μm retrievals were evaluated first through comparisons with the collocated DAOD10.6μm retrieved from the combined Imaging Infrared Radiometer (IIR) and CALIOP observations from our previous study (Zheng et al. 2022). The pixel-to-pixel comparison of the two retrievals indicates a good agreement (R~0.7) and a significant reduction of (~50 %) retrieval uncertainties largely thanks to the better constraint on dust size. In a climatological comparison, the seasonal and regional (5°×2°) mean DAOD10um retrievals based on our combined MODIS and CALIOP method are in good agreement with the two independent Infrared Atmospheric Sounding Interferometer (IASI) products over three dust transport regions (i.e., North Atlantic (NA; R = 0.9), Indian Ocean (IO; R = 0.8) and North Pacific (NP; R = 0.7)). Using the new retrievals from 2013 to 2017, we performed a climatological analysis of coarse mode dust Deff over global oceans. We found that dust Deff over IO and NP are up to 20 % smaller than that over NA. Over NA in summer, we found a ~50 % reduction of the number of retrievals with Deff > 5 μm from 15° W to 35° W and a stable trend of Deff average at 4.4 μm from 35° W throughout the Caribbean Sea (90° W). Over NP in spring, only ~5 % of retrieved pixels with Deff > 5 μm are found from 150° E to 180°, while the mean Deff remains stable at 4.0 μm throughout eastern NP. To our best knowledge, this study is the first to retrieve both DAOD and coarse-mode dust particle size over global oceans for multiple years. This retrieval dataset provides insightful information for evaluating dust long-wave radiative effects and coarse mode dust particle size in models.

Mineral dust aerosols impact the energy budget of Earth through interactions with radiation, clouds, atmospheric chemistry, the cryosphere and biogeochemistry. In this Review, we summarize these interactions and assess the resulting impacts of dust, and of changes in dust, on global climate and climate change. The total effect of dust interactions on the global energy budget of Earth — the dust effective radiative effect — is −0.2 ± 0.5 W m−2 (90% confidence interval), suggesting that dust net cools the climate. Global dust mass loading has increased 55 ± 30% since pre-industrial times, driven largely by increases in dust from Asia and North Africa, leading to changes in the energy budget of Earth. Indeed, this increase in dust has produced a global mean effective radiative forcing of −0.07 ± 0.18 W m−2, somewhat counteracting greenhouse warming. Current climate models and climate assessments do not represent the historical increase in dust and thus omit the resulting radiative forcing, biasing climate change projections and assessments of climate sensitivity. Climate model simulations of future changes in dust diverge widely and are very uncertain. Further work is thus needed to constrain the radiative effects of dust on climate and to improve the representation of dust in climate models.

The Community Atmosphere Model (CAM6.1), the atmospheric component of the Community Earth System Model (CESM; version 2.1), simulates the life cycle (emission, transport, and deposition) of mineral dust and its interactions with physio-chemical components to quantify the impacts of dust on climate and the Earth system. The accuracy of such quantifications relies on how well dust-related processes are represented in the model. Here we update the parameterizations for the dust module, including those on the dust emission scheme, the aerosol dry deposition scheme, the size distribution of transported dust, and the treatment of dust particle shape. Multiple simulations were undertaken to evaluate the model performance against diverse observations, and to understand how each update alters the modeled dust cycle and the simulated dust direct radiative effect. The model–observation comparisons suggest that substantially improved model representations of the dust cycle are achieved primarily through the new more physically-based dust emission scheme. In comparison, the other modifications induced small changes to the modeled dust cycle and model–observation comparisons, except the size distribution of dust in the coarse mode, which can be even more influential than that of replacing the dust emission scheme. We highlight which changes introduced here are important for which regions, shedding light on further dust model developments required for more accurately estimating interactions between dust and climate.

Desert dust accounts for most of the atmosphere’s aerosol burden by mass and produces numerous important impacts on the Earth system. However, current global climate models (GCMs) and land surface models (LSMs) struggle to accurately represent key dust emission processes, in part because of inadequate representations of soil particle sizes that affect the dust emission threshold, surface roughness elements that absorb wind momentum, and boundary-layer characteristics that control wind fluctuations. Furthermore, because dust emission is driven by small-scale (~1 km or smaller) processes, simulating the global cycle of desert dust in GCMs with coarse horizontal resolutions (~100 km) presents a fundamental challenge. This representation problem is exacerbated by dust emission fluxes scaling nonlinearly with wind speed above a threshold wind speed that is sensitive to land surface characteristics. Here, we address these fundamental problems underlying the simulation of dust emissions in GCMs and LSMs by developing improved descriptions of (1) the effect of soil texture on the dust emission threshold, (2) the effects of non-erodible roughness elements (both rocks and green vegetation) on the surface wind stress, and (3) the effects of boundary-layer turbulence on driving intermittent dust emissions. We then use the resulting revised dust emission parameterization to simulate global dust emissions in a standalone model forced by reanalysis meteorology and land surface fields. We further propose (4) a simple methodology to scale up high-resolution gridded dust emissions to the coarse resolution of GCMs. The resulting dust emission simulation shows substantially improved agreement against regional dust emissions observationally constrained by inverse modeling. We thus find that our revised dust emission parameterization can substantially improve dust emission simulations in GCMs and LSMs.

Mineral dust is the most abundant aerosol species by mass in the atmosphere, and it impacts global climate, biogeochemistry, and human health. Understanding these varied impacts on the Earth system requires accurate knowledge of dust abundance, size, and optical properties, and how they vary in space and time. However, current global models show substantial biases against measurements of these dust properties. For instance, recent studies suggest that atmospheric dust is substantially coarser and more aspherical than accounted for in models, leading to persistent biases in modelled impacts of dust on the Earth system. Here, we facilitate more accurate constraints on dust impacts by developing a new dataset: Dust Constraints from joint Observational-Modelling-experiMental analysis (DustCOMM). This dataset combines an ensemble of global model simulations with observational and experimental constraints on dust size distribution and shape to obtain more accurate constraints on three-dimensional (3-D) atmospheric dust properties than is possible from global model simulations alone. Specifically, we present annual and seasonal climatologies of the 3-D dust size distribution, 3-D dust mass extinction efficiency at 550 nm, and two-dimensional (2-D) atmospheric dust loading. Comparisons with independent measurements taken over several locations, heights, and seasons show that DustCOMM estimates consistently outperform conventional global model simulations. In particular, DustCOMM achieves a substantial reduction in the bias relative to measured dust size distributions in the 0.5–20 µm diameter range. Furthermore, DustCOMM reproduces measurements of dust mass extinction efficiency to almost within the experimental uncertainties, whereas global models generally overestimate the mass extinction efficiency. DustCOMM thus provides more accurate constraints on 3-D dust properties, and as such can be used to improve global models or serve as an alternative to global model simulations in constraining dust impacts on the Earth system.

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 (), 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 consistent uncertainties across spatiotemporal scales, which will enhance the use of the MIDAS dataset in a variety of future studies. Antonis Gkikas acknowledges support by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the 2nd Call for H.F.R.I. Research Projects to support Post-Doctoral Researchers (ATLANTAS, project number 544), as well as support from the European Union's Horizon 2020 research and innovation program under the Marie Skłodowska-Curie Actions (grant no. 749461; DUST-GLASS). Vassilis Amiridis acknowledges support from the European Research Council (grant no. 725698; D-TECT). Eleni Marinou was funded by a DLR VO-R young investigator group and the Deutscher Akademischer Austauschdienst (grant no. 57370121). Jasper F. Kok acknowledges support from National Science Foundation (NSF) grant 1552519. Carlos Pérez García-Pando acknowledges support from the European Research Council (grant no. 773051; FRAGMENT); the AXA Research Fund; and the Spanish Ministry of Science, Innovation and Universities (grant nos. RYC-2015-18690 and CGL2017-88911-R). The authors acknowledge support from the DustClim project as part of ERA4CS, an ERA-NET project initiated by JPI Climate and funded by FORMAS (SE), DLR (DE), BMWFW (AT), IFD (DK), MINECO (ES), and ANR (FR), with cofunding by the European Union (grant no. 690462). PRACE (Partnership for Advanced Computing in Europe) and RES (Red Española de Supercomputación) are acknowledged for awarding access to the MareNostrum Supercomputer in the Barcelona Supercomputing Center. We acknowledge support of this work by the PANhellenic infrastructure for Atmospheric Composition and climatE chAnge (PANACEA) project (grant no. MIS 5021516), which is implemented under the Horizon 2020 Action of Reinforcement of the Research and Innovation Infrastructure, funded by the Operational Programme Competitiveness, Entrepreneurship, and Innovation (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.

The role of mineral dust aerosol in the global radiative energy budget is often quantified by the dust direct radiative effect (DRE). The dust DRE strongly depends on dust aerosol optical depth (DAOD), therefore, DRE efficiency (DREE = DRE / DAOD) is widely compared across different studies to eliminate differences due to the various dust loads. Nevertheless, DREE is still influenced by the uncertainties associated with dust particle size distribution (PSD) and optical properties. In this study, we derive a global clear-sky size-resolved DREE dataset in both shortwave (SW) and longwave (LW) at top of the atmosphere (TOA) and surface based on satellite observations (i.e., satellite-retrieved dust extinction spatial and vertical distributions). In the DREE dataset, dust geometric diameter from 0.1 to 100 µm is divided into 10 bins and the corresponding monthly mean DREE (with respect to DAOD at 532 nm) for each size bin is derived by using the Rapid Radiative Transfer Model (RRTM). Three sets of state of the art dust refractive indices (RI) and two sets of dust shape models (sphere vs. spheroid) are adopted to investigate the sensitivity of dust DREE to dust absorption and shape. As a result, the size-resolved dust DREE dataset contains globally distributed monthly mean dust DREE at TOA and surface for each of 10 size bins with 5∘ (longitude) ×2∘ (latitude) resolution as well as for each dust RI and shape combination. The size-resolved dust DREE dataset can be used to readily calculate global dust DRE for any DAOD and dust PSD, including the uncertainty in the DRE induced by dust microphysical properties, (e.g., dust PSD, RI and shape). By calculating dust DRE based on DAOD climatology retrieved from different satellite sensors and based on different dust PSD, we find that uncertainty in the spatial pattern of DAOD induces more than 10 % of the uncertainty in SW dust DRE at TOA. The observation-based dust PSD induces around 15–20 % uncertainty in dust DRE at TOA and in the atmosphere. The sensitivity assessments of dust DRE to dust RI and shape further suggest that dust nonsphericity induces a negligible effect on dust DRE estimations, while dust RI turns out to be the most important factor in determining dust DRE, particularly in SW.

Mineral dust particles suspended in the atmosphere span more than three orders of magnitude in diameter, from

Aircraft measurement campaigns have revealed that super coarse dust (diameter >10μm) surprisingly accounts for approximately a quarter of aerosols by mass in the atmosphere. However, most global aerosol models either underestimate or do not include super coarse dust abundance. To address this problem, we use brittle fragmentation theory to develop a parameterization for the emitted dust size distribution that includes emission of super coarse dust. We implement this parameterization in the Community Earth System Model (CESM) and find that it brings the model in good agreement with aircraft measurements of super coarse dust close to dust source regions. However, the CESM still underestimates super coarse dust in dust outflow regions. Thus, we conclude that the model underestimation of super coarse atmospheric dust is in part due to the underestimation of super coarse dust emission and likely in part due to errors in deposition processes.

The presence of coarse mineral dust in the atmosphere has been substantiated in several recent measurement campaigns, which include observations of particles up to and above 100μm in diameter. Yet, atmospheric dust models either do not include particles larger than 20μm or severely underestimate their concentrations. One possibility for the underestimated concentrations is that models do not represent enhancements of particle transport due to subgrid-scale topography. Here, large-eddy simulations are used in combination with Lagrangian particle tracking to assess the impact of gentle two-dimensional topography with 50 and 100m elevation on the vertical transport of coarse dust in neutrally stratified conditions. The presence of topography significantly increases the likelihood that 5 and 20μm particles reach several hundred meters in altitude. Further, topography increases this likelihood by orders of magnitude for larger 60μm particles. Three mechanisms are observed to contribute to the increased vertical transport: a strong upward mean flow region on the uphill side of the topography, ejection of particles downwind of the topography crest, and enhanced vertical dispersion in the wake of the crest. The compounding effects of these mechanisms provide a pathway for coarse dust emitted from the surface to reach elevations where they can be further transported into the free atmosphere by large-scale motions such as convective plumes. While these findings are motivated by mineral dust observations, they are generally applicable to other heavy aerosols such as pollen.

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 MODIS-Aqua retrievals and 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 AOD – with about one 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 observationally based studies. Our work provides a simple, yet flexible method to estimate consistent uncertainties across spatiotemporal scales, which will enhance the use of the MIDAS dataset in future studies.

Natural wind-eroded soils contain a mixture of particle sizes. However, models for aeolian saltation are typically derived for sediment bed surfaces containing only a single particle size. To nonetheless treat natural mixed beds, models for saltation and associated dust aerosol emission have typically simplified aeolian transport either as a series of non-interacting single particle size beds or as a bed containing only the median or mean particle size. Here, we test these common assumptions underpinning aeolian transport models using measurements of size-resolved saltation fluxes at three natural field sites. We find that a wide range of sand size classes experience "equal susceptibility" to saltation at a single common threshold wind shear stress, contrary to the "selective susceptibility" expected for treatment of a mixed bed as multiple single particle size beds. Our observation of equal susceptibility refutes the common simplification of saltation as a series of non-interacting single particle sizes. Sand transport and dust emission models that use this incorrect assumption can be both simplified and improved by instead using a single particle size representative of the mixed bed.

Sand dunes and other active sands generally have a low content of fine grains and, therefore, are not considered to be major dust sources in current climate models. However, recent remote sensing studies have indicated that a surprisingly large fraction of dust storms are generated from regions covered by sand dunes, leading these studies to propose that sand dunes might be globally relevant sources of dust. To help understand dust emissions from sand dunes and other active sands, we present in situ field measurements of dust emission under natural saltation from a coastal sand sheet at Oceano Dunes in California. We find that saltation drives dust emissions from this setting that are on the low end of the range in emissions produced by non-sandy soils for similar wind speed. Laboratory analyses of sand samples suggest that these emissions are produced by aeolian abrasion of feldspars and removal of clay-mineral coatings on sand grain surfaces. We further find that this emitted dust is substantially finer than dust emitted from non-sandy soils, which could enhance its downwind impacts on human health, the hydrological cycle, and climate.

Mineral dust emission is a major process in determining the global dust cycle. At a global scale there is still uncertainty about the absolute and relative contribution of dust from different source areas and landforms. Dust sources are mainly considered as surfaces containing relatively high percentages of fine particles (e.g. Playas). Yet, active (dune) sands have been identified recently, by remote sensing studies, as dust sources in northern Africa, China, and elsewhere. Previous studies on dust emission from active sands suggested that dust can be generated by different aeolian mechanisms that are related to (i) re-emission of settled dust particles, (ii) clay coating removal, and (iii) abrasion of the sand grains. However, there is only limited information of the relative importance of the different mechanisms for producing this dust under natural aeolian (wind) conditions. This study integrates wind tunnel experiments and high resolution laboratory sand analyses to explore aeolian dust emission from active sands with conditions simulating the natural processes of saltation and explores the role of the different dust emission mechanisms. Sand samples from three sites with different characteristics of grain size, dust content, morphology, and mineralogy were used in the experiments. No dust emission was recorded for shear velocities below the saltation threshold. The initial content of dust-sized particles (>63 μm) in the sand sample was found to influence PM10 emission. PM10 concentrations were increased with the initial content of dust-sized particles in the sand. The experiments identify clay coatings removal as the dominant mechanism over time of dust emission in typical active sand dunes ( 40 μm). The dust emission observed in this study indicates that, in addition to the classic dust sources of non-sandy soils, sand bodies should also be taken into consideration in determining global dust emission.

Titan, the largest moon of Saturn, is characterized by gigantic linear dunes and an active dust cycle. Much like on Earth, these and other aeolian processes are caused by the wind-driven mobilization of surface grains, known as saltation. To date, very little is known about the conditions that allow for the occurrence of saltation on Titan. In fact, Titan saltation may be fundamentally different from Earth saltation given the denser atmosphere, the lower gravity, and the cohesion of its surface grains. Here, we draw on experiments, theory, and modeling to progress towards a comprehensive understanding of saltation on Titan. We find that aerodynamic lifting of surface grains requires strong wind speeds due to the high cohesion of the grains. However, saltation may be sustained through granular splash at wind speeds much smaller than those required to initiate grain motion. This suggests that most saltation transport on Titan is intermittent rather than continuous. We account for these insights by proposing a saltation mass flux parameterization specific for Titan conditions that accounts for transport intermittency, and use it to quantify yearly sediment transport with a general circulation model. The results show that Titan's prevailing atmospheric circulation is capable of generating highly intermittent yet significant saltation, yielding yearly transport rates similar to those on the most active dunes of Mars. Furthermore, we find that accounting for surface topography might be critical to answering open questions related to Titan's landscape evolution, including the formation of linear dunes in opposite direction to the prevailing circulation.

Mineral dust aerosols cool and warm the atmosphere by scattering and absorbing both solar (short-wave: SW) and thermal (long-wave: LW) radiation. However, large uncertainties remain in dust radiative effects, largely due to differences in the dust size distribution and optical properties simulated in Earth system models. Here, we improve the simulated dust properties with datasets that leverage measurements of size-resolved dust concentration and asphericity factor (improved simulation) in a coupled global chemical transport model (IMPACT) with a radiative transfer module (RRTMG) (default simulation). The global and annual average of dust aerosol optical depth at 550 nm (DAOD550) from the improved simulation (0.029) falls within the range of a semi-observation-based estimate (0.030 ± 0.005), in contrast to that (0.023) of the default simulation. Improved agreement against semi-observation-based estimate of the radiative effect efficiency was obtained using less absorptive SW and more absorptive LW dust refractive indices. Our sensitivity simulations reveal that the improved simulation leads to a similar net global dust radiative effect at the Top Of Atmosphere (TOA) on a global scale to the default simulation (−0.08 vs. −0.09 W ·m−2) but results in less cooling at the surface (−0.23 vs. −0.88 W ·m−2), because of enhanced LW warming by coarser aspherical dust. Our results thus suggest less atmospheric radiative heating due to aspherical dust with coarser size over the major source regions (0.15 vs. 0.79 W ·m−2 on a global scale).

A key challenge in accurate simulations of desert dust emission is the parameterization of the threshold wind speed above which dust emission occurs. However, the existing parameterizations yield a...

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

We present the dust module in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) version 2.0, a chemical weather prediction system that can be used for regional and global modeling at a range of resolutions. The representations of dust processes in MONARCH were upgraded with a focus on dust emission (emission parameterizations, entrainment thresholds, considerations of soil moisture and surface cover), lower boundary conditions (roughness, potential dust sources), and dust–radiation interactions. MONARCH now allows modeling of global and regional mineral dust cycles using fundamentally different paradigms, ranging from strongly simplified to physics-based parameterizations. We present a detailed description of these updates along with four global benchmark simulations, which use conceptually different dust emission parameterizations, and we evaluate the simulations against observations of dust optical depth. We determine key dust parameters, such as global annual emission/deposition flux, dust loading, dust optical depth, mass-extinction efficiency, single-scattering albedo, and direct radiative effects. For dust-particle diameters up to 20¿µm, the total annual dust emission and deposition fluxes obtained with our four experiments range between about 3500 and 6000¿Tg, which largely depend upon differences in the emitted size distribution. Considering ellipsoidal particle shapes and dust refractive indices that account for size-resolved mineralogy, we estimate the global total (longwave and shortwave) dust direct radiative effect (DRE) at the surface to range between about -0.90 and -0.63¿W¿m-2 and at the top of the atmosphere between -0.20 and -0.28¿W¿m-2. Our evaluation demonstrates that MONARCH is able to reproduce key features of the spatiotemporal variability of the global dust cycle with important and insightful differences between the different configurations. This research has been supported by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement no. 789630. Martina Klose has received funding through the Helmholtz Association’s Initiative and Networking Fund (grant agreement no. VH-NG-1533). Jeronimo Escribano and Matthew L. Dawson have received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreements no. 754433 (Jeronimo Escribano) and no. 747048 (Matthew Dawson). BSC co-authors acknowledge funding from the following: the European Research Council (FRAGMENT (grant no. 773051)); the AXA Research Fund; the Spanish Ministry of Science, Innovation and Universities (grant no. CGL2017- 88911-R); the EU H2020 project FORCES (grant no. 821205); the CMUG-CCI3-TECHPROP contract, an activity carried out under a program of and funded by the European Space Agency (ESA); and the DustClim project, which is part of ERA4CS, an ERA-NET initiated by JPI Climate and funded by FORMAS (SE), DLR (DE), BMWFW (AT), IFD (DK), MINECO (ES), and ANR (FR), with cofunding by the European Union (grant no. 690462). Ron L. Miller is funded by the NASA Modeling, Analysis and Prediction Program (NNG14HH42I). Jasper F. Kok acknowledges support from the National Science Foundation (NSF) under grant nos. 1552519 and 1856389 and the Army Research Office (grant no. W911NF-20- 2-0150). Yue Huang has received funding from the Columbia University Earth Institute Postdoctoral Research Fellowship and from NASA (grant no. 80NSSC19K1346) awarded under the Future Investigators in NASA Earth and Space Science and Technology (FINESST) program.

Desert dust aerosols affect Earth’s global energy balance through interactions with radiation(1,2), clouds(3,4), and ecosystems(5). But the magnitudes of these effects are so uncertain that it remains unclear whether atmospheric dust has a net warming or cooling effect on global climate(1,4,6). Consequently, it is still uncertain whether large changes in atmospheric dust loading over the past century have slowed or accelerated anthropogenic climate change(4,7–9), and the climate impact of possible future alterations in dust loading is similarly disputed(9,10). Here we use an integrative analysis of dust aerosol sizes and abundance to constrain the climatic impact of dust through direct interactions with radiation. Using a combination of observational, experimental, and model data, we find that atmospheric dust is substantially coarser than represented in current climate models. Since coarse dust warms global climate, the dust direct radiative effect (DRE) is likely less cooling than the ~0.4 W/m(2) estimated by models in a current ensemble(2,11–13). We constrain the dust DRE to - 0.20 (−0.48 to +0.20) W/m(2), which suggests that the dust DRE produces only about half the cooling that current models estimate, and raises the possibility that dust DRE is actually net warming the planet.