Non-stationary time series attribution for heatwaves over Europe
Pith reviewed 2026-05-25 07:34 UTC · model grok-4.3
The pith
A non-stationary Markov model attributes observed European heatwave time series to anthropogenic forcing with high confidence.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The likelihood ratio of the ERA5 observational time series, given the distributions estimated from CMIP6 historical natural-only and natural-plus-anthropogenic simulations via a non-stationary Markov process and the extremal pattern index, supplies very strong evidence for the anthropogenic scenario over Europe since the beginning of the 21st century; for central and southern Europe the influence of anthropogenic greenhouse gas emissions on heatwaves is already detectable in the 1970s, while no reliable signal appears beyond the general temperature increase in either the temporal dependence of extreme heat days or the shape of the extreme value distribution.
What carries the argument
The extremal pattern index (EPI) condenses spatial fields of extremes into a single compact time series; this series is then modeled by a non-stationary Markov process built on bivariate extreme value theory, allowing direct computation of the likelihood ratio between the two CMIP6 forcing ensembles.
If this is right
- Attribution remains strong even after the mean warming trend is removed from the data.
- For central and southern Europe the anthropogenic greenhouse-gas signal on heatwaves is already present in the 1970s.
- No additional signal beyond the overall temperature rise is found in the temporal clustering of extremes or in the shape of their distribution.
- The same time-series likelihood framework can be applied to other extreme-event types that exhibit clustering.
Where Pith is reading between the lines
- The method could be tested on other mid-latitude regions where heatwave clustering is pronounced to see whether the same early detection holds.
- If the CMIP6 natural-only runs systematically underestimate the variance of real-world natural variability, the reported likelihood ratios would overstate the strength of the anthropogenic signal.
- The absence of change in the shape of the extreme-value distribution suggests that future attribution work may need to focus on frequency and duration rather than on shifts in tail heaviness.
Load-bearing premise
The probability distributions fitted to the CMIP6 simulations accurately represent the true distributions of the extremal pattern index under natural-only and natural-plus-anthropogenic conditions.
What would settle it
A direct calculation in which the likelihood of the ERA5 time series under the natural-only CMIP6 distributions equals or exceeds its likelihood under the anthropogenic distributions would falsify the attribution result.
Figures
read the original abstract
The increasing occurrence of extreme weather events since the beginning of the 21st century has led to the development of new methods to attribute extreme events to anthropogenic climate change. The way in which the extreme event is defined has a major influence on the attribution result. A frequently disregarded or overlooked aspect concerns the temporal dependence and the clustering of extremes. This study presents an approach for attributing complete time series during extreme events to anthropogenic forcing. The approach is based on a non-stationary Markov process using bivariate extreme value theory to model the temporal dependence of the time series. We calculate the likelihood ratio of an observational time series from ERA5 given the distributions as estimated from CMIP6 simulations with historical natural-only and natural and anthropogenic forcing scenarios. The spatial fields are condensed by the extremal pattern index (EPI) as a compact description of spatial extremes. In addition, the study examines the extent to which attribution statements about the occurrence of extreme heat events change when the effect of the mean warming is eliminated. The resulting attribution statement provides very strong evidence for the scenario with anthropogenic drivers over Europe, especially since the beginning of the 21st century. For central and southern Europe, the influence of anthropogenic greenhouse gas emissions on heatwaves could already have been proven in the 1970s using today's methods. There is no reliable signal apart from a general increase in temperature, neither in terms of the temporal dependence of extreme heat days nor in terms of the shape of the extreme value distribution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a non-stationary Markov process based on bivariate extreme value theory to attribute complete time series of heatwaves (via the extremal pattern index EPI) over Europe. It computes the likelihood ratio of the ERA5 observational series against distributions fitted from CMIP6 historical natural-only versus natural-plus-anthropogenic ensembles, concluding very strong evidence for anthropogenic forcing (detectable since the 1970s in central/southern Europe) with no residual signal beyond mean warming in temporal dependence or extreme-value shape.
Significance. If the CMIP6-fitted distributions are faithful representations of the true conditional distributions of the EPI process, the method would offer a useful extension of attribution techniques that explicitly incorporates temporal clustering and non-stationarity. The early-detection claim and the finding of no additional signal beyond mean warming would be of broad interest to the detection-and-attribution community.
major comments (2)
- [Likelihood ratio computation] The attribution statement rests on the premise that the probability distributions estimated from the CMIP6 hist-nat and hist ensembles accurately represent the true distributions of the EPI under each forcing (abstract; § on likelihood ratio). No validation against held-out data, cross-validation, or sensitivity tests to model bias in tail shape or Markov dependence parameters is reported; any systematic discrepancy directly biases the likelihood ratio and therefore the central claim.
- [Results on mean-warming removal] The additional test that removes mean warming and reports no residual signal in temporal dependence or extreme-value shape (abstract) inherits the identical requirement that the CMIP6-fitted distributions faithfully capture the conditional process; without evidence on this point the conclusion that there is 'no reliable signal apart from a general increase in temperature' cannot be assessed.
minor comments (2)
- [Abstract] The abstract states 'very strong evidence' but supplies no numerical likelihood-ratio values, confidence intervals, or error bars; these should be reported in the main text and abstract.
- [Methods] Notation for the non-stationary Markov process and the bivariate extreme-value dependence function should be defined explicitly with equation numbers in the methods section.
Simulated Author's Rebuttal
We thank the referee for their insightful comments on our manuscript. We address the major concerns regarding the validation of the CMIP6-fitted distributions below and outline planned revisions to strengthen the paper.
read point-by-point responses
-
Referee: [Likelihood ratio computation] The attribution statement rests on the premise that the probability distributions estimated from the CMIP6 hist-nat and hist ensembles accurately represent the true distributions of the EPI under each forcing (abstract; § on likelihood ratio). No validation against held-out data, cross-validation, or sensitivity tests to model bias in tail shape or Markov dependence parameters is reported; any systematic discrepancy directly biases the likelihood ratio and therefore the central claim.
Authors: We recognize that the attribution results depend on the assumption that the CMIP6 ensembles provide faithful representations of the EPI distributions under the respective forcings. The manuscript does not include explicit held-out validation or cross-validation, as the full ensembles are used to estimate the parameters. However, the use of multi-model ensembles is standard in detection and attribution studies to account for model uncertainty. To address this, we will add sensitivity analyses to variations in tail shape and Markov dependence parameters in the revised version, including tests using subsets of models. This addition will be included as a new paragraph in the methods section and discussed in the results. revision: partial
-
Referee: [Results on mean-warming removal] The additional test that removes mean warming and reports no residual signal in temporal dependence or extreme-value shape (abstract) inherits the identical requirement that the CMIP6-fitted distributions faithfully capture the conditional process; without evidence on this point the conclusion that there is 'no reliable signal apart from a general increase in temperature' cannot be assessed.
Authors: The mean-warming removal analysis is indeed subject to the same modeling assumptions. We will revise the manuscript to explicitly state that this conclusion is conditional on the accuracy of the fitted distributions and incorporate the sensitivity tests described in response to the first comment. This will help assess the robustness of the finding that no residual signal exists beyond mean warming. revision: partial
Circularity Check
No significant circularity; derivation uses independent CMIP6 fits
full rationale
The paper estimates distributions from separate CMIP6 historical natural-only and natural-plus-anthropogenic ensembles, then forms the likelihood ratio of the ERA5 EPI series against those fitted distributions. This step does not reduce by the paper's own equations to any parameter fitted directly to the ERA5 observations themselves. No self-citation load-bearing steps, self-definitional constructions, or fitted-input-called-prediction patterns appear in the derivation chain. The central attribution result therefore remains independent of the target data and is self-contained against the external simulation benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption CMIP6 historical natural-only and natural-plus-anthropogenic ensembles furnish unbiased estimates of the true distributions of the extremal pattern index under each forcing scenario
- domain assumption The extremal pattern index condenses spatial extreme fields without discarding information relevant to the attribution question
Reference graph
Works this paper leans on
-
[1]
Barnes, C., Clarke, B., Rantanen, M., Skålevåg, A., Ødemark, K., Kjellström, E., Vahlberg, M., Singh, R., Otto, F., Zachariah, M., and et al.: Intense two-week heatwave in Fennoscandia hotter and more likely due to climate change, https://doi.org//10.25560/122924,
-
[2]
Chavez-Demoulin, V . and Davison, A.: Modelling Time Series Extremes, REVSTAT-Statistical Journal, 10, 109–133, https://doi.org/10.57805/revstat.v10i1.113,
-
[3]
Cooley, D. and Thibaud, E.: Decompositions of dependence for high-dimensional extremes, Biometrika, 106, 587–604, https://doi.org/10.1093/biomet/asz028,
-
[4]
Copernicus Climate Change Service: Heatwaves contribute to the warmest June on record in western Europe, https://climate.copernicus.eu/ heatwaves-contribute-warmest-june-record-western-europe, accessed on 08/10/2025.,
work page 2025
-
[5]
Faranda, D., Guinaldo, T., Pastor, J. F., Alberti, T., and Khodayar, S.: Attribution of the 2025 Mediterranean Marine Heatwave to Climate Change Using Analogues, https://hal.science/hal-05289765, preprint,
work page 2025
-
[6]
Fawcett, L. and Walshaw, D.: Markov chain models for extreme wind speeds, Environmetrics, 17, 795–809, https://doi.org/https://doi.org/10.1002/env.794,
-
[7]
Friederichs, P.: Statistical downscaling of extreme precipitation events using extreme value theory, Extremes, 13, 109–132, https://doi.org/10.1007/s10687-010-0107-5,
-
[8]
P., Shiogama, H., Funke, B., Hegerl, G., Knutti, R., Matthes, K., Santer, B
Gillett, N. P., Shiogama, H., Funke, B., Hegerl, G., Knutti, R., Matthes, K., Santer, B. D., Stone, D., and Tebaldi, C.: The Detection and Attribution Model Intercomparison Project (DAMIP v1.0) contribution to CMIP6, Geoscientific Model Development, 9, 3685–3697, https://doi.org/10.5194/gmd-9-3685-2016,
-
[9]
Gulev, S., Thorne, P., Ahn, J., Dentener, F., Domingues, C., Gerland, S., Gong, D., Kaufman, D., Nnamchi, H., Quaas, J., Rivera, J., Sathyen- dranath, S., Smith, S., Trewin, B., von Schuckmann, K., and V ose, R.: Changing State of the Climate System, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessme...
-
[10]
Hannart, A. and Naveau, P.: Probabilities of Causation of Climate Changes, Journal of Climate, 31, 5507 – 5524, https://doi.org/10.1175/JCLI-D-17-0304.1,
-
[11]
Hannart, A., Pearl, J., Otto, F. E. L., Naveau, P., and Ghil, M.: Causal Counterfactual Theory for the Attribution of Weather and Climate- Related Events, Bulletin of the American Meteorological Society, 97, 99 – 110, https://doi.org/10.1175/BAMS-D-14-00034.1,
-
[12]
Hegerl, G. C., Karl, T. R., Allen, M., Bindoff, N. L., Gillett, N., Karoly, D., Zhang, X., and Zwiers, F.: Climate Change Detection and Attribution: Beyond Mean Temperature Signals, Journal of Climate, 19, 5058–5077, http://www.jstor.org/stable/26259286,
-
[13]
Hersbach, H., Bell, B., Berrisford, P., Hirahara, S., Horányi, A., Muñoz-Sabater, J., Nicolas, J., Peubey, C., Radu, R., Schepers, D., et al.: The ERA5 global reanalysis, Quarterly Journal of the Royal Meteorological Society, 146, 1999–2049, https://doi.org/10.1002/qj.3803,
-
[14]
32 Hulme, M.: Attributing weather extremes to ‘climate change’: A review, Progress in Physical Geography: Earth and Environment, 38, 499– 511, https://doi.org/10.1177/0309133314538644,
-
[15]
Huser, R., Opitz, T., and Thibaud, E.: Max-infinitely divisible models and inference for spatial extremes, Scandinavian Journal of Statistics, 48, 321–348, https://doi.org/https://doi.org/10.1111/sjos.12491,
-
[16]
IEA: Global Energy Review: CO2 Emissions in 2021 Global emissions rebound sharply to highest ever level, https://www.iea.org/reports/ global-energy-review-co2-emissions-in-2021-2, accessed: 2025-07-30,
work page 2021
-
[17]
Iturbide, M., Gutiérrez, J. M., Alves, L. M., Bedia, J., Cerezo-Mota, R., Cimadevilla, E., Cofiño, A. S., Di Luca, A., Faria, S. H., Gorodet- skaya, I. V ., Hauser, M., Herrera, S., Hennessy, K., Hewitt, H. T., Jones, R. G., Krakovska, S., Manzanas, R., Martínez-Castro, D., Narisma, G. T., Nurhati, I. S., Pinto, I., Seneviratne, S. I., van den Hurk, B., a...
-
[18]
Jeon, S., Paciorek, C. J., and Wehner, M. F.: Quantile-based bias correction and uncertainty quantification of extreme event attribution statements, Weather and Climate Extremes, 12, 24–32, https://doi.org/10.1016/j.wace.2016.02.001,
-
[19]
F.: Principal Component Analysis for Extremes and Application to U.S
Jiang, Y ., Cooley, D., and Wehner, M. F.: Principal Component Analysis for Extremes and Application to U.S. Precipitation, Journal of Climate, 33, 6441–6451, https://doi.org/10.1175/JCLI-D-19-0413.1,
-
[20]
Kass, R. E. and Raftery, A. E.: Bayes Factors, Journal of the American Statistical Association, 90, 773–795, https://doi.org/10.1080/01621459.1995.10476572,
-
[21]
and Hense, A.: A Bayesian Assessment of Climate Change Using Multimodel Ensembles
Min, S.-K. and Hense, A.: A Bayesian Assessment of Climate Change Using Multimodel Ensembles. Part I: Global Mean Surface Tempera- ture, Journal of Climate, 19, 3237 – 3256, https://doi.org/10.1175/JCLI3784.1,
-
[22]
Otto, F. E. L., Barnes, C., Philip, S., Kew, S., van Oldenborgh, G. J., and Vautard, R.: Formally combining different lines of evidence in extreme-event attribution, Advances in Statistical Climatology, Meteorology and Oceanography, 10, 159–171, https://doi.org/10.5194/ascmo-10-159-2024,
-
[23]
Paciorek, C. J., Stone, D. A., and Wehner, M. F.: Quantifying statistical uncertainty in the attribution of human influence on severe weather, Weather and Climate Extremes, 20, 69–80, https://doi.org/10.1016/j.wace.2018.01.002,
-
[24]
Paule, R. C. and Mandel, J.: Consensus Values and Weighting Factors, Journal of Research of the National Bureau of Standards (1977), 87, 377–385, https://doi.org/10.6028/jres.087.022,
-
[25]
Perkins-Kirkpatrick, S., Alexander, L., King, A., Kew, S., Philip, S., Barnes, C., Maraun, D., Stuart-Smith, R., Jézéquel, A., Bevacqua, E., Burgess, S., Fischer, E., Hegerl, G., Kimutai, J., Koren, G., Lawal, K., Min, S.-K., New, M., Odoulami, R., Patricola, C., Pinto, I., Ribes, A., Shaw, T., Thiery, W., Trewin, B., Vautard, R., Wehner, M., and Zscheisc...
-
[26]
Philip, S., Kew, S., van Oldenborgh, G. J., Otto, F., Vautard, R., van der Wiel, K., King, A., Lott, F., Arrighi, J., Singh, R., and van Aalst, M.: A protocol for probabilistic extreme event attribution analyses, Advances in Statistical Climatology, Meteorology and Oceanography, 6, 177–203, https://doi.org/10.5194/ascmo-6-177-2020,
-
[27]
Seneviratne, S., Zhang, X., Adnan, M., Badi, W., Dereczynski, C., Di Luca, A., Ghosh, S., Iskandar, I., Kossin, J., Lewis, S., Otto, F., Pinto, I., Satoh, M., Vicente-Serrano, S., Wehner, M., and Zhou, B.: Weather and Climate Extreme Events in a Changing Climate, in: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the S...
-
[28]
Seong, M.-G., Min, S.-K., and Zhang, X.: A Bayesian Attribution Analysis of Extreme Temperature Changes at Global and Regional Scales, Journal of Climate, 35, 8189 – 8203, https://doi.org/10.1175/JCLI-D-22-0104.1,
-
[29]
Sippel, S., Mitchell, D., Black, M. T., Dittus, A. J., Harrington, L., Schaller, N., and Otto, F. E.: Combining large model en- sembles with extreme value statistics to improve attribution statements of rare events, Weather and Climate Extremes, 9, 25–35, https://doi.org/https://doi.org/10.1016/j.wace.2015.06.004, the World Climate Research Program Grand ...
-
[30]
Smith, R. L., Tawn, J. A., and Coles, S. G.: Markov chain models for threshold exceedances, Biometrika, 84, 249–268, https://doi.org/10.1093/biomet/84.2.249,
-
[31]
Stott, P. A., Stone, D. A., and Allen, M. R.: Human contribution to the European heatwave of 2003, Nature, 432, 610–614, https://doi.org/10.1038/nature03089,
-
[32]
Szemkus, S. and Friederichs, P.: Spatial patterns and indices for heat waves and droughts over Europe using a decomposition of extremal dependency, Advances in Statistical Climatology, Meteorology and Oceanography, 10, 29–49, https://doi.org/10.5194/ascmo-10-29-2024,
-
[33]
Wehner, M., Stone, D., Krishnan, H., AchutaRao, K., and Castillo, F.: The Deadly Combination of Heat and Humidity in India and Pakistan in Summer 2015, Bulletin of the American Meteorological Society, 97, S81 – S86, https://doi.org/10.1175/BAMS-D-16-0145.1,
-
[34]
Xoplaki, E., Ellsäßer, F., Grieger, J., Nissen, K. M., Pinto, J. G., Augenstein, M., Chen, T.-C., Feldmann, H., Friederichs, P., Gliksman, D., Goulier, L., Haustein, K., Heinke, J., Jach, L., Knutzen, F., Kollet, S., Luterbacher, J., Luther, N., Mohr, S., Mudersbach, C., Müller, C., Rousi, E., Simon, F., Suarez-Gutierrez, L., Szemkus, S., Vallejo-Bernal, ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.