Distribution Of Impact Craters Research Articles

Spherical multi-lacunarity reveals possible equatorial-polar differences in crater densities on the dwarf planet Ceres

The dwarf planet Ceres possesses a peculiar distribution of impact craters. It has been previously noted that the largest craters expected for Ceres are absent, while for smaller craters, the north polar region is the most heavily cratered. It thus appears that some process(es) have erased some of Ceres' craters, and the distribution of observed craters could point to the nature of these processes. For instance, a process tied to sunlight (e.g., relaxation or sublimation) could impart a latitudinal dependence, while a more regional distribution could point to a more endogenic process (e.g., cryovolcanism). Through a combination of spherical multi-lacunarity analysis and kernel density maps, we discover that the 20 largest craters (> 100 km across), though statistically indistinguishable from random, are stochastically concentrated in the south, while mid-size craters (20–70 km) show clustering at spatial scales >130 km, roughly the same scale that would be affected by emplacement of the 20 largest craters plus ejecta. A mask made from these largest craters plus ejecta reveals that a south polar region that matches the heavily cratered north is almost completely covered by the largest craters. Crater counts additionally reveal that this north polar region is older (i.e., more heavily cratered) than an equatorial region not masked by the largest craters, which in turn is about as cratered as the oldest of the 20 largest craters. Thus, it appears that there could be equator-polar differences to Ceres' crater distribution, or at the very least, that a latitudinal dependence cannot be discounted. Consequently, latitudinal variations in sunlight might be a controlling factor in the distribution of craters on Ceres.

The Effects of Terrain Properties Upon the Small Crater Population Distribution at Giordano Bruno: Implications for Lunar Chronology.

The distribution of impact craters on the ejecta of Giordano Bruno, a recent (<10 Ma) 22‐km diameter crater within the lunar highlands, exhibits substantial variations. We surveyed craters D ≥ 10 m across a 1,323 km2 area of Giordano Bruno's ejecta and compared the distribution of craters with variations in thermophysical properties derived from the Lunar Reconnaissance Orbiter Diviner instrument. We used Diviner‐derived rock abundance and nighttime regolith temperatures along with thermal model‐predicted surface temperatures for a diversity of terrains to identify and isolate areas of the ejecta based on thermophysical properties such as bulk density and thermal conductivity. We found that thermophysical properties of the ejecta vary considerably both laterally and vertically, and consistently differ from typical regolith, indicating the presence of higher thermal inertia materials. Crater‐size frequencies are significantly lower in areas with terrain properties exhibiting higher: rock abundance, nighttime temperatures, and/or modeled thermal inertia. This discrepancy in crater distribution increases for craters smaller than ∼25 m. These thermophysical variations indicate changes in the mechanical properties of the target materials. We suggest that these variations—specifically, terrain‐dependent crater scaling variations and impactor‐scale heterogeneities in material properties such as the presence or absence of large boulders—may influence crater diameters or inhibit crater production altogether in Giordano Bruno's ejecta; furthermore, these factors are size‐dependent.

Impact crater degradation, Oxia Planum, Mars

ABSTRACT The main goal of the European Space Agency (ESA) and Roscosmos ExoMars rover mission is to collect samples from the near-subsurface of Mars. The rover will look for any physical or chemical evidence of ancient life in the near subsurface. This map shows the distribution of impact craters at this proposed landing site in Oxia Planum on Mars. The map records 1199 impact craters > 500 m in diameter in a 5.0° × 2.5° region around Oxia Planum. The impact craters are symbolised based on the way different aspects of their morphology have degraded since their formation. The distribution of degradation and burial morphologies of impact craters can be used to determine where burial and erosion processes have occurred. Because the formation of impact craters is well constrained, occurs instantly and with a predictable flux, future studies could use this knowledge and our dataset to constrain when these events occurred.

Three-Dimensional Model of the Moon with Semantic Information of Craters Based on Chang'e Data.

China’s Chang’e lunar exploration project obtains digital orthophoto image (DOM) and digital elevation model (DEM) data covering the whole Moon, which are critical to lunar research. The DOM data have three resolutions (i.e., 7, 20 and 50 m), while the DEM has two resolutions (i.e., 20 and 50 m). Analysis and research on these image data effectively help humans to understand the Moon. In addition, impact craters are considered the most basic feature of the Moon’s surface. Statistics regarding the size and distribution of impact craters are essential for lunar geology. In existing works, however, the lunar surface has been reconstructed less accurately, and there is insufficient semantic information regarding the craters. In order to build a three-dimensional (3D) model of the Moon with crater information using Chang‘e data in the Chang‘e reference frame, we propose a four-step framework. First, software is implemented to annotate the lunar impact craters from Chang’e data by complying with our existing study on an auxiliary annotation method and open-source software LabelMe. Second, auxiliary annotation software is adopted to annotate six segments in the Chang’e data for an overall 25,250 impact crater targets. The existing but inaccurate craters are combined with our labeled data to generate a larger dataset of craters. This data set is analyzed and compared with the common detection data. Third, deep learning detection methods are employed to detect impact craters. To address the problem attributed to the resolution of Chang’e data being too high, a quadtree decomposition is conducted. Lastly, a geographic information system is used to map the DEM data to 3D space and annotate the semantic information of the impact craters. In brief, a 3D model of the Moon with crater information is implemented based on Chang’e data in the Chang‘e reference frame, which is of high significance.

The Physics of Changing Tectonic Regimes: Implications for the Temporal Evolution of Mantle Convection and the Thermal History of Venus

AbstractGiven similar sizes and compositions, Venus invites comparisons with Earth. While Earth convects in the plate tectonic regime, Venus' current and past tectonic states remain uncertain. Venus' impact crater distribution has led to competing models for its tectonic state, steady (e.g., stagnant lid) or catastrophic (e.g., episodic lid). Terrestrial geochemical evidence and geodynamic models suggest that planets may transition between tectonic states over time. Transitions can be triggered by increasing yield stress due to loss of water or by increasing surface temperature. Here, we revisit the classic endmember models of Venusian tectonics by modeling the transition from stable mobile lid convection into a stable stagnant lid in 3D spherical geometry. Transitions between states are governed by both regional and global scale instabilities, resulting in oscillations in heat flux, velocity, and magma production over Gyr time scales. The thermal and tectonic evolution of a convectively transitioning planet is neither catastrophic nor steady. Volcanism and tectonic yielding may be non‐global in nature. At any given time, portions of the planetary surface may reflect apparently different styles of convection, with some regions being highly active and other regions being sluggish and inactive. For Venus, the implications are that global melt production and resurfacing are a natural consequence of lid‐state evolution, without needing ad hoc resurfacing mechanisms which rely on a single governing lid‐state. Geologic evidence suggests that the transition from mobile lid to stagnant lid is still on‐going. Venus may have lost liquid surface water in the last third of its history.

The spatial distribution of impact craters on Ryugu

Asteroid 162173 Ryugu has numerous craters. The initial measurement of impact craters on Ryugu, by Sugita et al. (2019), is based on Hayabusa2 ONC images obtained during the first month after the arrival of Hayabusa2 in June 2018. Utilizing new images taken until February 2019, we constructed a global impact crater catalogue of Ryugu, which includes all craters larger than 20 m in diameter on the surface of Ryugu. As a result, we identified 77 craters on the surface of Ryugu. Ryugu shows variation in crater density which cannot be explained by the randomness of cratering; there are more craters at lower latitudes and fewer at higher latitudes, and fewer craters in the western bulge (160°E – 290°E) than in the region around the meridian (300°E – 30°E). This variation implies a complicated geologic history for Ryugu. It seems that the variation in crater density indicates that the equatorial ridge located in the western hemisphere is relatively young, while that located in the eastern hemisphere is a fossil structure formed during the short rotational period in the distant past.

Regolith mixing by impacts: Lateral diffusion of basin melt

Impact cratering has been the primary process to alter the distribution of lunar highland material since the formation of a crust. This impact history is recorded in the radiogenic clocks of impact melts which are accessible for study on lunar samples and meteorites. However, primary impact melt is exposed to a long-time gardening process (i.e. re-melting, excavation, burial, and re-excavation) by subsequent impacts resulting in a complex spatial distribution of materials representing specific impact events. To investigate the diffusion behavior of impact melt, a model tracing the evolving distribution of melt laterally and with depth was built using a Monte Carlo approach. Given scaling laws concerning melt production and ejecta distribution, the size-frequency distribution of impact craters, and the rate function for crater formation, we examine the evolution of melt component occurrence of different ages. Three mid- to late-forming basins (Serenitatis, Crisium, and Imbrium) are chosen as a case study for the diffusion of melt from major basin-forming events. The survival probability of basin melt occurrence at the Apollo and Luna sampling spots is derived. It is expected to find abundant Imbrium and Crisium melt at the Apollo and Luna sampling sites, consistent with the K-Ar radiometric dates of highland samples; whereas the older Serenitatis melt was subjected to the later long-term gardening, strongly influenced by later local impacts, and thus is less abundant. Understanding the diffusion of impact melt is helpful for interpretation of radiometric ages of lunar samples and can be used to predict the distribution of differently-aged melts at future landing/sampling sites such as the Chinese Chang'E-4 (CE-4).

Stereoscopic imaging determines space debris impact crater distribution and morphology

Stereoscopic imaging determines space debris impact crater distribution and morphology

高分解能衛星画像を用いた浅間山2004年9月1日噴火の噴石着弾痕分布図＜特集＞2004年浅間山噴火（3）

高分解能衛星画像を用いた浅間山2004年9月1日噴火の噴石着弾痕分布図＜特集＞2004年浅間山噴火（3）

Spatial distribution of impact craters on Deimos

Deimos, one of the Martian moons, has numerous impact craters. However, it is unclear whether crater saturation has been reached on this satellite. To address this issue, we apply a statistical test known as nearest-neighbor analysis to analyze the crater distribution of Deimos. When a planetary surface such as the Moon is saturated with impact craters, the spatial distribution of craters is generally changed from random to more ordered. We measured impact craters on Deimos from Viking and HiRISE images and found (1) that the power law of the size-frequency distribution of the craters is approximately −1.7, which is significantly shallower than those of potential impactors, and (2) that the spatial distribution of craters over 30m in diameter cannot be statistically distinguished from completely random distribution, which indicates that the surface of Deimos is inconsistent with a surface saturated with impact craters. Although a crater size-frequency distribution curve with a slope of −2 is generally interpreted as indicating saturation equilibrium, it is here proposed that two competing mechanisms, seismic shaking and ejecta emplacement, have played a major role in erasing craters on Deimos and are therefore responsible for the shallow slope of this curve. The observed crater density may have reached steady state owing to the obliterations induced by the two competing mechanisms. Such an occurrence indicates that the surface is saturated with impact craters despite the random distribution of craters on Deimos. Therefore, this work proposes that the age determined by the current craters on Deimos reflects neither the age of Deimos itself nor that of the formation of the large concavity centered at its south pole because craters should be removed by later impacts. However, a few of the largest craters on Deimos may be indicative of the age of the south pole event.

High-pressure minerals in eucrite suggest a small source crater on Vesta.

High-pressure minerals in meteorites are important records of shock events that have affected the surfaces of planets and asteroids. A widespread distribution of impact craters has been observed on the Vestan surface. However, very few high-pressure minerals have been discovered in Howardite-Eucrite-Diogenite (HED) meteorites. Here we present the first evidence of tissintite, vacancy-rich clinopyroxene, and super-silicic garnet in the eucrite Northwest Africa (NWA) 8003. Combined with coesite and stishovite, the presence of these high-pressure minerals and their chemical compositions reveal that solidification of melt veins in NWA 8003 began at a pressure of >~10 GPa and ceased when the pressure dropped to <~8.5 GPa. The shock temperature in the melt veins exceeded 1900 °C. Simulation results show that shock events that create impact craters of ~3 km in diameter (subject to a factor of 2 uncertainty) are associated with sufficiently high pressures to account for the occurrence of the high-pressure minerals observed in NWA 8003. This indicates that HED meteorites containing similar high-pressure minerals should be observed more frequently than previously thought.

Ejecta velocity distribution of impact craters formed on quartz sand: Effect of projectile density on crater scaling law

In order to clarify the effects of projectile density on ejecta velocity distributions for a granular target, impact cratering experiments on a quartz sand target were conducted by using eight types of projectiles with different densities ranging from 11gcm−3 to 1.1gcm−3, which were launched at about 200ms−1 from a vertical gas gun at Kobe University. The scaling law of crater size, the ejection angle of ejecta grains, and the angle of the ejecta curtain were also investigated. The ejecta velocity distribution obtained from each projectile was well described by the π-scaling theory of v0gR=k2x0R-1μ, where v0, g, R and x0 are the ejection velocity, gravitational acceleration, crater radius and ejection position, respectively, and k2 and μ are constants mostly depending on target material properties (Housen, K.R., Holsapple, K.A. [2011]. Icarus 211, 856–875). The value of k2 was found to be almost constant at 0.7 for all projectiles except for the nylon projectile, while μ increased with the projectile density, from 0.43 for the low-density projectile to 0.6–0.7 for the high-density projectile. On the other hand, the π-scaling theory for crater size gave a μ value of 0.57, which was close to the average of the μ values obtained from ejecta velocity distributions. The ejection angle, θ, of each grain decreased slightly with distance, from higher than 45° near the impact point to 30–40° at 0.6 R. The ejecta curtain angle is controlled by the two elementary processes of ejecta velocity distribution and ejection angle; it gradually increased from 52° to 63° with the increase of the projectile density. The comparison of our experimental results with the theoretical model of the crater excavation flow known as the Z-model revealed that the relationship between μ and θ obtained by our experiments could not be described by the Z-model (Maxwell, D.E. [1977]. In: Roddy, D.J., Pepin, R.O., Merrill, R.B. (Eds.), Impact and Explosion Cratering. Pergamon, NY, pp. 1003–1008). Therefore, we used the extended Z-model by Croft (Croft, S.K. [1980]. Proc. Lunar Sci. Conf. 11, 2347–2378), which could be applied to the crater excavation process when the point source was buried at the depth of d under the target surface, and then all the experimental results of μ and θ were reasonably explained by suitable Z and d values of the extended Z-model.

Rheological and petrological implications for a stagnant lid regime on Venus

Venus is physically similar to Earth but with no oceans and a hot dense atmosphere. Its near-random distribution of impact craters led to the inferences of episodic global resurfacing and a stagnant lid regime, and imply that it is not currently able to lose proportionately as much heat as Earth. This paper shows that a CO2-induced asthenosphere and decoupling of the mantle lid from the crust, caused by the elevated surface temperature, enables lid rejuvenation. Global hypsography implies a rate of 4·0±0·5km²a−1 and an implied heat loss rate of 32·8±3·6TW, ~90% of a scaled Earth-like rate of heat loss of 36TW. Estimates of the rate of lid rejuvenation by plume activity – 0·07 to 0·09km²a−1– imply that ten times the number of observed plumes are required to equal this rate of heat loss. However, lid rejuvenation by convection allows Venus to maintain a stable tectonic regime, with subcrustal horizontal extension (half-spreading) rates of between 25 and 50mma−1 determined from fits to topographic profiles across the principal rift systems. While the surface is largely detached from these processes, the association of rifting and other processes does imply that Venus is geologically active at the present day.

Watershed modeling in the Tyrrhena Terra region of Mars

Watershed analyses from high‐resolution image (Viking, Mars Orbiter Camera, and Thermal Emission Imaging System) and topographic (Mars Orbiter Laser Altimeter [MOLA]) data are used to qualitatively and quantitatively characterize highland fluvial systems and analyze the role of water in the evolution of Tyrrhena Terra (13°S–30°S, 265°W–280°W), Mars. In this study, Geographical Information System software is used in conjunction with MOLA Digital Elevation Models to delineate drainage basin divides, extract valley networks, and derive basin and network morphometric parameters (e.g., drainage density, stream order, bifurcation ratio, and relief morphometry) useful in characterizing the geologic and climatic conditions of watershed formation, as well as for evaluating basin “maturity” and processes of watershed development. Model‐predicted valley networks and watershed boundaries, which are dependent on the degree to which pixel sinks are filled in the topographic data set and a channelization threshold, are evaluated against image and topographic data, slope maps, and detailed maps of valley segments from photogeologic analyses. Valley morphologies, crater/valley relationships, and impact crater distributions show that valleys in Tyrrhena Terra are ancient. Based on geologic properties of the incised materials, valley and network morphologies, morphometric parameters, and the presence of many gullies heading at or near‐crater rim crests, surface runoff, derived from rainfall or snowmelt, was the dominant erosional process; sapping may have only played a secondary role in valley formation in Tyrrhena Terra. Furthermore, spatial and temporal relationships of dissected highland materials and impact craters, suggests widespread, but relatively short‐lived, erosion by runoff with most activity in the Noachian period.

Impact cratering records of the mid-sized, icy saturnian satellites

Resolution of Voyager 1 and 2 images of the mid-sized, icy saturnian satellites was generally not much better than 1 km per line pair, except for a few, isolated higher resolution images. Therefore, analyses of impact crater distributions were generally limited to diameters ( D) of tens of kilometers. Even with the limitation, however, these analyses demonstrated that studying impact crater distributions could expand understanding of the geology of the saturnian satellites and impact cratering in the outer Solar System. Thus to gain further insight into Saturn’s mid-sized satellites and impact cratering in the outer Solar System, we have compiled cratering records of these satellites using higher resolution Cassini ISS images. Images from Cassini of the satellites range in resolution from tens m/pixel to hundreds m/pixel. These high-resolution images provide a look at the impact cratering records of these satellites never seen before, expanding the observable craters down to diameters of hundreds of meters. The diameters and locations of all observable craters are recorded for regions of Mimas, Tethys, Dione, Rhea, Iapetus, and Phoebe. These impact crater data are then analyzed and compared using cumulative, differential and relative ( R) size-frequency distributions. Results indicate that the heavily cratered terrains on Rhea and Iapetus have similar distributions implying one common impactor population bombarded these two satellites. The distributions for Mimas and Dione, however, are different from Rhea and Iapetus, but are similar to one another, possibly implying another impactor population common to those two satellites. The difference between these two populations is a relative increase of craters with diameters between 10 and 30 km and a relative deficiency of craters with diameters between 30 and 80 km for Mimas and Dione compared with Rhea and Iapetus. This may support the result from Voyager images of two distinct impactor populations. One population was suggested to have a greater number of large impactors, most likely heliocentric comets (Saturn Population I in the Voyager literature), and the other a relative deficiency of large impactors and a greater number of small impactors, most likely planetocentric debris (Saturn Population II). Meanwhile, Tethys’ impact crater size-frequency distribution, which has some similarity to the distributions of Mimas, Dione, Rhea, and Iapetus, may be transitional between the two populations. Furthermore, when the impact crater distributions from these older cratered terrains are compared to younger ones like Dione’s smooth plains, the distributions have some similarities and differences. Therefore, it is uncertain whether the size-frequency distribution of the impactor population(s) changed over time. Finally, we find that Phoebe has a unique impact crater distribution. Phoebe appears to be lacking craters in a narrow diameter range around 1 km. The explanation for this confined “dip” at D = 1 km is not yet clear, but may have something to do with the interaction of Saturn’s irregular satellites or the capture of Phoebe.

Caloris impact basin: Exterior geomorphology, stratigraphy, morphometry, radial sculpture, and smooth plains deposits

The Mariner 10 flybys revealed the eastern third of the Caloris basin and demonstrated its importance as a chronostratigraphic marker for the geological history of Mercury. Uncertain after that mission were the size of the basin, the full distribution of materials deposited or modified during the impact process, and the stratigraphic and age relationships between Caloris and surrounding terrain. Images obtained during MESSENGER's January 2008 flyby of Mercury revealed that the Caloris basin is ∼ 15% larger than previously estimated and moderately elliptical (∼ 1525 by ∼ 1315 km). Basin-related sculpture and secondary craters are dispersed widely in areas surrounding the basin, confirming the widespread significance of this event as a stratigraphic marker. Units mapped around Caloris on the basis of Mariner 10 data are generally recognized in regions observed for the first time by MESSENGER, including most facies of the Caloris Group (the Caloris Montes, Odin Formation, and Van Eyck Formation). The only unit without obvious exposures in MESSENGER images is the Nervo Formation, which was originally interpreted as impact melt or fallback ejecta. New measurements of the size-frequency distribution of impact craters both within the Caloris basin and on surrounding units, particularly the Odin Formation, support the hypothesis that the vast majority of circum-Caloris plains postdate the basin and are likely volcanic. The MESSENGER observations re-emphasize the importance of understanding the exterior stratigraphy of the basin in order to utilize the Caloris event as a discrete time horizon in the geological history of Mercury.

Crater modification and geologic activity in Enceladus' heavily cratered plains: Evidence from the impact crater distribution

Although we can observe current activity on Saturn's satellite Enceladus with Cassini, insight into past activity is best achieved (for now) through studying the impact crater distributions. Furthermore, approximation of terrain ages can only be attained through calculations using crater densities and estimations of impact rates in the saturnian system. Here we focus on what the impact crater distribution in Enceladus' heavily cratered plains can tell us about Enceladus' geologic history. We use Cassini ISS images to count craters in the heavily cratered plains on Enceladus, along with Rhea, Dione, Tethys and Mimas as references, to develop and compare their size-frequency distributions. Comparisons of our counts show that Enceladus' cratered plains distribution is unique in that it appears to have a relative deficiency of craters for diameters ⩽ 2 km and ⩾ 6 km compared to the other satellites' heavily cratered plains. Our data also indicates that the impact crater density within the cratered plains changes with latitude. Specifically, both the north and south mid-latitude regions have approximately three times higher density than the equatorial region. We hypothesize that the “missing” small and large craters in Enceladus' cratered plains is due to a combination of viscous relaxation of the larger craters, and burial of the relaxed large craters and small craters by south polar plume and possibly E-ring material. We also conclude that the spatial density distribution is not consistent with recent polar wander.

On the negligible surface age of Triton

New mapping reveals 100 probable impact craters on Triton wider than 5 km diameter. All of the probable craters are within 90° of the apex of Triton's orbital motion (i.e., all are on the leading hemisphere) and have a cosine density distribution with respect to the apex. This spatial distribution is difficult to reconcile with a heliocentric (Sun-orbiting) source of impactors, be it ecliptic comets, the Kuiper Belt, the scattered disk, or tidally-disrupted temporary satellites in the style of Shoemaker–Levy 9, but it is consistent with head-on collisions, as would be produced if a prograde population of planetocentric (Neptune-orbiting) debris were swept up by retrograde Triton. Plausible sources include ejecta from impact on or disruption of inner/outer moons of Neptune. If Triton's small craters are mostly of planetocentric origin, Triton offers no evidence for or against the existence of small comets in the Kuiper Belt, and New Horizons observations of Pluto must fill this role. The possibility that the distribution of impact craters is an artifact caused by difficulty in identifying impact craters on the cantaloupe terrain is considered and rejected. The possibility that capricious resurfacing has mimicked the effect of head-on collisions is considered and shown to be unlikely given current geologic constraints, and is no more probable than planetocentrogenesis. The estimated cratering rate on Triton by ecliptic comets is used to put an upper limit of ∼50 Myr on the age of the more heavily cratered terrains, and of ∼6 Myr for the Neptune-facing cantaloupe terrain. If the vast majority of cratering is by planetocentric debris, as we propose, then the surface everywhere is probably less than 10 Myr old. Although the uncertainty in these cratering ages is at least a factor ten, it seems likely that Triton's is among the youngest surfaces in the Solar System, a candidate ocean moon, and an important target for future exploration.

Diapir-induced reorientation of Saturn's moon Enceladus

Enceladus is a small icy satellite of Saturn. Its south polar region consists of young, tectonically deformed terrain and has an anomalously high heat flux. This heat flux is probably due to localized tidal dissipation within either the ice shell or the underlying silicate core. The surface deformation is plausibly due to upwelling of low-density material (diapirism) as a result of this tidal heating. Here we show that the current polar location of the hotspot can be explained by reorientation of the satellite's rotation axis because of the presence of a low-density diapir. If the diapir is in the ice shell, then the shell must be relatively thick and maintain significant rigidity (elastic thickness greater than approximately 0.5 km); if the diapir is in the silicate core, then Enceladus cannot possess a global subsurface ocean, because the core must be coupled to the overlying ice for reorientation to occur. The reorientation generates large (approximately 10 MPa) tectonic stress patterns that are compatible with the observed deformation of the south polar region. We predict that the distribution of impact craters on the surface will not show the usual leading hemisphere-trailing hemisphere asymmetry. A low-density diapir also yields a potentially observable negative gravity anomaly.

Polar roller

The recent discovery of active water-venting from the surface of Saturn's icy moon Enceladus has caused a ripple of excitement as it suggests the possibility of subsurface liquid water. Geologically, Enceladus poses problems. One, the reason why the hotspot of geological activity is centred precisely on the south pole, has been tackled by Francis Nimmo and Robert T. Pappalardo. They identify reorientation of Enceladus's rotation axis as the likely cause, and an upwelling of low-density material to the surface as the cause of that reorientation. Comparison of the distribution of impact craters on the ‘leading’ and ‘trailing’ hemispheres will be a good test of the reorientation part of the theory.