What this study adds
Recent studies reveal that temperature may modify the ozone-related mortality risk. However, previous results are inconsistent and lack consensus on temperature-confounding control approaches in addition to including an interaction term. This study in eight European cities and 86 US cities found that inadequate control of temperature confounding by only including a categorical temperature term in time series models leads to residual confounding and an overestimation of the ozone-related mortality risk at high temperatures. Furthermore, this residual confounding and overestimation at high temperatures are mainly dominated by the incomplete control of heat effects rather than cold effects.
Epidemiological studies have demonstrated the adverse mortality effects of short-term ozone exposure, while temperature is generally conceptualized as a confounder.1–4 Recent studies revealed that temperature may also modify the ozone-related mortality risk.5 However, the confounding effect of temperature in previous studies has been controlled differently: some specified temperature only as a categorical variable,6,7 whereas some used both categorized temperature and smooth terms of continuous temperature to adjust for its nonlinear health effects.8–12 The rational for including only a categorical temperature term may be because of the scheme for identifying effect modification in epidemiology, where unadjusted strata-specific risk estimates are used to represent the modification effect if significant heterogeneity exists among the strata-specific risk estimates.13 However, the temperature–mortality relationship has been demonstrated to be nonlinear with increasing mortality risks for both high temperatures (heat) and low temperatures (cold).14,15 Thus, including only a categorical variable for temperature may not be adequate to fully rule out its nonlinear confounding effect and may lead to potential residual confounding.
Furthermore, the nonlinear temperature–mortality relationship has also been found to differ in lag structures between heat and cold effects. Heat effects are generally acute and occur within a few days, while cold effects are delayed and last up to 3 or 4 weeks.16–18 Previous studies tended to adjust only for temperature with a short lag (concurrent day or up to 3 days),6–10,12,19 which may lead to inadequate control for cold effects. When using short lags for temperature in the ozone–mortality association, inadequate control of cold effects could result in higher ozone effect estimates in winter than in summer in 21 East Asian cities.20
In this study, we aimed to evaluate the influence of varying adjustment for confounding by temperature on the modifying effect of temperature and potential residual confounding of the ozone-related mortality in eight European cities and 86 US cities, respectively.
We collected daily mortality, air temperature, and ozone data from 94 cities from two continents: eight European cities during 1999–2013 and 86 US cities during 1987–2000 (see eFigure 1; http://links.lww.com/EE/A3 for locations and eTable 1; http://links.lww.com/EE/A3 for summary statistics of study cities). These data have been previously used in the Ultrafine Particles and Health (UF&HEALTH) Study in Europe21 and the National Morbidity, Mortality, and Air Pollution Study (NMMAPS) in the United States.22,23 Details with regard to data collection can be found in the preceding studies.1,21,24
Daily counts from nonaccidental (International Classification of Diseases 9th revision 1–799 and International Statistical Classification of Diseases and Related Health Problems 10th Revision A00-R99) deaths for all ages were used to represent daily mortality. Consistent with the preceding studies, we used daily maximum 8-hour average of ozone in the European data21 and daily 24-hour average of ozone in the US data.1 In the European data, daily concentration in each city was calculated by averaging measurements from multiple local monitors with at least 75% daily data. In the US data, a 10% trimmed mean of air pollutant was used to protect against outliers and to average across monitors after correction for each monitor’s yearly average.25 We further collected daily concentration of particulate matter with aerodynamic diameter <10 μm (PM10) and nitrogen dioxide (NO2). Among the 86 US cities, 61 cities had at least 2 years of daily measurements for both PM10 and NO2. Data on daily mean temperature were obtained from local or airport meteorological stations.
We performed generalized additive models (GAMs) with an overdispersed Poisson distribution to evaluate the association between ozone and mortality. Consistent with previous studies,3,21 we used previous day’s ozone in Europe and the average of the same and previous day’s ozone in the United States while adjusting for long-term and seasonal trend and day of the week. Univariate penalized cubic regression splines were used to model the smooth terms of time trend. As in the preceding studies,7,21 6 degrees of freedom (df) per year for time trend was selected. We explored the sensitivity of the smoothness of time trend when using 6, 8, and 10 df per year by comparing generalized cross validation scores and mortality effect estimates.
To visually capture the potential interaction, we first applied the GAMs to fit a bivariate response surface.7,26 The ozone–temperature response surface was performed with a tensor product smooth, which was constructed using cubic regression splines with df chosen by a generalized cross-validation method. This data-driven approach has been used to illustrate the joint effects of ozone and temperature on daily mortality.8,10,27
To examine the effect modification by temperature, we categorized temperature into three levels: high (>75th city-specific percentile), medium (25th–75th city-specific percentile), and low (<25th city-specific percentile) in each city. In the GAMs, we included an interaction term between ozone and categorized temperature at the same lag structure. We first estimated the ozone-related mortality risks within each temperature category in each city and then pooled the city-level estimates using random-effects meta-analyses.28 We tested the statistical significance of differences between the temperature-specific estimates by calculating the 95% confidence interval (CI) according to Zeka et al.29 In sensitivity analysis, we also applied a different temperature stratification method by multiplying three binary temperature stratum indicators with ozone to allow ozone to have a different mortality effect within each temperature stratum.26
To investigate the influence of temperature-confounding control on temperature-stratified ozone-related mortality risks, we compared adjustments for temperature confounding by using only a categorical variable for temperature (i.e., the main effect of temperature strata in the interaction term) using additionally a continuous variable for temperature. In modeling the continuous temperature, we used two methods to account for differences in the lag structures of heat and cold effects: a distributed lag nonlinear model (DLNM) approach and an approach fitting separate terms for high and low temperatures. To explore which confounding (heat or cold) affects the temperature-stratified ozone–mortality associations, we applied two other approaches by adjusting only the heat or cold effect. Thus, we used a total of five approaches to control for confounding by temperature: (1) only the categorical term without adjustment for the smooth terms of temperature (“without sTemp”); (2) in addition to the categorical term, fitting a DLNM using a quadratic B-spline for temperature with three internal knots at 10th, 75th, and 90th percentiles, and a natural cubic B-spline for lag with 21 days and three internal knots at equally spaced log-scale values (“sTemp:DLNM”; see Gasparrini et al.14 for more information); (3) in addition to the categorical term, fitting separate natural spline smooth terms for low temperatures at lag 1–6 and high temperatures at lag 0–1 (“sTemp_heat+sTemp_cold”; see Stafoggia et al.21 for further details); (4) in addition to the categorical term, fitting a nature spline smooth term for high temperatures at lag 0–1 (“sTemp_heat”); and (5) in addition to the categorical term, fitting a nature spline smooth term for cold temperatures at lag 1–6 (“sTemp_cold”). The Akaike Information Criterion for quasi-Poisson was applied to compare the model fits among five approaches.
In addition, we fitted two-pollutant models to adjust for the potential confounding effect of co-pollutants (i.e., PM10 and NO2) using the same lag structure as ozone. As there were missing values for ozone and co-pollutants, we restricted the analyses of two-pollutant models to days with no missing values for either air pollutant.
To evaluate potential residual confounding in these five temperature control approaches, we applied a method proposed by Flanders et al.30: we added future ozone concentrations measured one day after death occurrence and tested its statistical significance. If future ozone concentrations are significantly associated with daily mortality, residual confounding or other potential bias is suggested. As Becher31 suggested that splitting the continuous confounding factor into four or five levels would remove most of the residual confounding, we also used a larger number of categories for the categorical temperature term (five and eight) defined by city-specific percentiles.
All analyses were performed with R software, version 3.2.1 (R Foundation for Statistical Computing, Vienna, Austria), using the packages mgcv,32 dlnm,33 and mvmeta.28
In the response surface model, the tensor product smooth terms of ozone and temperature were statistically significant in seven of the eight European cities and 29 of the 86 US cities, implying that an interaction between ozone and temperature might exist. Figure 1 graphically depicts the potential interactive effects of ozone and temperature on daily mortality in eight European cities and eight large US cities. The response surfaces show that temperature positively modify the O3 mortality effects in most European and US cities, but negative or no modification by temperature were also observed in a few cities such as Augsburg, Copenhagen, and Atlanta.
Figure 2 shows the pooled results of the temperature-modified ozone-related mortality risks using five temperature control methods for the European and US cities. If the nonlinear effect of temperature was not further controlled by smooth terms (without sTemp), the categorical temperature variable showed a significant effect modification of the ozone-related mortality risk for high temperatures in both the European and the US cities. However, ozone-related mortality risk estimates at high temperatures were significantly lower when adjusting for the nonlinear effect of temperature using approaches sTemp:DLNM, sTemp_heat+sTemp_cold, or sTemp_heat but remained similar when adjusting only for the cold effect (sTemp_cold). For example, when using the sTemp:DLNM approach, risk estimates at low, medium, and high temperatures were 0.23% (95% CI = −0.09%, 0.55%), 0.23% (95% CI = −0.06%, 0.53%), and 0.36% (95% CI = 0.04%, 0.68%) in European cities and 0.11% (95% CI = −0.31%, 0.54%), 0.17% (95% CI = −0.07%, 0.41%), and 0.59% (95% CI = 0.32%, 0.85%) in US cities. When using sTemp_heat+sTemp_cold or sTemp_heat approaches, the higher risk estimates at high temperatures relative to those at low or medium temperatures became nonsignificant for European cities but generally remained significant for US cities. Compared with without sTemp approach, other approaches controlling the smooth terms of temperature generally yielded much better model fits with lower Akaike Information Criterion for quasi-Poisson values (eTable 2; http://links.lww.com/EE/A3).
After adjustment for co-pollutants, the estimates of temperature-stratified ozone-related mortality risk remained similar for both European and US cities (Figure 3). Adjustment for PM10 or NO2 did not change the pattern among five temperature control methods. However, the confidence intervals in US cities became larger because the sample size decreased because of missing data for co-pollutants and fewer cities with available data.
Table 1 summarizes the regional acute effects of tomorrow’s ozone concentrations on daily mortality for the five temperature adjustment methods. In both regions, we observed significant associations between tomorrow’s ozone concentrations and current day’s mortality with adjustment approaches without sTemp and sTemp_cold, which implies residual confounding. Splitting of temperature into more categories (five or eight) yielded similar results as using three categories.
Our results remained robust when using 8 or 10 df per year for the time trend adjustment (eFigure 2; http://links.lww.com/EE/A3, eTable 3 and eTable 4; http://links.lww.com/EE/A3), when using different cutoff values for the temperature categories (20th/80th, 15th/85th, and 10th/90th; eFigure 3; http://links.lww.com/EE/A3 and eTable 5; http://links.lww.com/EE/A3), or when using an alternative temperature stratification method by fitting ozone effects within each temperature stratum (eFigure 4; http://links.lww.com/EE/A3).
Our findings show that inadequate control of temperature confounding by only including a categorical temperature term in time-series models leads to residual confounding and an overestimation of the ozone-related mortality risk at high temperatures. This residual confounding and overestimation at high temperatures are mainly dominated by the incomplete control of heat effects rather than cold effects.
Few time-series studies have investigated the modifying effect of the ozone-related mortality risk by temperature, and the results were inconsistent: some found a significant enhancement only for high temperatures,7,11,12,27,34 whereas others only for low temperatures,8–10 and some found higher ozone effects for both low and high temperatures.6 These inconsistencies may be attributable to different climatic and demographic conditions, population exposure patterns, air pollution characteristics of different countries, and different modeling and covariate control approaches.5,8 Our results suggest that inadequate control of temperature confounding may also contribute to these inconsistencies. In this study, the significant effect of ozone on mortality at high temperatures was largely reduced when including a nonlinear term for temperature effects, mainly attributable to the effect of heat (Figure 2). In previous studies, Ren et al.7 only controlled for temperature with a categorical term while Pascal et al.34 did not control for the heat effect, which may partially lead to the reported significant modifying effects of high temperatures. Similar to our sTemp_heat approach, Pattenden et al.11 adjusted for the nonlinear heat effect using a cubic spline of the 2-day average daily mean temperature and found an overall weak evidence of the ozone–heat interaction in 15 British conurbation.
We found that when the confounding effect of heat exposure was not fully removed, ozone-related mortality risks at high temperatures would be greatly overestimated. This bias remained after controlling for PM10 or NO2 in the two-pollutant models, suggesting that the residual confounding was not due to co-pollutants. This could be because air pollution were unlikely to confound studies of temperature effects35; thus, controlling for co-pollutants would not remove the residual confounding by temperature.
In this study, we used a city-specific percentile-based method to categorize temperature. This method used relative cutoff values to characterize differences in temperature distributions and population acclimatization to temperature changes in cities with different climate conditions. The city-specific temperature percentiles corresponded to different city-specific absolute temperatures. Even though temperature distributions differed across cities, the use of temperature percentiles rather than absolute values has been shown to yield greater consistency in observed mortality effects across cities.16 The minimum mortality temperature, at which mortality was the lowest, was found to be consistently around the 75th percentile across various climate patterns (from tropical to temperate), indicating a long-term population adaptation to their local climates.17 Thus, using the percentile-based method may be more appropriate to pool the temperature-stratified ozone effect estimates at each region.
In studies examining the effect modification by temperature on air pollution-related health impacts, bivariate response surface and stratification models were widely used.7–10,26 Bivariate response surface models are flexible without linear assumptions on association between independent and dependent variables, but is limited to visual inspection because of unavailable parametric estimates.7,19 In addition, response surface models also presented a challenge in including an additional lag dimension, which is critical in assessing interactions between temperature and air pollution with different lag structures. To quantify the effect modification by temperature, two types of temperature stratification models were applied: one with an interaction term between air pollutant and temperature levels7,8 and another one without including an interaction term and directly estimated the air pollution effects across temperature stratum.9,26 In our analyses, using these two stratification models generally yielded similar temperature-stratified ozone mortality risks and came to the same conclusion. However, using the later temperature stratification model without an interaction term resulted in negative ozone mortality effect estimates at low temperatures in US cities when cold effects were not controlled (eFigure 4; http://links.lww.com/EE/A3). Similarly, in a previous study using the temperature stratification model without an interaction term, negative PM10 mortality effects at low temperatures were observed when adjusting only for temperature at short lags (current day and average of previous 3 days).26
This study has some limitations. First, different metrics and units of ozone were used in European and US cities because the same metrics using the same units were not available. However, this difference is unlikely to influence our overall findings, as our focus was to evaluate the impact of temperature-confounding adjustment on the modifying effect of temperature in the ozone–mortality relationship, rather than to compare the difference of the ozone–mortality relationship between Europe and the United States. Moreover, while ozone might act as an intermediator in the heat–mortality relationship,36 mediation analysis is beyond the scope of our article.
In conclusion, our findings suggest that residual confounding of temperature can lead to an overestimation of its modifying effect on ozone–mortality associations. Future research investigating the modifying effect of temperature in the relationship between short-term air pollution exposure and adverse health outcomes shall be recommended to adequately adjust for the confounding by temperature.
Conflicts of interest statement
The authors declare that they have no conflicts of interest with regard to the content of this report.
K.C. was supported by the Alexander von Humboldt Foundation (Humboldt Research Fellowships for postdoctoral researchers).
The computing code is available upon request from the corresponding author. The data are not available for replication because the mortality data are confidential.
We thank the Instituto Nacional de Estadística and the Agència de Salut Pública de Barcelona for providing the mortality data and the Agencia Estatal de Meteorologia (Ministerio de Agricultura, Alimentación y Medio Ambiente) for providing the weather data for Spain. We thank the Institute of Environmental Assessment and Water Research (IDAEA-CSIC, Barcelona, Spain). Colleagues from IDAEA-CSIC were supported by the project PI15/00515, integrated in the National Plan for I+D+I and cofunded by the ISCIII-Directorate General for Evaluation and the European Regional Development Fund (FEDER). We thank “Information und Technik NRW, Düsseldorf, 2014” and “Landesamt für Natur, Umwelt und Verbraucherschutz Land NRW, Recklinghausen, www.lanuv.nrw” for providing, respectively, mortality and exposure data for the three cities of the Ruhr Area. We thank Dr. H. Ott from the Bavarian Environmental Agency (Bayerisches Landesamt für Umwelt) for providing the air pollution and meteorological data from Augsburg, Germany. We thank Helsinki Region Environmental Services Authority HSY for providing the air pollution data from Helsinki, Finland. We also thank Finnish Meteorological Institute for providing the weather data and Statistics Finland for providing the mortality data for Finland. The study has been conducted as a collaborative effort of the UF&HEALTH Study Group.
UF&HEALTH Study Group: S. Breitner, J. Cyrys, R. Hampel, F. Hennig, B. Hoffmann, T. Kuhlbusch; S. Lanzinger, A. Peters, U. Quass, A. Schneider, K. Wolf (Germany); E. Diapouli, K. Elefteriadis, K. Katsouyanni, E. Samoli, S. Vratolis (Greece); T. Ellermann, Z. Ivanovic-Andersen, S. Loft, A. Massling, C. Nordstrøm (Denmark); P. P. Aalto, M. Kulmala, T. Lanki, J. Pekkanen, P. Tiittanen, T. Yli-Tuomi (Finland); G. Cattani, A. Faustini, F. Forastiere, M. Inglessis, M. Renzi, M. Stafoggia (Italy); D. Agis, X. Basagaña, B. Jacquemin, N. Perez, J. Sunyer, A. Tobias (Spain); G. Bero-Bedada, T. Bellander (Sweden).
1. Bell ML, McDermott A, Zeger SL, Samet JM, Dominici F. Ozone
and short-term mortality
in 95 US urban communities, 1987–2000. JAMA 2004; 292:23722378.
2. Ito K, De Leon SF, Lippmann M. Associations between ozone
and daily mortality
: analysis and meta-analysis. Epidemiology 2005; 16:446457.
3. Peng RD, Samoli E, Pham L, Dominici F, Touloumi G, Ramsay T, et al. Acute effects of ambient ozone
in Europe and North America: results from the APHENA study. Air Qual Atmos Health 2013; 6:445453.
4. Wong C-M, Vichit-Vadakan N, Kan H, Qian Z. Public Health and Air Pollution in Asia (PAPA): a multicity study of short-term effects of air pollution on mortality
. Environ Health Perspect 2008; 116:11951202.
5. Li J, Woodward A, Hou X-Y, Zhu T, Zhang J, Brown H, et al. Modification of the effects of air pollutants on mortality
by temperature: a systematic review and meta-analysis. Sci Total Environ 2017; 575:15561570.
6. Jhun I, Fann N, Zanobetti A, Hubbell B. Effect modification of ozone
risks by temperature in 97 US cities. Environ Int 2014; 73:128134.
7. Ren C, Williams GM, Mengersen K, Morawska L, Tong S. Does temperature modify short-term effects of ozone
on total mortality
in 60 large eastern US communities?—An assessment using the NMMAPS data. Environ Int 2008; 34:451458.
8. Chen K, Yang HB, Ma ZW, Bi J, Huang L. Influence of temperature to the short-term effects of various ozone
metrics on daily mortality
in Suzhou, China. Atmos Environ 2013; 79:119128.
9. Cheng Y, Kan H. Effect of the interaction between outdoor air pollution and extreme temperature on daily mortality
in Shanghai, China. J Epidemiol 2012; 22:2836.
10. Liu T, Li TT, Zhang YH, Xu YJ, Lao XQ, Rutherford S, et al. The short-term effect of ambient ozone
is modified by temperature in Guangzhou, China. Atmos Environ 2013; 76:5967.
11. Pattenden S, Armstrong B, Milojevic A, Heal MR, Chalabi Z, Doherty R, et al. Ozone
, heat and mortality
: acute effects in 15 British conurbations. Occup Environ Med 2010; 67:699707.
12. Qian Z, He Q, Lin H-M, Kong L, Bentley CM, Liu W, et al. High temperatures enhanced acute mortality
effects of ambient particle pollution in the “oven” city of Wuhan, China. Environ Health Perspect 2008; 116:1172.
13. Webb P, Bain C. Essential Epidemiology: An Introduction for Students and Health Professionals 20112nd edNew York: Cambridge University Press.
14. Gasparrini A, Guo YM, Hashizume M, Lavigne E, Zanobetti A, Schwartz J, et al. Mortality
risk attributable to high and low ambient temperature: a multicountry observational study. Lancet 2015; 386:369375.
15. Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz JA. Temperature and mortality
in 11 cities of the Eastern United States. Am J Epidemiol 2002; 155:8087.
16. Anderson BG, Bell ML. Weather-related mortality
: how heat, cold, and heat waves affect mortality
in the United States. Epidemiology 2009; 20:205213.
17. Guo Y, Gasparrini A, Armstrong B, Li S, Tawatsupa B, Tobias A, et al. Global variation in the effects of ambient temperature on mortality
: a systematic evaluation. Epidemiology 2014; 25:781789.
18. Analitis A, Katsouyanni K, Biggeri A, Baccini M, Forsberg B, Bisanti L, et al. Effects of cold weather on mortality
: results from 15 European cities within the PHEWE Project. Am J Epidemiol 2008; 168:13971408.
19. Wilson A, Rappold AG, Neas LM, Reich BJ. Modeling the effect of temperature on ozone
. Ann Appl Stat 2014; 8:17281749.
20. Chen R, Cai J, Meng X, Kim H, Honda Y, Guo YL, et al. Ozone
and daily mortality
rate in 21 cities of East Asia: how does season modify the association? Am J Epidemiol 2014; 180:729736.
21. Stafoggia M, Schneider A, Cyrys J, Samoli E, Andersen ZJ, Bedada GB, et al. Association between short-term exposure to ultrafine particles and mortality
in eight European urban areas. Epidemiology 2017; 28:172180.
22. Peng RD, Dominici F. Statistical Methods for Environmental Epidemiology With R. R: A Case Study in Air Pollution and Health 2008; New York: Springer.
23. Peng RD, Welty LJ. The NMMAPSdata package. R news 2004; 4:1014.
24. Samet J, Zeger S, Dominici F, Curriero F, Coursac I, Dockery D, et al. The National Morbidity, Mortality
, and Air Pollution Study. Part II: morbidity and mortality
from air pollution in the United States. Res Rep Health Eff Inst 2000; 94:570.
25. Samet J, Dominici F, Zeger S, Schwartz J, Dockery D. The National Morbidity, Mortality
, and Air Pollution Study. Part I: methods and methodologic issues. Res Rep Health Eff Inst 2000; 94:514.
26. Roberts S. Interactions between particulate air pollution and temperature in air pollution mortality
time series studies. Environ Res 2004; 96:328337.
27. Burkart K, Canário P, Breitner S, Schneider A, Scherber K, Andrade H, et al. Interactive short-term effects of equivalent temperature and air pollution on human mortality
in Berlin and Lisbon. Environ Pollut 2013; 183:5463.
28. Gasparrini A, Armstrong B, Kenward MG. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat Med 2012; 31:38213839.
29. Zeka A, Zanobetti A, Schwartz J. Individual-level modifiers of the effects of particulate matter on daily mortality
. Am J Epidemiol 2006; 163:849859.
30. Flanders WD, Klein M, Darrow LA, Strickland MJ, Sarnat SE, Sarnat JA, et al. A method for detection of residual confounding in time-series and other observational studies. Epidemiology 2011; 22:5967.
31. Becher H. The concept of residual confounding in regression models and some applications. Stat Med 1992; 11:17471758.
32. Wood SN. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J R Stat Soc Ser C Appl Stat 2011; 73:336.
33. Gasparrini A. Distributed lag linear and non-linear models in R: the Package dlnm. J Stat Softw 2011; 43:120.
34. Pascal M, Wagner V, Chatignoux E, Falq G, Corso M, Blanchard M, et al. Ozone
and short-term mortality
in nine French cities: influence of temperature and season. Atmos Environ 2012; 62:566572.
35. Buckley JP, Samet JM, Richardson DB. Commentary: does air pollution confound studies of temperature? Epidemiology 2014; 25:242245.
36. Reid CE, Snowden JM, Kontgis C, Tager IB. The role of ambient ozone
in epidemiologic studies of heat-related mortality
. Environ Health Perspect 2012; 120:1627.