Does particulate matter modify the association between temperature and cardiorespiratory diseases?

BACKGROUND: A number of studies have shown that both temperature and air pollution are associated with health outcomes. In assessing air pollution effects, temperature is usually considered a confounder. However, only a few recent studies considered air pollution as confounders while assessing temperature effects. Few studies are available on whether or not air pollution modifies the temperature-disease relationship.

METHODS: In this study, we used three parallel Poisson generalized additive models to examine whether particulate matter < 10 [micro]m in aerodynamic diameter (P[M.sub.10]) modified the effects of minimum temperature on cardiorespiratory morbidity and mortality in Brisbane, Australia.

RESULTS: Results show that P[M.sub.10] statistically significantly modified the effects of temperature on respiratory and cardiovascular hospital admissions, all nonexternal-cause mortality, and cardiovascular mortality at different lags. The enhanced adverse temperature effects were found at higher levels of P[M.sub.10], but no clear evidence emerged for interactive effects on respiratory and cardiovascular emergency visits. Three parallel models produced similar results, which strengthened the validity of findings.

CONCLUSION: We conclude that it is important to evaluate the modification role of air pollution in the assessment of temperature-related health impacts.

KEY WORDS: air pollution, interaction, mortality, particulate matter, temperature. Environ Health Perspect 114:1690-1696 (2006). doi:10.1289/ehp.9266 available via http://dx.doi.org/ [Online 27 July 2006]

**********

The nature and magnitude of the association between temperature and human health has been increasingly recognized (Basu and Samet 2002; Martens 1997; Patz et al. 2000; Samet et al. 1998). Both hyperthermia and hypothermia are generally linked to cardiorespiratory morbidity or mortality (Braga et al. 2001; 2002; Kunst et al. 1993). The patterns of temperature-morbidity/mortality vary across regions, with J-, U-, or V-shapes most commonly observed (Basu and Samet 2002; Braga et al. 2002; Patz et al. 2000). In many regions of the world, death rates in winter are usually higher than those in summer, even though heat waves can cause excess deaths (Braga et al. 2002; McMichael et al. 2001). Seasonal variation in morbidity and mortality may also reflect factors beyond weather, including seasonal patterns of respiratory infections. Consequently, assessments of the effect of weather on human health have usually controlled for seasonality and sometimes for influenza epidemics (Schwartz et al. 2004).

Meanwhile, numerous studies have shown that air pollution is consistently associated with adverse health effect (Bell et al. 2004; Dominici et al. 2006; Samet et al. 2000). However, the role of air pollution is often ignored in assessing the health effects of temperature variability, except in a few recent studies adjusting for air pollution (O'Neill et al. 2003; Rainham and Smoyer-Tomic 2003). None of the previous studies have explored whether exposure to air pollution modifies the association between temperature and health outcomes. If substantial effect modification exists, an inappropriately specified model may result in bias. First, it may be inappropriate to consider air pollution only as a confounder in the assessment of the association between temperature and health outcomes, because air pollution may make people more vulnerable to the effects of temperature variability. Second, some studies have shown that temperature may modify the associations between air pollution and cardiorespiratory diseases (Choi et al. 1997; Katsouyanni et al. 1993; Ren and Tong 2006; Roberts 2004). There is often symmetry in modification--air pollution modifies temperature and then temperature modifies air pollution--but the magnitudes are likely to differ. Finally, the true magnitude of the association between temperature and health outcomes may be obscured if air pollution is an effect modifier of the relationship. In this study we used three parallel time-series models to explore whether particulate matter < 10 [micro]m in aerodynamic diameter (P[M.sub.10]) modified the effects of temperature on cardiorespiratory hospital admissions, emergency visits, and mortality in Brisbane, Australia, during the period 1996-2001.

Materials and Methods

Data collection. The data sets consisted of concurrent daily time series of health outcomes, weather, and air pollution collected in Brisbane City from 1 January 1996, to 31 December 2001. Brisbane City is the capital of Queensland, Australia, with a subtropical climate. In 2001, there were 0.89 million residents in Brisbane City (Brisbane City Council 2006).

Health outcome data in this study were provided by the Queensland Department of Health and comprised cardiovascular hospital admissions (CHA), cardiovascular emergency visits (CEV), cardiovascular mortality (CM), respiratory hospital admissions (RHA), respiratory emergency visits (REV), and all nonexternal-cause mortality (NECM). This analysis excluded respiratory mortality due to limited daily counts (mean, 1.5; range, 0-8). Discharge diagnosis was classified according to International Classification of Diseases, 9th Revision (ICD-9; World Health Organization 1975) (used until July 1999) or 10th Revision (ICD-10; World Health Organization 1992) codes: respiratory diseases (ICD-9: 460-519 or ICD-10: J00-J99), cardiovascular diseases (ICD-9: 390-448 or ICD-10: I00-I79), and external causes (ICD-9: E800-E999 or ICD-10: S00-U99). Influenza (ICD-9: 487.0-487.8 or ICD-10: J10-J11) was excluded from respiratory diseases, but occurrence of influenza outbreak was considered as a potential confounder in the data analysis. All cases were local residents of Brisbane City. The notification of incidence and mortality is a statuary requirement under the Health Act 1937 for all public and private hospitals, nursing homes and pathology services in Queensland. The Queensland Department of Health is responsible for data collection, management, and analysis (Queensland Government 2001).

Daily meteorologic data were supplied by the Australian Bureau of Meteorology (http://www.bom.gov.au/), including daily minimum temperature, relative humidity and rainfall for the period of this study. Air pollution data included ambient 24-hr average concentrations of P[M.sub.10] and ozone. All air pollution data were regularly recorded at a central monitoring site and provided by the Queensland Environmental Protection Agency (http://www.epa.qld.gov.au/).

Data analysis. Poisson generalized additive models (GAM) were employed to explore the associations of temperature and P[M.sub.10] with health outcomes. This assumed that the daily number of counts had an overdispersed Poisson distribution [E([Y.sub.t] = [[mu].sub.t]), var([Y.sub.t]) = [phi][[mu].sub.t]] (Dominici et al. 2004). GAM allows nonparametric smoothing functions to account for potentially nonlinear effects of confounding factors on the dependent variable, such as seasonal variation and weather conditions (Hastie and Tibshirani 1990). We used days of calendar time with a cubic smoothing function to control for the confounding effect of seasonality. We controlled for short-term fluctuation using day of the week as a factor. Other potential confounders, such as relative humidity and influenza outbreaks, were also adjusted for.

Before exploring effect modification of P[M.sub.10] on the temperature-health relationship, we used an independent model to explore the patterns of the relationship between temperature and health outcomes. The independent model is described below (Daniels et al. 2000; Hastie and Tibshirani 1990; Insightful Corporation 2001):

Log[E([Y.sub.t]|X)] = [alpha] + lo (tem[p.sub.t-i], span = 0.25) + lo ([pm.sub.t-i], span = 0.25) + s(seaso[n.sub.t], 7) + [gamma]Do[w.sub.t] + s(yea[r.sub.t], 3) + s(rai[n.sub.t-i], 4) + s(humi[d.sub.t-i], 4) + s(ozon[e.sub.t-i], 4) + [[beta].sub.f]fl[u.sub.t] + [[epsilon].sub.t] = [alpha] + lo(tem[p.sub.t-i], span = 0.25) + lo([pm.sub.t-i], span = 0.25) + COVs + [[epsilon].sub.t], [1]

where t refers to the day of the observation; i refers to lags; E([Y.sub.t]|X) denotes estimated daily case counts on day t; s(dot) and lo(dot) separately denote the cubic smoothing spline and LOESS smooth functions, respectively; [alpha] is the intercept term; tem[p.sub.t-i] is 24-hr minimum temperature on day t-i; [pm.sub.t-i] is P[M.sub.10] on day t-i; seaso[n.sub.t] denotes seasonality using days of calendar time. In accordance with the literature (Daniels et al. 2000), we used 7 degrees of freedom (df) per year for season so that little information from time scales longer than 2 months was included. Do[w.sub.t] is the day of week on day t, and [gamma] is a vector of coefficients. The variables rai[n.sub.t-i], humi[d.sub.t-i], and ozon[e.sub.t-i] refer to rainfall, relative humidity at 0900 hr and ozone on day t-i, respectively; fl[u.sub.t] represents the occurrence of influenza epidemics. Because > 99% of days only have 0 or 1 influenza cases, influenza was categorized as a dummy variable (0 cases, [greater than or equal to] 1 cases on day t). [[beta].sub.f] is the coefficient for influenza; [[epsilon].sub.t] is the residual. COVs represents all other covariates in the model.

Then we used three GAM models to assess whether P[M.sub.10] modified the association of temperature with health outcomes: a nonparametric bivariate response model, a nonstratification model, and a stratification model (Ren and Tong, in press; Roberts 2004). We used a bivariate model to explore visually the combining effects of both temperature and P[M.sub.10] with health outcomes. This was undertaken using a nonparametric smoothing function without linear assumptions that the two predictors linearly depend on outcomes. We used a nonstratification model quantitatively to examine the association of both above predictors with health outcomes with a linear assumption by including an interaction term of temperature and P[M.sub.10] as continuous functions. We used a stratification model quantitatively to assess the associations of temperature with health outcomes across P[M.sub.10] levels by including an interaction term of temperature and P[M.sub.10] in which P[M.sub.10] was categorized into two levels. The three models are described in detail below.

First, we used the nonparametric bivariate response model to identify the joint effects of minimum temperature and P[M.sub.10] on health outcomes. This can capture the relationship between independent and dependent variables without the need for strong assumptions (Hastie and Tibshirani 1990). This model provided a picture of the joint pattern of two predictors (temperature and P[M.sub.10]) on the dependent variable (each of cardiorespiratory morbidities and mortalities). Therefore, it can be used to observe whether or not there is an interactive effect of two continuous predictors on the dependent variable (Greenland 1993; Hastie and Tibshirani 1990). We modified Equation 1 to include a bivariate term for temperature and P[M.sub.10] as follows (Insightful Corporation 2001; Ren and Tong, in press; Roberts 2004):

Log[E([Y.sub.t]|X)] = [alpha] + lo(tem[p.sub.t-i], [pm.sub.t-i], span = 0.25) + COVs + [[epsilon].sub.t], [2]

where lo(tem[p.sub.t-i], [pm.sub.t-i]) means joint effect of temperature and P[M.sub.10] and COVs was the same as model 1.

Second, we used a nonstratification model, assuming a linear relationship, to estimate the interactive effects of P[M.sub.10] and minimum temperature on health outcomes. We added an interaction term to estimate increment in cardiorespiratory mortality/morbidity per unit change in ambient P[M.sub.10] and minimum temperature, as follows:

Log[E([Y.sub.t]|X)] = [alpha] + [[beta].sub.1][pm.sub.t-i] + [[beta].sub.2]tem[p.sub.t-i] + [[beta].sub.3]([pm.sub.t-i] : tem[p.sub.t-i]) + COVs + [[epsilon].sub.t], [3]

where [[beta].sub.1] denotes the increment in mortality/morbidity per unit increase in ambient P[M.sub.10] level, [[beta].sub.2] denotes the increment in mortality/morbidity per unit increase in temperature level, and [[beta].sub.3] estimates the interactive effect of P[M.sub.10] and temperature on health outcomes after adjustment for all other covariates. COVs was the same as in model 1.

Finally, we applied a stratification model to examine whether the effects of temperature on health outcomes were heterogeneous across different levels of P[M.sub.10]. We categorized P[M.sub.10] into two levels (low and high) and then examined whether temperature effects varied across levels of P[M.sub.10]. To assess effect modification in the high end of the temperature range, we used separate data sets to fit this model: one data set with the whole range of temperatures and another database with temperatures [greater than or equal to] 19.3 [degrees]C (75th percentile). We slightly modified Equation 3 as follows:

Log[E([Y.sub.t]|X)] = [alpha] + [[beta].sub.1]tem[p.sub.t-i] + [[beta].sub.2][pm.sub.kt-i] + [[beta].sub.3](tem[p.sub.t-i] : [pm.sub.kt-i]) + COVs + [[epsilon].sub.t], [4]

where [pm.sub.kt] represents levels of P[M.sub.10], tem[p.sub.t-i]:[pm.sub.kt-i] represents the interaction term of temperature and levels of P[M.sub.10], and other covariates were the same as in Equation 3. Because P[M.sub.10] was categorized into just two levels, each of [pm.sub.kt] and tem[p.sub.t-i] : [pm.sub.kt-i] has one coefficient denoted by [[beta].sub.2] and [[beta].sub.3], respectively. COVs was the same as in model 1.

S-plus software (version 6.2) was used in the data analyses (Chambers and Hastie 1993; Insightful Corporation 2001). To reduce potential bias caused by convergence, we used stricter criteria: 1.0 x [10.sup.-10] for both the local score algorithm and the backfitting algorithm (Dominici et al. 2002). We used the S-plus function gam.exact (Dominici et al. 2004; Internet-based Health and Air Pollution Surveillance System 2006) to correct the potential underestimation of the coefficient's standard error due to concurvity (Ramsay et al. 2003). Furthermore, analyses were restricted to days that contained values for all covariates in each model (> 87% of observations).

Results

We examined the distributions of each of the dependent variables, temperature, and P[M.sub.10] by time. The results show that there were strong seasonal patterns for RHA, CHA, REV, CEV, NECM, CM, and temperature, but the P[M.sub.10] pattern was less obvious (Figure 1). There were also apparent short-term fluctuations in health outcomes, minimum temperature, and the concentration of P[M.sub.10].

Table 1 provides summary statistics for individual health outcomes and explanatory variables. The results show considerable variation in each variable, ranging from 6 to 77 for RHA, 7 to 90 for CHA, 1 to 48 for REV, 2 to 38 for CEV, 5 to 42 for NECM, 1 to 31 for CM, 1.2 to 26 [degrees]C for temperature, and 2.5 to 60.0 [micro]g/[m.sup.3] for P[M.sub.10].

In the first model, there were inverse relationships between temperature and various measures of cardiorespiratory morbidity except for CHA, which showed a slight positive relationship (Figure 2). The patterns were similar at lags of 0, 1, or 2 days (RHA, CHA, REV, and CEV). Patterns of temperature effect on current day for morbidity and current day and lag 2 for mortality are presented (Figure 2). However, the relationships between temperature and cardiorespiratory mortality (NECM and CM) differed from the morbidity outcomes and varied by lengths of lag. For the current day, the associations of temperature with NECM and CM were relatively slight when the temperature was between 0 and 20[degrees]C and then increased quickly, but at lags of 1 and 2 days, the associations first decreased and then leveled off or slightly increased. Hence, the relationship between temperature and mortality forms a J- or U-shaped pattern, depending on the lag time (Basu and Samet 2002; Braga et al. 2002; Patz et al. 2000).

To explore potential effect modification of P[M.sub.10] and temperature on cardiorespiratory morbidity/mortality, we separately fitted bivariate response surface models (model 2) with individual health outcomes at each of three lags (0, 1, and 2 days). The results show interactive effects of P[M.sub.10] and temperature on RHA, REV, NECM, and CM at all time points, less so for CHA and CEV. Figure 3 illustrates the joint effects of P[M.sub.10] and temperature on each health outcome (RHA, CHA, REV, CEV, NECM, and CM) for the current day. Temperature effects were modified by levels of P[M.sub.10] for RHA, REV, NEMC, and CM, less so for CHA and CEV. Favorable temperature effects disappeared when P[M.sub.10] was above the mean or median (15.84 or 14.8 [micro]g/[m.sup.3]) for RHA and REV, but adverse temperature effects appeared for NECM and CM when P[M.sub.10] was above the mean or median. CHA was similar to RHA. The bivariate response surfaces differed from the independent model results, showing that the association between temperature and mortality changed with P[M.sub.10]. In fact, what at first appeared to be a J-shaped relationship in the independent model (model 1) became an approximately linear relationship when the joint effect of P[M.sub.10] and temperature was taken into account (model 2). There were inverse linear associations between minimum temperature and morbidity or mortality at low levels of P[M.sub.10] (< 20 [micro]g/[m.sup.3]). However, at higher levels of P[M.sub.10], the association between temperature and mortality was positively linear, whereas the associations with the various morbidity measures were weak.

Because no obvious J- or U-shaped patterns of the temperature-health relationship were observed in bivariate response models, we separately fitted nonstratification models (model 3) using each of the cardiorespiratory morbidity/mortality measures as a response variable with the same set of predictors at each of the lags (Table 2). The results indicate statistically significant interactive effects between temperature and P[M.sub.10] on RHA, CHA, NECM, and CM at different lags. For example, P[M.sub.10] modified the effects of temperature on RHA and CM at all lags, but modified the effects of temperature on NECM at lags of 0 and 2 days and CHA at lag 2, marginally at lag 0. No significant interactions were found for REV and CEV. The results were similar to those from model 2. Because the estimated effects of temperature variability differed with P[M.sub.10] levels, we present the estimated coefficients of model 3 instead of relative risks (Table 2).

To test sensitivity of changes in degrees of freedom related to the number of categories used for covariates, we refitted model 3 using 12 df for seaso[n.sub.t] (each month), 6 df for year (each year), and 8 df for rain, relative humidity, and ozone instead of the original df. Results show that increases in df changed the modeling outcomes only minimally.

Both the bivariate response surface and nonstratification models suggest that the effects of temperature on cardiorespiratory morbidity/mortality varied with levels of P[M.sub.10]. We then fitted the stratification model (model 4) to examine heterogeneity of temperature associations with health outcomes across different strata of P[M.sub.10], defined as greater than or less than the mean level (15.8 [micro]g/[m.sup.3]) of P[M.sub.10]. There were statistically significant interactions for RHA, CHA, NECM, and CM at different lags, but not for CEV and REV. Table 3 shows the percent changes in cardiorespiratory morbidity/mortality per 10[degrees]C increase in minimum temperature across the different levels of P[M.sub.10]. Temperature effects on cardiorespiratory morbidity/mortality varied across the different levels of P[M.sub.10]. For most lags and most health outcomes, the percent changes were higher when P[M.sub.10] levels were higher. For morbidity measures, this meant that the inverse association with temperature was less extreme for high P[M.sub.10]; for mortality measures, the association with minimum temperature became positive at high P[M.sub.10] levels. For example, when minimum temperature increased by 10[degrees]C (using data with the full range of temperature), there was a decrease in RHA of 7.2% and 1.0% on the current day, with low and high P[M.sub.10] levels, respectively. To examine the association at the high end of the temperature range with health outcomes, we also fitted model 4 using data sets constrained to the highest quartile ([greater than or equal to] 19.3[degrees]C) of temperature with the same cutoffs for temperature as the whole database. The pattern was even stronger when analyzed in the high-temperature data set (Table 3).

Discussion

In this study, we used three parallel time-series approaches to examine whether P[M.sub.10] modified the association between temperature and cardiorespiratory morbidity/mortality. Results show that P[M.sub.10] modified the effects of temperature on respiratory hospital admissions, cardiovascular hospital admissions, all non-external-cause mortality, and cardiovascular mortality in different lags. In particular, more adverse outcomes were evident with increasing temperature when P[M.sub.10] levels also increased. However, there were no significant interactive effects between temperature and P[M.sub.10] on respiratory and cardiovascular emergency visits. Three parallel models produced similar results.

In this study we used different health outcomes to examine consistency of findings. However, the findings from different health outcomes for the same observed groups varied somewhat. Reasons for this might include different age distributions, different events, and clinical features. For example, for respiratory hospital admissions and emergency visits, acute upper respiratory infection, pneumonia, and asthma were the dominant causes, and a high proportion of cases were identified in children. For cardiovascular hospital admission and emergency visits, angina pectoris, artrial fibrillation and flutter, and chronic ischemic heart diseases were the main causes, and elderly persons comprised most of theses cases. For MECN and CM, acute myocardial infarction, chronic ischemic heart diseases, and stroke were the main causes, and again elderly persons comprised most cases. Several studies have shown that age and preexisting diseases modify the air pollution-health association (Dubowsky et al. 2006; O'Neill et al. 2003). In considering the variation in findings for different health effects, we note that this study is designed to detect short-term effects of air pollution (within a few days). Mechanisms related to high temperature or P[M.sub.10] that precipitate acute illness may be different or have different magnitudes, depending on the underlying diseases.

Temperature and air pollution are generally highly correlated in many places (Holgate et al. 1999), and they may interact symmetrically to affect health outcomes. Although whether air pollution modifies temperature estimates has not been investigated so far, several studies have found evidence that temperature may modify the relationship between air pollution and morbidity or mortality (Choi et al. 1997; Katsouyanni et al. 1993; Ren and Tong, in press; Roberts 2004). For example, Katsouyanni et al. (1993) examined whether air pollution and ambient temperature had synergistic effects on excess mortality during the 1987 "heat wave" in Athens. They found a statistically significant modification of temperature on the association between exposure to sulphur dioxide and total excess mortality, although the main effect of this pollutant was not statistically significant. Roberts (2004) investigated the interaction between daily particulate air pollution and daily mean temperature on mortality in Cook County, Illinois, and Allegheny County, Pennsylvania, using data for 1987-1994. The study found that temperature modified the association between P[M.sub.10] and mortality, but the results were sensitive to the number of degrees of freedom. Our recent study also found that temperature significantly modified the association between P[M.sub.10] and health outcomes (Ren and Tong, in press). These findings support the hypothesis that P[M.sub.10] might modify the relationship between temperature and health outcomes.

It is biologically plausible that P[M.sub.10] modifies the effects of temperature on cardiorespiratory diseases. A range of studies have shown that P[M.sub.10] is consistently associated with health outcomes (Dominici et al. 2006; Samet et al. 2000). Exposure to P[M.sub.10] may directly affect airways through inhalation, including upper airways, bronchiole, and alveolus. The exposure could modulate the automatic nervous system and might further influence the cardiovascular system (Gordon 2003; Jeffery 1999). Some studies have shown that P[M.sub.10] is associated with decreased heart rate variation (Creason et al. 2001; Gold et al. 2000). Marked temperature changes also affect physiological and psychological stresses (Gordon 2003), which could aggravate preexisting diseases. Therefore, both high ambient temperature and high ambient P[M.sub.10] may interact to synergistically effect human morbidity/mortality.

Each of the three models used in this study has inherent advantages and disadvantages. The bivariate response surface model is a flexible approach to show the patterns of two continuous predictors on the dependent variable and explore whether potential interaction exists without a rigid assumption of linearing between predictors and the dependent variable (Greenland 1993; Hastie and Tibshirani 1990). However, this model can not provide parametric estimates for exposure effects; therefore, it may be difficult to judge whether interactive effects exist and also to compare the results from different studies. The nonstratification parametric model includes a pointwise product of two continuous variables. Parametric estimates for both predictors and their pointwise product can be obtained (Chambers and Hastie 1993). However, the linear assumption between the dependent variable and both continuous predictors is not necessarily met in all situations, especially for temperature and air pollution in a multisite study with variation in study populations. Moreover, the estimated coefficients of both predictors can not be simply interpreted as the main effects (Table 2) (Chambers and Hastie 1993). The stratification parametric model provides parametric estimates, which can be easily interpreted as main effects and interaction. The parametric estimates can be used in a meta-analysis. However, because the effect of one continuous predictor generally changed with another predictor level (Figure 3), the selection of cutoffs is still a challenge, especially in comparing the results from different studies.

We explored a marker of air pollution as a modifier of the relationship between temperature and cardiovascular morbidity/mortality. We used an independent model (model 1) and a bivariate response model (model 2) to examine the patterns of temperature with several health outcomes (RHA, CHA, REV, CEV, NECM, and CM). A J-shaped pattern was observed only for NECM and CM on the current day but not for other measures of cardiorespiratory morbidity. No obvious J-shaped pattern was observed in the bivariate response surface models, suggesting that the interaction between P[M.sub.10] and temperature may play an important role in model fit. Therefore, when modeling the health effects of air pollution and/or temperature, an interaction between these two factors should be carefully considered. Many studies have shown J-, U-, or V-shaped patterns of the temperature-health relationship (Basu and Samet 2002; Braga et al. 2002; Patz et al. 2000). The different patterns observed may be caused by different climate conditions across studies. For example, Braga et al. (2002) reported that greater variability of summer and winter temperature was associated with larger effects for hot and cold days, respectively, on respiratory deaths. However, Brisbane has a subtropical climate, with few extremely cold days (for example, the mean minimum temperature was 15.4[degrees]C and the lowest temperature was 1.2[degrees]C during the study period).

This study is an ecologic design, and misclassifications are possible for both health outcomes and exposure. Because a broad classification of diseases (cardiovascular, respiratory, and nonexternal classification of diseases) was used, we do not believe that misclassification for health outcomes is likely to be substantial. We used air pollution from one central monitoring site to represent individual exposure to P[M.sub.10] and this might result in misclassification. However, previous studies have shown that central fixed-site measurements may be treated as surrogates for personal exposure (Kim et al. 2005) and bias from the monitoring data might not be severe. Although some families in Brisbane have access to air conditioning, thus reducing exposure to high temperature when indoors, this effect is believed minimal due to the small proportion of houses with air conditioning as well as the outdoor lifestyle of Queensland residents.

There are two major strengths of this study. First, this is, to our knowledge, the first study to examine whether P[M.sub.10] modifies the association of temperature and a range of cardiorespiratory morbidity/mortality measures. Second, we performed three parallel statistical models with multiple health outcomes, and they produced similar findings, which strengthens the validity of findings.

However, this study also has two key limitations. First, caution is needed in interpreting any time-series study within a single location. This study was carried out in a single city with a subtropical climate, and 6 years of data are not extensive. Therefore, the results of this study may be difficult to generalize to other places. Second, this study is an ecologic design, in which bias from exposure measurement errors might occur to some degree due to lack of individual information.

Overall, we found statistically significant interactive effects of P[M.sub.10] and temperature on respiratory and cardiovascular hospital admission, all nonexternal-cause mortality, and cardiovascular mortality at different lags in Brisbane during the study period. The temperature effects were more adverse at high levels of P[M.sub.10]. These findings may have important implications in the assessment of health effects of temperature and the development of strategies and policies for controlling and preventing temperature-related deaths and diseases. However, it is necessary to determine whether a consistent finding could be found in other settings.

REFERENCES

Basu R, Samet JM. 2002. Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence. Epidemiol Rev 24:190-202.

Bell ML, McDermott A, Zeger SL, Samet JM, Dominici F. 2004. Ozone and short-term mortality in 95 US urbane communities, 1987-2000. JAMA 292:2372-2378.

Braga ALF, Zanobetti A, Schwartz J. 2001.The time course of weather related deaths. Epidemiology 12:662-667.

Braga ALF, Zanobetti A, Schwartz J. 2002. The effect of weather on respiratory and cardiovascular deaths in 12 U.S cities. Environ Health Perspect 110:859-863.

Brisbane City Council. 2006. Statistical Portrait of Brisbane. Available: http://svc189.bne146v.server-web.com/statistical_portrait/ [accessed 20 September 2006].

Chambers JM, Hastie TJ. 1993. Statistical Models in S. London:Chapman & Hall.

Choi K, Inou S, Shinozaki R. 1997. Air pollution, temperature, and regional differences in lung cancer mortality in Japan. Arch Environ Health 52:160-168.

Creason J, Neas L, Walsh D, Williams R, Sheldon L, Liao D, et al. 2001. Particulate matter and heart rate variability among elderly retirees: the Baltimore 1998 PM study. J Expo Anal Environ Epidemiol 11:116-122.

Daniels MJ, Dominici F, Samet JM, Zeger S. 2000. Estimating particulate matter-mortality dose-response curves and threshold levels: an analysis of daily time-series for the 20 largest US cities. Am J Epidemiol 152:397-406.

Dominici F, McDermott A, Hastie TJ. 2004. Improved semipara-metric time series models of air pollution and mortality. J Am Stat Assoc 99:938-948.

Dominici F, McDermott A, Zeger SL, Samet JM. 2002. On the use of generalized additive models in time-series studies of air pollution and health. Am J Epidemiol 156:193-203.

Dominici F, Peng RD, Bell ML, Pham L, McDermott A, Zeger SL, et al. 2006. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA 295:1127-1134.

Dubowsky SD, Suh Helen, Schwartz J, Coull BA, Gold DR. 2006. Diabetes, obesity, and hypertension may enhance associations between air pollution and markers of systemic inflammation. Environ Health Perspect 114:992-998.

Gold DR, Litonjua A, Schwartz J, Lovett E, Larson A, Nearing B, et al. 2000. Ambient pollution and heart rate variability. Circulation 101:1267-1273.

Gordon CJ. 2003. Role of environmental stress in the physiological response to chemical toxicants. Environ Res 92:1-7.

Greenland S. 1993. Basic problems in interaction assessment. Environ Health Perspect 101(suppl 4): 59-66.

Hastie T, Tibshirani R. 1990. Generalized Additive Models. New York:Chapman and Hall.

Holgate ST, Samet JM, Koren HS, Maynard RL. 1999. Air Pollution and Health. London:Academic Press.

Insightful Corporation. 2001. S-Plus 6 Guide to Statistics for windows. Seattle, WA:Insightful Corporation.

Internet-based Health & Air Pollution Surveillance System (iHAPSS). 2006. Software for Computing the Asymptotically Exact Standard Errors in GAM. Available:http://ihapss.biostat.jhsph.edu/software/gam.exact/gamDoc.htm [accessed 15 July 2006].

Jeffery P. 1999. Effects of cigarette and air pollutants on the lower respiratory tract. In: Air Pollution and Health (Holgate S, Samet J, Koren H, et al., eds). Sydney:Academic Press.

Katsouyanni K, Pantazopoulou A, Touloumi G, Tselepidaki I, Moustris K, Asimakopoulos D, et al. 1993. Evidence for interaction between air pollution and high temperature in the causation of excess mortality. Arch Environ Health 48:235-242.

Kim D, Sass-Kortsak A, Purdham JT, Dales RE, Brook JR. 2005. Associations between personal exposures and fixed-site ambient measurements of fine particulate matter, nitrogen dioxide, and carbon monoxide in Toronto, Canada. J Expo Anal Environ Epidemiol 15:1-12.

Kunst AE, Looman CW, Mackenbach JP. 1993. Outdoor air temperature and mortality in the Netherlands: a time-series analysis. Am J Epidemiol 137:331-341.

Martens WJM. 1997. Climate change, thermal stress and mortality changes. Soc Sci Med 46:331-344.

McMichael AJ, Githeko A, Akhtar R, Haines A, Kovats RS, Martens P, et al. 2001. Human health. In: Climate Change 2001; Impacts, Adaption, and Vulnerability (McCarthy J, Canziani O, Leary N, Dokken DJ, White KS. eds). New York:Cambridge University Press.

O'Neill MO, Zanobetti A, Schwartz J. 2003. Modifiers of temperature and mortality association in seven US cities. Am J Epidemiol 157:1074-1082.

Patz JA, Engelberg D, Last J. 2000. The effects of changing weather on public health. Annu Rev Public Health 21:271-303.

Queensland Government. 2001. Public Health Law and Indigenous Health Project--The Impact of the Health Act 1937 on Public Health in Queensland's Remote Indigenous Community. Available: http://www.health.qld.gov.au/phs/documents/osm/16395.pdf [accessed 15 July 2006].

Rainham DGC, Smoyer-Tomic KE. 2003. The role of air pollution in the relationship between a heat stress index and human mortality in Toronto. Environ Res 93:9-19.

Ramsay TO, Burnet RT, Krewski D. 2003. The effect of concurvity in generalized additive models linking mortality to ambient particulate matter. Epidemiology 14:18-23.

Ren C, Tong S. 2006. Temperature modifies the health effects of particulate matter in Brisbane, Australia. Int J Biometeorol doi: 10.1007/s00484-006-0054-7 [Online 12 September 2006].

Roberts S. 2004. Interaction between particulate air pollution and temperature in air pollution mortality time series studies. Environ Res 96:328-337.

Samet JM, Dominici F, Curriero FC, Coursac I, Zeger SLl. 2000. Fine particulate air pollution and mortality in 20 U.S. cities, 1987-1994. N Engl J Med 343:1742-1749.

Samet JM, Zeger S, Kelsall JE, Xu J, Kalkstein L. 1998. Does weather confound or modify the association of particulate air pollution with mortality? An analysis of the Philadelphia data, 1973-1980. Environ Res (Section A) 77:9-19.

Schwartz J, Samet JM, Patz JA. 2004. Hospital admissions for heart disease--the effects of temperature and humidity. Epidemiology 15:755-761.

World Health Organization. 1975. International Classification of Diseases, 9th Revision. Geneva:World Health Organization.

World Health Organization. 1992. International Classification of Diseases, 10th Revision. Geneva:World Health Organization.

Cizao Ren, (1) Gail M, Williams, (2) and Shilu Tong (1)

(1) School of Public Health, Queensland University of Technology, Kelvin Grove, Brisbane, Queensland, Australia; (2) School of Population Health, University of Queensland, Herston, Brisbane, Queensland, Australia

Address correspondence to C. Ren, School of Public Health, Queensland University of Technology, Kelvin Grove, Qld. 4059 Australia. Telephone: [phone]. Fax: [phone]. E-mail: [e-mail]

We thank J.M. Samet (Johns Hopkins University) and B. Newman (Queensland University of Technology) for their insightful comments, and the Queensland Environmental Protection Agency, Queensland Health Department, Australia Bureau of Meteorology and Australia Bureau of Statistics for providing data.

The study is funded partly by the Australian Research Council.

The authors declare they have no competing financial interests.

Received 14 April 2006; accepted 27 July 2006.

Table 1. Summary statistics for health outcomes, air pollutants, and
meteorologic conditions.

                                 25th                75th
Variable        Mean    Minimum  Percentile  Median  Percentile  Maximum

RHA (n)         33.16    6       23          33      42           77
CHA (n)         39.32    7       30          39      48           90
REV (n)         18.39    1       13          17      23           48
CEV (n)         17.17    2       14          17      20           38
NECM (n)        15.86    5       13          15      18           42
RM (n)           1.45    0        0           1       2            8
CM (n)           7.26    1        5           6       9           31
P[M.sub.10]     15.84    2.5     11.9        14.8    18.7         60
  ([micro]g/
  [m.sup.3])
Temperature     15.42    1.2     12.0        16.0    19.3         26
  ([degrees]C)
Humidity (%)    65.21   12.0     57          65      74          100
Ozone (ppb)     11.26    0        7          10      15           45
Rainfall (mm)    3.038   0        0           0.0     0.8        163

Table 2. Coefficients for main and interactive effects of minimum
temperature and P[M.sub.10] on different health outcomes using model 3.

Lag (a)  Variable     RHA          CHA            REV

Lag0     Temperature  -0.012019*   -0.003942      -0.009377*
         P[M.sub.10]  -0.004296**   0.000150      -0.000887
         Interaction   0.000471*    0.000232 (#)   0.000135
Lag1     Temperature  -0.009196*   -0.0003767     -0.010039**
         P[M.sub.10]  -0.002474**   0.000028      -0.004209
         Interaction   0.000339**   0.000124       0.000271
Lag2     Temperature  -0.008950*   -0.009764*     -0.009406*
         P[M.sub.10]  -0.004229**  -0.002946      -0.003440
         Interaction   0.000413*    0.000259**     0.000280 (#)

Lag (a)  Variable     CEV          NECM           CM

Lag0     Temperature  -0.003304    -0.008374**    -0.008193
         P[M.sub.10]   0.000737    -0.004022      -0.004565
         Interaction  -0.000015     0.000534*      0.000603**
Lag1     Temperature  -0.003636    -0.005321      -0.010081
         P[M.sub.10]  -0.001248    -0.002182      -0.003927
         Interaction   0.000116     0.000322 (#)   0.000556**
Lag2     Temperature  -0.008885**  -0.008556 (#)  -0.015128**
         P[M.sub.10]  -0.003383    -0.005177      -0.007141
         Interaction   0.000210     0.000441**     0.000809**

(a) Lag refers to 0, 1, or 2 days. *p < 0.01. **p < 0.05. (#) p < 0.10.

Table 3. Percent change (%) in cardiorespiratory morbidity/mortality per
10[degrees]C increase in temperature across the levels of P[M.sub.10].

Lag (a)                          Variable    RHA (95% CI)

Whole range of temp
Lag0                             PM Low (b)   -7.2 (-11.3 to -2.9)
                                 PM High      -1.0 (-5.0 to 3.2)
Lag1                             PM Low       -2.9 (-7.3 to 1.6)
                                 PM High      -2.4 (-6.5 to 1.8)
Lag2                             PM Low       -3.2 (-7.6 to 1.3)
                                 PM High      -0.4 (-4.5 to 3.8)
Temp [greater than or equal to]
  19.3[degrees]C
Lag0                             PM Low      -29.2 (-40.6 to -15.5)
                                 PM High       5.2 (-17.0 to 33.4)
Lag1                             PM Low       -7.6 (-23.4 to 11.5)
                                 PM High      -8.1 (-27.6 to 16.6)
Lag2                             PM Low      -12.8 (-27.4 to 4.8)
                                 PM High     -19.2 (-36.4 to 2.6)

Lag (a)                          Variable    CHA (95% CI)

Whole range of temp
Lag0                             PM Low (b)  -2.3 (-6.3 to 1.7)
                                 PM High      1.1 (-2.5 to 4.7)
Lag1                             PM Low      -2.6 (-6.7 to 1.4)
                                 PM High     -1.0 (-4.6 to 2.5)
Lag2                             PM Low      -8.2 (-12.2 to -4.2)
                                 PM High     -3.5 (-7.1 to 0.1)
Temp [greater than or equal to]
  19.3[degrees]C
Lag0                             PM Low      -4.7 (-19.2 to 10.0)
                                 PM High     -3.8 (-22.8 to 15.6)
Lag1                             PM Low       0.7 (-13.7 to 15.3)
                                 PM High     12.2 (-7.4 to 32.1)
Lag2                             PM Low      15.9 (1.1 to 30.8)
                                 PM High     12.8 (-6.7 to 32.7)

Lag (a)                          Variable    REV (95% CI)

Whole range of temp
Lag0                             PM Low (b)   -6.8 (-12.1 to -1.1)
                                 PM High      -6.6 (-11.5 to -1.4)
Lag1                             PM Low       -3.6 (-9.2 to 2.3)
                                 PM High      -5.3 (-10.3 to -0.1)
Lag2                             PM Low       -3.9 (-9.5 to 2.0)
                                 PM High      -4.4 (-9.4 to 1.9)
Temp [greater than or equal to]
  19.3[degrees]C
Lag0                             PM Low       -7.4 (-27.2 to 17.8)
                                 PM High       4.5 (-23.1 to 42.0)
Lag1                             PM Low        1.2 (-20.2 to 28.3)
                                 PM High      51.5 (11.1 to 106.6)
Lag2                             PM Low      -24.8 (-40.8 to -4.6)
                                 PM High      -8.8 (-33.3 to 24.7)

Lag (a)                          Variable    CEV (95% CI)

Whole range of temp
Lag0                             PM Low (b)   -2.2 (-7.8 to 3.5)
                                 PM High      -4.5 (-9.7 to 0.6)
Lag1                             PM Low       -2.1 (-7.8 to 3.7)
                                 PM High      -1.2 (-6.3 to 3.9)
Lag2                             PM Low       -6.7 (-12.4 to -1.0)
                                 PM High      -4.2 (-9.3 to 1.0)
Temp [greater than or equal to]
  19.3[degrees]C
Lag0                             PM Low      -13.2 (-33.8 to 7.8)
                                 PM High     -94.2 (-118.6 to 69.1)
Lag1                             PM Low      -14.9 (-35.3 to 6.0)
                                 PM High       1.5 (-25.7 to 29.6)
Lag2                             PM Low        8.0 (-13.0 to 29.4)
                                 PM High      29.4 (-51.8 to 117.5)

Lag (a)                          Variable    NECM (95%CI)

Whole range of temp
Lag0                             PM Low (b)  -1.4 (-7.3 to 4.8)
                                 PM High      2.8 (-2.7 to 8.7)
Lag1                             PM Low      -0.2 (-6.2 to 6.0)
                                 PM High      0.6 (-4.8 to 6.4)
Lag2                             PM Low      -3.9 (-9.6 to 2.2)
                                 PM High      1.1 (-4.3 to 6.9)
Temp [greater than or equal to]
  19.3[degrees]C
Lag0                             PM Low       9.9 (-12.9 to 38.5)
                                 PM High     14.0 (-15.4 to 53.7)
Lag1                             PM Low      35.1 (7.0 to 70.5)
                                 PM High     14.8 (-15.4 to 55.8)
Lag2                             PM Low      13.2 (-10.2 to 42.6)
                                 PM High     37.6 (1.9 to 85.8)

Lag (a)                          Variable    CM (95% CI)

Whole range of temp
Lag0                             PM Low (b)  -0.9 (-9.8 to 7.9)
                                 PM High      4.6 (-3.4 to 12.6)
Lag1                             PM Low      -1.4 (-10.2 to 7.5)
                                 PM High      0.0 (-7.9 to 8.0)
Lag2                             PM Low      -3.6 (-12.5 to 5.2)
                                 PM High      0.9 (-7.0 to 8.8)
Temp [greater than or equal to]
  19.3[degrees]C
Lag0                             PM Low       0.8 (-31.3 to 34.1)
                                 PM High     44.1 (0.8 to 89.2)
Lag1                             PM Low      18.9 (-13.6 to 52.4)
                                 PM High     26.2 (-17.3 to 71.6)
Lag2                             PM Low       6.7 (-25.4 to 39.8)
                                 PM High     23.2 (-19.1 to 67.4)

Abbreviations: CI, confidence interval; Temp, minimum temperature.
(a) Lag refers to 0, 1, or 2 days. (b) PM indicates P[M.sub.10],
categorized into low and high levels using mean of P[M.sub.10] as a
cutoff.