The Resilience of Life to Astrophysical Events

Much attention has been given in the literature to the effects of astrophysical events on human and land-based life. However, little has been discussed on the resilience of life itself. Here we instead explore the statistics of events that completely sterilise an Earth-like planet with planet radii in the range 0.5–1.5R⊕ and temperatures of ∼300 K, eradicating all forms of life. We consider the relative likelihood of complete global sterilisation events from three astrophysical sources – supernovae, gamma-ray bursts, large asteroid impacts, and passing-by stars. To assess such probabilities we consider what cataclysmic event could lead to the annihilation of not just human life, but also extremophiles, through the boiling of all water in Earth’s oceans. Surprisingly we find that although human life is somewhat fragile to nearby events, the resilience of Ecdysozoa such as Milnesium tardigradum renders global sterilisation an unlikely event.

oceans, so the threat of atmospheric removal is contained within that of oceans boiling. We are therefore led to consider death due to heat or radiation.

Analysis
To raise the entire ocean temperature by T requires a deposit of = E M T o  wherein M o is the ocean mass and  the specific heat capacity of water. In order to increase the temperature of the entirety of the Earth's oceans we need to introduce a large amount of thermal energy. The total mass of water in the oceans is around 1.35 × 10 21 kg. The specific heat capacity of water is 4184 J kg −1 °C −1 so we require 5.6 × 10 24 J to raise the ocean temperature by 1 °C. Thus the tardigrade with a tolerance of up to 100 °C would survive until around 5.6 × 10 26 J were deposited into the ocean. This is a lower bound-such heat would not be evenly distributed, being it most likely to be deposited in the upper ocean. To provide a conservative bound, we seek to minimise the depth of the deepest ocean on any planet-a uniform distribution of oceans across the planet's surface. When ocean mass is small compared to that of the planet, the depth of the ocean is approximately Here ρ = 3M p /(4πR p 3 ) is the average planet density, α the fraction of the mass in ocean (on Earth M o ≈ 2.3 × 10 −4 M ⊕ ) and ρ w the density of water. Most of Earth's water is contained within rocks. To remain conservative, we consider only the mass of liquid water in the oceans. There may exist planets that are almost entirely water (α ≈ 1), however for life as we know it, we focus on Earth-like planets with oceans on the surface of a rocky planet. We give these explicitly as we will assume they are broadly unchanged between planets. For the Earth, this implies that there must be an ocean of at least 2.5 km in depth. This is far shallower than the deepest points, however it will constitute a lower bound. The intensity of gamma rays is attenuated by interaction with matter by a factor exp(−μD), wherein D is the depth and μ the attenuation coefficient. This varies based on the material and the frequency of the incident radiation. The tardigrade is capable of withstanding over 6000 Gy (enough to endow every kilogram of material with 6000 J of energy). If the ocean depth is greater than µ log (700)/ (the latter figure being the ratio of the energy deposit per unit mass required to boil water to that to kill a tardigrade) the water above will be boiling. In fact, if we consider a sufficient radiative flux to kill a tardigrade at depth D, the total energy deposited upon the planet is at least . If our oceans are more than a few metres deep, this exceeds the threshold energy at which the oceans would boil before radiation would kill the tardigrade. We therefore consider temperature increase as the primary source of sterilisation.
Large asteroids are the leading candidate for causing of the Cretaceous-Tertiary extinction which took place 65 million years ago, annihilating approximately 75% of species on the planet leaving the Chicxulub crater. This event devastated larger land animals. Of those with masses over 25 kg only a few ectothermic species survived. However, around 90% of bony fish species survived 14 and deep ocean creatures were largely unaffected by the event. We estimate an upper bound for the energy deposited by an asteroid of mass M a as being its free-fall energy from infinity to the surface of the planet p p is the escape velocity of the planet (v e ≈ 11.2 kms −1 for Earth), and v ∞ is given by Öpik's close encounter theory 15 . In order to raise the ocean's temperature by T, we require an asteroid of mass To annihilate tardigrades on Earth we require a mass over ∼1.7 × 10 18 kg. The largest observed asteroids in the Solar System are Vesta and Pallas, with masses of 2.7 × 10 20 kg and 2.2 × 10 20 kg respectively. There are only 17 other known asteroids of sufficient mass, and a few dwarf planets, the most massive ones being Eris and Pluto, whose masses are 1.7 × 10 22 kg and 1.3 × 10 22 kg respectively. We reiterate that our estimate of the required energy is conservative-it is likely that it would take a significantly more massive impact as ocean heat would only be a fraction of the total energy. Since we consider Earth-like planets, the order of magnitude of this mass does not vary greatly between the largest and smallest planets-if oceans constitute an equal fraction of mass this changes by less than an order of magnitude.
In Fig. 1 we present a model for the impact rate of asteroids as a function of the mass. This is based on the extrapolation of relation between crater diameter and impact rate 16 . The mass of the object is related to crater diameter following ref. 17, assuming asteroids with density ρ ~ 5 g cm −3 entering the atmosphere with incidence angle of 90° with respect to the normal. In reality, most asteroids have ρ ≈ 2 g cm −3 , and in the case of comets this value is even lower, being the value here assumed a conservative assumption. This is highly dependent on the asteroid distribution in our Solar System-we assume in the absence of other evidence, that other systems are similar, however this remains to be verified.
The bulk of the energy output of a supernova is carried by the shock wave. To give an upper bound on the range at which a supernova would remove all life from the planet, we assume that the shock wave carries all the energy released. The fraction of energy incident on a planet of radius R p at a distance d from the supernova is given by the fraction of the sphere of radius d which is covered by the planet's surface π π R d /(4 ) p 2 2 . To raise the temperature of a planet by T, we would require a supernova within a distance d given: For the Earth, this sterilisation distance is around 0.04 pc, far closer than the closest stars, Proxima Centauri. Were a supernova to occur at that distance, the ocean temperature would only rise by about 0.1 °C. Furthermore, although there is a dependence on the mass of the planet, this dependence is quite weak. Note that none of the stars in the Alpha Centauri system are large enough to go supernova. The nearest potential supernova is the IK Pegasi system, approximately 45 pc away, which is three orders of magnitude farther than the estimated sterilisation radius.
To assess the relative risk faced by any planet in our galaxy, we approximate the odds of a close enough supernova happening over a timespan of 10 9 years. We find the expected number of stars of sufficient mass within the sterilisation distance of a planet, and the odds that one of these stars goes supernova. The galactic habitable zones, regions wherein complex life may evolve, depends on the occurrence rate of supernovae. A detailed simulation-based study was done by Lineweaver et al. 18 . We evaluate the rate of SN at a position (r, z) (cylindrical coordinates) as follows: Because the nearest star is about 1.3 pc away, we can conclude that Earth is located in a fortunate position. Near the galactic centre the density of stars increases and the probability of a SN sterilising life is higher. Nonetheless, this rate is almost insignificant even close to the galactic core, reaching only around 1% of planets being sterilised. The calculation for GRBs is similar to the one for SNe, but now we have to assume that the energy is collimated into in jets. As before, we assume the most pessimistic scenario-the smallest jet angle with the largest energy. The energy is typically the same as that of a supernova, 10 44 J, but the jet angles can be as low as 2°, hence the energy incident on a target of radius r which lies entirely within the beam angle at a distance d is: Hence for an increase of 100 °C in the ocean temperatures, we would need a GRB within about 13.8 pc; again, this is an upper limit. The rate of occurrence of short GRBs per volume in the universe is 0.04 Gpc −3 yr −1 , and long GRBs is 0.15 Gpc −3 yr −1 6 . We will restrict these occurrences to within galactic discs of stars, therefore we divide this by the product of the comoving number density of galaxies (≈10 7 Gpc −3 ) and the volume occupied of the galactic disc (10 11 pc 3 ), we find that the rate is around 2 × 10 −10 pc −3 Gyr −1 , and hence the probability of a GRB within the a distance at which it would sterilise a planet, aligned such that one of the beams hit the planet is 3.2 × 10 −10 Gyr −1 . This number is extremely small and we can conclude that such event is unlikely. is the mass of an asteroid whose impact crater is equal to the size of the planet, causing complete destruction. Here we find that even with the most conservative bound, the likelihood of complete sterilisation is lower than around 10 −5 over the lifetime of the planet.
Planetary systems can also be disrupted by passing-by stars. The rate, R, of such encounters for a given stellar density, n * , is where σ is the cross section to disrupt the orbit of the planet, and v the velocity of the star. Typically, v ≈ 40 km s −1 .
In the neighbourhood of the Solar System ∼ − − ⁎ n 10 pc 3 3 19 . In the case of Earth, σ ∼ − 10 pc 9 221 . Therefore, the rate of interactions of stars and the Earth-Sun system would be ∼ × − − R 3 10 Gyr 8 1 . If we repeat this calculation for the average stellar density in the galaxy, which is ∼ . − 0 1 pc 3 , the rate would be increased to ∼10 −6 Gyr −1 . This number is an upper bound-we expect that only a fraction of systems that experience disruption would eject a planet-yet it is still extremely small and we can conclude that ejection by this mechanism is a very rare event.

Discussion
Our analysis has focussed on providing an absolute upper bound for the rate of complete sterilisation of an Earth-like planet during its evolution, by considering the required events that would lead to the death of the hardiest species on Earth. With such assumptions, we find that the probability is less than 10 −7 per billion years. The overall likelihood of complete sterilisation is small even for planets which could exist around dwarf stars for ten trillion years, the most likely time for life to find itself 22 . For asteroids the impact rate of deadly objects is  − − 10 Gyr 5 1 . Eventually the evolution of the host star will render a planet sterile, either through expansion to the point of oceans boiling, or through a post-collapse freezing. There is a third scenario, where life continues around geothermal vents on a rogue planet until capture by a new host system, or the source of heat is extinguished. The time scale for the former is conservatively bounded by the rate of disruption of planetary systems by wandering stars, being therefore ∼30 Myr in regions with density of stars comparable to the vicinity of the Solar System; the time scale for the latter depends on properties of the planet, but are of the order of billion years. Consequently, life could perdure on a rogue planet long enough for it to be recaptured. We do not fully understand the mechanisms by which life started, but once it exists on an Earth-like planet, the complete removal of all life (other than through evolution of the host star) is a very unlikely event.
In this analysis we have considered only direct effects of astrophysical events such as the direct boiling of oceans; however such events could serve as a trigger for a second change. An example of this is the runaway greenhouse effect 23 in which a less energetic event could evaporate a fraction of the ocean, which in turn leads to increased temperatures. Such an occurrence would have a lower threshold energy, as total evaporation of the oceans is no longer necessary; however the precise details of such an event are not known. As an upper bound, we can presume that this would have to increase global temperatures by a larger amount than seasonal variations and changes in solar output. In essence this would reduce our energy threshold by about two orders of magnitude, still leaving such an occurrence very rare.
Throughout this analysis we have made several assumptions, the primary being that life elsewhere will be similar to that found on Earth. To justify this assumption, we consider that when searches are conducted for exoplanetary life, the template used is that of life on Earth. We make the further assumption that life will evolve to adapt to the extreme environments of exoplanets as it has to those on Earth. Again, we justify this by the ubiquity of life across environmental conditions. Finally, we note that the type of life we expect to survive all but the most extreme of events is that which could survive elsewhere in our own Solar System. The history of Mars indicates that it had an atmosphere 24 that could have supported life, albeit in conditions that on Earth would be considered inhospitable, with both the Mars 2020 and ExoMars programs aiming for experimental verification. Organisms with similar tolerances to radiation and temperature such as tardigrades would be the only kind of life that could survive long-terms in such conditions, Figure 2. The expected number of supernovae within the sterilisation sphere of radius 0.04 pc per Gyr as a function of galactic position. The black circle indicates the position of the Solar System. Closer to the galactic centre the stellar density is higher, and thus the likelihood of encountering a nearby supernova increases. However, this density is only sufficient to give a total rate of around 0.01 expected events per billion years, and thus total sterilisation through supernovae is still an improbable event. and even then they would have to be significantly below the surface. The subsurface oceans that are posited to exist on Europa and Enceladus would have conditions similar to the deep oceans of Earth where taridgrades are found -volcanic vents providing heat in an environment devoid of light. The presence of extremophiles in the locations would be a significant step in narrowing the question of life starting on exoplanets.