Earth ’ s atmosphere protects the biosphere from nearby supernovae

Geological evidence indicates that a supernova within 100 parsecs of Earth occurs around once per million years. Such nearby supernovas can produce an intense gamma-ray burst and a 100-fold increase of cosmic rays, lasting several centuries. We ﬁ nd that the effect of a short burst of gamma rays is small since they are strongly attenuated before reaching the lower stratosphere. Intense cosmic radiation affects stratospheric ozone but, due to compensating effects in catalytic chemical cycles, ozone depletion is moderate and comparable to that from current anthropogenic emissions. This also holds for the low-oxygen atmosphere during early evolution of terrestrial life. We estimate the increase in aerosol and clouds from a 100-fold increase of cosmic rays exerts a radiative forcing comparable in magnitude but opposite in sign to current anthropogenic climate forcing. We conclude that Earth ’ s atmosphere is effective at shielding the biosphere from nearby supernovae. Supernovae (SN) are exploding, massive stars caused by gravitational collapse during their ﬁ nal stage of evolution. They are the main source of galacticcosmicrays(GCRs;ionisingradiation,mainlyprotons)andcanalso generateintensegamma-raybursts.Detectionof 60 Fe (half-life 2.60 My) in oceansediments 1,2 ,andthepresenceofalocalgas-anddust-depletedbubble in the Milky Way, imply that several nearby core-collapse SN within about 100 parsec(pc,~300lightyears)ofEarthoccurredinthelast10 My 3 .Nearby SN were postulated to have the potential to cause mass extinctions 4,5 due to the

https://doi.org/10.1038/s43247-024-01490-9 Earth's atmosphere protects the biosphere from nearby supernovae Check for updates Theodoros Christoudias 1 , Jasper Kirkby 2,3 , Dominik Stolzenburg 4 , Andrea Pozzer 1,5 , Eva Sommer 6 , Guy P. Brasseur 7 , Markku Kulmala 8,9,10,11 & Jos Lelieveld 1,5 Geological evidence indicates that a supernova within 100 parsecs of Earth occurs around once per million years.Such nearby supernovas can produce an intense gamma-ray burst and a 100-fold increase of cosmic rays, lasting several centuries.We find that the effect of a short burst of gamma rays is small since they are strongly attenuated before reaching the lower stratosphere.Intense cosmic radiation affects stratospheric ozone but, due to compensating effects in catalytic chemical cycles, ozone depletion is moderate and comparable to that from current anthropogenic emissions.This also holds for the low-oxygen atmosphere during early evolution of terrestrial life.We estimate the increase in aerosol and clouds from a 100-fold increase of cosmic rays exerts a radiative forcing comparable in magnitude but opposite in sign to current anthropogenic climate forcing.We conclude that Earth's atmosphere is effective at shielding the biosphere from nearby supernovae.Supernovae (SN) are exploding, massive stars caused by gravitational collapse during their final stage of evolution.They are the main source of galactic cosmic rays (GCRs; ionising radiation, mainly protons) and can also generate intense gamma-ray bursts.Detection of 60 Fe (half-life 2.60 My) in ocean sediments 1,2 , and the presence of a local gas-and dust-depleted bubble in the Milky Way, imply that several nearby core-collapse SN within about 100 parsec (pc, ~300 light years) of Earth occurred in the last 10 My 3 .Nearby SN were postulated to have the potential to cause mass extinctions 4,5 due to the depletion of stratospheric ozone 6 and consequent high ultraviolet (UV) exposure, and from increased number of aerosol particles and cloud cover, resulting in climate cooling (Fig. 1).
Ruderman 7 was the first to suggest that increased cosmic radiation from a nearby SN may deplete Earth's atmospheric ozone.A number of studies 4,[8][9][10][11] have corroborated this conclusion.However, Crutzen and Brühl, and Gehrels et al. 12,13 , disputed catastrophic impacts, using two-dimensional chemistry-transport atmospheric models to infer that ozone depletion was significantly less severe than previously thought.They concluded that a supernova must occur within 8 pc to have a catastrophic effect on the biosphere.A signature of such a close supernova, however, has not been found in the geological record.Thomas 14 concluded that supernovainduced stratospheric ozone loss and consequent increased UV exposure can contribute to changes in species abundances, but probably not to mass extinctions.Nonetheless, Fields et al. 5 suggested that a nearby (~20 pc) supernova could intensify ionising radiation to such an extent that species extinctions can be expected.Furthermore, a connection between GCRs and cloud formation has been suggested to account for palaeoclimatic observations of solar-climate variability 15,16 .Svensmark et al. 17 argued that if a near-Earth supernova were to occur, the increase in ionising radiation would strongly increase cloud condensation nuclei (CCN) concentrations.Svensmark 18 proposed that SN strongly influences climate and, in turn, ocean circulation and marine biodiversity during the Phanerozoic.

Atmospheric chemical composition
Although low-dimensional atmospheric models have been used in previous studies, Earth system models (ESM), capturing the complex atmospheric circulation dynamics, chemistry, and process feedbacks, are needed to simulate stratospheric ozone loss in response to elevated ionisation 16 , leading to ion-induced nucleation and particle growth to CCN 19,20 .Here we use the ESM with Atmospheric Chemistry (EMAC) model 21 with ioninduced aerosol processes based on CERN CLOUD experiments to study the impacts of nearby SN on Earth's atmosphere and climate.We assume a representative nearby SN with GCR ionisation rates in the atmosphere that are 100 times present levels 4,7 .Though the energy spectrum of the incident GCRs from a nearby SN may change with time and can be modulated by the heliospheric magnetic field (solar wind) and the geomagnetic field, these influences on the GCR flux are relatively minor relative to that of the SN and galactic magnetic field alignment.The factor 100 represents the peak GCR flux estimated by Melott et al. 22 for a SN at 50 pc and with a turbulent galactic magnetic field.Our model incorporates comprehensive atmospheric chemistry, propagation of GCR ionisation, ion-induced nucleation of new particles and their growth to CCN, and aerosol-cloud interactions, between Earth's surface and 80 km altitude (see Methods for details).
Ionising radiation dissociates stratospheric nitrogen and oxygen molecules into ground-state and excited N and O atoms. .The NO 2 may then be photolysed, NO 2 + hν(<420 nm) → NO + O( 3 P), which is rapidly followed by O( 3 P) + O 2 + M → O 3 + M, thereby regenerating NO and causing no net ozone loss per cycle.However, the NO 2 may alternatively interact with an O atom, NO 2 + O → NO + O 2 , which represents a catalytic cycle that regenerates the NO radical and destroys two ozone molecules 26,27 .The oxidation of ground-state N( 4 S) is slow compared with the reaction N( 4 S) + NO → N 2 + O, which limits the amount of NO that can accumulate and, in turn, ozone depletion 12 .In fact, a large fraction of the odd nitrogen (N + NO + NO 2 + 2N 2 O 5 + HNO 2 + HNO 3 + HNO 4 ) is present as compounds other than chemically-active NO x (NO + NO 2 )-in particular HNO 3 -regardless of the initial formation in the atmosphere as NO 12 .
Ionising radiation also produces hydroxyl radicals (OH), from the reaction O( 1 D) + H 2 O → 2OH.The yield is two OH molecules per ion pair 16,28 .In the lower and middle atmosphere, OH and HO 2 radicals rapidly cycle, establishing the HO x family and, in the absence of NO, causing net catalytic destruction of ozone 29 .The reaction cycle is OH +O 3 → HO 2 + O 2 followed by HO 2 + O 3 → OH + O 2 + O 2 or HO 2 + O → OH + O 2 .However, in the presence of NO, the hydroperoxyl radical, HO 2 , may regenerate OH via HO 2 + NO → NO 2 + OH.The NO 2 can then photolyse to produce an oxygen atom which then forms O 3 , resulting in a null cycle for ozone removal 30 .GCRs therefore produce reactive radicals that contribute to both the destruction and the formation of ozone, leading to a relatively smaller net change.
The overall relative change of NO x and HO x throughout the atmosphere is shown in Fig. 2a and b, respectively.While NO x increases throughout the troposphere and stratosphere, HO x decreases in the stratosphere and moderately increases in the troposphere.Our model simulations reveal substantially less ozone loss by cosmic rays from a nearby SN than previous estimates 31 ; Fig. 2c shows the relative change of ozone versus latitude and altitude, and Fig. 3a shows the change in ozone column thickness versus time-of-year and latitude.
Even though the increase of GCR flux incident on Earth is isotropic, the impact is largely localised over the polar regions due to the weak shielding (magnetic deflection) by the geomagnetic field, together with atmospheric circulations that dynamically isolate high-from middle latitudes in the lower stratosphere.The maximum ozone depletion over the poles is less than the present-day anthropogenic ozone hole over Antarctica, which amounts to an ozone column loss of 60-70% 32 .On the other hand, there is an increase of ozone in the troposphere, but it is well within the levels resulting from recent anthropogenic pollution 33 .
We find that the maximum mean stratospheric O 3 depletion, occurring at the end of the boreal summer, is of the order of 10 11 kg, equivalent to around 10% globally, which is close to the observed global ozone depletion due to anthropogenic halogen compounds of about 5-6% 32 .Although significant, it is unlikely that such ozone changes would have a major impact on the biosphere, especially because most of the ozone loss is found to occur at high latitudes.
Low oxygen atmosphere.We have also studied the impact of nearby SN on a low-oxygen atmosphere, to assess their possible effect on the evolution of life.We considered atmospheric conditions that prevailed during the Cambrian between 540 and 485 million years ago.During this period, the diversity of life forms sharply increased, and plants and animals began to populate the land, abandoning the UV protection of the oceans 34 .The emerging vegetation used photosynthesis to assimilate CO 2 and release O 2 , raising its atmospheric fraction from about 2% to 10% and facilitating the evolution of terrestrial animals 35 .
We simulated a 2% oxygen atmosphere since this would likely represent conditions where the emerging biosphere on land would still be particularly sensitive to ozone depletion.We have conducted a sensitivity test with reduced atmospheric oxygen levels to be directly comparable to the nominal simulation with SN-enhanced GCR.This does not purport to be a full paleoclimate assessment, but reducing the oxygen content of the atmosphere to a few percent allows us to quantify the nearby SN impacts on the ozone layer.In the absence of SN, the global average ozone column thickness under 2% O 2 conditions was about 100 DU (black contours in Fig. 3b), i.e.only a factor of three less than today's atmospheric average of roughly 300 DU, as observed by the WMO Global Ozone Observing System, despite the presence of only one-tenth the oxygen.A reduced amount of stratospheric ozone allows deeper penetration of short-wave UV radiation that photodissociates oxygen (hv < 242 nm).Therefore, O 3 formation relocates to lower altitudes with higher air and oxygen density, causing a vertical redistribution of ozone, which moderates ozone reduction in the vertical column.The relatively thinner 100 DU ozone layer was already able to absorb sufficient UV to allow evolution of life on land.The effect on this atmosphere of a  Bottom: relative increase factor of (d) bisulfate ions (HSO À 4 ), (e) new particle formation rates (>1.7 nm), and (f) total aerosol concentrations, due to a 100-fold increase in GCR intensity from a nearby supernova.100-fold increase in GCR ionisation from a nearby supernova is shown in Fig. 3b).Ozone loss is about 10-25% at mid-latitudes and an order of magnitude lower in the tropics (Fig. 3).At high latitudes, the ozone loss brings the column thickness to similar values as existed in the tropics at that time.At minimum ozone levels near the poles during Antarctic spring, an SN perturbation may even have increased the column ozone, since the short-wave UV solar radiation could penetrate to low altitudes where NO x -catalysed photochemistry enhances rather than decreases the formation of O 3 .We conclude that these changes of atmospheric ozone are unlikely to have had a major impact on the emerging biosphere on land during the Cambrian.
We also considered the effects of an intense γ-ray burst from an SN at 100 pc 31 .We assume a γ-ray intensity of 10 kJ m −2 , lasting for a few seconds.This is comparable to a brief factor-ten increase of the solar electromagnetic energy incident on the Earth.Despite receiving this huge electromagnetic pulse at short wavelengths, the planetary impact is remarkably small due to the protection of the atmosphere, which is around 28 radiation lengths in depth.The γ-ray burst is efficiently absorbed by ionisation of stratospheric air above 10 hPa (~30 km altitude), where the incident gamma rays are efficiently attenuated by two orders of magnitude.The resulting NO produced at high altitudes is then dispersed by the global stratospheric (Brewer-Dobson) circulation.Our simulations show that the atmospheric NO x content returns to the unperturbed level within a few weeks after the burst.

Aerosols, clouds, and climate
We have also simulated the effects of a nearby supernova on CCN and clouds due to increased ion-induced nucleation from a 100-fold increase of GCR intensity (Fig. 4a and b).The global mean increase of the number concentration of CCN (0.2% supersaturation) at cloud base level is about 10-20%.However, in regions where the unperturbed CCN concentrations are low, such as the Pacific and Southern Oceans, the increase of CCN is larger and can reach up to 100%.These changes, while climatically relevant, are comparable to the contrast between the pristine pre-industrial atmosphere and the polluted present-day atmosphere 36 .
The muted response of CCN to the intense ionisation from a nearby SN results both from the shorter than present-day ion lifetimes and an effective buffering of new particles to reach sizes above around 50 nm, where they constitute CCN.Nucleation and growth compete for the same condensable vapours so an increased nucleation rate implies slower growth rates for the new particles and a higher susceptibility to be lost by coagulation with preexisting aerosol particles.The effect of the shorter ion lifetimes is as follows.
At high ionisation rates, the ion lifetime is dominated by ion-ion recombination.The ion pair production rate ip , where α = 1.6 × 10 −6 cm 3 s −1 throughout the troposphere is the ion-ion recombination coefficient 37 and n ip [cm −3 ] is the ion-pair concentration.The ion lifetime is τ . So, at 100 times present ionisation rates, ion lifetimes are a factor 10 shorter than in the present-day atmosphere.During the supernova irradiation, the ion lifetime is τ = 56 s at ground (Q = 200 i.p. cm −3 s −1 ), and 18 s at 10 km (Q = 4000 i.p. cm −3 s −1 ).These lifetimes are much shorter than the mean time for a primary ion to form a bisulfate (HSO À 4 ) ion, marking the preliminary step for ion-induced nucleation to take place.The collision rate of each negative ion with an H 2 SO 4 molecule is k where k −0 = 2 × 10 −9 cm 3 s −1 is the charged-neutral collision rate coefficient 38 and [H 2 SO 4 ] [cm −3 ] is the sulphuric acid vapour concentration.For example, at a representative sulphuric acid concentration of 10 6 cm −3 , the collision rate of each negative ion with a sulphuric acid molecule is k = 2 × 10 −3 s −1 , equivalent to a mean time of 500 s.This implies that only around 10% or less of the ions will form an HSO À 4 ion before they are neutralised.Moreover, a further 500 s is required to nucleate the first molecular cluster, H 2 SO 4 Á HSO À 4 , during which time the HSO À 4 ion undergoes neutralization by ion-ion recombination.Consequently, only a small fraction of the ions deposited by GCR from a nearby SN are available to drive ion-induced nucleation of sulphuric acid particles.
These SN effects are incorporated into our model and lead to relatively modest changes in HSO À 4 concentration, new particle (>1.7 nm) formation rates, and aerosol particle concentrations, as shown in Fig. 2d-f, respectively.The global tropospheric changes in bisulfate ion concentration (Fig. 2d) are small, despite 100-fold increasing GCR intensity and a tenfold increase in ion pair concentrations.The largest increases occur in the upper troposphere, where the ion-ion recombination rate is around a factor ten slower than at lower altitudes.The rate of formation of new particles is enhanced above the tropopause in high latitudes and in the tropical lower stratosphere, as well as close to the surface (Fig. 2f).However, the overall increase in nucleation-mode aerosol particles (Fig. 2e) from two orders of magnitude increase in GCR flux is relatively modest, at about a factor 2-5 in the troposphere and mid-stratosphere.This increases particles at CCN sizes by less than a factor of 2 (Fig. 4b).
To further investigate the effect of increased ionisation on particle growth rates, we performed highly detailed box-model simulations using the ion-UHMA model 39 .This showed growth rate enhancements for newlyformed particles of less than a factor 1.5, which is found for upper tropospheric conditions, despite a factor of 100 higher ionisation rates.This supports our EMAC model simulations and will not result in significantly increased survival rates for new particles to reach CCN sizes.
Our simulation of the outgoing shortwave solar radiation at the top of the atmosphere during a nearby SN is shown in Fig. 4c, and the fractional change in outgoing solar radiation relative to the unperturbed atmosphere is shown in Fig. 4d.The increased CCN concentrations at high latitudes enhance the cloud albedo and, consequently, exert a cooling effect.Figure 4e shows the resultant zonal average radiative forcing (net change in shortwave plus longwave energy flux at the top of the atmosphere).There is a small positive climate forcing in the tropics (due to raised cloud top heights) and a sustained negative forcing (cloud brightening) in most of the rest of the world.We estimate the global mean radiative forcing from a 100-fold increase in GCR intensity is −2.5 W/m 2 (a cooling).This is of opposite sign but comparable in magnitude to anthropogenic radiative forcing since the pre-industrial age 40 , though the timescales are very different.

Conclusions
Using an Earth system model with comprehensive atmospheric chemistry and with ion and aerosol nucleation processes based on the CERN CLOUD experiments, we have simulated the impacts on the atmosphere and climate of a nearby supernova that produces a 100-fold increase in GCR intensity.Our study does not consider the direct health risks to humans and animals resulting from exposure to elevated ionising radiation.The considered scenarios contain some simplifying, but justified assumptions, in particular regarding the low-oxygen atmosphere, which aims to represent the early Cambrian ozone layer but not the climate of the period.While the general result that the atmosphere and climate are not strongly perturbed by SN is robust, this work cannot serve as a model of specific SN impacts on the Earth system in their full complexity.The current mean annual dose equivalent from all-natural sources of ionising radiation is 2.4 mSv (mostly Rn), of which cosmic rays contribute 0.35 mSv 41 .A factor 100 increase in exposure from ground-level cosmic rays (35 mSv/y) is equivalent to a lifetime dose of around 2 Sv in the tropics and mid-latitudes.Higher increases in cosmic ray flux could occur in the unlikely case of an ordered interstellar magnetic field aligned along the line of sight to the supernova 22 .Although 2 Sv would be a very high acute exposure, there is little epidemiological data to indicate the effects of chronic (long-term) exposure reaching this level.For comparison, the 50%-mortality acute (short-term) dose is 1-10 Sv for mammals and birds, 10-100 Sv for reptiles and fish, and 10 3 -10 4 Sv for protozoa, bacteria, and insects 42 .A lifetime exposure that is comparable to the lethal acute dose may have some consequences for life on Earth and deserves further study.
We find the seasonal local maximum depletion of stratospheric ozone is about 30%, located over the polar regions.This is of similar magnitude to the present-day polar ozone loss caused by anthropogenic emissions.The average clear sky (no cloud) UV-index change near the equator and in middle latitudes, where the majority of life resides and has evolved, is less than two (equivalent to ~15% increase), with more pronounced differences in polar latitudes.However, the change in solar radiation reaching the surface, when considering the total atmospheric column albedo and extinction, is only a few percent, with the highest increase limited to the polar regions.Although such ozone depletion may exert UV stress on the biota at high latitudes, the effect on the rest of the biosphere is minor 43 .The same applies to the low-oxygen atmosphere that existed during the Cambrian explosion of terrestrial life on Earth.
Interestingly, in an atmosphere with ten times less oxygen than today, the ozone layer is estimated to be only three times thinner.We similarly find that γ-ray bursts from nearby SN do not cause substantial ozone loss or otherwise impact the biosphere.Concerning increases in atmospheric aerosol from ion-induced nucleation, we find relatively muted responses due to the limited availability of condensable vapours and to reduced ion lifetimes.We estimate the increase in CCN abundance exerts a global mean radiative forcing of −2.5 W/m 2 (a cooling), which is comparable in magnitude but of opposite sign to current anthropogenic climate forcing.
Overall, we find that nearby SN are unlikely to have caused mass extinctions on Earth.We conclude that our planet's atmosphere and geomagnetic field effectively shield the biosphere from the effects of nearby SN, which has allowed life to evolve on land over the last hundreds of million years.

Methods
The EMAC model is a numerical chemistry and climate simulation system that includes sub-models describing tropospheric and middle atmosphere processes and their interaction with oceans, land, vegetation, and human influences 21 .It uses the second version of the Modular Earth Submodel System (MESSy2) to link multi-institutional computer codes.The core atmospheric model is the 5th generation European Centre Hamburg general circulation model (ECHAM5, 44,45 ).For the present study, we applied EMAC (ECHAM5 version 5.3.02,MESSy version 2.55.0) in the T42L47 resolution 46 , i.e. with a spherical truncation of T42 (corresponding to a quadratic Gaussian grid of approximately 2.8 degrees in latitude and longitude) with 47 vertical hybrid pressure levels up to 0.01 hPa (~80 km altitude).
EMAC uses a modal representation of aerosols dynamics (GMXe) 47 that describes the aerosol size distribution as seven interacting log-normal distributions, of which four modes are soluble and three are insoluble.New particles are added directly to the nucleation mode.We employ the Mainz Chemistry mechanism developed by 48 , evaluated by 49 , and described previously for use in global climate-chemistry simulations 50,51 .The EMAC model reproduces the observed middle atmosphere trace gas, ozone, aerosol, and water cycles, including the characteristic "tape recorder" signal in stratospheric water vapour 52 .The aerosol representation in EMAC, including mineral dust particles, has been documented by ref. 53.The activation of aerosol particles as CCN and the formation of cloud droplets are described by ref. 54, and the activation of ice nuclei in mixed-phase and cirrus clouds follow the work of 55 and references therein.Our simulations employ the MSBM submodel 56,57 to account for the formation of polar stratospheric clouds detailed heterogeneous chemistry reactions of Ozone-reducing compounds 58 .
The applied model setup comprises the aerosol nucleation submodel detailed in ref. 20 that implements the parameterisation to calculate aerosol particle formation rates from CERN CLOUD measurements, including ioninduced nucleation (ref.19 Supplementary Materials Sec. 8).The model calculation of ion-pair production rate that mediates the ion-induced or ionmediated nucleation is also described in ref. 20.The model includes a GCR ionisation scheme 59 as a function of atmospheric depth, GCR modulation, and geomagnetic cut-off rigidity.Values between the tabulated points are calculated by linear interpolation, similar to 19 .The geomagnetic cut-off rigidity is calculated following the method of ref. 60, using the IONS EMAC submodel 20 , based on the International Geomagnetic Reference Field (IGRF) 61 for the year 2010.The heliospheric modulation and geomagnetic field coefficients used for the simulations can be found in Table 1.
For the supernova simulations, the IONS submodel ion-pair production is increased by a factor of 100 relative to the control simulations with present-day rates, throughout the atmosphere 22 .We simulated year-long periods under present-day (control) and supernova conditions (with 100fold increase in GCR), after a spin-up period of one year.Each GCR ion-pair produced in the atmosphere is associated with the production of 1.25 NO 23 and 2 OH molecules 16,28 .
Ionisation in the atmosphere from GCRs primarily results in charged oxygen, O À 2 , and nitrogen, N þ 2 , ion pairs.Our model then propagates the primary negative ions to vapours having progressively lower proton affinity.The subsequent ion-transfer reactions are complex, but tropospheric and stratospheric negative ion composition is generally dominated by NO À 3 and its various hydrates 62 .However, NO À 3 ions will transfer their charge upon collision with a sulphuric acid molecule to form HSO À 4 and its hydrates 63 .It is not known how HSO À 4 concentrations respond to an increasing primary ionisation.Here we use a parameterised scheme to approximate a possible moderation effect due to the enhanced ion-ion recombination.We model the ion cascade as a set of coupled ordinary differential equations describing the evolution of atmospheric ion composition, focusing on the major negative species in a confined air parcel.We formulated an initial value problem and derived an analytical expression for the steady-state solution.
We follow the reactive chain of ref. 64 for propagating negative ions.We used the charged-neutral collision rate coefficient of 2 × 10 9 cm 3 s −1 measured by ref. 38 as the reaction rate constant for the ion-neutral molecule collisions.When a negative ion collides with a neutral monomer of higher gas phase acidity (lower proton affinity) the negative charge is transferred with unit probability, i.e. a single collision is sufficient.The negative ions are neutralised at the rate where the positive ion number concentration and Q is the ion pair production rate, such that k ¼ ffiffiffiffiffiffi αQ p .The model includes the altitude dependence of the ion-ion recombination coefficient, α 37 .We have used the Tropospheric Ultraviolet and Visible Radiation Model 65 to calculate the difference in UV index during clear-sky (no cloud) conditions, taking into account the Ozone layer and aerosol optical thickness differences between the SN and nominal GCR simulations.
The prompt SN γ-ray spectrum is non-thermal and is approximated using 66 evenly spaced logarithmic energy bins in the range 0.001 MeV ≤ E ≤ 10 MeV, following the Band spectrum 13 .We scale the top of the atmosphere incident photon flux to the estimated γ-ray flux of 10 kJ/ m 2 .We note that a conservative estimate by Melott et al. 66 suggested that the probable nearest burst pointed at the Earth in the last billion years was ~2 kpc (perhaps 1 kpc).We interpolate the values of the photon attenuation coefficient from the Particle Data Group 67 to calculate mass lengths for the GRB spectrum for each of the 66 bins and estimate the ion-pair production at each atmospheric depth assuming 43.6 eV per conversion.We account for the atmospheric air mass at oblique incidence relative to that at zenith using the plane-parallel atmosphere approximation, scaling the atmospheric path length by the secant of the γ-ray zenith angle up to 85°.
For growth of nucleated particles to CCN sizes, we additionally used the University of Helsinki Multicomponent Aerosol model for neutral and charged particles (ion-UHMA).Ion-UHMA is a sectional (60 sections between 1.8 and 1000 nm) box model that simulates aerosol dynamics (condensation, coagulation, and deposition) and ion dynamics, ion-aerosol interaction, and ion-ion recombination during NPF, where the nucleation rates (neutral, positive and negative are treated as input and use the parametrisations from 19 .Sub-1.8 nm charged clusters are treated dynamically in the model, with an ion-pair production rate of 3 cm −3 s −1 or 300 cm −3 s −1 , and the aerosol population is modelled with three sub-populations (neutral and charged in both polarities) where the collision efficiencies also consider charge and dipole effects.More details can be found in ref. 39.We performed two sets of simulations (each with a high and low ion-pair production rate and across different sulphuric acid concentrations 0.1-5 × 10 7 molecules/cm 3 ), representative of the lower troposphere where organics participate in the growth, and of the upper troposphere with no background aerosol and only H 2 SO 4 .
Excited N( 2 D) atoms may react with O 2 to produce nitric oxide radicals, N( 2 D) + O 2 → NO + O.The yield is 1.25 NO molecule per ion pair 23 .Nitric oxide rapidly reacts with ozone to produce nitrogen dioxide, NO + O 3 → NO 2 + O 2 24,25

Fig. 1 |
Fig. 1 | Supernova impacts.Potential atmospheric and climate impacts of a nearby supernova.

Fig. 2 |
Fig. 2 | Relative annual mean change factor in atmospheric composition from a nearby supernova.Top: relative change of (a) NO x (NO + NO 2 ), (b) HO x (OH + HO 2 ), and (c) O 3 .The dashed line shows the mean tropopause height.

Fig. 3 |
Fig. 3 | Change in atmospheric ozone from a nearby supernova.Atmospheric ozone column percentage decrease (colour) due to a 100-fold increase in GCR intensity over nominal versus latitude and time of year.Black contours with numerical contour level magnitudes indicate the annual average O 3 column thickness in Dobson Units (1 Dobson Unit corresponds to 2.7 × 10 16 molecules cm −2 ) for nominal GCR intensities: (a) Present-day atmospheric composition, without ozone-holeinducing anthropogenic emissions.(b) Sensitivity test for 2% oxygen atmosphere, similar to conditions on Earth during the early Cambrian, a period renowned for its remarkable surge in biodiversity and evolutionary pace.The figure shows the important shielding effect of the geomagnetic field, which deflects a large fraction of the GCRs at low latitudes.

Fig. 4 |
Fig. 4 | Changes in cloud condensation nuclei (CCN) and radiative forcing from a nearby supernova.a Number concentration of CCN (0.2% supersaturation) at cloud base level (top of the boundary layer) after a 100-fold increase in GCR intensity and (b) fractional change in CCN relative to the present day.c Outgoing shortwave solar radiation at the top of the atmosphere with a 100-fold increase in GCR intensity and (d) fractional change in outgoing solar radiation relative to the present day, due to increased cloud albedo.e Zonal average radiative forcing (net change in shortwave plus longwave energy fluxes at the top of the atmosphere) due to a nearby supernova, and percentage change in solar radiation reaching the surface of the Earth.