Influence of Saharan Dust Intrusions on Respiratory Medication Dispensing

1. Introduction

The exposome, first introduced by Wild [1], refers to the totality of environmental exposures encountered across the lifespan and their biological interactions. This concept complements genomic research by offering a broader framework to understand complex, multifactorial diseases [2,3,4]. Among the diverse external exposures, those related to air quality are particularly critical. In this context, desert dust intrusions—such as those originating from the Sahara—represent episodic yet substantial components of the atmospheric exposome. These events contribute to sharp increases in particulate matter and transport not only mineral particles but also microorganisms, allergens, and toxic compounds, with significant implications for human health [5,6].

Desert dust intrusions are large-scale atmospheric phenomena driven by strong winds and arid conditions, capable of lofting fine and coarse particles into the atmosphere and carrying them thousands of kilometers. As a result, air quality can deteriorate far from the emission source, affecting populations across entire regions [7,8]. While dust has traditionally been considered a natural element of Earth’s geochemical cycles, the rising frequency and intensity of dust storms—exacerbated by land degradation, desertification, and climate change—have transformed them into an emerging public health challenge [9].

One of the most relevant consequences of dust exposure is its impact on respiratory health. Fine particulate matter (PM₁₀ and PM₂.₅), abundant during dust events, can penetrate the lower respiratory tract and trigger inflammation, oxidative stress, and bronchial hyperresponsiveness [10,11]. Epidemiological studies have consistently linked dust episodes with higher rates of asthma and chronic obstructive pulmonary disease (COPD) exacerbations, emergency department visits, hospital admissions, and mortality, particularly among vulnerable populations [12,13,14]. Other natural hazards such as wildfires and volcanic eruptions also produce airborne pollutants—fine particulates, carbon monoxide, irritant gases—that contribute to respiratory morbidity. In recent years, wildfires have become more frequent and intense, causing episodes of severe air pollution across wide geographic areas [15,16]. In addition, desert dust can carry biological material including allergens, bacteria, and fungi, which may further aggravate symptoms in patients with asthma or immunologic susceptibility [17,18].

Real-world data increasingly support the notion that these natural airborne pollution events are followed by short-term surges in inhaled medication use, particularly short-acting beta-agonists (SABAs), inhaled corticosteroids (ICS), and combination maintenance therapies such as ICS–LABA and long-acting antimuscarinics [19,20]. In regions affected by Saharan dust—such as Southern Europe and the Middle East—dispensing records have shown a rise in rescue and maintenance inhaler use within days of exposure peaks [21,22]. Similar trends have been observed during wildfire seasons in North America and Australia, where pharmacy sales of SABAs increase in parallel with acute respiratory distress [23,24]. These patterns highlight both the clinical burden of environmental episodes and the utility of inhaler dispensing as a proxy for population-level respiratory impact. Moreover, inhaler utilization trends reveal disparities in access and preparedness: during periods of elevated pollutant exposure, patients with asthma or COPD may rely more heavily on rescue inhalers, especially when access to routine care is disrupted. This raises important questions regarding the adequacy of preventive strategies, the timeliness of public health alerts, and the resilience of healthcare systems during environmental emergencies.

Previous research suggests that the respiratory impact of air pollution may be more accurately captured through patterns of medication use than by symptom reporting or exacerbation counts [25,26]. Using Spain’s universal prescription registry, we investigated Saharan dust intrusions (SDI) in the Canary Islands, a setting uniquely suited to this question given its frequent exposure to SDI [27,28] and elevated prevalence of asthma and COPD [29,30,31]. The objective was to determine whether SDI are associated with changes in dispensing of SABA reliever therapy and ICS–LABA maintenance therapy.

2. Materials and Methods

2.1. Definition of ICS-LABA and SABA Sales

Sales data (number of units dispensed) of Inhaled Corticosteroids plus Long-Acting Beta-Agonists (ICS-LABA) and Short-Acting Beta-Agonists (SABA) were retrieved from 29 and 31 pharmacies located in the two provinces of the Canary Islands (Spain), Santa Cruz de Tenerife (TF) and Las Palmas de Gran Canaria (GC), respectively. The data timeline was established from June 2017 to May 2022 (60 months of follow-up). The monthly sum of sales for each drug was calculated, resulting in 60 monthly observations for each province. Sales were assigned to the month following dispensing to better reflect the effective consumption period.

2.2. Definition of SDI

Days with SDI were identified considering the daily concentration of Particulate Matter with a diameter ≤ 10 µm (PM₁₀), retrieved from the Canary Islands Air Quality Control and Monitoring Network [32]. Daily mean PM₁₀ concentrations were calculated from 30 and 21 monitoring stations in TF and GC, respectively, and the timeline included the same period of 60 months defined above. The presence of SDI was declared when the mean PM₁₀ concentration exceeded 40 µg/m³, a threshold at which air quality is classified as “regular” according to the criteria defined by the Canary Islands Air Quality Control and Monitoring Network [32].

Since different SDI atmospheric scenarios take place in the Canary Islands depending on the season, the effect of seasonality was also considered. During the winter (October 1^st to March 31^st), Saharan dust is mainly carried by low-level continental African trade winds (Harmattan winds), reaching altitudes of less than 700 meters above sea level and mostly affecting urban populated areas. However, in the summer (April 1^st to September 30^th), dust is mainly transported by the northern branch of the high-altitude Saharan Air Layer, reaching higher altitudes where urban-populated regions are almost absent [33,34,35].

2.3. Statistical Analysis

Data analyses were performed using R software (version 4.3.2; R Core Team, 2023) within RStudio (version 2023.09.1+494; Posit Software, PBC). Outliers were identified using the interquartile range (IQR) criteria (datapoints 1.5 × IQR above the third quartile or below the first quartile were excluded). Outliers were identified independently for each analysis group and removed prior to statistical testing. Data normality was assessed using the Kolmogorov-Smirnov test (n > 50). Linear regression models were applied to evaluate the association between ICS-LABA or SABA sales with the presence/absence of SDI, the total number of days with SDI per month, and the monthly average PM₁₀ concentration. Additionally, linear regression models were adjusted for seasonality (winter or summer), defined as described above. Statistical significance was declared for those comparisons with p-value < 0.05, and correlations were evaluated with the Pearson’s coefficient (r).

2.4. Ethical Approval and Consent to Participate

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Ethics Committee of CEIm Hospital Universitario de Canarias, Tenerife, Spain with the reference number P.I.-2017/72 on October 30, 2017. The study was conducted using aggregated pharmaceutical sales records and publicly available air quality data therefore, no identifiable personal information was collected, and all data were anonymized at the pharmacy level. The requirement for informed consent was waived, as the study involved only secondary, non-identifiable data and posed no risk to individuals.

3. Results

3.1. Effect of SDI over SABA and ICS-LABA Sales in the Canary Islands

The monthly sales of SABA and ICS-LABA during the 60-months follow-up was studied in the two provinces of the Canary Islands (Table 1).

Results showed that SABA sales were 14.8% lower in TF, compared with GC (Figure 1A), while ICS-LABA sales were 10.9% higher (Figure 1B).

These results suggest a greater adherence on rescue medication in GC, while TF showed slightly higher use of maintenance therapy, and possibly better adherence to preventive treatment, pointing to differences in disease management practices, or patient adherence. Therefore, the effect of SDI over SABA and ICS-LABA sales was independently studied in both provinces, comparing the sales that took place in months where PM₁₀ reached 40 µg/m³ in at least one day, with the sales in months with absence of SDI. SABA sales were similar between the two groups, both in TF and GC, meaning that SDI did not affect the dispensation of rescue medication in the two provinces studied (Figure 2A). Interestingly, months with presence of SDI showed significant higher ICS-LABA sales in the two provinces (Figure 2B). In the case of GC, the increment reaches 10.2%, while in the province of TF, ICS-LABA sales were increased by 5.7%. These results indicate that SDI exposure is not translated in a higher adherence of rescue medication but is associated with and slightly higher use of maintenance therapy in both provinces.

3.2. Effect of SDI Frequency and Intensity over SABA and ICS-LABA Sales

First, the effect of SDI frequency was evaluated by testing the associations between the number of days with SDI per month and the drugs sales, using linear regression models. No associations were found in the case of SABA (not shown), while ICS-LABA sales showed a significant association with the number of days with SDI per month, both in GC and TF (Figure 3A). To investigate the effect of SDI intensity, associations between the monthly average of PM₁₀ concentration and the drug sales were also evaluated. Again, SABA sales showed no associations with PM₁₀ average concentration (not shown), while significant associations were confirmed for ICS-LABA in both provinces (Figure 3B). Nevertheless, although associations were significantly supported, linear correlations were certainly low for both provinces, when evaluated by the Pearson’s r coefficient (Figure 3).

Overall, these results support a significant but weak correlation between the frequency and intensity of SDI and ICS-LABA sales, but not for SABA. This tendency was observed in both provinces independently, suggesting that the more the population is exposed to SDI, the greater the demand for maintenance therapy, but not of rescue medication.

3.3. Effect of SDI Seasonality over SABA and ICS-LABA Sales

Seasonality strongly modulates SDI in the Canary Islands. In winter, SDI mostly affects populated areas, while in summer it reaches more elevated areas, where urban centres are absent. Therefore, drug sales were compared only for the months with SDI presence, between the winter and summer seasons (Figure 4). Results showed that ICS-LABA sales were significantly higher during winter months with SDI presence than in summer months with SDI (Figure 4A). This effect was observed in both provinces independently, reaching 9.6% and 14.3% in GC and TF, respectively (Figure 4A). Moreover, when the comparison between seasons were repeated, but including those months without SDI only, no differences were observed (p_GC = 0.078; p_TF = 0.257). These results indicate that the increase of ICS-LABA sales detected in the Canary Islands as result of SDI is mostly associated to the events that take place during the winter.

When restricted to months with SDI, SABA sales remained significantly higher in winter compared with summer in both provinces (18.8% in GC and 18.9% in TF; Figure 4B). Comparable seasonal differences were also observed during months without SDI (pGC = 0.035; pTF = 0.039). These findings indicate that the increase in SABA sales is attributable to seasonality and not to the presence of SDI.

Finally, based on these results, linear regression models showed in Figure 2 were adjusted to include seasonality as a covariate. Findings showed that the increase in ICS-LABA sales during SDI remained significant even after adjustment by season in both provinces (p_GC = 0.022, p_TF = 0.027). This analysis showed an increase of 1,014.0 ± 431.6 and 1,181.7 ± 520.1 units of ICS-LABA in those months with SDI presence for GC and TF, respectively.

4. Discussion

Although desert dust is an inherent component of the Earth’s climate system, its intensification in recent decades—fueled by desertification, land degradation, and climate change—has emerged as a major public health concern [36]. The Sahara, which alone accounts for more than half of global dust emissions, is the predominant source, with recurrent plumes that often reach Southern Europe and particularly affect the Canary Islands [37,38].

Desert dust intrusions expose the respiratory epithelium to a heterogeneous mixture of mineral particles, metals, organic compounds, and microbial agents that compromise barrier integrity and trigger inflammatory cascades. Experimental models show that dust particles induce oxidative stress and activate the NLRP3 inflammasome, promoting release of IL-1β and epithelial alarmins such as IL-33 and TSLP [39,40]. These responses amplify type-2 and neutrophilic inflammation, enhance bronchial hyperresponsiveness, and impair mucociliary clearance, thereby worsening airway obstruction in both asthma and COPD [41,42]. Repeated or sustained epithelial injury can lower the threshold for symptoms such as wheezing, cough, and dyspnea, prompting greater reliance on inhaler therapies. Epidemiological studies linking dust exposure to exacerbations and healthcare use are consistent with these mechanistic insights and support the interpretation of inhaler utilization as a clinical marker of dust-related exposure burden within the exposome framework [43,44].

In this study, SDI in the Canary Islands were associated with a measurable increase in ICS–LABA dispensations but not in SABA use. This association persisted after adjustment for seasonal variation, indicating that dust episodes exert an independent effect on maintenance therapy utilization. By contrast, the winter increase in SABA sales was attributable to seasonal trends rather than SDI exposure. These results align with previous evidence linking natural air pollution events to heightened medication demand yet the differentiation between reliever and maintenance therapy is notable [45,46,47], While SABA use primarily reflects acute symptom relief, the rise in ICS–LABA sales suggests escalation of preventive therapy, whether initiated by patients or reinforced by clinicians. The modest but significant correlation between ICS–LABA sales and both frequency and intensity of SDI further supports this interpretation.

The absence of a parallel rise in SABA use invites two complementary explanations. First, SDI events in the Canary Islands have consistently been associated with increased morbidity from asthma and COPD, including higher admissions and exacerbations, often with a temporal lag suggestive of subacute inflammation rather than immediate bronchospasm [48]. This delayed effect may encourage greater reliance on maintenance therapy to stabilize symptom burden. Second, the growing adoption of Maintenance and Reliever Therapy (MART/SMART), in which ICS–formoterol serves as both controller and reliever, could shift symptom-driven medication demand into ICS–LABA dispensations rather than SABA, in line with current guideline recommendations [49,50]. Together, these mechanisms suggest that SDI exposure increases both the need for intensified anti-inflammatory therapy and the demand for symptom relief covered under ICS–LABA prescriptions, while leaving SABA utilization largely unaffected.

From a clinical and public health perspective, these findings carry several implications. Inhaler dispensing data represent a valuable real-world indicator of the population impact of environmental exposures, complementing hospital admissions and patient-reported outcomes. The stronger effect observed in winter highlights the role of seasonal atmospheric dynamics in shaping exposure risk and enhances the value of integrating environmental alerts into preventive respiratory care strategies. Moreover, the reliance on maintenance rather than rescue therapy raises important questions regarding adherence, prescribing practices, and whether SDI events prompt sustained adjustments in asthma and COPD management.

Limitations of the present study must be acknowledged. Dispensing data indicate medication supply rather than actual use and may be affected by prescribing practices, stockpiling, or pharmacy availability [51,52,53]. The absence of individual-level clinical outcomes prevents assessing whether increased ICS–LABA use improves control or reduces exacerbations. Dust composition—including allergens, bacteria, and fungi—was also not characterized, despite its potential impact on respiratory outcomes in atopic populations [54,55]. Future research should apply advanced time-series models (e.g., distributed lag or autoregressive) to capture delayed and cumulative effects of SDI, and link dispensing patterns with clinical outcomes such as exacerbations, emergency visits, and lung function. Incorporating chemical and biological dust analyses may further clarify exposure–response mechanisms and support targeted preventive strategies in regions recurrently affected by desert dust.

5. Conclusions

This study provides novel evidence that Saharan dust intrusions significantly increase the use of ICS–LABA maintenance therapy in the Canary Islands, independent of seasonal variation. In contrast, SABA utilization patterns reflected only seasonal trends, not dust exposure. These findings highlight important differences in patient and prescriber responses to environmental events and underscore the utility of pharmaceutical dispensing records as a sensitive, real-world proxy for respiratory burden. Public health systems in regions affected by desert dust should anticipate temporary increases in maintenance therapy demand, particularly during winter episodes when populated areas are most exposed. Integrating dust forecasts into early-warning systems and patient education programs may improve disease control, reduce exacerbations, and strengthen healthcare preparedness for future environmental challenges.