The smoke product is off.
Satellite and kriged products can miss diffuse, aged smoke, especially in the East where plumes often sit aloft.
Look for: disagreement between smoke products and ground monitors that varies by region.Use economics to trace each step between a fire burning in Quebec or Manitoba and a death in a U.S. county—then match each step to evidence that public data can actually provide.
A smoky day in a U.S. county can come from U.S., Mexican, or Canadian fires. Exposure and attribution are different measurements.
Public mortality data are county-month; fire and wind data are daily and gridded, so every source is aggregated to one panel.
The per-µg/m³ effect comes from the research design; the total death count also relies on an attribution model.
The project’s umbrella question is:
How many deaths in the United States does smoke from Canadian wildfires cause, and what are those deaths worth?
The answer is not one jump from a fire to a death. It is a chain of physical and economic links, and each link is measured by a different dataset.
The 2023 Canadian fire season sent heavy smoke into the eastern United States, including the June 7–8 episode over New York. Published estimates of the U.S. death toll exist, but they are health-impact assessments: they multiply modeled smoke exposure by a dose–response function estimated in other populations (Zhang et al. 2025; Pye et al. 2026). No study yet estimates U.S. mortality from Canadian smoke using variation that is specific to Canadian fires. That gap is the project’s opening.
Four different quantities are often called “the death toll from smoke.” Name the one you mean before comparing studies.
Zhang et al. report about 4,100 acute U.S. deaths from 2023 Canadian smoke (Zhang et al. 2025).
Identified by a research design, not assumed from another setting.
Needs a model of how much of each county’s smoke came from Canada.
Converts deaths into an externality measured in dollars.
| Comparison | What differs | What it means | Can the project match it? |
|---|---|---|---|
| Acute vs. chronic | Short-run deaths vs. deaths from long-term exposure | Zhang et al.’s 4,100 is acute; their ~33,000 is chronic, a different estimand | Acute, yes; a county-month design identifies short-run mortality |
| Same-day vs. monthly | Lag 0 vs. deaths within a month (plus later months) | A monthly slope includes within-month lags | Approximately, with the lag difference stated |
| Thresholded vs. all smoke | Only “Canada smoke days” vs. all Canadian smoke | A narrower day filter mechanically lowers the count | Yes, by rebuilding their day filter on our grid |
| National vs. sample | Whole U.S. vs. counties large enough for stable counts | Excluded small counties must be added back | Yes, through state-remainder units |
Interpretation rule: a number above or below 4,100 is not evidence that the health-impact assessment is wrong until lag window, exposure definition, and geography are matched. The project builds that comparison one change at a time.
A naive regression of deaths on smoke is biased because smoke does not arrive at random.
Deathsct depend on smokect but smoke co-moves with heat, drought, stagnant air, and the weather patterns that also grow fires
Heat alone kills, so a smoke–death correlation can overstate the effect of smoke. Measurement error in smoke products and people avoiding smoke push the other way. The design must supply variation in smoke that is unrelated to local conditions.
Causal studies of U.S. smoke find that smoke raises mortality and emergency visits, often nonlinearly (Miller et al. 2024; Heft-Neal et al. 2023). National panels link smoke PM2.5 to mortality and project large future burdens (Qiu et al. 2025). Smoke also lowers earnings and test scores and changes behavior (Borgschulte et al. 2024; Wen and Burke 2022; Burke et al. 2022). These studies mostly rely on smoke from U.S. fires. None isolates smoke from Canada.
The 2023 estimates apply a pooled multi-country wildfire relative risk for acute deaths (Chen et al. 2021) and a meta-analytic relative risk for chronic deaths (Chen and Hoek 2020). If the project estimates a different slope for the eastern United States, at least four explanations compete:
Satellite and kriged products can miss diffuse, aged smoke, especially in the East where plumes often sit aloft.
Look for: disagreement between smoke products and ground monitors that varies by region.Smoke that crosses 1,000+ km ages chemically, so its toxicity per µg/m³ may differ from fresh smoke.
Look for: slopes that change with transport distance, holding the population fixed.Air purifiers, alerts, and habits may cut harm where smoke is common; eastern residents have had less practice.
Look for: smaller per-unit effects in counties with more past smoke, within the same region.Age, baseline health, air conditioning, and housing change how much outdoor smoke reaches lungs.
Look for: slopes that vary with age structure and air-conditioning prevalence (Janssen et al. 2002).Chinese pollution carried by dust storms raised mortality in South Korea (Jia and Ku 2019).
Project use: a design that combines transport events with source-country pollution.Trajectory-based exposure to Chinese pollution and Korean mortality, valued as an international externality (Heo et al. 2025).
Project use: the template for turning a spillover estimate into a dollar cost.Fire radiative power from Indonesian fires instruments for Singapore’s pollution and health outcomes (Sheldon and Sankaran 2017).
Project use: the closest precedent for a foreign-fire instrument.Smoke from Indonesian fires increased early-life mortality (Jayachandran 2009).
Project use: evidence that distant fire smoke has measurable mortality effects.Changes in wind direction move pollution across places for reasons unrelated to local health, which gives a credible instrument for pollution (Deryugina et al. 2019; Anderson 2020; Schlenker and Walker 2016). Fires upwind vs. downwind of a place provide similar variation (Rangel and Vogl 2019). This project combines both ideas: Canadian fire activity located upwind of a county, carried by the day’s winds.
Because such an index is built from a formula, some counties are always more exposed simply by geography. Recentering subtracts the exposure a county would expect from other fire seasons under the same winds, so that only the unusual part of each season’s fires remains (Borusyak and Hull 2023).
People adapt to environmental hazards. Heat–mortality relationships weakened as air conditioning spread and are smaller where heat is common (Barreca et al. 2016; Heutel et al. 2021). Pollution alerts change behavior (Neidell 2009), and households pay to defend themselves with air purifiers (Ito and Zhang 2020). These results make adaptation a plausible reason why eastern U.S. counties, which rarely faced heavy smoke before 2023, might suffer more per unit.
identified slope × smoke attributed to Canada × population, summed over counties and months
deaths × value of a statistical life, with a life-years version because smoke deaths skew old (Viscusi and Aldy 2003; Aldy and Viscusi 2008)
The table treats each paper as design → estimand → limitation → project use.
| Study | Context and design | What it credibly contributes | Transfer to this project |
|---|---|---|---|
| Zhang et al. (Zhang et al. 2025) | 2023 Canadian fires; transport modeling + borrowed dose–response | North American acute and chronic burden; 4,100 acute U.S. deaths | The comparator; match its lag, day filter, and geography before comparing |
| Pye et al. (Pye et al. 2026) | EPA air-quality modeling of 2023 smoke, PM2.5, ozone, air toxics | Official assessment including ozone and toxics | Reminds us smoke is a mixture; the project’s slope is per µg/m³ of the Canadian smoke mixture |
| Miller, Molitor & Zou (Miller et al. 2024) | U.S. smoke shocks; instrumental variables; elderly mortality | Causal, concave dose–response for smoke | A benchmark slope; concavity matters for monthly averages |
| Qiu et al. (Qiu et al. 2025) | U.S. county panel of smoke PM2.5 and mortality | National smoke–mortality function and future burden | Nearest national estimate; not identified from Canadian variation |
| Deryugina et al. (Deryugina et al. 2019) | Wind-direction instrument; Medicare mortality | Credible PM2.5–mortality effect from wind shifts | Design template for wind-driven variation |
| Sheldon & Sankaran (Sheldon and Sankaran 2017) | Indonesian fire radiative power → Singapore | Foreign fire intensity as an instrument | Precedent for using fire radiative power across a border |
| Jia & Ku (Jia and Ku 2019) | Dust-transport days × Chinese pollution → Korea | Transboundary mortality effect | Shows how transport timing isolates a foreign source |
| Borusyak & Hull (Borusyak and Hull 2023) | Econometric theory for formula instruments | Recentering removes bias from non-random exposure | Why the project subtracts expected exposure from other seasons |
| Heutel, Miller & Molitor (Heutel et al. 2021) | U.S. temperature–mortality by climate region | Harm per hot day is smaller where heat is common | Motivates testing whether smoke harm falls with past smoke |
| Childs et al. (Childs et al. 2022) | Machine-learning daily smoke PM2.5 for U.S. counties | Standard smoke measure for 2006–2023 | Robustness smoke product; the main product covers through 2025 |
Several studies examine the 2023 episode directly. Health-impact assessments estimate the North American and U.S. burden (Zhang et al. 2025; Pye et al. 2026). In New York City, asthma-related emergency visits rose during the June smoke (Chen et al. 2023). Among about 53,000 hemodialysis patients in 22 states, days under the 2023 smoke plume were associated with higher same-day mortality and hospitalization (Song et al. 2025). Satellite smoke climatologies show that smoke is a regular summer presence over U.S. cities: it is present on about one in five days that exceed the 24-hour PM2.5 standard, and on 10–20% of ozone exceedances in places such as the Northeast corridor (Kaulfus et al. 2017; Brey and Fischer 2016). The record 2023 Canadian season followed early snowmelt, drought, and unusual heat (Jain et al. 2024).
A new study must add something these do not:
Estimate the slope from Canadian fires upwind of each county, not from a borrowed function.
Rebuild the published estimate step by step so a difference is interpretable.
Cover 2006–2025, not only 2023, so identification does not rest on one episode.
Convert attributed deaths into dollars that no market or treaty currently prices.
| Source | Role | Useful measures | Appropriate join | Cannot establish |
|---|---|---|---|---|
| NASA FIRMS | Fire activity (instrument) | fire location, date, radiative power, satellite, collection | 2° source cell × local day | Smoke amount, smoke height, or where the smoke went |
| ERA5 daily statistics | Transport winds; weather regimes | u/v winds at 850 and 700 hPa; 500 hPa height; surface weather | grid → county population weights | Whether smoke actually arrived at the surface |
| CSU kriged smoke PM2.5 (Dryad) | Main smoke measure (2006–2025) | total PM2.5, background, HMS flag; smoke = flag × (total − background) | 15 km grid → county-day | The source country of the smoke |
| Childs et al. smoke PM2.5 | Robustness smoke measure (2006–2023) | daily smoke PM2.5 by county | county-day | 2024–2025; source attribution |
| NOAA HMS smoke polygons | Smoke extent and episode flags | analyst-drawn plume outlines and density | polygon → county-day | Surface concentration (plumes are often aloft in the East) |
| EPA EMBER | Canadian share of smoke (2023) | base, zero-fire, and zero-Canadian-fire model runs | 36 km grid → county-day | An independent check (it shares fire inputs with the instrument) |
| EPA AQS | Ground truth | daily monitor PM2.5 and ozone | monitor → county-day | Concentrations where there is no monitor |
| CDC WONDER | Deaths (county × year, county × month) | resident deaths by cause and age | county-month | Counts of 1–9 (suppressed); individual records |
| NCHS public-use deaths (NBER mirror) | Pre-sample monthly deaths, 1999–2004 | county × month deaths, counties ≥100,000 people | fixed county geography | Small counties or any year after 2004 |
| SEER populations | Denominators and weights | county population by age | county-year → month | Deaths or exposure |
Aggregate once, carefully. Fires are daily and gridded, winds are daily and gridded, deaths are county-month. Document every conversion—local solar day vs. UTC day, land vs. water area, population vs. area weights—because each one can shift the answer.
Question. How many U.S. deaths does smoke from Canadian wildfires cause, identified from Canadian-specific variation, and what are they worth?
Unit. County × month, 2006–2025, contiguous United States; excluding the COVID window (March 2020–December 2022) from estimation.
Instrument. Canadian fire radiative power carried toward each county by 850–700 hPa winds, recentered on 23 permuted fire seasons (2003–2025).
Outcome. All-cause deaths per 100,000 over the exposure month and the next three months, which nets out deaths merely moved forward by a few weeks.
Estimands. Keep three objects separate:
Interpretation. The slope is identified by the design. The total also depends on how much smoke is attributed to Canada, so the project reports several attribution models and an instrument-only version beside the headline.
| Hypothesis | Observable implication | Rival explanation | Evidence that would move the claim forward |
|---|---|---|---|
| H1 · Level | Unusual upwind Canadian fire raises smoke and then deaths | Fire seasons coincide with heat waves | Balance tests on weather that smoke cannot change; heat controls; placebo downwind fires |
| H2 · Scope | The matched estimate differs from the borrowed dose–response | Lag, day filter, and geography differ | A bridge that changes one scope element at a time |
| H3 · Adaptation | Per-unit harm is larger where past smoke was rare | Age, air conditioning, income, smoke chemistry | Within-region variation in past exposure, holding those rivals fixed |
| H4 · Persistence | Effects appear in later months too | Deaths moved forward (harvesting) | Cumulative effects over 0–3 months and beyond |
| H5 · Mechanism | Alerts and searches rise with Canadian smoke | Attention driven by news, not exposure | NWS alerts and search data tied to measured smoke |
When and where did Canadian smoke reach U.S. counties?
FIRMS + HMS + smoke PM2.5 + populationHow much of each county’s smoke came from Canada?
Add EMBER model runs and the transport indexHow many deaths did Canadian smoke cause per µg/m³?
Add the recentered instrument and outcome-blind gatesWhat is that externality worth, and to whom?
Add the value of a statistical life and age profilesMove upward only when the data support the next claim. Levels 1 and 2 are open to student projects now.
Student projects stay at Levels 1–2 and never use 2006–2026 mortality.
Rebuild a county-day series (smoke, fires, or weather) from the raw files and check it against the project panel.
Choose a smoke episode (for example June 2023 or July 2021) and document it: dates, source fires, plume maps, surface PM2.5.
Put HMS plumes, CSU smoke, Childs smoke, and AQS monitors side by side for that episode.
Map population-weighted smoke by county; test how the map changes with area vs. population weights.
State the claim, the unit, the comparison, the rivals, and the claim ceiling.
Change one construction choice (threshold, day definition, weights) and report what moves.
Document each Canadian-smoke episode reaching the East, 2006–2025, with dates, source provinces, and affected population.
Quantify where HMS, CSU, Childs, and monitors agree or disagree, and why the East is harder.
Compare EMBER’s Canadian share with a simple wind-direction share for June 2023.
Link NWS air-quality alerts and search interest in air purifiers to measured smoke.
Trace adaptation, smoke chemistry, or air conditioning from theory to a testable prediction.
Record exceptional-event filings, alert rules, and U.S.–Canada disputes over smoke, with dates.
For every result, record:
Report the identified per-µg/m³ effect first, then each attribution and scope choice that turns it into deaths and dollars—so a reader can see exactly what the number depends on.