Comparison of spatiotemporal prediction models of daily exposure of individuals to ambient nitrogen dioxide and ozone in Montreal, Canada

Spatio-temporal evaluation of air pollution using ground-based and satellite data during COVID-19 in Ecuador

The concentration of gases in the atmosphere is a topic of growing concern due to its effects on health, ecosystems etc. Its monitoring is commonly carried out through ground stations which offer high precision and temporal resolution. However, in countries with few stations, such as Ecuador, these data fail to adequately describe the spatial variability of pollutant concentrations. Remote sensing data have great potential to solve this complication. This study evaluates the spatiotemporal distribution of nitrogen dioxide (NO2) and ozone (O3) concentrations in Quito and Cuenca, using data obtained from ground-based and Sentinel-5 Precursor mission sources during the years 2019 and 2020. Moreover, a Linear Regression Model (LRM) was employed to analyze the correlation between ground-based and satellite datasets, revealing positive associations for O3 (R2 = 0.83, RMSE = 0.18) and NO2 (R2 = 0.83, RMSE = 0.25) in Quito; and O3 (R2 = 0.74, RMSE = 0.23) and NO2, (R2 = 0.73, RMSE = 0.23) for Cuenca. The agreement between ground-based and satellite datasets was analyzed by employing the intra-class correlation coefficient (ICC), reflecting good agreement between them (ICC ≥0.57); and using Bland and Altman coefficients, which showed low bias and that more than 95% of the differences are within the limits of agreement. Furthermore, the study investigated the impact of COVID-19 pandemic-related restrictions, such as social distancing and isolation, on atmospheric conditions. This was categorized into three periods for 2019 and 2020: before (from January 1st to March 15th), during (from March 16th to May 17th), and after (from March 18th to December 31st). A 51% decrease in NO2 concentrations was recorded for Cuenca, while Quito experienced a 14.7% decrease. The tropospheric column decreased by 27.3% in Cuenca and 15.1% in Quito. O3 showed an increasing trend, with tropospheric concentrations rising by 0.42% and 0.11% for Cuenca and Quito respectively, while the concentration in Cuenca decreased by 14.4%. Quito experienced an increase of 10.5%. Finally, the reduction of chemical species in the atmosphere as a consequence of mobility restrictions is highlighted. This study compared satellite and ground station data for NO2 and O3 concentrations. Despite differing units preventing data validation, it verified the Sentinel-5P satellite's effectiveness in anomaly detection. Our research's value lies in its applicability to developing countries, which may lack extensive monitoring networks, demonstrating the potential use of satellite technology in urban planning.

Comparison of PM2.5 Air Pollution Exposures and Health Effects Associations Using 11 Different Modeling Approaches in the Women's Health Initiative Memory Study (WHIMS)

BACKGROUND

Many approaches to quantifying air pollution exposures have been developed. However, the impact of choice of approach on air pollution estimates and health-effects associations remains unclear.

OBJECTIVES

Our objective is to compare particulate matter with aerodynamic diameter ≤ 2.5 μ m (PM 2.5 ) concentrations and resulting health effects associations using multiple estimation approaches previously used in epidemiologic analyses.

METHODS

We assigned annual PM 2.5 exposure estimates from 1999 to 2004 derived from 11 different approaches to Women's Health Initiative Memory Study (WHIMS) participant addresses within the contiguous US. Approaches included geostatistical interpolation approaches, land-use regression or spatiotemporal models, satellite-derived approaches, air dispersion and chemical transport models, and hybrid models. We used descriptive statistics and plots to assess relative and absolute agreement among exposure estimates and examined the impact of approach on associations between PM 2.5 and death due to natural causes, cardiovascular disease (CVD) mortality, and incident CVD events, adjusting for individual-level covariates and climate-based region.

RESULTS

With a few exceptions, relative agreement of approach-specific PM 2.5 exposure estimates was high for PM 2.5 concentrations across the contiguous US. Agreement among approach-specific exposure estimates was stronger near PM 2.5 monitors, in certain regions of the country, and in 2004 vs. 1999. Collectively, our results suggest but do not quantify lower agreement at local spatial scales for PM 2.5 . There was no evidence of large differences in health effects associations with PM 2.5 among estimation approaches in analyses adjusted for climate region.

CONCLUSIONS

Different estimation approaches produced similar spatial patterns of PM 2.5 concentrations across the contiguous US and in areas with dense monitoring data, and PM 2.5 -health effects associations were similar among estimation approaches. PM 2.5 estimates and PM 2.5 -health effects associations may differ more in samples drawn from smaller areas or areas without substantial monitoring data, or in analyses with finer adjustment for participant location. Our results can inform decisions about PM 2.5 estimation approach in epidemiologic studies, as investigators balance concerns about bias, efficiency, and resource allocation. Future work is needed to understand whether these conclusions also apply in the context of other air pollutants of interest. https://doi.org/10.1289/EHP12995.

Immortal Time Bias With Time-Varying Exposures in Environmental Epidemiology: A Case Study in Lung Cancer Survival

Immortal time bias is a well-recognized bias in clinical epidemiology but is rarely discussed in environmental epidemiology. Under the target trial framework, this bias is formally conceptualized as a misalignment between the start of study follow-up (time 0) and treatment assignment. This misalignment can occur when attained duration of follow-up is encoded into treatment assignment using minimums, maximums, or averages. The bias can be exacerbated in the presence of time trends commonly found in environmental exposures. Using lung cancer cases from the California Cancer Registry (2000-2010) linked with estimated concentrations of particulate matter less than or equal to 2.5 μm in aerodynamic diameter (PM2.5), we replicated previous studies that averaged PM2.5 exposure over follow-up in a time-to-event model. We compared this approach with one that ensures alignment between time 0 and treatment assignment, a discrete-time approach. In the former approach, the estimated overall hazard ratio for a 5-μg/m3 increase in PM2.5 was 1.38 (95% confidence interval: 1.36, 1.40). Under the discrete-time approach, the estimated pooled odds ratio was 0.99 (95% confidence interval: 0.98, 1.00). We conclude that the strong estimated effect in the former approach was likely driven by immortal time bias, due to misalignment at time 0. Our findings highlight the importance of appropriately conceptualizing a time-varying environmental exposure under the target trial framework to avoid introducing preventable systematic errors.

A performance comparison study on PM2.5 prediction at industrial areas using different training algorithms of feedforward-backpropagation neural network (FBNN)

Presence of particulate matters with aerodynamic diameter of less than 2.5 μm (PM_2.5) in the atmosphere is fast increasing in Malaysia due to industrialization and urbanization. Prolonged exposure of PM_2.5 can cause serious health effects to human. This research is aimed to identify the most reliable model to predict the PM_2.5 pollution using multi-layered feedforward-backpropagation neural network (FBNN). Air quality and meteorological data were collected from Department of Environment (DOE) Malaysia. Six different training algorithms consisting of thirteen various training functions were trained and compared. FBNN model with the highest coefficient correlation (R²) and lowest root mean square error (RMSE), mean absolute error (MAE) and mean absolute percentage error (MAPE) were selected as the best performing model. Levenberg Marquardt (trainlm) is the best performing algorithms compared to other algorithms with R² value of 0.9834 and the lowest error values for RMSE (2.3981), MAE (1.7843) and MAPE (0.1063).

Trends and characteristics of ozone and nitrogen dioxide related health impacts in Chinese cities

Ambient ozone pollution has been becoming severe and attributed to considerable health impacts in China. Nitrogen dioxide (NO₂) is involved in atmospheric ozone production while also affecting public health directly. Joint control ozone and NO₂ pollution would be of significance. This study quantitatively assessed the health impact attributed to ambient ozone and NO₂ pollution in 338 Chinese cities from 2015 to 2020. The results reveal the generally opposite trends of ozone- and NO₂-related health impacts in China. From 2015-2020, respiratory and chronic obstructive pulmonary disease (COPD) health impacts attributed to ozone in 338 cities increased by 65.30% and 63.98%. The NO₂-attributed health impacts decreased by 24.80% and 24.62%. In 2020, the ozone- and NO₂-related respiratory health impacts were 3.96 million DALYs (disability-adjusted life years) and 1.47 million DALYs. High health impacts are concentrated in big cities and city clusters. In 2020, the sum of ozone- and NO₂-related respiratory health impacts in the top 20 cities was 0.98 million DALYs and 0.44 million DALYs, accounting for 24.70% and 30.24% of the 338 cities. The population attribution fraction analysis identified the increasing distributional consistency of ozone and NO₂-related health impacts, emphasizing the necessity and possible efficiency of ozone-NO₂ joint control. Emission source analysis based on gridded data provided a reference for understanding health impacts and developing targeted strategies.

Impact of different exposure models and spatial resolution on the long-term effects of air pollution

Long-term exposure to air pollution has been related to mortality in several epidemiological studies. The investigations have assessed exposure using various methods achieving different accuracy in predicting air pollutants concentrations. The comparison of the health effects estimates are therefore challenging. This paper aims to compare the effect estimates of the long-term effects of air pollutants (particulate matter with aerodynamic diameter less than 10 μm, PM₁₀, and nitrogen dioxide, NO₂) on cause-specific mortality in the Rome Longitudinal Study, using exposure estimates obtained with different models and spatial resolutions. Annual averages of NO₂ and PM₁₀ were estimated for the year 2015 in a large portion of the Rome urban area (12 × 12 km²) applying three modelling techniques available at increasing spatial resolution: 1) a chemical transport model (CTM) at 1km resolution; 2) a land-use random forest (LURF) approach at 200m resolution; 3) a micro-scale Lagrangian particle dispersion model (PMSS) taking into account the effect of buildings structure at 4 m resolution with results post processed at different buffer sizes (12, 24, 52, 100 and 200 m). All the exposures were assigned at the residential addresses of 482,259 citizens of Rome 30+ years of age who were enrolled on 2001 and followed-up till 2015. The association between annual exposures and natural-cause, cardiovascular (CVD) and respiratory (RESP) mortality were estimated using Cox proportional hazards models adjusted for individual and area-level confounders. We found different distributions of both NO₂ and PM₁₀ concentrations, across models and spatial resolutions. Natural cause and CVD mortality outcomes were all positively associated with NO₂ and PM₁₀ regardless of the model and spatial resolution when using a relative scale of the exposure such as the interquartile range (IQR): adjusted Hazard Ratios (HR), and 95% confidence intervals (CI), of natural cause mortality, per IQR increments in the two pollutants, ranged between 1.012 (1.004, 1.021) and 1.018 (1.007, 1.028) for the different NO₂ estimates, and between 1.010 (1.000, 1.020) and 1.020 (1.008, 1.031) for PM₁₀, with a tendency of larger effect for lower resolution exposures. The latter was even stronger when a fixed value of 10 μg/m³ is used to calculate HRs. Long-term effects of air pollution on mortality in Rome were consistent across different models for exposure assessment, and different spatial resolutions.

Ambient Fine Particulate Matter and Preterm Birth in California: Identification of Critical Exposure Windows

Exposure to ambient fine particulate matter (particulate matter ≤2.5 μm in aerodynamic diameter (PM2.5)) during pregnancy is associated with preterm birth (PTB), a leading cause of infant morbidity and mortality. Results from studies attempting to identify etiologically relevant exposure periods of vulnerability have been inconsistent, possibly because of failure to consider the time-to-event nature of the outcome and lagged exposure effects of PM2.5. In this study, we aimed to identify critical exposure windows for weekly PM2.5 exposure and PTB in California using California birth cohort data from 2005-2010. Associations were assessed using distributed-lag Cox proportional hazards models. We assessed effect-measure modification by race/ethnicity by calculating the weekly relative excess risk due to interaction. For a 10-μg/m3 increase in PM2.5 exposure over the entire period of gestation, PTB risk increased by 11% (hazard ratio = 1.11, 95% confidence interval: 1.09, 1.14). Gestational weeks 17-24 and 36 were associated with increased vulnerability to PM2.5 exposure. We find that non-Hispanic black mothers may be more susceptible to effects of PM2.5 exposure than non-Hispanic white mothers, particularly at the end of pregnancy. These findings extend our knowledge about the existence of specific exposure periods during pregnancy that have the greatest impact on preterm birth.

Comparison of model estimates from an intra-city land use regression model with a national satellite-LUR and a regional Bayesian Maximum Entropy model, in estimating NO2 for a birth cohort in Sydney, Australia

BACKGROUND

Methods for estimating air pollutant exposures for epidemiological studies are becoming more complex in an effort to minimise exposure error and its associated bias. While land use regression (LUR) modelling is now an established method, there has been little comparison between LUR and other recent, more complex estimation methods. Our aim was to develop a LUR model to estimate intra-city exposures to nitrogen dioxide (NO₂) for a Sydney cohort, and to compare those with estimates from a national satellite-based LUR model (Sat-LUR) and a regional Bayesian Maximum Entropy (BME) model.

METHODS

Satellite-based LUR and BME estimates were obtained using existing models. We used methods consistent with the European Study of Cohorts for Air Pollution Effects (ESCAPE) methodology to develop LUR models for NO₂ and NOx. We deployed 46 Ogawa passive samplers across western Sydney during 2013/2014 and acquired data on land use, population density, and traffic volumes for the study area. Annual average NO₂ concentrations for 2013 were estimated for 947 addresses in the study area using the three models: standard LUR, Sat-LUR and a BME model. Agreement between the estimates from the three models was assessed using interclass correlation coefficient (ICC), Bland-Altman methods and correlation analysis (CC).

RESULTS

The NO₂ LUR model predicted 84% of spatial variability in annual mean NO₂ (RMSE: 1.2 ppb; cross-validated R²: 0.82) with predictors of major roads, population and dwelling density, heavy traffic and commercial land use. A separate model was developed that captured 92% of variability in NOx (RMSE 2.3 ppb; cross-validated R²: 0.90). The annual average NO₂ concentrations were 7.31 ppb (SD: 1.91), 7.01 ppb (SD: 1.92) and 7.90 ppb (SD: 1.85), for the LUR, Sat-LUR and BME models respectively. Comparing the standard LUR with Sat-LUR NO₂ cohort estimates, the mean estimates from the LUR were 4% higher than the Sat-LUR estimates, and the ICC was 0.73. The Pearson's correlation coefficients (CC) for the LUR vs Sat-LUR values were r = 0.73 (log-transformed data) and r = 0.69 (untransformed data). Comparison of the NO₂ cohort estimates from the LUR model with the BME blended model indicated that the LUR mean estimates were 8% lower than the BME estimates. The ICC for the LUR vs BME estimates was 0.73. The CC for the logged LUR vs BME estimates was r = 0.73 and for the unlogged estimates was r = 0.69.

CONCLUSIONS

Our LUR models explained a high degree of spatial variability in annual mean NO₂ and NOx in western Sydney. The results indicate very good agreement between the intra-city LUR, national-scale sat-LUR, and regional BME models for estimating NO₂ for a cohort of children residing in Sydney, despite the different data inputs and differences in spatial scales of the models, providing confidence in their use in epidemiological studies.

Does the air pollution model influence the evidence of socio-economic disparities in exposure and susceptibility?

Studies assessing socio-economic disparities in air pollution exposure and susceptibility are usually based on a single air pollution model. A time stratified case-crossover study was designed to assess the impact of the type of model on differential exposure and on the differential susceptibility in the relationship between ozone exposure and daily mortality by socio-economic strata (SES) in Montreal. Non-accidental deaths along with deaths from cardiovascular and respiratory causes on the island of Montreal for the period 1991-2002 were included as cases. Daily ozone concentration estimates at partictaipants' residence were obtained from the five following air pollution models: Average value (AV), Nearest station model (NS), Inverse-distance weighting interpolation (IDW), Land-use regression model with back-extrapolation (LUR-BE) and Bayesian maximum entropy model combined with a land-use regression (BME-LUR). The prevalence of a low household income (< 20,000/year) was used as socio-economic variable, divided into two categories as a proxy for deprivation. Multivariable conditional logistic regressions were used considering 3-day average concentrations. Multiplicative and additive interactions (using Relative Excess Risk due to Interaction) as well as Cochran's tests were calculated and results were compared across the different air pollution models. Heterogeneity of susceptibility and exposure according to socio-economic status (SES) were found. Ratio of exposure across SES groups means ranged from 0.75 [0.74-0.76] to 1.01 [1.00-1.02], respectively for the LUR-BE and the BME-LUR models. Ratio of mortality odds ratios ranged from 1.01 [0.96-1.05] to 1.02 [0.97-1.08], respectively for the IDW and LUR-BE models. Cochran's test of heterogeneity between the air pollution models showed important heterogeneity regarding the differential exposure by SES, but the air pollution model was not found to influence heterogeneity regarding the differential susceptibility. The study showed air pollution models can influence the assessment of disparities in exposure according to SES in Montreal but not that of disparities in susceptibility.

Associations between ambient air pollution and daily mortality in a cohort of congestive heart failure: Case-crossover and nested case-control analyses using a distributed lag nonlinear model

BACKGROUND

Persons with congestive heart failure may be at higher risk of the acute effects related to daily fluctuations in ambient air pollution. To meet some of the limitations of previous studies using grouped-analysis, we developed a cohort study of persons with congestive heart failure to estimate whether daily non-accidental mortality were associated with spatially-resolved, daily exposures to ambient nitrogen dioxide (NO₂) and ozone (O₃), and whether these associations were modified according to a series of indicators potentially reflecting complications or worsening of health.

METHODS

We constructed the cohort from the linkage of administrative health databases. Daily exposure was assigned from different methods we developed previously to predict spatially-resolved, time-dependent concentrations of ambient NO₂ (all year) and O₃ (warm season) at participants' residences. We performed two distinct types of analyses: a case-crossover that contrasts the same person at different times, and a nested case-control that contrasts different persons at similar times. We modelled the effects of air pollution and weather (case-crossover only) on mortality using distributed lag nonlinear models over lags 0 to 3 days. We developed from administrative health data a series of indicators that may reflect the underlying construct of "declining health", and used interactions between these indicators and the cross-basis function for air pollutant to assess potential effect modification.

RESULTS

The magnitude of the cumulative as well as the lag-specific estimates of association differed in many instances according to the metric of exposure. Using the back-extrapolation method, which is our preferred exposure model, we found for the case-crossover design a cumulative mean percentage changes (MPC) in daily mortality per interquartile increment in NO₂ (8.8 ppb) of 3.0% (95% CI: -0.4, 6.6%) and for O₃ (16.5 ppb) 3.5% (95% CI: -4.5, 12.1). For O₃ there was strong confounding by weather (unadjusted MPC = 7.1%; 95% CI: 1.7, 12.7%). For the nested case-control approach the cumulative MPC for NO₂ in daily mortality was 2.9% (95% CI: -0.9, 6.9%) and for O₃ 7.3% (95% CI: 3.0, 11.9%). We found evidence of effect modification between daily mortality and cumulative NO₂ and O₃ according to the prescribed dose of furosemide in the nested case-control analysis, but not in the case-crossover analysis.

CONCLUSIONS

Mortality in congestive heart failure was associated with exposure to daily ambient NO₂ and O₃ predicted from a back-extrapolation method using a land use regression model from dense sampling surveys. The methods used to assess exposure can have considerable influence on the estimated acute health effects of the two air pollutants.

Blending Multiple Nitrogen Dioxide Data Sources for Neighborhood Estimates of Long-Term Exposure for Health Research

Exposure to traffic related nitrogen dioxide (NO₂) air pollution is associated with adverse health outcomes. Average pollutant concentrations for fixed monitoring sites are often used to estimate exposures for health studies, however these can be imprecise due to difficulty and cost of spatial modeling at the resolution of neighborhoods (e.g., a scale of tens of meters) rather than at a coarse scale (around several kilometers). The objective of this study was to derive improved estimates of neighborhood NO₂ concentrations by blending measurements with modeled predictions in Sydney, Australia (a low pollution environment). We implemented the Bayesian maximum entropy approach to blend data with uncertainty defined using informative priors. We compiled NO₂ data from fixed-site monitors, chemical transport models, and satellite-based land use regression models to estimate neighborhood annual average NO₂. The spatial model produced a posterior probability density function of estimated annual average concentrations that spanned an order of magnitude from 3 to 35 ppb. Validation using independent data showed improvement, with root mean squared error improvement of 6% compared with the land use regression model and 16% over the chemical transport model. These estimates will be used in studies of health effects and should minimize misclassification bias.