Articles | Volume 20, issue 23
Atmos. Chem. Phys., 20, 15285–15295, 2020

Special issue: The Modular Earth Submodel System (MESSy) (ACP/GMD inter-journal...

Atmos. Chem. Phys., 20, 15285–15295, 2020

Research article 09 Dec 2020

Research article | 09 Dec 2020

Weaker cooling by aerosols due to dust–pollution interactions

Weaker cooling by aerosols due to dust–pollution interactions
Klaus Klingmüller1, Vlassis A. Karydis2, Sara Bacer3, Georgiy L. Stenchikov4, and Jos Lelieveld1,5 Klaus Klingmüller et al.
  • 1Max Planck Institute for Chemistry, Hahn-Meitner-Weg 1, 55128 Mainz, Germany
  • 2Forschungszentrum Jülich GmbH, IEK-8, 52425 Jülich, Germany
  • 3LEGI, Université Grenoble Alpes, CNRS, Grenoble INP, Grenoble, France
  • 4King Abdullah University of Science and Technology, Thuwal 23955-6900, Saudi Arabia
  • 5The Cyprus Institute, P.O. Box 27456, 1645 Nicosia, Cyprus

Correspondence: Klaus Klingmüller (


The interactions between aeolian dust and anthropogenic air pollution, notably chemical ageing of mineral dust and coagulation of dust and pollution particles, modify the atmospheric aerosol composition and burden. Since the aerosol particles can act as cloud condensation nuclei, this affects the radiative transfer not only directly via aerosol–radiation interactions, but also indirectly through cloud adjustments. We study both radiative effects using the global ECHAM/MESSy atmospheric chemistry-climate model (EMAC) which combines the Modular Earth Submodel System (MESSy) with the European Centre/Hamburg (ECHAM) climate model. Our simulations show that dust–pollution–cloud interactions reduce the condensed water path and hence the reflection of solar radiation. The associated climate warming outweighs the cooling that the dust–pollution interactions exert through the direct radiative effect. In total, this results in a net warming by dust–pollution interactions which moderates the negative global anthropogenic aerosol forcing at the top of the atmosphere by (0.2 ± 0.1) W m−2.

1 Introduction

A prime objective of current atmospheric and climate science is the deeper understanding of ambient aerosols and their interactions with clouds. This is motivated by their central role in two areas of societal importance, public health and climate change. The inhalation of aerosols allows fine particles to enter deep into the respiratory system or even translocate through the lungs into the cardiovascular system, causing a multitude of health challenges and making exposure to fine particulate air pollution one of the main public health risks worldwide (Lelieveld et al.2015, 2019a, b; Chowdhury et al.2020). On the other hand, aerosols modify the albedo of the Earth, predominantly increasing the reflection of solar radiation and thus cooling the planet. Since the emissions of anthropogenic greenhouse gases are accompanied by those of aerosols, to a large extent through common source categories, the greenhouse warming has been partially masked by the aerosol effects on climate (IPCC2014; Lelieveld et al.2019a).

The planetary albedo can be increased both directly by interactions of the anthropogenic aerosol particles with solar radiation and indirectly by enhanced cloudiness or cloud brightness caused by aerosol particles acting as cloud condensation nuclei. These direct and indirect effects are estimated to contribute a negative effective radiative forcing (ERF) of 0.45 W m−2 each, adding up to about 0.9 W m−2 (IPCC2014).

Since not all aerosols in the atmosphere are of anthropogenic origin – in fact natural aerosols including aeolian dust and sea salt are the most abundant components by mass – the anthropogenic pollutants form a mixture with natural aerosols. On the one hand particulate pollution coagulates with natural particles and on the other hand natural particles are exposed to chemical ageing. Klingmüller et al. (2019) showed that the interactions between natural mineral dust and anthropogenic pollution enhance the global net cooling through the direct radiative effects and have a significant impact on regional radiative transfer. Here we extend the analysis to include the indirect radiative effects. The abundant atmospheric water vapour represents a vast source of cloud water so that cloud optical depths are typically much larger than aerosol optical depths. Therefore, cloud adjustments potentially leverage the aerosol radiative effect, and we may expect the indirect radiative effect of the dust–pollution interactions to be even more significant than the direct effect.

We use the global ECHAM/MESSy atmospheric chemistry-climate model (EMAC) which combines the Modular Earth Submodel System (MESSy) with the European Centre/Hamburg (ECHAM) climate model. It includes implementations of an extensive set of relevant physical and chemical processes, including detailed parametrisations of mineral dust ageing, cloud droplet activation and ice crystal formation in cirrus and mixed-phase clouds.

The model and its configuration are described in Sect. 2 followed by an outline of the methodology of our analysis in Sect. 3. Results for the dust–pollution interaction effect on the cloud condensate are presented in Sect. 4 and the resulting effects on radiative transfer in Sect. 5. Conclusions are presented in Sect. 6.

2 Model description

The EMAC model version and configuration used in the present study are largely identical to those used by Klingmüller et al. (2018, 2019), combining ECHAM 5.3.02 and MESSy 2.52. However, to allow decadal simulations, the horizontal resolution has been reduced to a Gaussian T63 grid with a grid spacing of 1.875 along latitudes and about 1.86 along longitudes, corresponding to an edge length of the individual grid cells of around 200 km or less. The number of vertical levels remains at 31. Moreover, the present study uses the EDGARv4.3 (Emissions Database for Global Atmospheric Research) database for anthropogenic emissions (Crippa et al.2016) and a back port of the CLOUD submodel from MESSy 2.54 to benefit from recent improvements of the cloud parametrisations.

As in the previous studies, the GFEDv3.1 (Global Fire Emissions Database) (Randerson et al.2013) and AeroCom (Aerosol Comparisons between Observations and Models) (Dentener et al.2006) databases provide biomass burning and sea salt emissions, respectively. Mineral dust emissions are calculated online by the submodel ONEMIS (Kerkweg et al.2006b) using the dust emission scheme presented by Klingmüller et al. (2018) which is based on Astitha et al. (2012). It differentiates the Ca++, K+, Mg++ and Na+ fractions in mineral particles originating from different deserts (Karydis et al.2016).

The MESSy submodels most relevant for aerosols include the Global Modal Aerosol Extension (GMXe) (Pringle et al.2010a, b). It simulates the microphysics of four soluble (nucleation, Aitken, accumulation, coarse) and three insoluble (Aitken, accumulation, coarse) aerosol log-normal modes with fixed geometric standard deviations (σg=2 for the coarse modes, σg=1.59 for all others). The count median dry radius of each mode can vary between fixed boundaries at 6, 60 nm and 1 µm. Super coarse mineral dust particles are therefore only included as part of the coarse modes with mean radius larger than 1 µm, and their mass is probably underrepresented (Adebiyi and Kok2020). However, their role in the dust–pollution–cloud interactions is limited by their low number concentration, corresponding to a low probability of pollution particles coagulating with them, and a comparably short atmospheric residence time, leaving less time for chemical ageing.

Within GMXe, the gas–aerosol partitioning can be computed by ISORROPIA II (Fountoukis and Nenes2007) or EQSAM4clim (Equilibrium Simplified Aerosol Model V4 for climate simulations) (Metzger et al.2016); here we use ISORROPIA II. Assuming diffusion-limited condensation, it calculates the amount of gas kinetically able to condense using the accommodation coefficients in Table S1 in the Supplement. Subsequently the mass is re-distributed between the gas and aerosol phases to obtain the amount of condensed material (Pringle et al.2010a, b). This means that gaseous compounds from anthropogenic pollution, including sulfuric acid, nitric acid, hydrochloric acid and ammonia, can condense on mineral dust particles and initiate their chemical ageing, which is the primary interaction between mineral dust and gaseous pollution, primarily through reactions of acids with mineral cations.

Insoluble particles are transferred to the soluble modes if sufficient hydrophilic material has accumulated to cover the particles with 10 molecular monolayers or if they coagulate with soluble particles (Vignati et al.2004; Stier et al.2005; Pringle et al.2010a, b). In particular, freshly emitted mineral dust is assumed to be hydrophobic and thus emitted into the insoluble aerosol modes (at approximately 89 %, the majority of the mass is emitted into the coarse mode and the remainder into the accumulation mode), but chemical ageing can transfer the mineral dust particles to soluble modes. The chemical ageing and partitioning of organic aerosol compounds are implemented in the submodel ORACLE (Organic Aerosol Composition and Evolution) (Tsimpidi et al.2014, 2018). GMXe and ORACLE interact with the gas phase, for which the chemistry is simulated by the submodel MECCA (Module Efficiently Calculating the Chemistry of the Atmosphere) (Sander et al.2019).

In addition to the condensation of hydrophilic compounds, the coagulation with soluble particles also transfers insoluble particles to the soluble modes. Within GMXe, coagulation is implemented following Vignati et al. (2004) using the coagulation coefficient equation from Fuchs (1964). All aerosol components are affected by coagulation irrespective of their sources and their chemical composition, including components represented by “bulk” tracers, which are treated as chemically inert, and the major particulate pollutants black carbon, organic compounds, sulfates, nitrate and ammonium. This makes coagulation the primary interaction between mineral dust and particulate pollution. Aside from modifying the composition and hygroscopicity of dust particles, it has a significant effect on the burden of particulate pollution. Because typically mineral dust particles are coarser than pollution particles like soot or sulfate particles, coagulation with dust transfers fine particulate pollution to coarser modes, decreasing the number concentration, especially in the fine modes. Once in the coarse mode, the pollution is affected by the shorter atmospheric residence time of coarse particles, which reduces the mass concentration of particulate pollutants. After being transferred to the hydrophilic modes, mineral dust particles grow by taking up water and act as cloud condensation nuclei. The hygroscopic growth increases the deposition rate and affects the optical properties.

Table 1Globally averaged annual mean cloud properties and contributions thereto, based on the SST simulations. “Total” represents the simulation with all emissions, “Mineral dust” and “Anthropogenic pollution” include the effect of dust–pollution interactions (Δdustxintx and Δpolxintx in Eq. 1), and “Dust–pollution interactions” are given by the interaction term Δintx, Eq. (2). The corresponding results from the nudged simulations are provided in Table S4 in the Supplement.

Download Print Version | Download XLSX

The AEROPT (AERosol OPTical properties) submodel (Lauer et al.2007; Klingmüller et al.2014) calculates the aerosol optical properties assuming the aerosol components within each mode to be well mixed in spherical particles with a volume-averaged refractive index. The refractive indices considered by AEROPT for the individual components are compiled from the OPAC 3.1 database (Hess et al.1998) (black carbon, mineral dust), the HITRAN 2004 database (Rothman et al.2005) (organic carbon, sea salt, ammonium sulfate, water), Kirchstetter et al. (2004) (organic carbon for λ<0.7µm) and additional mineral dust values for λ>2.5µm. The full dataset is specified in the Supplement of Klingmüller et al. (2014). The imaginary part of the dust refractive index provided by the OPAC dataset attains a minimum of 4×10-3 at visible and near-infrared wavelengths. This is lower than the former recommendation by the World Meteorological Organisation of 8×10-3 (Deepak et al.1983), but even smaller and regionally varying values are found in more recent literature (Kaufman et al.2001; Müller et al.2011; Di Biagio et al.2019). Even though a larger imaginary part of the refractive index corresponds to stronger absorption, we obtain a distinctive negative climate forcing attributed to mineral dust (Table 2). In our simulations the modelled dust is usually internally mixed with other components, and especially water, so that the effective imaginary refractive index of the entire particles is often lower than the value assumed for pure dust. Using a smaller value for pure dust would further enhance the negative direct forcing of dust and to some extent the direct forcing through the dust–pollution interactions. However, the dominant indirect effect of the interactions would not be affected.

The aerosol optical properties are considered by the radiative transfer submodel RAD (Dietmüller et al.2016) to account for the aerosol–radiation coupling. In the solar spectrum, absorption and scattering are computed using extinction coefficient, single scattering albedo and asymmetry parameter, whereas in the terrestrial spectrum scattering is neglected. The latter approximation is valid for particles much smaller than the wavelength and is therefore largely justified for the long terrestrial wavelengths but might be inaccurate in the presence of super coarse particles (Di Biagio et al.2020). However, this only affects the direct radiative effect. In the context of the present study, the indirect radiative effect turns out to be much more relevant.

Aerosol removal by wet deposition is calculated by the scavenging submodel SCAV (Tost et al.2006a), dry deposition and sedimentation by the submodels DDEP and SEDI (Kerkweg et al.2006a). The aerosol and in particular the mineral dust representation in EMAC have a proven track record (e.g. Abdelkader et al.2015, 2017; Metzger et al.2016; Klingmüller et al.2018; Brühl et al.2018; Metzger et al.2018; Ma et al.2019): the dust aerosol optical depth is consistent with observations not only at visible wavelengths, but also in the infrared at 10 µm (Klingmüller et al.2018). This is an indication of a realistic particle size distribution, provided that the ratio of the extinction efficiencies at visible and infrared wavelengths is realistic. The uncertainty in this ratio is expected to be small compared to other uncertainties, given that the spectral extinction efficiency is calculated consistently throughout the spectrum and, unlike the single scattering albedo, is hardly sensitive to the aforementioned uncertainties in the imaginary part of the refractive index.

Large-scale clouds are simulated by the submodel CLOUD (Jöckel et al.2006), where different parametrisations of cloud droplet formation and ice nucleation are implemented. We use a two-moment stratiform cloud microphysics scheme (Lohmann et al.1999, 2007; Lohmann and Kärcher2002) in combination with the UAF (Unified Activation Framework) cloud droplet activation parametrisation (Kumar et al.2011; Karydis et al.2011, 2017). For the ice crystal formation we use the comprehensive parametrisation for cirrus and mixed-phase clouds implemented by Bacer et al. (2018) based on Barahona and Nenes (2009). Convective clouds are calculated by the CONVECT submodel (Jöckel et al.2006), where interactions with aerosols are not taken into account. CONVECT provides a choice of convection schemes (Tost et al.2006b), and here we use the scheme of Tiedtke (1989) including modifications by Nordeng (1994). The optical properties of clouds which serve as input for the radiative transfer submodel RAD are computed by the submodel CLOUDOPT (Dietmüller et al.2016). The model yields a global annual mean cloud liquid water path around 80 g m−2 (Tables 1 and S3 in the Supplement), which is well within the range of other climate model results (32 to 125 g m−2Lebsock and Su2014) and observations (30 to 90 g m−2Lohmann and Neubauer2018). Likewise, the modelled annual mean global cloud ice water path of about 15 g m−2 (Tables 1 and S3 in the Supplement) is consistent with results from other models (e.g. 14.8 g m−2Lohmann and Neubauer2018) and close to observed values (e.g. (25 ± 7) g m−2Li et al.2012).

Table 2Globally averaged annual mean TOA ERFs in W m−2, based on the SST simulations. “Mineral dust” and “Anthropogenic pollution” include the effect of dust–pollution interactions (Δdustxintx and Δpolxintx in Eq. 1); “Dust–pollution interactions” are given by the interaction term Δintx, Eq. (2). The corresponding forcings obtained from the nudged simulations are provided in Table S4 in the Supplement.

Download Print Version | Download XLSX

A complete list of the MESSy submodels used in our simulations is provided in Table S2 in the Supplement. Descriptions of each submodel and further references can be found online in the MESSy submodel list (MESSy2020).

3 Methodology

We apply a similar analysis to Klingmüller et al. (2019), which is based on simulations with four different emission set-ups: a baseline simulation with neither dust nor anthropogenic emissions (“0”), a simulation with dust but without anthropogenic emissions (“dust”), a simulation with anthropogenic pollution but without dust emissions (“pol”) and a full simulation considering all emissions.

In the anthropogenic pollution-free simulations (“0”, “dust”) we disable the EDGAR emissions including SO2, NH3, NOx, and black and organic carbon emissions but retain the greenhouse gases. We attribute 90 % of the GFED biomass burning emissions to human activities (Levine2014) and reduce them accordingly, whereas we do not consider anthropogenic factors in dust emissions such as land use and climate change (Klingmüller et al.2016), assuming all dust emissions to be natural.

A result x from the full simulation (e.g. the annual global mean cloud liquid water content) is related to the corresponding result from the baseline simulation x0 by

(1) x = x 0 + Δ dust x + Δ pol x + Δ int x ,

where Δdustx=xdust-x0, Δpolx=xpol-x0 and

(2) Δ int x = x - x dust - x pol + x 0 ,

which represents the effect of the dust–pollution interactions. In the absence of such interactions, the term Δintx vanishes. We apply Eq. (2) to the annual mean cloud liquid and ice water paths and radiative fluxes.

To quantify the effects of the different emission set-ups and the dust–pollution interactions on radiation, we consider the ERF, which is defined as the change in net TOA downward radiative flux after allowing for atmospheric temperatures, water vapour and clouds to adjust but with sea surface temperature (SST) and sea ice cover fixed at climatological values (IPCC2014). Note that positive downward fluxes correspond to downward (incoming) radiation; negative values correspond to upward (outgoing) radiation. The ERF accounts for rapid adjustments by radiative and dynamical feedbacks, whereas it excludes long-term climate responses involving the much slower thermal equilibration of the oceans. Due to the limited constraints on the atmospheric dynamics in SST simulations, the meteorological variability is large, and hence a sufficient number of years has to be simulated to obtain statistically significant results. We perform SST simulations long enough to yield significant globally averaged results; however, detailed regional analysis would require much longer SST simulations. In order to nevertheless gain insights from regional evaluation, we additionally use simulations where the model dynamics above the boundary layer is nudged to meteorological analyses of the European Centre for Medium-Range Weather Forecasts (ECMWF). Within the boundary layer, in the topmost layers and to some extent in between, nudged quantities like the temperature may still respond to other variables such as radiative fluxes (soft nudging). The nudging greatly reduces the influence of inter-annual variability on statistical analysis. The results from the nudged simulations turn out to be largely consistent with those of the SST simulations (Tables 1, 2 vs. S3, S4 in the Supplement); in particular, the estimates for the total global radiative effect of the dust–pollution interactions agree within the error bounds, so that the use of nudged simulations for the regional analysis is reasonable and helpful.

With prescribed SST we run ensembles of 16 simulations, each covering 1 year. As there is one ensemble for each of the four emission set-ups, in total this amounts to 64 SST simulations. The ensemble members are obtained by perturbing temperature and humidity in the fourth year of a common spin-up simulation, followed by an additional spin-up of the individual ensemble members to attain a total of 5 spin-up years. The perturbation is implemented by adding a uniformly distributed random variable ranging from 0.1 ‰ to 0.1 ‰ of the perturbed quantity, so that the perturbation is numerically but not meteorologically relevant. Emission data for 2010 are used for all simulations. The nudged simulations cover 10 years from 2006 to 2015, and 2 simulation years prior to that period were used for the model spin-up. To estimate the uncertainties of the 10-year mean values for the nudged simulations and the ensemble mean values for the SST simulations, we compute the standard error of the mean (SEM) of the annual values.

In the analysis of the SST results, we substitute the variable x in Eq. (2) with global annual mean values, and for the nudged results we skip the global averaging and apply the equation to the annual mean for each grid cell separately to obtain the spatial distribution of the interaction term.

Substituting x for the global annual mean net flux F at the top of the atmosphere (TOA) in the SST simulations, ΔdustF=Fdust-F0 corresponds to the total ERF of mineral dust including all rapid adjustments, analogous to ΔpolF to the anthropogenic aerosol ERF (both excluding the dust–pollution interactions) and ΔintF to the ERF of dust–pollution interactions. In the case of the nudged simulations, the possible adjustments are constrained, so that the resulting forcings are in between the ERF and the radiative forcing RF as defined by IPCC (2014), where only the stratospheric temperature is allowed to adjust. For this reason the forcings from the nudged simulations are not directly comparable to RF and ERF results but, as mentioned above, provide valuable information about the regional effects.

To compute the direct radiative effect of aerosols, the radiative transfer code is called twice for every model time step. The first call considers scattering and absorption by aerosols and is used to calculate the heating rates affecting the temperature; the second call ignores scattering and absorption by aerosols and computes the radiative fluxes and heating rates only for diagnostic output. The difference of the radiative fluxes from both calls yields the instantaneous forcing (IRF) due to the direct radiative effect of aerosols Fari. Since both calls are performed with identical clouds, the cloud forcing is excluded, and only little statistical noise is introduced by the strong variability of clouds. Nevertheless, in this way we obtain the direct radiative forcing in the presence of clouds, which is typically smaller than the clear-sky forcing. The difference of the instantaneous aerosol forcings in the SST simulations with and without mineral dust ΔdustFari=Fari,dust-Fari,0 yields the aerosol–radiation interaction contribution to the ERF of dust, i.e. the direct radiative effect of dust. Analogously ΔpolFari and ΔintFari represent the direct radiative effect of particulate pollution and the dust–pollution interactions. The direct radiative forcings are subtracted from the corresponding total aerosol radiative forcings to extract the indirect radiative forcings; e.g. the indirect contribution to the dust–pollution interaction forcing is ΔintF−ΔintFari.

Figure 1Annual mean effect of the dust–pollution interactions on the liquid (a) and ice (b) cloud water, calculated by applying Eq. (2) to the results of the nudged simulations. Over stippled regions the results differ from zero by less than 2 times the standard error of the mean (SEM) of the annual values.

4 Effects on the cloud condensate

Hydrophilic particulate anthropogenic pollution enhances the cloud droplet formation and thus the liquid water content (Table 1). However, in the presence of mineral dust particles this effect is reduced because fine pollution particles coagulate with coarse dust particles, decreasing the particle number and virtually cleaning the atmosphere of fine particulate pollution. Moreover, the adsorption activation of mineral dust particles occurs early on in the cloud formation process (Kumar et al.2011), reducing the maximum supersaturation and inhibiting the activation of small pollution particles. These effects reduce the number of cloud condensation nuclei (Karydis et al.2011, 2017) and decrease the cloud liquid water path as shown in Fig. 1a. Especially over East and South Asia, where strong pollution emissions mix with aeolian dust from the Taklamakan, Gobi and Thar deserts, the reduction is substantial and regionally exceeds 40 g m−2. Even over polluted regions in Europe and the USA which are only occasionally exposed to dust intrusions, we obtain a small but significant reduction. This negative impact of the dust–pollution interactions over large parts of the Northern Hemisphere leads to a reduction of the global mean cloud liquid water path in Fig. 1a by (1.10 ± 0.03) g m−2. A comparable reduction by (1.5 ± 0.2) g m−2 is obtained in the SST simulations (Table 1). Relative to the mean liquid water path in the SST simulation considering all emissions (85.5 ± 0.1) g m−2, these reductions appear to be rather moderate. The reason is that the transport time periods between most of the major dust sources, especially the Sahara and the Middle East, and major pollution sources like North America and Europe to a large degree exceed the dust aerosol lifetime. In Asia these sources are less distant, while pollution emissions are generally larger. Thus, the strong effects over Asia might provide an outlook for regions with emerging pollution sources close to dust sources in Africa and the Middle East. But already today, due to the critical influence of clouds on radiative transfer, the relatively small changes in the water paths cause substantial radiative forcings, as will be discussed in the next section.

The dust–pollution interaction effect on the cloud ice water path, shown in Fig. 1b, is less distinct. A negative impact is obtained over the Sahel. The direct radiative effect of mineral dust over the Sahara warms the atmosphere by absorption of solar radiation (Fig. S1a in the Supplement). This increases the atmospheric capacity to hold moisture and the vertical water vapour transport (Fig. S1b in the Supplement). As a result, more moisture is available for ice cloud formation (Fig. S1c in the Supplement). Since the net direct radiative effect of the dust–pollution interactions cools the atmosphere over the Sahara (Klingmüller et al.2019), it moderates the enhancement of ice cloud formation. A similar net cooling effect is found over the region around the Taklamakan and Gobi deserts. In this region with generally high ice water content, anthropogenic pollution enhances the ice water path, but adding dust reduces the number of anthropogenic ice nucleation particles via coagulation. In contrast, a positive impact is obtained over coastal regions of Canada and Greenland around 60 N, probably due to aerosol and cloud feedbacks on the polar and Ferrel cell circulations and associated vertical moisture transport. However, because of the comparably small radiative fluxes at these latitudes, this has relatively little impact on radiative fluxes from a global perspective. Due to the regionally varying sign of the dust–pollution interaction effect on cloud ice, the global mean in Fig. 1b is close to zero, (0.027 ± 0.003) g m−2. The corresponding value in the SST simulations, (0.02 ± 0.03) g m−2, is consistent with this result, being several orders of magnitude smaller than the global mean ice water path in the SST simulation considering all emissions, (14.70 ± 0.01) g m−2 (Table 1).

Figure 2Annual mean indirect effect of the dust–pollution interactions on the solar (a) and terrestrial (b) radiative forcing at the top of the atmosphere, calculated by applying Eq. (2) to the results of the nudged simulations. Over stippled regions the results differ from zero by less than 2 times the SEM of the annual values.

5 Radiative effects

The reduction of the cloud water content by dust–pollution interactions has a significant impact on the transfer of solar radiation (“shortwave”, SW), which is shown in Fig. 2a. With reduced liquid cloud water, less solar radiation is reflected back to space; i.e. the outgoing radiation and the associated negative contribution to the net flux decrease, corresponding to a net positive forcing at the TOA. Comparing Figs. 2a and 1a reveals the one-to-one correspondence of the dust–pollution interaction effect on the liquid cloud water and solar radiation. Over the polluted regions of the Northern Hemisphere, i.e. Asia, Europe and North America, and over the Atlantic Ocean along the northern African coast in the Saharan dust outflow, the positive forcing can exceed 2 W m−2. Globally averaged, the net forcing in the solar spectrum shown in Fig. 2a is (0.23 ± 0.01) W m−2, and the SST simulations yield an ERF of (0.3 ± 0.1) W m−2 (Table 2).

Figure 3Total (direct and indirect, SW and LW) annual mean effect of the dust–pollution interactions on the radiative forcing at the top of the atmosphere, calculated by applying Eq. (2) to the results of the nudged simulations. Over stippled regions the results differ from zero by less than 2 times the SEM of the annual values.

On the other hand, the dust–pollution interaction effect on the terrestrial spectrum (“longwave”, LW), Fig. 2b, is directly related to the effects on ice clouds, Fig. 1b. This is most distinct over the Sahel but also apparent over the East Asian deserts. The reduced cloud ice water path over these regions traps less outgoing terrestrial radiation, resulting in a net cooling from the dust–pollution interactions. Over the Sahel the terrestrial TOA forcing reaches 2 W m−2. With regard to the radiative energy budget, the regions with a significant dust–pollution interaction effect on cloud ice in Fig. 1b are of different relevance. The Sahel, where the dust–pollution interactions reduce cloud ice, is relatively close to the Equator, and accordingly stronger radiative fluxes are affected by the cloud ice changes than in the other regions, hence the global net radiative effect related to cloud ice is more relevant than the global net effect on cloud ice itself. Globally averaged, the net forcing in the terrestrial spectrum shown in Fig. 2b is (0.05 ± 0.01) W m−2, and the SST simulations yield an ERF of (0.08 ± 0.09) W m−2 (Table 2).

Thus, a substantial positive forcing in the solar spectrum is partially compensated by a negative forcing in the terrestrial spectrum to yield a still considerable, positive net forcing associated with the effect of dust–pollution interactions on clouds. The global distribution of the total net forcing at the TOA including the direct radiative effect is shown in Fig. 3. The regional forcing ranges from below 2 W m−2 over the Sahel to above 2 W m−2 over Asia. Even though overall these contributions partially counterbalance, at (0.15 ± 0.02) W m−2 the corresponding global mean forcing in Fig. 3 is significantly positive. Consistently, the ERF in the SST simulations is (0.2 ± 0.1) W m−2 (Table 2).

Figure 4Estimates of the global anthropogenic aerosol forcings at the top of the atmosphere (TOA) in the presence (“With dust”, FFdust) or absence (“Without dust”, FpolF0) of aeolian dust, based on the SST simulations. The total forcings comprise the direct (green, Fari-Fari,dust and Fari,pol-Fari,0) and indirect (blue) forcings. The change caused by including mineral dust corresponds to the positive forcing of dust–pollution interactions ΔintF (red). Darker colours represent the SEM.


Figure 4 summarises the direct and indirect global TOA ERF of anthropogenic aerosol interacting with mineral dust and in the absence of mineral dust, obtained from the SST simulations. Despite the more negative anthropogenic aerosol direct radiative forcing in the presence of mineral dust, already reported by Klingmüller et al. (2019), the effect of mineral dust on the total forcing is clearly dominated by the moderation of the indirect forcing. The figure highlights the importance of the dust–pollution interactions for assessing the cooling effect of anthropogenic aerosol: the cooling is substantially reduced by the interactions from (0.81 ± 0.06) W m−2 ERF, which is close to 0.9 W m−2 estimated by IPCC (2014), down to (0.6 ± 0.1) W m−2.

6 Conclusions

We have studied the effects of interactions between mineral dust and anthropogenic pollution on clouds and radiation by analysing comprehensive global simulations performed with the atmospheric chemistry-climate model EMAC. Four different emission configurations representing all possible combinations of including and excluding dust and pollution were considered. Comparing the results for these four scenarios allowed us to isolate the effect of the dust–pollution interactions from the individual effects of dust and pollution. Several aspects make this analysis challenging, should be considered when interpreting the results, and may leave room for refinements in future studies. Naturally, clouds are subject to strong variability; hence, although we performed ensemble simulations, there is considerable statistical uncertainty in the present results. This adds to the need to evaluate differences in differences of results from a number of simulations to obtain the interaction effect, which increases the relative error. Moreover, a wide range of physical and chemical processes is involved in the dust–pollution interactions, and accordingly many submodels and parametrisations within EMAC contribute to our final result and uncertainty. Even though the parametrisations are well established and tested, the analysis might be sensitive to systematic errors of some of them.

The analysis reveals that the cloud water path is reduced by the dust–pollution interactions as they moderate the cloud water path increase caused by anthropogenic pollution. The reason for this moderation is that mineral dust particles decrease the number of anthropogenic cloud condensation nuclei by coagulation and additionally limit the activation of the fine hydrophilic anthropogenic particles by lowering the maximum supersaturation through adsorption activation. Dust–pollution interaction effects on the cloud ice content are noticeable as well but less relevant.

The atmospheric radiative transfer is very sensitive to the reduction of the cloud water path. Generally, dust–pollution interactions affect the radiative transfer at all wavelengths (solar and terrestrial) by modifying both the direct aerosol–radiation interactions and the indirect radiative effect of aerosols via cloud adjustments. However, the total radiative effect of the dust–pollution interactions is dominated by the impact through the indirect effect which, in contrast to the direct effect, exerts an overall positive TOA net forcing. The impact on the indirect radiative effect in turn is dominated by that on solar radiation fluxes. In this case, the aforementioned decrease in the cloud water path reduces the cloud albedo and the reflection of solar radiation, resulting in a positive contribution to the radiative net flux.

We estimate that dust–pollution interactions contribute (0.2 ± 0.1) W m−2 to the global mean anthropogenic aerosol effective radiative forcing, significantly reducing the climate cooling effect of atmospheric aerosols. In view of this considerable contribution to the atmospheric energy balance, it is recommended to account for the dust–pollution interactions in assessments of climate change, especially because on a regional scale effects can be even larger. The net global effect partially depends on regionally counteracting positive and negative radiative forcings. This study emphasises the importance of continued efforts to improve the understanding and parametrisations of the processes involved in order to reduce the uncertainty of future climate simulations.

Code and data availability

The Modular Earth Submodel System (MESSy) is continuously further developed and applied by a consortium of institutions. The usage of MESSy and access to the source code are licensed to all affiliates of institutions which are members of the MESSy Consortium. Institutions can become a member of the MESSy Consortium by signing the MESSy Memorandum of Understanding. More information can be found on the MESSy Consortium website (, last access: 30 October 2020, MESSy2020). The ECHAM climate model is available to the scientific community under the MPI-M Software License Agreement (, last access: 30 October 2020, MPI-C2020). The simulation results analysed in this study are available at (last access: 30 October 2020, Klingmüller2020).


The supplement related to this article is available online at:

Author contributions

KK performed the simulations assisted by VAK and SB, analysed the model results and wrote the article supported by JL, SB, VAK and GLS. All the authors discussed the results and contributed to the final manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “The Modular Earth Submodel System (MESSy) (ACP/GMD inter-journal SI)”. It is not associated with a conference.


The research reported in this publication has received funding from the MaxWater initiative of the Max Planck Society and the King Abdullah University of Science and Technology project “Combined Radiative and Air Quality Effects of Anthropogenic Air Pollution and Dust over the Arabian Peninsula”.

Financial support

This research has been supported by the King Abdullah University of Science and Technology (grant CRG3, grant no. URF/1/2180-01-01).

The article processing charges for this open-access
publication were covered by the Max Planck Society.

Review statement

This paper was edited by Yves Balkanski and reviewed by two anonymous referees.


Abdelkader, M., Metzger, S., Mamouri, R. E., Astitha, M., Barrie, L., Levin, Z., and Lelieveld, J.: Dust–air pollution dynamics over the eastern Mediterranean, Atmos. Chem. Phys., 15, 9173–9189,, 2015. a

Abdelkader, M., Metzger, S., Steil, B., Klingmüller, K., Tost, H., Pozzer, A., Stenchikov, G., Barrie, L., and Lelieveld, J.: Sensitivity of transatlantic dust transport to chemical aging and related atmospheric processes, Atmos. Chem. Phys., 17, 3799–3821,, 2017. a

Adebiyi, A. A. and Kok, J. F.: Climate models miss most of the coarse dust in the atmosphere, Sci. Adv., 6, eaaz9507,, 2020. a

Astitha, M., Lelieveld, J., Abdel Kader, M., Pozzer, A., and de Meij, A.: Parameterization of dust emissions in the global atmospheric chemistry-climate model EMAC: impact of nudging and soil properties, Atmos. Chem. Phys., 12, 11057–11083,, 2012. a

Bacer, S., Sullivan, S. C., Karydis, V. A., Barahona, D., Krämer, M., Nenes, A., Tost, H., Tsimpidi, A. P., Lelieveld, J., and Pozzer, A.: Implementation of a comprehensive ice crystal formation parameterization for cirrus and mixed-phase clouds in the EMAC model (based on MESSy 2.53), Geosci. Model Dev., 11, 4021–4041,, 2018. a

Barahona, D. and Nenes, A.: Parameterizing the competition between homogeneous and heterogeneous freezing in ice cloud formation – polydisperse ice nuclei, Atmos. Chem. Phys., 9, 5933–5948,, 2009. a

Brühl, C., Schallock, J., Klingmüller, K., Robert, C., Bingen, C., Clarisse, L., Heckel, A., North, P., and Rieger, L.: Stratospheric aerosol radiative forcing simulated by the chemistry climate model EMAC using Aerosol CCI satellite data, Atmos. Chem. Phys., 18, 12845–12857,, 2018. a

Chowdhury, S., Pozzer, A., Dey, S., Klingmueller, K., and Lelieveld, J.: Changing risk factors that contribute to premature mortality from ambient air pollution between 2000 and 2015, Environ. Res. Lett., 15, 074010,, 2020. a

Crippa, M., Janssens-Maenhout, G., Dentener, F., Guizzardi, D., Sindelarova, K., Muntean, M., Van Dingenen, R., and Granier, C.: Forty years of improvements in European air quality: regional policy-industry interactions with global impacts, Atmos. Chem. Phys., 16, 3825–3841,, 2016. a

Deepak, A. and Gerber, H. E.: Report of the Experts Meeting on Aerosols and Their Climatic Effects, World Meteorological Organization, Geneva, Switzerland, 1983. a

Dentener, F., Kinne, S., Bond, T., Boucher, O., Cofala, J., Generoso, S., Ginoux, P., Gong, S., Hoelzemann, J. J., Ito, A., Marelli, L., Penner, J. E., Putaud, J.-P., Textor, C., Schulz, M., van der Werf, G. R., and Wilson, J.: Emissions of primary aerosol and precursor gases in the years 2000 and 1750 prescribed data-sets for AeroCom, Atmos. Chem. Phys., 6, 4321–4344,, 2006. a

Di Biagio, C., Formenti, P., Balkanski, Y., Caponi, L., Cazaunau, M., Pangui, E., Journet, E., Nowak, S., Andreae, M. O., Kandler, K., Saeed, T., Piketh, S., Seibert, D., Williams, E., and Doussin, J.-F.: Complex refractive indices and single-scattering albedo of global dust aerosols in the shortwave spectrum and relationship to size and iron content, Atmos. Chem. Phys., 19, 15503–15531,, 2019. a

Di Biagio, C., Balkanski, Y., Albani, S., Boucher, O., and Formenti, P.: Direct Radiative Effect by Mineral Dust Aerosols Constrained by New Microphysical and Spectral Optical Data, Geophys. Res. Lett., 47, e2019GL086186,, 2020. a

Dietmüller, S., Jöckel, P., Tost, H., Kunze, M., Gellhorn, C., Brinkop, S., Frömming, C., Ponater, M., Steil, B., Lauer, A., and Hendricks, J.: A new radiation infrastructure for the Modular Earth Submodel System (MESSy, based on version 2.51), Geosci. Model Dev., 9, 2209–2222,, 2016. a, b

Fountoukis, C. and Nenes, A.: ISORROPIA II: a computationally efficient thermodynamic equilibrium model for K+-Ca2+-Mg2+-NH4+-Na+-SO42−-NO3-Cl-H2O aerosols, Atmos. Chem. Phys., 7, 4639–4659,, 2007. a

Fuchs, N. A.: The Mechanics of Aerosols, Pergamon Press, Oxford, 1964. a

Hess, M., Koepke, P., and Schult, I.: Optical Properties of Aerosols and Clouds: The Software Package OPAC, B. Am. Meteorol. Soc., 79, 831–844,<0831:OPOAAC>2.0.CO;2, 1998. a

IPCC: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, USA,, 2013. a, b, c, d, e

Jöckel, P., Tost, H., Pozzer, A., Brühl, C., Buchholz, J., Ganzeveld, L., Hoor, P., Kerkweg, A., Lawrence, M. G., Sander, R., Steil, B., Stiller, G., Tanarhte, M., Taraborrelli, D., van Aardenne, J., and Lelieveld, J.: The atmospheric chemistry general circulation model ECHAM5/MESSy1: consistent simulation of ozone from the surface to the mesosphere, Atmos. Chem. Phys., 6, 5067–5104,, 2006. a, b

Karydis, V. A., Kumar, P., Barahona, D., Sokolik, I. N., and Nenes, A.: On the effect of dust particles on global cloud condensation nuclei and cloud droplet number, J. Geophys. Res.-Atmos., 116, D23204,, 2011. a, b

Karydis, V. A., Tsimpidi, A. P., Pozzer, A., Astitha, M., and Lelieveld, J.: Effects of mineral dust on global atmospheric nitrate concentrations, Atmos. Chem. Phys., 16, 1491–1509,, 2016. a

Karydis, V. A., Tsimpidi, A. P., Bacer, S., Pozzer, A., Nenes, A., and Lelieveld, J.: Global impact of mineral dust on cloud droplet number concentration, Atmos. Chem. Phys., 17, 5601–5621,, 2017. a, b

Kaufman, Y. J., Tanré, D., Dubovik, O., Karnieli, A., and Remer, L. A.: Absorption of sunlight by dust as inferred from satellite and ground-based remote sensing, Geophys. Res. Lett., 28, 1479–1482,, 2001. a

Kerkweg, A., Buchholz, J., Ganzeveld, L., Pozzer, A., Tost, H., and Jöckel, P.: Technical Note: An implementation of the dry removal processes DRY DEPosition and SEDImentation in the Modular Earth Submodel System (MESSy), Atmos. Chem. Phys., 6, 4617–4632,, 2006a. a

Kerkweg, A., Sander, R., Tost, H., and Jöckel, P.: Technical note: Implementation of prescribed (OFFLEM), calculated (ONLEM), and pseudo-emissions (TNUDGE) of chemical species in the Modular Earth Submodel System (MESSy), Atmos. Chem. Phys., 6, 3603–3609,, 2006b. a

Kirchstetter, T. W., Novakov, T., and Hobbs, P. V.: Evidence that the spectral dependence of light absorption by aerosols is affected by organic carbon, J. Geophys. Res.-Atmos., 109, 21208,, 2004. a

Klingmüller, K.: Weaker cooling by aerosols due to dust-pollution interactions – simulation data, available at: (last access: 30 October 2020), 2020. a

Klingmüller, K., Steil, B., Brühl, C., Tost, H., and Lelieveld, J.: Sensitivity of aerosol radiative effects to different mixing assumptions in the AEROPT 1.0 submodel of the EMAC atmospheric-chemistry–climate model, Geosci. Model Dev., 7, 2503–2516,, 2014. a, b

Klingmüller, K., Pozzer, A., Metzger, S., Stenchikov, G. L., and Lelieveld, J.: Aerosol optical depth trend over the Middle East, Atmos. Chem. Phys., 16, 5063–5073,, 2016. a

Klingmüller, K., Metzger, S., Abdelkader, M., Karydis, V. A., Stenchikov, G. L., Pozzer, A., and Lelieveld, J.: Revised mineral dust emissions in the atmospheric chemistry–climate model EMAC (MESSy 2.52 DU_Astitha1 KKDU2017 patch), Geosci. Model Dev., 11, 989–1008,, 2018. a, b, c, d

Klingmüller, K., Lelieveld, J., Karydis, V. A., and Stenchikov, G. L.: Direct radiative effect of dust–pollution interactions, Atmos. Chem. Phys., 19, 7397–7408,, 2019. a, b, c, d, e

Kumar, P., Sokolik, I. N., and Nenes, A.: Cloud condensation nuclei activity and droplet activation kinetics of wet processed regional dust samples and minerals, Atmos. Chem. Phys., 11, 8661–8676,, 2011. a, b

Lauer, A., Eyring, V., Hendricks, J., Jöckel, P., and Lohmann, U.: Global model simulations of the impact of ocean-going ships on aerosols, clouds, and the radiation budget, Atmos. Chem. Phys., 7, 5061–5079,, 2007. a

Lebsock, M. and Su, H.: Application of active spaceborne remote sensing for understanding biases between passive cloud water path retrievals, J. Geophys. Res.-Atmos., 119, 8962–8979,, 2014. a

Lelieveld, J., Evans, J. S., Fnais, M., Giannadaki, D., and Pozzer, A.: The contribution of outdoor air pollution sources to premature mortality on a global scale, Nature, 525, 367–371,, 2015. a

Lelieveld, J., Klingmüller, K., Pozzer, A., Burnett, R. T., Haines, A., and Ramanathan, V.: Effects of fossil fuel and total anthropogenic emission removal on public health and climate, P. Natl. Acad. Sci., 116, 7192–7197,, 2019a. a, b

Lelieveld, J., Klingmüller, K., Pozzer, A., Pöschl, U., Fnais, M., Daiber, A., and Münzel, T.: Cardiovascular disease burden from ambient air pollution in Europe reassessed using novel hazard ratio functions, Eur. Heart J., 40, 1590–1596,, 2019b. a

Levine, J.: 5.5-Biomass Burning: The Cycling of Gases and Particulates from the Biosphere to the Atmosphere, in: Treatise on Geochemistry (Second Edition), edited by: Holland, H. D. and Turekian, K. K., Elsevier, Oxford, UK, 139–150,, 2014. a

Li, J.-L. F., Waliser, D. E., Chen, W.-T., Guan, B., Kubar, T., Stephens, G., Ma, H.-Y., Deng, M., Donner, L., Seman, C., and Horowitz, L.: An observationally based evaluation of cloud ice water in CMIP3 and CMIP5 GCMs and contemporary reanalyses using contemporary satellite data, J. Geophys. Res.-Atmos., 117, D16105,, 2012. a

Lohmann, U. and Kärcher, B.: First interactive simulations of cirrus clouds formed by homogeneous freezing in the ECHAM general circulation model, J. Geophys. Res.-Atmos., 107, D104105,, 2002. a

Lohmann, U. and Neubauer, D.: The importance of mixed-phase and ice clouds for climate sensitivity in the global aerosol–climate model ECHAM6-HAM2, Atmos. Chem. Phys., 18, 8807–8828,, 2018. a, b

Lohmann, U., Feichter, J., Chuang, C. C., and Penner, J. E.: Prediction of the number of cloud droplets in the ECHAM GCM, J. Geophys. Res.-Atmos., 104, 9169–9198,, 1999. a

Lohmann, U., Stier, P., Hoose, C., Ferrachat, S., Kloster, S., Roeckner, E., and Zhang, J.: Cloud microphysics and aerosol indirect effects in the global climate model ECHAM5-HAM, Atmos. Chem. Phys., 7, 3425–3446,, 2007. a

Ma, J., Brühl, C., He, Q., Steil, B., Karydis, V. A., Klingmüller, K., Tost, H., Chen, B., Jin, Y., Liu, N., Xu, X., Yan, P., Zhou, X., Abdelrahman, K., Pozzer, A., and Lelieveld, J.: Modeling the aerosol chemical composition of the tropopause over the Tibetan Plateau during the Asian summer monsoon, Atmos. Chem. Phys., 19, 11587–11612,, 2019. a

MESSy: The Modular Earth Submodel System, available at:, last access: 30 October 2020. a, b

Metzger, S., Steil, B., Abdelkader, M., Klingmüller, K., Xu, L., Penner, J. E., Fountoukis, C., Nenes, A., and Lelieveld, J.: Aerosol water parameterisation: a single parameter framework, Atmos. Chem. Phys., 16, 7213–7237,, 2016. a, b

Metzger, S., Abdelkader, M., Steil, B., and Klingmüller, K.: Aerosol water parameterization: long-term evaluation and importance for climate studies, Atmos. Chem. Phys., 18, 16747–16774,, 2018. a

Müller, T., Schladitz, A., Kandler, K., and Wiedensohler, A.: Spectral particle absorption coefficients, single scattering albedos and imaginary parts of refractive indices from ground based in situ measurements at Cape Verde Island during SAMUM-2, Tellus B, 63, 573–588,, 2011. a

MPI-C: Code availability, available at:, last access: 30 October 2020. a

Nordeng, T.-E.: Extended versions of the convective parametrization scheme at ECMWF and their impact on the mean and transient activity of the model in the tropics, ECMWF, Reading, UK,, 1994. a

Pringle, K. J., Tost, H., Message, S., Steil, B., Giannadaki, D., Nenes, A., Fountoukis, C., Stier, P., Vignati, E., and Lelieveld, J.: Description and evaluation of GMXe: a new aerosol submodel for global simulations (v1), Geosci. Model Dev., 3, 391–412,, 2010a. a, b, c

Pringle, K. J., Tost, H., Metzger, S., Steil, B., Giannadaki, D., Nenes, A., Fountoukis, C., Stier, P., Vignati, E., and Lelieveld, J.: Corrigendum to “Description and evaluation of GMXe: a new aerosol submodel for global simulations (v1)” published in Geosci. Model Dev., 3, 391–412, 2010, Geosci. Model Dev., 3, 413–413,, 2010b. a, b, c

Randerson, J., van der Werf, G., Giglio, L., Collatz, G., and Kasibhatla, P.: Global Fire Emissions Database, Version 3 (GFEDv3.1), ORNL DAAC,, 2013. a

Rothman, L. S., Jacquemart, D., Barbe, A., Chris Benner, D., Birk, M., Brown, L. R., Carleer, M. R., Chackerian, C., Chance, K., Coudert, L. H., Dana, V., Devi, V. M., Flaud, J.-M., Gamache, R. R., Goldman, A., Hartmann, J.-M., Jucks, K. W., Maki, A. G., Mandin, J.-Y., Massie, S. T., Orphal, J., Perrin, A., Rinsland, C. P., Smith, M. A. H., Tennyson, J., Tolchenov, R. N., Toth, R. A., Vander Auwera, J., Varanasi, P., and Wagner, G.: The HITRAN 2004 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 96, 139,, 2005.  a

Sander, R., Baumgaertner, A., Cabrera-Perez, D., Frank, F., Gromov, S., Grooß, J.-U., Harder, H., Huijnen, V., Jöckel, P., Karydis, V. A., Niemeyer, K. E., Pozzer, A., Riede, H., Schultz, M. G., Taraborrelli, D., and Tauer, S.: The community atmospheric chemistry box model CAABA/MECCA-4.0, Geosci. Model Dev., 12, 1365–1385,, 2019. a

Stier, P., Feichter, J., Kinne, S., Kloster, S., Vignati, E., Wilson, J., Ganzeveld, L., Tegen, I., Werner, M., Balkanski, Y., Schulz, M., Boucher, O., Minikin, A., and Petzold, A.: The aerosol-climate model ECHAM5-HAM, Atmos. Chem. Phys., 5, 1125–1156,, 2005. a

Tiedtke, M.: A Comprehensive Mass Flux Scheme for Cumulus Parameterization in Large-Scale Models, Mon. Weather Rev., 117, 1779–1800,<1779:ACMFSF>2.0.CO;2, 1989. a

Tost, H., Jöckel, P., Kerkweg, A., Sander, R., and Lelieveld, J.: Technical note: A new comprehensive SCAVenging submodel for global atmospheric chemistry modelling, Atmos. Chem. Phys., 6, 565–574,, 2006a. a

Tost, H., Jöckel, P., and Lelieveld, J.: Influence of different convection parameterisations in a GCM, Atmos. Chem. Phys., 6, 5475–5493,, 2006b. a

Tsimpidi, A. P., Karydis, V. A., Pozzer, A., Pandis, S. N., and Lelieveld, J.: ORACLE (v1.0): module to simulate the organic aerosol composition and evolution in the atmosphere, Geosci. Model Dev., 7, 3153–3172,, 2014 a

Tsimpidi, A. P., Karydis, V. A., Pozzer, A., Pandis, S. N., and Lelieveld, J.: ORACLE 2-D (v2.0): an efficient module to compute the volatility and oxygen content of organic aerosol with a global chemistry–climate model, Geosci. Model Dev., 11, 3369–3389,, 2018. a

Vignati, E., Wilson, J., and Stier, P.: M7: An efficient size-resolved aerosol microphysics module for large-scale aerosol transport models, J. Geophys. Res.-Atmos., 109, D22202,, 2004. a, b

Short summary
Particulate air pollution cools the climate and partially masks the greenhouse warming by reflecting sunlight and enhancing the reflection by clouds. The intensity of this cooling depends on interactions between pollution and desert dust within the atmosphere. Our simulations with a global atmospheric chemistry-climate model indicate that these interactions significantly weaken the cooling.
Final-revised paper