Preliminary studies concerning subsurface probes for the exploration of icy planetary bodies

It is well-known that the permanent terrestrial ice sheets (glaciers and polar caps) contain a lot of information about the recent geological history and in particular about climatic changes. Extrapolating this fact to other ice sheets in the solar system (e.g. the Mars polar regions, the icy moons of the outer planets, etc.), we may expect a similar wealth of information. To obtain this information it is possible to drill holes or melt the ice by a heated probe, which in this way is able to penetrate the surface and investigate the deeper layers in situ. In the latter case the driving agent is the heating power and the weight of the probe. In this paper we consider the application of such “melting probes” for exploring the structure of ice sheets in extraterrestrial environments. We describe several laboratory experiments with simple melting probes performed under cryo-vacuum conditions and compare the results with tests in a terrestrial environment. The experiments revealed that under space conditions the downward motion of a heated probe in an ice sheet is characterized by intermittent periods of sublimation and melting of the surrounding ice, sometimes interrupted by periods where a part of the probe's outer surface is frozen to the surrounding ice. This leads to a temporary blocking of the probe's downward motion. A similar situation can occur when the trailing tether is frozen in behind the probe. During the periods of ice sublimation the penetration process is significantly more power consuming, due to the large difference between the latent heat of sublimation and the latent heat of melting for water ice.     

Study of a thermal drill head for the exploration of subsurface planetary ice layers

The recently discovered water vapor plumes on Saturn's moon Enceladus, the polar caps of planet Mars and the possible ice volcanism on the Jovian satellites call for suitable techniques to explore deep ice layers of the solar system bodies. This paper presents a novel approach to deliver scientific probes into deeper layers of planetary ice. Several existing locomotion concepts and techniques for such probes are presented. After studying the mathematical framework of the melting locomotion process, melting tests with different head forms were done to evaluate the influence of the head's geometry on the melting process. This work led to a novel concept of a thermal drill head, using heat and mechanical drill in combination to penetrate the ice. We compare the performance of such a hybrid concept versus the melting penetration alone by a mathematical model and tests in ice with a prototype of the melting drill head.     

Recent rheologic processes on dark polar dunes of Mars: Driven by interfacial water?

In springtime on HiRISE images of the Southern polar terrain of Mars flow-like or rheologic features were observed. Their dark color is interpreted as partly defrosted surface where the temperature is too high for CO₂ but low enough for H₂O ice to be present there. These branching streaks grow in size and can move by an average velocity of up to about 1 m/day and could terminate in pond-like accumulation features. The phenomenon may be the result of interfacial water driven rheologic processes. Liquid interfacial water can in the presence of water ice exist well below the melting point of bulk water, by melting in course of interfacial attractive pressure by intermolecular forces (van der Waals forces e.g.), curvature of water film surfaces, and e.g. by macroscopic weight, acting upon ice. This melting phenomenon can be described in terms of “premelting of ice”. It is a challenging consequence, that liquid interfacial water unavoidably must in form of nanometric layers be present in water ice containing soil in the subsurface of Mars. It is the aim of this paper to study possible rheologic consequences in relation to observations, which seem to happen at sites of dark polar dunes on Mars at present. The model in this work assumes that interfacial water accumulates at the bottom of a translucent water-ice layer above a dark and insolated ground. This is warmed up towards the melting point of water. The evolving layer of liquid interfacial water between the covering ice sheet and the heated ground is assumed to drive downward directed flow-like features on slopes, and it can, at least partially, infiltrate (seep) into a porous ground. There, in at least temporarily cooler subsurface layers, the infiltrated liquid water refreezes and forms ice. The related stress built-up is shown to be sufficient to cause destructive erosive processes. The above-mentioned processes may cause change in the structure and thickness of the covering ice and/or may cause the movement of dune grains. All these processes may explain the observed springtime growing and downward extension of the slope streaks analyzed here.     

Seas under ice: Stability of liquid-water oceans within icy worlds

The present-day existence of internal oceans under the outer ice shell of several icy satellites of the Solar System has been recently proposed. The presence of antifreeze substances decreasing ice’s melting point (and tidal heating in Europa’s case) has been generally believed to allow the stability of such oceans; limited cooling of the water (ice plus liquid) layer, due to stability against convection or to stagnant lid convection in the icy shell, have been also considered. Here we propose that even pure liquid-water oceans could survive today within several icy worlds, and we consider some factors affecting thermal modeling in these bodies. So, the existence of such oceans would be a natural consequence of the physical properties of water ice, independently from the addition of antifreeze substances or any other special conditions. The inclusion of these substances would contribute to expand the conditions for water to stay liquid and to increase ocean’s volume.     

Combined modeling of thermal evolution and accretion of trans-neptunian objects—Occurrence of high temperatures and liquid water

We have calculated the early thermal evolution of trans-neptunian objects by means of a thermal evolution code that takes into account simultaneous accretion. The set of coupled partial differential equations for ²⁶Al radioactive heating, transformation of amorphous to crystalline ice and melting of water ice was solved numerically for small porous icy (cometary-like) bodies growing to final radii between 2 and 32 km and accreting between 20 and 44 AU. Accretion within a swarm of gravitationally interacting small bodies was calculated self-consistently with a simple accretion algorithm and thermal evolution of a typical member of the swarm was tracked in a parameter-space survey. We find that including accretion in numerical modeling of thermal evolution leads to a broad variety of thermally processed icy bodies and that the early occurrence of liquid water and extended crystalline ice interiors may be a very common phenomenon. The pristine nature of small icy bodies becomes thus restricted to a particular set of initial conditions. Generally, long-period comets should be more thermally affected than short-period ones.     

Efficiency of neutrino-induced radio measurements to inspect local areas of Enceladus

The Cassini probe observed a young and smooth surface around the south pole of Enceladus, while around the north pole the surface was found to be relatively old and inactive (Porco, C.C. et al. [2006]. Science 311, 1393–1401). This heterogeneous surface implies that the ice thickness of Enceladus is not uniform between the north and south polar regions. Determining the thickness of the icy layer is important to confirm the existence of an internal ocean as well as to reveal the heating mechanism of Enceladus. We show that the measurement of radio waves induced by cosmic neutrinos can be an effective method to constrain the ice thickness of a localized area where conventional gravity or electromagnetic field measurements cannot be used. This method could be used to constrain the thickness of the icy layer on Enceladus even if the ice is a few tens of kilometers thick, measuring over a period of several years, which greatly exceeds the ability of radar sounding, and hence could be used in future orbiter missions.     

Impact crater formation in icy layered terrains on Mars

Abstract— We present numerical simulations of crater formation under Martian conditions with a single near‐surface icy layer to investigate changes in crater morphology between glacial and interglacial periods. The ice fraction, thickness, and depth to the icy layer are varied to understand the systematic effects on observable crater features. To accurately model impact cratering into ice, a new equation of state table and strength model parameters for H₂O are fitted to laboratory data. The presence of an icy layer significantly modifies the cratering mechanics. Observable features demonstrated by the modeling include variations in crater morphometry (depth and rim height) and icy infill of the crater floor during the late stages of crater formation. In addition, an icy layer modifies the velocities, angles, and volumes of ejecta, leading to deviations of ejecta blanket thickness from the predicted power law. The dramatic changes in crater excavation are a result of both the shock impedance and the strength mismatch between layers of icy and rocky materials. Our simulations suggest that many of the unusual features of Martian craters may be explained by the presence of icy layers, including shallow craters with well‐preserved ejecta blankets, icy flow related features, some layered ejecta structures, and crater lakes. Therefore, the cratering record implies that near‐surface icy layers are widespread on Mars.     

Earth analog for high-latitude landforms and recent flows on Mars: Icy debris fans in the Wrangell Volcanic Field, Alaska

New Mars Reconnaissance Orbiter HiRise and CRISM imagery of polar layered terrain of Mars reveals striking similarities to icy debris fans along the base of steep escarpments in Alaska formed in high-latitude periglacial environments. Process and morphologic observations of a deglaciating site in the Wrangell Mountains reveal a complex suite of supraglacial processes involved in the construction of icy debris fans. Snow, ice, and sediment are delivered to the fans from degradation of an upper-level icecap. Alaskan icy debris fans were studied during an 8-day reconnaissance mission in July 2006. We directly observed 289 major depositional events dominated by dry snow/ice avalanches, but also including icy debris flows, rockfalls, small jokulhlaups, and glacial calving. Small fans with larger catchments receive episodic icy debris flows triggered by outburst flows that mobilize rockfall sediment temporarily stored in catchments above the fan apex. Large fans with smaller catchments have better linkage to the upper icecap, providing a direct pathway for frequent large avalanches. The large, avalanche-dominated fans thicken rapidly from an overabundance of snow/ice supply to the point where they become hybrid fan-glaciers. Surficial geology evolves rapidly in this high-latitude environment through both depositional events and solar-driven albedo changes that occur daily. Ground penetrating radar surveys show that subsurface sedimentary architecture and fan evolution is similar to the active surface processes and deposits observed on the fans. Direct field observations of active geomorphic processes provide unique insights on the pace and nature of high-latitude landscape evolution during climate changes on both planets.     

Growth and evolution of small porous icy bodies with an adaptive-grid thermal evolution code: I. Application to Kuiper belt objects and Enceladus

We present a new 1-dimensional thermal evolution code suited for small icy bodies of the Solar System, based on modern adaptive grid numerical techniques, and suited for multiphase flow through a porous medium. The code is used for evolutionary calculations spanning 4.6×¹⁰⁹ yr of a growing body made of ice and rock, starting with a 10 km radius seed and ending with an object 250 km in radius. Initial conditions are chosen to match two different classes of objects: a Kuiper belt object, and Saturn's moon Enceladus. Heating by the decay of ²⁶Al, as well as long-lived radionuclides is taken into account. Several values of the thermal conductivity and accretion laws are tested. We find that in all cases the melting point of ice is reached in a central core. Evaporation and flow of water and vapor gradually remove the water from the core and the final (present) structure is differentiated, with a rocky, highly porous core of 80 km radius (and up to 160 km for very low conductivities). Outside the core, due to refreezing of water and vapor, a compact, ice-rich layer forms, a few tens of km thick (except in the case of very high conductivity). If the ice is initially amorphous, as expected in the Kuiper belt, the amorphous ice is preserved in an outer layer about 20 km thick. We conclude by suggesting various ways in which the code may be extended.     

Thermal evolution of trans‐Neptunian objects,icy satellites,and minor icy planets in the early solar system

Numerical simulations are performed to understand the early thermal evolution and planetary scale differentiation of icy bodies with the radii in the range of 100–2500 km. These icy bodies include trans‐Neptunian objects, minor icy planets (e.g., Ceres, Pluto); the icy satellites of Jupiter, Saturn, Uranus, and Neptune; and probably the icy‐rocky cores of these planets. The decay energy of the radionuclides, ²⁶Al, ⁶⁰Fe, ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th, along with the impact‐induced heating during the accretion of icy bodies were taken into account to thermally evolve these planetary bodies. The simulations were performed for a wide range of initial ice and rock (dust) mass fractions of the icy bodies. Three distinct accretion scenarios were used. The sinking of the rock mass fraction in primitive water oceans produced by the substantial melting of ice could lead to planetary scale differentiation with the formation of a rocky core that is surrounded by a water ocean and an icy crust within the initial tens of millions of years of the solar system in case the planetary bodies accreted prior to the substantial decay of ²⁶Al. However, over the course of billions of years, the heat produced due to ⁴⁰K, ²³⁵U, ²³⁸U, and ²³²Th could have raised the temperature of the interiors of the icy bodies to the melting point of iron and silicates, thereby leading to the formation of an iron core. Our simulations indicate the presence of an iron core even at the center of icy bodies with radii ≥500 km for different ice mass fractions.     

Experimental study on collisional disruption of highly porous icy bodies

Knowing the collisional process among small porous icy bodies in the outer solar system is a key to understanding the formation of EKBOs and the evolution of icy planetesimals. Impact experiments of sintered porous ice spheres with 40%, 50%, 60% and 70% porosity were conducted by using three types of projectiles at the impact velocity from 2.4 to 489 m/s, and we studied the effects of porosity on the collisional processes. Projectile sticking occurred at the impact velocity higher than 44 m/s for 60% porosity targets and higher than 13 m/s for 70% porosity targets. The antipodal velocity of the porous ice target increased with the increase of energy density, Q, and it increased slightly with the increase of porosity, although it was exceptionally high in cases when the projectile penetrated the target. The shattering strength of porous ice targets was found to decrease from 100 to 31 J/kg with the increase of porosity from 40% to 70%. The cumulative fragment mass distribution was found to depend on the energy density and the target porosity, and the slopes of the distribution in the small fragment region were almost flat for more porous targets. We reanalyzed the cumulative fragment mass distribution and first obtained the empirical equation showing the fragment mass distribution of porous ice targets as a function of the energy density and the porosity.     

Ice grain size and the rheology of the martian polar deposits

The history and dynamics of the martian polar deposits (MPD), the largest known water reservoirs on Mars, are of great interest, but estimates of ice grain size are required before detailed modeling can be performed. We clarify the microphysical processes that may control grain size in the MPD. If the MPD are ∼2% dust by mass, the maximum ice grain size is ∼1 mm due to grain boundary pinning by silicate microparticles. Relatively dusty layers in the MPD will have smaller grain sizes. If MPD ice has a very low impurity content and has experienced a significant amount of strain, grains may reach a steady state size of ∼1.5 to 3 mm due to dynamic recrystallization, wherein a steady state grain size is maintained due to the balance of grain growth and destruction during flow. If the near-bed ice in the MPD is warmed close to its melting point and has been extensively sheared, grain sizes at its base may be between 10 and 40 mm, by analogy with warm, dirty, near-bed ice in terrestrial ice sheets.     

The Use of Geophysical Data in Studies of the Historical Solar-Terrestrial Environment

Recent studies of the solar-terrestrial environment for the past 500 years have necessitated the use of a variety of historical databases: nitrates in ice cores, knowledge of large volcanic eruptions, sunspot numbers, mid-latitude aurora and geomagnetic records. The nitrate data are being used to identify large solar proton fluence events. The volcanic record helps to provide time markers for the ice core. The records of major geomagnetic storms and mid-latitude aurora have been used for additional identification. We also know that the Earth’s magnetic field is evolving with a present rapid decrease in magnitude. In addition the wandering magnetic pole must be considered in ascertaining what was “mid latitude” in historic times versus “mid latitude” in 2000. We illustrate how these databases are being used in recent studies of historic solar proton events.     

Constraints on the composition of the martian south polar cap from gravity and topography

The polar caps of Mars have long been acknowledged to be composed of unknown proportions of water ice, solid CO₂ (dry ice), and dust. Gravity and topography data are here analyzed over the southern cap to place constraints on its density, and hence composition. Using a localized spectral analysis combined with a lithospheric flexure model of ice cap loading, the best fit density of the volatile-rich south polar layered deposits is found to be 1271 kg m⁻³ with 1-σ limits of 1166 and 1391 kg m⁻³. The best fit elastic thickness of this geologically young deposit is 140 km, though any value greater than 102 km can fit the observations. The best fit density implies that about 55% dry ice by volume could be sequestered in these deposits if they were completely dust free. Alternatively, if these deposits were completely free of solid CO₂, the dust content would be constrained to lie between about 14 and 28% by volume. The bulk thermal conductivity of the polar cap is not significantly affected by these maximum allowable concentrations of dust. However, even if a moderate quantity of solid CO₂ were present as horizontal layers, the bulk thermal conductivity of the polar cap would be significantly reduced. Reasonable estimates of the present day heat flow of Mars predict that dry ice beneath the thicker portions of the south polar cap would have melted. Depending on the quantity of solid CO₂ in these deposits today, it is even possible that water ice could melt where the cap is thickest. If independent estimates for either the dust or CO₂ content of the south polar cap could be obtained, and if radar sounding data could determine whether this polar cap is presently experiencing basal melting or not, it would be possible to use these observations to place tight constraints on the present day heat flow of Mars.     

Water ice polymorphs and their significance on planetary surfaces

Impacts into an icy surface could produce significant amounts of high pressure forms of water ice. Due to the relatively low ambient surface temperatures on satellites in the outer solar system and the modest temperature rises accompanying the impact pressures required for water ice metamorphism, high-pressure polymorphs will be created by and may remain after large cratering events. If so, those high-pressure ices should be ubiquitous. Low-pressure cubic ice may be abundant as well. Impacts into an icy regolith may both produce high-pressure polymorphs from ice I and destroy high-pressure polymorphs already present. The result will be an (unknown) equilibrium concentration of high pressure polymorphs in the regolith. Polymorphs may be detectable and mappable by reflection spectroscopy at vacuum ultraviolet and mid-infrared wavelengths.     

Radiative habitable zones in martian polar environments

The biologically damaging solar ultraviolet (UV) radiation (quantified by the DNA-weighted dose) reaches the martian surface in extremely high levels. Searching for potentially habitable UV-protected environments on Mars, we considered the polar ice caps that consist of a seasonally varying CO2 ice cover and a permanent H2O ice layer. It was found that, though the CO2 ice is insufficient by itself to screen the UV radiation, at approximately 1 m depth within the perennial H2O ice the DNA-weighted dose is reduced to terrestrial levels. This depth depends strongly on the optical properties of the H2O ice layers (for instance snow-like layers). The Earth-like DNA-weighted dose and Photosynthetically Active Radiation (PAR) requirements were used to define the upper and lower limits of the northern and southern polar Radiative Habitable Zone (RHZ) for which a temporal and spatial mapping was performed. Based on these studies we conclude that photosynthetic life might be possible within the ice layers of the polar regions. The thickness varies along each martian polar spring and summer between approximately 1.5 and 2.4 m for H2O ice-like layers, and a few centimeters for snow-like covers. These martian Earth-like radiative habitable environments may be primary targets for future martian astrobiological missions. Special attention should be paid to planetary protection, since the polar RHZ may also be subject to terrestrial contamination by probes.     

Topography on satellite surfaces and the shape of asteroids

Calculations of the topography and shape of planetary bodies are presented for two sets of models. One set of models deals with the effects of static loading on bodies, taking into account strengths of materials, density, and size. The other set considers the effects of creep deformation on model bodies of differing composition, size and temperature. Application of these models to asteroids and satellites of the major planets indicates that model, even the largest asteroids could retain highly nonspherical shapes, and the four large satellites of Jupiter could sustain statically loaded topography on the order of 10 km. (2) If silicate asteroids have not been heated to near the melting temperature of silicates, initial topography should survive for at least 109 yr under creep deformation. Topography on an insulated icy asteroid will be rapidly reduced if it is of larger scale than the insulating layer, no matter what the thermal history. (3) Of the Galilean satellites of Jupiter, J1 and J2 should retain topography created on silicate surfaces since their formation (or since the surfaces were near the silicate melting temperature. If ice layers of any significant thickness exist, topography on a scale smaller than the layer's thickness will be reduced rapidly. (4) J4 and J3 probably fit an icy model throughout and topography of all scales may be reduced with relaxation times < 10⁶yr. These satellites are thus likely to preserve only very recent features on their surfaces, in contrast to the other Galilean satellites. If melting has taken place since formation, these conclusions become even stronger. (5) Of the satellites of the other planets, only Titan appears likely to have undergone topographic reduction by creep, under the models presented. However, if ices other than water are present in large proportion on these satellites relaxation times for topography may be shorter than calculated from the water ice models.     

Origin of ice diapirism, true polar wander, subsurface ocean, and tiger stripes of Enceladus driven by compositional convection

We consider the scenario in which the presence of ammonia in the bulk composition of Enceladus plays a pivotal role in its thermochemical evolution. Because ammonia reduces the melting temperature of the ice shell by 100 K below that of pure water ice, small amounts of tidal dissipation can power an “ammonia feedback” mechanism that leads to secondary differentiation of Enceladus within the ice shell. This leads to compositionally distinct zones at the base of the ice shell arranged such that a layer of lower density (and compositionally buoyant) pure water ice underlies the undifferentiated ammonia-dihydrate ice layer above. We then consider a large scale instability arising from the pure water ice layer, and use a numerical model to explore the dynamics of compositional convection within the ice shell of Enceladus. The instability of the layer can easily account for a diapir that is hemispherical in scale. As it rises to the surface, it co-advects the warm internal temperatures towards the outer layers of the satellite. This advected heat facilitates the generation of a subsurface ocean within the ice shell of Enceladus. This scenario can simultaneously account for the origin of asymmetry in surface deformation observed on Enceladus as well as two global features inferred to exist: a large density anomaly within the interior and a subsurface ocean underneath the south polar region.     

Laboratory simulation experiments on the solid-state greenhouse effect in planetary ices

Icy surfaces like the polar caps of Mars, comets, Edgeworth-Kuiper belt objects or the surface areas of many moons in the outer Solar System behave different than rock and soil surfaces when irradiated by solar light. The latter ones absorb and reflect incoming solar radiation immediately at the surface. In contrast, ices are partially transparent in the visible spectral range and opaque in the infrared. Due to this fact it is possible for the solar radiation to reach a certain depth and increase the temperature of the sub-surface layers directly. This internal temperature rise is called “solid-state greenhouse effect,” in analogy to the classical greenhouse effect in an atmosphere. It may play an important role in the energy balance of various icy bodies in the Solar System. Within the scope of a project conducted at the Space Research Institute of the Austrian Academy of Sciences in Graz the solid-state greenhouse effect was investigated experimentally and theoretically. A number of experiments with diverse materials, focussing mainly on layered samples with a surface cover consisting of transparent H₂O-ice, were performed. The samples were irradiated under cryo-vacuum conditions by a solar simulator. The temperature distributions inside the samples were measured and compared with the results of numerical model calculations. We found that the predicted sub-surface temperature maximum is very clearly measurable in glass beads samples with various particle size distributions, but can also be detected in transparent compact surface ice layers. However, in the latter case it is less distinct than originally expected. Measuring the effect by laboratory methods turned out to be a difficult task due to the shallow depth where the temperature maximum occurs.