Modeling the morphological diversity of impact craters on icy satellites

Impact craters on icy satellites display a wide range of morphologies, some of which have no counterpart on rocky bodies. Numerical simulation studies have struggled to reproduce the diversity of features, such as central pits and transitions in crater depth with increasing diameter, observed on the icy Galilean satellites. The transitions in crater depth (at diameters of about 26 and 150 km on Ganymede and Callisto) have been interpreted as reflecting subsurface structure. Using the CTH shock physics code, we model the formation of craters with diameters between 400 m and about 200 km on Ganymede using different subsurface temperature profiles. Our calculations include recent improvements in the model equation of state for H₂O and quasi-static strength parameters for ice. We find that the shock-induced formation of dense high-pressure polymorphs (ices VI and VII) creates a gap in the crater excavation flow, which we call discontinuous excavation. For craters larger than about 20 km, discontinuous excavation concentrates a hot plug of material (>270 K and mostly on the melting curve) in the center of the crater floor. The size and occurrence of the hot plug are in good agreement with the observed characteristics of central pit craters, and we propose that a genetic link exists between them. We also derive depth versus diameter curves for different internal temperature profiles. In a 120 K isothermal crust, calculated craters larger than about 30 km diameter are deeper than observed and do not reproduce the transition at about 26 km diameter. Calculated crater depths are shallower and in good agreement with observations between about 30 and 150 km diameter using a warm thermal gradient representing a convective interior. Hence, the depth-to-diameter transition at about 26 km reflects thermal weakening of ice. Finally, simulation results generally support the hypothesis that the anomalous interior morphologies for craters larger than 100 km are related to the presence of a subsurface ocean.     

Experimental impact shock chemistry on planetary icy satellites

We present new experimental results on impact shock chemistry into icy satellites of the outer planets. Icy mixtures of pure water ice with CO₂, Na₂CO₃, CH₃OH, and CH₃OH/(NH₄)₂SO₄ at 77 K were ablated with a powerful pulsed laser—a new technique used to simulate shock processes which can occur during impacts. New products were identified by GC-MS and FTIR analyses after laser ablation. Our results show that hydrogen peroxide is formed in irradiated H₂O/CO₂ ices with a final concentration of 0.23%. CO and CH₃OH were also detected as main products. The laser ablation of frozen H₂O/Na₂CO₃ generates only CO and CO₂ as destruction products from the salt. Pulsed irradiation of water ice containing methanol leads also to the formation of CO and CO₂, generates methane and more complex molecules containing carbonyl groups like acetaldehyde, acetone, methyl formate, and a diether, dimethyl formal. The last three compounds are also produced when adding ammonium sulfate to H₂O/CH₃OH ice, but acetone is more abundant. The formation of two hydrocarbons, CH₄ and C₂H₆ is observed as well as the production of three nitrogen compounds, nitrous oxide, hydrogen cyanide, and acetonitrile.     

Diameter to depth dependence of impact craters

Laboratory impact experiments in the micron to millimeter projectile size range in silicate and metal targets have been performed in order to clarify the still ambigously interpreted velocity dependence of the crater diameter to depth ratios $\text{(}\text{D}\text{T}\text{)}$. The experimental results clearly show the independence of the $\text{D}\text{T}$ ratio of velocities above a threshold velocity of 3–4 km s^?1. The $\text{D}\text{T}$ ratio is a function of target properties and of projectile density ?. For a given target, the resulting approximate relation is $\text{D}\text{T}\text{～ ?}{}^{\mathrm{?\alpha }}$ with α varying between $\text{1}\text{2}\text{and}\text{1}\text{5}$.     

Viscous relaxation of craters within the martian south polar layered deposits

The South Polar Layered Deposits (PLD) are of fundamental importance to martian climatology, as they may comprise the largest reservoir of near-surface water on Mars. The South PLD exhibit relatively young crater retention surface ages, which are widely attributed to recent resurfacing. However, we show that both constructional and destructional resurfacing mechanisms (such as dust deposition and water ice sublimation, respectively) are inconsistent with the size, depth, and spatial distributions of South PLD craters. We demonstrate that another process—viscous creep relaxation of dusty water ice—is more compatible with the observed cratering of South PLD surfaces. The results of our finite element relaxation simulations suggest that, despite their apparent youthfulness, the PLD have been stable for at least several hundred million and perhaps even over a billion years. Consequently, our modeling implies that the time scales for the formation and preservation of the layers characteristic of the South (and possibly North) PLD are much longer than generally assumed.     

Shock remagnetization associated with meteorite impact at planetary surfaces

Shock remagnetization is a significant mode of alteration of the intensity and direction of magnetization in planetary crustal rocks subjected to the dynamic and thermochemical effects associated with meteorite impact. Shock remagnetization will take place almost instantaneously during and following the transient shock episode, and over longer times depending on residual temperature effects associated with shock heating and the production of impact melt. Remagnetization will follow certain demagnetization effects. The following transitions and residual effects will result in remagnetization of planetary crustal material:

1. First order reversible crystallographic transitions in bodycentered cubic iron-nickel alloys.
2. Second order Curie temperature transitions in face-centered cubic iron-nickel alloys.
3. Shock induced uniaxial anisotropy due to magnetoelasstic coupling of magnetic vectors to the shock wave.
4. Shock melting of iron containing silicates.
5. Subsolidus reduction and FeO decomposition.
6. Partial ther moremanence due to post-shock temperature.
7. Total thermoremanence due to post-shock temperature.
8. Production of a superparamagnetic distribution of iron which is sensitive to surface temperature fluctuation.
9. Thermal effects in metal and alloy phases.

Lunar breccia and soil samples are generally more reduced than crystalline rocks and some of th's reduction is subsolidus probably associated with the transient thermal effects due to meteorite impact in teh porous reglith.     

Phobos: Surface density of impact craters

The preliminary conclusion of the Mariner 9 Television Team that the surfaces of Phobos and Deimos are saturated with craters larger than 0.2 km in diameter is reconsidered on the basis of more extensive and uniform crater counts. For Phobos, it is verified that the surface appears saturated with craters larger than 1 km in diameter. For craters smaller than 1 km, the data points fall below the saturation curve, and it is not clear that all the departure can be explained in terms of loss of resolution. For Deimos, because of the paucity of craters visible in the Mariner 9 images, a statistically meaningful crater density curve cannot be constructed. Definitive crater density curves for subkilometer craters can only be established once additional imagery at a resolution of better than 100 m is obtained. Such imagery will be provided by the 1976 Viking Orbiters.     

Machine cataloging of impact craters on Mars

This study presents an automated system for cataloging impact craters using the MOLA 128 pixels/degree digital elevation model of Mars. Craters are detected by a two-step algorithm that first identifies round and symmetric topographic depressions as crater candidates and then selects craters using a machine-learning technique. The system is robust with respect to surface types; craters are identified with similar accuracy from all different types of martian surfaces without adjusting input parameters. By using a large training set in its final selection step, the system produces virtually no false detections. Finally, the system provides a seamless integration of crater detection with its characterization. Of particular interest is the ability of our algorithm to calculate crater depths. The system is described and its application is demonstrated on eight large sites representing all major types of martian surfaces. An evaluation of its performance and prospects for its utilization for global surveys are given by means of detailed comparison of obtained results to the manually-derived Catalog of Large Martian Impact Craters. We use the results from the test sites to construct local depth-diameter relationships based on a large number of craters. In general, obtained relationships are in agreement with what was inferred on the basis of manual measurements. However, we have found that, in Terra Cimmeria, the depth/diameter ratio has an abrupt decrease at ∼38°S regardless of crater size. If shallowing of craters is attributed to presence of sub-surface ice, a sudden change in its spatial distribution is suggested by our findings.     

Variation of albedo with solar incidence angle on planetary surfaces

While it has been known for a long time that in general the albedo of a surface depends on incidence angle, this fact is commonly neglected in many calculations of planetary surface temperatures. We show that the effect is especially pronounced for bright surfaces. For objects such as Ganymede and Io, the effect produces substantially cooler temperatures near the poles and terminators that would be calculated under the assumption of a constant albedo—a factor which may be important in determining the stability of frosts on such surfaces.     

Radar topography of domes on planetary surfaces

We investigate the possibility of measuring the heights and morphology of viscously emplaced domes using radar imagery. We accurately reproduce the known height and shape of a terrestrial salt dome, and estimate the heights of several venusian pancake domes to within a factor of two. The terrestrial salt dome is consistent with a Bingham flow, while the much larger venusian pancake domes are consistent with a Newtonian flow. Applying the same techniques to Ganesa Macula, a potential cryovolcanic dome on Titan, we estimate a height between 2.0-4.9 km. Additional factors such as variable roughness and composition might account for some of the discrepancies observed.     

Morphological characteristics of lunar craters with small depth/diameter ratio. II

The morphological characteristics of craters with relatively small ratio depth/diameter are discussed. It is observed that many morphological similarities exist among craters possessing ratiosd/D which do not differ considerably. The distribution of 1933 craters with respect to diameter, depth and the depth/diameter ratio is presented.     

Morphological characteristics of lunar craters with small depth/diameter ratio. I

The value of the ratio depth/diameter is examined as a criterion of morphological classification of lunar craters. Craters of small depth/diameter ratio are studied in detail.     

Modeling global impact effects on middle-sized icy bodies: applications to Saturn's moons

Disrupted terrains that form as a consequence of giant impacts may help constrain the internal structures of planets, asteroids, comets and satellites. As shock waves and powerful seismic stress waves propagate through a body, they interact with the internal structure in ways that may leave a characteristic impression upon the surface. Variations in peak surface velocity and tensile stress, related to landform degradation and surface rupture, may be controlled by variations in core size, shape and density. Caloris Basin on Mercury and Imbrium Basin on the Moon have disturbed terrain at their antipodes, where focusing is most intense for an approximately symmetric spheroid. Although, the icy saturnian satellites Tethys, Mimas, and Rhea possess giant impact structures, it is not clear whether these structures have correlated disrupted terrains, antipodal or elsewhere. In anticipation of high-resolution imagery from Cassini, we investigate antipodal focusing during giant impacts using a 3D SPH impact model. We first investigate giant impacts into a fiducial 1000 km diameter icy satellite with a variety of core radii and compositions. We find that antipodal disruption depends more on core radius than on core density, suggesting that core geometry may express a surface signature in global impacts on partially differentiated targets. We model Tethys, Mimas, and Rhea according to their image-derived shapes (triaxial for Tethys and Mimas and spherical for Rhea), varying core radii and densities to give the proper bulk densities. Tethys shows greater antipodal values of peak surface velocity and peak surface tensile stress, indicating more surface damage, than either Mimas or Rhea. Results for antipodal and global fragmentation and terrain rupture are inconclusive, with the hydrocode itself producing global disruption at the limits of model resolution but with peak fracture stresses never exceeding the strength of laboratory ice.     

Geology of large impact craters on Ganymede: Implications on thermal and tectonic histories

The geometry of the furrows of Galileo Regio indicates that they are not of impact origin, and irrelevant to discussion about large impact effects. The detailed study of three large impact basins indicates that their transient cavity radii are different from previously reported values. Because of the relations between crater's size and lithospheric thickness, these new values of basins radii would constrain further models of Ganymede's thermal evolution. The geometry of lineaments around these three basins, which occurred on grooved terrains, indicates that these impacts induced tectonic motions along a preexisting planetary wide grid pattern. This pattern influenced also the formation of the furrows on Galileo Regio. That would indicate that the grooved terrains are only superficial layers and that they were formed without destruction or rotation of their basement.     

Venusian impact craters on Magellan images: View from Venera 15/16

The review of the data of the analysis of images of Venera 15/16 mission is done. These data are used to make some predictions of how many impact craters the MAGELLAN space probe would see on Venus and how they would look. The possible trends of crater areal distribution and characteristic features of the size-frequency relationship are discussed.'Geology and Tectonics of Venus', special issue edited by Alexander T. Basilevsky (USSR Acad. of Sci. Moscow), James W. Head (Brown University, Providence), Gordon H. Pettengill (MIT, Cambridge, Massachusetts) and R. S. Saunders (J.P.L., Pasadena).     

Transitional impact craters on the Moon: Insight into the effect of target lithology on the impact cratering process

We studied a data set of 28 well‐preserved lunar craters in the transitional (simple‐to‐complex) regime with the aim of investigating the underlying cause(s) for morphological differences of these craters in mare versus highland terrains. These transitional craters range from 15 to 42 km in diameter, demonstrating that the transition from simple to complex craters is not abrupt and occurs over a broad diameter range. We examined and measured the following crater attributes: depth (d), diameter (D), floor diameter (D_f), rim height (h), and wall width (w), as well as the number and onset of terraces and rock slides. The number of terraces increases with increasing crater size and, in general, mare craters possess more terraces than highland craters of the same diameter. There are also clear differences in the d/D ratio of mare versus highland craters, with transitional craters in mare targets being noticeably shallower than similarly sized highland craters. We propose that layering in mare targets is a major driver for these differences. Layering provides pre‐existing planes of weakness that facilitate crater collapse, thus explaining the overall shallower depths of mare craters and the onset of crater collapse (i.e., the transition from simple to complex crater morphology) at smaller diameters as compared to highland craters. This suggests that layering and its interplay with target strength and porosity may play a more significant role than previously considered.     

Displaced mass,depth, diameter,and effects of oblique trajectories for impact craters formed in dense crystalline rocks

Empirical formulae are presented for calculating the displaced mass, depth, diameter, and effects of oblique trajectories for impact craters formed in dense crystalline rocks. The formulae are applicable to craters with diameters from approximately 10^?3–10³ cm that require, respectively, impact kinetic energies of approximately 10 to 10¹⁶ ergs for their formation. The experimental results are in poor agreement with Öpik’s theoretical calculations and raise questions on the validity of his theoretical model.     

Radial transpression ridges: A new structural feature of complex impact craters

Abstract— The gravity‐driven collapse of complex impact craters induces the mass transfer of large rock volumes. In distal parts of a crater, inward movements dominate; whereas in the centre of an impact structure, outward movements occur if the central uplift collapses. The particle trajectory field is centro‐symmetric, which signifies that the conditions for plane‐strain deformation are not fullfilled. Converging particle trajectories can be compensated either by a bulk thickening of inward sliding masses (folding, repetition of rock units along thrust faults, plastic flow) or by the formation of localised radial transpression ridges (RTR) at the edges of individual landslides. In these transpression ridges, material is uplifted to accommodate the converging mass flow. Different modes of uplift are possible including radial folding, lateral overthrusting, and the formation of positive flower structures. A simple geometric model is used to estimate the amount of transpression thickening and bulk thickening of inward sliding masses on the base of volumetric considerations. The existence of RTRs is confirmed by structural investigations at the Siljan impact structure, Sweden, as well as at other complex impact craters on Earth.     

Photometry of rough planetary surfaces: The role of multiple scattering

Topographic features affect the scattering properties of planetary surfaces by casting shadows and altering the local incidence and emission angles. Measurements of this phenomenon were obtained on the Cornell goniometer for both high and low albedo surfaces. For the low albedo surface, the decrease in reflected radiation due to topography increases sharply with increasing phase angle, whereas for the high albedo sample the effects are approximately constant between phase angles of 30 and 70°. The observations are in good agreement with a theoretical model in the case of the dark surface. However, for the high albedo surface the model overestimates the effects by about a factor of 2, since it does not include the partial illumination of shadows by multiple scattering. For both high and low albedo surfaces, the effects of topography do not become significant until a phase angle of 30–40°.     

Shock metamorphism and impact melting in small impact craters on Earth: Evidence from Kamil crater,Egypt

Kamil is a 45 m diameter impact crater identified in 2008 in southern Egypt. It was generated by the hypervelocity impact of the Gebel Kamil iron meteorite on a sedimentary target, namely layered sandstones with subhorizontal bedding. We have carried out a petrographic study of samples from the crater wall and ejecta deposits collected during our first geophysical campaign (February 2010) in order to investigate shock effects recorded in these rocks. Ejecta samples reveal a wide range of shock features common in quartz‐rich target rocks. They have been divided into two categories, as a function of their abundance at thin section scale: (1) pervasive shock features (the most abundant), including fracturing, planar deformation features, and impact melt lapilli and bombs, and (2) localized shock features (the least abundant) including high‐pressure phases and localized impact melting in the form of intergranular melt, melt veins, and melt films in shatter cones. In particular, Kamil crater is the smallest impact crater where shatter cones, coesite, stishovite, diamond, and melt veins have been reported. Based on experimental calibrations reported in the literature, pervasive shock features suggest that the maximum shock pressure was between 30 and 60 GPa. Using the planar impact approximation, we calculate a vertical component of the impact velocity of at least 3.5 km s^?1. The wide range of shock features and their freshness make Kamil a natural laboratory for studying impact cratering and shock deformation processes in small impact structures.     

The rayed crater Zunil and interpretations of small impact craters on Mars

A 10-km diameter crater named Zunil in the Cerberus Plains of Mars created ∼¹⁰⁷ secondary craters 10 to 200 m in diameter. Many of these secondary craters are concentrated in radial streaks that extend up to 1600 km from the primary crater, identical to lunar rays. Most of the larger Zunil secondaries are distinctive in both visible and thermal infrared imaging. MOC images of the secondary craters show sharp rims and bright ejecta and rays, but the craters are shallow and often noncircular, as expected for relatively low-velocity impacts. About 80% of the impact craters superimposed over the youngest surfaces in the Cerberus Plains, such as Athabasca Valles, have the distinctive characteristics of Zunil secondaries. We have not identified any other large (?10 km diameter) impact crater on Mars with such distinctive rays of young secondary craters, so the age of the crater may be less than a few Ma. Zunil formed in the apparently youngest (least cratered) large-scale lava plains on Mars, and may be an excellent example of how spallation of a competent surface layer can produce high-velocity ejecta (Melosh, 1984, Impact ejection, spallation, and the origin of meteorites, Icarus 59, 234-260). It could be the source crater for some of the basaltic shergottites, consistent with their crystallization and ejection ages, composition, and the fact that Zunil produced abundant high-velocity ejecta fragments. A 3D hydrodynamic simulation of the impact event produced 10¹⁰ rock fragments ?10 cm diameter, leading to up to 10⁹ secondary craters ?10 m diameter. Nearly all of the simulated secondary craters larger than 50 m are within 800 km of the impact site but the more abundant smaller (10-50 m) craters extend out to 3500 km. If Zunil is representative of large impact events on Mars, then secondaries should be more abundant than primaries at diameters a factor of ∼1000 smaller than that of the largest primary crater that contributed secondaries. As a result, most small craters on Mars could be secondaries. Depth/diameter ratios of 1300 small craters (10-500 m diameter) in Isidis Planitia and Gusev crater have a mean value of 0.08; the freshest of these craters give a ratio of 0.11, identical to that of fresh secondary craters on the Moon (Pike and Wilhelms, 1978, Secondary-impact craters on the Moon: topographic form and geologic process, Lunar Planet. Sci. IX, 907-909) and significantly less than the value of ∼0.2 or more expected for fresh primary craters of this size range. Several observations suggest that the production functions of Hartmann and Neukum (2001, Cratering chronology and the evolution of Mars, Space Sci. Rev. 96, 165-194) predict too many primary craters smaller than a few hundred meters in diameter. Fewer small, high-velocity impacts may explain why there appears to be little impact regolith over Amazonian terrains. Martian terrains dated by small craters could be older than reported in recent publications.