Ukraine–Russia conflict: year 1 after detonation
After detonation, we explore the spatiotemporal evolution of the aerosol particles and the solar-radiative heating of the lower stratosphere (at 20 km) for the first 20 days (Fig. 1a–f), which sets the stage for the longer-term climate response. The spatial dispersion of the BC aerosol particles is illustrated in Fig. 1a–c, which shows the evolution of total aerosol optical depth (AOD) at 550 nm on days 1, 5, and 20. On day 1 (Fig. 1a), the BC particles remain confined to Eastern Europe. By day 5 (Fig. 1b), they have been lofted further into the stratosphere and been transported eastward by prevailing westerly winds across the Northern Hemisphere (NH). By day 20 (Fig. 1c), the particles have dispersed extensively throughout the NH, spreading towards the pole and equator. Corresponding stratospheric temperature changes at 20 km altitude are shown in panels d–f for days 10, 15 and 20. A modest hemispheric-mean stratospheric warming of ~2 °C in the Northern Hemisphere (NH) and <1 °C in the Southern Hemisphere (SH) is evident by day 10 (Fig. 1d). By day 15 (Fig. 1e) the NH mean stratospheric temperatures have increased by ~5 °C (SH ~1.8 °C), and by day 20 (Fig. 1f) it exceeds ~8 °C (SH ~ 1.6 °C), reflecting both the dispersion and solar-radiative heating of the BC particles. These patterns illustrate how even a moderate (5 Tg) BC atmospheric perturbation from a ‘limited’ nuclear conflict can rapidly disperse BC aerosols, globally blocking out sunlight and driving stratospheric heating.
Fig. 1: Evolution of aerosol optical depth and stratospheric temperature anomalies (ensemble mean) following 5 Tg BC aerosol release into the stratosphere (9–13 km) from 100 × 15 kt detonations at the Ukraine-Russia border (32.1° E–45° E, 46.4–52.3° N; ≈ 6 × 10⁵ km²).
The alternative text for this image may have been generated using AI.
a–c Total aerosol optical depth anomaly at 550 nm on days 1, 5 and 20 after detonation, illustrating aerosol spatial dispersion in the days after detonation. d–f Stratospheric temperature anomalies at 20 km on days 10, 15 and 20, illustrating a delayed but accelerating atmospheric thermal response. Hemisphere-mean anomalies are annotated beneath each temperature panel. Stippling in d–f marks grid points where the anomaly is not statistically significant (two-tailed Student’s t test, p ≥ 0.05).
Building on the hemispheric dispersion and solar-radiative heating shown in Fig. 1, we examine the vertical distribution of BC mass concentration and temperature changes at the start (Month 1, January) and end (Month 12, December) of the first year (Fig. 2a–d). Figure 2a, b shows the zonal-mean latitude–altitude cross-sections of BC mass concentration (ensemble mean) in first and last month of the year. In the first month (Fig. 2a), the mean BC concentration maxima occurs at ~15 km with concentrations > 0.1 µg m−3 in the NH. By contrast, by end of the year (Fig. 2b, month 12) elevated BC concentrations are observed throughout the stratospheric column across both hemispheres, illustrating the interhemispheric transport of BC. Figure 2c, d shows the corresponding zonal mean temperature anomalies (ensemble mean). In the first month (Fig. 2c), peak zonal mean warming of ~15 °C occurs near 20 km altitude. By year-end (Fig. 2d), the stratospheric temperature increase becomes more vertically extensive throughout the atmospheric column, with zonal mean temperature anomalies as high as ~40 °C. The evolution of the vertical and spatial dispersion of BC, and the associated thermal response highlights the critical role of aerosol-chemistry-dynamics feedbacks in modulating both the timing and magnitude of the Earth’s climate response.
Fig. 2: Zonal-mean vertical profiles of BC mass concentration and temperature response (ensemble mean) in month 1 (January) and 12 (December) after detonation (BC stratospheric emission occurring between Jan 1st and 7th).
The alternative text for this image may have been generated using AI.
Black carbon mass concentration (μg m−3) in Month 1 (January) (a) and Month 12 (December) (b). Corresponding zonal-mean temperature anomalies (°C) in Month 1 (January) (c) and Month 12 (December) (d).
Figure 3 illustrates both the seasonal cycle and annual-mean climate responses in year 1 post-detonation: panels a, b the change in downward solar radiation, c, d surface temperature anomalies, e, f precipitation anomalies, and d–h stratospheric wind perturbations. We also highlight anomalies across the land regions of Russia and USA, because they are the world’s two largest nuclear-armed states, and their climates serve as high-latitude and mid-latitude benchmarks for assessing how even a regionally confined conflict in Eastern Europe could impose severe climate stress on the principal nuclear powers.
Fig. 3: Seasonal and annual-mean surface and stratospheric circulation response in the first-year post detonation.
The alternative text for this image may have been generated using AI.
Surface solar radiation anomaly: a Monthly mean surface solar radiation anomalies (shading represents the ensemble variability) for the Northern Hemisphere (blue), Southern Hemisphere (red), Russia (green) and USA (purple). b Annual-mean surface solar flux anomaly. Surface temperature anomalies: c Monthly mean temperature anomalies d Annual-mean temperature anomaly. Precipitation anomaly: e Monthly rainfall anomalies f Annual mean precipitation percentage change. Stratospheric wind anomalies: g Time series of maximum zonal (blue), meridional (orange) and vertical (green) wind anomalies at 40 km. h Annual-mean spatial map of wind-speed anomalies (shading) showing strengthened subtropical jets and altered polar circulations. Stippling indicates grid points where the anomaly is not statistically significant (two-tailed Student’s t test, p ≥ 0.05). Together, these panels demonstrate how the solar-radiative properties of BC drive substantial changes to climate variables in the 1st year after nuclear detonation.
As the BC disperses across the atmosphere globally, it induces a significant negative anomaly in the incoming surface solar radiation in year 1. In the NH, the downward solar radiation anomaly is ≈ −20 W m−2 by April (Fig. 3a). As the BC spreads into the SH stratosphere, the solar flux anomaly is ≈ −15 W m−2 by end of year (SH summer). Mid-latitude continental zones in North America and Russia show anomalies between –10 and –50 W m−2 throughout the year, underscoring the vulnerability of key agricultural regions. This is reflected in annual-mean reductions in net primary productivity (NPP) of up to 0.2 kgC m−2 across large parts of the NH, particularly over regions of North America and Asia (Supplementary Section 2 and Fig. S1). The annual-mean downward solar radiation anomalies (Fig. 3b) are substantial, ≈−13.4 W m−2 in the NH and ≈−6.3 W m−2 in the SH. The spatial pattern of solar flux anomaly closely follows the changing stratospheric BC mass concentration through the year (Fig. 2), with the largest springtime reductions occurring over the United States, Russia, and much of Eurasia, regions responsible for a significant share of global crop production.
Surface temperature cooling develops rapidly during the first year after detonation (Fig. 3c, d). The seasonal cycle (Fig. 3c) shows that surface temperature anomalies in the NH develop rapidly within the first few months after detonation, reaching ≈−1.5 °C by the end of the first year. By contrast, the SH shows negligible temperature anomalies year-round, owing to the thermal buffering influence of the Southern Ocean. The annual-mean temperature response (Fig. 3d) confirms a hemispheric-average cooling of ≈ −1 °C in the NH compared to a slight warming of +0.01 °C in the SH. Strong regional cooling (>−4 °C) occurs over mid and high-latitude continental regions. Russia experiences the strongest anomaly (≈–5 °C by end of year), owed to its high-latitude location, minimal buffering from the ocean, and the low heat capacity of its continental land surface. The United States also experiences substantial cooling (≈−3 °C in spring), with ensemble members showing year-end anomalies as extreme as ≈−10 °C.
Significant shifts in global precipitation are observed in year 1 as a direct response to the strong NH surface cooling and associated changes in large-scale circulation (Figs. 3e, f and 4a, b). Annual-mean precipitation decreases by 20–40% across much of the NH mid-latitudes, with locally larger decreases (up to 80%) over densely populated and agricultural regions of Asia and West Africa (Fig. 3f), alongside land regions of the United States and Russia experiencing precipitation decreases of up to 20 mm month−1 in year 1 (Fig. 3e). These patterns are also consistent with the seasonal cycle shown in Fig. 4a, where India and Niger experience pronounced monsoon season decreases of ~ 40–100 mm month−1 in July–August. The spatial pattern of drying corresponds closely to changes in the meridional overturning circulation and vertical winds. Meridional Stream-function anomalies and vertical winds over longitudinal sectors (64°E–91°E and 16°W–13°E, respectively) (Supplementary Section 2 and Figs. S2 and S3), indicate a supressed assent over India and West Africa during the NH summer (JJA) resulting in dryer conditions relative to the control simulation. In contrast, the SH mid-latitudes experience substantial precipitation increases, especially over Southern Africa and Australia, where annual precipitation rises by up to 100%. The seasonal cycle shows that both Australia and Namibia exhibit consistent precipitation increases from February-April, ~20–60 mm month−1 more than the control simulation (Fig. 4b). Stream-function and vertical wind anomalies over longitudinal sectors (9°E–40°E and 112°E–154°E) (Supplementary Section 2, Figs. S4 and S5) confirm the presence of a strengthened anomalous counterclockwise cell and increased vertical winds which indicates enhanced low-level ascent, which explains the increases in precipitation. These hemispheric contrasts are also influenced by the southward displacement of the Intertropical Convergence Zone (ITCZ) by approximately 2–6° between March and June caused by the asymmetric NH cooling. The weakened, shifted ITCZ produces lower equatorial precipitation by ~0.5–1.0 mm day−1 in year 1 (Supplementary Section 2 and Fig. S6). Collectively, the strongly altered overturning circulation and displaced Intertropical tropical convergence zones explain the hemispheric dipole: widespread dryer conditions across the NH and enhanced rainfall over the SH subtropical regions.
Fig. 4: Regional Precipitation impacts in the first-year post detonation.
The alternative text for this image may have been generated using AI.
Simulated year-1 monthly precipitation (mm month-1) for a India and Niger (NH), and b Australia and Namibia (SH). Comparison of Control (blue) and 5 Tg BC UKRRUS (red) scenarios. Solid and dashed lines represent the ensemble means for specific countries (indicated in the legend), with shaded bands indicating the ensemble range.
The stratospheric wind dynamics changes driven by atmospheric temperature changes are prominently illustrated through wind anomalies at ~40 km altitude (Fig. 3g–h), with the column vertical profile of the zonal, meridional and vertical winds shown in Fig. 5a–c. The time series (Fig. 3g) highlights sustained zonal jet anomalies with velocities ranging between 50 and 200 m s−1, meridional wind anomalies between 20 and 70 m s−1, along with smaller but significant vertical wind anomalies between 0.02 and 0.1 m s−1, indicative of pronounced changes in stratospheric vertical transport and mixing processes, peaking notably in the spring and summer months. The annual-mean wind anomalies at 40 km (Fig. 3h) show strengthening of the subtropical and polar jets in both hemispheres, alongside perturbations to the vertical winds at high latitudes. In comparison, previous studies19 modelling a 5 Tg BC India–Pakistan conflict reported lower maximum stratospheric zonal wind anomalies of ~40 m s−1 during the NH-winter, whereas our modelled scenario shows peak zonal wind anomaly between 100 and 200 m s−1 during the same time period, reflecting differences in injection latitude, soot heating distribution, and structural inter-model differences.
Fig. 5: Zonal mean vertical profile of the zonal, meridional and vertical winds (ensemble mean) in the first-year post detonation.
The alternative text for this image may have been generated using AI.
a Zonal-wind anomaly, Δu (m s−1). b Meridional-wind anomaly, Δv (m s−1). c Vertical-wind anomaly, Δw (cm s−1).
To evaluate how the geographic origin of the BC alters the subsequent Earth system response, we compare two war scenarios, each with 5 Tg BC released into stratosphere: one originating from the Ukraine–Russia border region (“UKRRUS”) and the other from the widely studied India–Pakistan regional conflict scenario (“INDPAK”)6,7,9,16,17. Previous modelling studies of the 5 Tg INDPAK scenario consistently show declines in surface temperature and global mean precipitation6,7,9. While there are notable differences in model setups/configurations in past modelling studies—which includes date of detonation, horizontal and vertical resolution, and aerosol microphysics; our INDPAK simulations produce quantitatively similar results. Reported global mean surface cooling ranges from 0.5 to 2 °C6,7,9,20, with reductions in global mean precipitation during the first year of 0.2–0.5 mm day−1. Figure 6a shows the zonal mean AOD anomaly for days 8 and 20 for both scenarios. The UKRRUS aerosol dispersion (blue/green) remains more confined to mid and high latitudes, whereas in the INDPAK scenario (red/orange), there is more aerosol dispersion into the SH and less so towards the NH high latitudes. Figure 6b illustrates the corresponding annual-mean AOD anomaly, highlighting stronger aerosol loading across Northern Eurasia and the Arctic in the UKRRUS case relative to INDPAK. This spatial pattern reflects the more persistent poleward transport of soot in the higher-latitude UKRRUS injection compared to the tropical-latitude INDPAK source. Consistent with these spatial contrasts in aerosol loading, Figure S7 (Supplementary Section 2) shows that both scenarios produce a comparable global-mean cooling of approximately 0.5 °C, but the UKRRUS case has an earlier peak NH cooling (end of year 1) compared with the INDPAK scenario (end of year 2). These differences in aerosol dispersion translate into distinct spatial patterns of climate disruption (Fig. 6c–e). Consequently, surface temperatures are lower at northern latitudes for UKRRUS (up to –2 °C) than for the INDPAK scenario, but somewhat warmer at lower latitudes (≈ + 1 °C over parts of South Asia) (Fig. 6c). These responses follow the spatial distribution of the change in the surface solar fluxes (Fig. 6d). Solar dimming is more pronounced in mid-latitude continental regions for UKRRUS (–0.7 W m−2 hemispheric mean) but yields a net positive anomaly in the SH (+2.2 W m−2) because the INDPAK BC plume disperses more across the tropics and into the SH. Relative to INDPAK, the UKRRUS scenario drives a hemispheric asymmetry in the precipitation response (Fig. 6e), characterised by widespread suppression in the NH (up to 20 mm month-1) and localised enhancement over Southeast Asia (up to 50 mm month-1). These anomalies are primarily driven by a southward displacement of the ITCZ, which shifts up to 4° further south by month 6 in the UKRRUS scenario compared to INDPAK (Supplementary Section 2 and Fig. S6). Consistent with this shift, analysis of the meridional stream function and vertical winds confirm strengthened vertical ascent over Southeast Asia driving enhanced rainfall during the NH summer (Supplementary Section 2 and Fig. S8).
Fig. 6: Comparison between a Ukraine-Russia and an India-Pakistan Conflict to assess the influence of detonation location on climate.
The alternative text for this image may have been generated using AI.
a Zonal-mean aerosol optical depth (AOD) anomalies at days 8 (solid) and 20 (dashed) for the Ukraine–Russia (blue/green) and India–Pakistan (red/orange) scenarios. b Annual-mean AOD anomaly at 550 nm (UKRRUS – INDPAK), c Annual-mean surface temperature difference (UKRRUS – INDPAK), d Annual-mean change in downward solar radiation (UKRRUS – INDPAK). e Annual-mean change in precipitation (UKRRUS – INDPAK). Stippling in (b–d) marks grid points where the anomaly is not statistically significant (two-tailed Student’s t test, p ≥ 0.05).
Ukraine-Russia conflict: multi-year climate response
To evaluate the long-term climate impacts of the nuclear conflict, Fig. 7 illustrates the temporal evolution of key surface climate variables over the subsequent 10-year period. In each panel, solid lines denote the ensemble-mean anomalies, while the shaded envelopes indicate the ensemble spread for the Northern Hemisphere (blue), Southern Hemisphere (red), Russia (green) and USA (purple).
Fig. 7: Ten-year evolution of climate anomalies post detonation.
The alternative text for this image may have been generated using AI.
a Monthly surface air temperature anomalies (°C). b Downward solar radiation anomalies (W m−2). c Precipitation change (%). a–c Illustrate the multi-year climate disruption and recovery which lasts ~6 years. The bold line indicates the ensemble mean and shaded bands representing the 6-member ensemble spread, for the Northern Hemisphere (blue), Southern Hemisphere (red), Russia (green), and USA (purple).
Land regions of Russia experience the largest negative temperature anomalies (Fig. 7a), as low as ~–6 °C (ensemble member minimum). The NH mean surface temperature also falls sharply (approximately –1 °C), whereas the SH temperature response is negligible owing to the buffering effect of the Southern Ocean. Between years 2 and 4, surface temperatures gradually recover as atmospheric BC concentrations decrease, and by year 6, anomalies approach baseline conditions. The temperature response over the US is similar to Russia with the ensemble minimum as low as ~−4 °C.
Incoming surface solar radiation anomalies (Fig. 7b) closely tracks the temperature changes. In the first year, the NH experiences a maximum decrease of approximately –13 W m−2, with the US showing the largest decrease (≈–30 W m−2) and Russia around –20 W m−2. As the aerosol dispersed globally into the SH, the SH anomalies go as low as approximately −10 W m−2 between years 1 and 2. From years 1 through 8, solar flux anomalies globally recover and approach baseline conditions.
Precipitation anomalies (Fig. 7c) closely track the temporal evolution of surface temperature and surface solar radiation anomalies. During the first 2 years, NH mean precipitation declines by roughly 10–20%, equivalent to about −6 to −10 mm month−1, with particularly strong reductions over Russia (up to −25%, or ≈ −15 to −20 mm month−1). The USA experiences more modest decreases of ~5–10% in comparison to Russia, albeit with substantial ensemble spread. These suppressed precipitation rates persist through years 3–6 before gradually returning toward baseline conditions. Precipitation anomalies in the SH remain minimal, constrained to within ±5% over the entire simulation period.
Together, these multi-year time series demonstrate that a nuclear war in Eastern Europe produces rapid, hemisphere-wide changes to surface temperature, solar radiation, and precipitation. This is followed by a gradual climate recovery phase, with our model projecting a return to near-baseline conditions within ~6 years.
Near-term radioactive fallout—48 h post detonation
Near-term fallout within the first 48 h after detonation is estimated using the simplified empirical model of Glasstone and Dolan (1977)11. Due to the inherent complexity and sensitivity of fallout dynamics to environmental parameters such as wind speed, precipitation, and conflict locations, this approach provides a practical framework for estimating early radioactive contamination. This approach is derived from historical nuclear weapons test data and has been applied to nuclear war impact modelling studies in the recent past16, and enables us to estimate areas exposed to radioactive contamination within the first 48 h following detonation. The 5 Tg BC that is emitted into the upper troposphere-lower stratosphere could result from a combination of both groundbursts and airbursts. Our radiological fallout estimates (methodology documented in the Supplementary Section 3) are carried out for 100 simultaneous 15-kiloton (kt) surface nuclear detonations scenario, evenly distributed across approximately 613,000 km² region along the Ukraine–Russia border (46.4°–52.3°N, 32.1°–45.0°E) (Fig. 8a).
Fig. 8: Global fallout analysis post detonation at the Ukraine-Russia border.
The alternative text for this image may have been generated using AI.
a Location of the target region at the Ukraine-Russia border where 100 nuclear detonations (each 15 kt yield) are uniformly distributed. b Temporal evolution of cumulative global and hemispheric BC deposition (in Tg) over a 10-year period, highlighting the e-folding timescale (~3.5 years) for global BC removal. Shaded areas represent ensemble spread. c 10-year cumulative effective dose (mSv) from Cesium-137(Cs-137) and Strontium-90 (Sr-90) deposition. d Top 10 countries by mean cumulative effective dose (mSv) e, Top 10 countries by Collective radiation dose (person–Sv)
The mean surface winds over the region in the first month is ~ 5.3 m s−1 (19 km h−1) towards the northeast (~40°). Climatologically, radioactive debris from surface detonations would primarily be expected to disperse towards the northeast shortly after the explosions. By adapting the Glasstone–Dolan fallout model to these specific wind conditions, we determine that each 15-kt detonation generates plumes that extend downwind approximately 62.7 km for radiation doses exceeding 1 Sv, 26.1 km for doses exceeding 5 Sv, and approximately 17.9 km for the highest lethal threshold of 10 Sv (Supplementary Section 3 and Table S1) in 48 h.
The fallout plume width ranges from ~0.9 km for ≥10 Sv to ~3.2 km for ≥1 Sv. Collectively, the 100 detonations produce a cumulative contaminated area of approximately ~20,000 km² (in 48 h) at doses exceeding 1 Sv, ~3500 km² above 5 Sv, and ~1600 km² above 10 Sv. For perspective, the 48-h fallout zone which experience doses ≥5 Sv surpasses the size of the Chernobyl exclusion zone (~2600 km²)21, although the isotopic composition would differ, with weapon fallout dominated by short-lived fission products (e.g. I-131, Ba-140, Te-132) rather than the longer-lived Cs-137 and Sr-90 typical of Chernobyl22.
Considering the average population density of the region (~49 persons/km², with a total population of ~30.2 million), we estimate significant human exposure and health impacts. Approximately a million people would experience radiation doses exceeding 1 Sv, a threshold typically sufficient to induce acute radiation sickness. Within zones receiving ≥5 Sv, roughly 170,000 people would encounter severe radiation sickness. The most critically impacted zone, receiving ≥10 Sv, would affect approximately an estimated 80,000 people, facing virtually certain fatal outcomes without advanced medical care. For perspective, the regulatory upper limit for artificial public exposure is only 1 mSv per year23, highlighting the extreme magnitude of these doses.
Beyond immediate health implications, the long-term consequences would persistently affect the region. Areas receiving ≥1 Sv would require extended evacuation, prolonged exclusion, or intensive remediation measures. Residual soil contamination, disruption of agriculture, and extensive infrastructure damage would contribute to prolonged socioeconomic instability. This analysis underscores that even a limited, regional nuclear conflict involving relatively low-yield weapons can produce extensive, long-lasting humanitarian and environmental devastation24.
Long term radioactive fallout from global BC transport and deposition
To assess the long-term fallout arising from global BC transport and deposition, we consider the release of radionuclides during nuclear detonation that can adhere to BC particles. Radionuclide attachment to aerosol particles has been well documented, including in observations following the Fukushima nuclear accident12,13. The self-lofting of BC into the stratosphere, combined with nuclear war-induced precipitation reductions delays wet removal of BC aerosol, prolonging its residence time. These BC-driven atmospheric changes allow radioactive debris to spread broadly before deposition. The atmospheric circulation and gradual surface deposition of BC facilitates widespread global dispersal of these radionuclides. Long-lived radionuclides such as Cesium-137 (Cs-137; T₁/₂ ≈ 30.2 yr) and Strontium-90 (Sr-90; T₁/₂ ≈ 28.8 yr) are routinely monitored in long-term fallout assessments owing to their persistent environmental residence, propensity for biological uptake, and well-documented health impacts over decadal timescales25. Assuming the weapons used are Uranium-235 (U-235) fission bombs, we estimate the yields of Cs-137 and Sr-90 and their adherence to BC, with assumptions and parameters detailed in Supplementary Section 4 and Table S2. We recognise that the extent to which radionuclides from surface detonations reach the stratosphere depends on their proximity to large fires and associated convective plumes; therefore, these estimates are an idealised upper-bound case for the distribution of radionuclides with soot.
A decade after detonation, deposition patterns exhibit a hemispheric asymmetry: fallout is initially confined largely to the NH, followed by substantial cross-equatorial transport into the SH, and ultimately dispersing radioactive contamination across the globe (Fig. 8b–e). Figure 8b shows that approximately two-thirds of the injected BC is deposited from the atmosphere in the first 4 years, corresponding to an atmospheric e-folding time of ~3.5 years (red dashed line). By year 10 nearly all the 5 Tg of BC has deposited globally. The surface deposition is initially concentrated in the NH, which receives most of the fallout in the first 1–2 years, but gradual cross-equatorial transport results in ~40% of the BC being deposited in the SH by year 10.
Maps of long-lived radionuclide fallout (Supplementary Section 2 and Fig. S9a, b) for Cs-137 and Sr-90 show a widespread global deposition (methodology and assumptions are documented in Supplementary Section 4), with the highest surface radiation levels up to ~10 Bq m−2 in the NH by year 10. Both Cs-137 and Sr-90 show a very similar spatial pattern (as both isotopes are transported with the same BC aerosol particles). In contrast, large parts of the SH areas receive negligible fallout (<0.1 Bq m−2), but some areas close of 45S receive >0.4 Bq m−2. These results indicate that the NH mid-latitudes, including parts of Central Asia, Europe and the Middle East accumulate the bulk of the global radioactive fallout.
Despite the broad geographic extent of radioactive contamination, the cumulative dose map (Fig. 8c) shows that resulting radiation exposures remain very low. The 50-year cumulative effective dose from Cs137 + Sr90 deposition is at most ~0.9 mSv in the most affected regions (e.g. parts of Central Asia). Most of the NH land area receives ~ 0.1–0.3 mSv, and virtually all populated areas in the SH stay below 0.07 mSv. Such doses are far below the natural background (~2.4 mSv yr−1) and by themselves would pose little direct health risk26. However, Fig. 8d indicates discernible differences in the radiation dose at the country level. The highest national mean 50-year cumulative doses (≈0.25–0.48 mSv per person) occur in smaller countries situated under the primary BC deposition areas: for example, Tajikistan (0.48 mSv) and Bhutan (0.38 mSv) top the list, followed by countries in Europe and Central Asia. Larger countries such as Russia have much higher localised deposition near the detonation zone, but their country average dose is diluted by vast areas with minimal fallout.
Collective radiation dose metrics26 (Fig. 8e) underscore how population size modulates impact. China and India’s large population, combined with its subtropical latitude, leads to a high collective radiation dose (China ~400,000 person–Sv and India ~200,000 person-Sv), followed by countries in Central Asia, Africa and Europe. Figure 7b–e demonstrate that BC-driven atmospheric transport can distribute radioactive fallout across hemispheres, leading to low-level contamination and associated risks far from the conflict region.
For perspective, the doses estimated here are several orders of magnitude smaller than those produced by atmospheric nuclear weapons testing prior to the 1963 Partial Test Ban Treaty. By that time, approximately 189 Mt of fission yield had been released into the atmosphere, roughly ten times the total fission yield assumed in our scenario, resulting in global mean annual effective doses peaking near 0.1–0.15 mSv in just 1 year (1963)27. In contrast, our modelled Cs-137 + Sr-90 fallout yields peak regional cumulative doses of ~0.8 mSv over a period of 50 years, underscoring how much smaller the long-term radiological burden would be under our simulated scenario.