Next Article in Journal
Experimental Assessment of Recycled Diesel Spill-Contaminated Domestic Wastewater Treated by Reed Beds for Irrigation of Sweet Peppers
Next Article in Special Issue
Flood-Exposure is Associated with Higher Prevalence of Child Undernutrition in Rural Eastern India
Previous Article in Journal
Modified Ideal Cardiovascular Health Status is Associated with Lower Prevalence of Stroke in Rural Northeast China
Previous Article in Special Issue
Impacts of Climatic Variability on Vibrio parahaemolyticus Outbreaks in Taiwan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of an Early-Warning System for Heat Wave-Related Mortality in Europe: Implications for Sub-seasonal to Seasonal Forecasting and Climate Services

1
Catalan Institute of Climate Sciences (IC3), Carrer Doctor Trueta, 203, 3a, 08005 Barcelona, Spain
2
World Health Organization (WHO) Regional Office for Europe, European Centre for Environment and Health, Platz der Vereinten Nationen 1, 53113 Bonn, Germany
3
National Institute of Health and Medical Research, INSERM U988 and U1198, University of Montpelier, Building 24, Place Eugène Bataillon-CC105, 34095 Montpellier, Cedex 05, France
4
Ecole Pratique des Hautes Etudes, 75014 Paris, France
5
Division of Geriatrics, Department of Internal Medicine, Rehabilitation and Geriatrics, Geneva University Hospitals, University of Geneva, Chemin du., Pont-Bochet, 1226 Thônex, Switzerland
6
Catalan Institution for Research and Advanced Studies (ICREA), Passeig de Lluís Companys, 23, 08010 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2016, 13(2), 206; https://doi.org/10.3390/ijerph13020206
Submission received: 30 November 2015 / Accepted: 1 February 2016 / Published: 6 February 2016
(This article belongs to the Special Issue Climate Change and Human Health)

Abstract

:
Heat waves have been responsible for more fatalities in Europe over the past decades than any other extreme weather event. However, temperature-related illnesses and deaths are largely preventable. Reliable sub-seasonal-to-seasonal (S2S) climate forecasts of extreme temperatures could allow for better short-to-medium-term resource management within heat-health action plans, to protect vulnerable populations and ensure access to preventive measures well in advance. The objective of this study is to assess the extent to which S2S climate forecasts could be incorporated into heat-health action plans, to support timely public health decision-making ahead of imminent heat wave events in Europe. Forecasts of apparent temperature at different lead times (e.g., 1 day, 4 days, 8 days, up to 3 months) were used in a mortality model to produce probabilistic mortality forecasts up to several months ahead of the 2003 heat wave event in Europe. Results were compared to mortality predictions, inferred using observed apparent temperature data in the mortality model. In general, we found a decreasing transition in skill between excellent predictions when using observed temperature, to predictions with no skill when using forecast temperature with lead times greater than one week. However, even at lead-times up to three months, there were some regions in Spain and the United Kingdom where excess mortality was detected with some certainty. This suggests that in some areas of Europe, there is potential for S2S climate forecasts to be incorporated in localised heat–health action plans. In general, these results show that the performance of this climate service framework is not limited by the mortality model itself, but rather by the predictability of the climate variables, at S2S time scales, over Europe.

1. Introduction

Temperature-related illness and death is a continuing public health concern, and in many regions is presenting an increasing burden on public health systems due to climate change [1,2]. The Intergovernmental Panel on Climate Change (IPCC) concluded that there is medium confidence that the observed warming in the climate system has increased heat-related human mortality in some regions [3], and heat waves have been responsible for more fatalities in Europe over the past decades than any other extreme weather event [4].
Europe emerges as an especially responsive area to temperature rise under climate change, particularly during the warm season [5]. The anomalous and persistent heat of the summer of 2003 caused over 70,000 additional deaths across twelve countries in Western Europe [6]. This event revealed the lack of reactivity of society and health services to such an extreme event. Since then, some European countries have implemented strict measures and protocols to minimize the negative effects of heat waves. Some of these efforts have successfully reduced the impact of the following extreme events, such as in France, where mortality rates were reduced by 70% lower than expected during the 2006 episode with regard to the expected mortality [7]. Certain adverse health effects could easily be avoided if informed decisions were made prior to heat waves, to protect vulnerable populations, such as children and the elderly, and ensure access to preventive measures well in advance [8,9].
A WHO study [10] estimated an increase in annual heat-related deaths of 15,280 additional deaths in 2030, and 32,152 additional deaths in 2050, for the over 65 year-old population in the WHO European Region (including central Asia), assuming no adaptation. However, this increase might be offset by a similar reduction in cold-related human mortality before 2050, whereas such compensation would not take place afterwards, with a resulting average of 15,000 excess deaths per year [11]. Moreover, heat waves are relatively infrequent compared with the total number of hot days, thus the effect of heat waves is likely to represent only a fraction of the total heat impact [12]. Deaths attributable to extreme heat are roughly as frequent as those attributable to moderate heat [13]. The effects of exposure can be directly heat-related (heat stroke, heat fatigue and dehydration) or heat stress, which can contribute to a worsening of respiratory and cardiovascular diseases, electrolyte disorders and kidney problems [8]. A myriad of other indirect health effects should also be considered, including food shortages and malnutrition, social conflicts and aggravation of chronic illnesses.
Temperature-related illnesses and deaths are largely preventable [7,14,15]. One of the most effective health system preparations for this emergency is the development and implementation of action plans for preparedness and response, such as heat-health action plans (HHAP). Development, implementation and subsequent evaluation and improvement of these preparedness and response plans would lead to a reduction in temperature-related mortality and stronger climate-resilient health systems. A HHAP relies on early-warning systems for timely activation and to allow for longer-term resource planning. Forecast data and projections for temperature are primarily used by health systems within the decision-making frameworks of these plans. Improved climate forecasts would allow for better short-to-medium-term resource management within health systems and would help authorities prepare and respond ahead of heat waves.
HHAPs can be evaluated based on inclusion of nine core elements [16,17,18]. These elements focus on short-term activation and implementation during the event, defining specific roles and actors. Three of these core elements are of specific interest for the development of climate services for temperature and health: a timely alert system, preparedness of the health/social care system, and long-term urban planning. Due to their speed of implementation, these three core elements have potentially different lead-times for forecasting: alert systems triggering plan activation are typically days in advance. However, if a country had not experienced such an extreme temperature event for several years, earlier warning could provide sufficient time to review and update existing action plans. Longer-term health system preparedness and better resource management would benefit from sub-seasonal-to-seasonal (S2S) forecasts to allow for increasing capacity in the health system and improving the resilience of the primary health services, especially in the case of prolonged and/or recurrent heat waves where mortality significantly increases with intensity and duration of heat waves. Long-term urban planning for extreme weather occurs on annual-to-decadal (or even longer timescales) and includes cross-sectoral issues such as increasing green and blue spaces, changes in building design and improvement in housing stock, changes in land-use decisions, reductions in energy consumption, and improved transport policies.
In an assessment of the 18 countries of the WHO European Region with known heat–health action plans, 16 had a clearly defined alert system and a health system preparedness component, all included an information plan, but only four included long-term urban planning within the HHAP itself [16]. None of these HHAPs routinely incorporated S2S climate forecasts. Thus, there exists potential for these forecasts to be incorporated into pre-existing plans in the European Region, whilst countries yet to develop a HHAP could incorporate such information at the design stage.
In a previous study [19], the performance of a climate-driven mortality model to provide probabilistic predictions of exceeding emergency mortality thresholds for a heat wave and a cold spell scenario was evaluated. The mortality model was formulated using observed (reanalysis) apparent temperature data. This observed temperature data was then used to produce spatio-temporal probabilistic mortality estimates for a heat wave and cold spell scenario. The model showed considerable skill, particularly for the heat wave scenario (1–15 August 2003), successfully anticipating the occurrence or non-occurrence of mortality rates exceeding the emergency threshold (75th percentile of the mortality distribution) for 89% of the 54 regions, given a probability decision threshold of 70%. The use of observed apparent temperature in the mortality model represented an upper bound to forecast skill, given a perfect climate forecast. In this study, we replace the observed apparent temperature with forecast apparent temperature at 8 different lead times, from 1 day to 3 months, to drive the model and produce probabilistic mortality forecasts ahead of the 2003 heat wave event in Europe. This allows the assessment of the extent to which S2S climate forecasts could be incorporated into HHAP, to support timely public health decision-making ahead of imminent extreme temperature events in Europe.

2. Experimental Section

2.1. Temperature-Mortality Model

Daily mortality data corresponding to 187 NUTS2 regions across 16 countries in Europe were obtained from 1998 to 2003 (see [6] for details). The NUTS2 regions (i.e., second level of the Nomenclature of Territorial Units for Statistics) are basic regions of the economic territory of the European Union for the application of regional policies. Due to the extremely large differences in population between these regions (from 18 million in North Rhine-Westphalia to less than 100,000 in some Slovenian regions), regions were grouped in 54 larger and more homogeneously populated areas. Several factors were taken into account in this process, such as geographical proximity, state borders (e.g., similar warning systems or adaptation measures might be applied in each country), and similarities in regional temperature/mortality dependencies or comfort temperatures [11]. Location-specific average mortality rates, at given temperature intervals over the entire time period, were modelled to account for the increased mortality observed during both high and low temperature extremes and differing comfort temperatures between the 54 aggregations [11,19]. The temperature-mortality dependency for each aggregation was estimated as follows:
Apparent temperature, defined by the following equation, is the climatological input to the mortality model:
T a p p = 2.653 + 0.994 T a i r + 0.0153 T d e w p t 2
where Tapp is the apparent temperature, Tair the air temperature and Tdewpt the dew point temperatures (at 2 metres), in degrees Celsius, from the reanalysis ERA-Interim dataset [20]. We refer to this reanalysis data as “observed” apparent temperature data herein. To relate the climate variables to mortality data for the 54 aggregated regions, the apparent temperature data at gridpoints found inside the region were identified and averaged for each time step. In case the region was smaller than the gridsquares, the value of the nearest neighbour grid to the centroid of the region was used.
The range of apparent temperatures was divided in equally spaced intervals. Days belonging to each interval were grouped, and daily temperature and mortality data within each interval were averaged. Interval mean mortality was smoothed using a centred 31-term filter, corresponding to nearly 3 °C intervals. The lowest value defines the interval of comfort temperature. This threshold divides the range of temperatures into “warm” and “cold” tails [11]. The model used to fit the temperature-mortality curves was formulated as follows:
y i k N ( α j + β 1 j x i k + β 2 j x 2 i k + β 3 j x 3 i k , σ j 2 )
where yik is the logarithm of the average mortality rate (per million population) in region i, and at temperature interval k. Then, for each region i, the log mortality rate was formulated as a non-linear function of temperature, xik, (a third order polynomial) with location specific intercept αj. This formulation directly corresponds to the linear predictor of a Poisson count model, with the population in each region as a model offset. When comparing observed to predicted mortality in the subsequent analysis, the log mortality incidence rates are transformed to mortality counts. Note, parameters are fitted separately for the warm tail (j = w) and cold tail (j = c), depending on whether the temperature is greater than (xikxim) or less than (xik < xim) the comfort temperature (i.e., the temperature of minimum mortality), xim. The comfort temperature occurs twice a year around June and September, defining a summer season of warm tail temperatures with non-linear mortality sensitivity and a long three-season period (autumn, winter and spring) of cold tail temperatures with near-linear mortality increases with decreasing temperature [11]. This model framework distinguishes between warm and cold periods, without restricting the model to calendar months, thus allowing temperature extremes, which might occur outside the typical range of “warm” or “cold” months, to be detected.
The model was fitted in a Bayesian probabilistic framework, with 1000 samples of mortality rates generated for each time step [19]. This allows the simulation of probabilistic predictions of daily mortality in space and time, given any user-defined emergency and probability decision thresholds. In order to simulate mortality predictions for a heat wave scenario, spatio-temporal observed apparent temperature data, xit, where t is the time step (daily), are combined with 1000 samples of the parameters estimated from the warm tail model (j = w):
y i t N ( α w + β 1 w x i t + β 2 w x 2 i t + β 3 w x 3 i t , σ w 2 ) ,   if   x i t x i m
This gives 1000 samples of daily mortality rates for each region and each day, yit. Daily mortality rates can then be averaged for the climatological events of interest (e.g., the heat wave in August 2003).

2.2. Climate Forecasts

Climate forecasts from two different datasets were used: The European Centre for Medium-Range Weather Forecasts (ECMWF) seasonal forecast system version 4 (System4) and the ECMWF sub-seasonal forecast system. Hindcasts (i.e., retrospective forecasts) for the 2003 event were obtained using the most recent state-of-the-art forecast system.
Seasonal forecasts seek to take advantage of aspects of the climate system with long-term memory, such as the oceans, to predict climate anomalies one or more months ahead of a given season. To estimate uncertainty, each forecast consists of an ensemble of forecasts, obtained by perturbing the initial conditions. Therefore, the information produced can be summarised in probabilistic terms. The ECMWF System4 [21] consists of a global coupled ocean-atmosphere model initiated from reanalysis. The available hindcasts cover 30 years (1981–2010). In the present study two seasonal forecasts, with 15 ensemble members each, were used, one starting 1 July 2003 (one month before the event), the other starting 1 May 2003 (three months before the event).
Sub-seasonal forecasts (or extended-range forecasts) lie between medium-range weather forecasts and seasonal forecasts, and take advantage of the predictability of phenomena like the Madden-Julian oscillation or jet stream blocking events [22]. The climate model used for the sub-seasonal forecasts is similar to the model used to produce seasonal forecasts, but with a higher resolution. The model is run twice a week, out to 46 days ahead. The available hindcasts cover 20 years (1995–2015), with 10 ensemble members. In the present study, 6 sub-seasonal forecasts were used with the following start dates: 30-07-2003, 27-07-2003, 23-07-2003, 20-07-2003, 16-07-2003 and 13-07-2003, which correspond to lead times of 1, 4, 8, 11, 15 and 18 days, respectively.
Both sets of forecasts are produced using multiple climate model runs, which produce high frequency (up to hourly) data. The difference between sub-seasonal and seasonal forecasts is the complexity of the model and the lead time: the longer the lead time, the longer the period over which variables should be averaged, in order to find predictability. However, forecast providers usually supply data for every day in the forecast period, to enable the study of indicators, such as the number of heat wave days in a given month or season. In our case, the mortality model needs daily frequency data as input. For this study, daily mortality forecasts are aggregated over a longer period of 15 days (see below).
As with the observed temperature data, both sub-seasonal and seasonal forecasts were aggregated to the 54 regions, as described above. Therefore, 54 time series with a daily resolution, for each ensemble member and forecast lead time, were computed. Both series for 2-metre temperature and dew point temperature were generated and then combined, following Equation (1), to produce regional sub-seasonal and seasonal forecasts of apparent temperature.
The regional apparent temperature forecasts were then bias-corrected, following [23], before using them as an input to the mortality model. This technique removes “model drift” [24] and aligns forecasts with past observations. Both sub-seasonal and seasonal apparent temperature forecast ensembles data were then used to run the mortality model, generating 1000 samples of daily mortality rates for each region, each day and each ensemble member (i.e., 10,000 samples for the sub-seasonal forecasts and 15,000 samples for the seasonal forecasts, for each day and region). By incorporating the whole ensemble of temperature forecasts in the probabilistic Bayesian framework, the additional source of uncertainty generated by using forecast instead of observed temperatures is accounted for, as well as uncertainty in the mortality estimates themselves.
For comparison, observed apparent temperature (from ERA-Interim reanalysis) was used to run the mortality model, producing 1000 samples of daily mortality rates for each region and each day. This allows an assessment of skill loss when replacing observed with forecast climate data to run the mortality model. As a case study, daily mortality rate samples were averaged for the heat wave period of interest, 1–15 August 2003, using both observed and S2S forecast climate data as inputs to the mortality model.
Following [19], region specific emergency thresholds of mortality rates were set at the 75th percentile (3rd quartile) of the mortality distribution, for the days in the time period in which temperatures were greater than the comfort temperature, i.e., in the warm tail distribution. This allowed us to assess the ability of the temperature-driven mortality model to determine mortality rates exceeding this emergency threshold. The probability of exceeding this threshold in each region was determined using the mortality model driven by (i) observed apparent temperature data (i.e., the proportion of 1000 samples that exceeded the location specific emergency threshold); (ii) the sub-seasonal apparent temperature forecast data, with lead times of 1, 4, 8, 11, 15 and 18 days (i.e., the proportion of 10,000 samples that exceeded the location specific emergency threshold) and (iii) the seasonal apparent temperature forecast data, with lead times of 1 and 3 months (i.e., the proportion of 15,000 samples that exceeded the location specific emergency threshold).

3. Results

3.1. Excess Mortality Probability Maps Using Climate Forecasts at Increasing Lead Times

Figure 1a shows the predicted probability of mortality rates exceeding the 75th percentile of the mortality distribution for the heat wave period 1–15 August 2003. The corresponding observations (i.e., whether the mortality rate exceeded the threshold or not) are displayed in Figure 1b. There is considerable agreement between predictions and observations, with the model correctly predicting with high confidence that mortality would exceed the emergency threshold across most of Spain, France and Northern Italy (see [19] for further details). We use this example of a mortality prediction, using observed temperature, as a reference (or benchmark). This allows us to compare mortality predictions, using forecasts of apparent temperature at increasing lead times from sub-seasonal (1 to 18 days, see Figure 2) to seasonal (1 and 3 months, see Figure 3) time scales.
Figure 1. (a) Probabilistic map of exceeding emergency daily mortality threshold (75th percentile of daily mortality distribution in the warm tail) using a mortality model driven with observed apparent temperature data during the heat wave scenario (115 August 2003); (b) Corresponding observations for the same period. The graduated colour bar represents the probability of exceeding the mortality threshold (ranging from 0%, pale colours, to 100%, deep colours). Source: [19].
Figure 1. (a) Probabilistic map of exceeding emergency daily mortality threshold (75th percentile of daily mortality distribution in the warm tail) using a mortality model driven with observed apparent temperature data during the heat wave scenario (115 August 2003); (b) Corresponding observations for the same period. The graduated colour bar represents the probability of exceeding the mortality threshold (ranging from 0%, pale colours, to 100%, deep colours). Source: [19].
Ijerph 13 00206 g001
Figure 2. Probabilistic maps of exceeding emergency daily mortality threshold (75th percentile of daily mortality distribution in the warm tail) during the heat wave scenario (1–15 August 2003) using sub-seasonal forecasts of apparent temperature at lead times of (a) 1 day; (b) 4 days; (c) 8 days; (d) 11 days; (e) 15 days; (f) 18 days. The graduated colour bar represents the probability of exceeding the mortality threshold (ranging from 0%, pale colours, to 100%, deep colours).
Figure 2. Probabilistic maps of exceeding emergency daily mortality threshold (75th percentile of daily mortality distribution in the warm tail) during the heat wave scenario (1–15 August 2003) using sub-seasonal forecasts of apparent temperature at lead times of (a) 1 day; (b) 4 days; (c) 8 days; (d) 11 days; (e) 15 days; (f) 18 days. The graduated colour bar represents the probability of exceeding the mortality threshold (ranging from 0%, pale colours, to 100%, deep colours).
Ijerph 13 00206 g002aIjerph 13 00206 g002b
Figure 3. Probabilistic maps of exceeding emergency daily mortality threshold (75th percentile of daily mortality distribution in the warm tail) during the heat wave scenario (1–15 August 2003) using seasonal forecast apparent temperature at lead times of (a) 1 month and (b) 3 months. The graduated colour bar represents the probability of exceeding the mortality threshold (ranging from 0%, pale colours, to 100%, deep colours).
Figure 3. Probabilistic maps of exceeding emergency daily mortality threshold (75th percentile of daily mortality distribution in the warm tail) during the heat wave scenario (1–15 August 2003) using seasonal forecast apparent temperature at lead times of (a) 1 month and (b) 3 months. The graduated colour bar represents the probability of exceeding the mortality threshold (ranging from 0%, pale colours, to 100%, deep colours).
Ijerph 13 00206 g003aIjerph 13 00206 g003b
At lead times of 1 and 4 days, excess mortality is reasonably well predicted, with high certainty of exceeding the emergency threshold for much of western Europe (compare Figure 2a,b to Figure 1b). As lead time increases, the mortality forecasts become less certain and the sub-seasonal climate forecasts begin to erroneously detect excess mortality in southern Italy (compare Figure 2c–f) with Figure 1b). Interestingly, even at seasonal time scales, excess mortality is correctly predicted with some degree of certainty for some regions, including central and eastern Spain and southwestern United Kingdom (compare Figure 3a,b to Figure 1b). These results are promising for the development of more localised early warning systems.
Figure 4. ROC curves for the binary event of exceeding the emergency mortality threshold of the 75th percentile in each of the 54 regions for the heat wave scenario (115 August 2003), using the probabilistic mortality model driven by forecast apparent temperature data at lead times ranging from 1 day, 4 days, 8 days, 11 days, 15 days, 18 days, 1 month and 3 months. The ROC curve for the mortality model driven by observed apparent temperature data is shown for reference (black curve).
Figure 4. ROC curves for the binary event of exceeding the emergency mortality threshold of the 75th percentile in each of the 54 regions for the heat wave scenario (115 August 2003), using the probabilistic mortality model driven by forecast apparent temperature data at lead times ranging from 1 day, 4 days, 8 days, 11 days, 15 days, 18 days, 1 month and 3 months. The ROC curve for the mortality model driven by observed apparent temperature data is shown for reference (black curve).
Ijerph 13 00206 g004

3.2. Skill Assessment for Increasing Forecast Lead Time

Relative operating characteristic (ROC) curve can be used to indicate the hit rates and false alarm rates that would result from using different probability decision thresholds to determine a binary event (in this case, exceeding or not the emergency mortality threshold in regions across Europe). The ROC score or AUC (equivalent to the Area Under the modelled ROC curves), is a widely used measure of skill [25]. The ROC score characterises the quality of a forecast system by describing the system’s ability to anticipate correctly the occurrence or non-occurrence of pre-defined events. A ROC score value of 0.5 indicates zero skill while a value of 1 repr in this case exceeding esents perfect skill.
Figure 5. The ROC score or area under the ROC curve (AUC) and associated 95% confidence intervals for the binary event of exceeding the emergency mortality threshold of the 75th percentile in each of the 54 regions for the heat wave scenario (115 August 2003), using the probabilistic mortality model driven by forecast apparent temperature data at lead times ranging from 1 day, 4 days, 8 days, 11 days, 15 days, 18 days (sub-seasonal forecasts), 1 month and 3 months (seasonal forecasts). The AUC for the mortality model driven by observed apparent temperature data is shown for reference. The dotted horizontal line indicates the skill threshold of AUC = 0.5. The red vertical line distinguishes lead-times at sub-seasonal (less than a month) and seasonal (greater than a month) time scales. Results are also displayed in the accompanying table.
Figure 5. The ROC score or area under the ROC curve (AUC) and associated 95% confidence intervals for the binary event of exceeding the emergency mortality threshold of the 75th percentile in each of the 54 regions for the heat wave scenario (115 August 2003), using the probabilistic mortality model driven by forecast apparent temperature data at lead times ranging from 1 day, 4 days, 8 days, 11 days, 15 days, 18 days (sub-seasonal forecasts), 1 month and 3 months (seasonal forecasts). The AUC for the mortality model driven by observed apparent temperature data is shown for reference. The dotted horizontal line indicates the skill threshold of AUC = 0.5. The red vertical line distinguishes lead-times at sub-seasonal (less than a month) and seasonal (greater than a month) time scales. Results are also displayed in the accompanying table.
Ijerph 13 00206 g005

Climate Forecast Lead TimeMean AUC (95% Confidence Intervals)
Observed0.97 (0.93, 1.00)
1 day0.81 (0.73, 0.88)
4 days0.76 (0.67, 0.84)
8 days0.43 (0.37, 0.5)
11 days0.21 (0.14, 0.29)
15 days0.47 (0.41, 0.52)
18 days0.21 (0.13, 0.29)
1 month0.22 (0.14, 0.29)
3 months0.36 (0.29, 0.43)
Figure 4 shows ROC curves for the probabilistic mortality predictions for the 2003 heat wave scenario, using forecast temperatures at increasing lead times, ranging from 1 day to 3 months. The ROC curve obtained using the observed temperature data is also shown. Figure 5 shows the corresponding ROC score or AUC (i.e., the area under the modelled ROC curve) for each lead time, and associated confidence intervals, computed with 2000 stratified bootstrap replicates [26]. Results show that significant overall skill exists for the two shorter lead times (1 and 4 days). This is an important finding, since the non-linear nature of the mortality model makes it very sensitive to errors in temperature, especially in the warm tail. These mortality forecasts are more informative than simple temperature forecasts because they use empirical temperature-mortality relationships that take into account the different degrees of adaptation of the regions to warm temperatures. AUC for lead times longer than 4 days are below the reference value of 0.5. In the case of the ROC score, this implies that the forecast is no better than the mortality prediction derived from a climatological or random climate forecast. The “negative” skill found using longer lead times is the result of the climate forecasts not being able to reproduce the record-breaking summer heat wave, and instead, simulating a local warm event in southern Italy. Therefore, excess mortality is incorrectly predicted in southern Italy, leading to potential false alarms.

4. Discussion

In this study, we have used S2S climate forecast data at different lead times (i.e., 1, 4, 8, 11, 15 and 18 days, and 1 and 3 months) to drive a mortality model [11,19], and compared results to the predictions inferred from observational climate data [19]. Results have shown that there is a smooth and monotonically decreasing transition in skill between excellent predictions with observational data and predictions with almost no skill a week in advance and beyond. These results show that the performance of the scheme is not limited by the mortality model itself, but is instead constrained by the skill of the climate forecasting system. However, even at longer lead times, the system is able to forecast excess mortality with some confidence in certain regions of Spain (see Figure 2 and Figure 3). These regions coincide with the regions were skill is found in the seasonal forecast of summer months [21]. While this may be a coincidence, results suggest that forecasts with longer lead times are possible in these regions. This needs to be confirmed in further studies, using temperature and mortality data over a longer time period.
This simplistic mortality model is formulated using only one thermal variable. Other socio-economic factors, which can strongly influence mortality rates during heat waves, are not explicitly considered in the model. However, the temperature-mortality relationships and the emergency thresholds are region dependent. Therefore, they implicitly include different sensitivity to temperatures in each region. Despite its simplicity, results show that the model successfully predicts excess mortality for the 2003 heat wave event when observed temperature data is used. Predictability is only lost when replacing observed with forecast temperature. Thus, the lack of skill at increasing lead times is due to the quality of the climate forecasts, rather than the temperature-mortality model itself.
This work indicates that the predictability of temperature-related mortality is strictly determined by the chaotic nature of the atmosphere. Indeed, temperature-related mortality is predictable as far as climate forecasts provide useful skill for the prediction of temperature variables. Beyond this window of around one week, the limited skill in climate forecasts prevents the accurate prediction of excess mortality at lead times up to a season ahead. Therefore, it is the lack of climate forecast skill in Europe that limits the predictability of temperature-related mortality for the majority of regions. However, there appears to be some predictability of temperature-related mortality in parts of Spain and the United Kingdom, even at lead times of 3 months. These areas deserve further attention to investigate the potential for more localised early-warning systems.
The spatial coverage of the mortality dataset allowed a Europe-wide assessment of the capability of state-of-the-art seasonal climate forecast to predict excess mortality in the region. However, this study is restricted by the relatively short time period for which the mortality data was available (1998–2003), thus providing only one example of an extreme heat wave for analysis. As more mortality data becomes available, an assessment of the loss of skill in predicting excess mortality with increasing forecast lead-time, for other more recent European heat waves, can be compared to corroborate the findings presented here.
Despite the limitations of this study, the results have strong implications for current research initiatives in Europe. As it is widely known, the Horizon 2020 programme of the European Commission is highlighting the importance of climate services for Europe and globally within the framework of European research and excellence. Particularly, its specific challenge is the provision of trustworthy science-based climate information to government, public and private decision-makers as a fundamental prerequisite for both properly managing the risks society is facing and seizing the opportunities this implies. However, if climate information is to improve the decision-making process and the resilience of society in Europe, climate services need to be based on skillful climate forecasts that can provide useful information as input data for the impact models. In that regard, our work highlights the current potential weaknesses of the climate service strategy for Europe, given that seasonal forecasts in this region are limited for informing certain impact models. Therefore, if climate services initiatives are to succeed in achieving their challenging goals, efforts should in parallel be directed towards the improvement of climate forecasting in this region of the planet. The limited skill of seasonal forecasts in Europe is linked to the low predictability of the planetary waves developing in the polar jet stream. However, several encouraging results have been published pointing to near-future improvements in seasonal forecasts [27,28], while sub-seasonal forecasts are consistently being improved each year [29].

5. Conclusions

In certain areas of Europe, there is potential for longer lead-times to be incorporated into pre-existing HHAP in the European Region, whilst countries yet to develop a HHAP could incorporate such information at the design stage. However, our results indicate that a compromise will have to be reached between user needs and the capabilities of seasonal climate forecasts over Europe, to provide skillful mortality predictions in advance of imminent extreme temperature events.

Acknowledgments

The research leading to these results has received funding from the EUPORIAS project (grant agreement No. 308291) funded by the European Commission’s Seventh Framework Research Programme. The mortality data collection was supported by the EU Community Action Programme for Public Health (grant agreement No. 2005114). Joan Ballester gratefully acknowledges funding from the European Commission through a Marie Curie International Outgoing Fellowship (project MEMENTO from the FP7-PEOPLE-2011-IOF call), and from the European Commission and the Catalan Government through a Marie Curie—Beatriu de Pinós Fellowship (project 00068 from the BP-DGR-2014-B call). We acknowledge ECMWF for the ERA-Interim and sub-seasonal forecast data, which were downloaded from the public data portal (http://apps.ecmwf.int/datasets/). System4 seasonal forecasts were downloaded from the ECOMS-UDG portal (https://meteo.unican.es/trac/wiki/udg/ecoms). We are grateful to colleagues at the Federal Office of Meteorology and Climatology MeteoSwiss in the EUPORIAS project for helpful feedback on an earlier version of this study.

Author Contributions

Rachel Lowe took the leading role in the design of the study, collation of the data and devised the Bayesian model framework. Markel García-Díez obtained the seasonal climate forecast data and performed bias corrections, helped design the study and interpret the results. Joan Ballester formulated the temperature-mortality curves, which were used to fit the models, helped design the study and interpret the results. James Creswick provided expert advice on public health impact and user perspective and contributed to the introduction and discussion. Jean-Marie Robine and François R. Herrmann provided the mortality data and expert advice on its interpretation. Xavier Rodó participated in the design of the study and interpretation of results. All authors wrote and edited the manuscript. James Creswick is a staff member of the World Health Organization (WHO) Regional Office for Europe. The authors alone are responsible for the views expressed in this publication and they do not necessarily represent the decision or stated policy of the World Health Organization.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Patz, J.A.; Campbell-Lendrum, D.; Holloway, T.; Foley, J.A. Impact of regional climate change on human health. Nature 2005, 438, 310–317. [Google Scholar] [CrossRef] [PubMed]
  2. Ostro, B.; Barrera-Gómez, J.; Ballester, J.; Basagaña, X.; Sunyer, J. The impact of future summer temperature on public health in Barcelona and Catalonia, Spain. Int. J. Biometeorol. 2012, 56, 1135–1144. [Google Scholar] [CrossRef] [PubMed]
  3. Pachauri, R.K.; Meyer, L. IPCC Climate Change 2014: Synthesis Report; Intergovernmental Panel on Climate Change: Geneva, Switzerland, 2014. [Google Scholar]
  4. European Environment Agency (EEA). Climate Change, Impacts and Vulnerability in Europe 2012; European Environment Agency: Copenhagen, Denmark, 2012. [Google Scholar]
  5. Giorgi, F. Climate change hot-spots. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef]
  6. Robine, J.-M.; Cheung, S.L.K.; Roy, S.L.; Oyen, H.V.; Griffiths, C.; Michel, J.-P.; Herrmann, F.R. Death toll exceeded 70,000 in Europe during the summer of 2003. C. R. Biol. 2008, 331, 171–178. [Google Scholar] [CrossRef] [PubMed]
  7. Fouillet, A.; Rey, G.; Wagner, V.; Laaidi, K.; Empereur-Bissonnet, P.; Le Tertre, A.; Frayssinet, P.; Bessemoulin, P.; Laurent, F.; De Crouy-Chanel, P.; et al. Has the impact of heat waves on mortality changed in France since the European heat wave of summer 2003? A study of the 2006 heat wave. Int. J. Epidemiol. 2008, 37, 309–317. [Google Scholar] [CrossRef] [PubMed]
  8. Aström, C.; Orru, H.; Rocklöv, J.; Strandberg, G.; Ebi, K.L.; Forsberg, B. Heat-related respiratory hospital admissions in Europe in a changing climate: A health impact assessment. BMJ Open 2013, 3, e001842. [Google Scholar] [CrossRef] [PubMed]
  9. Hajat, S.; Barnard, L.T.; Butler, C. Heat-related and cold-related mortality and morbidity. In Climate Change and Global Health; CABI International: Wallingford, UK, 2014; p. 21. [Google Scholar]
  10. World Health Organization. Quantitative Risk Assessment of the Effects of Climate Change on Selected Causes of Death, 2030s and 2050s; WHO: Geneva, Switzerland, 2014. [Google Scholar]
  11. Ballester, J.; Robine, J.-M.; Herrmann, F.R.; Rodó, X. Long-term projections and acclimatization scenarios of temperature-related mortality in Europe. Nat. Commun. 2011. [Google Scholar] [CrossRef] [PubMed]
  12. Smith, K.R.; Woodward, A.; Campbell-Lendrum, D.; Chadee, D.; Honda, Y.; Liu, Q.; Olwoch, J.M.; Revich, B.; Sauerborn, R. Human health: Impacts, adaptation, and co-benefits. In Climate Change 2014: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel of Climate Change; Field, C.B., Barros, V.R., Dokken, D.J., Mach, K.J., Mastrandrea, M.D., Bilir, T.E., Chatterjee, M., Ebi, K.L., Estrada, Y.O., Genova, R.C., Girma, B., Kissel, E.S., Levy, A.N., MacCracken, S., Mastrandrea, P.R., White, L.L., Eds.; Cambridge University Press: Cambridge, UK, 2014; pp. 709–754. [Google Scholar]
  13. Gasparrini, A.; Guo, Y.; Hashizume, M.; Lavigne, E.; Zanobetti, A.; Schwartz, J.; Tobias, A.; Tong, S.; Rocklöv, J.; Forsberg, B.; et al. Mortality risk attributable to high and low ambient temperature: A multicountry observational study. Lancet 2015, 386, 369–375. [Google Scholar] [CrossRef]
  14. Donaldson, G.; Ermakov, S.; Komarov, Y.; McDonald, C.; Keatinge, W. Cold related mortalities and protection against cold in Yakutsk, eastern Siberia: Observation and interview study. BMJ 1998, 317, 978–982. [Google Scholar] [CrossRef] [PubMed]
  15. Bobb, J.F.; Peng, R.D.; Bell, M.L.; Dominici, F. Heat-related mortality and adaptation to heat in the United States. Environ. Health Perspect. 2014, 122. [Google Scholar] [CrossRef] [PubMed]
  16. Bittner, M.I.; Matthies, E.F.; Dalbokova, D.; Menne, B. Are European countries prepared for the next big heat-wave? Eur. J. Public Health 2014, 24, 615–619. [Google Scholar] [CrossRef] [PubMed]
  17. McGregor, G.R.; Bessemoulin, P.; Ebi, K.L.; Menne, B. Heatwaves and Health: Guidance on Warning-System Development; World Meteorological Organization and World Health Organization: Geneva, Switzerland, 2015. [Google Scholar]
  18. Matthies, F.; Bickler, G.; Marin, N.; Hales, S. WHO Europe, Heat–Health Action Plans; Matthies, F., Bickler, G., Marin, N., Hales, S., Eds.; WHO Regional Office for Europe: Copenhagen, Denmark, 2008. [Google Scholar]
  19. Lowe, R.; Ballester, J.; Creswick, J.; Robine, J.-M.; Herrmann, F.R.; Rodó, X. Evaluating the performance of a climate-driven mortality model during heat waves and cold spells in Europe. Int. J. Environ. Res. Public Health 2015, 12, 1279–1294. [Google Scholar] [CrossRef] [PubMed]
  20. Dee, D.; Uppala, S.; Simmons, A.; Berrisford, P.; Poli, P.; Kobayashi, S.; Andrae, U.; Balmaseda, M.; Balsamo, G.; Bauer, P.; et al. The ERA-Interim reanalysis: Configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc. 2011, 137, 553–597. [Google Scholar] [CrossRef]
  21. Molteni, F.; Stockdale, T.; Balmaseda, M.; Balsamo, G.; Buizza, R.; Ferranti, L.; Magnusson, L.; Mogensen, K.; Palmer, T.; Vitart, F. The New ECMWF Seasonal Forecast System (System 4); European Centre for Medium-Range Weather Forecasts: Reading, UK, 2011. [Google Scholar]
  22. Vitart, F. Monthly forecasting at ECMWF. Mon. Weather Rev. 2004, 132, 2761–2779. [Google Scholar] [CrossRef]
  23. Mahlstein, I.; Spirig, C.; Liniger, M.A.; Appenzeller, C. Estimating daily climatologies for climate indices derived from climate model data and observations. J. Geophys. Res. Atmos. 2015, 120, 2808–2818. [Google Scholar] [CrossRef] [PubMed]
  24. Manzanas, R.; Fernández, J.; Magariño, M.E.; Gutiérrez, J.M.; José Doblas-Reyes, F.; Nikulin, G.; Buontempo, C. Assessing the Drift of Seasonal Forecasts. In EGU General Assembly Conference Abstracts; Copernicus Publications: Göttingen, Germany, 2014; Volume 16, p. 15360. [Google Scholar]
  25. Mason, S.J.; Graham, N.E. Areas beneath the relative operating characteristics (ROC) and relative operating levels (ROL) curves: Statistical significance and interpretation. Q. J. R. Meteorol. Soc. 2002, 128, 2145–2166. [Google Scholar] [CrossRef]
  26. Carpenter, J.; Bithell, J. Bootstrap confidence intervals: When, which, what? A practical guide for medical statisticians. Stat. Med. 2000, 19, 1141–1164. [Google Scholar] [CrossRef]
  27. Scaife, A.; Arribas, A.; Blockley, E.; Brookshaw, A.; Clark, R.; Dunstone, N.; Eade, R.; Fereday, D.; Folland, C.; Gordon, M.; et al. Skillful long-range prediction of European and North American winters. Geophys. Res. Lett. 2014, 41, 2514–2519. [Google Scholar] [CrossRef]
  28. Prodhomme, C.; Doblas-Reyes, F.; Bellprat, O.; Dutra, E. Impact of land-surface initialization on sub-seasonal to seasonal forecasts over Europe. Clim. Dyn. 2015, 1–17. [Google Scholar] [CrossRef] [Green Version]
  29. Vitart, F. Evolution of ECMWF sub-seasonal forecast skill scores. Q. J. R. Meteorol. Soc. 2014, 140, 1889–1899. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Lowe, R.; García-Díez, M.; Ballester, J.; Creswick, J.; Robine, J.-M.; Herrmann, F.R.; Rodó, X. Evaluation of an Early-Warning System for Heat Wave-Related Mortality in Europe: Implications for Sub-seasonal to Seasonal Forecasting and Climate Services. Int. J. Environ. Res. Public Health 2016, 13, 206. https://doi.org/10.3390/ijerph13020206

AMA Style

Lowe R, García-Díez M, Ballester J, Creswick J, Robine J-M, Herrmann FR, Rodó X. Evaluation of an Early-Warning System for Heat Wave-Related Mortality in Europe: Implications for Sub-seasonal to Seasonal Forecasting and Climate Services. International Journal of Environmental Research and Public Health. 2016; 13(2):206. https://doi.org/10.3390/ijerph13020206

Chicago/Turabian Style

Lowe, Rachel, Markel García-Díez, Joan Ballester, James Creswick, Jean-Marie Robine, François R. Herrmann, and Xavier Rodó. 2016. "Evaluation of an Early-Warning System for Heat Wave-Related Mortality in Europe: Implications for Sub-seasonal to Seasonal Forecasting and Climate Services" International Journal of Environmental Research and Public Health 13, no. 2: 206. https://doi.org/10.3390/ijerph13020206

APA Style

Lowe, R., García-Díez, M., Ballester, J., Creswick, J., Robine, J. -M., Herrmann, F. R., & Rodó, X. (2016). Evaluation of an Early-Warning System for Heat Wave-Related Mortality in Europe: Implications for Sub-seasonal to Seasonal Forecasting and Climate Services. International Journal of Environmental Research and Public Health, 13(2), 206. https://doi.org/10.3390/ijerph13020206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop