Chapter 7: Assessing Habitability Beyond Earth

All known life on Earth inhabits environments that maintain conditions between certain extremes of temperature, chemical composition, energy availability, and so on (Chapter 6). Life may have emerged in similar environments elsewhere in the Solar System and beyond. The ongoing search for life elsewhere mainly focuses on those environments most likely to support life, now or in the past-that is, potentially habitable environments. Discussion of habitability is necessarily based on what we know about life on Earth, as it is our only example. This chapter gives an overview of the known and presumed requirements for life on Earth and discusses how these requirements can be used to assess the potential habitability of planetary bodies across the Solar System and beyond. We first consider the chemical requirements of life and potential feedback effects that the presence of life can have on habitable conditions, and then the planetary, stellar, and temporal requirements for habitability. We then review the state of knowledge on the potential habitability of bodies across the Solar System and exoplanets, with a particular focus on Mars, Venus, Europa, and Enceladus. While reviewing the case for the potential habitability of each body, we summarize the most prominent and impactful studies that have informed the perspective on where habitable environments are likely to be found.

Nonsynchronous rotation of icy moon ice shells: The thermal wind perspective

The ice shells of icy satellites have been hypothesized to undergo nonsynchronous rotation (NSR) under the influence of tidal torques and/or ocean currents. In this work, the author proposes that the thermal wind relationship can be combined with geostrophic turbulence theory to predict ocean stress onto the ice shell inside the tangent cylinder. High-resolution numerical simulations validate the prediction within a factor of 2. For the prediction to be valid, the rotation effect must dominate (Rossby number < 1), and the upper ocean should be stratified. The latter can be achieved with sufficiently large ice thickness variations [the threshold for Europa is O(100) m]. Using this framework, once the ice rheology, thickness variations and NSR rate are determined, one may be able to estimate the ocean overturn timescale and put constraints on the ocean vertical diffusivity or the heat flux originating from the silicate core.

Ice Transit and Performance Analysis for Cryorobotic Subglacial Access Missions on Earth and Europa

Ice-covered ocean worlds, such as the Jovian moon Europa, are some of the prime targets for planetary exploration due to their high astrobiological potential. While upcoming space exploration missions, such as the Europa Clipper and JUICE missions, will give us further insight into the local cryoenvironment, any conclusive life detection investigation requires the capability to penetrate and transit the icy shell and access the subglacial ocean directly. Developing robust, autonomous cryorobotic technology for such a mission constitutes an extremely demanding multistakeholder challenge and requires a concentrated interdisciplinary effort between engineers, geoscientists, and astrobiologists. An important tool with which to foster cross-disciplinary work at an early stage of mission preparation is the virtual testbed. In this article, we report on recent progress in the development of an ice transit and performance model for later integration in such a virtual testbed. We introduce a trajectory model that, for the first time, allows for the evaluation of mission-critical parameters, such as transit time and average/overall power supply. Our workflow is applied to selected, existing cryobot designs while taking into consideration different terrestrial, as well as extraterrestrial, deployment scenarios. Specific analyses presented in this study show the tradeoff minimum transit time and maximum efficiency of a cryobot and allow for quantification of different sources of uncertainty to cryobot's trajectory models.

Exploring the Interior of Europa with the Europa Clipper

The Galileo mission to Jupiter revealed that Europa is an ocean world. The Galileo magnetometer experiment in particular provided strong evidence for a salty subsurface ocean beneath the ice shell, likely in contact with the rocky core. Within the ice shell and ocean, a number of tectonic and geodynamic processes may operate today or have operated at some point in the past, including solid ice convection, diapirism, subsumption, and interstitial lake formation. The science objectives of the Europa Clipper mission include the characterization of Europa's interior; confirmation of the presence of a subsurface ocean; identification of constraints on the depth to this ocean, and on its salinity and thickness; and determination of processes of material exchange between the surface, ice shell, and ocean. Three broad categories of investigation are planned to interrogate different aspects of the subsurface structure and properties of the ice shell and ocean: magnetic induction, subsurface radar sounding, and tidal deformation. These investigations are supplemented by several auxiliary measurements. Alone, each of these investigations will reveal unique information. Together, the synergy between these investigations will expose the secrets of the Europan interior in unprecedented detail, an essential step in evaluating the habitability of this ocean world.

The Impact of Extreme Low Temperatures on Raman Spectra of Amino Acids Relevant for the Search for Life on Europa

Raman spectroscopy, an emerging technology for in situ space exploration, has been suggested for life detection for the Europa Lander Mission. However, obtaining spectra of samples from the europan icy shell requires measurements at temperatures down to -233°C, which will affect the Raman spectra of any potential biosignatures. In this study, we obtained Raman spectra of amino acids using a 785 nm Raman system at temperatures ranging down to -196°C, analogous to Europa's surface and near subsurface. Significant Raman band width narrowing and decreasing variance were observed at lower temperatures leading to higher-precision Raman measurements, which required higher spectral resolution that could be as high as 2 cm^-1 for full identification of amino acids. Such spectral resolution is much higher than the resolution of contemporary Raman instruments for planetary exploration and may be particularly problematic for miniaturized instruments. Shifting of Raman bands to both higher and lower frequencies by as much as ∼25 cm^-1 together with changes in the Raman band intensity were recorded. The emergence of new bands and diminishing of the original bands also occurred for some amino acids. A significantly increased fluorescence background was observed in spectra of fluorescent molecules (i.e., tryptophan). A link between the type of vibrational modes associated with Raman bands and the change in their Raman shift at extreme low temperatures was identified and described. This link offers an exciting new method of molecule identification solely based on the comparison of spectra collected at two different temperatures and could greatly improve the identification capabilities in Raman spectroscopy for a wide array of applications.

How does salinity shape ocean circulation and ice geometry on Enceladus and other icy satellites?

Of profound astrobiological interest, Enceladus appears to have a global saline subsurface ocean, indicating water-rock reaction at present or in the past, an important mechanism in the moon's potential habitability. Here, we investigate how salinity and the partition of heat production between the silicate core and the ice shell affect ocean dynamics and the associated heat transport-a key factor determining equilibrium ice shell geometry. Assuming steady-state conditions, we show that the meridional overturning circulation of the ocean, driven by heat and salt exchange with the poleward-thinning ice shell, has opposing signs at very low and very high salinities. Regardless of these differing circulations, heat and fresh water converge toward the equator, where the ice is thick, acting to homogenize thickness variations. Among scenarios explored here, the pronounced ice thickness variations observed on Enceladus are most consistent with heating that is predominantly in the ice shell and a salinity of intermediate range.

A perturbation method for evaluating the magnetic field induced from an arbitrary, asymmetric ocean world analytically

Magnetic investigations of icy moons have provided some of the most compelling evidence available confirming the presence of subsurface, liquid water oceans. In the exploration of ocean moons, especially Europa, there is a need for mathematical models capable of predicting the magnetic fields induced under a variety of conditions, including in the case of asymmetric oceans. Existing models are limited to either spherical symmetry or assume an ocean with infinite conductivity. In this work, we use a perturbation method to derive a semi-analytic result capable of determining the induced magnetic moments for an arbitrary layered body, provided each layer is nearly spherical. Crucially, we find that degree-2 tidal deformation results in changes to the induced dipole moments. We demonstrate application of our results to models of plausible asymmetry from the literature within the oceans of Europa and Miranda and the ionospheres of Callisto and Triton. For the models we consider, we find that in the asymmetric case, the induced magnetic field differs by more than 2 nT near the surface of Europa, 0.25-0.5 nT at 1 R above Miranda and Triton, and is essentially unchanged for Callisto. For Miranda and Triton, this difference is as much as 20%-30% of the induced field magnitude. If measurements near the moons can be made precisely to better than a few tenths of a nT, these values may be used by future spacecraft investigations to characterize asymmetry within the interior of icy moons.

Silicate Volcanism on Europa's Seafloor and Implications for Habitability

Habitable ocean environments on Europa require an influx of reactants to maintain chemical disequilibrium. One possible source of reactants is seafloor volcanism. Modeling has shown that dissipation of tidal energy in Europa's asthenosphere can generate melt, but melt formation cannot be equated with volcanism. Melt must also be transported through Europa's cold lithosphere to erupt at the seafloor. Here, we use two models of dike propagation to show that dikes can only traverse the lithosphere if either the fracture toughness of the lithosphere or the flux into the dike is large (>500 MPa m1/2 or ∼1 m2 s-1, respectively). We conclude that cyclic volcanic episodes might provide reactants to Europa's ocean if magma accumulates at the base of the lithosphere for several thousand years. However, if dikes form too frequently, or are too numerous, the magma flux into each will be insufficient, and volcanism cannot support a habitable ocean environment.

Dynamic Europa ocean shows transient Taylor columns and convection driven by ice melting and salinity

The deep (~100 km) ocean of Europa, Jupiter's moon, covered by a thick icy shell, is one of the most probable places in the solar system to find extraterrestrial life. Yet, its ocean dynamics and its interaction with the ice cover have received little attention. Previous studies suggested that Europa's ocean is turbulent using a global model and taking into account non-hydrostatic effects and the full Coriolis force. Here we add critical elements, including consistent top and bottom heating boundary conditions and the effects of icy shell melting and freezing on ocean salinity. We find weak stratification that is dominated by salinity variations. The ocean exhibits strong transient convection, eddies, and zonal jets. Transient motions organize in Taylor columns parallel to Europa's axis of rotation, are static inside of the tangent cylinder and propagate equatorward outside the cylinder. The meridional oceanic heat transport is intense enough to result in a nearly uniform ice thickness, that is expected to be observable in future missions.

A Metamorphic Origin for Europa's Ocean

Europa likely contains an iron-rich metal core. For it to have formed, temperatures within Europa reached ≳ 1250 K. Going up to that temperature, accreted chondritic minerals - for example, carbonates and phyllosilicates - would partially devolatilize. Here, we compute the amounts and compositions of exsolved volatiles. We find that volatiles released from the interior would have carried solutes, redox-sensitive species, and could have generated a carbonic ocean in excess of Europa's present-day hydrosphere, and potentially an early CO 2 atmosphere. No late delivery of cometary water was necessary. Contrasting with prior work, CO 2 could be the most abundant solute in the ocean, followed by Ca 2 + , SO 4 2 - , and HCO 3 - . However, gypsum precipitation going from the seafloor to the ice shell decreases the dissolved S/Cl ratio, such that Cl > S at the shallowest depths, consistent with recently inferred endogenous chlorides at Europa's surface. Gypsum would form a 3-10 km thick sedimentary layer at the seafloor.

Induced magnetic moments from a nearly spherical ocean

The five largest planets all have strong intrinsic magnetic fields that interact with their satellites, many of which contain electrically conducting materials on global scales. Conducting bodies exposed to time-varying magnetic fields induce secondary magnetic fields from movement of eddy currents. In the case of spherically symmetric conducting bodies, matching magnetic solutions at the boundary results in relatively simple relations between the excitation field and the induced field. In this work, we determine the more complicated induced magnetic field from a near-spherical conductor, where the outer boundary is expanded in spherical harmonics. Under the approximations that the excitation field is uniform at a single frequency, the product of wavenumber and radius for the body is large, and the average radius of the body is large compared to the perturbation from spherical symmetry, we find that each spherical harmonic in the shape expansion induces discrete magnetic moments that are independent from the other harmonics in the expansion. That is, simple superposition applies to the magnetic moments induced by each perturbation harmonic. We present a table of the magnetic moments induced by each spherical harmonic up to degree 2 in the perturbed shape. We also present a simple formula by which the induced magnetic field may be evaluated for any arbitrary shape described by expanding the radius of the conducting body in spherical harmonics. Unlike the Earth, many moons in the Solar System are tidally locked to their parent bodies, and many also contain saline, subsurface oceans. Conductive material in these moons is therefore expected to be non-spherical. Accounting for the boundary shape of Europa's ocean will be critical for interpretation of Europa Clipper magnetic measurements near the moon, where the effects of quadrupole-and-higher magnetic moments will be most apparent. The results of this work permit magnetic studies considering non-spherical oceans of satellites for the first time.

Spontaneous formation of geysers at only one pole on Enceladus's ice shell

The ice shell on Enceladus, an icy moon of Saturn, exhibits strong asymmetry between the northern and southern hemispheres, with all known geysers concentrated over the south pole, even though the expected pattern of tidal forced deformation should be symmetric between the north and south poles. Using an idealized ice-evolution model, we demonstrate that this asymmetry may form spontaneously, without any noticeable a priori asymmetry (such as a giant impact or a monopole structure of geological activity), in contrast to previous studies. Infinitesimal asymmetry in the ice shell thickness due to random perturbations are found to be able to grow indefinitely, ending up significantly thinning the ice shell at one of the poles, thereby allowing fracture formation there. Necessary conditions to trigger this hemispheric symmetry-breaking mechanism are found analytically. A rule of thumb we find is that, for Galilean and Saturnian icy moons, the ice shell can undergo hemispheric symmetry breaking only if the mean shell thickness is around 10 to 30 km.

Near/Far Side Asymmetry in the Tidally Heated Moon

Using viscoelastic mass spring model simulations to track heat distribution inside a tidally perturbed body, we measure the near/far side asymmetry of heating in the crust of a spin synchronous Moon in eccentric orbit about the Earth. With the young Moon within. 8 Earth radii of the Earth, we find that tidal heating per unit area in a lunar crustal shell is asymmetric due to the octupole order moment in the Earth's tidal field and is 10 to 20% higher on its near side than on its far side. Tidal heating reduces the crustal basal heat flux and the rate of magma ocean crystallization. Assuming that the local crustal growth rate depends on the local basal heat flux and the distribution of tidal heating in latitude and longitude, a heat conductivity model illustrates that a moderately asymmetric and growing lunar crust could maintain its near/far side thickness asymmetry but only while the Moon is near the Earth.

The surface temperature of Europa

Previous estimates of the annual mean surface temperature of Jupiter's moon, Europa, neglected the effect of the eccentricity of Jupiter's orbit around the Sun, the effect of the emissivity and heat capacity of Europa's ice, the effect of the eclipse of Europa (i.e., the relative time that Europa is within the shadow of Jupiter), the effect of Jupiter's radiation, and the effect of Europa's internal heating. Other studies concentrated on the diurnal cycle but neglected some of the above factors. In addition, to our knowledge, the seasonal cycle of the surface temperature of Europa was not estimated. Here we systematically estimate the diurnal, seasonal and annual mean surface temperature of Europa, when Europa's obliquity, emissivity, heat capacity, and eclipse, as well as Jupiter's radiation, internal heating, and eccentricity, are all taken into account. For a typical internal heating rate of 0.05 W m - 2 , the equator, pole, and the global and mean annual mean surface temperatures are 96 K, 46 K, and 90 K, respectively. We found that the temperature at the high latitudes is significantly affected by the internal heating, especially during the winter solstice, suggesting that measurements of high latitude surface temperatures can be used to constrain the internal heating. We also estimate the incoming solar radiation to Enceladus, the moon of Saturn.

Self-Assembling Ice Membranes on Europa: Brinicle Properties, Field Examples, and Possible Energetic Systems in Icy Ocean Worlds

Brinicles are self-assembling tubular ice membrane structures, centimeters to meters in length, found beneath sea ice in the polar regions of Earth. We discuss how the properties of brinicles make them of possible importance for chemistry in cold environments-including that of life's emergence-and we consider their formation in icy ocean worlds. We argue that the non-ice composition of the ice on Europa and Enceladus will vary spatially due to thermodynamic and mechanical properties that serve to separate and fractionate brines and solid materials. The specifics of the composition and dynamics of both the ice and the ocean in these worlds remain poorly constrained. We demonstrate through calculations using FREZCHEM that sulfate likely fractionates out of accreting ice in Europa and Enceladus, and thus that an exogenous origin of sulfate observed on Europa's surface need not preclude additional endogenous sulfate in Europa's ocean. We suggest that, like hydrothermal vents on Earth, brinicles in icy ocean worlds constitute ideal places where ecosystems of organisms might be found.

Vital Signs: Seismology of Icy Ocean Worlds

Ice-covered ocean worlds possess diverse energy sources and associated mechanisms that are capable of driving significant seismic activity, but to date no measurements of their seismic activity have been obtained. Such investigations could reveal the transport properties and radial structures, with possibilities for locating and characterizing trapped liquids that may host life and yielding critical constraints on redox fluxes and thus on habitability. Modeling efforts have examined seismic sources from tectonic fracturing and impacts. Here, we describe other possible seismic sources, their associations with science questions constraining habitability, and the feasibility of implementing such investigations. We argue, by analogy with the Moon, that detectable seismic activity should occur frequently on tidally flexed ocean worlds. Their ices fracture more easily than rocks and dissipate more tidal energy than the <1 GW of the Moon and Mars. Icy ocean worlds also should create less thermal noise due to their greater distance and consequently smaller diurnal temperature variations. They also lack substantial atmospheres (except in the case of Titan) that would create additional noise. Thus, seismic experiments could be less complex and less susceptible to noise than prior or planned planetary seismology investigations of the Moon or Mars. Key Words: Seismology-Redox-Ocean worlds-Europa-Ice-Hydrothermal. Astrobiology 18, 37-53.

Plume Activity and Tidal Deformation on Enceladus Influenced by Faults and Variable Ice Shell Thickness

We investigated the effect of variations in ice shell thickness and of the tiger stripe fractures crossing Enceladus' south polar terrain on the moon's tidal deformation by performing finite element calculations in three-dimensional geometry. The combination of thinning in the polar region and the presence of faults has a synergistic effect that leads to an increase of both the displacement and stress in the south polar terrain by an order of magnitude compared to that of the traditional model with a uniform shell thickness and without faults. Assuming a simplified conductive heat transfer and neglecting the heat sources below the ice shell, we computed the global heat budget of the ice shell. For the inelastic properties of the shell described by a Maxwell viscoelastic model, we show that unrealistically low average viscosity of the order of 10¹³ Pa s is necessary for preserving the volume of the ocean, suggesting the important role of the heat sources in the deep interior. Similarly, low viscosity is required to predict the observed delay of the plume activity, which hints at other delaying mechanisms than just the viscoelasticity of the ice shell. The presence of faults results in large spatial and temporal heterogeneity of geysering activity compared to the traditional models without faults. Our model contributes to understanding the physical mechanisms that control the fault activity, and it provides potentially useful information for future missions that will sample the plume for evidence of life. Key Words: Enceladus-Tidal deformation-Faults-Variable ice shell thickness-Tidal heating-Plume activity and timing. Astrobiology 17, 941-954.

Science potential from a Europa lander

The prospect of a future soft landing on the surface of Europa is enticing, as it would create science opportunities that could not be achieved through flyby or orbital remote sensing, with direct relevance to Europa's potential habitability. Here, we summarize the science of a Europa lander concept, as developed by our NASA-commissioned Science Definition Team. The science concept concentrates on observations that can best be achieved by in situ examination of Europa from its surface. We discuss the suggested science objectives and investigations for a Europa lander mission, along with a model planning payload of instruments that could address these objectives. The highest priority is active sampling of Europa's non-ice material from at least two different depths (0.5-2 cm and 5-10 cm) to understand its detailed composition and chemistry and the specific nature of salts, any organic materials, and other contaminants. A secondary focus is geophysical prospecting of Europa, through seismology and magnetometry, to probe the satellite's ice shell and ocean. Finally, the surface geology can be characterized in situ at a human scale. A Europa lander could take advantage of the complex radiation environment of the satellite, landing where modeling suggests that radiation is about an order of magnitude less intense than in other regions. However, to choose a landing site that is safe and would yield the maximum science return, thorough reconnaissance of Europa would be required prior to selecting a scientifically optimized landing site.

Thermal diapirism and the habitability of the icy shell of Europa

Europa's chaos and lenticulae features may have originated by thermal diapirs related to convective plumes. Warm ice plumes could be habitable, since their temperature is close to the ice melting temperature. Moreover, thermal plumes intruding into the lower stagnant lid warm several kilometers of country ice above 230 K for periods of 10(5) years, and hundreds of meters above 240 K for periods of 10(4) years. Diapir coalescence generating chaos areas should provide a large zone with temperature above approximately 240 K for thousands of years. A temperature above approximately 230 K is potentially interesting for astrobiology, since it corresponds to the lowest temperature at which microbial metabolic activity in Antarctic ice has been reported. So, the warming by thermal plumes could cause an aureole of biological activation/reactivation in the country ice. Adaptation of life to either high salinity or low temperature is similar: it requires the synthesis of compatible solutes, like trehalose or glycerol, which are efficient cryoprotectants. We therefore propose that the future astrobiological exploration of Europa should include the search for compatible solutes in chaos and lenticulae features.