Next Article in Journal
Stress-Induced Endocrine and Immune Dysfunctions in Caregivers of People with Eating Disorders
Previous Article in Journal
Teacher Rated School Ethos and Student Reported Bullying—A Multilevel Study of Upper Secondary Schools in Stockholm, Sweden
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 14
Issue 12
10.3390/ijerph14121562
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impacts of the 2015 Heat Waves on Mortality in the Czech Republic—A Comparison with Previous Heat Waves

by
Aleš Urban
1,2,*,
Hana Hanzlíková
1,3,
Jan Kyselý
1,4,5 and
Eva Plavcová
1
1
Institute of Atmospheric Physics, Czech Academy of Sciences, Boční II 1401, 4 14131 Prague, Czech Republic
2
Faculty of Science, Charles University, Albertov 6, 2 12843 Prague, Czech Republic
3
Institute of Geophysics, Czech Academy of Sciences, Boční II 1401, 4 14131 Prague, Czech Republic
4
Faculty of Environmental Sciences, Czech University of Life Sciences, Kamýcká 129, 6 16521 Prague, Czech Republic
5
Global Change Research Centre, Czech Academy of Sciences, Bělidla 986, 60300 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2017, 14(12), 1562; https://doi.org/10.3390/ijerph14121562
Submission received: 29 September 2017 / Revised: 10 November 2017 / Accepted: 8 December 2017 / Published: 13 December 2017
(This article belongs to the Section Environmental Health)
Browse Figures
Versions Notes

Abstract

:
This study aimed to assess the impacts of heat waves during the summer of 2015 on mortality in the Czech Republic and to compare them with those of heat waves back to the previous record-breaking summer of 1994. We analyzed daily natural-cause mortality across the country’s entire population. A mortality baseline was determined using generalized additive models adjusted for long-term trends, seasonal and weekly cycles, and identified heat waves. Mortality deviations from the baseline were calculated to quantify excess mortality during heat waves, defined as periods of at least three consecutive days with mean daily temperature higher than the 95th percentile of annual distribution. The summer of 2015 was record-breaking in the total duration of heat waves as well as their total heat load. Consequently, the impact of the major heat wave in 2015 on the increase in excess mortality relative to the baseline was greater than during the previous record-breaking heat wave in 1994 (265% vs. 240%). Excess mortality was comparable among the younger age group (0–64 years) and the elderly (65+ years) in the 1994 major heat wave while it was significantly larger among the elderly in 2015. The results suggest that the total heat load of a heat wave needs to be considered when assessing its impact on mortality, as the cumulative excess heat factor explains the magnitude of excess mortality during a heat wave better than other characteristics such as duration or average daily mean temperature during the heat wave. Comparison of the mortality impacts of the 2015 and 1994 major heat waves suggests that the recently reported decline in overall heat-related mortality in Central Europe has abated and simple extrapolation of the trend would lead to biased conclusions even for the near future. Further research is needed toward understanding the additional mitigation measures required to prevent heat-related mortality in the Czech Republic and elsewhere.

1. Introduction

Many studies have documented seasonal and meteorological effects on human mortality in various geographical areas [1,2,3]. In mid-latitudes, mortality shows a typical pattern with a peak occurring during the cold season and the lowest mortality rates in summer [4,5]. If mortality anomalies relative to the seasonal pattern are considered, however, periods of extremely high temperatures (heat waves) have a significant impact on mortality such as those in 1995 in the U.S. [6], 2003 in southwestern Europe [7], or 2010 in Russia [8]. As summer mean temperatures are expected to increase during the 21st century [9], the frequency and intensity of extreme heat waves are likely to be significantly greater in comparison to recent climate [10,11,12].
As a result of climate change, increased heat-related mortality is projected in the future [13,14,15]. Meanwhile, its compensation by reduced cold-related mortality has been questioned [9,16,17]. Despite increasing temperature in recent decades, a significantly reduced risk of heat-related mortality has been observed [18,19]. This has been documented in North America [20,21,22,23,24], Europe [23,25,26,27,28], as well as Asia and Australia [23,29,30]. The magnitude and significance of heat-related mortality change may differ, however, according to region, length of the examined period, number and timing of heat waves, and the methodology used [23,26,31,32,33]. Although many studies have assessed temporal changes in relative risk of increased mortality per unit temperature rise [21,24], much less has been done to compare effects of heat events on mortality with respect to their meteorological and spatial-temporal characteristics [34,35,36].
Physiological (acclimatization), behavioural (clothing, working time), and technological (urban planning, building design, warning systems) adaptations have been often hypothesized as the reasons for declining vulnerability to extreme temperatures [21,24,25,26]. Although increased use of air conditioning has been considered one of the main causes of reduced susceptibility to heat [19,24,37], larger temporal declines in heat-related mortality were not significantly linked to larger increases in air conditioning prevalence when different cities and regions were compared [21,30]. Limited understanding as to which specific mechanisms and adaptation measures have caused the decrease in heat-related mortality [18,19,38] brings uncertainty to projections for the effect of future climate change on temperature-related mortality [13,14,15,39,40]. Therefore, further research is required to improve the effectiveness of such planned adaptation measures as heat warning systems [38,41]. This includes seeking to better understand heat-wave characteristics that affect human health.
In 1994, Central Europe experienced its hottest summer since at least the 1950s [42]. An exceptional impact on excess mortality was accordingly documented in several countries, including the Czech Republic [43,44], Germany [45], Poland [46], and The Netherlands [47]. Since 1994, several notable heat waves have occurred in the Central European area, including those in 2003 [44], 2006 [48,49], and 2013 [50]. Globally the hottest year to that date was 2015 [51], and Central European countries experienced intense, record-breaking heat waves [52,53]. In Switzerland, the summer of 2015 had as great impact on excess mortality as did the climatologically comparable summer of 2003 [54]. In England, on the other hand, the impact of the 2013 heat wave was considerably less than seen in 2003 and 2006 even though the temperature magnitude of the heat waves was comparable [55]. These two examples demonstrate that the impact of a given heat wave on mortality is affected by more factors than simple temperature characteristics. Duration and intensity of a heat stress period [25,36], timing of a heat wave within a season [44,56,57], air quality in a given area [49], and/or the magnitude of mortality in the previous winter due to influenza epidemics [44,58] can substantially influence the overall impact. Additionally, the structure of a given population, including its demographics and socioeconomic status, can significantly affect the impact of heat on mortality in a specific region [59,60,61,62]. Although the impact of hot weather on mortality is generally largest in urban environments [63], it may be significant also in rural areas [62,64].
Several studies have evaluated mean effects of hot spells on mortality [65] and morbidity [66] in the Czech Republic, as well as their spatial differences [67], using data up to 2009. Hůnová et al. [49] examined the effect of heat waves in 2003 and 2006 on mortality in Prague, but to date the impacts of the heat waves since 2010 have not been evaluated. This study aims to assess the impact of the 2015 heat waves on excess natural-cause mortality in the Czech Republic and to compare it with previous heat waves back to 1994 while focusing on the previously record-breaking summer of 1994. Meteorological characteristics of heat waves, considering their intensity, duration, and temperature change with respect to the previous conditions, are compared and linked to the magnitude of heat-related mortality during heat waves. The decreasing trend in heat-related mortality documented by previous studies is evaluated in the context of changes in impacts of heat waves, and demographic and meteorological factors modifying the observed impacts are discussed.

2. Data and Methods

2.1. Mortality Data

Daily natural-cause mortality (A00-R99 according to the International Classification of Diseases, version 10 (ICD-10)) was analysed in the Czech Republic during 1994–2015 (population about 10.5 million inhabitants as of 2015, having changed only little since 1994). The data were collected and processed by the Czech Statistical Office and the Institute of Health Information and Statistics of the Czech Republic. The data were stratified by gender and age into four subgroups—a total of 2,253,788 deaths were recorded in the national registry during the examined period, out of which 1,110,755 were males, 1,143,033 females, 486,789 younger than 65 years (0–64 years), and 1,766,999 deceased were 65 years and older (65+ years). As the demographic structure of the population has shifted toward a higher proportion of the elderly during the examined period (see Supplementary Materials Figure S1), we calculated daily mortality rates (per 10,000,000 inhabitants) in order to account for these changes. We applied the direct standardization procedure [67] using the mid-year population in each year and the standard WHO European population as a reference [68].

2.2. Meteorological Data

Daily air temperature was calculated as the mean from 12 high-quality meteorological stations operated by the Czech Hydrometeorological Institute, representing the spatial population distribution in the Czech Republic (Figure 1), and used as a proxy variable of thermal conditions. While some studies have shown a theoretical advantage of indices combining multiple variables in order to assess human thermal comfort [69] as well as their better performance on mortality data [70,71], the differences between various temperature measures in evaluating heat’s impact on mortality are generally insignificant [36,72,73,74].

2.3. Definition of a Heat Wave

To select only those heat waves affecting the majority of the area, we used one temperature time series calculated from the 12 stations distributed across the whole area. Heat waves are defined as periods of at least 3 days with daily mean temperature (Tmean) reaching or exceeding the 95th percentile (21.2 °C) of the annual distribution (estimated from data over 1994–2015), with at least 1 day exceeding the 98th percentile (23.2 °C). Although there is no universally accepted heat wave definition in public health studies, the 95th percentile of temperature distribution is generally understood as a threshold of a heat wave according to the WHO and WMO Guidance on Warning-System Development [75]. The second threshold (the 98th percentile) is employed in order to ensure that a heat wave involves at least one day with very high temperature; definitions of heat waves with two thresholds are relatively frequent in biometeorological and climatological studies [44,47,76].

2.4. Excess Heat Evaluation

Koppe and Jendritzky [56] pointed out the importance of short-term (within-season) acclimatization to heat for the overall impact of a heat wave on human mortality, and accordingly they introduced the HeRATE approach to heat-related mortality assessment. A similar method by Nairn and Fawcett [77] employs the excess heat factor (EHF) for evaluating the heat wave intensity with respect to previous thermal conditions. EHF consists of two excess heat indices (EHIs): the significance index (EHIsig), a measure of heat wave intensity, and the acclimatization index (EHIaccl), taking into account short-term acclimatization to ambient temperature over the recent past. Daily EHIsig and EHIaccl were calculated using Equations (1) and (2), respectively:
EHIsigi = (Ti + Ti–1 + Ti–2)/3 − T95, and
EHIaccli = (Ti + Ti–1 + Ti–2)/3 − (Ti–3 + … + Ti–32)/30,
where Ti denotes daily mean temperature on day i, T95 represents the 95th percentile of daily mean temperature, EHIsigi is defined as the exceedance of the previous three-day mean temperature (starting on day i) above T95, and EHIaccli on day i is calculated as the difference between the three-day mean temperature and a mean of the prior 30 days. It should be noted that we slightly modified the original Nairn and Fawcett’s (2015) equations by calculating the three-day mean temperature over days i to i − 2 instead of days i to i + 2, since human mortality on day i can be affected by thermal conditions on day i or previous days, but not on upcoming days. Consequently, daily EHF was calculated using Equation (3):
EHFi = max(0,EHIsigi) × max(1,EHIaccli)
Finally, values of EHIsig, EHIaccl, and EHF on all days of a heat wave were summed, in order to evaluate its total heat load.

2.5. Statistical Analyses

2.5.1. Excess Mortality Calculation

To assess the mortality increase during heat waves, daily baseline for natural-cause mortality was determined using generalized additive models (“mgcv” package in R (version 2.15.2) [78]). Specifically, we employed Equation (4):
Log[E(M)] = α + s(time; df = 6 × n(y)) + factor(dow) + factor(hwd),
where E(M) is the predicted daily mortality, α is the model intercept, time is a variable to account for long-term trends and seasonality, n(y) is the number of years in the time series, s() represents thin plate smoothing spline [78], and factor variables define categorical variables for days of week (dow) and heat-wave days (hwd). A smoothing spline with six degrees of freedom (df) per year (total df = 132) was applied after sensitivity tests evaluated by a combination of procedures, including partial autocorrelation plots, models’ unbiased risk estimator criterion [78], percentage of deviation explained, and visual examination of model-predicted baseline plots. The final number of df used was also in accordance with an estimate (137 df) made via generalized cross-validation using the “bruto” function in the “mda” R package (version 0.4-9 [79]) according to Peng and Dominici [80]. The weekly cycle in mortality data, with maximum on Monday and minimum on Sunday was adjusted by the categorical variable dow. Finally, the mortality baseline was adjusted by a categorical variable representing days of heat waves identified in this study (hwd). The hwd variable reaches values 0/1, where 1 is used for the heat wave days. The coefficient of the hwd variable was subtracted from the model-predicted baseline on heat-wave days and therefore the excess mortality during the heat waves was calculated with respect to the baseline adjusted for the effect of heat waves [55].
Daily numbers of deaths per standardized 10 million inhabitants and their relative deviations from the mortality baseline (excess mortality) are shown in Figure S2. Cumulative excess mortality per standardized 10 million inhabitants during the identified heat waves was calculated for the population as a whole and individually for population groups showing usually different response to heat: males and females, and age groups younger and older than 65 years (0–64 and 65+ years, respectively). The significance of the mortality excess was evaluated by 95% confidence intervals of the cumulative excess mortality, calculated as the lower and upper limit factors for a Poisson-distributed variable according to Schoenberg [81].

2.5.2. Associations between Heat Waves and Excess Mortality

Characteristics of heat waves such as average daily mean temperature (average Tmean), intensity (EHIsig), duration, and timing within a year have been documented as important factors modifying the eventual impact on mortality [34,36]. Using the cumulative excess heat factor (∑EHF) allows us to consider all these important meteorological components defining the total heat load of a heat wave together—its intensity represented by the temperature anomaly during the heat wave (EHIsig), timing within a season represented by the temperature change with respect to the previous conditions (EHIaccl), and duration (∑). To evaluate associations between heat waves and mortality with respect to the baseline mortality in a given year, the sum of relative mortality deviations (∑RMD) during a heat wave was calculated and linked to the above mentioned heat wave characteristics, using linear regression models.

2.5.3. Mortality Displacement

Although heat waves are defined by the temperature thresholds presented above, the actual periods of heat stress are likely to be slightly longer (e.g., a hot period may have started to develop several days before a heat wave) and the total effect of a heat wave may be slightly greater than that of the identified heat wave period. On the other hand, the effect of a short-term mortality displacement (harvesting) after the end of a heat wave may substantially reduce the net effect of a heat wave on excess mortality [82,83]. To evaluate the magnitude of these two effects, we defined a so-called “extended heat wave period” (EHP) [43] with several days of positive mortality deviations before a heat wave and several days with negative mortality deviations after a heat wave. A three-day centered moving average of mortality deviations was calculated to evaluate the mortality displacement. The beginning of EHP was delimited as the first three-day period with a positive average mortality deviation before a heat wave while the end of EHP was defined as the last three-day period with a negative average mortality deviation after a heat wave [43]. A positive phase of EHP lasted till the central day of the last three-day period with positive averaged mortality deviations and the rest of EHP was defined as a negative phase. Multiple heat events separated by days without a drop of the three-day average of mortality deviations below zero were considered as a single EHP. The final magnitude of the mortality displacement was calculated as the absolute value of the ratio of the sum of the three-day-averaged mortality deviations during the negative phase of EHP to the sum of the three-day-averaged mortality deviations during the positive phase.

3. Results

3.1. General Characteristics of Heat Waves

A total of 49 heat waves were identified between 1994 and 2015 in the Czech Republic, with maximum of four heat waves per year (Table 1). Although we did not strictly specify the season of our analyses, all identified heat waves (except one in May 2005) occurred in summer (June–August). Most of the heat waves were associated with positive cumulative mortality deviation and the positive anomalies were significantly (p = 0.05) different from the mortality baseline in 40 events (see Supplementary Materials Tables S1 and S2A,B). The mean relative rise in natural-cause mortality during the heat waves was 10.7% per day, while the maximum daily excess mortality was greater than 30% on the hottest days. Differences in the mean impact of heat waves were observed between population groups. While the mean relative mortality deviation during heat waves was not significantly (t-test at p = 0.05) different between women and men (11.7% vs. 10.1%, respectively), it was significantly larger in the elderly (65+ years) than the younger (0–64 years) age group (11.6% vs. 8.3%, respectively).
2015 was the record-breaking year during the examined period regarding the total duration of heat waves (36 days) as well as their cumulative excess heat factor (∑EHF = 606 °C; Figure 2). Table 1 shows meteorological characteristics of heat waves in six years with the total heat wave duration longer than 15 days. Although the longest heat wave during the examined period occurred in 1994 (18 days; 22 July–8 August; denoted 1994_4), the heat wave 2015_3, lasting 13 days (3–15 August), was the hottest in the study period according to the average daily mean temperature (avgTmean = 26.0 °C, compared to 24.1 °C in 1994_4) as well as ∑EHF (292.4 °C and 176.6 °C, respectively).

3.2. Effects of the Major Heat Waves on Excess Mortality

The two persistent heat waves—1994_4 and 2015_3 (hereafter termed major heat waves)—were remarkable not only because of their meteorological characteristics, but also due to their impact on excess mortality, which was substantially larger than during any other heat wave (Table 2 and Table S2A,B). The 1994_4 major heat wave was associated with 1197 excess deaths (95% confidence interval (95% CI: 1002–1397) per standardized 10 million inhabitants, in comparison to 847 (95% CI: 711–988) deaths per 10 million inhabitants during 2015_3. However, these values represented 240% and 265% cumulative relative mortality deviation (∑RMD), respectively, thus indicating that the impact on mortality relative to the mortality baseline in the given year was larger in 2015. A substantial difference in the impact of the two major heat waves was observed between the two age groups. While ∑RMD was comparable among the younger and the elderly age groups (254% and 234%, respectively) in the 1994_4 heat wave, it was significantly larger (p = 0.05) among the elderly in 2015_3 (132% and 313%, respectively).
Although the fourth heat wave in 2006 (18–28 July; 2006_4) was the third longest heat wave during the examined period (11 days) and avgTmean was equal to 1994_4 (24.1 °C), its impact on ∑RMD (93%) ranked as the 11th largest among all heat waves. Similarly, the longest heat wave in 2003 had a relatively small impact on mortality (∑RMD = 60%), despite its 10-day duration (1–10 August; 2003_3). In both cases, the values of ∑EHF (72.7 °C in 2003_3 and 99.7 °C in 2006_4, respectively) were much lower than during the major heat waves (1994_4 and 2015_3). On the contrary, other heat waves, such as 1994_1, 2007_1 (Tables S1 and S2B), 2013_1, and 2015_1 had quite a strong impact on excess mortality despite their relatively short duration (4, 8, 5, and 7 days, respectively). These heat waves occurred as the first heat waves in a given season, accompanied by a substantial increase of temperature with respect to previous conditions (Table 1 and Figure 3). Therefore, their total heat load was also relatively large (∑EHF = 99.4, 171.3, 158.7, and 143.1 °C, respectively).

3.3. Associations between Heat Wave Characteristics and Excess Mortality

The above presented results highlight the importance of considering the total heat load of a heat wave when evaluating its impact on human health. This finding was also supported by the linear regression analysis performed among all 49 heat waves. As shown in Figure 4, increasing ∑EHF during a heat wave was significantly linked to increasing ∑RMD (R2 = 0.75). The ∑RMD of a heat wave depended more strongly on ∑EHF than other indicators of heat wave intensity such as average (avgTmean, R2 = 0.31) and maximum daily mean temperature (maxTmean, R2 = 0.31) during a heat wave, as well as its duration (R2 = 0.59). Table 3 shows that relationships between ∑RMD and independent variables were generally stronger among the elderly (65+ years) than the younger (0–64 years) age group.

3.4. Mortality Displacement

Generally, the order of heat waves with the largest impact on mortality was maintained also if the mortality displacement effect during the extended heat wave period (EHP) was considered (Table 2). The net cumulative mortality deviation was positive for most EHPs (see Table S3). In several cases, the net cumulative mortality deviation during EHPs was larger than during the heat waves alone due to an extraordinarily long period of consecutive positive daily excess mortality preceding or following the heat wave. The proportion of the harvesting effect on excess mortality during the positive phase of EHPs was estimated to be on average 23% during the most extreme heat waves (with ∑EHF ≥ 100 °C), while it was larger, on average 54%, during heat waves with ∑EHF < 100 °C. The magnitude of the mortality displacement effect varied substantially between individual cases, however. While the magnitude of the mortality displacement exceeded 30% during EHPs of the intense but relatively short heat waves in 2007 (2007_1) and 2013 (2013_1) (38% and 37%, respectively), comparable heat waves in 2010 (2010_3) and 2015 (2015_1) revealed a substantially lower harvesting effect (5% and 10%, respectively). In the latter two cases, the harvesting effect of the heat waves was weakened by a subsequent hot period occurring shortly thereafter (Figure 3). This was true especially in 2015, when the 2015_1 heat wave (7 days; ∑EHF = 143.1 °C) was followed by the even more severe 2015_2 heat wave (10 days; ∑EHF = 145.1 °C), but the impact of the second heat wave on excess mortality was substantially weaker. When only the two major heat waves (1994_4 and 2015_3) are compared, the magnitude of the mortality displacement was slightly larger in 1994 (27%) than in 2015 (22%).

4. Discussion

We analyzed effects of heat waves on natural-cause mortality in the Czech Republic during summer 2015 and compared them with the effects of previous heat waves back to 1994. Two major heat waves, occurring on 22 July–8 August 1994 and 3–15 August 2015, respectively, were exceptional in both the climatological extremity and their impact on excess mortality. Although the excess mortality (deaths per standardized 10 million inhabitants) was higher in 1994, the increase in mortality relative to the baseline in a given year was larger in 2015. These findings suggest that the health response to the major heat-stress events has not improved significantly during the examined period and that the total heat load of a heat wave, given by its intensity, duration, and temperature change with respect to the previous conditions, constitutes one of the main factors affecting the magnitude of the mortality increase.

4.1. Effect of the Total Heat Load

Studies reporting significant decrease in heat-related mortality over time have mostly assessed temporal changes in relative risk of increased mortality per unit of temperature increase, aggregating all days above a certain threshold [21,24]. This approach may improve our knowledge of temperature–mortality relationships in general, but it does not tell much about effects of individual heat events with respect to their meteorological and spatial-temporal characteristics [34,35,36]. The extraordinary impact of summer 2015 on heat-related mortality has been already documented in several Central European countries [54,77]. The relative mortality increase during the 2015 major heat wave in the Czech Republic was larger than in 1994, even though the 2015 event was shorter by five days. When comparing the overall magnitude of the two major heat waves, however, the 2015 one was the most extreme during the examined period as measured by average daily mean temperature during the heat wave (avgTmean), the cumulative temperature anomaly above the 95th percentile threshold (∑EHIsig), cumulative temperature increase with respect to the previous conditions (∑EHIaccl), as well as cumulative excess heat factor (∑EHF). These findings are in accordance with those of Hoy et al. [53], who have reported summer 2015 in the Central European region as record-breaking since the beginning of the 20th century in terms of several heat-wave characteristics.
On the other hand, while the summer of 2003 (June–August) was the second warmest in the Czech Republic during the examined period (according to mean temperature), temperature increases during heat waves were generally small compared to the non-heat wave conditions, thus resulting in low cumulative excess heat factor of the heat waves in 2003. Therefore, despite its length comparable to the 2015 major heat wave, the total heat load of the August 2003 heat wave was substantially lower. The importance of the total heat load of a heat wave (considered as a function of its absolute temperature anomaly, temperature change, and duration [34,36]) in the final effect on heat-related mortality was demonstrated in this study by the strong correlation between the cumulative relative mortality deviation (∑RMD) and ∑EHF during a heat wave. The difference between the low total heat load during the 2003 heat wave in Central Europe demonstrated here as well as in other studies [31,42,45,77] and, on the other hand, its strong intensity in Western Europe [50,52], contributed significantly to the spatial differences in the mortality response across Europe [44].
Despite the declared importance of the total heat load in the final impact of a heat wave on excess mortality, heat waves with comparable ∑EHF in our study had in some cases substantially different impact on ∑RMD. These variations were related mainly to different conditions preceding the heat waves. For example, while in 2007 a mild summer without obvious temperature extremes prevailed in Central Europe prior to the heat wave with the third largest ∑RMD during the examined period, a comparably intense heat wave in 2010 (2010_3) had a weaker impact on excess mortality as it was preceded by two other heat waves earlier in the season. On the other hand, the cumulative excess mortality during both heat waves was comparable when the extended heat wave periods (EHPs) including the harvesting effect were considered. This finding emphasizes limitations of the regression analyses of associations between ∑RMD and heat wave characteristics, assuming that all heat-related excess deaths occurred only during the identified heat waves. On the other hand, it is also not very likely that all excess deaths that occurred during EHPs were associated to the heat stress. Although we attempted to estimate mortality displacement during individual EHPs, the resulting values have to be perceived very carefully, since quantifying the real proportion of deaths that would have occurred regardless of the heat wave is still not a completely resolved issue [82,83].
In addition, the study did not take into account other potential confounding variables or effect modifiers, such as air pollution [84] and mortality in previous winters due to influenza epidemics [85,86]. As air pollution concentrations vary considerably across space and causal relationships among temperature, air pollution, and human health are associated with large uncertainty [87], we did not consider the effect of air pollution in this study. Possible effects of the influenza epidemic in the winter season 2002/2003 on reduced excess mortality during the following summer have been discussed by Kyselý and Kříž [44]. A study from Slovakia [76] suggests that the impact of summer 2015 on excess mortality could have been even stronger, if no influenza epidemic occurred earlier in the year.

4.2. Changes in Population Structure and Socioeconomic Factors

In addition to meteorological factors, demographic changes over time may significantly modify the effect of heat. Age-standardization of mortality is commonly applied in epidemiological studies in order to adjust the time series for changes in population structure over time and space [68]. Due to age-standardization, lower heat wave-related mortality (deaths per 10 million inhabitants) was estimated in 2015 than in 1994. However, when the relative mortality deviations from the baseline were considered, the impact of the 2015 major heat wave was even greater than that of the 1994 heat wave. This finding seems to be in contradiction with Kyselý and Plavcová [26], who found a significantly decreasing trend in relative excess heat-related mortality in the Czech Republic over 1986–2009. The significant decrease was related primarily to substantial socioeconomic development, improved quality of life (including more efficient health care and lifestyle changes), and significantly decreased environmental deprivation due to rapid transition from a developing to developed economy in the Czech Republic (and other Central European countries) in the 1990s, following the so-called “Velvet Revolution” and end of the communist regime in 1989 [26]. Socioeconomic and lifestyle changes were much less pronounced in the most recent past, which may have contributed to abating the trend in heat-related mortality. A third factor contributing to significantly decreasing heat-related mortality over 1986–2009 and reported by Kyselý and Plavcová [26] might have been essentially the absence of a major heat wave with an impact on mortality comparable to that of the 1994 event in the second half of their examined period. Ha and Kim [32], for example, excluded summer 1994 with an unprecedented heat wave in South Korea from their analysis, thus resulting in predominantly insignificant changes in temperature-related mortality over time (1993–2009).
Significant decreases in heat-related mortality risk during recent decades have been observed also in the U.S. [20] and other developed regions of the world [23,25,28,30]. These have been associated mainly with improved health care, better public awareness of heat-related risks, and increased prevalence of air conditioning [20,24,88]. The decline has abated since the mid-1990s in the U.S., however, and a further substantial decrease may be unlikely, as the air conditioning availability is already near its saturation in some regions [24,88,89]. The results of our study suggest that a similar abatement of this favourable development may be observed or expected in other regions, too, albeit for different reasons (air conditioning is much less widespread in Central Europe), and that major or unprecedented heat waves may substantially influence the estimated trends. As their frequency, intensity, and duration are projected to increase in a warmer climate, a simple extrapolation of observed trends in heat-related mortality would lead to biased conclusions even for the near future.
Together with the favourable socioeconomic development over recent decades, the demographic structure of the population has shifted toward a higher proportion of the elderly during the same period. While in 1994, 13% of the Czech population was older than 65 years of age, this figure was already 18% in 2015 [90]. Our results showed that the relative increase in mortality was comparable between the younger age group (0–64 years) and the elderly (65+ years) during the 1994 major heat wave, while it was substantially larger among the elderly in 2015. These findings suggest that the 0–64 years age group has become less vulnerable to heat stress over time, while the total impact of the 2015 major heat wave was comparable with that of the 1994 heat wave due to increased proportion of the elderly (65+ years) among the affected population. An analysis at a finer detail of age groups would be needed to understand our findings in the context of changing population structure and increasing life expectancy. However, as the elderly with chronic cardiovascular and respiratory diseases constitute the population group at the highest risk of heat-related death [57,65,66] and their proportion in the total population has been increasing over time, they should be in the particular interest when developing new adaptation measures reducing impacts of the extreme heat on public health. In order to quantify the effects of population ageing on changes in the heat-vulnerability properly, a measure such as years of life lost [16]—considering the life expectancy at the time of death—should be analyzed in follow-up research.

4.3. Heat-and-Health Warning Systems

Heat-and-health warning systems (HHWSs) have been often presented as an effective tool of heat-related mortality alleviation [9,41,91]. Reduced impact of heat waves on mortality has been documented in southwestern European countries after the 2003 heat wave, when local HHWSs were set up [55,92,93]. On the other hand, a study by de’Donato et al. [33] comparing changes in heat-related mortality risk between two seven-year periods before and after 2003 showed that the beneficial effect of HHWSs has not been consistent across Europe. In spite of many methodological similarities among HHWSs in various countries, significant differences in key characteristics of the plans exist [94], especially regarding steps following the heat alert issuance. Therefore, it is necessary to evaluate the effectiveness of these systems in order to understand which actions actually increase resilience of populations [9,38,41,94].
The Czech Republic does not so far have an HHWS based on temperature–health relationships and including steps to follow after an alert’s issuance [94]. As demonstrated by recent studies [33,95], North and Central European cities may be particularly vulnerable to increasing intensity and frequency of heat waves in the near future. The findings of this study suggest that meteorological characteristics of a predicted heat period should be considered within the context of the previous temperature and mortality development when its possible impact on excess mortality is estimated [44,55]. Otherwise, the heat wave’s severity in relation to impacts may be misinterpreted. Excess heat factor [77] or a similar method considering the total heat load of a heat wave with respect to previous conditions may be one of the ways to improve HHWSs (as already documented by [96]), as better understanding of weather events is essential for properly predicting the possible human health impact [24,97].
For better understanding of associations between weather events and human health when developing an effective HHWS, the WMO/WHO guideline [75] recommends consideration of location-specific relationships between ambient temperature and mortality. Therefore, it needs to be acknowledged that associations between heat waves and excess mortality found in this study for the nation-wide population may not necessarily represent associations existing in individual cities and regions. As previous study has shown [67], however, modelling temperature-mortality relationships in relatively small cities and regions is difficult due to small sample sizes that do not bring robust results. Moreover, as the aim of this study was to compare the impact of the extraordinary summer 2015 on excess mortality with the previous important heat waves back to 1994 and not to define a HHWS methodology for a specific city, we followed the methodological approach applied in previous studies assessing the impact of heat waves on mortality in the Czech Republic and its changes over time [26,43,44].

5. Conclusions

Impacts of heat waves in summer of 2015 on natural-cause mortality in the Czech Republic were evaluated and compared with those of previous heat waves back to the record-breaking summer of 1994. The summer 2015 exceeded the 1994 one in both the intensity of heat waves and magnitude of heat-related mortality. Despite the smaller excess mortality associated with the major heat wave in 2015 in comparison to 1994, the relative increase from the baseline was larger in 2015. The effect of the major 2015 heat wave was more pronounced in particular due to stronger impact on excess mortality among the elderly (65+ years) than the younger age group (0–64 years) in 2015 compared to 1994. The results suggest that the cumulative excess heat factor of a heat wave, considered as a function of absolute temperature anomaly during a heat wave, heat wave duration, and temperature increase relative to the previous conditions, explains better the magnitude of excess mortality during a heat wave than other characteristics such as duration and average daily mean temperature during the heat wave alone. Further research is needed to improve the effectiveness of measures to prevent heat-related mortality while considering changes in the intensity, duration, and frequency of heat waves due to ongoing climate change as well as changes in the demographic, epidemiological, and socioeconomic structure of the population over space and time. The comparison of mortality impacts of the 2015 and 1994 major heat waves suggests that the recently reported decline in overall heat-related mortality in Central Europe has abated and a simple extrapolation of the trend would lead to biased conclusions even for the near future.

Supplementary Materials

The following are available online at www.mdpi.com/1660-4601/14/12/1562/s1, Table S1: Characteristics of heat waves during 1994–2015 in the Czech Republic, Table S2A: Impact of heat waves on mortality due to natural causes for the whole population, males and females, and the younger (0–64 years) and the elderly (65+ years) population groups, during 1994–2004 in the Czech Republic, Table S2B: Impact of heat waves on mortality due to natural causes for the whole population, males and females, and the younger (0–64 years) and the elderly (65+ years) population groups, during 2005–2015 in the Czech Republic, Table S3: Net mortality change over extended heat wave periods during 1994–2015 in the Czech Republic, Figure S1: Population structure in the Czech Republic in 1994 and 2015, Figure S2: Top: Daily numbers of natural-cause deaths per standardized 10 million inhabitants in the Czech Republic during 1994–2015. The black line denotes the mortality baseline adjusted for long-term trend, seasonality, weekly cycle, and the effect of heat waves. Bottom: Relative mortality deviations (excess mortality) from the baseline in the same population and period.

Acknowledgments

Aleš Urban was supported by the Czech Academy of Sciences Programme for Research and Mobility Support of Starting Researchers, project NO. MSM100421604, and by the Charles University Grant Agency, student project NO. 3300-243-227 249. Jan Kyselý and Eva Plavcová were supported by the Czech Science Foundation, project NO. 16-22000S. Thanks are due to Bohumír Kříž and Jan Kynčl, National Institute of Public Health, for preparing epidemiological datasets as well as for useful discussions concerning their interpretation. This study was partly elaborated during Aleš Urban’s research stay at the School of Geographical Science & Urban Planning, Arizona State University. We acknowledge the excellent working conditions and support provided by the institute. We are particularly grateful to David M. Hondula, Arizona State University, for reviewing this manuscript and providing helpful comments. Thanks go also to all anonymous reviewers that helped improve the original manuscript.

Author Contributions

Aleš Urban, Hana Hanzlíková and Jan Kyselý defined the study design. Hana Hanzlíková was responsible for gathering, cleaning and standardization of the mortality data. Aleš Urban was responsible for gathering and cleaning of the meteorological data. Aleš Urban and Eva Plavcová carried out the statistical analyses and prepared the figures. Aleš Urban drafted the first manuscript. Hana Hanzlíková, Jan Kyselý and Eva Plavcová provided critical revisions of the manuscript. All authors contributed to discussion on the study design and the final interpretation of results, and read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gosling, S.N.; Lowe, J.A.; McGregor, G.R.; Pelling, M.; Malamud, B.D. Associations between elevated atmospheric temperature and human mortality: A critical review of the literature. Clim. Chang. 2009, 92, 299–341. [Google Scholar] [CrossRef]
  2. Hajat, S.; Kosatky, T. Heat-related mortality: A review and exploration of heterogeneity. J. Epidemiol. Community Health 2010, 64, 753–760. [Google Scholar] [CrossRef] [PubMed]
  3. Burkart, K.; Khan, M.M.H.; Schneider, A.; Breitner, S.; Langner, M.; Krämer, A.; Endlicher, W. The effects of season and meteorology on human mortality in tropical climates: A systematic review. Trans. R. Soc. Trop. Med. Hyg. 2014, 108, 393–401. [Google Scholar] [CrossRef] [PubMed]
  4. Bhaskaran, K.; Gasparrini, A.; Hajat, S.; Smeeth, L.; Armstrong, B. Time series regression studies in environmental epidemiology. Int. J. Epidemiol. 2013, 42, 1187–1195. [Google Scholar] [CrossRef] [PubMed]
  5. Huang, C.; Chu, C.; Wang, X.; Barnett, A.G. Unusually cold and dry winters increase mortality in Australia. Environ. Res. 2015, 136, 1–7. [Google Scholar] [CrossRef] [PubMed]
  6. Semenza, J.C.; Mccullough, J.E.; Flanders, W.D.; Mcgeehin, M.A.; Lumpkin, J.R. Excess hospital admissions during the July 1995 Heat Wave in Chicago. Am. J. Prev. Med. 1999, 16, 269–277. [Google Scholar] [CrossRef]
  7. Larsen, J. Setting the Record Straight: More than 52,000 Europeans Died from Heat in Summer 2003. Earth Policy Institute, 2006. Available online: http://www.earth-policy.org/plan_b_updates/2006/update56 (accessed on 20 March 2017).
  8. Barriopedro, D.; Fischer, E.M.; Luterbacher, J.; Trigo, R.M.; García-Herrera, R. The hot summer of 2010: Redrawing the temperature record map of Europe. Science 2011, 332, 220–224. [Google Scholar] [CrossRef] [PubMed]
  9. Intergovernmental Panel on Climate Change (IPCC). Climate Change 2014: Impacts, Adaptation, and Vulnerability; Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  10. Ballester, J.; Rodó, X.; Giorgi, F. Future changes in Central Europe heat waves expected to mostly follow summer mean warming. Clim. Dyn. 2009, 35, 1191–1205. [Google Scholar] [CrossRef]
  11. Fischer, E.M.; Schär, C. Consistent geographical patterns of changes in high-impact European heatwaves. Nat. Geosci. 2010, 3, 398–403. [Google Scholar] [CrossRef]
  12. Lhotka, O.; Kyselý, J.; Farda, A. Climate change scenarios of heat waves in Central Europe and their uncertainties. Theor. Appl. Climatol. 2017. [Google Scholar] [CrossRef]
  13. Huang, C.; Barnett, A.G.; Wang, X.; Vaneckova, P.; Fitzgerald, G.; Tong, S. Projecting Future heat-related mortality under climate change scenarios: A systematic review. Environ. Health Perspect. 2011, 119, 1681–1690. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Li, T.; Ban, J.; Horton, R.M.; Bader, D.A.; Huang, G.; Sun, Q.; Kinney, P.L. Heat-related mortality projections for cardiovascular and respiratory disease under the changing climate in Beijing, China. Sci. Rep. 2015, 5, 4–11. [Google Scholar] [CrossRef] [PubMed]
  15. Kim, D.-W.; Ravinesh, D.K.; Chung, C.D.J. Projection of heat wave mortality related to climate change in Korea. Nat. Hazards 2016, 80, 623–637. [Google Scholar] [CrossRef]
  16. Huang, C.R.; Barnett, A.G.; Wang, X.M.; Tong, S. Effects of extreme temperatures on years of life lost for cardiovascular deaths: A time series study in Brisbane, Australia. Circ. Cardiovasc. Qual. Outcomes 2012, 5, 609–614. [Google Scholar] [CrossRef] [PubMed]
  17. Li, T.; Horton, R.M.; Kinney, P.L. Projections of seasonal patterns in temperature-related deaths for Manhattan, New York. Nat. Clim. Chang. 2013, 3, 717–721. [Google Scholar] [CrossRef] [PubMed]
  18. Hondula, D.M.; Balling, R.C.; Vanos, J.K.; Georgescu, M. Rising temperatures, human health, and the role of adaptation. Curr. Clim. Chang. Rep. 2015, 1, 144–154. [Google Scholar] [CrossRef]
  19. Arbuthnott, K.; Hajat, S.; Heaviside, C.; Vardoulakis, S. Changes in population susceptibility to heat and cold over time: Assessing adaptation to climate change. Environ. Health 2016, 15, S33. [Google Scholar] [CrossRef] [PubMed]
  20. Davis, R.E.; Knappenberger, P.C.; Novicoff, W.M.; Michaels, P.J. Decadal changes in summer mortality in U.S. cities. Int. J. Biometeorol. 2003, 47, 166–175. [Google Scholar] [CrossRef] [PubMed]
  21. 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, 811–816. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Petkova, E.P.; Gasparrini, A.; Kinney, P.L. Heat and mortality in New York City since the beginning of the 20th century. Epidemiology 2014, 25, 554–560. [Google Scholar] [CrossRef] [PubMed]
  23. Gasparrini, A.; Guo, Y.; Hashizume, M.; Kinney, P.L.; Petkova, E.P.; Lavigne, E. Temporal variation in heat-mortality associations: A multicountry study. Environ. Health Perspect. 2015, 123, 1200–1207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Sheridan, S.C.; Dixon, P.G. Spatiotemporal trends in human vulnerability and adaptation to heat across the United States. Anthropocene 2016. [Google Scholar] [CrossRef]
  25. Matzarakis, A.; Muthers, S.; Koch, E. Human biometeorological evaluation of heat-related mortality in Vienna. Theor. Appl. Climatol. 2011, 105, 1–10. [Google Scholar] [CrossRef]
  26. Kyselý, J.; Plavcová, E. Declining impacts of hot spells on mortality in the Czech Republic, 1986–2009: Adaptation to climate change? Clim. Chang. 2012, 113, 437–453. [Google Scholar] [CrossRef]
  27. Schifano, P.; Leone, M.; De Sario, M.; de’Donato, F.; Bargagli, A.M.; D’Ippoliti, D. Changes in the effects of heat on mortality among the elderly from 1998–2010: Results from a multicenter time series study in Italy. Environ. Health 2012, 11, 58. [Google Scholar] [CrossRef] [PubMed]
  28. Åström, D.O.; Forsberg, B.; Edvinsson, S.; Rocklöv, J. Acute fatal effects of short-lasting extreme temperatures in Stockholm, Sweden: Evidence across a century of change. Epidemiology 2013, 24, 820–829. [Google Scholar] [CrossRef] [PubMed]
  29. Coates, L.; Haynes, K.; O’Brien, J.; McAneney, J.; de Oliveira, F.D. Exploring 167 years of vulnerability: An examination of extreme heat events in Australia 1844–2010. Environ. Sci. Policy 2014, 42, 33–44. [Google Scholar] [CrossRef]
  30. Ng, S.; Boeckmann, M.; Ueda, K.; Zeeb, H.; Nitta, H.; Watanabe, C.; Honda, Y. Heat-related mortality: Effect modification and adaptation in Japan from 1972 to 2010. Glob. Environ. Chang. 2016, 39, 234–243. [Google Scholar] [CrossRef]
  31. Výberči, D.; Švec, M.; Faško, P.; Savinová, H.; Trizna, M.; Mičietová, E. The effects of the 1996–2012 summer heat events on human mortality in Slovakia. Morav. Geogr. Rep. 2015, 23, 57–70. [Google Scholar] [CrossRef]
  32. Ha, J.; Kim, H. Changes in the association between summer temperature and mortality in Seoul, South Korea. Int. J. Biometeorol. 2013, 57, 535–544. [Google Scholar] [CrossRef] [PubMed]
  33. De’Donato, F.K.; Leone, M.; Scortichini, M.; De Sario, M.; Katsouyanni, K.; Lanki, T.; Basagaña, X.; Ballester, F.; Åström, C.; Paldy, A.; et al. Changes in the effect of heat on mortality in the last 20 years in nine European cities. Results from the PHASE project. Int. J. Environ. Res. Public Health 2015, 12, 15567–15583. [Google Scholar] [CrossRef] [PubMed]
  34. Anderson, G.B.; Bell, M.L. Heat waves in the United States: Mortality risk during heat waves and effect modification by heat wave characteristics in 43 U.S. communities. Environ. Health Perspect. 2011, 119, 210–218. [Google Scholar] [CrossRef] [PubMed]
  35. Gasparrini, A.; Armstrong, B. The impact of heat waves on mortality. Epidemiology 2011, 22, 68–73. [Google Scholar] [CrossRef] [PubMed]
  36. Barnett, A.G.; Hajat, S.; Gasparrini, A.; Rocklöv, J. Cold and heat waves in the United States. Environ. Res. 2012, 112, 218–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ostro, B.; Rauch, S.; Green, R.; Malig, B.; Basu, R. The effects of temperature and use of air conditioning on hospitalizations. Am. J. Epidemiol. 2010, 172, 1053–1061. [Google Scholar] [CrossRef] [PubMed]
  38. Boeckmann, M.; Rohn, I. Is planned adaptation to heat reducing heat-related mortality and illness? A systematic review. BMC Public Health 2014, 14, 1112. [Google Scholar] [CrossRef] [PubMed]
  39. Hondula, D.M.; Georgescu, M.; Balling, R.C. Challenges associated with projecting urbanization-induced heat-related mortality. Sci. Total Environ. 2014, 490, 538–544. [Google Scholar] [CrossRef] [PubMed]
  40. Petkova, E.P.; Vink, J.K.; Horton, R.M.; Gasparrini, A.; Bader, D.A.; Francis, J.D.; Kinney, P.L. Towards more comprehensive projections of urban heat-related mortality: Estimates for New York City under multiple population, adaptation, and climate scenarios, (January). Environ. Health Perspect. 2016. [Google Scholar] [CrossRef] [PubMed]
  41. Toloo, G.S.; Fitzgerald, G.; Aitken, P.; Verrall, K.; Tong, S. Are heat warning systems effective? Environ. Health 2013, 12, 27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Lhotka, O.; Kyselý, J. Characterizing joint effects of spatial extent, temperature magnitude and duration of heat waves and cold spells over Central Europe. Int. J. Climatol. 2015, 35, 1232–1244. [Google Scholar] [CrossRef]
  43. Kyselý, J. Mortality and displaced mortality during heat waves in the Czech Republic. Int. J. Biometeorol. 2004, 49, 91–97. [Google Scholar] [CrossRef] [PubMed]
  44. Kyselý, J.; Kříž, B. Decreased impacts of the 2003 heat waves on mortality in the Czech Republic: An improved response? Int. J. Biometeorol. 2008, 52, 733–745. [Google Scholar] [CrossRef] [PubMed]
  45. Gabriel, K.M.A.; Endlicher, W.R. Urban and rural mortality rates during heat waves in Berlin and Brandenburg, Germany. Environ. Pollut. 2011, 159, 2044–2050. [Google Scholar] [CrossRef] [PubMed]
  46. Kuchcik, M.; Degórski, M. Heat- and cold-related mortality in the north-east of Poland as an example of the socio-economic effects of extreme hydrometeorological events in the Polish Lowland. Geogr. Pol. 2009, 82, 69–78. [Google Scholar] [CrossRef]
  47. Huynen, M.M.; Martens, P.; Schram, D.; Weijenberg, M.P.; Kunst, A.E. The impact of heat waves and cold spells on mortality rates in the Dutch population. Environ. Health Perspect. 2001, 109, 463–470. [Google Scholar] [CrossRef] [PubMed]
  48. Kyselý, J. Recent severe heat waves in central Europe: How to view them in a long-term prospect? Int. J. Climatol. 2010, 30, 89–109. [Google Scholar] [CrossRef]
  49. Hůnová, I.; Brabec, M.; Malý, M.; Knobová, V.; Braniš, M. Major heat waves of 2003 and 2006 and health outcomes in Prague. Air Qual. Atmos. Health 2016. [Google Scholar] [CrossRef]
  50. Lhotka, O.; Kyselý, J. Hot Central-European summer of 2013 in a long-term context. Int. J. Climatol. 2015, 35, 4399–4407. [Google Scholar] [CrossRef]
  51. Brown, D.; Cabbage, M.; McCarthy, L. NASA, NOAA Analyses Reveal Record-Shattering Global Warm Temperatures in 2015. 2016. Available online: www.nasa.gov/press-release/nasa-noaa-analyses-reveal-record-shattering-global-warm-temperatures-in-2015 (accessed on 20 March 2017).
  52. Russo, S.; Sillmann, J.; Fischer, E.M. Top ten European heatwaves since 1950 and their occurrence in the future. Environ. Res. Lett. 2015, 10, 124003. [Google Scholar] [CrossRef]
  53. Hoy, A.; Hänsel, S.; Skalák, P.; Ustrnul, Z.; Bochníček, O. The extreme European summer of 2015 in a long-term perspective. Int. J. Climatol. 2016. [Google Scholar] [CrossRef]
  54. Vicedo-Cabrera, A.; Ragettli, M.; Schindler, C.; Röösli, M. Excess mortality during the warm summer of 2015 in Switzerland. Swiss Med. Wkly. 2016, 146, w14379. [Google Scholar] [CrossRef] [PubMed]
  55. Green, H.K.; Andrews, N.; Armstrong, B.; Bickler, G.; Pebody, R. Mortality during the 2013 heatwave in England—How did it compare to previous heatwaves? A retrospective observational study. Environ. Res. 2016, 147, 343–349. [Google Scholar] [CrossRef] [PubMed]
  56. Koppe, C.; Jendritzky, G. Inclusion of short-term adaptation to thermal stresses in a heat load warning procedure. Meteorol. Z. 2005, 14, 271–278. [Google Scholar] [CrossRef]
  57. Rocklöv, J.; Ebi, K.; Forsberg, B. Mortality related to temperature and persistent extreme temperatures: A study of cause-specific and age-stratified mortality. Occup. Environ. Med. 2011, 68, 531–536. [Google Scholar] [CrossRef] [PubMed]
  58. Rocklöv, J.; Forsberg, B.; Meister, K. Winter mortality modifies the heat-mortality association the following summer. Eur. Respir. J. 2009, 33, 245–251. [Google Scholar] [CrossRef] [PubMed]
  59. Hajat, S.; Kovats, R.S.; Lachowycz, K. Heat-related and cold-related deaths in England and Wales: Who is at risk? Occup. Environ. Med. 2007, 64, 93–100. [Google Scholar] [CrossRef] [PubMed]
  60. Hondula, D.M.; Davis, R.E.; Leisten, M.J.; Saha, M.V.; Veazey, L.M.; Wegner, C.R. Fine-scale spatial variability of heat-related mortality in Philadelphia County, USA, from 1983–2008: A case-series analysis. Environ. Health 2012, 11, 16. [Google Scholar] [CrossRef] [PubMed]
  61. Maier, G.; Grundstein, A.; Jang, W.; Li, C.; Naeher, L.P.; Shepherd, M. Assessing the performance of a vulnerability index during oppressive heat across Georgia, United States. Weather Clim. Soc. 2014, 6, 253–263. [Google Scholar] [CrossRef]
  62. Kovach, M.M.; Konrad, C.E., II; Fuhrmann, C.M. Area-level risk factors for heat-related illness in rural and urban locations across North Carolina, USA. Appl. Geogr. 2015, 60, 175–183. [Google Scholar] [CrossRef]
  63. O’Neill, M.S.; Carter, R.; Kish, J.K.; Gronlund, C.J.; White-Newsome, J.L.; Manarolla, X.; Zanobetti, A.; Schwartz, J.D. Preventing heat-related morbidity and mortality: New approaches in a changing climate. Maturitas 2009, 64, 98–103. [Google Scholar] [CrossRef] [PubMed]
  64. Urban, A.; Davídkovová, H.; Kyselý, J. Heat- and cold-stress effects on cardiovascular mortality and morbidity among urban and rural populations in the Czech Republic. Int. J. Biometeorol. 2014, 58, 1057–1068. [Google Scholar] [CrossRef] [PubMed]
  65. Kyselý, J.; Plavcová, E.; Davídkovová, H.; Kynčl, J. Comparison of hot and cold spell effects on cardiovascular mortality in individual population groups in the Czech Republic. Clim. Res. 2011, 49, 113–129. [Google Scholar] [CrossRef]
  66. Hanzlíková, H.; Plavcová, E.; Kynčl, J.; Kříž, B.; Kyselý, J. Contrasting patterns of hot spell effects on morbidity and mortality for cardiovascular diseases in the Czech Republic, 1994–2009. Int. J. Biometeorol. 2015, 59, 1673–1684. [Google Scholar] [CrossRef] [PubMed]
  67. Urban, A.; Burkart, K.; Kyselý, J.; Schuster, C.; Plavcová, E.; Hanzlíková, H.; Štepanek, P.; Lakes, T. Spatial patterns of heat-related cardiovascular mortality in the Czech Republic. Int. J. Environ. Res. Public Health 2016, 13, 284. [Google Scholar] [CrossRef] [PubMed]
  68. Ahmad, O.; Boschi-Pinto, C.; Lopez, A.D.; Murray, C.J.L.; Lozano, R.; Inoue, M. Age Standardization of Rates: A New WHO Standard; World Health Organization: Geneva, Switzerland, 2001; Available online: http://www.who.int/healthinfo/paper31.pdf (accessed on 20 March 2017).
  69. Jendritzky, G.; de Dear, R.; Havenith, G. UTCI—Why another thermal index? Int. J. Biometeorol. 2012, 56, 421–428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Morabito, M.; Crisci, A.; Messeri, A.; Capecchi, V.; Modesti, P.A.; Gensini, G.F.; Orlandini, S. Environmental temperature and thermal indices: What is the most effective predictor of heat-related mortality in different geographical contexts? Sci. World J. 2014, 961750. [Google Scholar] [CrossRef] [PubMed]
  71. Ruuhela, R.; Jylhä, K.; Lanki, T.; Tiittanen, P.; Matzarakis, A. Biometeorological assessment of mortality related to extreme temperatures in Helsinki region, Finland, 1972–2014. Int. J. Environ. Res. Public Health 2017, 14, 944. [Google Scholar] [CrossRef] [PubMed]
  72. Burkart, K.; Schneider, A.; Breitner, S.; Khan, M.H.; Krämer, A.; Endlicher, W. The effect of atmospheric thermal conditions and urban thermal pollution on all-cause and cardiovascular mortality in Bangladesh. Environ. Pollut. 2011, 159, 2035–2043. [Google Scholar] [CrossRef] [PubMed]
  73. Vaneckova, P.; Neville, G.; Tippett, V.; Aitken, P.; FitzGerald, G.; Tong, S. Do biometeorological indices improve modeling outcomes of heat-related mortality? J. Appl. Meteorol. Climatol. 2011, 50, 1165–1176. [Google Scholar] [CrossRef]
  74. Urban, A.; Kyselý, J. Comparison of UTCI with other thermal indices in the assessment of heat and cold effects on cardiovascular mortality in the Czech Republic. Int. J. Environ. Res. Public Health 2014, 11, 952–967. [Google Scholar] [CrossRef] [PubMed]
  75. McGregor, G.R.; Bessemoulin, P.; Ebi, K.L.; Menne, B. Heatwaves and Health: Guidance on Warning-System Development; World Meteorological Organization: Geneva, Switzerland; World Health Organization: Geneva, Switzerland, 2015; ISBN 978-92-63-11142-5. [Google Scholar]
  76. Výberči, D.; Labudová, L.; Eštóková, M.; Faško, P.; Trizna, M. Human mortality impacts of the 2015 summer heat spells in Slovakia. Theor. Appl. Climatol. 2017. [Google Scholar] [CrossRef]
  77. Nairn, J.R.; Fawcett, R.J.B. The Excess Heat Factor: A metric for heatwave intensity and its use in classifying heatwave severity. Int. J. Environ. Res. Public Health 2015, 12, 227–253. [Google Scholar] [CrossRef] [PubMed]
  78. Wood, S. Generalized Additive Models: An Introduction with R; Chapman & Hall/CRC: London, UK, 2006; p. 392. ISBN 1-58488-474-6. [Google Scholar]
  79. Hastie, T.; Tibshirani, R.; Leisch, F.; Hornik, K.; Ripley, B.D. Mda: Mixture and Flexible Discriminant Analysis. R Package Version 0.4-9. Stanford University: Stanford, CA, USA, 2016. Available online: https://cran.r-project.org/web/packages/mda/index.html (accessed on 20 March 2017).
  80. Peng, R.D.; Dominici, F. Statistical Methods for Environmental Epidemiology with R; Springer: New York, NY, USA, 2008; p. 144. ISBN 978-0-387-78166-2. [Google Scholar]
  81. Schoenberg, B.S. Calculating confidence intervals for rates and ratios. Neuroepidemiology 1983, 2, 257–265. [Google Scholar] [CrossRef]
  82. Saha, M.V.; Davis, R.E.; Hondula, D.M. Mortality displacement as a function of heat event strength in 7 US cities. Am. J. Epidemiol. 2014, 179, 467–474. [Google Scholar] [CrossRef] [PubMed]
  83. Armstrong, B.; Gasparrini, A.; Hajat, S. Estimating mortality displacement during and after heat waves. Am. J. Epidemiol. 2014, 179, 1405–1406. [Google Scholar] [CrossRef] [PubMed]
  84. Burkart, K.; Canário, P.; Breitner, S.; Schneider, A.; Scherber, K.; Andrade, H.; Alcoforado, M.J.; Endlicher, W. Interactive short-term effects of equivalent temperature and air pollution on human mortality in Berlin and Lisbon. Environ. Pollut. 2013, 183, 54–63. [Google Scholar] [CrossRef] [PubMed]
  85. Ha, J.; Kim, H.; Hajat, S. Effect of previous-winter mortality on the association between summer temperature and mortality in South Korea. Environ. Health Perspect. 2011, 119, 542–546. [Google Scholar] [CrossRef] [PubMed]
  86. Davis, R.E.; Dougherty, E.; McArthur, C.; Huang, Q.S.; Baker, M.G. Cold, Dry Air is Associated with Influenza and Pneumonia Mortality in Auckland, New Zealand. Influenza Other Respir. Viruses 2015, 10, 310–313. [Google Scholar] [CrossRef] [PubMed]
  87. Buckley, J.P.; Samet, J.M.; Richardson, D.B. Commentary: Does air pollution confound studies of temperature? Epidemiology 2014, 25, 242–245. [Google Scholar] [CrossRef] [PubMed]
  88. Sheridan, S.C.; Kalkstein, A.J.; Kalkstein, L.S. Trends in heat-related mortality in the United States, 1975–2004. Nat. Hazards 2009, 50, 145–160. [Google Scholar] [CrossRef]
  89. U.S. Energy Information Administration: Residential Energy Consumption Survey (RECS 2009). Available online: https://www.eia.gov/consumption/residential/reports/2009/air-conditioning.php (accessed on 28 October 2017).
  90. CZSO—Czech Statistical Office. 2017. Available online: https://www.czso.cz/staticke/animgraf/cz/ (accessed on 20 March 2017).
  91. Ebi, K.L.; Exuzides, K.A.; Lau, E.; Kelsh, M.; Barnston, A. Weather changes associated with hospitalizations for cardiovascular diseases and stroke in California, 1983–1998. Int. J. Biometeorol. 2004, 49, 48–58. [Google Scholar] [CrossRef] [PubMed]
  92. 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]
  93. Morabito, M.; Profili, F.; Crisci, A.; Francesconi, P.; Gensini, G.F.; Orlandini, S. Heat-related mortality in the Florentine area (Italy) before and after the exceptional 2003 heat wave in Europe: An improved public health response? Int. J. Biometeorol. 2012, 56, 801–810. [Google Scholar] [CrossRef] [PubMed]
  94. Lowe, D.; Ebi, K.L.; Forsberg, B. Heatwave early warning systems and adaptation advice to reduce human health consequences of heatwaves. Int. J. Environ. Res. Public Health 2011, 8, 4623–4648. [Google Scholar] [CrossRef] [PubMed]
  95. Ward, K.; Lauf, S.; Kleinschmit, B.; Endlicher, W. Heat waves and urban heat islands in Europe: A review of relevant drivers. Sci. Total Environ. 2016, 569–570, 527–539. [Google Scholar] [CrossRef] [PubMed]
  96. Scalley, B.D.; Spicer, T.; Jian, L.; Xiao, J.; Nairn, J.; Robertson, A.; Weeramanthri, T. Responding to heatwave intensity: Excess Heat Factor is a superior predictor of health service utilisation and a trigger for heatwave plans. Aust. N. Z. J. Public Health 2015, 39, 582–587. [Google Scholar] [CrossRef] [PubMed]
  97. Montero, J.C.; Mirón, I.J.; Criado, J.J.; Linares, C.; Díaz, J. Comparison between two methods of defining heat waves: A retrospective study in Castile-La Mancha (Spain). Sci. Total Environ. 2010, 408, 1544–1550. [Google Scholar] [CrossRef] [PubMed]
Ijerph 14 01562 g001 550
Figure 1. Location of the Czech Republic and meteorological stations used to calculate mean temperature (black dots).
Figure 1. Location of the Czech Republic and meteorological stations used to calculate mean temperature (black dots).
Ijerph 14 01562 g001
Ijerph 14 01562 g002 550
Figure 2. Total duration of heat waves in each year during 1994–2015 in the Czech Republic, and their total heat load according to the cumulative excess heat factor (∑EHF).
Figure 2. Total duration of heat waves in each year during 1994–2015 in the Czech Republic, and their total heat load according to the cumulative excess heat factor (∑EHF).
Ijerph 14 01562 g002
Ijerph 14 01562 g003 550
Figure 3. Relative excess mortality during summers (June–August) in six years with the total heat wave (HW) duration greater than 15 days. Tmean indicates daily mean temperature for the Czech Republic. A horizontal dotted line shows the heat wave threshold temperature—95th percentile of the annual daily mean temperature distribution during 1994–2015.
Figure 3. Relative excess mortality during summers (June–August) in six years with the total heat wave (HW) duration greater than 15 days. Tmean indicates daily mean temperature for the Czech Republic. A horizontal dotted line shows the heat wave threshold temperature—95th percentile of the annual daily mean temperature distribution during 1994–2015.
Ijerph 14 01562 g003
Ijerph 14 01562 g004 550
Figure 4. Relationship between cumulative excess mortality for the whole population and (a) cumulative excess heat factor (∑EHF), (b) average daily mean temperature (average Tmean), (c) maximum daily mean temperature (maximum Tmean), and (d) duration of heat waves, during 1994–2015 in the Czech Republic. Dash-dotted lines indicate 95% confidence interval of the regression line. Slopes of the fitted regression lines are presented in Table 3.
Figure 4. Relationship between cumulative excess mortality for the whole population and (a) cumulative excess heat factor (∑EHF), (b) average daily mean temperature (average Tmean), (c) maximum daily mean temperature (maximum Tmean), and (d) duration of heat waves, during 1994–2015 in the Czech Republic. Dash-dotted lines indicate 95% confidence interval of the regression line. Slopes of the fitted regression lines are presented in Table 3.
Ijerph 14 01562 g004
Table
Table 1. Characteristics of heat waves within six years with the longest total heat-wave duration during 1994–2015 in the Czech Republic.
Table 1. Characteristics of heat waves within six years with the longest total heat-wave duration during 1994–2015 in the Czech Republic.
Heat WaveStartEndDaysavgTmean (°C)maxTmean (°C)∑EHIsig (°C)∑EHIaccl (°C)∑EHF (°C)
1994_126 June29 June424.926.09.134.999.4
_22 July4 July322.824.20.714.210.0
_312 July17 July622.023.22.815.412.2
_422 July8 August1824.126.951.549.3176.6
2003_110 June12 June323.224.33.616.921.1
_219 July22 July422.825.03.115.220.0
_31 August10 August1023.325.417.130.772.7
_412 August14 August323.525.85.14.48.2
2006_119 June22 June422.723.83.332.234.4
_224 June26 June322.923.71.019.212.7
_310 July13 July423.124.76.29.615.9
_418 July28 July1124.125.526.530.499.7
2010_19 June12 June423.124.25.838.359.9
_229 June4 July622.323.53.729.425.4
_39 July17 July924.526.324.148.3163.7
2013_117 June21 June524.626.513.552.8158.7
_226 July29 July425.428.513.224.590.0
_31 August8 August824.126.717.722.565.2
2015_11 July7 July724.426.916.046.3143.1
_216 July25 July1024.227.626.048.2145.1
_33 August15 August1326.028.357.655.3292.4
_427 August1 September622.525.48.14.926.2
Table
Table 2. Impact of heat waves on mortality for the whole population, and the younger (0–64 years) and the elderly (65+ years) age groups, within six years with the longest total heat-wave duration during 1994–2015. The variables represent sum of excess deaths per standardized 10,000,000 inhabitants (excess mortality) and its 95% confidence intervals (CI), mean relative mortality deviation (ØRMD), and cumulative relative mortality deviation (∑RMD) during heat waves. Displaced mortality represents the absolute value of the ratio of the sum of three-day-averaged mortality deviations during the negative phase of the extended heat wave period (EHP) to the sum of three-day-averaged mortality deviations during the positive phase (see Section 2.5.3). --- indicates that the heat wave was considered as a single EHP together with the following one. * denotes statistically significant excess mortality.
Table 2. Impact of heat waves on mortality for the whole population, and the younger (0–64 years) and the elderly (65+ years) age groups, within six years with the longest total heat-wave duration during 1994–2015. The variables represent sum of excess deaths per standardized 10,000,000 inhabitants (excess mortality) and its 95% confidence intervals (CI), mean relative mortality deviation (ØRMD), and cumulative relative mortality deviation (∑RMD) during heat waves. Displaced mortality represents the absolute value of the ratio of the sum of three-day-averaged mortality deviations during the negative phase of the extended heat wave period (EHP) to the sum of three-day-averaged mortality deviations during the positive phase (see Section 2.5.3). --- indicates that the heat wave was considered as a single EHP together with the following one. * denotes statistically significant excess mortality.
Heat WaveWhole Population0–64 Years65+ Years
Excess Mortality (n)ØRMD (%)∑RMD (%)Displaced Mortality (%)Excess Mortality (n)ØRMD (%)∑RMD (%)Excess Mortality (n)ØRMD (%)∑RMD (%)
1994_1552 (455; 654) *26.8107.2---63 (17; 115) *10.943.7487 (403; 577) *32.7130.8
_289 (12; 170) *5.817.318.453 (12; 98) *12.336.936 (−28; 105)3.39.8
_3144 (36; 257) *4.828.928.321 (−35; 80)2.414.3123 (32; 220) *5.734.1
_41197 (1002; 1397) *13.4240.326.6352 (250; 459) *14.1254.4844 (678; 1014) *13.0234.1
2003_140 (−31; 115)3.08.921.819 (−17; 59)5.416.121 (−40; 85)2.06.1
_2265 (181; 355) *15.561.88.494 (51; 142) *20.782.9171 (99; 248) *13.554.0
_3254 (124; 388) *6.059.7---47 (−18; 116)4.342.9203 (92; 319) *6.564.5
_421 (−49; 94)1.64.743.512 (−23; 51)3.310.010 (−49; 74)1.03.1
2006_1249 (168; 336) *15.662.42.460 (19; 105) *14.156.3190 (120; 265) *16.164.5
_2−11 (−76; 60)−1.0−3.0200.9−27 (−57; 9)−8.8−26.315 (−42; 75)1.54.6
_3191 (111; 276) *12.048.117.957 (17; 101) *13.855.2135 (66; 208) *11.445.7
_4369 (237; 506) *8.593.436.766 (2; 136) *5.965.3308 (194; 427) *9.5105.0
2010_1211 (133; 293) *14.457.41.835 (−3; 78)9.236.9178 (110; 250) *16.365.1
_2101 (10; 197) *4.627.71.9152 (6; 103) *9.355.551 (−28; 133)3.118.5
_3471 (354; 593) *14.4129.95.2105 (47; 168) *12.4111.9366 (265; 472) *15.1136.2
2013_1418 (330; 511) *24.3121.737.396 (54; 143) *23.0115.1325 (248; 407) *25.0125.0
_2259 (184; 340) *19.577.8---81 (43; 123) *23.794.7176 (112; 246) *17.770.8
_3116 (15; 222) *4.334.633.5−25 (−73; 29)−3.6−29.0139 (51; 231) *7.055.6
2015_1304 (208; 405) *13.695.09.591 (43; 144) *16.1112.9217 (134; 304) *12.990.3
_2160 (49; 276) *5.050.324.212 (−42; 70)1.615.6148 (52; 249) *6.261.7
_3847 (711; 988) *20.4265.422.4105 (42; 174) *10.2132.2750 (630; 874) *24.1313.3
_4211 (117; 311) *9.465.854.955 (9; 105) *9.767.8158 (76; 245) *9.365.2
Table
Table 3. Relationships of cumulative excess mortality for the whole population, the younger (0–64 years) and the elderly (65+ years) age groups to heat wave characteristics (independent variables x), namely cumulative excess heat factor (∑EHF), cumulative significance index (∑EHIsig), cumulative acclimatization index (∑EHIaccl), average daily mean temperature (avgTmean), maximum daily mean temperature (maxTmean), and length of heat waves (duration). Regression coefficients are reported along with their associated R2. All relationships were statistically significant at p < 0.001.
Table 3. Relationships of cumulative excess mortality for the whole population, the younger (0–64 years) and the elderly (65+ years) age groups to heat wave characteristics (independent variables x), namely cumulative excess heat factor (∑EHF), cumulative significance index (∑EHIsig), cumulative acclimatization index (∑EHIaccl), average daily mean temperature (avgTmean), maximum daily mean temperature (maxTmean), and length of heat waves (duration). Regression coefficients are reported along with their associated R2. All relationships were statistically significant at p < 0.001.
VariableWhole Population0–64 Years65+ Years
SlopeR2SlopeR2SlopeR2
∑EHF (°C)0.810.760.540.380.910.80
∑EHIsig (°C)3.970.722.760.394.390.74
∑EHIaccl (°C)2.990.572.330.393.230.56
avgTmean (°C)31.90.3116.30.0937.50.36
maxTmean (°C)21.20.3111.80.1124.60.39
duration (days)13.90.5911.40.4514.70.56

Share and Cite

MDPI and ACS Style

Urban, A.; Hanzlíková, H.; Kyselý, J.; Plavcová, E. Impacts of the 2015 Heat Waves on Mortality in the Czech Republic—A Comparison with Previous Heat Waves. Int. J. Environ. Res. Public Health 2017, 14, 1562. https://doi.org/10.3390/ijerph14121562

AMA Style

Urban A, Hanzlíková H, Kyselý J, Plavcová E. Impacts of the 2015 Heat Waves on Mortality in the Czech Republic—A Comparison with Previous Heat Waves. International Journal of Environmental Research and Public Health. 2017; 14(12):1562. https://doi.org/10.3390/ijerph14121562

Chicago/Turabian Style

Urban, Aleš, Hana Hanzlíková, Jan Kyselý, and Eva Plavcová. 2017. "Impacts of the 2015 Heat Waves on Mortality in the Czech Republic—A Comparison with Previous Heat Waves" International Journal of Environmental Research and Public Health 14, no. 12: 1562. https://doi.org/10.3390/ijerph14121562

APA Style

Urban, A., Hanzlíková, H., Kyselý, J., & Plavcová, E. (2017). Impacts of the 2015 Heat Waves on Mortality in the Czech Republic—A Comparison with Previous Heat Waves. International Journal of Environmental Research and Public Health, 14(12), 1562. https://doi.org/10.3390/ijerph14121562

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