# Stratospheric impacts on dust transport and air pollution in West Africa and the Eastern Mediterranean

Ying Dai, Peter Hitchcock, Natalie Mahowald, Daniela Domeisen, Douglas Hamilton, Longlei Li, Beatrice Marticorena, Maria Kanakidou, Nikolaos Mihalopoulos, Adwoa Aboagye-Okyere

► To cite this version:

Ying Dai, Peter Hitchcock, Natalie Mahowald, Daniela Domeisen, Douglas Hamilton, et al. Stratospheric impacts on dust transport and air pollution in West Africa and the Eastern Mediterranean. Nature Communications, 2022, 13 (1), pp.7744. 10.1038/s41467-022-35403-1. hal-04238785

HAL Id: hal-04238785

https://hal.science/hal-04238785v1

Submitted on 13 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

# Stratospheric impacts on dust transport and air pollution in West Africa and the Eastern Mediterranean

Received: 18 May 2022

Accepted: 25 November 2022

Published online: 14 December 2022

<figure><img src="image_3.png" /></figure>

Check for updates

Ying Dai©<fnref n="1" />, Peter Hitchcock©<fnref n="1" />✉, Natalie M. Mahowald©<fnref n="1" />, Daniela I. V. Domeisen©²,³, Douglas S. Hamilton©<fnref n="1" />,⁹, Longlei Li<fnref n="1" />, Beatrice Marticorena⁴, Maria Kanakidou©<fnref n="5" />,⁶,⁷, Nikolaos Mihalopoulos<fnref n="5" />,⁸ & Adwoa Aboagye-Okyere<fnref n="1" />

Saharan dust intrusions strongly impact Atlantic and Mediterranean coastal regions. Today, most operational dust forecasts extend only 2–5 days. Here we show that on timescales of weeks to months, North African dust emission and transport are impacted by sudden stratospheric warmings (SSWs), which establish a negative North Atlantic Oscillation-like surface signal. Chemical transport models show a large-scale dipolar dust response to SSWs, with the burden in the Eastern Mediterranean enhanced up to 30% and a corresponding reduction in West Africa. Observations of inhalable particulate (PM₁₀) concentrations and aerosol optical depth confirm this dipole. On average, a single SSW causes 680–2460 additional premature deaths in the Eastern Mediterranean and prevents 1180–2040 premature deaths in West Africa from exposure to dust-source fine particulate (PM₂.₅). Currently, SSWs are predictable 1–2 weeks in advance. Altogether, the stratosphere represents an important source of subseasonal predictability for air quality over West Africa and the Eastern Mediterranean.

The world's largest dust source is North Africa<fnref n="1" />, where strong winds can lift large amounts of dust from bare, dry soils into the atmosphere. Once entrained into the atmosphere, North African dust is transported towards the Atlantic Ocean and/or towards the Mediterranean Sea²,³. Dust aerosols therefore can have vast health, visibility, environmental, and economic impacts on large population centers and industrial areas concentrated along the coasts of the Atlantic and the Mediterranean⁴⁻⁷. During episodes of Saharan dust storms that affect Europe, there is substantial evidence of the windblown desert dust's association with mortality and morbidity⁸⁻¹¹. Apart from devastating

health impacts, dust also impacts the environment, transport, and infrastructure¹²⁻¹⁴. Monetizing these impacts can translate to hundreds of million dollars just from a single dust storm¹³. According to the UN¹⁵, in the Middle East and North Africa (MENA) region, about 13 billion USD in Gross Domestic Product (GDP) are lost every year due to dust storms alone. According to the World Bank¹³, dust pollution and dust storms cost MENA over 150 billion USD annually and over 2.5% of GDP for most countries in the region. Forecasts of hazardous air pollution caused by dust aerosols are crucial to help reduce personal exposure to outdoor air pollution, and would help decision makers to restrict

<figure><img src="image_4.png" /></figure>

<footnote n="1">¹Department of Earth and Atmospheric Sciences, Cornell University, Ithaca, NY 14853, USA. ²University of Lausanne, Lausanne, Switzerland. ³ETH Zurich, Zurich, Switzerland. ⁴Laboratoire Interuniversitaire des Systèmes Atmosphériques, Universités Paris Est-Paris Diderot-Paris 7, UMR CNRS 7583, Créteil, France.</footnote>

<footnote n="5">⁵Environmental Chemical Processes Laboratory (ECPL), Department of Chemistry, University of Crete, Heraklion, Greece. ⁶Center of Studies of Air quality and Climate Change, Institute for Chemical Engineering Sciences, Foundation for Research and Technology Hellas, Patras, Greece. ⁷Excellence Chair, Institute of Environmental Physics, University of Bremen, Bremen, Germany. ⁸Institute for Environmental Research and Sustainable Development, National Observatory of Athens, Pendeli, Greece. ⁹Present address: Department of Marine, Earth, and Atmospheric Science, NC State University, Raleigh, NC, USA.</footnote>

<footnote>e-mail: aph28@cornell.edu</footnote>
anthropogenic emissions of domestic, traffic, and industrial pollutants during predicted periods of high pollution¹⁶. However, across much of Africa, dust still arrives without accurate intensity warnings¹⁷. In many parts of Europe, today's dust models are only capable of providing operational dust forecasts for several days in advance¹⁸,¹⁹. Case studies of North African dust storms with strong impact over southern Europe found a common large-scale upper-level precursor (a double Rossby wave breaking process) developing 5–10 days prior to dust storm formation. Using this precursor signal as an indicator might provide early warnings of high-impact North African dust outbreaks over southern Europe at lead times of 5–10 days²⁰,²¹. In general, the predictive lead time is limited to several days. Extending this lead time to weeks and even months would provide medical, agricultural, and nautical stakeholders more time to prepare for risks associated with airborne dust¹⁷.

One potential approach towards filling this capability-need gap is to consider the impact of the stratosphere. Stratospheric variability is well-recognized as an important source of predictability for surface weather and climate on subseasonal to seasonal (S2S) timescales, corresponding to timescales between 2 weeks and 2 months²²⁻²⁴. A particularly extreme phenomenon of the winter stratosphere that contributes to S2S predictability at the surface is the so-called sudden stratospheric warming (SSW). An SSW corresponds to a breakdown of the stratospheric polar vortex, in which the polar stratosphere warms rapidly and the westerly circumpolar winds weaken dramatically and reverse to easterlies²⁵. Currently, individual SSW events can be predicted about 1–2 weeks in advance²⁶ and can be detected early on with satellite observations²⁷. Subsequently, SSWs tend to influence surface weather for up to a few weeks²⁸. This influence allows for timescales of weeks to months to study the surface weather impacts of stratospheric forcing, pointing to their predictive value for surface weather on S2S timescales²⁹.

The tropospheric circulation response to SSW events is commonly characterized as an equatorward shift of the North Atlantic eddy-driven jet, corresponding to a negative phase of the North Atlantic Oscillation (NAO)³⁰⁻³³. The negative phase of the NAO is associated with precipitation extremes over Southern Europe³⁴ and cold air outbreaks over Northern Eurasia and the eastern United States³⁵⁻³⁷. Beyond its impact on surface meteorological fields such as temperature and precipitation, the NAO exerts a large-scale control on the dust transport out of North Africa², which makes up almost half of the dust source on Earth³⁸. Despite the reported link between the NAO and African dust export, the connection between SSWs and desert dust export has not been explored, and could contribute to long-term predictions of African dust export.

In this article, we present evidence of the impact of SSWs on African dust export and the consequent impact on air quality over S2S timescales, with chemical transport models and station observations finding up to 30% changes in aerosol concentrations over West Africa and the Eastern Mediterranean caused by SSWs. In particular, SSWs induce meteorological conditions that increase northward African dust transport to the Eastern Mediterranean, while at the same time reduce African dust emission and decrease westward dust transport to West Africa. Consequently, the former leads to enhanced dust burden in the Eastern Mediterranean and the latter results in a reduction in dust burden within West Africa. These changes caused by SSWs worsen air pollution in Southern European countries such as Greece, while they help reduce air pollution in West African countries such as Senegal.

## Results

### Detecting a large-scale dust response to sudden stratospheric warming events

As described in the introduction, SSWs can lead to a tropospheric circulation response that resembles the negative phase of the NAO. In

terms of the sea level pressure (SLP) field, a negative NAO corresponds to an anticyclonic anomaly at Stykkisholmur, Iceland and a cyclonic anomaly at Lisbon, Portugal³⁹. Composite averages over SSW episodes (defined as the 30-day periods after the onset of SSW events; see Methods) indeed reveal a cyclonic SLP anomaly over southern Europe/northern Africa (Fig. 1a, black contours and shading), similar to that associated with a negative NAO (Fig. 1a, yellow contours). On the one hand, this cyclonic SLP anomaly corresponds to a decrease in the strength of the subtropical ridge in northern Africa and the strong pressure gradient associated with it (Fig. 1b). In this region, a weakening of the south-north pressure gradient decreases the speed of northeasterly trade winds. On the other hand, the cyclonic SLP anomaly causes strong southwesterly winds from the Sahara towards the Eastern Mediterranean (Fig. 1a, black arrows). Note that the downward influence of SSWs is more intense and far reaching than its projection on to the NAO would imply. In particular, the cyclonic SLP anomaly caused by SSWs (Fig. 1a, black contours and shading) is more strongly negative and extended further eastward than the component that can be associated with the NAO (Fig. 1a, yellow contours). This NAO-SSW comparison suggests that the downward influence of SSWs is not only longer-lasting⁴⁰ but also stronger and broader than a negative NAO. Understanding their difference–potentially due to additional remote influences that affect both the stratosphere and surface weather–warrants further investigation, but is outside the scope of the present work. Given the reported link between the NAO and desert dust export², the downward influence of SSWs, which is closely linked to a negative NAO (Fig. 1a), points to a potentially important role of SSWs in controlling the North African dust export.

To detect if SSWs exert a large-scale control on North African dust export, two chemical transport models with specified meteorology are used (see Methods). Owing to their spatio-temporal continuity, the simulated surface dust concentrations enable us to obtain the large-scale dust response to SSWs. The surface dust response in CESM2 features a meridionally dipolar pattern with a positive response over the Eastern Mediterranean and a negative response over wide areas of West Africa (Fig. 1c). The east-westward aligned negative response over West Africa matches geographically well with a regional belt where surface dust concentrations are among the highest during the extended winter season from November to March (Fig. 1d). The MERRA2-surface dust concentrations show a very similar dipolar response (Supplementary Fig. 2a), demonstrating the robustness of the large-scale surface dust response to SSWs.

To further quantify the dust response to SSWs, SSW-caused changes in dust pollution are obtained by comparing the dust pollution during SSW episodes with those during non-SSW episodes. For this purpose, 1000 sets of non-SSW episodes are identified (see Methods). The dust pollution is measured by the level of surface dust concentrations and the number of high dust pollution days. These days correspond to those on which the specified daily thresholds 250, 500, and 1000 µg m⁻³ are exceeded (Supplementary Table 1). In the Eastern Mediterranean, SSWs lead to a 20–30% increase in the level of surface dust concentrations and a 20–40% increase in the number of high dust pollution days (Supplementary Table 1). In West Africa, decreases are simulated. The level of surface dust concentrations is reduced by 10–20%, and the number of high dust pollution days decreases between 10 and 30%, depending on the model and threshold used (Supplementary Table 1). All these changes are relative to the non-SSW mean, which is defined as an average over 1000 sets of non-SSW episodes (see Methods). The quantitative analysis shows that dust pollution at the surface is influenced to a major extent by SSWs.

The dust represents one of the main natural contributions to atmospheric particulate matter (PM), which has direct effects on people's health when inhaled. We next present results of how SSWs, via influencing dust pollution at the surface, impact mortality from exposure to dust-source fine particulate matter less than 2.5 µm in
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

**Fig. 1** | Large-scale meteorological and dust response to sudden stratospheric warming events (SSWs). **a** Composites of anomalous sea level pressure (SLP) (contours) and 10 m wind (arrows) during SSW episodes over the 1980/81-2013/14 extended winters (from November to March, NDJFM). The SLP and 10 m wind fields come from the MERRA2 reanalysis. Shadings indicate SLP anomalies that are statistically significant at the *p* < 0.10 level based on a two-tailed Monte Carlo test. Wind vectors shown are statistically significant at least for one component. Scaling for the wind vectors is given in the bottom-right corner (units: ms⁻¹). Yellow contours indicate the regressed SLP anomalies onto the extended wintertime (NDJFM) North Atlantic Oscillation (NAO) index (see Methods). To make them comparable

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

to the composite SLP anomalies following SSWs, the regressed SLP anomalies are scaled by the average value of the NAO index during SSW episodes, which is slightly greater than -0.3 (a histogram is shown in Supplementary Fig. 1). **b** The corresponding extended wintertime (NDJFM) climatologies of SLP and 10 m wind fields. **c, d** The same as **a**, **b** but for CESM2 surface dust concentrations (units: µg m⁻³). The two red boxes in panel **c** indicate the Eastern Mediterranean (30°–33.75°N, 10°–30°E) and West Africa (16.25°–23.75°N, 12.5°–23.75°E) regions with large and significant dust response to SSWs. Made with Natural Earth. Free vector and raster map data @ naturalearthdata.com.

diameter (PM₂.₅) (see Methods). Generally, SSWs lead to additional premature deaths in the Eastern Mediterranean via enhancing PM₂.₅ concentrations and reduced premature deaths in West Africa by decreasing PM₂.₅ concentrations (Supplementary Fig. 3). On average, a single SSW can cause 2460 (±930, 1σ) premature deaths in the Eastern Mediterranean (the red box in Supplementary Fig. 3a) and prevent 2040 (±660, 1σ) premature deaths in West Africa (the blue box in Supplementary Fig. 3a), as estimated from the CESM2 output of dust-source PM₂.₅. MERRA2 suggests a smaller but comparable signal, that is, a single SSW can lead to 680 (±210, 1σ) additional premature deaths in the Eastern Mediterranean and 1180 (±750, 1σ) reduced premature deaths in West Africa. The difference in the magnitude of the response in different models is consistent with the uncertainty in the dust cycle across different dust models⁴¹. The impact of SSWs on dust pollution mortality suggests that skillful forecasts of dust pollution during SSWs may alleviate pressure on health and social care systems.

## Meteorological causes of the dipolar dust response

Understanding the meteorological causes of the dipolar dust response is important because it provides theoretical support for the model results. As we show next, during SSW episodes, the enhanced dust pollution over the Eastern Mediterranean is linked to increased northward African dust transport, while the reduction in dust pollution within West Africa arises from a combination of

suppressed African dust emission and decreased westward dust transport.

In the Eastern Mediterranean, the southwesterly winds, which facilitate northward Sahara desert dust transport, create a potential for high dust concentrations (Fig. 2a). For example, the maximum particulate matter less than 10 µm in diameter (PM₁₀) concentrations level for 2015 recorded at Finokalia (Crete, southern Greece) was related to a desert dust outbreak event caused by a strong southwesterly airflow⁴². During SSW episodes, for the Eastern Mediterranean, a significant increase in the fraction of time and area that the winds are southwesterly is seen (Fig. 2b), consistent with the fact that the SSW-induced cyclonic SLP anomaly is associated with southwesterly winds from the Sahara towards the Eastern Mediterranean (Fig. 1a, black arrows). In this sense, the enhanced dust burden over the Eastern Mediterranean during SSW episodes (Fig. 1c) is likely related to the increased northward Saharan dust transport induced by an increase in the fraction of southwesterly winds (Fig. 2b). In terms of the surface wind speeds, however, there are no significant changes due to the SSW (Supplementary Fig. 4a). Note that the Eastern Mediterranean is also exposed to dust originating from western Asia (mainly the Arabian Peninsula). Therefore, Arabian Desert dust imported under an easterly flow can also lead to dust intrusions into the Eastern Mediterranean⁴³⁻⁴⁵. However, during SSW episodes, a significant decrease in the fraction of easterly wind is
seen (Fig. 2b), indicating that the enhanced dust burden over the Eastern Mediterranean during SSW episodes is unlikely contributed by dust from the Arabian Desert. There is no evidence for changes in dust burden due to the SSW in the Western Mediterranean basin (Fig. 1c).

In West Africa, a decrease in surface wind velocity is associated with a lower probability of high dust concentrations (Fig. 2c), as the slower winds pick up and transport less dust from the Sahara to West Africa⁴⁶. In particular, extremely intense PM₁₀ concentrations at Cinzana (Mali) are recorded only for higher than 5 ms⁻¹ local wind velocities⁴⁷. During SSW episodes, for West Africa, a weakening of surface wind velocities is seen (Fig. 2d), as can be inferred from the

decreased probability of the high wind velocities (higher than 5 ms⁻¹) and the increased probability of the low wind velocities (lower than 4 ms⁻¹) (Fig. 2d). In this case, the reduced dust burden within West Africa during SSW episodes (Fig. 1c) may be linked to the suppressed dust emission and reduced westward dust transport in response to a reduction in surface wind velocities (Fig. 2d). In terms of the surface wind directions, however, there are no significant changes due to the SSW (Supplementary Fig. 4b).

To sum up, during the extended boreal winter season (NDJFM) (Fig. 3a), the subtropical ridge tends to be located over northern Africa (Fig. 3a, red contour line; Supplementary Fig. 5a, thick contour line), and is associated with strong northeasterly trade winds (Fig. 3a, black

<figure><img src="image_9.png" /></figure>

**Fig. 2 |** Observed changes in surface wind fields within the Eastern Mediterranean and West Africa. The Eastern Mediterranean and West Africa correspond to the two red boxes shown in Fig. 1c. **a** Bivariate histogram of surface dust concentrations (from the CESM2 simulation) and surface wind directions (from the MERRA2 reanalysis) at all grid points within the Eastern Mediterranean during 1980/81-2013/14 extended winters (from November to March, NDJFM), with colors indicating the probability of occurrences falls within each bin. For surface wind directions, "NE", "SE", "SW", and "NW" indicate north-easterlies, south-easterlies, south-westerlies, and north-westerlies, respectively. **b** Histogram of surface wind directions at all grid points within the Eastern Mediterranean during SSW episodes (gray filled bars) and during 1000 sets of non-SSW episodes (see Methods). For

1000 sets of non-SSW episodes, black unfilled bars show the mean and black whiskers show the spread. The boxes (filled in black, red, and blue colors) indicate the 5th and 95th percentiles. The boxes are filled in red and blue when SSW-caused increases and decreases are statistically significant at the *p* < 0.10 level based on a two-tailed Monte Carlo test, respectively. **c, d** The same as **a, b** but for surface wind velocities (from the MERRA2 reanalysis) at all grid points within West Africa. To generate the bivariate histogram plots, we partition the surface dust concentrations into 500 bins of width 20 µg m⁻³ within [0, 10,000], the surface wind direction into 720 bins of width 0.5° within [0, 360], and the surface wind velocity into 700 bins of width 0.02 ms⁻¹ within [0, 14], respectively.

<figure><img src="image_10.png" /><figcaption>**a** Wintertime Climatology</figcaption></figure>

**Fig. 3 |** Schematic presentation of the stratospheric impacts on tropospheric circulation and North African dust. **a** The extended wintertime (from November to March, NDJFM) climatology of the tropospheric circulation and North African dust. **b** The same as panel a but for the sudden stratospheric warming (SSW) episodes. Various quantities used to illustrate the dust response include the sub-tropical ridge (red contour line; see Supplementary Fig. 5 for a description), surface winds (black and blue arrows), dust emission (brown cloud-shaped pattern), and dust transport (brown arrows). As can be seen, during the extended winter season (NDJFM) (panel a), the subtropical ridge tends to be located over northern Africa.

<figure><img src="image_11.png" /><figcaption>**b** Days [1,30] after SSWs</figcaption></figure>

The associated northeasterly trade winds pick up and transport much of the dust aerosols southwestward. Meanwhile, the Eastern Mediterranean is under the influence of midlatitude westerlies. During SSW episodes (panel b), the subtropical ridge retreats to the North Atlantic Ocean. Weaker northeasterly trade winds pick up and transport less dust aerosols southwestward, while southwesterly winds are more likely to occur over the Eastern Mediterranean, increasing northward dust transport. Made with Natural Earth. Free vector and raster map data @ naturalearthdata.com.
<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

**Fig. 4 | Observed air quality response to sudden stratospheric warming events (SSWs).** a Regression of anomalous sea level pressure (SLP) (from the MERRA2 reanalysis) onto the station measurements of daily time-series of surface inhalable particulate (PM₁₀) concentrations from Finokalia, Greece. Warm and cold shadings indicate positive and negative anomalies that are statistically significant at the *p* < 0.10 level based on a Student's *t*-test, respectively. Histogram of **b** the level of PM₁₀ concentrations (units: µg m⁻³) at Finokalia, Greece and **c**, **d** the expected number of poor air quality days per 30-day period for two air quality standards

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

(AQS), during 1000 sets of non-SSW episodes (colored bars, see Methods). Thick vertical colored line indicates the non-SSW mean; gray shading indicates the 5th–95th percentile confidence interval derived from the non-SSW spread. Thick vertical black line indicates the SSW composite. In panels **c**, **d**, poor air quality days are defined when the daily-mean PM₁₀ concentrations exceed the daily limit of **c** 50 µg m⁻³ and **d** 260 µg m⁻³, respectively. **e–h** The same as **a–d** but for station observations of surface PM₁₀ concentrations from M'Bour, Senegal. Made with Natural Earth. Free vector and raster map data @ naturlearthdata.com.

arrows). In arid regions such as North Africa, strong trade winds favor dust emissions to the atmosphere (Fig. 3a, brown cloud-shaped pattern) and help transport much of the dust aerosols southwestward (Fig. 3a, brown arrow). Meanwhile, the Eastern Mediterranean is under the influence of midlatitude westerlies (Fig. 3a, blue arrow). During SSW episodes, a cyclonic SLP anomaly is observed in the region of southern Europe/northern Africa, which is a typical pressure center of action of the negative NAO phase (Fig. 1a). Over the Eastern Mediterranean, the cyclonic SLP anomaly brings southwesterly winds (Fig. 3b, blue arrow), which favor northward dust transport (Fig. 3b, brown arrow, pointing northeast) and thus enhance dust burden in the Eastern Mediterranean (Fig. 1c). Over the continent to the south, in the presence of the cyclonic SLP anomaly, the subtropical ridge retreats to the North Atlantic Ocean (Fig. 3b, red contour line; Supplementary Fig. 5b, thick contour line). This corresponds to a weakening of the northeasterly trade winds (Fig. 3b, black arrows), which results in a decrease in African dust emission (Fig. 3b, brown cloud-shaped pattern) and a reduction in westward dust transport (Fig. 3b, brown arrow,

pointing southwest). As a result, reduced dust burden is seen in West Africa (Fig. 1c).

## Observational evidence of the dipolar dust response

Changes in dust transport and emission will shape changes in human environment. Of particular importance will be those changes in air quality impacting the local population of the affected regions. Here we use a range of observational datasets to quantify the impact of SSWs on the air quality of countries in Southern Europe and West Africa.

First, we use the daily-mean PM₁₀ concentrations collected at background stations in Greece and Senegal (see Methods). Overall, the regressed SLP anomaly onto the PM₁₀ time-series at the two stations suggests that, higher-than-normal PM₁₀ concentrations in Greece tend to occur under the influence of a cyclonic SLP anomaly (Fig. 4a), while higher-than-normal PM₁₀ concentrations in Senegal tend to occur in the presence of an anticyclonic SLP anomaly (Fig. 4e). Given the fact that a cyclonic SLP anomaly is seen following SSWs (Fig. 1a), above-normal PM₁₀ concentrations (worse air quality) in Greece, and below-
normal PM₁₀ concentrations (better air quality) in Senegal can be expected during SSW episodes.

To examine whether this expectation holds, two air quality indicators are considered: (a) the level of PM₁₀ concentrations, and (b) the expected number of poor air quality days when PM₁₀ concentrations exceed the given thresholds (two thresholds are considered, see Methods). In Greece, all air quality indicators during SSW episodes (black vertical lines in Fig. 4b-d) lie entirely beyond the upper end of the non-SSW spread (color histograms in Fig. 4b-d), indicating that the air quality in Greece is worsened significantly during SSW episodes. In particular, the level of PM₁₀ concentrations reaches 31 µg m⁻³ (black vertical line in Fig. 4b), an increase of 64% relative to the non-SSW mean (19 µg m⁻³; purple vertical line in Fig. 4b); the expected number of poor air quality days per 30-day period exceeding 50 µg m⁻³ reaches 3.8 days (black vertical line in Fig. 4c), roughly four times the non-SSW mean (one day; green vertical line in Fig. 4c); the expected number of poor air quality days per 30-day period exceeding the higher threshold of 260 µg m⁻³ also increases by a factor of four (Fig. 4d), though PM₁₀ concentrations at this station rarely reach this level. In Senegal, the level of PM₁₀ concentrations and the expected number of poor air quality days exceeding the threshold of 260 µg m⁻³ decrease significantly during SSW episodes (Fig. 4f, h), indicating that the air quality in Senegal is improved considerably due to the SSW. In particular, the level of PM₁₀ concentrations is 136 µg m⁻³ (black vertical line in Fig. 4f), a decrease of 34% relative to the non-SSW mean (207 µg m⁻³; purple vertical lines in Fig. 4f); the expected number of poor air quality days exceeding the threshold of 260 µg m⁻³ decreases to 2.4 days (black vertical line in Fig. 4h), about 34% of the non-SSW mean (7.0 days; green vertical line in Fig. 4h). The lower threshold of 50 µg m⁻³ is nearly always exceeded at this station, and no significant changes are seen during SSW episodes in this case (Fig. 4g). We note nonetheless that changes in the expected number of poor air quality days during SSW episodes are robust across a wide range of thresholds (Supplementary Fig. 6), and that exposure to air pollutants at any level has observable health impacts⁴⁸⁻⁵⁰.

To sum up, SSWs lead to a significant increase in the level of PM₁₀ concentrations in Greece and a significant decrease in the level of PM₁₀ concentrations in Senegal. The oppositely signed air quality changes in response to SSWs (worsened in Greece and improved in Senegal) detected in station observations of PM₁₀ concentrations show good agreement with the dipolar dust response derived from the modeled datasets (Fig. 1c; Supplementary Fig. 2a).

Using the ground-based Aerosol Robotic Network (AERONET) direct measurements of aerosol optical depth (AOD) yields responses similar to the above results, which consist of above-normal AOD (increased by 10%) in Greece and below-normal AOD (reduced by 11%) in Senegal (Supplementary Table 2). The consistency of model results with station observations lends strong support to the impacts of SSWs on air pollution in Southern Europe and West Africa.

## Discussion

Saharan dust intrusions are a recurrent problem threatening public health and economic well-being of a large population in Southern Europe and West Africa. Forecasting such high impact but short duration events poses profound scientific challenges⁵¹. Major progress in the basic understanding and predictive capability can take decades¹⁷. Today, many developing countries across Africa do not have the computing power to forecast air quality¹⁷. In many parts of Europe, most forecasts of air quality only cover 2–5 days. Extending this to longer term forecasts would give communities more time to prepare for risks associated with air pollution.

Using a wide range of datasets, our results demonstrate the existence of large-scale stratospheric impacts on North African dust emission and transport on S2S timescales with the consequent impact on air quality, implying the potential for long-term dust and air

pollution predictions. In particular, chemical transport models consistently suggest a large-scale dipole dust response to SSWs, which consists of an enhanced dust signal in the Eastern Mediterranean and a reduction in West Africa. The dust pollution changes caused by SSWs can further impact premature deaths from exposure to PM₂.₅. According to our estimate, a single SSW, on average, can cause 680–2460 premature deaths in the Eastern Mediterranean via dust-source PM₂.₅ enhancement and prevent 1180–2040 premature deaths in West Africa through dust-source PM₂.₅ reduction. The changes from an SSW event in dust pollution mortality are comparable to the changes in cold weather mortality in the UK caused by an SSW, which is 620 additional deaths per event²². Using station observations of PM₁₀ concentrations and AERONET AOD yields very similar air pollution responses to SSWs, including worsened air quality in Greece and improved air quality in Senegal, in line with the dipolar dust response derived from models. The stratospheric impacts on dust transport and air pollution at the surface are driven by the role of SSWs in establishing a negative NAO-like signal at the surface. This NAO phase is typically associated with a weakened subtropical ridge, which, through increasing southwesterly winds in the Eastern Mediterranean and weakening northeasterly trade winds in West Africa, creates large-scale meteorological conditions favorable for the dipolar dust response. The agreement between models and the consistency with station observations and the meteorological causes lend strong support to the stratospheric impacts on air pollution at the surface. Overall, the changes in air pollution in response to SSWs range between 10 and 30% (with a couple of outliers reaching 60%) compared to the non-SSW mean, depending on the region of concern and on the dataset.

The oppositely signed air quality changes between the Eastern Mediterranean and West Africa in response to SSWs suggest that the occurrence of SSWs is not always a tale of destruction. While SSWs act to worsen air pollution in Southern European countries such as Greece, for West African countries such as Senegal, air quality improves in response to SSWs although PM₁₀ concentrations remain at very high levels. In addition to weak vortex events such as SSWs, there are also strong vortex events that correspond to tropospheric weather regimes of opposite sign to those following SSWs²⁸. Compared to SSWs, the strong vortex events might have the opposite effect on the African dust export and strongly worsen air quality in West African countries instead. Furthermore, so-called final warming events⁵³ in spring can lead to same-signed surface anomalies as SSW events⁵⁴, and may therefore extend the surface impacts of the stratosphere into spring, with similar effects on dust transport as SSW events. We note that stratospheric extreme events do not necessarily cause dust storms, rather they act over a period of weeks to months to increase the probability of meteorological conditions favorable for high dust pollution in certain regions. This knowledge helps assess the likelihood that high dust pollution may occur and thus can provide an advance warning for adverse conditions.

Currently, extreme stratospheric polar vortex events are predictable at 1–2 weeks lead times in operational S2S prediction systems²⁶. Furthermore, once an extreme stratospheric event occurs, it can lead to changes at surface air pollution that can last for weeks to months. In this context, the stratosphere represents a non-negligible source of predictability for surface air pollution on S2S timescales, particularly for Southern Europe and West Africa. These findings provide a window of opportunity for advanced air quality warnings at lead times of weeks to months based on the forecast of extreme stratospheric events, and call for consideration of the stratospheric variability in air quality forecasting systems.

## Methods

### Datasets

The meteorological fields examined include daily-mean sea level pressure (SLP), 10 m zonal wind (U10 m), and 10 m meridional wind
(V10 m) from MERRA2 (the Modern-Era Retrospective Analysis for Research and Applications, version 2)⁵⁵. These meteorological fields have a 1.25° × 1.25° latitude-longitude resolution.

For air quality related fields, we use daily data from a wide range of sources, including model outputs, station observations, and aerosol reanalysis. Firstly, we examine simulated daily-mean surface dust concentrations from CESM2 (Community Earth System Model, version 2)⁵⁶ and MERRA2 (The MERRA2 Aerosol Reanalysis)⁵⁷. Note that the dust field from MERRA2 is not directly constrained by the aerosol assimilation included in MERRA2 although the overall optical depth is constrained⁵⁷. The dust field from MERRA2 is thus to a large extent a model-based product, just like that from CESM2⁵⁶. Dust in CESM2 and MERRA2 simulations is modeled via the chemical transport models DEAD (Dust Entrainment And Deposition model)⁵⁸ and GOCART (Goddard Chemistry, Aerosol, Radiation, and Transport model)⁵⁹, respectively. Both CESM2 and MERRA2 simulations are driven by MERRA2 meteorology. The CESM2 output covers the years 1980–2014 with a latitude-longitude resolution at -0.9° × 1.25°. The MERRA2 output covers the years 1980–2020 with a latitude-longitude resolution at 1.25° × 1.25°. Owing to their spatio-temporal continuity, the model outputs enable us to obtain the large-scale dust response to stratospheric extreme events.

Secondly, we present a verification of the findings derived from the model output with station observations of surface PM₁₀ concentrations and aerosol optical depth (AOD) (see Supplementary Fig. 7 for the station location and Supplementary Table 3 for the temporal coverage). The PM₁₀ concentrations are measured at four stations, including three sites in West Africa (M'Bour, Senegal; Cinzana, Mali; Banizoumbou, Niger) from INDAAF (International Network to study Deposition and Atmospheric composition in Africa) and one site at Finokalia (southern Greece) from Finokalia Atmospheric Observatory. The station AOD at the aforementioned four stations comes from AERONET (Aerosol Robotic Network)⁶⁰. Unlike the model outputs, the observational data are not recorded everyday due to network and/or power failures, instrument maintenance, etc. In this study, the days with missing values are discarded. This is unlikely to affect the representativeness of the measurements given their high recovery rate⁴⁷. Note that the potential impacts of sampling error from short observational records and limited background stations may represent a potential source of uncertainty in observed air quality responses (Supplementary Note 1).

Thirdly, the gridded AOD from the MERRA2 aerosol reanalysis⁵⁷ is analyzed. Unlike the MERRA2 dust mentioned above, the MERRA2 AOD is directly constrained by the aerosol assimilation and therefore is a blend of observations with model hindcasts⁵⁷. To enable evaluation against the station AOD from AERONET, the MERRA2 AOD time-series at the grid point closest to each AERONET station is evaluated.

For the gridded datasets, daily anomalies at each grid point are obtained by subtracting the seasonal cycle on that calendar day. The seasonal cycle is defined as the mean and first three Fourier harmonics of the daily climatology to remove the high-frequency sampling variability, following previous studies⁶¹⁻⁶³.

## Sudden stratospheric warnings

Sudden stratospheric warming events (SSWs) are defined by the reversal of the zonal-mean westerlies at 60°N and 10 hPa¹⁰. The first day on which the daily-mean zonal-mean zonal wind at 60°N and 10 hPa is easterly, is defined as the central date of the warming. Using the daily-mean zonal-mean zonal wind at 60°N and 10 hPa from the MERRA2 reanalysis, 24 SSW events are identified in the 40 extended winters (November–March, NDJFM) during 1980/81-2019/20 (Supplementary Table 4) where the MERRA2 dust is available. Of these, 22 SSWs occur during the 1980/81-2013/14 extended winters (NDJFM) when the CESM2 dust is available. The definition of SSWs is insensitive to the reanalysis dataset because the central dates for SSWs in each reanalysis

are in good agreement for the analyzed time period⁶⁴. The impact of SSWs on the near-surface circulation is often described as resembling a negative phase of the North Atlantic Oscillation (NAO)³⁰⁻³³. In this study, the NAO index is defined from the difference between normalized SLP between Lisbon, Portugal, and Stykkisholmur, Iceland³⁹. This definition has been used in a previous study to demonstrate the role of the NAO in controlling the North African dust export².

## SSW episodes and non-SSW episodes

To examine the tropospheric response to SSWs, we define “SSW episodes” as the 30-day periods after the onset of SSW events, because this is when the downward influence of SSWs on the troposphere can be most clearly identified²⁸,³¹. For comparison, we define “non-SSW episodes” as the 30-day periods after the onset of non-SSW events. Each set of non-SSW events have the equivalent central dates as the observed SSWs but are randomly taken from extended winters (NDJFM) without an SSW. To account for the uncertainty in random sampling, 1000 sets of non-SSW events are generated. As reference, the mean of 1000 sets of non-SSW episodes (hereafter referred to as “the non-SSW mean”) is used to detect by how much the air pollution during SSW episodes differs from that during non-SSW episodes, and the spread of the 1000 sets of non-SSW episodes (hereafter referred to as “the non-SSW spread”) is used to determine to what extent the changes caused by SSWs are statistically significant.

## Premature mortality estimation

SSW events are associated with changes in surface dust concentrations, which impact air quality and human health. After the occurrence of SSWs, changes in mortality from exposure to PM₂.₅ are estimated as⁶⁵: $ΔMort = M_b × N_{pop} × [\exp(βΔX) - 1]$, where M_b is the baseline mortality obtained from the Global Burden of Disease (GBD) study, N pop is the exposed population obtained from the Center for International Earth Science Information Network, β is the concentration-response factor for mortality due to short-term exposure to PM₂.₅ and is based on a 2-day averaged anomaly in PM₂.₅ concentrations (a 10 µg m⁻³ increase in 2-day averaged PM₂.₅ was found to be associated with a 0.98% increase in total mortality⁶⁶), ΔX is the 2-day averaged anomaly in PM₂.₅ concentrations over days [1, 30] following SSW events (22 SSWs for CESM2 and 24 SSWs for MERRA2). Here, the model output of surface PM₂.₅ concentrations specific to dust aerosols is used. In this way, the estimated mortality changes are attributable to SSWs given their impact on surface dust aerosols. There are uncertainties in the β deduced as the concentration-response factor, as well as what baseline mortality to use in the exposure by gender, age or comorbidities, but here we use a standard relationship and standard mortality values used for air pollution impacts⁶⁶.

Mortality changes are calculated at each grid point for each of the 15, 2-day periods over days [1, 30] following the occurrence of SSWs. The 15 here indicates the number of 2-day periods over days [1, 30] after the onset of SSWs. For example, the 1st 2-day period corresponds to days [1, 2], the 2nd 2-day period corresponds to days [3, 4], the 3rd 2-day period corresponds to days [5, 6], and the 15th 2-day period corresponds to days [29, 30]. For each SSW event, the mortality changes at each of the 15, 2-day periods are added together to obtain an estimate of the excess deaths that occurred during the SSW event. The excess deaths for each SSW event are further used to calculate the average excess deaths per event (Supplementary Fig. 3a, c). Similar to the dipolar dust-source PM₂.₅ response pattern (Supplementary Fig. 3b, d), the mortality changes also exhibit a dipolar pattern but with the centers moving towards large population centers (Supplementary Fig. 3a, c). The mortality changes estimated from CESM2 (Supplementary Fig. 3a) and MERRA2 (Supplementary Fig. 3c) are generally comparable although CESM2 suggests larger changes than MERRA2. In fact, across the models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6), there is a wide range in dust
emissions wherein CESM2 is at the higher end while MERRA2 is at the lower end⁴¹. In this sense, our results obtained from CESM2 and MERRA2 represent the diversity in dust emission and therefore provide a robust evaluation of the possible health impacts of SSWs.

For each of the designated regions (the red and blue boxes in Supplementary Fig. 3a, c), the mortality changes at each grid point within the region and each of the 15, 2-day periods for each SSW event are added together to get an estimate of the total excess deaths that occurred during each SSW event. This gives 22 individual numbers for CESM2, one per SSW event (or 24 individual numbers for MERRA2). These 22 numbers obtained from CESM2 (or 24 numbers from MERRA2) are further used to calculate the average excess deaths per event, as well as a standard deviation across events.

## Air quality indicators

While the smaller particles (PM<sub>2.5</sub>) have been used in the previous section, because observations of these are not available we focus on observational records of PM<sub>10</sub> in this section as an indicator of air quality because long-term observations of PM<sub>10</sub> are available for the studied areas. Using the daily-mean surface PM<sub>10</sub> concentrations, two air quality indicators are defined: (a) the level of PM<sub>10</sub> concentrations, and (b) the expected number of poor air quality days per 30-day period. Poor air quality days are defined when the daily-mean PM<sub>10</sub> concentrations exceed the given air quality standard. In this study, for each station, two air quality standards are used. These include the European Union air quality standard (daily-mean PM<sub>10</sub> concentrations of 50 µg m<sup>-3</sup>) and a higher threshold (daily-mean PM<sub>10</sub> concentrations of 260 µg m<sup>-3</sup>). The latter is the National Air Quality Standard for Senegal⁶⁷; it has also been used in a previous study as the air quality standard for Tunisia, which is located close to desert dust source⁶⁸. Note that exposure to permissible concentrations of air pollutants still has observable health impacts⁴⁸⁻⁵⁰, indicating that the air quality standards should be regularly reviewed and revised as new scientific evidence emerges on adverse effects on public health and the environment.

Here is an example of how we compute the two air quality indicators. At Finokalia (Greece), the daily-mean PM<sub>10</sub> concentrations are collected between years 2004–2019, wherein 9 SSW events took place (Supplementary Table 3). Regarding the first air quality indicator, the level of PM<sub>10</sub> concentrations during SSW episodes is computed as the average of daily-mean PM<sub>10</sub> concentrations over the 270 days following the 9 SSW events (30 days per SSW event). The level of PM<sub>10</sub> concentrations during each set of non-SSW episodes is computed in the same way. To obtain the second air quality indicator, we first count the number of poor air quality days during SSW episodes and during each set of non-SSW episodes. However, a direct comparison of the number of poor air quality days between SSW and non-SSW episodes would not be fair if the number of days with available data differs between SSW and non-SSW episodes, given the fact that the observational data are not recorded everyday due to network and/or power failures, instruments maintenance, etc. To account for this uncertainty, the number of poor air quality days is divided by the number of available data days to obtain the rate of poor air quality days. To facilitate understanding, the rate of poor air quality days is multiplied by 30 to obtain the expected number of poor air quality days per 30-day period. Note that the expected number of poor air quality days obtained in this way is not necessarily integer.

To assess whether SSWs are able to cause changes in the air quality, the air quality indicators during SSW episodes are compared to the distribution of those during non-SSW episodes. The differences, if any, are expected to arise from the downward influence of SSWs.

## Statistical significance assessment

Statistical significance of the regression analysis used in this study is evaluated using a two-tailed Student's t-test, where the effective degrees of freedom are calculated following equation (1) of ref. 69.

Statistical significance of the composite analysis based on SSWs is assessed with a two-tailed Monte Carlo test. For anomalous fields (the seasonal cycle is subtracted), one thousand sets of composites are computed around *N* randomly chosen 'central dates' from the pool of days during all extended winters (NDJFM) with available data, where *N* is equal to the number of SSWs involved in the composite analysis. For full fields (the seasonal cycle is not subtracted), one thousand sets of composites based on non-SSW events are computed, so as to account for the influence of seasonal cycle in the full fields. The composite fields based on SSWs are considered significant (at the *p* < 0.10 level by a two-sided test) if it lies outside of the 5th–95th percentile range derived from the randomly generated composites.

## Data availability

The MERRA2 data used in this study come from https://gmao.gsfc.nasa.gov/reanalysis/MERRA-2/. The CESM2 data are available from the Cornell eCommons repository (https://doi.org/10.7298/xqqj-qk90). The PM<sub>10</sub> concentrations measured at three stations in West Africa (M'Bour, Senegal; Cinzana, Mali; Banizoumbou, Niger) come from INDAAF. The PM<sub>10</sub> concentrations measured at Finokalia (southern Greece) come from Finokalia Atmospheric Observatory (https://finokalia.chemistry.uoc.gr/). The station AOD at the aforementioned four stations comes from AERONET. The baseline mortality used in this study is obtained from the Global Burden of Disease (GBD) study (http://ghdx.healthdata.org/gbd-results-tool), the exposed population is obtained from the Center for International Earth Science Information Network (https://sedac.ciesin.columbia.edu/data/set/gpw-v4-population-count-rev11).

## Code availability

The CESM2 model is freely available to the scientific community. The code can be accessed through a public GitHub repository (https://github.com/ESCOMP/CESM).

## References

1. Goudie, A. S. & Middleton, N. J. Saharan dust storms: nature and consequences. *Earth Sci. Rev.* **56**, 179–204 (2001).

2. Moulin, C., Lambert, C. E., Dulac, F. & Dayan, U. Control of atmospheric export of dust from North Africa by the North Atlantic Oscillation. *Nature* **387**, 691–694 (1997).

3. Tsamalis, C., Chédin, A., Pelon, J. & Capelle, V. The seasonal vertical distribution of the Saharan Air Layer and its modulation by the wind. *Atmos. Chem. Phys.* **13**, 11235–11257 (2013).

4. Alastuey, A. et al. Characterisation of TSP and PM2.5 at Izaña and Sta. Cruz de Tenerife (Canary Islands, Spain) during a Saharan Dust Episode (July 2002). *Atmos. Environ.* **39**, 4715–4728 (2005).

5. Kishcha, P. et al. Vertical distribution of Saharan dust over Rome (Italy): Comparison between 3-year model predictions and lidar soundings. *J. Geophys. Res. Atmos.* **110**, 6 (2005).

6. Perez, C. et al. A long Saharan dust event over the western Mediterranean: lidar, Sun photometer observations, and regional dust modeling. *J. Geophys. Res. Atmos.* **111**, 15 (2006).

7. Escudero, M. et al. A methodology for the quantification of the net African dust load in air quality monitoring networks. *Atmos. Environ.* **41**, 5516–5524 (2007).

8. Jiménez, E., Linares, C., Martínez, D. & Díaz, J. Role of Saharan dust in the relationship between particulate matter and short-term daily mortality among the elderly in Madrid (Spain). *Sci. Total Environ.* **408**, 5729–5736 (2010).

9. Mallone, S. et al. Saharan dust and associations between particulate matter and daily mortality in Rome, Italy. *Environ. Health Perspect.* **119**, 1409–1414 (2011).

10. Samoli, E., Kougea, E., Kassomenos, P., Analitis, A. & Katsouyanni, K. Does the presence of desert dust modify the effect of PM₁₀ on
mortality in Athens, Greece? *Sci. Total Environ.* **409**,
2049–2054 (2011).

11. Karanasiou, A. et al. Health effects from Sahara dust episodes in Europe: literature review and research gaps. *Environ. Int.* **47**, 107–114 (2012).

12. Goudie, A. S. & Middleton, N. J. *Desert Dust in the Global System* (Springer-Verlag, 2006)

13. World Bank. *Sand and Dust Storms in the Middle East and North Africa Region-Sources, Costs, and Solutions*. https://www.worldbank.org/en/region/mena/publication/sand-and-dust-storms-in-the-middle-east-and-north-africa-mena-region-sources-costs-and-solutions. (2019)

14. Monteiro, A. et al. Multi-sectoral impact assessment of an extreme African dust episode in the Eastern Mediterranean in March 2018. *Sci. Total Environ.* **843**, 156861 (2022).

15. United Nations Environment Assembly. UNEA-2 FACT SHEET: SAND AND DUST STORMS. https://www.unep.org/resources/factsheet/sand-and-dust-storms. (2016).

16. Baklanov, A. & Zhang, Y. Advances in air quality modeling and forecasting. *Glob. Transitions* **2**, 261–270 (2020).

17. Kumar, R., Peuch, V.-H., Crawford, J. & Brasseur, G. Five steps to improve air-quality forecasts. *Nature* **561**, 27–29 (2018).

18. Pérez, C. et al. Atmospheric dust modeling from meso to global scales with the online NMMB/BSC-Dust model - Part 1: Model description, annual simulations and evaluation. *Atmos. Chem. Phys.* **11**, 13001–13027 (2011).

19. Kukkonen, J. et al. A review of operational, regional-scale, chemical weather forecasting models in Europe. *Atmos. Chem. Phys.* **12**, 1–87 (2012).

20. Dhital, S., Kaplan, M. L., Orza, J. A. G. & Fiedler, S. Atmospheric Dynamics of a Saharan Dust Outbreak Over Mindelo, Cape Verde Islands, Preceded by Rossby Wave Breaking: Multiscale Observational Analyses and Simulations. *J. Geophys. Res. Atmos.* **125**, 2020–032975 (2020).

21. Orza, J. A. G., Dhital, S., Fiedler, S. & Kaplan, M. L. Large scale upper-level precursors for dust storm formation over North Africa and poleward transport to the Iberian Peninsula. Part I: An observational analysis. *Atmos. Environ.* **237**, 117688 (2020).

22. Butler, A. et al. in *Sub-Seasonal to Seasonal Prediction Robertson*, (eds. A.W., Vitart, F. eds.) Ch. 11, 223–241. (Elsevier, 2019).

23. Domeisen, D. I. V. & Butler, A. H. Stratospheric drivers of extreme events at the Earth's surface. *Commun. Earth Environ.* **1**, 59 (2020).

24. Domeisen, D. I. V. et al. The role of the stratosphere in subseasonal to seasonal prediction: 2. predictability arising from stratosphere-troposphere coupling. *J. Geophys. Res. Atmos.* **125**, 2019–030923 (2020).

25. Baldwin, M. P. et al. Sudden stratospheric warmings. *Rev. Geophys.* **59**, 2020–000708 (2021).

26. Domeisen, D. I. V. et al. The role of the stratosphere in subseasonal to seasonal prediction: 1. predictability of the stratosphere. *J. Geophys. Res. Atmos.* **125**, 2019–030920 (2020).

27. Li, Y., Kirchengast, G., Schwärz, M., Ladstädter, F. & Yuan, Y. Monitoring sudden stratospheric warmings using radio occultation: a new approach demonstrated based on the 2009 event. *Atmos. Measure. Tech.* **14**, 2327–2343 (2021).

28. Baldwin, M. P. & Dunkerton, T. J. Stratospheric harbingers of anomalous weather regimes. *Science* **294**, 581–584 (2001).

29. Scaife, A. A. et al. Long-range prediction and the stratosphere. *Atmos. Chem. Phys.* **22**, 2601–2623 (2022).

30. Charlton, A. & Polvani, L. A new look at stratospheric sudden warmings. Part I: Climatology and modeling benchmarks. *J. Clim.* **20**, 449–469 (2007).

31. Hitchcock, P. & Simpson, I. R. The downward influence of stratospheric sudden warmings. *J. Atmos. Sci.* **71**, 3856–3876 (2014).

32. Maycock, A. C., Masukwedza, G. I. T., Hitchcock, P. & Simpson, I. R. A regime perspective on the North Atlantic Eddy-driven jet response to sudden stratospheric warmings. *J. Clim.* **33**, 3901–3917 (2020).

33. Dai, Y. & Hitchcock, P. Understanding the basin asymmetry in surface response to sudden stratospheric warmings from an ocean-atmosphere coupled perspective. *J. Clim.* **34**, 8683–8698 (2021).

34. Ayarzagüena, B. et al. Stratospheric connection to the abrupt end of the 2016/2017 Iberian Drought. *Geophys. Res. Lett.* **45**, 12639–12646 (2018).

35. Kolstad, E. W., Breiteig, T. & Scaife, A. A. The association between stratospheric weak polar vortex events and cold air outbreaks in the Northern Hemisphere. *Quart. J. R. Meteorol. Soc.* **136**, 886–893 (2010).

36. King, A. D., Butler, A. H., Jucker, M., Earl, N. O. & Rudeva, I. Observed relationships between sudden stratospheric warmings and European climate extremes. *J. Geophys. Res. Atmos.* **124**, 13943–13961 (2019).

37. Lehtonen, I. & Karpechko, A. Y. Observed and modeled tropospheric cold anomalies associated with sudden stratospheric warmings. *J. Geophys. Res. Atmos.* **121**, 1591–1610 (2016).

38. Kok, U. F. et al. Contribution of the world's main dust source regions to the global cycle of desert dust. *Atmos. Chem. Phys.* **21**, 8169–8193 (2021).

39. Hurrell, J. W. Decadal trends in the north atlantic oscillation: regional temperatures and precipitation. *Science* **269**, 676–679 (1995).

40. Domeisen, D. I. V. Estimating the frequency of sudden stratospheric warming events from surface observations of the North Atlantic Oscillation. *J. Geophys. Res. Atmos.* **124**, 3180–3194 (2019).

41. Zhao, A., Ryder, C. L. & Wilcox, L. J. How well do the CMIP6 models simulate dust aerosols? *Atmos. Chem. Phys.* **22**, 2095–2119 (2022).

42. Gkikas, A., Giannaros, T., Kotroni, V. & Lagouvardos, K. Assessing the radiative impacts of an extreme desert dust outbreak and the potential improvements on short-term weather forecasts: The case of February 2015. *Atmos. Res.* **226**, 152–170 (2019).

43. Dayan, U., Heffter, J., Miller, J. & Gutman, G. Dust intrusion events into the Mediterranean Basin. *J. Appl. Meteorol. Climatol.* **30**, 1185–1199 (1991).

44. Dayan, U. & Levy, I. The influence of meteorological conditions and atmospheric circulation types on PM₁₀ and visibility in Tel Aviv. *J. Appl. Meteorol.* **44**, 606–619 (2005).

45. Bodenheimer, S., Lensky, I. M. & Dayan, U. Characterization of Eastern Mediterranean dust storms by area of origin; North Africa vs. Arabian Peninsula. *Atmos. Environ.* **198**, 158–165 (2019).

46. Helgren, D. M. & Prospero, J. M. Wind velocities associated with dust deflation events in the Western Sahara. *J. Clim. Appl. Meteorol.* **26**, 1147–1151 (1987).

47. Marticorena, B. et al. Temporal variability of mineral dust concentrations over West Africa: analyses of a pluriannual monitoring from the AMMA Sahelian Dust Transect. *Atmos. Chem. Phys.* **10**, 8899–8915 (2010).

48. Papadogeorgou, G., Kioumourtzoglou, M.-A., Braun, D. & Zanobetti, A. Low levels of air pollution and health: effect estimates, methodological challenges, and future directions. *Curr. Environ. Health Rep.* **6**, 105–115 (2019).

49. Mendoza, D. L. et al. Impact of low-level fine particulate matter and ozone exposure on absences in K-12 students and economic consequences. *Environ. Res. Lett.* **15**, 114052 (2020).

50. Yazdi, M. D. et al. Long-term effect of exposure to lower concentrations of air pollution on mortality among US Medicare participants and vulnerable subgroups: a doubly-robust approach. *Lancet Planet. Health* **5**, 689–697 (2021).

51. Menut, L., Chiapello, I. & Moulin, C. Predictability of mineral dust concentrations: The African Monsoon Multidisciplinary Analysis
first short observation period forecasted with CHIMERE-DUST. *J. Geophys. Res. Atmos.* **114**, 7 (2009).

52. Charlton-Perez, A. J., Huang, W. T. K. & Lee, S. H. Impact of sudden stratospheric warmings on United Kingdom mortality. *Atmos. Sci. Lett.* **22**, 1013 (2021).

53. Black, R. & McDaniel, B. A. The dynamics of Northern Hemisphere stratospheric final warming events. *J. Atmos. Sci.* **64**, 2932–2946 (2007).

54. Butler, A. H., Perez, A. C., Domeisen, D. I. V., Simpson, I. R. & Sjoberg, J. Predictability of northern hemisphere final stratospheric warmings and their surface impacts. *Geophys. Res. Lett.* **43**, 23 (2019).

55. Gelaro, R. et al. The modern-era retrospective analysis for research and applications, version 2 (MERRA-2). *J. Clim.* **30**, 5419–5454 (2017).

56. Hamilton, D. S. et al. Recent (1980 to 2015) trends and variability in daily-to-interannual soluble iron deposition from dust, fire, and anthropogenic sources. *Geophys. Res. Lett.* **47**, 2020–089688 (2020).

57. Randles, C. A. et al. The MERRA-2 Aerosol Reanalysis, 1980 Onward. Part I: System Description and Data Assimilation Evaluation. *J. Clim.* **30**, 6823–6850 (2017).

58. Zender, C. S., Bian, H. & Newman, D. Mineral Dust Entrainment and Deposition (DEAD) model: Description and 1990s dust climatology. *J. Geophys. Res. Atmos.* **108**, 14 (2003).

59. Chin, M. et al. Tropospheric aerosol optical thickness from the GOCART model and comparisons with satellite and sun photometer measurements. *J. Atmos. Sci.* **59**, 461–483 (2002).

60. Holben, B. N. et al. AERONET-A federated instrument network and data archive for aerosol characterization. *Remote Sens. Environ.* **66**, 1–16 (1998).

61. Compo, G. P. & Sardeshmukh, P. D. Storm track predictability on seasonal and decadal scales. *J. Clim.* **17**, 3701–3720 (2004).

62. Simpson, I. R., Hitchcock, P., Shepherd, T. G. & Scinocca, J. F. Southern Annular Mode Dynamics in Observations and Models. Part I: The Influence of Climatological Zonal Wind Biases in a Comprehensive GCM. *J. Clim.* **26**, 3953–3967 (2013).

63. Ogrosky, H. R. & Stechmann, S. N. Assessing the equatorial long-wave approximation: asymptotics and observational data analysis. *J. Atmos. Sci.* **72**, 4821–4843 (2015).

64. Butler, A. H., Sjoberg, J. P., Seidel, D. J. & Rosenlof, K. H. A sudden stratospheric warming compendium. *Earth Syst. Sci. Data* **9**, 63–76 (2017).

65. Borchers Arriagada, N. et al. Unprecedented smoke-related health burden associated with the 2019-20 bushfires in eastern Australia. *Med. J. Aust.* **213**, 282–283 (2020).

66. Zanobetti, A. & Schwartz, J. The effect of fine and coarse particulate air pollution on mortality: a national analysis. *Environ. Health Perspect.* **117**, 898–903 (2009).

67. Joss, M. K., Eeftens, M., Gintowt, E., Kappeler, R. & Künzli, N. Time to harmonize national ambient air quality standards. *Int. J. Public Health* **62**, 453–462 (2017).

68. Bouet, C. et al. Impact of desert dust on air quality: what is the meaningfulness of daily PM standards in regions close to the sources? The Example of Southern Tunisia. *Atmosphere* **10**, 452 (2019).

69. Kosaka, Y., Chowdary, J. S., Xie, S.-P., Min, Y.-M. & Lee, J.-Y. Limitations of seasonal predictability for summer climate over east Asia and the northwestern Pacific. *J. Clim.* **25**, 7574–7589 (2012).

## Acknowledgements

Y.D. acknowledges support from Cornell University. Support from the Swiss National Science Foundation through projects PPOOP2_170523 and PPOOP2_198896 to D.I.V.D. is gratefully acknowledged. M.K. and N.M. acknowledge support through the project PANACEA (MIS 5021516) funded by the Operational Programme "Competitiveness, Entrepreneurship and Innovation" (NSRF 2014-2020) and co-financed by Greece and the European Union Regional Development Fund. The French National Observatory Service INDAAF is supported by the INSU/CNRS, the IRD (Institut de Recherche pour le Développement), and the Observatoires des Sciences de l'Univers EFLUVE and Observatoire Midi-Pyrénées. We would like to thank the PIs and operators from France, Niger, Mali, and Senegal for maintaining the INDAAF stations (https://indaafobs-mip.fr).

## Author contributions

P.H. and D.I.V.D. initially conceived of the study. Y.D., P.H., N.M.M., and D.S.H. jointly designed the study. Y.D. performed the data analysis, created the figures, and wrote the manuscript. D.S.H. provided the CESM2 simulation and suggested the mortality calculation. L.L. suggested the examination of MERRA2 aerosol reanalysis. B.M., M.K., and N.M. provided the observational PM₁₀ concentrations. A.A.-O. calculated the mortality changes. All authors discussed the results and edited the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-35403-1.

**Correspondence** and requests for materials should be addressed to Peter Hitchcock.

**Peer review information** *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022