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⁻¹, 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.

Aerosols from biomass burning (BB) emissions are poorly constrained in global and regional models, resulting in a high level of uncertainty in understanding their impacts. In this study, we compared six BB aerosol emission datasets for 2008 globally as well as in 14 regions. The six BB emission datasets are (1) GFED3.1 (Global Fire Emissions Database version 3.1), (2) GFED4s (GFED version 4 with small fires), (3) FINN1.5 (FIre INventory from NCAR version 1.5), (4) GFAS1.2 (Global Fire Assimilation System version 1.2), (5) FEER1.0 (Fire Energetics and Emissions Research version 1.0), and (6) QFED2.4 (Quick Fire Emissions Dataset version 2.4). The global total emission amounts from these six BB emission datasets differed by a factor of 3.8, ranging from 13.76 to 51.93 Tg for organic carbon and from 1.65 to 5.54 Tg for black carbon. In most of the regions, QFED2.4 and FEER1.0, which are based on satellite observations of fire radiative power (FRP) and constrained by aerosol optical depth (AOD) data from the Moderate Resolution Imaging Spectroradiometer (MODIS), yielded higher BB aerosol emissions than the rest by a factor of 2–4. By comparison, the BB aerosol emissions estimated from GFED4s and GFED3.1, which are based on satellite burned-area data, without AOD constraints, were at the low end of the range. In order to examine the sensitivity of model-simulated AOD to the different BB emission datasets, we ingested these six BB emission datasets separately into the same global model, the NASA Goddard Earth Observing System (GEOS) model, and compared the simulated AOD with observed AOD from the AErosol RObotic NETwork (AERONET) and the Multiangle Imaging SpectroRadiometer (MISR) in the 14 regions during 2008. In Southern Hemisphere Africa (SHAF) and South America (SHSA), where aerosols tend to be clearly dominated by smoke in September, the simulated AOD values were underestimated in almost all experiments compared to MISR, except for the QFED2.4 run in SHSA. The model-simulated AOD values based on FEER1.0 and QFED2.4 were the closest to the corresponding AERONET data, being, respectively, about 73 % and 100 % of the AERONET observed AOD at Alta Floresta in SHSA and about 49 % and 46 % at Mongu in SHAF. The simulated AOD based on the other four BB emission datasets accounted for only ∼50 % of the AERONET AOD at Alta Floresta and ∼20 % at Mongu. Overall, during the biomass burning peak seasons, at most of the selected AERONET sites in each region, the AOD values simulated with QFED2.4 were the highest and closest to AERONET and MISR observations, followed closely by FEER1.0. However, the QFED2.4 run tends to overestimate AOD in the region of SHSA, and the QFED2.4 BB emission dataset is tuned with the GEOS model. In contrast, the FEER1.0 BB emission dataset is derived in a more model-independent fashion and is more physically based since its emission coefficients are independently derived at each grid box. Therefore, we recommend the FEER1.0 BB emission dataset for aerosol-focused hindcast experiments in the two biomass-burning-dominated regions in the Southern Hemisphere, SHAF, and SHSA (as well as in other regions but with lower confidence). The differences between these six BB emission datasets are attributable to the approaches and input data used to derive BB emissions, such as whether AOD from satellite observations is used as a constraint, whether the approaches to parameterize the fire activities are based on burned area, FRP, or active fire count, and which set of emission factors is chosen.

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 (PM₂₀) is approximately 5000 Tg yr⁻¹, which is greater than most models account for. This larger PM₂₀ 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.

A radiative transfer interface has been developed to simulate the UV aerosol index (AI) from the NASA Goddard Earth Observing System version 5 (GEOS-5) aerosol assimilated fields. The purpose of this work is to use the AI and aerosol absorption optical depth (AAOD) derived from the Ozone Monitoring Instrument (OMI) measurements as independent validation for the Modern Era Retrospective analysis for Research and Applications Aerosol Reanalysis (MERRAero). MERRAero is based on a version of the GEOS-5 model that is radiatively coupled to the Goddard Chemistry, Aerosol, Radiation, and Transport (GOCART) aerosol module and includes assimilation of aerosol optical depth (AOD) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. Since AI is dependent on aerosol concentration, optical properties and altitude of the aerosol layer, we make use of complementary observations to fully diagnose the model, including AOD from the Multi-angle Imaging SpectroRadiometer (MISR), aerosol retrievals from the AErosol RObotic NETwork (AERONET) and attenuated backscatter coefficients from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) mission to ascertain potential misplacement of plume height by the model. By sampling dust, biomass burning and pollution events in 2007 we have compared model-produced AI and AAOD with the corresponding OMI products, identifying regions where the model representation of absorbing aerosols was deficient. As a result of this study over the Saharan dust region, we have obtained a new set of dust aerosol optical properties that retains consistency with the MODIS AOD data that were assimilated, while resulting in better agreement with aerosol absorption measurements from OMI. The analysis conducted over the southern African and South American biomass burning regions indicates that revising the spectrally dependent aerosol absorption properties in the near-UV region improves the modeled-observed AI comparisons. Finally, during a period where the Asian region was mainly dominated by anthropogenic aerosols, we have performed a qualitative analysis in which the specification of anthropogenic emissions in GEOS-5 is adjusted to provide insight into discrepancies observed in AI comparisons.

The radiative effects of Saharan dust aerosols are investigated in the NASA GEOS‐5 atmospheric general circulation model. A sectional aerosol microphysics model (CARMA) is run online in GEOS‐5. CARMA treats the dust aerosol lifecycle, and its tracers are radiatively coupled to GEOS‐5. A series of AMIP‐style simulations are performed, in which input dust optical properties (particle shape and refractive index) are varied. Simulated dust distributions for summertime Saharan dust compare well to observations, with best results found when the most absorbing dust optical properties are assumed. Dust absorption leads to a strengthening of the summertime Hadley cell circulation, increased dust lofting to higher altitudes, and a strengthening of the African easterly jet, resulting in increased dust atmospheric lifetime and farther northward and westward transport. We find a positive feedback of dust radiative forcing on emissions, in contrast with previous studies, which we attribute to our having a relatively strong longwave forcing caused by our simulating larger effective particle sizes. This longwave forcing reduces the magnitude of midday net surface cooling relative to other studies, and leads to a nighttime warming that results in higher nighttime wind speeds and dust emissions. The radiative effects of dust particle shape have only minor impact on transport and emissions, with small (~5%) impact on top of atmosphere shortwave forcing, in line with previous studies, but relatively more pronounced effects on shortwave atmospheric heating and surface forcing (~20% increase in atmospheric forcing for spheroids). Shape effects on longwave heating terms are of order ~10%.