Impacts into oceans and seas

Impacts of cosmic bodies into oceans and seas lead to the formation of very high waves. Numerical simulations of 3-km and 1-km comets impacting into a 4 km depth ocean with a velocity of 20 km/sec have been conducted. For a 1-km body, depth of the interim crater in the sea bed is about 8 km below ocean level, and the height of the water wave is 10 m at a distance of 2000 km from the impact point. As the water wave runs into shallows, a huge tsunami hits the coast. The height of the wave strongly depends on the coastal and sea bed topography. If the impact occurred near the shore, the huge mass of water strikes the cliffs and the near shore mountain ridges and can cause displacement of the rocks, initiate landslides, and change the relief. Thus, impact into oceans and seas is an important geological factor. Cosmic bodies of small sizes are disrupted by aerodynamic forces. Fragments of a 100-m radius comet striking the water surface create an unstable cavity in the water of about 1 km radius. Its collapse also creates tsunami. A simple estimate has been made using the light curves from recent atmosphere explosions detected by satellites. The results of our assessment of the characteristics of meteoroids which caused these intense light flashes suggests that fragments of a 25-m stony body with initial impact velocity 15 to 20 km/sec will hit the surface. For a 75-m iron body striking the sea with a depth of 600 m, the height of the wave is 10 m at 200–300 km distance from the impact.     

Impacts into weak dust targets under microgravity and the formation of planetesimals

We carried out 16 collision experiments in the drop tower in Bremen, Germany. Dust projectiles and solid projectiles of several mm in size impacted a dust target 5 cm in depth and width at velocities between 3.5 and 21.5 m/s. For solid impactors we found significant mass loss on the front (impact) side of the target. Mass loss depended on the impact velocity and projectile type (solid sphere or dust) and was up to 35 times the projectile mass for targets of the lowest tensile strength. Typical fragment velocities on the front side of the target ranged from 3 to 12 cm/s. The ejecta velocity was independent of the impact velocity but it increased with projectile mass. On the back side of the target (opposite to the impact side) mass was ejected from the target above a certain threshold impact velocity. Ejection velocity on the back side increased with impact velocity and is larger for solid projectiles than for dust projectiles. In one case a slightly stronger target gained mass in a slow dust-dust collision. We verified that collisions of dust projectiles with compact, very strong dust targets lead to a more massive target accreting part of the projectile. Applied to planetesimal formation, the experiments suggest that the maximum possible ejecta velocity from a body of several cm in size after a collision is small. Ejecta were slow enough that they were reaccreted by means of gas flow if large pores were part of the body's morphology. While very weak bodies cannot grow in the primary collision at the given velocities, this can lead to growth by secondary collisions. Slight compression, which could result from preceding collisions, might lead to immediate growth of a body in slow collisions by adding projectile mass.     

Liquid water oceans in ice giants

Aptly named, ice giants such as Uranus and Neptune contain significant amounts of water. While this water cannot be present near the cloud tops, it must be abundant in the deep interior. We investigate the likelihood of a liquid water ocean existing in the hydrogen-rich region between the cloud tops and deep interior. Starting from an assumed temperature at a given upper tropospheric pressure (the photosphere), we follow a moist adiabat downward. The mixing ratio of water to hydrogen in the gas phase is small in the photosphere and increases with depth. The mixing ratio in the condensed phase is near unity in the photosphere and decreases with depth; this gives two possible outcomes. If at some pressure level the mixing ratio of water in the gas phase is equal to that in the deep interior, then that level is the cloud base. The gas below the cloud base has constant mixing ratio. Alternately, if the mixing ratio of water in the condensed phase reaches that in the deep interior, then the surface of a liquid ocean will occur. Below this ocean surface, the mixing ratio of water will be constant. A cloud base occurs when the photospheric temperature is high. For a family of ice giants with different photospheric temperatures, the cooler ice giants will have warmer cloud bases. For an ice giant with a cool enough photospheric temperature, the cloud base will exist at the critical temperature. For still cooler ice giants, ocean surfaces will result. A high mixing ratio of water in the deep interior favors a liquid ocean. We find that Neptune is both too warm (photospheric temperature too high) and too dry (mixing ratio of water in the deep interior too low) for liquid oceans to exist at present. To have a liquid ocean, Neptune's deep interior water to gas ratio would have to be higher than current models allow, and the density at 19 kbar would have to be ≈0.8 g/cm³. Such a high density is inconsistent with gravitational data obtained during the Voyager flyby. In our model, Neptune's water cloud base occurs around 660 K and 11 kbar, and the density there is consistent with Voyager gravitational data. As Neptune cools, the probability of a liquid ocean increases. Extrasolar “hot Neptunes,” which presumably migrate inward toward their parent stars, cannot harbor liquid water oceans unless they have lost almost all of the hydrogen and helium from their deep interiors.     

Fluid dynamics of local martian magma oceans

The evolution of a melt region produced by a large impact during Mars formation is addressed. While some impact induced melt is redistributed during crater excavation, sufficiently large impacts (much larger than basin forming impacts) generate an intact melt region which is retained beneath the excavation zone, i.e., a local magma ocean. Local magma ocean evolution depends on the effective rheology controlling large scale deformation of the solid part of the planet, mechanism of crystallization, and melt region size. Within the uncertainties of various parameters, two scenarios are possible. For sufficiently weak rheology or large melt region size, evolution is characterized by rapid extrusion and formation of a global magma ocean. For sufficiently strong rheology or small melt region size, in situ crystallization to a partially molten solid state occurs prior to isostatic adjustment. Subsequent to in situ crystallization, local magma ocean evolution depends on melt region size and efficiency of lateral redistribution compared to bulk conductive cooling. For large melt regions, lateral spreading occurs via plastic deformation and results in an asymmetric, global, partial melt layer. For small melt region size, viscous spreading viscous can result in bulk cooling below the solidus prior to formation of a global layer. A hypothesis for the origin of the hemispherical crustal dichotomy and Tharsis rise is suggested. The dichotomy is associated with a global partial melt layer produced by evolution of a large, local magma ocean. After dichotomy formation, evolution of a second, smaller, local magma ocean is related to Tharsis development.     

Formation of early water oceans on rocky planets

Terrestrial planets, with silicate mantles and metallic cores, are likely to obtain water and carbon compounds during accretion. Here I examine the conditions that allow early formation of a surface water ocean (simultaneous with cooling to clement surface conditions), and the timeline of degassing the planetary interior into the atmosphere. The greatest fraction of a planet’s initial volatile budget is degassed into the atmosphere during the end of magma ocean solidification, leaving only a small fraction of the original volatiles to be released into the atmosphere through later volcanism. Rocky planets that accrete with water in their bulk mantle have two mechanisms for producing an early water ocean: First, if they accrete with at least 1 to 3 mass% of water in their bulk composition, liquid water may be extruded onto the planetary surface at the end of magma ocean solidification. Second, at initial water contents as low as 0.01 mass% or lower, during solidification a massive supercritical fluid and steam atmosphere is produced that collapses into a water ocean upon cooling. The low water contents required for this process indicate that rocky super-Earth exoplanets may be expected to commonly produce water oceans within tens to hundreds of millions of years of their last major accretionary impact, through collapse of their atmosphere.     

Sub-ice volcanoes and ancient oceans/lakes: a Martian challenge

New instruments on board the Mars Global Surveyor (MGS) spacecraft began providing accurate, high-resolution image and topography data from the planet in 1997. Though data from the Mars Orbiter Laser Altimeter (MOLA) are consistent with hypotheses that suggest large standing bodies of water/ice in the northern lowlands in the planet's past history, Mars Orbiter Camera (MOC) images acquired to test these hypotheses have provided negative or ambiguous results. In the absence of classic coastal features to test the paleo-ocean hypothesis, other indicators need to be examined. Tuyas and hyaloclastic ridges are sub-ice volcanoes of unique appearance that form in ponded water conditions on Earth. Features with similar characteristics occur on Mars. MOLA analyses of these Martian features provide estimates of the height of putative ice/water columns at the edge of the Utopia Planitia basin and within Ophir Chasma of Valles Marineris, and support the hypotheses of a northern ocean on Mars.     

Impacts into porous foam targets: possible implications for the disruption of comet nuclei

We have conducted a series of impact experiments to examine the response of very porous foam targets to various impacts. Under near-vacuum conditions, closed-pore and open-pore foam targets were subjected to ∼1 km s⁻¹ impacts from aluminum and foam projectiles. We found that open-pore targets absorbed the impacts with little or no global fragmentation or noticeable cratering, exhibiting only local damage along the path of the projectile, which tunneled through the target. Closed-pore targets exhibited nearly explosive disruption, apparently resulting from stresses built up within the target due to internal pressurization from air that could not escape the target interior during evacuation of the impact chamber. These results suggest that build-up of internal volatile pressure within the nuclei of collisionally or dynamically unevolved comets could allow comparatively small impacts onto their surfaces to result in disproportionately disruptive outcomes.     

The tidal response of Ganymede and Callisto with and without liquid water oceans

Calculations of the tidal responses of Ganymede and Callisto reveal that tidal amplitudes on these bodies may be as large as a few meters if a liquid ocean exists to decouple the surface ice from the interior. Tides on Ganymede's surface can exceed 7 m peak-to-peak variation, while on Callisto the tidal amplitude can exceed 5 m in the presence of a liquid ocean. Without an ocean, tidal amplitudes are less than 0.5 m on Ganymede and less than 0.3 m on Callisto. An orbiting spacecraft using an altimeter for crossover analysis and Doppler tracking from Earth should be able to achieve sufficient accuracy to identify the tidal amplitude to within about a meter over the course of a few months (observing tens of tidal cycles).     

Subsurface oceans and deep interiors of medium-sized outer planet satellites and large trans-neptunian objects

The detection of induced magnetic fields in the vicinity of the jovian satellites Europa, Ganymede, and Callisto is one of the most surprising findings of the Galileo mission to Jupiter. The observed magnetic signature cannot be generated in solid ice or in silicate rock. It rather suggests the existence of electrically conducting reservoirs of liquid water beneath the satellites' outermost icy shells that may contain even more water than all terrestrial oceans combined. The maintenance of liquid water layers is closely related to the internal structure, composition, and thermal state of the corresponding satellite interior. In this study we investigate the possibility of subsurface oceans in the medium-sized icy satellites and the largest trans-neptunian objects (TNO's). Controlling parameters for subsurface ocean formation are the radiogenic heating rate of the silicate component and the effectiveness of the heat transfer to the surface. Furthermore, the melting temperature of ice will be significantly reduced by small amounts of salts and/or incorporated volatiles such as methane and ammonia that are highly abundant in the outer Solar System. Based on the assumption that the satellites are differentiated and using an equilibrium condition between the heat production rate in the rocky cores and the heat loss through the ice shell, we find that subsurface oceans are possible on Rhea, Titania, Oberon, Triton, and Pluto and on the largest TNO's 2003 UB₃₁₃, Sedna, and 2004 DW. Subsurface oceans can even exist if only small amounts of ammonia are available. The liquid subsurface reservoirs are located deeply underneath an ice-I shell of more than 100 km thickness. However, they may be indirectly detectable by their interaction with the surrounding magnetic fields and charged particles and by the magnitude of a satellite's response to tides exerted by the primary. The latter is strongly dependent on the occurrence of a subsurface ocean which provides greater flexibility to a satellite's rigid outer ice shell.     

Physics and Chemistry of Impacts

Earth, Moon, and Planets -     

Impacts into quartz sand: Crater formation,shock metamorphism,and ejecta distribution in laboratory experiments and numerical models

We investigated the ejection mechanics by a complementary approach of cratering experiments, including the microscopic analysis of material sampled from these experiments, and 2‐D numerical modeling of vertical impacts. The study is based on cratering experiments in quartz sand targets performed at the NASA Ames Vertical Gun Range. In these experiments, the preimpact location in the target and the final position of ejecta was determined by using color‐coded sand and a catcher system for the ejecta. The results were compared with numerical simulations of the cratering and ejection process to validate the iSALE shock physics code. In turn the models provide further details on the ejection velocities and angles. We quantify the general assumption that ejecta thickness decreases with distance according to a power‐law and that the relative proportion of shocked material in the ejecta increase with distance. We distinguish three types of shock metamorphic particles (1) melt particles, (2) shock lithified aggregates, and (3) shock‐comminuted grains. The agreement between experiment and model was excellent, which provides confidence that the models can predict ejection angles, velocities, and the degree of shock loading of material expelled from a crater accurately if impact parameters such as impact velocity, impactor size, and gravity are varied beyond the experimental limitations. This study is relevant for a quantitative assessment of impact gardening on planetary surfaces and the evolution of regolith layers on atmosphereless bodies.     

High-Resolution Simulations of the Impacts of Asteroids into the Venusian Atmosphere II: 3D Models

We compare high-resolution 2D and 3D numerical hydrocode simulations of asteroids striking the atmosphere of Venus. Our focus is on aerobraking and its effect on the size of impact craters. We consider impacts both by spheres and by the real asteroid 4769 Castalia, a severely nonspherical body in a Venus-crossing orbit. We compute mass and momentum fluxes as functions of altitude as global measures of the asteroid's progress. We find that, on average, the 2D and 3D simulations are in broad agreement over how quickly an asteroid slows down, but that the scatter about the average is much larger for the 2D models than for the 3D models. The 2D models appear to be rather strongly susceptible to the “butterfly effect,” in which tiny changes in initial conditions (e.g., 0.05% change in the impact velocity) produce quite different chaotic evolutions. By contrast, the global properties of the 3D models appear more reproducible despite seemingly large differences in initial conditions. We argue that this difference between 2D and 3D models has its root in the greater geometrical constraints present in any 2D model, and in particular in the global conservation of enstrophy in 2D that forces energy to pool in large-scale structures. It is the interaction of these artificial large-scale structures that causes slightly different 2D models to diverge so greatly. These constraints do not apply in 3D and large scale structures are not observed to form. A one-parameter modified pancake model reproduces the expected crater diameters of the 3D Castalias reasonably well.     

The axial angular momentum of the general circulation of the oceans

With the Hamburg geostrophic ocean general circulation model, the temporal variation of the angular momentum of the oceans is studied. The input wind stresses consist of climatological winds plus the deviations of the actual winds 1981–1989. The motion part of the angular momenta of the two hemispheres resembles the (counter-phase) behaviour of the corresponding atmospheric angular momentum. The dominating matter part however shows no phase shift between the hemispheres. The sum of the angular momenta of the solid Earth and the atmosphere is deduced from observations; in principle it should be inversely correlated to the oceanic angular momentum. A first attempt to recognize the correlation between these empirical data and the theoretical values is described.     

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.     

Detecting the glint of starlight on the oceans of distant planets

We propose that astronomers will be eventually be able to discriminate between extrasolar Earth-like planets with surface oceans and those without using the shape of phase light curves in the visible and near-IR spectrum. We model the visible light curves of planets having Earth-like surfaces, seasons, and optically-thin atmospheres with idealized diffuse-scattering clouds. We show that planets partially covered by water will appear measurably brighter near crescent phase (relative to Lambertian planets) because of the efficient specular reflection (“glint”) of starlight incident on their surfaces at a highly oblique angle. Planets on orbits within 30° of edge-on orientation (50% of all planets) will show pronounced glint over a sizeable range of orbital longitudes, from quadrature to crescent, all outside the glare of their parent stars. Also, water-covered planets will appear darker than a Lambertian disk near full illumination. Finally, we show that planets with a mixed land/water surface will polarize the reflected signal by as much as 30-70%. These results suggest several new ways of directly identifying water on distant planets.     

NEOs and Impacts: The Galactic Connection

The orbit of the solar system within the Galaxy is subject both to adiabatic tidal forces and to irregular buffeting caused by encounters with molecular clouds and spiral arms, while the long-period comet flux varies with the differential of these forces. Quasi-periodic episodes of Earth bombardment by Oort cloud comets are expected, with the half-period of the solar vertical oscillations being prominent. Various terrestrial datasets are examined and show evidence of a periodicity ∼27±1 Myr at a confidence level ≳ 99.5%, consistently with the expectations. Bombardment episodes thus appear to exert a controlling influence on global terrestrial processes. Conversely, the long-term terrestrial record provides an empirical constraint on the source regions of NEOs as well as on the dark matter content of the Galactic disc. During such episodes, climatic effects due to stratospheric dusting by very large comets are likely to be a major factor in geological and biological change. This revised version was published online in July 2006 with corrections to the Cover Date.     

Impacts into non-polar ice-rich paleodeposits on Mars: Excess ejecta craters, perched craters and pedestal craters as clues to Amazonian climate history

We compare three previously independently studied crater morphologies - excess ejecta craters, perched craters, and pedestal craters - each of which has been proposed to form from impacts into an ice-rich surface layer. Our analysis identifies the specific similarities and differences between the crater types; the commonalities provide significant evidence for a genetic relationship among the morphologies. We use new surveys of excess ejecta and perched craters in the southern hemisphere in conjunction with prior studies of all of the morphologies to create a comprehensive overview of their geographic distributions and physical characteristics. From these analyses, we conclude that excess ejecta craters and perched craters are likely to have formed from the same mechanism, with excess ejecta craters appearing fresh while perched craters have experienced post-impact modification and infilling. Impacts that led to these two morphologies overwhelmed the ice-rich layer, penetrating into the underlying martian regolith, resulting in the excavation of rock that formed the blocky ejecta necessary to armor the surface and preserve the ice-rich deposits. Pedestal craters, which tend to be smaller in diameter, have the same average deposit thickness as excess ejecta and perched craters, and form in the same geographic regions. They rarely have ejecta around their crater rims, instead exhibiting a smooth pedestal surface. We interpret this to mean that they form from impacts into the same type of ice-rich paleodeposit, but that they do not penetrate through the icy surface layer, and thus do not generate a blocky ejecta covering. Instead, a process related to the impact event appears to produce a thin, indurated surface lag deposit that serves to preserve the ice-rich material. These results provide a new basis to identify the presence of Amazonian non-polar ice-rich deposits, to map their distribution in space and time, and to assess Amazonian climate history. Specifically, the ages, distribution and physical attributes of the crater types suggest that tens to hundreds of meters of ice-rich material has been episodically emplaced at mid latitudes in both hemispheres throughout the Amazonian due to obliquity-driven climate variations. These deposits likely accumulated more frequently in the northern lowlands, resulting in a larger population of all three crater morphologies in the northern hemisphere.     

Radarclinometry of the sand seas of Africa’s Namibia and Saturn’s moon Titan

Radarclinometry is a powerful technique for estimating heights of landforms in synthetic aperture radar (SAR) images of planetary surfaces. In particular, it has been used to estimate heights of dunes in the sand seas of Saturn’s moon Titan (Lorenz, R.D., and 39 colleagues [2006]. Science 312, 724-727). In this work, we verify the technique by comparing dune heights derived from radarclinometry to known topography of dune fields in the Namib sand sea of western Africa. We compared results from three different image grid spacings, and found that 350 m/pixel (the same spacing at which the Cassini RADAR data was processed) is sufficient to determine dune height for dunes of similar morphometry to those of the Namib sand sea. At this grid spacing, height estimates derived from radarclinometry are largely representative of, though may underestimate by as much as 30%, or overestimate by as much as 40%, true dune height. Applying the technique to three regions on Titan, we estimate dune heights of 45-180 m, and dune spacings of 2.3-3.3 km. Obtaining accurate heights of Titan’s dunes will help to constrain the total organic inventory on Titan.     

Numerical Modeling of Marine Target Impacts

The principles of the numerical modeling of marine impacts of large cosmic bodies are described. Three underwater impact structures, MjØlnir, Lockne, and Eltanin, are considered with the aim of studying the characteristics of the crater formation at varying sea depths; the distinctions between the underwater and continental craters are discussed. The mechanisms for tsunami-wave generation are studied at different ratios of sea depth to impactor size. The calculation results are compared to the experimental data obtained during underwater nuclear explosions.