Desert dust is assumed to have substantial adverse effects on human health. However, the epidemiologic evidence is still inconsistent, mainly because previous studies used different metrics for dust ...exposure and its corresponding epidemiologic analysis. We aim to provide a standardized approach to the methodology for evaluating the short-term health effects of desert dust.
We reviewed the methods commonly used for dust exposure assessment, from use of a binary metric for the occurrence of desert dust advections to a continuous one for quantifying particulate matter attributable to desert dust. We presented alternative time-series Poisson regression models to evaluate the dust exposure-mortality association, from the underlying epidemiological and policy-relevant questions. A set of practical examples, using a real dataset from Rome, Italy, illustrate the different modeling approaches.
We estimate substantial effects of desert dust episodes and particulate matter with diameter <10 μm (PM10) on daily mortality. The estimated effect of non-desert PM10 was 1.8% (95% confidence interval CI = 0.4, 3.2) for a 10 μg/m rise of PM10 at lag 0 for dust days, 0.4% (95% CI = -0.1, 0.8) for non-dust days, and 0.6% (95% CI = -0.5, 2.1) for desert PM10.
The standardized modeling approach we propose could be applicable elsewhere, in and near hot spots, which could lead to more consistent evidence on the health effects of desert dust from future studies.
Particulate matter (PM) air pollution is one of the major causes of death worldwide, with demonstrated adverse effects from both short-term and long-term exposure. Most of the epidemiological studies ...have been conducted in cities because of the lack of reliable spatiotemporal estimates of particles exposure in nonurban settings. The objective of this study is to estimate daily PM10 (PM < 10 μm), fine (PM < 2.5 μm, PM2.5) and coarse particles (PM between 2.5 and 10 μm, PM2.5–10) at 1-km2 grid for 2013–2015 using a machine learning approach, the Random Forest (RF). Separate RF models were defined to: predict PM2.5 and PM2.5–10 concentrations in monitors where only PM10 data were available (stage 1); impute missing satellite Aerosol Optical Depth (AOD) data using estimates from atmospheric ensemble models (stage 2); establish a relationship between measured PM and satellite, land use and meteorological parameters (stage 3); predict stage 3 model over each 1-km2 grid cell of Italy (stage 4); and improve stage 3 predictions by using small-scale predictors computed at the monitor locations or within a small buffer (stage 5). Our models were able to capture most of PM variability, with mean cross-validation (CV) R2 of 0.75 and 0.80 (stage 3) and 0.84 and 0.86 (stage 5) for PM10 and PM2.5, respectively. Model fitting was less optimal for PM2.5–10, in summer months and in southern Italy. Finally, predictions were equally good in capturing annual and daily PM variability, therefore they can be used as reliable exposure estimates for investigating long-term and short-term health effects.
•Estimates of fine and coarse particles at fine spatiotemporal scale are lacking in Italy•We applied a multistage random forest model combining PM data with satellite, land-use and meteorology•We imputed missing satellite AOD data using ensemble atmospheric models•We estimated daily PM10, PM2.5 and PM2.5-10 at a 1-km2 grid over Italy for the years 2013-2015•Our model displayed good CV fitting (R2=0.75 for PM10, R2=0.80 for PM2.5, R2=0.64 for PM2.5-10) and negligible bias
Few European studies have investigated the effects of long-term exposure to both fine particulate matter (≤ 2.5 µm; PM2.5) and nitrogen dioxide (NO2) on mortality.
We studied the association of ...exposure to NO2, PM2.5, and traffic indicators on cause-specific mortality to evaluate the form of the concentration-response relationship.
We analyzed a population-based cohort enrolled at the 2001 Italian census with 9 years of follow-up. We selected all 1,265,058 subjects ≥ 30 years of age who had been living in Rome for at least 5 years at baseline. Residential exposures included annual NO2 (from a land use regression model) and annual PM2.5 (from a Eulerian dispersion model), as well as distance to roads with > 10,000 vehicles/day and traffic intensity. We used Cox regression models to estimate associations with cause-specific mortality adjusted for individual (sex, age, place of birth, residential history, marital status, education, occupation) and area (socioeconomic status, clustering) characteristics.
Long-term exposures to both NO2 and PM2.5 were associated with an increase in nonaccidental mortality hazard ratio (HR) = 1.03 (95% CI: 1.02, 1.03) per 10-µg/m3 NO2; HR = 1.04 (95% CI: 1.03, 1.05) per 10-µg/m3 PM2.5. The strongest association was found for ischemic heart diseases (IHD) HR = 1.10 (95% CI: 1.06, 1.13) per 10-µg/m3 PM2.5, followed by cardiovascular diseases and lung cancer. The only association showing some deviation from linearity was that between NO2 and IHD. In a bi-pollutant model, the estimated effect of NO2 on mortality was independent of PM2.5.
This large study strongly supports an effect of long-term exposure to NO2 and PM2.5 on mortality, especially from cardiovascular causes. The results are relevant for the next European policy decisions regarding air quality.
Display omitted
•Long-term road traffic and railway noise are associated with most CVD causes of death.•Risk increases often start well below the WHO Environmental Noise guideline ...levels.•Associations are independent of air pollution.•Higher levels of noise intermittency are independently associated with each outcome.•Relative and absolute risk are higher in males compared to females.
Death from cardiovascular diseases (CVD) has been associated with transportation noise. This nationwide cohort, with state-of-the-art exposure assessment, evaluates these associations by noise source.
Road traffic, railway and aircraft noise for 2001 and 2011 were linked to 4.1 million adults in the Swiss National Cohort, accounting for address history. Mean noise exposure in 5-year periods was calculated. Time-varying Cox regression models, with age as timescale, were applied to all and cause-specific cardiovascular causes of death. Models included all three noise sources plus PM2.5, adjusted for individual and spatial covariates. Nighttime noise events for all sources combined (expressed as intermittency ratio or number of events) were considered in sensitivity analyses. Absolute excess risk was calculated by multiplying deaths/100,000 person-years by the excess risk (hazard ratio-1) within each age/sex group.
During a 15-year follow-up, there were 277,506 CVD and 34,200 myocardial infarction (MI) deaths. Associations (hazard ratio; 95%-CIs) for road traffic, railway and aircraft noise and CVD mortality were 1.029 (1.024–1.034), 1.013 (1.010–1.017), and 1.003 (0.996–1.010) per 10 dB Lden, respectively. Associations for MI mortality were a respective 1.043 (1.029–1.058), 1.020 (1.010–1.030) and 1.040 (1.020–1.060) per 10 dB Lden. Blood pressure-related, ischemic heart disease, and all stroke mortality were significantly associated with road traffic and railway noise, while ischemic stroke mortality was associated with aircraft noise. Associations were mostly linear, often starting below 40 dB Lden for road traffic and railway noise. Higher levels of noise intermittency were also independently associated with each outcome. While the absolute number of deaths attributed to noise increased with age, the hazard ratios declined with age. Relative and absolute risk was higher in males compared to females.
Independent of air pollution, transportation noise exposure is associated with all and cause-specific CVD mortality, with effects starting below current guideline limits.
In order to investigate associations between air pollution and adverse health effects consistent fine spatial air pollution surfaces are needed across large areas to provide cohorts with comparable ...exposures. The aim of this paper is to develop and evaluate fine spatial scale land use regression models for four major health relevant air pollutants (PM2.5, NO2, BC, O3) across Europe.
We developed West-European land use regression models (LUR) for 2010 estimating annual mean PM2.5, NO2, BC and O3 concentrations (including cold and warm season estimates for O3). The models were based on AirBase routine monitoring data (PM2.5, NO2 and O3) and ESCAPE monitoring data (BC), and incorporated satellite observations, dispersion model estimates, land use and traffic data. Kriging was performed on the residual spatial variation from the LUR models and added to the exposure estimates. One model was developed using all sites (100%). Robustness of the models was evaluated by performing a five-fold hold-out validation and for PM2.5 and NO2 additionally with independent comparison at ESCAPE measurements. To evaluate the stability of each model's spatial structure over time, separate models were developed for different years (NO2 and O3: 2000 and 2005; PM2.5: 2013).
The PM2.5, BC, NO2, O3 annual, O3 warm season and O3 cold season models explained respectively 72%, 54%, 59%, 65%, 69% and 83% of spatial variation in the measured concentrations. Kriging proved an efficient technique to explain a part of residual spatial variation for the pollutants with a strong regional component explaining respectively 10%, 24% and 16% of the R2 in the PM2.5, O3 warm and O3 cold models. Explained variance at fully independent sites vs the internal hold-out validation was slightly lower for PM2.5 (65% vs 66%) and lower for NO2 (49% vs 57%). Predictions from the 2010 model correlated highly with models developed in other years at the overall European scale.
We developed robust PM2.5, NO2, O3 and BC hybrid LUR models. At the West-European scale models were robust in time, becoming less robust at smaller spatial scales. Models were applied to 100 × 100 m surfaces across Western Europe to allow for exposure assignment for 35 million participants from 18 European cohorts participating in the ELAPSE study.
•Robust PM2.5, NO2, BC and O3 hybrid LUR models at a 100x100 m resolution for Western Europe were developed•Models included large scale SAT and CTM estimates and fine scale traffic and land use and were further improved with kriging•Models were robust in time at European scale, becoming less robust at smaller spatial scales.
Spatiotemporally resolved particulate matter (PM) estimates are essential for reconstructing long and short-term exposures in epidemiological research. Improved estimates of PM2.5 and PM10 ...concentrations were produced over Italy for 2013–2015 using satellite remote-sensing data and an ensemble modeling approach. The following modeling stages were used: (1) missing values of the satellite-based aerosol optical depth (AOD) product were imputed using a spatiotemporal land-use random-forest (RF) model incorporating AOD data from atmospheric ensemble models; (2) daily PM estimations were produced using four modeling approaches: linear mixed effects, RF, extreme gradient boosting, and a chemical transport model, the flexible air quality regional model. The filled-in MAIAC AOD together with additional spatial and temporal predictors were used as inputs in the three first models; (3) a geographically weighted generalized additive model (GAM) ensemble model was used to fuse the estimations from the four models by allowing the weights of each model to vary over space and time. The GAM ensemble model outperformed the four separate models, decreasing the cross-validated root mean squared error by 1–42%, depending on the model. The spatiotemporally resolved PM estimations produced by the suggested model can be applied in future epidemiological studies across Italy.
Empirical spatial air pollution models have been applied extensively to assess exposure in epidemiological studies with increasingly sophisticated and complex statistical algorithms beyond ordinary ...linear regression. However, different algorithms have rarely been compared in terms of their predictive ability.
This study compared 16 algorithms to predict annual average fine particle (PM2.5) and nitrogen dioxide (NO2) concentrations across Europe. The evaluated algorithms included linear stepwise regression, regularization techniques and machine learning methods. Air pollution models were developed based on the 2010 routine monitoring data from the AIRBASE dataset maintained by the European Environmental Agency (543 sites for PM2.5 and 2399 sites for NO2), using satellite observations, dispersion model estimates and land use variables as predictors. We compared the models by performing five-fold cross-validation (CV) and by external validation (EV) using annual average concentrations measured at 416 (PM2.5) and 1396 sites (NO2) from the ESCAPE study. We further assessed the correlations between predictions by each pair of algorithms at the ESCAPE sites.
For PM2.5, the models performed similarly across algorithms with a mean CV R2 of 0.59 and a mean EV R2 of 0.53. Generalized boosted machine, random forest and bagging performed best (CV R2~0.63; EV R2 0.58–0.61), while backward stepwise linear regression, support vector regression and artificial neural network performed less well (CV R2 0.48–0.57; EV R2 0.39–0.46). Most of the PM2.5 model predictions at ESCAPE sites were highly correlated (R2 > 0.85, with the exception of predictions from the artificial neural network). For NO2, the models performed even more similarly across different algorithms, with CV R2s ranging from 0.57 to 0.62, and EV R2s ranging from 0.49 to 0.51. The predicted concentrations from all algorithms at ESCAPE sites were highly correlated (R2 > 0.9). For both pollutants, biases were low for all models except the artificial neural network. Dispersion model estimates and satellite observations were two of the most important predictors for PM2.5 models whilst dispersion model estimates and traffic variables were most important for NO2 models in all algorithms that allow assessment of the importance of variables.
Different statistical algorithms performed similarly when modelling spatial variation in annual average air pollution concentrations using a large number of training sites.
•Multiple statistical algorithms with very different assumptions were compared.•Despite the difference in modeling frameworks, predictions among the models exhibit generally good agreement.•The use of an external evaluation dataset strengthens evaluation by cross-validation.
Standardized mortality surveillance data, capable of detecting variations in total mortality at population level and not only among the infected, provide an unbiased insight into the impact of ...epidemics, like COVID-19 (Coronavirus disease). We analysed the temporal trend in total excess mortality and deaths among positive cases of SARS-CoV-2 by geographical area (north and centre-south), age and sex, taking into account the deficit in mortality in previous months.
Data from the Italian rapid mortality surveillance system was used to quantify excess deaths during the epidemic, to estimate the mortality deficit during the previous months and to compare total excess mortality with deaths among positive cases of SARS-CoV-2. Data were stratified by geographical area (north vs centre and south), age and sex.
COVID-19 had a greater impact in northern Italian cities among subjects aged 75-84 and 85+ years. COVID-19 deaths accounted for half of total excess mortality in both areas, with differences by age: almost all excess deaths were from COVID-19 among adults, while among the elderly only one third of the excess was coded as COVID-19. When taking into account the mortality deficit in the pre-pandemic period, different trends were observed by area: all excess mortality during COVID-19 was explained by deficit mortality in the centre and south, while only a 16% overlap was estimated in northern cities, with quotas decreasing by age, from 67% in the 15-64 years old to 1% only among subjects 85+ years old.
An underestimation of COVID-19 deaths is particularly evident among the elderly. When quantifying the burden in mortality related to COVID-19, it is important to consider seasonal dynamics in mortality. Surveillance data provides an impartial indicator for monitoring the following phases of the epidemic, and may help in the evaluation of mitigation measures adopted.
AbstractObjectiveTo assess short term mortality risks and excess mortality associated with exposure to ozone in several cities worldwide.DesignTwo stage time series analysis.Setting406 cities in 20 ...countries, with overlapping periods between 1985 and 2015, collected from the database of Multi-City Multi-Country Collaborative Research Network.PopulationDeaths for all causes or for external causes only registered in each city within the study period.Main outcome measuresDaily total mortality (all or non-external causes only).ResultsA total of 45 165 171 deaths were analysed in the 406 cities. On average, a 10 µg/m3 increase in ozone during the current and previous day was associated with an overall relative risk of mortality of 1.0018 (95% confidence interval 1.0012 to 1.0024). Some heterogeneity was found across countries, with estimates ranging from greater than 1.0020 in the United Kingdom, South Africa, Estonia, and Canada to less than 1.0008 in Mexico and Spain. Short term excess mortality in association with exposure to ozone higher than maximum background levels (70 µg/m3) was 0.26% (95% confidence interval 0.24% to 0.28%), corresponding to 8203 annual excess deaths (95% confidence interval 3525 to 12 840) across the 406 cities studied. The excess remained at 0.20% (0.18% to 0.22%) when restricting to days above the WHO guideline (100 µg/m3), corresponding to 6262 annual excess deaths (1413 to 11 065). Above more lenient thresholds for air quality standards in Europe, America, and China, excess mortality was 0.14%, 0.09%, and 0.05%, respectively.ConclusionsResults suggest that ozone related mortality could be potentially reduced under stricter air quality standards. These findings have relevance for the implementation of efficient clean air interventions and mitigation strategies designed within national and international climate policies.
Display omitted
•Desert dust episodes were more intense and frequent in Ahvaz than in Tehran.•The effect of PM10 in mortality was higher for dust than non-dust days in Ahvaz.•In Tehran the PM10 ...effect was slightly higher during non-dust days.•Effects were evidenced for PM2.5 only in Tehran for non-dust days at lags 2 and 3.•Middle East desert dust is an important risk factor to human health in this region.
Increased atmospheric particulate matter (PM) concentrations are commonly observed during desert dust days in Iran, but there is still no evidence of their effects on human health. We aimed to evaluate the association between daily mortality and exposure to PM10 and PM2.5 during dust and non-dust days in Tehran and Ahvaz, two major Middle Eastern cities with different sources, intensity, and frequency of desert dust days.
We identified desert dust days based on exceeding a daily PM10 concentration threshold of 150 µg/m3 between 2014 and 2017, checking for low PM2.5/PM10 ratio typical of dust days. We used a time-stratified case-crossover design to estimate the short-term effects of PM10 and PM2.5 concentrations on daily mortality during dust and non-dust days. Data was analyzed using conditional Poisson regression models.
Higher concentrations of PM and frequency of desert dust days were observed in Ahvaz rather than Tehran. In Ahvaz, the effect of PM10 at lag 0 was much higher during dust days, an increment of 10 μg/m3 was associated with 3.28% (95%CI = 2.42, 4.15) increase of daily mortality, than non-dust days, 1.03% (95%CI = −0.02, 2.08), while in Tehran, was slightly higher during non-dust days, 0.72% (95%CI = 0.23, 1.23), than in dust days, 0.49% (95%CI = −0.22, 1.20). No statistically significant associations were observed between PM2.5 and daily mortality in Ahvaz, while in Teheran the effect of PM2.5 increased significantly during non-dust days at lag 2, 1.89% (95%CI = 0.83, 1.2.95 and lag 3, 1.88% (95%CI = 0.83, 1.2.95).
The study provides evidence that exposure to PM during Middle East dust days is an important risk factor to human health in arid regions and areas affected by desert dust events.