Jacob Buffo

and 3 more

Compositional heterogeneities within Europa’s ice shell likely impact the dynamics and habitability of the ice and subsurface ocean, but the total inventory and distribution of impurities within the shell is unknown. In sea ice on Earth, the thermochemical environment at the ice-ocean interface governs impurity entrainment into the ice. Here, we simulate Europa’s ice-ocean interface and bound the impurity load (1.053-14.72 g/kg (parts per thousand weight percent, or ppt) bulk ice shell salinity) and bulk salinity profile of the ice shell. We derive constitutive equations that predict ice composition as a function of the ice shell thermal gradient and ocean composition. We show that evolving solidification rates of the ocean and hydrologic features within the shell produce compositional variations (ice bulk salinities of 5-50% of the ocean salinity) that can affect the material properties of the ice. As the shell thickens, less salt is entrained at the ice-ocean interface, which implies Europa’s ice shell is compositionally homogeneous below ~ 1 km. Conversely, the solidification of water filled fractures or lenses introduces substantial compositional variations within the ice shell, creating gradients in mechanical and thermal properties within the ice shell that could help initiate and sustain geological activity. Our results suggest that ocean materials entrained within Europa’s ice shell affect the formation of geologic terrain and that these structures could be confirmed by planned spacecraft observations.

Jacob Buffo

and 2 more

Ocean worlds have been identified as high-priority astrobiology targets due to the link between life and liquid water. Young surface terrain on many icy bodies indicates they support active geophysical cycles that may facilitate ocean-surface transport that could provide observables for upcoming missions. Accurately interpreting spacecraft observations requires constraining the relationship between ice shell characteristics and interior dynamics. On Earth, the composition, physical characteristics, and bioburden of ocean-derived ices are related to their formation history and parent fluid composition. In such systems the ice-ocean interface, which exists as a multiphase mushy layer, dictates the overlying ice’s properties and evolution. Inclusion of the physics governing these boundaries is a novel strategy in modeling planetary ices, and thus far has been limited to 1D approaches. Here we present results from 2D simulations of an archetypal ice-ocean world. We track the evolution of temperature, salinity, porosity, and brine velocity within a thickening ice shell enabling us to place improved constraints on ice-ocean world properties, including: the composition of planetary ice shells, the thickness and hydraulic connectivity of ice-ocean interfaces, and heterogeneous dynamics/structures in the interfacial mushy layer. We show that stable eutectic horizons are likely a common feature of ice-ocean worlds and that ocean composition plays an important role in governing the structure and dynamics of the interface, including the formation of chemical gradient-rich regions within the mushy layer. We discuss the geophysical and astrobiological implications of our results and highlight how they can be validated by instrument specific measurements.

Jacob Buffo

and 11 more

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 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 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 which 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”.

Jacob Buffo

and 5 more

Non-ice impurities within the ice shells of ocean worlds (e.g., Europa, Enceladus, Titan) are believed to play a fundamental role in their geophysics and habitability and may become a surface expression of subsurface ocean properties. Heterogeneous entrainment and distribution of impurities within planetary ice shells have been proposed as mechanisms that can drive ice shell overturn, generate diverse geological features, and facilitate ocean-surface material transport critical for maintaining a habitable subsurface ocean. However, current models of ice shell composition suggest that impurity rejection at the ice-ocean interface of thick contemporary ice shells will be exceptionally efficient, resulting in relatively pure, homogeneous ice. As such, additional mechanisms capable of facilitating enhanced and heterogeneous impurity entrainment are needed to reconcile the observed physicochemical diversity of planetary ice shells. Here we investigate the potential for hydrologic features within planetary ice shells (sills and basal fractures), and the unique freezing geometries they promote, to provide such a mechanism. By simulating the two-dimensional thermal and physicochemical evolution of these hydrological features as they solidify, we demonstrate that bottom-up solidification at sill floors and horizontal solidification at fracture walls generate distinct ice compositions and provide mechanisms for both enhanced and heterogeneous impurity entrainment. We compare our results with magmatic and metallurgic analogs that exhibit similar micro- and macroscale chemical zonation patterns during solidification. Our results suggest variations in ice-ocean/brine interface geometry could play a fundamental role in introducing compositional heterogeneities into planetary ice shells and cryoconcentrating impurities in (re)frozen hydrologic features.