The study, published in the journal Advances in Statistical Climatology, Meteorology and Oceanography by Pascal Meurer, Sebastian Buschow, Svenja Szemkus and Petra Friederichs, tackles a blind spot in conventional attribution. Most attribution analyses compress a sprawling, spatially extended heatwave into a single number, such as a regional average temperature, and then apply univariate extreme value statistics. In doing so they discard the temporal structure of the event, the fact that extreme heat persists for days, that heatwaves unfold as sequences of hot days and nights, and that this persistence is central to their impact on human health, agriculture and infrastructure. The Bonn team set out to build an attribution framework that keeps that temporal dependence intact and can even attribute summers that mix heat extremes with cold spells, such as the varied European summer of 2025.

The first challenge is dimensionality. Daily maximum temperature over Europe is a vast spatial field, and an attribution method that must evaluate probabilities for extremes cannot handle thousands of grid points at once. The researchers therefore condense the spatial information using the extremal pattern index, or EPI, a tool developed by Szemkus and Friederichs in 2024. The method is a cousin of principal component analysis, but instead of decomposing ordinary correlations it decomposes the dependence between extremes. Temperature anomalies are first transformed to a standard Fréchet distribution, a mathematical step that deliberately gives most weight to large positive anomalies and little weight to small or negative ones. A tail pairwise dependence matrix is then built from all pairs of land grid points, and its eigenvectors reveal spatial patterns where extremes tend to occur together. Projecting each day’s temperature field onto the ten leading extremal patterns yields a single daily number, the EPI, which summarises how strongly a large-scale extreme heat pattern is active. In southern Europe, these ten patterns capture roughly seventy-five percent of the variance of the Fréchet-standardised extremes.

With the spatial field compressed into a daily time series, the team models the sequence of EPI values across each summer as a first-order Markov process, meaning that the state of any given day depends statistically on the previous day. The crucial innovation lies in how the extremes are handled. Rather than the traditional practice of declustering, in which consecutive threshold exceedances are collapsed into a single cluster maximum to manufacture independence, the researchers model the dependence directly using bivariate extreme value theory. Pairs of consecutive days are described by a joint distribution whose margins follow a generalised Pareto distribution above a high threshold, the standard peaks-over-threshold framework, while the dependence between the two days is captured by a stable tail dependence function, here the simple logistic model with a single parameter. A censored likelihood scheme divides the two-dimensional plane into four regions depending on whether each day exceeds the threshold, so that non-extreme observations contribute only the probability of not exceeding, and the estimation of extremal dependence remains untainted by ordinary weather.

Because the climate is not stationary, the model must also evolve. The parameters of the marginal distribution, the exceedance probability, the scale and the shape of the generalised Pareto distribution, and the day-to-day dependence parameter are all allowed to change slowly over time, described by low-order Legendre polynomials that are orthogonal, smooth and independent of the data. The team implements two complementary formulations. In the first, the threshold is held fixed at the ninety-fifth percentile and the exceedance probability is modelled with logistic regression, so that a warming world shows up as a rising probability of crossing the same line. In the second, the exceedance probability is fixed and the threshold itself is allowed to drift, estimated by quantile regression, so that the general warming trend is absorbed into a moving baseline. This second formulation answers a subtly different question: is there any climate change signal in the tail behaviour of heatwaves beyond the simple shift of the whole temperature distribution?

Attribution then proceeds through a likelihood ratio grounded in causal counterfactual theory. The researchers fit the non-stationary Markov model separately to two families of climate simulations from the CMIP6 archive: the historical scenario with both natural and anthropogenic forcing, extended after 2014 by the SSP2-4.5 scenario, and the historical-natural scenario in which greenhouse gases and aerosols from human activity are removed and only solar and volcanic drivers remain. The likelihood of the observed EPI time series from the ERA5 reanalysis, which provides a consistent picture of the atmosphere since 1940, is then evaluated under each scenario. Because summers are treated as independent, the seasonal likelihood ratios multiply into a Bayes factor, and the accumulated evidence over any span of summers can be tracked from 1940 to the present. Uncertainty is handled by bootstrapping ensemble members within each climate model, and the log-likelihood ratios from different models are combined with a random-effects meta-analytic model that separates natural variability from model uncertainty.

The results are unambiguous for the first attribution question. The full record of European summers from 1940 to 2025 is attributed to the scenario with anthropogenic emissions with decisive evidence in every region examined, northern, central and southern Europe. In the southern European region, the accumulated Bayes factor after the summer of 2025 reaches roughly 2.25 times ten to the power of twenty-six, meaning the observed summers are that many times more probable under the world with human emissions than under the counterfactual natural-only world. More remarkable still is the timing: with a fixed threshold, the level of decisive evidence was already crossed in 1962 for central Europe and 1964 for southern Europe. Individual recent summers stand out sharply, with 2010, 2018 and 2022 providing strong evidence even on their own in northern Europe, and the scorching southern European summer of 2025, in which nearly half the days exceeded the ninety-fifth percentile of the EPI, falling into the decisive category. Interestingly, the summer of 2025 did not top the list despite the exceptional Scandinavian heatwave, because the method attributes the entire season, and a relatively cool June diluted the evidence compared with attribution studies that examined only the hot July phase.

The answer to the second question is more nuanced and, in a sense, reassuring in one narrow respect. When the warming trend is absorbed into a time-varying threshold, the likelihood ratios of most recent summers hover around one, indicating no detectable signal beyond the shift of the distribution. In the anthropogenic simulations, the shape parameter of the exceedance distribution trends negative after 2010, which would imply a bounded upper limit on temperatures relative to the shifted baseline, but this tendency is not reproduced in the ERA5 observations and may be an artefact of the climate models. In other words, relative to a moving baseline, the tail behaviour of heatwave extremes has not changed in a consistent, detectable way. The authors are careful to stress that this statement does not contradict the first finding: the increasing frequency and intensity of European heatwaves remains a robust and decisive result, and the dependence between consecutive hot days has clearly strengthened over time, reflecting longer and more persistent heatwaves.

The broader significance of the work extends beyond the headline numbers. By adapting the censored threshold model for extremes to non-stationary conditions, the Bonn team has created a framework that can attribute any stretch of time, from a single summer to a multi-decade record, and that is less dependent on how an event is defined than conventional approaches. The method also doubles as a test of climate models themselves: comparing the non-stationary parameter estimates from CMIP6 simulations with those from ERA5 reveals how faithfully the models reproduce the historical evolution of temperature extremes, and the historical simulations generally track the observed trends well. As warming accelerates and non-stationarity intensifies, the authors argue, sophisticated statistical machinery of this kind will become essential for evaluating extremes in both past and future climates, complementing the rapid-response analyses of the World Weather Attribution with a perspective that captures the full temporal texture of a heatwave rather than a single compressed snapshot.

Subject of Research: Non-stationary time series attribution of European heatwaves to anthropogenic climate forcing

Article Title: Non-stationary time series attribution for heatwaves over Europe

Article References: Meurer, P., Buschow, S., Szemkus, S., & Friederichs, P. (2026). Non-stationary time series attribution for heatwaves over Europe. Advances in Statistical Climatology, Meteorology and Oceanography, 12(2), 243-274. https://doi.org/10.5194/ascmo-12-243-2026

Image Credits: AI Generated

DOI: 10.5194/ascmo-12-243-2026

Keywords: heatwaves, extreme event attribution, climate change, extreme value theory, Markov process, CMIP6, ERA5, Europe, Bayes factor, non-stationarity, greenhouse gases, statistical climatology

The methodological innovation at the heart of the study lies in what the authors call health impact event attribution. Traditional extreme event attribution, which has matured rapidly over the past decade, typically asks how climate change altered the probability of a meteorological anomaly such as a heatwave, using metrics like regional temperature maxima. But a hotter day does not translate directly into a death, and the relationship between temperature and mortality varies enormously with geography, climate, age structure and adaptation. The new approach closes that gap by coupling attribution statistics to epidemiological exposure-response functions calibrated for each of 34 European countries. Rather than asking how much more likely the 2022 heatwave was, the team asked how much more likely the mortality event itself was, the actual human toll, in a climate warmed by roughly 1.2 degrees Celsius compared with a pre-industrial baseline.

Technically, the researchers constructed generalized extreme value distributions for both temperature and heat-related mortality, allowing them to estimate return periods, how often an event of a given magnitude would be expected to occur, under current and counterfactual pre-industrial conditions. The ratio of these probabilities, known as the probability ratio, then expresses how much climate change has multiplied the odds of an extreme. Crucially, the mortality side of the analysis was built on weekly death counts from Eurostat, publicly available data that are updated in near real time. Temperature inputs came from the ERA5-Land reanalysis produced by the European Centre for Medium-Range Weather Forecasts, while global mean surface temperature changes were drawn from NASA’s GISTEMP record. The reliance on weekly rather than daily mortality data was validated using recent methodological work showing that temporally aggregated health data can yield unbiased temperature-mortality estimates, a finding that dramatically expands where such analyses can be performed.

The results reveal a continent divided. Across Europe as a whole, 2022-like heat-mortality events became substantially more probable, but the amplification in Southern Europe dwarfed the continental average. The Mediterranean region, already identified in the scientific literature as a climate change hotspot, experienced nonlinear increases in the likelihood of deadly heat, meaning that each additional fraction of a degree of warming multiplies the risk faster than the last. The study also projected forward, estimating probability ratios at 1.5 and 2.0 degrees Celsius of global warming, the benchmark thresholds of the Paris Agreement. Those projections show the risk continuing to climb steeply, with Southern Europe again bearing the steepest increases. The nonlinear character of these changes is particularly troubling because it implies that modest additional warming will produce disproportionately large jumps in mortality risk, compressing the time available for societies to adapt.

Equally significant are the demographic disparities the analysis uncovered. When the researchers stratified their models by age, they found that the elderly, particularly those aged 80 and older, faced sharply elevated probability ratios compared with younger populations. This pattern reflects well-documented physiological vulnerabilities, including diminished thermoregulatory capacity, higher prevalence of chronic disease, and greater use of medications that impair the body’s ability to cope with heat stress. Prior research has also documented sex-specific differences in heat vulnerability and adaptation, and the new framework’s ability to dissect such heterogeneity is one of its principal strengths. The authors emphasize that these findings make the case for population-specific analyses: a national average can conceal the fact that a small, highly vulnerable subgroup absorbs a wildly disproportionate share of the climate-driven mortality increase.

The 2022 summer serves as the study’s anchor case because it was, by most measures, the hottest summer ever recorded in Europe. Spain registered its warmest summer in its historical series, the United Kingdom breached 40 degrees Celsius for the first time, and the Copernicus Climate Change Service documented the season as a continental anomaly. Previous work by several of the same authors estimated that the summer of 2022 caused more than 60,000 heat-related deaths across Europe, while companion studies attributed a substantial fraction of that toll directly to anthropogenic warming. The new attribution analysis reframes those mortality counts in probabilistic terms: without human influence on the climate, an event of that lethality would have been extraordinarily rare, whereas in today’s climate it has become a recurring hazard with a return period measured in years rather than centuries.

What distinguishes this work from earlier impact attribution efforts is not only its scope, spanning 34 countries, but its speed and reproducibility. Because the mortality data are public and updated weekly, and because the analytical code has been released openly on GitHub, the framework can in principle be deployed rapidly after any extreme heat episode, transforming attribution from a retrospective academic exercise into a near-real-time public health tool. The authors argue that this capability could feed directly into health emergency forecasting systems, several of which are already being developed and tested in Europe, allowing authorities to anticipate and respond to deadly heat before the peak toll is counted. Rapid attribution of health impacts also carries weight beyond epidemiology, providing quantified evidence relevant to climate litigation, loss-and-damage negotiations and the legal and moral accounting of emissions.

The study arrives amid a broader scientific consensus that heat has become Europe’s leading weather-related killer. The Lancet Countdown and the World Meteorological Organization have both documented rising heat mortality, and modeling studies covering hundreds of European cities project that the burden will grow substantially without aggressive adaptation. Yet adaptation itself remains uneven. Research from Spain and the Netherlands shows that populations can shift their minimum mortality temperature over time, a signature of acclimatization, but the pace of such physiological and behavioral adaptation lags far behind the pace of warming. Heat-health warning systems exist across much of Europe but vary widely in coverage, trigger thresholds and effectiveness. The new findings sharpen the argument that adaptation investments, from cooling centers and urban greening to occupational heat protections and targeted outreach to the elderly, should be prioritized precisely where and for whom the probability ratios are highest.

Ultimately, the study delivers a double message. First, the fingerprint of climate change is now legible not just in thermometers but in mortality statistics: the deaths of thousands of Europeans in recent summers are statistically attributable to emissions, with Southern Europeans facing risks an order of magnitude beyond the continental average. Second, the tools needed to see that fingerprint are now fast, transparent and grounded in freely available data, meaning that the era of waiting years to understand the human cost of an extreme summer is over. As global temperatures continue to rise toward and potentially past the 1.5-degree threshold, the nonlinear escalation of heat-mortality risk documented in this analysis offers both a warning and a roadmap. The regions and populations identified as most vulnerable are known, the data streams needed to monitor the threat in real time exist, and the remaining variable is the speed with which societies choose to act on the evidence.

Subject of Research: Attribution of extreme heat-related mortality events in Europe to anthropogenic climate change

Article Title: Extreme event attribution for heat-related mortality due to anthropogenic climate change across Europe

Article References: Beck, T. M., Gudmundsson, L., Schumacher, D. L., Seneviratne, S. I., Achebak, H., & Ballester, J. (2026). Extreme event attribution for heat-related mortality due to anthropogenic climate change across Europe. Nature Health. https://doi.org/10.1038/s44360-026-00193-z

Image Credits: AI Generated

DOI: 10.1038/s44360-026-00193-z

Keywords: extreme event attribution, heat-related mortality, climate change, Europe, heatwaves, Nature Health, epidemiological models, Southern Europe, public health, adaptation, probability ratio, 2022 heatwave