According to the European Environment Agency, Europe is the fastest-warming continent in the world, and heat stress is one of the most severe climate risks.1 In 2023, Europe recorded one of its warmest summers in history. Globally, yearly average temperatures continue to rise as well. Many countries are experiencing an increase in their average temperature and temperature volatility, leading to extreme heat or cold. Thus, people are being exposed to more heatwaves, which have a significant impact on human health, as demonstrated by the high number of excess deaths during the European heatwaves in 2003, 2018 and 2022.2,3,4
To analyze the mortality impact of heatwaves, one should first define what a heatwave is – but due to regional variations in temperature, there is no standardized definition. Typically, heatwaves are defined by the number of days (duration) with the maximum temperature exceeding a certain threshold. The threshold can be either a fixed value or a percentile of the local temperature distribution.5 It is also important to consider the possibility of delayed effects of extreme heat on mortality.6 For example, heat stress or dehydration may worsen chronic conditions and thereby lead to higher mortality in the weeks following a heatwave.
In this article, we analyze the impact of heatwaves on mortality in selected European countries from 2000–2019 and in South Korea from 2010–2019. We estimate quasi-Poisson regression models7 with actual over expected (A/E) deaths as the target variable and include the number of days in a week with maximum temperatures at or above a fixed threshold to capture the effects of extreme heat duration. We also include one- and two-week lags of extreme heat duration to explore its delayed effects on the mean A/E.
Data Sources
We considered the following groups of countries:
- The DACH countries (Germany, Austria and Switzerland)
- The Benelux countries (Belgium, Netherlands and Luxembourg)
- The Baltic countries (Estonia, Latvia and Lithuania)
- South Korea
Since the land area of Germany is larger than the areas of the other selected countries, the analysis on a national level for Germany would have been less informative and was conducted at the state level instead.
We took the actual weekly deaths from the Short-term Mortality Fluctuations (STMF) data series of the Human Mortality Database (HMD) for all countries except Germany,8 for which we used weekly state-level deaths from the German Federal Statistical Office.9 Thereby, we obtained more data points for Germany and accounted for state variations in temperature. For the European countries, we included weekly mortality data over the period 2000–2019, preceding the COVID‑19 pandemic. For South Korea, the STMF weekly mortality data was available only for 2010–2019.
For all countries except Germany, we used the HMD population data at the beginning of a year.10 For Germany, population size at the end of a year for each state was taken from the German Federal Statistical Office,11 and we assumed that population size at the end of a year was equal to population size at the beginning of the next year. Further, we used life tables from the HMD12 for all countries, including Germany, to compute expected deaths for each country or state in each year with the population size at the beginning of the year and life table of the previous year, which was the newest life table at that time. We divided the expected yearly deaths by 52, removing week 53 where applicable for consistency purposes. Weeks in these data sources are defined according to the ISO 8601 date and time standard.
We used CPC Global Unified Temperature data,13 WorldPop population counts14 and country15 (or German state16) polygons to compute the weighted mean of maximum temperature for each country (or German state) for each day over the period 2000–2019.17 Population counts were used as weights, since high temperatures are expected to have a stronger effect on mortality in densely populated areas, as more people are at risk there.
Modeling Approach
We estimated quasi-Poisson regression models to investigate the impact of extreme heat duration on the mean value of actual deaths while accounting for expected deaths in different countries or groups of countries. To this end, the A/E ratio was taken as the target variable and expected deaths as the weight variable. The number of days with maximum temperatures at or above a certain threshold in a given week, one week prior and two weeks prior were included as categorical predictor variables alongside country (for Europe), year, age group and gender. We computed the 95% confidence intervals for parameter estimates using the delta method.18,19
Since extreme heat events typically occur during the summer months, we restricted the data we used for modeling to weeks 18–40. Higher levels of the extreme heat duration were combined into one category (3+, 4+ or 5+) due to a low occurrence of those levels, and we chose the level of 0 days as the base level in all models.
Extreme heat in Europe
Figure 1 presents modeling results for the DACH, Benelux and Baltic countries. The maximum temperature threshold used to define extreme heat was set to 30°C, similar to studies analyzing heat-related mortality in Germany and Austria.20,21 The results in the rest of the article are based on the state-level data for Germany and country-level data for the other countries.
Figure 1 – Effects of Extreme Heat Duration in the DACH, Benelux and Baltic Countries