Skip to content

Environmental Health Perspectives

Facebook Page EHP Twitter Feed Open Access icon  

Research December 2014 | Volume 122 | Issue 12

Email this to someoneShare on FacebookTweet about this on TwitterShare on LinkedInShare on Google+Share on StumbleUpon
Environ Health Perspect; DOI:10.1289/ehp.1306954

Variability in Temperature-Related Mortality Projections under Climate Change

Tarik Benmarhnia,1,2,3 Marie-France Sottile,4,5 Céline Plante,6 Allan Brand,3,7 Barbara Casati,4 Michel Fournier,6 and Audrey Smargiassi1,3,7

Author Affiliations open
1Département de santé environnementale et santé au travail (DSEST), Université de Montréal, Montréal, Québec, Canada; 2EHESP School of Public Health (École des hautes études en santé publique), Rennes, Sorbonne-Paris Cité, Rennes, France; 3Chaire sur la pollution de l’air, les changements climatiques et la santé, DSEST, Université de Montréal, Montréal, Québec, Canada; 4Consortium Ouranos, Montréal, Québec, Canada; 5Ministère du Développement durable, de l’Environnement et des Parcs, Québec, Canada; 6Direction de santé publique de l’Agence de la santé et des services sociaux de Montréal, Québec, Canada; 7Institut National de Santé Publique du Québec, Montréal, Québec, Canada

PDF icon PDF Version (254 KB)

  • Background: Most studies that have assessed impacts on mortality of future temperature increases have relied on a small number of simulations and have not addressed the variability and sources of uncertainty in their mortality projections.

    Objectives: We assessed the variability of temperature projections and dependent future mortality distributions, using a large panel of temperature simulations based on different climate models and emission scenarios.

    Methods: We used historical data from 1990 through 2007 for Montreal, Quebec, Canada, and Poisson regression models to estimate relative risks (RR) for daily nonaccidental mortality in association with three different daily temperature metrics (mean, minimum, and maximum temperature) during June through August. To estimate future numbers of deaths attributable to ambient temperatures and the uncertainty of the estimates, we used 32 different simulations of daily temperatures for June–August 2020–2037 derived from three global climate models (GCMs) and a Canadian regional climate model with three sets of RRs (one based on the observed historical data, and two on bootstrap samples that generated the 95% CI of the attributable number (AN) of deaths). We then used analysis of covariance to evaluate the influence of the simulation, the projected year, and the sets of RRs used to derive the attributable numbers of deaths.

    Results: We found that < 1% of the variability in the distributions of simulated temperature for June–August of 2020–2037 was explained by differences among the simulations. Estimated ANs for 2020–2037 ranged from 34 to 174 per summer (i.e., June–August). Most of the variability in mortality projections (38%) was related to the temperature–mortality RR used to estimate the ANs.

    Conclusions: The choice of the RR estimate for the association between temperature and mortality may be important to reduce uncertainty in mortality projections.

  • Citation: Benmarhnia T, Sottile MF, Plante C, Brand A, Casati B, Fournier M, Smargiassi A. 2014. Variability in temperature-related mortality projections under climate change. Environ Health Perspect 122:1293–1298;

    Address correspondence to A. Smargiassi, Direction de Santé Publique de Montréal, 1301 Sherbrooke East, Montreal, Québec, H2L 1M3 Canada. Telephone: (514) 528-2400. E-mail:

    We are grateful to T. Kosatsky for his useful comments. We also acknowledge the Data Access Integration (DA) Team for providing the data and technical support.

    This project was supported by the Québec Government Fonds vert of the Action 21 of the Plan d’action 2006–2012 sur les changements climatiques (PACC) and by Health Canada. The DAI Portal (​index.php/pages/dai) is made possible by a collaboration among the Global Environmental and Climate Change Centre (GEC3), the Adaptation and Impacts Research Division (AIRD) of Environment Canada, and the Drought Research Initiative (DRI). The Ouranos Consortium (in Quebec) provides access of the Canadian Regional Climate Model (CRCM).

    The authors declare they have no actual or potential competing financial interests.

    Received: 15 April 2013
    Accepted: 16 July 2014
    Advance Publication: 18 July 2014
    Final Publication: 1 December 2014

  • PDF icon Supplemental Material (1.2 MB) PDF


Among environmental determinants of health, weather and climate have received increasing attention related to awareness of climate change and the documentation of both usual and catastrophic heat-related mortality (Campbell-Lendrum and Woodruff 2007; Ebi 2008; Kovats and Hajat 2008; O’Neill and Ebi 2009; Patz et al. 2005).

Increases in ambient summer temperatures over city-specific thresholds have been associated with an increase in mortality [Hajat and Kosatsky 2010; Intergovernmental Panel on Climate Change (IPCC) 2007a, 2007b]. Although present-day health effects of summer temperatures have been well characterized (Armstrong 2006), the extent to which future changes in summer temperatures will affect human health has received relatively little attention (Campbell-Lendrum and Woodruff 2007; Costello et al. 2009; Ebi and Gamble 2005; Huang et al. 2011; McMichael et al. 2006; Peng et al. 2011). Global average ambient temperatures have been projected to increase under any scenario of increasing greenhouse gas (GHG) concentrations (IPCC 2001, 2007a). Various models have been developed in climate science (Caya et al. 1995; McFarlane et al. 1992) to estimate future temperatures according to different GHG emission scenarios [IPCC Special Report on Emissions Scenarios (SRES); IPCC 2000].

Studies published since 2008 that have estimated the impacts of future temperatures on mortality have mainly been conducted in Europe and North America, and the time periods used for baseline data and projections have varied among them. Most studies have evaluated projected temperatures based on only a small number of climate change simulations, with the exception of Li et al. (2013). For a review of most of these studies, see Huang et al. (2011) and Gosling et al. (2009).

Divergences in temperature projections, for example, due to model structure and GHG emissions scenarios, may occur; to capture the maximal range of possible future temperatures and health impacts, it may thus be important to consider a large number of simulations when assessing health impacts.

We assessed the variability of temperature projections and future mortality distributions, using a large panel of temperature simulations based on climate models and emission scenarios for the period 2020–2037 on the island of Montreal, which includes the city of Montreal, the largest city of the province of Quebec, Canada.


To predict mortality attributable to past (1990–2007) and future temperatures (2020–2037) in Montreal, we first used historical data for 1990–2007 to estimate three sets of relative risks (RRs) for associations between mortality and temperatures during June–August on the island of Montreal with Poisson models. We refer to “sets” of RRs, given that there is one RR per degree temperature due to the nonlinear relation between temperature and mortality. One set of RRs was based on the observed historical data, and two were based on individual bootstrap samples constructed from the observed historical data that generated the 95% confidence interval (CI) of the attributable number (AN) of deaths. We then used the three sets of RRs to predict mortality attributable to past (1990–2007) and future temperatures (2020–2037), using the observed historical temperatures and 32 different temperature simulations (for 1990–2007 and 2020–2037) based on three general circulation models (GCMs) and the Canadian Regional Climate Model (RCM) (Caya et al. 1995). Predicted temperature distributions were corrected by applying the daily translation method (Mpelasoka and Chiew 2009) based on observed versus simulated temperatures for 1990–2007. Finally, using analysis of covariance (ANCOVA), we examined factors that influence the variability in future ANs, including the three sets of RRs and the 32 temperature simulations used to estimate future ANs (2020–2037). The project was carried out in the context of the Quebec ministerial health surveillance plan, which obtained ethics approval from the Quebec Public Health Ethical Health Surveillance Committee.

Observations. Mortality data. The mortality file comprises data on residents of Montreal who died on the island during June, July, and August of 1990 through 2007. We included all underlying nonaccidental causes of death and excluded deaths for International Classification of Diseases, Revision 9 (ICD-9) codes 800–999 (injury and poisoning) and Revision 10 (ICD-10) codes S00-T98 (injury, poisoning, and other consequences of external causes).

Observed temperatures and ozone (O3) levels. We computed daily (0000–2300 hours) mean, minimum, and maximum outdoor temperatures using data for June–August 1990–2007 obtained from the Environment Canada meteorological observation station at the Montréal Pierre Elliott Trudeau International Airport (Environment Canada 2012), located about 20 km from the city core (Doyon et al. 2008). We obtained hourly measurements of O3 at seven fixed-site monitoring stations from the Environment Canada National Air Pollution Surveillance Network (Environment Canada 2012). We averaged hourly concentrations over all stations and computed daily mean concentrations of O3 from these values for the summers (June–August) of 1990–2007.

Projection of temperatures and methods applied to correct simulated temperatures. Models and simulations. The 32 simulations used in the present study provided climate data for 1990–2007 and 2020–2037, as well as daily temperature estimates (Table 1). We selected the near future period 2020–2037 for projections to provide a climate change signal out of the climate natural variability that still corresponds to a near future for public health consideration. For each climate model, we considered the temperature time series for the grid point nearest to Montreal. Twenty-two of the simulations used in this study were based on the GCM simulations of the World Climate Research Programme (WCRP) Coupled Model Intercomparison Project, phase 3 (CMIP3) multimodel data set (Meehl et al. 2007), which used three SRES GHG emissions scenarios: B1, A1B, and A2 representing mild, medium, and strong future emissions of greenhouse gases and aerosols, respectively (IPCC 2000). GCMs have a resolution of 200–300 km, and the RCM has a resolution of approximately 50 km. The ensemble of the CMIP3 climate models is the multimodel data set used for the IPCC Fourth Assessment Report (IPCC 2007a). The CMIP3 simulations are available from the Program for Climate Model Diagnostic and Intercomparison archive ( Most of the CMIP3 GCMs covered three discrete time periods (1961–2000, 2046–2065, and 2081–2100), whereas the time periods selected for the present study (1990–2007 and 2020–2037) limited our choices to three GCMs with continuous simulations. Therefore, we used 15 simulations (with varying initial conditions) of the Canadian Global Circulation Model, version 3.1 (CGCM3.1/T47; Flato et al. 2000), developed at the Canadian Centre for Climate Modeling & Analysis (CCCMA) and three simulations of the same model truncated at a finer resolution (CGCM3.1/T63); three simulations of the Mark 3.5 (Mk3.5) climate model developed at the CSIRO Atmospheric Research, Australia (Cai et al. 2005); and one simulation of the German Coupled Global Climate Model (ECHAM5; Jungclaus et al. 2006) developed at the Max Planck Institute (MPI) of Meteorology (Germany). See Supplemental Material, Figure S1, for grids representing the GCM models.

Table 1 - See HTML for full tableTable 1 – Simulations of temperatures and climate models used.

View Table (HTML Version)
View larger image (TIF File)

In addition to the GCM simulations, we used 10 simulations of the Canadian Regional Climate Model (CRCM; Caya et al. 1995; Laprise 2008; Plummer et al. 2006), versions 4.1 and 4.2 (de Elía and Côté 2010; Paquin 2010) that mirror the recent CRCM evolution (within version 4) and include some minor modifications in parameters associated with surface processes and O3 data. RCMs use GCMs to produce climate projections at a higher resolution on a regional (usually continental) scale (see Christensen et al. 2007 and references therein). CRCM simulations over two domains were considered: a domain covering North America (200 × 192 grid points) and a domain centered over Quebec (111 × 87 grid points), both with a horizontal grid-size mesh of 45 km (true at 60°N). Grids representing the RCM models are shown in Supplemental Material, Figure S2.

There are different sources of uncertainty that affect temperature projections (de Elía and Côté 2010; Déqué et al. 2007). These include a) the effects introduced by the RCM downscaling on temperature projections (here termed resolution); b) the GHG emission scenarios (SRES) based on different states of future human activities; c) the climate model (GCM) itself (each model is conceived with a different structural design, numerical scheme, and physical parameterizations); d) the initial conditions, which—because of the chaotic nature of the climate system—generate an intrinsic “natural variability” in the climate model response; and e) the domain (the area covered by the RCM simulation). In the present study, we assessed the influence of different simulations on the variability of temperature projections and thus on the variability of future mortality distributions.

To ensure that the 32 simulations covered most of the range of possible future temperature (minimum, maximum, and mean) simulations, we computed the differences between the average of future (2046–2065) and historical simulated temperature distributions (1971–2000) for the three temperature metrics (daily mean, minimum, and maximum temperature) based on the 32 simulations versus all simulations provided by the Ouranos Consortium on Regional Climatology and Adaptation to Climate Change ( (n = 127 simulations for mean temperatures, n = 111 for maximum and minimum temperatures). This analysis confirmed that the 32 simulations included in the present analysis covered most of the temperature variability distribution from all simulations available (see Supplemental Material, Figure S3).

Correction for simulated temperatures. GCMs and RCMs project future temperatures at specific locations with errors related to the scale of the model predictions and other factors (Mpelasoka and Chiew 2009). Correction methods originally developed for hydroclimatology studies (Teng et al. 2011) include notably constant scaling, daily scaling, and daily translation (DT) methods (Mpelasoka and Chiew 2009). For the present analysis, we used the DT method to derive correction factors based on differences between observed historical data for 1990–2007 from a local meteorological station and simulated data for the same location and time period, thereby accounting for errors related to scaling as well as other sources of error.

Deaths attributable to historical and future temperatures. Relationships between observed historical temperatures and mortality. We used separate generalized linear Poisson models to estimate associations between daily death counts [relative risks (RRs) compared with the daily mean death count over the entire period] and daily mean, maximum, and minimum observed temperatures during June, July, and August of 1990–2007 (Armstrong 2006). We used cubic B-splines of time with 5 degrees of freedom (splines package in R; to control for secular trends in the mortality series (Hajat et al. 2007) and modeled the day of the season (1–92) using a spline with three knots (10th, 50th, and 90th percentiles) to control for seasonal patterns. To account for the nonlinear relationship between mortality and temperature, we modeled each temperature variable as a cubic spline with 5 knots (corresponding to 0, 25th, 50th, 75th, and 100th percentiles). Thus, a different RR was estimated for each degree of temperature, creating a set of RRs based on observed historical data. We also included daily mean levels of O3 in the regression as a simple continuous variable. To ensure that no autocorrelation remained in the residuals, we visually inspected partial autocorrelation plots and used the white noise statistical test (null hypothesis not rejected at p > 0.05; Lobato and Velasco 2004). We also visually inspected the plots of modeled Pearson residuals against the predicted values to verify that there was no important overdispersion.

Generation of 95% CIs of the ANs. There is no standard analytical way to estimate the 95% CI bounds of the distribution of a set of RRs (i.e., when there is one RR per degree temperature due to the nonlinear relation). We thus estimated the statistical uncertainty of the ANs calculated using the observed historical temperatures with 1,000 bootstrap samples, from which we selected the 2.5% and the 97.5% of the ANs based on the observed data, and the corresponding sets of RRs that produced them (i.e., corresponding point estimates per degree temperature produced with the same parameters as the Poisson model developed with the historical data).

The bootstrap samples were developed from the observed daily data as follows. For each of the 1,000 samples, we drew with replacement 18 times from day ones, the day twos, and so on (bootstrap samples were stratified by day of “summer”). Thus each of the 1,000 samples contained 1,656 days (92 days × 18 years).

Calculation of deaths attributable to temperatures. We first calculated the attributable fraction of daily deaths (AF) for each daily temperature metric value Ti (mean, maximum, or minimum) using the three sets of RRs described above (i.e., one based on the observed historical data, and two based on individual bootstrap samples). We used Equation 1 to calculate the AF only when the RR was > 1 and the daily minimum, mean, and maximum temperatures were above 15°C, 20°C, and 20°C, respectively.

AF(Ti) = [RR(Ti) – 1]/RR(Ti). [1]

The total number of ANs per year for the 1990–2007 and 2020–2037 periods for a given temperature metric (mean, maximum, or minimum) was then estimated using Equation 2:

AN = Σ[AF(Ti) × MDC × ND(Ti)], [2]

where MDC is the mean observed daily death count for the period 1990–2007 and ND(Ti) is the number of days with the temperature metric (observed or simulated) of Ti; values are summed from the minimum value of Ti (i.e., 15°C, 20°C, and 20°C for daily minimum, mean, and maximum temperatures, respectively) for which the estimated RR for the temperature metric and mortality based on historical data was > 1 up to the maximum value of Ti.

In our calculation of deaths attributable to future temperatures, we assumed that there would be no change in the mean daily death count in the future, no demographical change, no change in O3 levels, and no adaptation to heat from populations.

Variability analysis. We studied the influence of the different simulations on the future temperature and mortality distributions (2020–2037) focusing on daily mean temperature (not on daily maximum or minimum).

Temperature projection variability. We performed an ANCOVA in which the variable to be explained was the simulated daily mean temperature and the predictors were the simulation (n = 32) and the year modeled (n = 18).

Variability in ANs. We also used an ANCOVA to evaluate the influence of the set of RRs used to represent the temperature–mortality association, the simulation used to project future temperatures, and the year of the simulation on the estimated number of deaths attributable to temperature during each future year of 2020–2037. The three sets of RRs used included the set based on the observed historical data and the two sets of RRs based on the bootstrap data samples that generated the 95% CIs of ANs.


Observations. For the 18-year period 1990–2007 in Montreal, 61,356 nonaccidental deaths occurred during June–August, of which 79.9% were among people > 65 years of age. The average observed values for daily mean, maximum, and minimum temperatures were 20.4°C, 24.9°C, and 15.6°C for June, July, and August of 1990–2007. For the distributions of observed daily temperatures and O3 levels for the summers (i.e., June–August) of 1990–2007, see Supplemental Material, Table S1. The average daily O3 concentration for the same period was 25.3 μg/m3 (range, 0.83–76.1 μg/m3) (see Supplemental Material, Table S1).

Simulated temperatures. Historical temperatures. Average values of simulated daily mean temperatures for 1990–2007 over the 32 simulations were lower than observed values before the DT correction was applied (16.6°C compared with 20.4°C), whereas simulated mean daily temperatures were closer to the observed values after correction (e.g., mean 20.1°C) (Table 2). We found similar results with daily maximum and minimum temperatures (data not shown).

Table 2 - See HTML for full tableTable 2 – Distribution of observed (1990–2007), historical simulated (1990–2007), and future simulated (2020–2037) average daily mean temperatures (°C) (= 1,656 days).

View Table (HTML Version)
View larger image (TIF File)

Future temperatures. All DT-corrected simulations for June–August 2020–2037 suggested an increase in daily mean temperatures in Montreal compared with observed temperatures (Table 2): The average daily mean temperature from the 32 future simulations was 20.9°C (range: 20.3–21.3°C), compared with 20.4°C for the observed temperatures for June–August 1990–2007. We observed notable differences between the 32 DT-corrected and uncorrected values (uncorrected average daily mean, 17.6°C). However, the effect of the correction was not constant at all percentiles; for the future simulations, we observed greater differences in the ranges at higher and lower percentiles. For example, the average (range) of the corrected and uncorrected future simulated daily mean temperatures for the lowest percentile were 8.2°C (5.2–10.5°C) and 5.5°C (3.8–8.7°C), respectively, and for the 50th percentile were 21.1°C (20.3–21.6°C) and 17.7°C (17.1–18.1°C), respectively. Results were similar for daily maximum and minimum temperatures (data not shown).

Variability in temperature projections. Table 3 presents the percent of variance in daily temperature projections explained by the simulation and the year. Less than one percent of the daily temperature variation was explained by the choice of simulation model and year simulated.

Table 3 - See HTML for full tableTable 3 – Effect of simulation and year on daily mean temperature projections from an ANCOVA model (n = 52,992).

View Table (HTML Version)
View larger image (TIF File)

Deaths attributable to historical and future temperatures. For the three sets of RRs for each temperature unit used to calculate deaths attributable to temperature, see Supplemental Material, Table S2. Applying the RR (> 1) of the temperature–mortality relationships, we estimated a mean of 62 (95% CI: 32, 86) attributable deaths per summer (i.e., June–August) for the years 1990–2007. For maximum and minimum historical daily temperatures, we estimated 55 (95% CI: 32, 79) and 38 (95% CI: 9, 61) attributable deaths, respectively. Deaths attributable to temperatures and 95% CIs calculated by bootstrapping are presented in Supplemental Material, Table S3.

The estimated numbers of deaths attributable to daily mean, maximum, and minimum temperatures during June–August 2020–2037 based on the 32 DT-corrected simulations are presented in Figure 1, along with estimated numbers of deaths attributable to temperatures during 1990–2007 based on observed and DT-corrected simulated temperatures. Average numbers of deaths attributable to daily mean temperature during each year for 2020–2037 based on all 32 simulations were higher than the average number based on observed mean daily temperatures for 1990–2007. For maximum and minimum daily temperatures, 100% and 72% of the projections, respectively, produced estimated ANs of daily deaths above the average estimated for observed daily temperatures in 1990–2007. Average ANs of deaths based on simulated daily mean temperature during each year for 1990–2007 were similar to the average numbers based on observed mean daily temperatures for that time period. However, we observed differences in ANs between simulated minimum and maximum temperatures for this period.

Figure 1 - Annual average deaths attributable to temperature during June–August (y-axis) estimated based on observed data for 1990–2007, simulated data for 1990–2007, and simulated data for 2020–2037 and daily mean, maximum, or minimum daily temperature values (x-axis).Figure 1 – Estimated average annual deaths attributable to temperature (daily mean, daily maximum, or daily minimum) during June–August based on observed data for 1990–2007 and simulated data for 1990–2007 and 2020–2037. Simulated data are based on 32 simulations from RCMs and GCMs and corrected by the DT method. Deaths attributable to simulated temperatures were estimated with a set of RRs based on the observed historical data. The upper lines represent the ANs calculated with the set of RRs that generated the upper 95% CI bound of the ANs for the highest attributable number of deaths; the lower lines represent the ANs calculated with the set of RRs that generated the lower 95% CI bound of the ANs for the lowest attributable number of deaths. The three sets of RRs were also used to generate the attributable numbers for the observed data.

View larger image (TIF File)

Variability in mortality projections. Using the RR based on the observed historical data and the 32 simulations to estimate ANs for 2020–2037, we noted a high degree of variability in summer deaths attributable to daily temperatures (range, 65–129 summer deaths for mean daily temperatures; 78–161 for maximum daily temperatures; and 30–53 for minimum daily temperatures) (see Supplemental Material, Table S3).

Table 4 presents the percentage of variance in future yearly deaths attributable to temperature projections that is explained by the simulation, the year, and the set of RRs used to calculate the ANs. A higher part of the variability in mortality projections than in temperature projections was explained by the simulation: Six percent of the yearly variation of deaths attributable to future temperatures was explained by the choice of the simulation, compared with < 1% for the projections of temperatures. Nonetheless, most of the variability in mortality projections (38%) was related to the temperature–mortality RRs used to estimate the attributable fraction of heat-related deaths.

Table 4 - See HTML for full tableTable 4 – Effect of simulation, year, and set of RRs on deaths attributable to future temperatures by season from an ANCOVA model (n = 1,728).

View Table (HTML Version)
View larger image (TIF File)

Figure 2 shows the yearly average estimated number of attributable deaths associated with each daily mean temperature value with RR > 1; the numbers for each distribution sum to the average total estimated number of attributable deaths. Most of the estimated attributable deaths occur on days with daily mean temperatures between 24°C and 28°C; this reflects the number of days with temperatures in this range.

Figure 2 - Graph of total number deaths attributable to temperature during June–August (y-axis) according to mean daily temperature (x-axis) and average values for 1990–2007, average values for 1990–2007, estimated average values for 2020–2037, and estimated minimum and maximum values for 2020–2037.Figure 2 – Yearly estimates of average, minimum, and maximum numbers of deaths during June–August attributable to observed temperatures in 1990–2007 and predicted temperatures in 2020–2037 based on 32 simulations (corrected using the DT method) according to mean daily temperature (°C).

View larger image (TIF File)


In the present study, we estimated the variability in projections of deaths attributable to temperatures in Montreal during June–August 2020–2037 using 32 RCM and GCM temperature simulations (with different climate models, SRES, domains, versions, and members). Using DT-corrected simulated June–August temperatures, we found an increase in estimated numbers of deaths attributable to daily mean ambient temperatures during 2020–2037, with a large variability ranging from 34 to 174 deaths per summer, compared with 62 deaths attributable to daily observed mean temperatures in 1990–2007. We found that a small portion of the estimated variability in mortality projections was due to the different simulations (i.e., variability due to characteristics of the simulations); most of the variability was associated with the RR used to calculate deaths attributable to temperature.

The uncertainty related to the temperature–mortality relationship had much more impact on heat-related mortality projections than did uncertainty due to climate models. This may be a result of the chosen temperature simulations. Although the simulations chosen covered an important part of the temperature variability from all simulations available, we did not include extreme simulations in our study (see Supplemental Material, Figure S3). Our choice of simulations may thus contribute to underestimating the contribution of the climate model projections. Furthermore, the uncertainty related to the temperature–mortality relationship may be large because of the propagation of the errors associated with the repeated use of the same mortality risk for a given temperature that occurs frequently; small changes in the risk function are magnified when applied repeatedly over numerous summer days. Future work is thus needed to better assess the uncertainty of the temperature–mortality relationship and climate model projections. Finally, we used a crude way to apportion the variation in the future ANs to different sources of variations, using only the sets of RRs that generated the lower and upper bounds of the 95% CI of the ANs and not the shape of the distribution of probable RRs. The explained variance by the sets of RRs thus represents an upper estimate.

Uncorrected simulated values of historical temperatures were quite different from the observed data. Therefore, to examine the variability in mortality projections, we corrected the simulated temperatures using the DT method. We used this method because it allows the validation of historical simulations on the basis of observed temperatures. Future work should seek to validate and improve different methods of correcting simulated temperatures in mortality projections, as is done in hydroclimatology studies (Teng et al. 2011).

To date, studies published on the near-future impact of summer temperatures on mortality in Montreal have reported an increase in deaths attributable to heat relative to current numbers of deaths in the summer (Cheng et al. 2008; Doyon et al. 2008; Martin et al. 2011). Although these studies used different methods to estimate future mortality, it is still possible to compare their results with our ours. Martin et al. (2011) calculated a predicted change in annual heat-related mortality rate per 100,000 population, and they estimated an increase of 152 heat-related attributable deaths per summer for the period 2031–2050, compared with 1981–2000. Doyon et al. (2008) calculated the equivalent of 79 summer heat-related deaths for 2020 (2% increase) and 81 deaths, for 2050 (6% increase). Cheng et al. (2008) estimated 96.3 summer heat-related deaths per year during 2040–2059. These results correspond to the range of our results. Regarding our results on factors affecting the variability in mortality projections, some differences are likely due to the RRs used to estimate the ANs of heat-related deaths. Lower numbers of summer heat-related deaths are likely in southern countries because heat thresholds are generally higher in communities closer to the equator (Hajat and Kosatsky 2010).

Methods for estimating mortality projections inevitably rely on assumptions concerning heat-related mortality relationships. In the present study, we chose to not take into account lag effects (Goldberg et al. 2011), harvesting, or mortality displacement (Zanobetti et al. 2002), nor did we consider intraurban variations of temperatures and risks (Smargiassi et al. 2009). Future studies should consider these aspects. We estimated only nonaccidental causes of death, and we did not conduct specific analyses for cardiovascular or respiratory causes of death (Goldberg et al. 2011; Halonen et al. 2011), even though the distribution of specific causes of death may vary in the future. Thus, further work should address these aspects. Our study also had other limitations: We assumed the same mean daily death count in the future, no demographic changes, and no population adaptation to heat, such as through access to air conditioning (Rogot et al. 1992). It is difficult to conclude what the impacts of these assumptions might be. On one hand, with demographic changes there are likely to be more vulnerable populations (elderly populations in particular); on the other hand, adaptations and mitigation measures may reduce climate impacts (Patz et al. 2008). New climate change impact studies taking into account these specific adaptation and mitigation measures should be performed. Humidity or dew point simulation for future years are needed to address the effect of these other climate factors.

Effects of climate change on health will affect most populations in the next decades, and put the lives and well-being of billions of people at increased risk. Our results suggest that the choice of the RR estimate for the association between temperature and mortality may be important in reducing uncertainty in mortality projections.


Armstrong B. 2006. Models for the relationship between ambient temperature and daily mortality. Epidemiology 17:624–631.

Cai W, Shi G, Cowan T, Bi D, Ribbe J. 2005. The response of the Southern Annular Mode, the East Australian Current, and the southern mid-latitude ocean circulation to global warming. Geoph Res Lett 32(23):L23706; doi: 10.1029/2005GL024701.

Campbell-Lendrum D, Woodruff R. 2007. Climate Change: Quantifying the Health Impact at National and Local Levels. WHO Environmental Burden of Disease Series No. 14. Geneva:World Health Organization. Available:​07/9789241595674_eng.pdf [accessed 23 October 2014].

Caya D, Laprise R, Giguère M, Bergeron G, Blanchet JP, Stocks B, et al. 1995. Description of the Canadian regional climate model. Water Air Soil Pollut 82:477–482.

Cheng CS, Campbell M, Li Q, Li G, Auld H, Day N, et al. 2008. Differential and combined impacts of extreme temperatures and air pollution on human mortality in south–central Canada. Part II: future estimates. Air Qual Atmos Health 1:223–235.

Christensen JH, Hewitson B, Busuioc A, Chen A, Gao X, Held I, et al. 2007. Regional climate projections. In: Climate Change, 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, et al., eds). Cambridge, UK:Cambridge University Press, 847–940. Available:​/ar4/wg1/ar4-wg1-chapter11.pdf [accessed 21 October 2014].

Costello A, Abbas M, Allen A, Ball S, Bell S, Bellamy R, et al. 2009. Managing the health effects of climate change: Lancet and University College London Institute for Global Health Commission. Lancet 373(9676):1693–1733.

de Elía R, Côté H. 2010. Climate and climate change sensitivity to model configuration in the Canadian RCM over North America. Meteorologische Zeitschrift 19:325–339.

Dequé M. 2007. Frequency of precipitation and temperature extremes over France in an anthropogenic scenario: model results and statistical correction according to observed values. Global Planetary Change 57:16–26.

Doyon B, Bélanger D, Gosselin P. 2008. The potential impact of climate change on annual and seasonal mortality for three cities in Quebec, Canada. Int J Health Geogr 7:23; doi: 10.1186/1476-072X-7-23.

Ebi KL. 2008. Healthy people 2100: modeling population health impacts of climate change. Clim Change 88:5–19.

Ebi KL, Gamble JL. 2005. Summary of a workshop on the development of health models and scenarios: strategies for the future. Environ Health Perspect 113:335–338; doi: 10.1289/ehp.7380.

Environment Canada. 2012. Climate. Available: [accessed 1 May 2012].

Flato GM, Boer GJ, Lee WG, McFarlane NA, Ramsden D, Reader MC, et al. 2000. The Canadian Centre for Climate Modelling and Analysis global coupled model and its climate. Clim Dyn 16:451–467.

Goldberg MS, Gasparrini A, Armstrong B, Valois MF. 2011. The short-term influence of temperature on daily mortality in the temperate climate of Montreal, Canada. Environ Res 111:853–860.

Gosling SN, Lowe JA, McGregor GR, Pelling M, Malamud BD. 2009. Associations between elevated atmospheric temperature and human mortality: a critical review of the literature. Clim Change 92:299–341.

Hajat S, Kosatky T. 2010. Heat-related mortality: a review and exploration of heterogeneity. J Epidemiol Community Health 64(9):753–760.

Hajat S, Kovats RS, Lachowycz K. 2007. Heat-related and cold-related deaths in England and Wales: who is at risk? Occup Environ Med 64:93–100.

Halonen JI, Zanobetti A, Sparrow D, Vokonas PS, Schwartz J. 2011. Relationship between outdoor temperature and blood pressure. Occup Environ Med 68:296–301.

Huang C, Barnett AG, Wang X, Vaneckova P, FitzGerald G, Tong S. 2011. Projecting future heat-related mortality under climate change scenarios: a systematic review. Environ Health Perspect 119:1681–1690; doi: 10.1289/ehp.1103456.

IPCC (Intergovernmental Panel on Climate Change). 2000. Special Report on Emissions Scenarios. Cambridge, UK:Cambridge University Press. Available:​sion/index.php?idp=0 [accessed 14 October 2014].

IPCC (Intergovernmental Panel on Climate Change). 2001. Climate Change 2001: Third Assessment Report. Available:​pcc_tar/ [accessed 14 October 2014].

IPCC (Intergovernmental Panel on Climate Change). 2007a. Climate Change 2007: Impacts, Adaptation and Vulnerability. Contribution of Working Group II to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change (Parry ML, Canziani OF, Palutikof JP, van der Linden PJ, Hanson CE, eds). Cambridge, UK:Cambridge University Press. Available:​/publications_ipcc_fourth_assessment_rep​ort_wg2_report_impacts_adaptation_and_vu​lnerability.htm [accessed 10 October 2014].

IPCC (Intergovernmental Panel on Climate Change). 2007b. Climate Change 2007: Synthesis Report: Geneva:IPCC. Available:​/ar4/syr/ar4_syr.pdf [accessed 10 October 2014].

Jungclaus JH, Keenlyside N, Botzet M, Haak H, Luo JJ, Latif M, et al. 2006. Ocean circulation and tropical variability in the coupled model ECHAM5/MPI-OM. J Climate 19:3952–3972.

Kovats RS, Hajat S. 2008. Heat stress and public health: a critical review. Annu Rev Public Health 29:41–55.

Laprise R. 2008. Regional climate modelling. J Comput Phys 227:3641–3666.

Li T, Horton RM, Kinney PL. 2013. Projections of seasonal patterns in temperature-related deaths for Manhattan, New York. Nat Clim Change 3:717–721.

Lobato I, Velasco C. 2004. A simple test of normality of time series. Econometric Theory 20:671–689.

Martin SL, Cakmak S, Hebbern CA, Avramescu ML, Tremblay N. 2011. Climate change and future temperature-related mortality in 15 Canadian cities. Int J Biometeorol 56(4):605–619.

McFarlane NA, Boer GJ, Blanchet JP, Lazare M. 1992. The Canadian Climate Centre second-generation general circulation model and its equilibrium climate. J Clim 5:1013–1044.

McMichael AJ, Woodruff RE, Hales S. 2006. Climate change and human health: present and future risks. Lancet 367:859–869.

Meehl GA, Covey C, Taylor KE, Delworth T, Stouffer RJ, Latif M, et al. 2007. The WCRP CMIP3 multimodel dataset: a new era in climate change research. Bull Am Meteor Soc 88:1383–1394.

Mpelasoka FS, Chiew FHS. 2009. Influence of rainfall scenario construction methods on runoff projections. J Hydrometeor 10:1168–1183.

O’Neill MS, Ebi KL. 2009. Temperature extremes and health: impacts of climate variability and change in the United States. J Occup Environ Med 51:13–25.

Paquin D. 2010. Évaluation du MRCC4 en passé récent (1961–1999). Climate Simulation Group, Rapport interne no 15. Montreal:Consortium Ouranos.

Patz J, Campbell-Lendrum D, Gibbs H, Woodruff R. 2008. Health impact assessment of global climate change: expanding on comparative risk assessment approaches for policy making. Annu Rev Public Health 29:27–39.

Patz JA, Campbell-Lendrum D, Holloway T, Foley JA. 2005. Impact of regional climate change on human health. Nature 438:310–317.

Peng RD, Bobb JF, Tebaldi C, McDaniel L, Bell ML, Dominici F. 2011. Toward a quantitative estimate of future heat wave mortality under global climate change. Environ Health Perspect 119:701–706; doi: 10.1289/ehp.1002430.

Plummer DA, Caya D, Frigon A, Côté H, Giguère M, Paquin D, et al. 2006. Climate and climate change over North America as simulated by the Canadian RCM. J Climate 19:3112–3132.

Rogot E, Sorlie PD, Backlund E. 1992. Air-conditioning and mortality in hot weather. Am J Epidemiol 136:106–116.

Smargiassi A, Goldberg MS, Plante C, Fournier M, Baudouin Y, Kosatsky T. 2009. Variation of daily warm season mortality as a function of micro-urban heat islands. J Epidemiol Community Health 63:659–664.

Teng J, Chiew F, Timbal B, Vaze J, Wang Y, Wang B, et al. 2011. Assessing historical and future runoff modelled using rainfall from the analogue downscaling method. In: MODSIM11, 19th International Congress on Modelling and Simulation, Perth, Australia, 12–16 December 2011 (Chan F, Marinova D, Andersson R, eds). Perth, Australia:Modelling and Simulation Society of Australia and New Zealand, 2789–2795. Available:​eng.pdf [accessed 14 October 2014].

Zanobetti A, Schwartz J, Samoli E, Gryparis A, Touloumi G, Atkinson R, et al. 2002. The temporal pattern of mortality responses to air pollution: a multicity assessment of mortality displacement. Epidemiology 13:87–93.

WP-Backgrounds Lite by InoPlugs Web Design and Juwelier Schönmann 1010 Wien