Thickness of a Europan ice shell from impact crater simulations

Geologic Constraints on the Formation and Evolution of Saturn's Mid-Sized Moons

Saturn's mid-sized icy moons have complex relationships with Saturn's interior, the rings, and with each other, which can be expressed in their shapes, interiors, and geology. Observations of their physical states can, thus, provide important constraints on the ages and formation mechanism(s) of the moons, which in turn informs our understanding of the formation and evolution of Saturn and its rings. Here, we describe the cratering records of the mid-sized moons and the value and limitations of their use for constraining the histories of the moons. We also discuss observational constraints on the interior structures of the moons and geologically-derived inferences on their thermal budgets through time. Overall, the geologic records of the moons (with the exception of Mimas) include evidence of epochs of high heat flows, short- and long-lived subsurface oceans, extensional tectonics, and considerable cratering. Curiously, Mimas presents no clear evidence of an ocean within its surface geology, but its rotation and orbit indicate a present-day ocean. While the moons need not be primordial to produce the observed levels of interior evolution and geologic activity, there is likely a minimum age associated with their development that has yet to be determined. Uncertainties in the populations impacting the moons makes it challenging to further constrain their formation timeframes using craters, whereas the characteristics of their cores and other geologic inferences of their thermal evolutions may help narrow down their potential histories. Disruptive collisions may have also played an important role in the formation and evolution of Saturn's mid-sized moons, and even the rings of Saturn, although more sophisticated modeling is needed to determine the collision conditions that produce rings and moons that fit the observational constraints. Overall, the existence and physical characteristics of Saturn's mid-sized moons provide critical benchmarks for the development of formation theories.

Multiring basin formation constrains Europa's ice shell thickness

Jupiter's moon Europa hosts a subsurface ocean under an ice shell of uncertain thickness. Europa has two multiring basins that exhibit several concentric rings. The formation of these multiring basins is thought to be sensitive to the thickness and thermal structure of the ice shell. Here, we simulate multiring basin forming impacts on Europa finding that a total ice shell greater than 20 kilometers thick is required to reproduce observed ring structures. Thin ice shells (<15 kilometers thick) result in compressional tectonics inconsistent with observed ring structures. Our simulations are also sensitive to the thermal structure of the ice shell and indicate that Europa's at least 20-kilometer ice shell is composed of a 6- to 8-kilometer-thick conductive lid overlying warm convecting ice. The constraints on Europa's ice shell structure resulting from this work are directly relevant to our understanding of the potential habitability of Europa.

Ice Shell Structure and Composition of Ocean Worlds: Insights from Accreted Ice on Earth

Accreted ice retains and preserves traces of the ocean from which it formed. In this work, we study two classes of accreted ice found on Earth-frazil ice, which forms through crystallization within a supercooled water column, and congelation ice, which forms through directional freezing at an existing interface-and discuss where each might be found in the ice shells of ocean worlds. We focus our study on terrestrial ice formed in low temperature gradient environments (e.g., beneath ice shelves), consistent with conditions expected at the ice-ocean interfaces of Europa and Enceladus, and we highlight the juxtaposition of compositional trends in relation to ice formed in higher temperature gradient environments (e.g., at the ocean surface). Observations from Antarctic sub-ice-shelf congelation ice and marine ice show that the purity of frazil ice can be nearly two orders of magnitude higher than congelation ice formed in the same low temperature gradient environment (∼0.1% vs. ∼10% of the ocean salinity). In addition, where congelation ice can maintain a planar ice-water interface on a microstructural scale, the efficiency of salt rejection is enhanced (∼1% of the ocean salinity) and lattice soluble impurities such as chloride are preferentially incorporated. We conclude that an ice shell that forms by gradual thickening as its interior cools would be composed of congelation ice, whereas frazil ice will accumulate where the ice shell thins on local (rifts and basal fractures) or regional (latitudinal gradients) scales through the operation of an "ice pump."

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.

Ultrafast visualization of incipient plasticity in dynamically compressed matter

Plasticity is ubiquitous and plays a critical role in material deformation and damage; it inherently involves the atomistic length scale and picosecond time scale. A fundamental understanding of the elastic-plastic deformation transition, in particular, incipient plasticity, has been a grand challenge in high-pressure and high-strain-rate environments, impeded largely by experimental limitations on spatial and temporal resolution. Here, we report femtosecond MeV electron diffraction measurements visualizing the three-dimensional (3D) response of single-crystal aluminum to the ultrafast laser-induced compression. We capture lattice transitioning from a purely elastic to a plastically relaxed state within 5 ps, after reaching an elastic limit of ~25 GPa. Our results allow the direct determination of dislocation nucleation and transport that constitute the underlying defect kinetics of incipient plasticity. Large-scale molecular dynamics simulations show good agreement with the experiment and provide an atomic-level description of the dislocation-mediated plasticity.

Understanding incipient plasticity has been experimentally limited by spatial and temporal resolution. Here the authors report ultra-fast, in situ electron diffraction measurement of dislocation defect dynamics in the early stage of plastic deformation in Al under laser-driven compression.

How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa

Jupiter’s moon Europa harbors one of the most likely environments for extant extraterrestrial life. Determining whether Europa is truly habitable requires understanding the structure and thickness of its ice shell, including the existence of perched water or brines. Stereo-derived topography from images acquired by NASA Galileo’s Solid State Imager (SSI) of Europa are often used as a constraint on ice shell structure and heat flow, but the uncertainty in such topography has, to date, not been rigorously assessed. To evaluate the current uncertainty in Europa’s topography we generated and compared digital terrain models (DTMs) of Europa from SSI images using both the open-source Ames Stereo Pipeline (ASP) software and the commercial SOCET SET® software. After first describing the criteria for assessing stereo quality in detail, we qualitatively and quantitatively describe both the horizontal resolution and vertical precision of the DTMs. We find that the horizontal resolution of the SOCET SET® DTMs is typically 8–11× the root mean square (RMS) pixel scale of the images, whereas the resolution of the ASP DTMs is 9–13× the maximum pixel scale of the images. We calculate the RMS difference between the ASP and SOCET SET® DTMs as a proxy for the expected vertical precision (EP), which is a function of the matching accuracy and stereo geometry. We consistently find that the matching accuracy is ~0.5 pixels, which is larger than well-established “rules of thumb” that state that the matching accuracy is 0.2–0.3 pixels. The true EP is therefore ~1.7× larger than might otherwise be assumed. In most cases, DTM errors are approximately normally distributed, and errors that are several times the derived EP occur as expected. However, in two DTMs, larger errors (differences) occur and correlate with real topography. These differences primarily result from manual editing of the SOCET SET® DTMs. The product of the DTM error and the resolution is typically 4–8 pixel2 if calculated using the RMS image scale for SOCET SET® DTMs and the maximum images scale for the ASP DTMs, which is consistent with recent work using martian data sets and suggests that the relationship applies more broadly. We evaluate how ASP parameters affect DTM quality and find that using a smaller subpixel refinement kernel results in DTMs with smaller (better) resolution but, in some cases, larger gaps, which are sometimes reduced by increasing the size of the correlation kernel. We conclude that users of ASP should always systematically evaluate the choice of parameters for a given dataset.

Star-shaped crack pattern of broken windows

Broken thin brittle plates like windows and windshields are ubiquitous in our environment. When impacted locally, they typically present a pattern of cracks extending radially outward from the impact point. We study the variation of the pattern of cracks by performing controlled transverse impacts on brittle plates over a broad range of impact speed, plate thickness, and material properties, and we establish from experiments a global scaling law for the number of radial cracks incorporating all these parameters. A model based on Griffith's theory of fracture combining bending elastic energy and fracture energy accounts for our observations. These findings indicate how the postmortem shape of broken samples are related to material properties and impact parameters, a procedure relevant to forensic science, archaeology, or astrophysics.

Exploration and protection of Europa's biosphere: implications of permeable ice

Europa has become a high-priority objective for exploration because it may harbor life. Strategic planning for its exploration has been predicated on an extreme model in which the expected oceanic biosphere lies under a thick ice crust, buried too deep to be reached in the foreseeable future, which would beg the question of whether other active satellites might be more realistic objectives. However, Europa's ice may in fact be permeable, with very different implications for the possibilities for life and for mission planning. A biosphere may extend up to near the surface, making life far more readily accessible to exploration while at the same time making it vulnerable to contamination. The chances of finding life on Europa are substantially improved while the need for planetary protection becomes essential. The new National Research Council planetary protection study will need to go beyond its current mandate if meaningful standards are to be put in place.

Shock deformation of face-centred-cubic metals on subnanosecond timescales

Despite its fundamental importance for a broad range of applications, little is understood about the behaviour of metals during the initial phase of shock compression. Here, we present molecular dynamics (MD) simulations of shock-wave propagation through a metal allowing a detailed analysis of the dynamics of high strain-rate plasticity. Previous MD simulations have not seen the evolution of the strain from one- to three-dimensional compression that is observed in diffraction experiments. Our large-scale MD simulations of up to 352 million atoms resolve this important discrepancy through a detailed understanding of dislocation flow at high strain rates. The stress relaxes to an approximately hydrostatic state and the dislocation velocity drops to nearly zero. The dislocation velocity drop leads to a steady state with no further relaxation of the lattice, as revealed by simulated X-ray diffraction.

Clathrate hydrates of oxidants in the ice shell of Europa

Europa's icy surface is radiolytically modified by high-energy electrons and ions, and photolytically modified by solar ultraviolet photons. Observations from the Galileo Near Infrared Mapping Spectrometer, ground-based telescopes, the International Ultraviolet Explorer, and the Hubble Space Telescope, along with laboratory experiment results, indicate that the production of oxidants, such as H2O2, O2, CO2, and SO2, is a consequence of the surface radiolytic chemistry. Once created, some of the products may be entrained deeper into the ice shell through impact gardening or other resurfacing processes. The temperature and pressure environments of regions within the europan hydrosphere are expected to permit the formation of mixed clathrate compounds. The formation of carbon dioxide and sulfur dioxide clathrates has been examined in some detail. Here we add to this analysis by considering oxidants produced radiolytically on the surface of Europa. Our results indicate that the bulk ice shell could have a approximately 1.7-7.6% by number contamination of oxidants resulting from radiolysis at the surface. Oxidant-hosting clathrates would consequently make up approximately 12-53% of the ice shell by number relative to ice, if oxidants were entrained throughout. We examine, in brief, the consequences of such contamination on bulk ice shell thickness and find that clathrate formation could lead to substantially thinner ice shells on Europa than otherwise expected. Finally, we propose that double occupancy of clathrate cages by O2 molecules could serve as an explanation for the observation of condensed-phase O2 on Europa. Clathrate-sealed, gas-filled bubbles in the near surface ice could also provide an effective trapping mechanism, though they cannot explain the 5771 A (O2)2 absorption.

Locating potential biosignatures on Europa from surface geology observations

We evaluated the astrobiological potential of the major classes of geologic units on Europa with respect to possible biosignatures preservation on the basis of surface geology observations. These observations are independent of any formational model and therefore provide an objective, though preliminary, evaluation. The assessment criteria include high mobility of material, surface concentration of non-ice components, relative youth, textural roughness, and environmental stability. Our review determined that, as feature classes, low-albedo smooth plains, smooth bands, and chaos hold the highest potential, primarily because of their relative young age, the emplacement of low-viscosity material, and indications of material exchange with the subsurface. Some lineaments and impact craters may be promising sites for closer study despite the comparatively lower astrobiological potential of their classes. This assessment will be expanded by multidisciplinary examination of the potential for habitability of specific features.

Thickness constraints on the icy shells of the galilean satellites from a comparison of crater shapes

A thin outer ice shell on Jupiter's large moon Europa would imply easy exchange between the surface and any organic or biotic material in its putative subsurface ocean. The thickness of the outer ice shell is poorly constrained, however, with model-dependent estimates ranging from a few kilometres to ten or more kilometres. Here I present measurements of depths of impact craters on Europa, Ganymede and Callisto that reveal two anomalous transitions in crater shape with diameter. The first transition is probably related to temperature-dependent ductility of the crust at shallow depths (7 8 km on Europa). The second transition is attributed to the influence of subsurface oceans on all three satellites, which constrains Europa's icy shell to be at least 19 km thick. The icy lithospheres of Ganymede and Callisto are equally ice-rich, but Europa's icy shell has a thermal structure about 0.25 0.5 times the thicknesses of Ganymede's or Callisto's shells, depending on epoch. The appearances of the craters on Europa are inconsistent with thin-ice-shell models and indicate that exchange of oceanic and surface material could be difficult.

Energy cycling and hypothetical organisms in Europa's ocean

While Europa has emerged as a leading candidate for harboring extraterrestrial life, the apparent lack of a source of free energy for sustaining living systems has been argued. In this theoretical analysis, we have quantified the amount of energy that could in principle be obtained from chemical cycling, heat, osmotic gradients, kinetic motion, magnetic fields, and gravity in Europa's subsurface ocean. Using reasonable assumptions based on known organisms on Earth, our calculations suggest that chemical oxidation-reduction cycles in Europa's subsurface ocean could support life. Osmotic and thermal gradients, as well as the kinetic energy of convection currents, also represent plausible alternative sources of energy for living systems at Europa. Organisms thriving on these gradients could interact with each other to form the complex energy cycling necessary for establishing a stable ecosystem.