Search

Desert dust is an important atmospheric aerosol that affects the Earth's climate, biogeochemistry, and air quality. However, current Earth system models (ESMs) struggle to accurately capture the impact of dust on the Earth's climate and ecosystems, in part because these models lack several essential aeolian processes that couple dust with climate and land surface processes. In this study, we address this issue by implementing several new parameterizations of aeolian processes detailed in our companion paper in the Community Earth System Model version 2 (CESM2). These processes include (1) incorporating a simplified soil particle size representation to calculate the dust emission threshold friction velocity, (2) accounting for the drag partition effect of rocks and vegetation in reducing wind stress on erodible soils, (3) accounting for the intermittency of dust emissions due to unresolved turbulent wind fluctuations, and (4) correcting the spatial variability of simulated dust emissions from native to higher spatial resolutions on spatiotemporal dust variability. Our results show that the modified dust emission scheme significantly reduces the model bias against observations compared with the default scheme and improves the correlation against observations of multiple key dust variables such as dust aerosol optical depth (DAOD), surface particulate matter (PM) concentration, and deposition flux. Our scheme's dust also correlates strongly with various meteorological and land surface variables, implying higher sensitivity of dust to future climate change than other schemes' dust. These findings highlight the importance of including additional aeolian processes for improving the performance of ESM aerosol simulations and potentially enhancing model assessments of how dust impacts climate and ecosystem changes.

Estimating past aerosol radiative effects and their uncertainties is an important topic in climate science. Aerosol radiative effects propagate into large uncertainties in estimates of how present and future climate evolves with changing greenhouse gas emissions. A deeper understanding of how aerosols affected the atmospheric energy budget under past climates is hindered in part by a lack of relevant paleo-observations and in part because less attention has been paid to the problem. Because of the lack of information we do not seek here to determine the change in the radiative forcing due to aerosol changes but rather to estimate the uncertainties in those changes. Here we argue that current uncertainties from emission uncertainties (90 % confidence interval range spanning 2.8 W m−2) are just as large as model spread uncertainties (2.8 W m−2) in calculating preindustrial to present-day aerosol radiative effects. There are no estimates of radiative forcing for important aerosols such as wildfire and dust aerosols in most paleoclimate time periods. However, qualitative analysis of paleoclimate proxies suggests that changes in aerosols between different past climates are similar in magnitude to changes in aerosols between the preindustrial and present day; plus, there is the added uncertainty from the variability in aerosols and fires in the preindustrial. From the limited literature we crudely estimate a paleoclimate aerosol uncertainty for the Last Glacial Maximum relative to preindustrial of 4.8 W m−2, and we estimate the uncertainty in the aerosol feedback in the natural Earth system over the paleoclimate (Last Glacial Maximum to preindustrial) to be about 3.2 W m−2 K−1. In order to more accurately assess the uncertainty in historical aerosol radiative effects, we propose a new model intercomparison project, which would include multiple plausible emission scenarios tested across a range of state-of-the-art climate models over the historical period. These emission scenarios would then be compared to the available independent aerosol observations to constrain which are most probable. In addition, future efforts should work to characterize and constrain paleo-aerosol forcings and uncertainties. Careful propagation of aerosol uncertainties in the literature is required to ensure an accurate quantification of uncertainties in projections of future climate changes.

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 the dust lognormal volume particle size distribution (PSD) corresponding to retrieved Deff with the in situ-measured dust PSDs from the AERosol Properties – Dust (AER-D), Saharan Mineral Dust Experiment (SAMUM-2), and Saharan Aerosol Long-Range Transport and Aerosol–Cloud-Interaction Experiment (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 DAOD retrievals indicates a good agreement (R∼0.7) and a significant reduction in (∼50 %) retrieval uncertainties largely thanks to the better constraint on dust size. In a climatological comparison, the seasonal and regional (2∘×5∘) mean DAOD10 µm 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 is up to 20 % smaller than that over NA. Over NA in summer, we found a ∼50 % reduction in the number of retrievals with Deff>5 µm from 15 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 to 180∘ E, while the mean Deff remains stable at 4.0 µm throughout eastern NP. To the best of our 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 longwave radiative effects and coarse-mode dust particle size in models.

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 nonerodible 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 rescale lower-resolution dust emission simulations to match the spatial variability of higher-resolution emission simulations in 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.

Most global aerosol models approximate dust as spherical particles, whereas most remote sensing retrieval algorithms approximate dust as spheroidal particles with a shape distribution that conflicts with measurements. These inconsistent and inaccurate shape assumptions generate biases in dust single-scattering properties. Here, we obtain dust single-scattering properties by approximating dust as triaxial ellipsoidal particles with observationally constrained shape distributions. We find that, relative to the ellipsoidal dust optics obtained here, the spherical dust optics used in most aerosol models underestimate dust single-scattering albedo, mass extinction efficiency, and asymmetry parameter for almost all dust sizes in both the shortwave and longwave spectra. We further find that the ellipsoidal dust optics are in substantially better agreement with observations of the scattering matrix and linear depolarization ratio than the spheroidal dust optics used in most retrieval algorithms. However, relative to observations, the ellipsoidal dust optics overestimate the lidar ratio by underestimating the backscattering intensity by a factor of ∼2. This occurs largely because the computational method used to simulate ellipsoidal dust optics (i.e., the improved geometric optics method) underestimates the backscattering intensity by a factor of ∼2 relative to other computational methods (e.g., the physical geometric optics method). We conclude that the ellipsoidal dust optics with observationally constrained shape distributions can help improve global aerosol models and possibly remote sensing retrieval algorithms that do not use the backscattering signal.

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.

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 consistent uncertainties across spatiotemporal scales, which will enhance the use of the MIDAS dataset in a variety of future studies.

Mineral dust aerosols cool and warm the atmosphere by scattering and absorbing solar (shortwave: SW) and thermal (longwave: LW) radiation. However, significant uncertainties remain in dust radiative effects, largely due to differences in the dust size distribution and spectral optical properties simulated in Earth system models. Dust models typically underestimate the coarse dust load (more than 2.5 µm in diameter) and assume a spherical shape, which leads to an overestimate of the fine dust load (less than 2.5 µm) after the dust emissions in the models are scaled to match observed dust aerosol optical depth at 550 nm (DAOD550). Here, we improve the simulated dust properties with data sets that leverage measurements of size-resolved dust concentration, asphericity factor, and refractive index in a coupled global chemical transport model with a radiative transfer module. After the adjustment of size-resolved dust concentration and spectral optical properties, the global and annual average of DAOD550 from the simulation increases from 0.023 to 0.029 and falls within the range of a semi-observationally based estimate (0.030 ± 0.005). The reduction of fine dust load after the adjustment leads to a reduction of the SW cooling at the top of the atmosphere (TOA). To improve agreement against a semi-observationally based estimate of the radiative effect efficiency at TOA, we find that a less absorptive SW dust refractive index is required for coarser aspherical dust. Thus, only a minor difference is estimated for the net global dust radiative effect at TOA (−0.08 vs. −0.00 W m−2 on a global scale). Conversely, our sensitivity simulations reveal that the surface warming is substantially enhanced near the strong dust source regions (less cooling to −0.23 from −0.60 W m−2 on a global scale). Thus, less atmospheric radiative heating is estimated near the major source regions (less heating to 0.15 from 0.59 W m−2 on a global scale), because of enhanced LW warming at the surface by the synergy of coarser size and aspherical shape.

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

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.

The large uncertainty in the mineral dust direct radiative effect (DRE) hinders projections of future climate change due to anthropogenic activity. Resolving modeled dust mineral speciation allows for spatially and temporally varying refractive indices consistent with dust aerosol composition. Here, for the first time, we quantify the range in dust DRE at the top of the atmosphere (TOA) due to current uncertainties in the surface soil mineralogical content using a dust mineral-resolving climate model. We propagate observed uncertainties in soil mineral abundances from two soil mineralogy atlases along with the optical properties of each mineral into the DRE and compare the resultant range with other sources of uncertainty across six climate models. The shortwave DRE responds region-specifically to the dust burden depending on the mineral speciation and underlying shortwave surface albedo: positively when the regionally averaged annual surface albedo is larger than 0.28 and negatively otherwise. Among all minerals examined, the shortwave TOA DRE and single scattering albedo at the 0.44–0.63 µm band are most sensitive to the fractional contribution of iron oxides to the total dust composition. The global net (shortwave plus longwave) TOA DRE is estimated to be within −0.23 to +0.35 W m−2. Approximately 97 % of this range relates to uncertainty in the soil abundance of iron oxides. Representing iron oxide with solely hematite optical properties leads to an overestimation of shortwave DRE by +0.10 W m−2 at the TOA, as goethite is not as absorbing as hematite in the shortwave spectrum range. Our study highlights the importance of iron oxides to the shortwave DRE: they have a disproportionally large impact on climate considering their small atmospheric mineral mass fractional burden (∼2 %). An improved description of iron oxides, such as those planned in the Earth Surface Mineral Dust Source Investigation (EMIT), is thus essential for more accurate estimates of the dust DRE.

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.

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.

The role of mineral dust in climate and ecosystems has been largely quantified using global climate and chemistry model simulations of dust emission, transport, and deposition. However, differences between these model simulations are substantial, with estimates of global dust aerosol optical depth (AOD) that vary by over a factor of 5. Here we develop an observationally based estimate of the global dust AOD, using multiple satellite platforms, in situ AOD observations and four state-of-the-science global models over 2004–2008. We estimate that the global dust AOD at 550 nm is 0.030 ± 0.005 (1σ), higher than the AeroCom model median (0.023) and substantially narrowing the uncertainty. The methodology used provides regional, seasonal dust AOD and the associated statistical uncertainty for key dust regions around the globe with which model dust schemes can be evaluated. Exploring the regional and seasonal differences in dust AOD between our observationally based estimate and the four models in this study, we find that emissions in Africa are often overrepresented at the expense of Asian and Middle Eastern emissions and that dust removal appears to be too rapid in most models.

The mineralogy of desert dust is important due to its effect on radiation, clouds and biogeochemical cycling of trace nutrients. This study presents the simulation of dust radiative forcing as a function of both mineral composition and size at the global scale, using mineral soil maps for estimating emissions. Externally mixed mineral aerosols in the bulk aerosol module in the Community Atmosphere Model version 4 (CAM4) and internally mixed mineral aerosols in the modal aerosol module in the Community Atmosphere Model version 5.1 (CAM5) embedded in the Community Earth System Model version 1.0.5 (CESM) are speciated into common mineral components in place of total dust. The simulations with mineralogy are compared to available observations of mineral atmospheric distribution and deposition along with observations of clear-sky radiative forcing efficiency. Based on these simulations, we estimate the all-sky direct radiative forcing at the top of the atmosphere as + 0.05 Wm−2 for both CAM4 and CAM5 simulations with mineralogy. We compare this to the radiative forcing from simulations of dust in release versions of CAM4 and CAM5 (+0.08 and +0.17 Wm−2) and of dust with optimized optical properties, wet scavenging and particle size distribution in CAM4 and CAM5, −0.05 and −0.17 Wm−2, respectively. The ability to correctly include the mineralogy of dust in climate models is hindered by its spatial and temporal variability as well as insufficient global in situ observations, incomplete and uncertain source mineralogies and the uncertainties associated with data retrieved from remote sensing methods.

The complex nature of mineral dust aerosol emission makes it a difficult process to represent accurately in weather and climate models. Indeed, results in the companion paper indicate that many large-scale models underestimate the dust flux's sensitivity to the soil's threshold friction velocity for erosion. We hypothesize that this finding explains why many dust cycle simulations are improved by using an empirical dust source function that shifts emissions towards the world's most erodible regions. Here, we both test this hypothesis and evaluate the performance of the new dust emission parameterization presented in the companion paper. We do so by implementing the new emission scheme into the Community Earth System Model (CESM) and comparing the resulting dust cycle simulations against an array of measurements. We find that the new scheme shifts emissions towards the world's most erodible regions in a manner that is strikingly similar to the effect of implementing a widely used source function based on satellite observations of dust source regions. Furthermore, model comparisons against aerosol optical depth measurements show that the new physically based scheme produces a statistically significant improvement in CESM's representation of dust emission, which exceeds the improvement produced by implementing a source function. These results indicate that the need to use an empirical source function is eliminated, at least in CESM, by the additional physics in the new scheme, and in particular by its increased sensitivity to the soil's threshold friction velocity. Since the threshold friction velocity is affected by climate changes, our results further suggest that many large-scale models underestimate the global dust cycle's climate sensitivity.

Simulations of the dust cycle and its interactions with the changing Earth system are hindered by the empirical nature of dust emission parameterizations in weather and climate models. Here we take a step towards improving dust cycle simulations by using a combination of theory and numerical simulations to derive a physically based dust emission parameterization. Our parameterization is straightforward to implement into large-scale models, as it depends only on the wind friction velocity and the soil's threshold friction velocity. Moreover, it accounts for two processes missing from most existing parameterizations: a soil's increased ability to produce dust under saltation bombardment as it becomes more erodible, and the increased scaling of the dust flux with wind speed as a soil becomes less erodible. Our treatment of both these processes is supported by a compilation of quality-controlled vertical dust flux measurements. Furthermore, our scheme reproduces this measurement compilation with substantially less error than the existing dust flux parameterizations we were able to compare against. A critical insight from both our theory and the measurement compilation is that dust fluxes are substantially more sensitive to the soil's threshold friction velocity than most current schemes account for.

This study examines the uncertainties in simulating mass balance and radiative forcing of mineral dust due to biases in the dust size parameterization. Simulations are conducted quasi-globally (180° W–180° E and 60° S–70° N) using the WRF-Chem model with three different approaches to represent dust size distribution (8-bin, 4-bin, and 3-mode). The biases in the 3-mode or 4-bin approaches against a relatively more accurate 8-bin approach in simulating dust mass balance and radiative forcing are identified. Compared to the 8-bin approach, the 4-bin approach simulates similar but coarser size distributions of dust particles in the atmosphere, while the 3-mode approach retains more fine dust particles but fewer coarse dust particles due to its prescribed σg of each mode. Although the 3-mode approach yields up to 10 days of longer dust mass lifetime over the remote oceanic regions than the 8-bin approach, the three size approaches produce a similar dust mass lifetime (3.2 days to 3.5 days) on quasi-global average, reflecting that the global dust mass lifetime is mainly determined by the dust mass lifetime near the dust source regions. With the same global dust emission (~4600 Tg yr−1), the 8-bin approach produces a dust mass loading of 39 Tg, while the 4-bin and 3-mode approaches produce 3% (40.2 Tg) and 25% (49.1 Tg) higher dust mass loading, respectively. The difference in dust mass loading between the 8-bin approach and the 4-bin or 3-mode approaches has large spatial variations, with generally smaller relative difference (−2) and atmospheric warming (0.39~0.96 W m−2) and in a tremendous difference of a factor of ~10 in dust TOA (top of atmosphere) cooling (−0.24~−2.20 W m−2). The impact of different size representations on dust radiative forcing efficiency is smaller. An uncertainty of a factor of 2 is quantified in dust emission estimation due to the different size parameterizations. This study also highlights the uncertainties in modeling dust mass and number loading, deposition fluxes, and radiative forcing resulting from different size parameterizations, and motivates further investigation of the impact of size parameterizations on modeling dust impacts on air quality, climate, and ecosystems.

The present study investigates the dust emission and load over the Mediterranean basin using the coupled chemistry–aerosol–regional climate model RegCM-4. The first step of this work focuses on dust particle emission size distribution modeling. We compare a parameterization in which the emission is based on the individual kinetic energy of the aggregates striking the surface to a recent parameterization based on an analogy with the fragmentation of brittle materials. The main difference between the two dust schemes concerns the mass proportion of fine aerosol that is reduced in the case of the new dust parameterization, with consequences for optical properties. At the episodic scale, comparisons between RegCM-4 simulations, satellite and ground-based data show a clear improvement using the new dust distribution in terms of aerosol optical depth (AOD) values and geographic gradients. These results are confirmed at the seasonal scale for the investigated year 2008. This change of dust distribution has sensitive impacts on the simulated regional dust budget, notably dry dust deposition and the regional direct aerosol radiative forcing over the Mediterranean basin. In particular, we find that the new size distribution produces a higher dust deposition flux, and smaller top of atmosphere (TOA) dust radiative cooling. A multi-annual simulation is finally carried out using the new dust distribution over the period 2000–2009. The average SW radiative forcing over the Mediterranean Sea reaches −13.6 W m−2 at the surface, and −5.5 W m−2 at TOA. The LW radiative forcing is positive over the basin: 1.7 W m−2 on average over the Mediterranean Sea at the surface, and 0.6 W m−2 at TOA.

The size distribution of mineral dust aerosols partially determines their interactions with clouds, radiation, ecosystems, and other components of the Earth system. Several theoretical models predict that the dust size distribution depends on the wind speed at emission, with larger wind speeds predicted to produce smaller aerosols. The present study investigates this prediction using a compilation of published measurements of the size-resolved vertical dust flux emitted by eroding soils. Surprisingly, these measurements indicate that the size distribution of naturally emitted dust aerosols is independent of the wind speed. The recently formulated brittle fragmentation theory of dust emission is consistent with this finding, whereas other theoretical models are not. The independence of the emitted dust size distribution with wind speed simplifies both the interpretation of geological records of dust deposition and the parameterization of dust emission in atmospheric circulation models.