Saharan dust and association between particulate matter and case-specific mortality: a case-crossover analysis in Madrid (Spain)
© Díaz et al; licensee BioMed Central Ltd. 2012
Received: 11 October 2011
Accepted: 8 March 2012
Published: 8 March 2012
Saharan dust intrusions are a common phenomenon in the Madrid atmosphere, leading induce exceedances of the 50 μg/m3- EU 24 h standard for PM10.
We investigated the effects of exposure to PM10 between January 2003 and December 2005 in Madrid (Spain) on daily case-specific mortality; changes of effects between Saharan and non-Saharan dust days were assessed using a time-stratified case-crossover design.
Saharan dust affected 20% of days in the city of Madrid. Mean concentration of PM10 was higher during dust days (47.7 μg/m3) than non-dust days (31.4 μg/m3). The rise of mortality per 10 μg/m3 PM10 concentration were always largely for Saharan dust-days. When stratifying by season risks of PM10, at lag 1, during Saharan dust days were stronger for respiratory causes during cold season (IR% = 3.34% (95% CI: 0.36, 6.41) versus 2.87% (95% CI: 1.30, 4.47)) while for circulatory causes effects were stronger during warm season (IR% = 4.19% (95% CI: 1.34, 7.13) versus 2.65% (95% CI: 0.12, 5.23)). No effects were found for cerebrovascular causes.
We found evidence of strongest effects of particulate matter during Saharan dust days, providing a suggestion of effect modification, even though interaction terms were not statistically significant. Further investigation is needed to understand the mechanism by which Saharan dust increases mortality.
The increased levels of particulate matter in a given geographical area are mainly influenced by the intrusion of natural origin such as those related to the advection of dust from the desert . Obviously, the areas which are closer to the great deserts are the most affected by this type of events. In the case of the Saharan desert, Southern European countries are largely influenced by dust intrusions . Dust particles can also be transported over long distances by atmospheric circulation, reaching the Near East end the Americas .
Saharan dust events can contribute to exceedances of PM10 daily European Union Limit of 50 μg/m3 [2, 4]. Moreover, these particles can carry biological material that makes them potentially harmful to health . Among the diseases that would be affected by dust intrusions are those related to exacerbations of respiratory diseases, especially paediatric asthma [6–8].
Recently, this subject is rising from studies conducted in Southern Europe, although they have not given consistent results. In Spain, a study conducted in Barcelona  reported strongest effects of coarse particles (PM10-2.5) on total daily mortality during Saharan dust intrusions. In Italy, a study conducted in Rome  also found strongest effects of PM10 and coarse particles on daily mortality for cardiovascular causes, although others in Emilia-Romagna region  did not found an increase of risk of total and case-specific mortality due to PM10 during Saharan dust days. Lastly, studies in Athens did not found effect modification of PM10 on total and case-specific mortality , although they did found an increase of risk on paediatric asthma exacerbation .
The city of Madrid has high levels of pollution of anthropogenic origin, which are strongly increased those days with Saharan dust intrusions, and its geographical and meteorological conditions that did not facilitate the dispersion of pollutants. The aim of this study is to evaluate the effect of Saharan dust on the association between particulate matter and daily case-specific mortality in Madrid.
The city of Madrid constitutes a dense metropolitan area located on the central region of Spain. The main emission source of atmospheric pollutants is road traffic. Its geographical setting has climatological conditions directly linked to emission levels, with frequent anticyclonic situations in summer and winter, which impede the dispersion of pollutants. Saharan dust events on the central region of Spain generally occur 50 to 100 days per year, mainly between spring and autumn, mostly caused by an atmospheric depression West or Southwest of Portugal or by anticyclonic conditions over Algeria inducing Southern winds over the Iberian Peninsula .
Daily mortality was obtained from the Madrid Regional Inland Revenue Department (Consejería de Hacienda), which is the department responsible for mortality registry, from 1st January 2003 to 31st December 2005. The outcomes of interest were organic causes except accidents (ICD-10: A00-R99), respiratory causes (ICD-10: J00-J99), circulatory causes (ICD-10: I00-I99) and cerebrovascular mortality (ICD-10: I60-I69).
Daily mean concentrations of PM10 were obtained from the automated network of the Madrid's City Comprehensive Air-Pollution Monitoring, Forecasting and Information System (Red de Control de Contaminación Atmosférica del Ayuntamiento de Madrid). Daily mean levels of gaseous pollutants (NO2, SO2) and O3 were also collected.
Saharan dust intrusions, a dichotomous indicator variable used to show whether or not there had been a Saharan dust intrusion in the Madrid atmosphere on the date in question, with the pertinent data being obtained from the Directorate-General for Environmental Quality & Assessment (Dirección General de Calidad and Evaluación Ambiental) at the Ministry for the Environment and Rural & Marine Habitats (Ministerio de Medio Ambiente y Medio Rural y Marino) (available at: http://www.calima.ws/, last accessed February 17th, 2012). Briefly, an integrated methodological approach was used to identify days on which air masses from the Sahara-Sahel region were transported to the central region of Spain through back-trajectory analysis (Hysplit model), information from NRL SKIRON and BSCDREAM dust maps, and satellite images provided by the NASA SeaWiFS .
Daily mean temperature and relative humidity registered at the Madrid-Retiro Observatory, situated at city centre, were recorded by the National Meteorology Agency (Agencia Nacional de Meteorología).
Design and statistical analysis
The association of daily concentrations of PM10 with daily mortality was investigated using a case-crossover design . This uses the day on which the outcome of interest (mortality) occurs as a case day. Exposure on case days is compared with exposure on days on which the outcome of interest does not occur (control days) . A time-stratified approach used to represent exposure on control days from the same day of the week, month, and year as case days, minimising bias from time trends in the exposure series and from other short-term time-varying confounders .
We fitted a basic Poisson regression model to the daily mortality data that included potential confounders: weather, influenza, epidemics and time trends and seasonality. A natural cubic spline with 3 degrees of freedom was built to adjust for the potential confounding of mean temperature (4-day average) and a linear term to adjust for relative humidity (4-day average). We also controlled for gaseous pollutants and O3, by means of a 4-day average linear term, and for public holidays and influenza epidemics by means of categorical variables. Finally, we included a three-way interaction term between day of the week, month, and year to control for both seasonality and time trends. This choice was motivated by the need to replicate the adjustment made by the case-crossover design with the time-stratified approach for the selection of control days . To take into account possible overdispersion of daily counts of deaths, we used quasi-likelihood estimation.
The effects of exposure to PM10 were examined for the same day (lag 0) to 4 days after exposure (lag 4). Estimated effects are reported as IR% and 95% CI associated with a 10 μg/m3 increase of PM10. A dummy variable was created for the presence or absence of Saharan dust outbreaks at exposure days and its interaction with PM10 was fitted to test for effect modification by Saharan dust outbreaks. Because of most of Saharan dust intrusions occurs in Spain during the hot season, as previously reported , separate estimates were undertaken for the period May-September (warm season) and October-April (cold season). All analyses were done using Stata, release 11, statistical software (StataCorp, College Station, TX, 2010).
Summary statistics for daily case-specific mortality, PM10 and weather variables during Saharan dust and non-dust days in Madrid, 2003-2005.
Non-Saharan dust days (n = 878)
Saharan dust days (n = 218)
Association between particulate matter and daily case-specific mortality during Saharan dust and non-dust days, by season*
Non-Saharan dust days
Saharan dust days
p-value for interaction
Table 2 shows results also stratified by season at strongest lag previously found (lag 1). Risk of PM10 for organic causes was larger, and statistically significant, during cold season for non-dust days (1.32% [95% CI: 0.67, 1.97]) whilst during warm season risk was larger for Saharan dust days (2.53% [95% CI: 1.01, 4.08]) also providing evidence of effect modification with interaction terms marginally significant (p = 0.150 and p = 0.066 for cold and warm season, respectively). Moreover, during cold season effects of PM10 were stronger, and statistically significant, for respiratory causes being the larger risk for Saharan dust days (3.34% [95% CI: 0.36, 6.41] versus 2.87% [95% CI: 1.30, 4.47]) while for warm season effects were stronger for circulatory causes and risk was larger for Saharan dust days (4.19% [95% CI: 1.34, 7.13] versus 2.65% [95% CI: 0.12, 5.23]). However, none of the interaction terms were statistically significant. No effects were found for cerebrovascular causes for both seasons.
Saharan dust transport events reaching the city of Madrid showed strongest effects of daily concentrations of PM10 on case-specific mortality due to organic, respiratory and circulatory causes. Therefore, during Saharan dust-days these effects were larger in the cold season for respiratory causes and in the warm season for circulatory ones.
The distribution of daily concentrations of PM10 reported in our city, with higher levels on Saharan dust intrusion days, are consistent with the increased contribution of particulate matter from the Sahara because of the natural advections of this. Saharan dust intrusions tend to occur more frequently in the warmer months, which agrees with the frequency of synoptic-scale weather situations , mainly characterized by higher temperatures and low concentrations in relative humidity. Similar synoptic behaviour and frequency of Saharan dust intrusions have been described in the study conducted in the Emilia-Romagna region , while in Athens are less frequent and have no clear seasonal pattern .
The effects of PM10 on case specific mortality found in our study agreed with those previously reported in other European cities. The APHEA-2 study reported that an increase in PM10 by 10 μg/m3 was associated with increases of 0.76% in cardiovascular deaths and 0.58% in respiratory deaths . However, only few studies have assessed the effect of Saharan dust events on the association between PM10 and daily mortality, reporting inconsistent findings. Our results, showing an increase in mortality due to organic, respiratory and circulatory causes during Saharan dust days, are similar to those recently found in Rome . However, another study conducted in the Emilia-Romagna region  concluded that Saharan dust outbreaks are an independent risk factor that increases respiratory mortality. On the contrary, a latest study in Athens  did not found effect modification of PM10 on total and case-specific mortality due to dust events. In addition, studies focused on daily mortality related to dust of non-Saharan dust origin, like those in Washington  and Seoul , also fail to find any effect on Saharan dust events.
Currently, some studies have been focus on the effect of coarse-fraction (PM10-2.5). A study conducted in Barcelona  found a stronger effect of PM10-2.5 on daily total mortality during Saharan dust days, as well as the previously cited in Rome  that also found evidence of stronger effects on cardiovascular mortality. However, the health effects of coarse particles are still controversial, since most of the available studies were not able to isolate an independent effect of coarse particles from fine particles . Furthermore, previous studies of African dust size distributions showed that long distant transport reduces the amount of heavier and larger particles in the Saharan air masses, thus increasing the relative contribution of smaller particles . These data suggested than the usual definition of coarse fraction (PM10-2.5) usually adopted in air pollution monitoring also in epidemiological studies should be use with caution as specific marker of Saharan dust transports, which are probably better overall described by PM10-1 . However, we did not consider other PM fractions since from a previous study we did not found differences between Saharan and non-Saharan dust days for PM2.5 , and no PM1 levels were monitored during the study period. Nevertheless, the issue of harmful effects of Saharan dust has recently been addressed by taking PM10 into account, as recently did in the Emilia-Romagna region  and in Athens . Thus, our study is entirely comparable with those in terms of outcome and exposure evaluated, increasing the ongoing evidence on the mortality induced by the influx of Saharan dust in southern Europe.
A major strength of our study is that seasonal dependency, only conducted previously in the study conducted in the Emilia-Romagna region , has also been addressed. We found strongest effects of PM10 in daily mortality by circulatory causes during the warm season, and on respiratory causes during cold season. This confirms the seasonal pattern of case-specific mortality in Madrid, described in earlier studies more than two decades ago . This seasonal dependency was not found in the Emilia-Romagna region . This could indicate that the activity of biological agents in Saharan air masses differs depending on the season and the environmental conditions experienced by air masses travelling from the Sahara to Southern Europe [4, 11]. An interaction with other environmental factors not included in the analyses could be another feasible explanation.
As limitation, we have not addressed for contextual socioeconomic status and other potential modifiers, such as age and sex. However, the study design that would have meant a spatial adjustment had been very different to that conducted in our study . From a pathophysiological point of view there it is well known that the elderly group is most vulnerable to the short-term effects of air pollution on health , we have not considered age groups due to it was not available. Furthermore, as well as in previous studies [10–12], we do not have available data on chemical and biological components or on the copresence of other pollutants in particulate matter during Saharan dust days, so cannot determine specific components that could account for stronger estimated effects on Saharan dust days, or confirm the natural or anthropogenic origin of particulate matter. Only the study in Barcelona  reported data on chemical composition of fine and coarse particles for Saharan dust and non-dust days, concluding that although coarse particles seem to be more hazardous during Saharan dust days, differences in chemical composition did not explain these observations. Inconsistencies in the existing evidence suggests that characteristics associated with the short-term effects of desert dust should be further studied, with special attention to their composition, biological properties, and anthropogenic components that may be associated with the dust particles .
However, the results found in our study jointly with the knowledge of a possible biological mechanism that indicates that particulate material from the desert carry biological material which makes them particularly harmful to health [5, 26, 27], as well as the findings from other Mediterranean cities, such as Barcelona  and Rome , with similar settings like Madrid, make obvious the need to conduct further studies with a similar methodology in Southern Europe.
Our study adds to existing evidence of short-term effects of PM10 on the risk of daily case-specific mortality, especially for organic, respiratory and circulatory causes. During Saharan dust intrusion we also observed effects of PM10 on mortality due to respiratory causes in the cold season and to circulatory causes in the warm one. These shows the need that air quality standards must consider the potential harmful effects of dust from natural sources. Furthermore, additional research is needed to understand the mechanism by which particulate matter from desert sources can increase mortality.
International Classification of Diseases 10th revision
Percentage increases in risk of death
The study was supported by the Health Research Fund, Subdirectorate-General for Assessment and Development of Research (Project PI080354) and by the Carlos III Institute of Health (Project SEPY1453/07).
- Rodríguez S, Querol X, Alastuey A, Viana MM, Mantilla E: Events affecting levels and seasonal evolution of airborne matter concentrations in the Western Mediterranean. Environ Sci Technol. 2003, 37: 216-222. 10.1021/es020106p.View Article
- Rodriguez S, Querol X, Alastuey A, Kallos G, Kakaliagou O: Saharan dust contributions to PM10 and TSP levels in Southern and Eastern Spain. Atmos Environ. 2001, 35: 2433-2447. 10.1016/S1352-2310(00)00496-9.View Article
- Goudie AS, Middleton NJ: Saharan dust storms: nature and consequences. Earth Sci Rev. 2001, 56: 179-204. 10.1016/S0012-8252(01)00067-8.View Article
- Matassoni L, Pratesi G, Centioli D, Cadoni F, Malesani P, Caricchia AM, di Bucchianico AD: Saharan dust episodes in Italy: influence on PM10 daily limit value (DLV) exceedances and the related synoptic. J Environ Monit. 2009, 11: 1586-1594. 10.1039/b903822a.View Article
- Griffin DW: Atmospheric movement of microorganisms in clouds of desert dust and implications for human health. Clin Microbiol Rev. 2007, 20: 459-477. 10.1128/CMR.00039-06.View Article
- Gyan K, Henry W, Lacaille S, Laloo A, Lamsee-Ebanks C, McKay S, Antoine RM, Monteil MA: African dust clouds are associated with increased paediatric asthma accident and emergency admissions on the Caribbean island of Trinidad. Int J Biometeorol. 2005, 49: 371-376. 10.1007/s00484-005-0257-3.View Article
- Kanatini KT, Ito I, Al-Delaimy WK, Adachi Y, Mathews WC, Ramsdell JW, Toyama Asian Desert Dust and Asthma Study Team: Desert dust exposure is associated with increased risk of asthma hospitalization in children. Am J Resp Crit Care Med. 2010, 182: 1475-1481. 10.1164/rccm.201002-0296OC.View Article
- Samoli E, Nastos PT, Paliatsos AG, Katsouyanni K, Priftis KN: Acute effects of air pollution on pediatric asthma exacerbation: Evidence of association and effect modification. Environ Res. 2011, 111: 418-424. 10.1016/j.envres.2011.01.014.View Article
- Pérez L, Tobías A, Querol X, Künzly N, Pey J, Alastuey A, Viana M, Valero N, González-Cabré M, Sunyer J: Coarse particles from Saharan dust and daily mortality. Epidemiol. 2008, 19: 800-807. 10.1097/EDE.0b013e31818131cf.View Article
- Mallone S, Stafoggia M, Faustini A, Gobbi GP, Marconi A, Forastiere F: Saharan Dust and Associations between Particulate Matter and Daily Mortality in Rome, Italy. Environ Health Perspect. 2011, 119: 1409-1414.View Article
- Sajani SZ, Miglio R, Bonasoni P, Cristofanelli P, Marinoni A, Sartini C, Goldoni CA, De Girolamo G, Lauriola P: Saharan dust and daily mortality in Emilia-Romagna (Italy). Occup Environ Med. 2011, 68: 446-451. 10.1136/oem.2010.058156.View Article
- Samoli E, Kougea E, Kassomenos P, Analitis A, Katsouyanni K: Does the presence of desert dust modify the effect of PM10 on mortality in Athens, Greece?. Sci Total Environ. 2011, 409: 2049-2054. 10.1016/j.scitotenv.2011.02.031.View Article
- Maclure M: The case-crossover design: a method for studying transient effects on the risk of acute events. Am J Epidemiol. 1991, 133: 144-153.
- Janes H, Sheppard L, Lumley T: Case-crossover analyses of air pollution exposure data: referent selection strategies and their implications for bias. Epidemiol. 2005, 16: 717-726. 10.1097/01.ede.0000181315.18836.9d.View Article
- Lu Y, Symons JM, Geyh AS, Zeger SL: An approach to checking case-crossover analyses based on equivalence with time-series methods. Epidemiol. 2008, 19: 169-175. 10.1097/EDE.0b013e3181632c24.View Article
- Trigo RM, García-Herrera R, Díaz J, Trigo IF, Valente MA: How exceptional was the early August 2003 heat wave in France?. Geophys Res Lett. 2005, 32: L10701-View Article
- Analitis A, Katsouyanni K, Dimakopoulou K, Samoli E, Nikoloulopoulos AK, Petasakis Y, Touloumi G, Schwartz J, Anderson HR, Cambra K, Forastiere F, Zmirou D, Vonk JM, Clancy L, Kriz B, Bobvos J, Pekkanen J: Short-term effects of ambient particles on cardiovascular and respiratory mortality. Epidemiology. 2006, 17: 230-233. 10.1097/01.ede.0000199439.57655.6b.View Article
- Schwartz J, Norris G, Larson T, Sheppard L, Claiborne C, Koenig J: Episodes of high coarse particles concentrations are not associated with increased mortality. Environ Health Perspect. 1999, 107: 339-342. 10.1289/ehp.99107339.View Article
- Lee JT, Son JY, Cho YS: A Comparison of mortality related to urban air particles between periods with Asian dust days and without Asian dust days in Seoul Korea, 2000-2004. Environ Res. 2007, 105: 409-413. 10.1016/j.envres.2007.06.004.View Article
- Brunekreef B, Forsberg B: Epidemiological evidence of effects of coarse airborne particles on health. Eur Respir J. 2005, 26: 309-318. 10.1183/09031936.05.00001805.View Article
- Prospero JM, Olmez I, Ames M: Al and Fe in PM2.5 and PM10 in suspended particles in south-central Florida: the impact of the long range transport of African mineral dust. Water Air Soil Pollut. 2001, 125: 291-317. 10.1023/A:1005277214288.View Article
- Tobías A, Pérez L, Díaz J, Linares C, Pey J, Alastruey A, Querol X: Short-term effects of particulate matter on total mortality during Saharan dust outbreaks: A case-crossover analysis in Madrid (Spain). Sci Total Environ. 2011, 386-389. 412-413
- Alberdi JC, Díaz J: Modelización de la mortalidad diaria en la Comunidad de Madrid de 1986 a 1991. Gac Sanit. 1997, 11: 9-15.View Article
- Barceló MA, Saez M, Saurina C: Spatial variability in mortality inequalities, socioeconomic deprivation, and air pollution in small areas of the Barcelona Metropolitan Region, Spain. Sci Total Environ. 2009, 407: 5501-5523. 10.1016/j.scitotenv.2009.07.028.View Article
- Aga E, Samoli E, Touloumi G, Anderson HR, Cadum E, Forsberg B, Goodman P, Goren A, Kotesovec F, Kriz B, Macarol-Hiti M, Medina S, Paldy A, Schindler C, Sunyer J, Tittanen P, Wojtyniak B, Zmirou D, Schwartz J, Katsouyanni K: Short-term effects of ambient particles on mortality in the elderly: results from 28 cities in the APHEA2 project. Eur Respir J Suppl. 2003, 40: 28s-33s.View Article
- McCarthy M: Dust clouds implicated in spread of infection. Lancet. 2001, 358: 478-10.1016/S0140-6736(01)05677-X.View Article
- Kellogg C, Griffin D: Aerobiology and the global transport of desert dust. Trends Ecol Evol. 2006, 21: 638e44-View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.