
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39227653
69914
10.1038/s41598-024-69914-2
Article
Giant impact on early Ganymede and its subsequent reorientation
Hirata Naoyuki hirata@tiger.kobe-u.ac.jp

https://ror.org/03tgsfw79 grid.31432.37 0000 0001 1092 3077 Graduate School of Science, Kobe University, Rokkodai 1-1 657-8501, Kobe, Japan
3 9 2024
3 9 2024
2024
14 199825 4 2024
9 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Ganymede has an ancient impact structure called a furrow system. The furrow system is the largest impact structure in the outer solar system, and the impact should have significantly affected Ganymede’s early history; however, its effects are poorly understood. No attention has been given to the center of the furrow system coinciding with Ganymede's tidal axis, indicating that mass redistribution induced by the furrow-forming impact caused a reorientation (true polar wander) of Ganymede. We propose that the impact ejecta created a mass anomaly that reoriented the impact site toward the tidal axis. We found that an impactor with a radius of 150 km and an incidence angle between 60° and 90° most accurately reproduces the current location of the furrow system. We predict that future explorations would reveal remnant topographic profiles or gravity anomalies associated with the furrow-forming impact and reorientation. Additionally, various possible explanations for the reorientation of Ganymede, such as an impactor-origin mascon beneath the basin or a thickness variation in the lithosphere, should be studied.

Subject terms

Planetary science
Rings and moons
http://dx.doi.org/10.13039/501100001691 Japan Society for the Promotion of Science 20K14538 Hirata Naoyuki http://dx.doi.org/10.13039/100021128 Hyogo Science and Technology Association issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Ganymede is the largest satellite in the solar system and has many unique features, including tectonic troughs known as furrows1–5. Furrows are the oldest surface features recognized on Ganymede because they are crosscut by any impact craters with diameters exceeding 10 km3. Therefore, furrows can provide a window into the early history of Ganymede. Furrows have been proposed to be fragments of multiring impact basin structures, similar to those of the Valhalla or Asgard basins on Callisto1,2,6,7. The largest furrow system is present across Galileo and Marius Regios (the so-called Galileo‒Marius furrow system), and it is the remnant of an ancient giant impact, which extends concentrically from a single point of Ganymede, 21° S 179° W8 (Fig. 1); however, estimating the size of the impactor is complicated because of the absence of an identifiable clear rim8–11. Although there are a few small furrow systems on Ganymede, only the Galileo–Marius furrow system is examined in this study.Figure 1 Distribution of furrows and location of the center of the furrow system shown in the hemisphere that always faces away from Jupiter (top) and the cylindrical projection map of Ganymede (bottom). The distribution of furrows was obtained from Collins et al.10. The gray regions represent geologically young terrain without furrows. Furrows (green lines) exist only on geologically old terrains (black regions).

We propose that the ejecta mass of the furrow-forming impact created a large positive gravity anomaly around the impact center, which subsequently led to the reorientation of Ganymede. No attention has been given to the fact that the center of the furrow system coincides with the longitude of Ganymede's tidal axis; however, this coincidence implies that Ganymede has experienced significant reorientation. Interestingly, the center of the furrow system has a very similar geometric location to that of the center of the Sputnik Planitia (18° N 178° E12), which is the largest impact basin on Pluto. Both centers are located along the tidal longitude and deviate poleward by 20° from the tidal axis. Note that the coordinates of Ganymede and Pluto are defined such that the tidal axes are at a latitude of 0° and longitudes of 0°/180°. Previous studies have proposed that the Sputnik Planitia has a large positive gravity anomaly that caused significant reorientation of Pluto12–14.

In general, a positive gravity anomaly on the surface of a tidally locked satellite leads to the reorientation of the satellite wherein the anomaly approaches the tidal axis; alternatively, a negative gravity anomaly leads to a reorientation wherein the anomaly approaches the rotation axis15–18. The crater hole and ejecta blanket act as negative and positive gravity anomalies, respectively. Both anomalies mostly cancel each other out, but the overall gravity anomaly becomes slightly negative, leading to a net effect of the reorientation of the crater center toward the rotation axis19. In fact, some of the largest impact basins, such as the South Pole Aitken basin (Moon) and Rheasilvia basin (Vesta), are located close to the polar region. In contrast, the Sputnik Planitia basin is located near the tidal axis. Previous studies12,13 have proposed that isostasy can be achieved within the basin floor if Pluto has a global ocean, which, together with the nitrogen deposits within the basin, creates a large positive gravity anomaly. Because the isostatically compensated topography (e.g., supported by the buoyancy of the asthenosphere or subsurface ocean) exhibits no free-air gravity anomaly, the negative gravity anomaly of the basin depression vanishes even if it has a negative topographic profile. In short, the basin depression of the Sputnik Planitia is likely to be compensated for and does not exhibit a negative gravity anomaly, whereas the ejecta blanket is not compensated for and thus exhibits a positive gravity anomaly.

A similar mechanism is appropriate for the case of the furrow-forming impact on Ganymede because the formation of multiring basins on the Jovian icy satellites has been explained by the collapse of a crater basin owing to fluid or warm convective ice beneath a thin ice shell6,20–23. Although it is difficult to argue whether Ganymede had a differentiated interior with an ocean or an undifferentiated interior without an ocean at the time of the formation of the furrow system, numerical simulations24–27 have shown that both a thin ice shell over the ocean and warm convective ice without an ocean can reproduce the shallowing of the basin floor and widening of the crater rim of the multiring basin. Isostasy would have been achieved within the basin floor in both subsurface models (see also Methods). The ejecta blanket alone then played the role of a gravity anomaly, leading to reorientation toward the tidal axis.

We used two different techniques to investigate the reorientation of Ganymede: the first technique balances the degree-2 gravity anomaly caused by a simplified ejecta blanket model against the remnant rotational bulge, and the second technique calculates the perturbation of the moment of inertia from a more complex ejecta blanket model without taking the remnant bulge into account.

Evaluation on the basis of comparison with the remnant bulge

We evaluated the extent to which ejecta mass would have been sufficient to induce significant reorientation in Ganymede. Many significant studies have been conducted on the reorientation of solid bodies in the solar system16–18,28–33. For a rigid spinning tidally locked body, mass redistribution, such as cratering or subsurface diapirism, leads to a new principal axis rotational state (reorientation), while the remnant rotational and tidal bulge stabilizes the rotation axis. Therefore, the true polar-wander solutions are determined by the balance between the remnant bulge and a load (i.e., the mass creating the gravity anomaly)34,35. The dimensionless parameter Q, defined as the ratio of the degree-2 gravitational potential perturbation of the load to the fossil part of the rotational bulge, is useful in determining solutions34,35; this parameter is expressed as follows:1 Q≡35GMg20Ω2R3k2T∗-k2T

where g20 is the degree-2 order-0 coefficient of the gravitational potential perturbation of the mass redistribution at the object surface; G is the gravitational constant; M, R, and Ω are the mass, radius, and rotation angular frequency of Ganymede, respectively; and k2T∗ and k2T are the degree-2 fluid and actual Love numbers, respectively. If the load is too small compared with the existing bulge (i.e., Q≪1), reorientation does not occur because the bulge plays a role in stabilizing the current rotation axis. Here, Q>0 (Q<0) indicates that the load plays the role of a positive (negative) gravity anomaly, and Q=1 indicates that the gravity anomaly of the load is comparable with that of the rotational bulge. As an example, the load of Q=1 leads to a reorientation of 7° if its latitude is 45°16.

If the ejecta blanket with a total volume of Vejecta has a uniform thickness and simple annulus shape with inner and outer radii of Rϕ and Rθ (θ=2ϕ) and isostasy is achieved within the basin radius (Rϕ), g20 is given by.2 g20=ρc(1-Cn)Vejecta25Mcos2ϕ+cosϕcosθ-sin2θ

where ϕ and θ are the angles from the object center, Cn is the degree of isostatic compensation, and ρc is the density of the ejecta blanket (Methods). Note that this simple annulus provides a good approximation of the ejecta distribution because the wavelength of the degree-2 gravitational potential is sufficiently longer than the basin radius. On the basis of the Z model (e.g., ref.36), we assumed that the total volume ejected from a crater (Vejecta) was a fraction of the total volume displaced from the transient crater (Vtc), Vejecta=14Vtc, and that the size of the transient crater was a simple hemisphere, Vtc=23πrtc3. Additionally, we used a relation of Rϕ=1.3rtc (e.g., ref.36). As a result, we obtain the value of Q as a function of the transient crater radius, rtc (Fig. 2). The values of Cn and Δk2=k2T∗-k2T, which depend on the interior structure of Ganymede 4 billion years ago, are highly uncertain. Various estimations exist regarding the thickness of the lithosphere at the time of furrow formation: (i) the thickness of the lithosphere was estimated as 6–10 km on the basis of the width and spacing of furrows37–39, and (ii) the effective elastic thickness was estimated as 0.5 km on the basis of flexural uplift around furrows, although flexural uplift may have formed during a geologically active period of Ganymede long after furrow formation40. If the elastic thickness was t= 10 km or 0.5 km, we can assume that Cn=0.87 or 0.99, respectively (Methods). When t= 10 km, we can assume two cases: Δk2=0.6 if Ganymede had no ocean and its mantle had a high viscosity of η> 1014 Pa s or Δk2=0.03 if Ganymede had an ocean or its mantle had a low viscosity of η< 1012 Pa s (Methods). When t= 0.5 km, we can assume two cases: Δk2=0.6 if Ganymede had no ocean and its mantle had a high viscosity of η> 1014 Pa s, or Δk2=0.005 if Ganymede had an ocean or its mantle had a low viscosity of η< 1012 Pa s (Methods). Figure 2 presents the four cases: t,Δk2=(10 km, 0.6), (10 km, 0.03), (0.5 km, 0.6), and (0.5 km, 0.005). As a result, the volume of the ejecta blanket that creates a gravity anomaly with Q=1 corresponds to 100% of the ejecta created by a transient crater with radii of 480, 160, and 220 km when t,Δk2=(10 km, 0.6), (10 km, 0.03), and (0.5 km, 0.005), respectively (Fig. 2). Therefore, a transient crater larger than 160 ~ 480 km in radius can lead to significant reorientation. When t,Δk2= (0.5 km, 0.6), we can obtain Q≪1 for any transient crater radius.Figure 2 Parameter Q is a function of the transient crater radius, where t represents the lithospheric thickness and where Δk2=k2T∗-k2T represents the difference between the degree-2 fluid and actual Love numbers. Here, Q≫1 (Q≪1) means that the load of the ejecta leads to significant (small) reorientation of Ganymede. Here, the gravity anomaly inside the basin radius is assumed to be zero.

Global distribution of the ejecta blanket and the most stable orientation

The location of the Sputnik Planitia, which deviates 20° from the tidal axis, was explained by the balance between the bulge and mass anomaly of the basin12. However, the center of the furrow system cannot be explained in the same manner if the remnant bulge of Ganymede is readjusted on a short timescale. In general, some fraction of the remnant bulge eventually vanishes with the occurrence of reorientation; a new bulge is immediately formed in response to the new rotation state, and the gravity anomaly ultimately migrates to the tidal axis15. We calculated the center of the gravity anomaly (i.e., the minimum principal axis of the load) of the mass distribution of the ejecta blanket created by the furrow-forming impact. This approach ignores the stabilizing effect of any remnant bulge and would thus be expected to underestimate the size of the impactor compared with the results derived above for a simplified ejecta blanket. Numerical simulations and experimental studies have shown that although the shape of the crater depression is not sensitive to the impact incidence angle unless the angle is very shallow, the shape of the ejecta blanket is strongly sensitive to the impact incidence angle, and most impacts produce asymmetric ejecta patterns (Methods). Unless an ejecta blanket is perfectly symmetric, the center of the load does not match that of the crater. Therefore, we considered the impact incidence angle, utilized the methods described by Hirata et al.41 and Hirata42, calculated the trajectories of ejecta particles, and reproduced the thickness distribution of the ejecta blanket.

In this ejecta model, the initial launch position, velocity, and volume of ejecta particles were determined via the scaling law developed by Housen and Holsapple43 and Raducan et al.44. The ejecta model was obtained by experimental studies on Earth and has been used in studies of secondary craters and ejecta on icy satellites, including Ganymede45,46. Because the furrow-forming impact is in a gravity regime that does not depend on the strength of the target material, the interior structure of Ganymede would hardly affect the ejecta model. For example, an experimental study47 demonstrated that the shattering strength in impact experiments of icy clay samples with water contents of up to 35 wt% is not different from that of basalt, although the static tensile strength of basalt is approximately 8 times greater than that of icy clay samples. Housen and Holsapple43 provided three ejecta models in the gravity regime (water, dry sand, and glass microspheres as target surfaces), although the three ejecta models do not differ markedly from each other. The ejecta model of Housen and Holsapple43 was based on an experiment with a normal impact (incidence angle of 90°). Raducan et al.44 updated the ejecta model by including the incidence angle of the impactor as a parameter. The trajectories of the ejecta particles were solved via Hill’s equation48. We defined ejecta particles that reached below the object surface as colliding with the object and ejecta particles that reached an altitude greater than the Hill radius of the object as escaping from the object. Additionally, we removed the particles that land within the radius of a basin rim to assume an isostatically compensated basin floor.

The left plates in Figs. 3, 4, and 5 show fifteen examples of the global distribution of an ejecta blanket in the cases of an ejecta launch angle of 45°, C4 (dry sand as a target surface) as the set of scaling constants, a basin rim radius of Rϕ=1000 km, and other parameters presented in Supplementary Table 1. We calculated the moment of inertia tensor of the mass distribution of each ejecta blanket and the eigenvalues and eigenvectors of the tensor to obtain the principal moments of inertia and principal axes of each ejecta blanket. We remapped each ejecta blanket in body-fixed coordinates after reorientation to ensure that Ganymede had the most stable rotation; the maximum and minimum principal axes of the mass distribution matched the rotation and tidal axes, respectively (right plates in Figs. 3, 4, and 5). In many cases, the current location of the center of the furrow system can be reproduced by this ejecta model (Figs. 3 and 4). However, this could not be reproduced in some cases (Fig. 5). Although reproducing the current location of the furrow system involves many parameters and is complicated, our results suggest that (i) the impact incidence angle is between 60° and 90° or (ii) the uprange direction is roughly on the west side. This is almost independent of the initial location of the impact site (Fig. 4).Figure 3 (left) Global distribution of ejecta blankets before reorientation using impact parameters #1–5 presented in Supplementary Table 1. The white, gray, and black “ + ” symbols represent the minimum, intermediate, and maximum principal axes, respectively. (right) Global ejecta thickness after reorientation to ensure that Ganymede had the most stable rotation. The vertical and horizontal axes in each plate indicate the latitude and east longitude, respectively, where the top left corner is 90° N and 0° E, and the bottom right corner is 90° S and 360° E. The color bar is represented on a logarithmic scale.

Figure 4 Global distribution of ejecta blankets before reorientation (left) and after reorientation (right) using impact parameters #6–10 presented in Supplementary Table 1. The white, gray, and black “ + ” symbols represent the minimum, intermediate, and maximum principal axes, respectively. The vertical and horizontal axes in each plate indicate the latitude and east longitude, respectively, where the top left corner is 90° N and 0° E, and the bottom right corner is 90° S and 360° E. The color bar is represented on a logarithmic scale.

Figure 5 Global distribution of ejecta blankets before reorientation (left) and after reorientation (right) using impact parameters #11–15 presented in Supplementary Table 1.

In this ejecta model, a small mound with a height of approximately 10 km and a width of 300 ~ 600 km owing to the accumulation of ejecta appears near the sub-Jovian point (Figs. 3d, e, 4a, b). This may agree with the small mound with a height of 3 km located near the sub-Jovian point (1.5° E, 0.5° N) discovered by the Juno flyby and legacy data reanalysis49,50. When the uprange direction is roughly on the west side, the mound tends to appear near the sub-Jovian point after reorientation because the small mound shifted slightly from the exact opposite point of the impact owing to the effect of Ganymede's rotation. Viscous relaxation, isostatic readjustment, and/or uncertainty of the ejecta model may be responsible for the difference in height and shape between the discovered mound and our ejecta mound.

Furthermore, via the Monte Carlo method, we examined the probability that the center of the furrow system after the most stable reorientation matched the current location. The initial longitude, latitude, uprange direction, and incidence angle of the impactor follow the functions of 2πx1,acos2x2-1,2πx3,and,acos(x4), respectively, where 0≤x1,x2,x3,x4<1 are uniformly distributed random numbers51. We examined cases for an impact incidence angle between 60° and 90° or between 30° and 60°; a target material of C4 or C1 (water); a transient crater radius of 500, 600, 700, 800, or 900 km; and an ejecta launch angle of 25°, 30°, or 45° (angle between the initial launch velocity vector and the object surface). Note that an ejecta launch angle of ~ 28° may be the most appropriate for a giant impact because of the curvature of the target surface52,53. For each case, 1000 trials were performed. Figure 6 shows the initial and final locations of the crater centers in the case of an impact incidence angle between 60° and 90°, a target material of C4, and an ejecta launch angle of 30° or 45°. This figure shows that the basin center tends to shift poleward from the tidal axis by the basin radius along the tidal longitude. Figure 7 shows the relationship between the transient crater and the probability that the center of the furrow system after the most stable reorientation moved within 5° or 10° from the point deviated poleward by 20° from the tidal axis.Figure 6 (top) Initial location of the center of the crater, (middle) final location of the center of the crater after reorientation, and (bottom) final location in the coordinate centered at the tidal axis. We have displayed cases with ejecta launch angles of 45° (left) and 30° (right). Note that because the true polar wander solution does not distinguish between the southern and northern hemispheres, points plotted in the Northern Hemisphere are shown with the opposite sign on the bottom plate. Owing to satellite rotation, it is not symmetrical with respect to the east‒west direction.

Figure 7 Probability that the center of the furrow system after the most stable reorientation moves within 5° (left) or 10° (right) from the point that deviates poleward by 20° from the tidal axis as a function of the transient crater radius. Here, we assume an impact incidence angle (inc) between 30° and 60° or between 60° and 90°; an ejecta launch angle (ej) of 25°, 30°, or 45°; and a target material (tar) of C1 or C4. The legend is common on both plates.

For example, the probability that the impact center falls within 5° from the point that deviates poleward by 20° from the tidal axis under the conditions of an impact incidence angle between 60° and 90°, a transient crater radius of 700 km, and an ejecta launch angle of 30° is 57.2% (Fig. 7). This transient crater radius creates a gravity anomaly with Q>1, as shown in Fig. 2. If reorientation does not occur, the probability is 0.76%. Note that there are four locations on Ganymede that deviate poleward by 20° from the tidal axis. The probability tends to be maximum when the size of the transient crater radius is 700 or 800 km, which is consistent with the above arguments for the Q value. A transient crater radius of 700 or 800 km is produced by an impactor with a radius of 150 km11. In conclusion, an impactor with a radius of 150 km and an incidence angle between 60° and 90° most likely reproduces the current location of the furrow system.

Discussion

Although we demonstrated that the load of the ejecta blanket can reorient the impact site toward the tidal axis, a variety of other possible explanations for the reorientation of Ganymede cannot be ruled out because of the high uncertainty of the conditions of Ganymede 4 billion years ago. First, because inhomogeneous lithospheric thickness leads to an inhomogeneous degree of isostatic compensation, even if the ejecta blanket is perfectly symmetric, the ejecta rim supported by a thicker lithosphere would have greater gravity anomalies. Therefore, variations in lithospheric thickness, rather than an asymmetric ejecta distribution, may be responsible for the deviation of the center of the furrow system from the tidal axis by 20°. Second, it is proposed that the variation in the thickness of the icy shell driven by the latitudinal variation in solar radiation leads to the reorientation of Ganymede32 or that a captured rocky impactor core beneath the Sputnik Planitia basin leads to the reorientation of Pluto33, although we do not discuss these hypotheses in this paper. Third, the mass anomaly of the ejecta blanket and bulge supported by the lithosphere may vanish on a long timescale, whereas the presence of localized gravity anomalies (e.g., ref.54) indicates that relaxation does not always occur and that long-term mass anomalies can persist. However, it is difficult to estimate the relaxation time of the ejecta blanket and the remnant bulge (Method). If the gravity anomaly of the ejecta blanket vanishes before that of the remnant bulge, Ganymede should reorient back to its initial position. Alternatively, the cooling of Ganymede during the timeframe possibly increased the thickness of the lithosphere and should have fixed the orientation of Ganymede. These issues also need to be discussed in further studies.

The reorientation of Ganymede provides a window into the early history of the Jovian satellites. It is proposed that the 10-km-thick (0.5-km-thick) lithosphere of Ganymede indicates a surface heat flow of approximately 40 mW/m2 (60–80 mW/m2) at the time of furrow formation39,40, which exceeds the maximum plausible heat fluxes expected from radiogenic heating alone 4.5 billion years ago (nearly 27 mW/m255). Therefore, additional energy sources, such as tidal heating, gravity segregation, and/or loss of accretional heat40, may have been present at the time of furrow formation. In particular, large-amplitude librations and/or nonsynchronous rotations can themselves cause significant heating due to tidal dissipation. We estimate that Ganymede should have attained its new orientation approximately 1000 years after the furrow-forming impact as a result of tidal dissipation (Methods). In addition, a 150 km-radius impactor was proposed to pump early Ganymede's free eccentricity56.

Heliocentric impactors preferentially hit the leading hemisphere of a tidally locked satellite, and the impact crater density of Ganymede decreases from the apex (center of the leading hemisphere) to the antapex (center of the trailing hemisphere) of its motion; however, the degree of apex–antapex asymmetry in various tidally locked satellites is considerably lower than that of theoretical estimates for ecliptic comets51,57,58. Previous studies51,57,58 have proposed various possible explanations for this difference, such as crater saturation, nonsynchronous rotation, reorientation, and nearly isotropic comets as the dominant impactor population. If the load of the ejecta has fixed the orientation of Ganymede, further reorientation would be unlikely in the case of Ganymede. This could be an area that needs further study.

The amount of ejecta in the Valhalla basin on Callisto is equivalent to Q= 10, which is sufficiently large to cause a reorientation of Callisto (Methods); however, Valhalla is far from both the pole and tidal axes. Possible explanations include the following: (i) Callisto's rotation at the time of formation of the Valhalla basin was faster than that at present, (ii) isostasy is not archived on the basin floor, (iii) other surface loads exist on Callisto, or (iv) Callisto has a large remnant bulge; however, we cannot argue these explanations conclusively.

Many areas of Ganymede still have not been imaged with sufficient resolution5, and further data from future explorations are required for further discussions of the tectonic landforms formed by reorientation and the age of furrow formation. For example, tectonic patterns (other than furrows) resulting from Ganymede’s reorientation have not been discovered, although tectonic patterns resulting from the reorientation of Pluto and the Jovian satellites have been investigated13,28,31. New fractures may not have formed if the reactivation of existing fractures, or furrows, absorbed reorientation-induced stresses. The current accuracy of the gravity and topography measurements of Ganymede59–61 are insufficient for comparisons between individual topographic features owing to the lack of a global shape model of Ganymede and a local topographic map around the center of the furrow system. Although it is not clear that the gravity anomaly has vanished, the topographic profiles of the ejecta blankets should still be visible in the topography even if they are mostly viscously relaxed. Future explorations would reveal such a remnant of topographic profiles associated with the furrow-forming impact and the reorientation of Ganymede, which would provide insights into Ganymede’s early history and highlight its differences compared with those of other Jovian satellites.

Methods

The degree of isostatic compensation and the value of g20

Following Melosh19, we assumed the following conditions for a crater basin: (i) the ejecta blanket of excavated material has a uniform thickness, te, and a simple annulus shape with inner and outer radii of Rϕ and Rθ, respectively, where ϕ and θ are the angles from the center of the object, respectively, and the outer radius is twice the inner radius, θ=2ϕ; and (ii) the basin floor has a uniform depth, h, and a simple circular shape with a radius of Rϕ. In this case, a simple geometric solution gives the volume of the ejecta blanket and hole as follows:3 Vejecta=2πR2te(cosϕ-cosθ),and

4 Vhole=2πR2h(1-cosϕ)

In general, both volumes are equal because of the conservation of mass19, even if the very shallow basin was created after the collapse of a transient crater; subsequently, we can use the simple relation of h=te(cosϕ-cosθ)/(1-cosϕ). We considered the variations in the density of a nearly spherical object expressed by variations over a series of concentric spherical shells, and the density perturbations were azimuthally symmetric. The density perturbation in this model of the ejecta blanket model is expressed as follows:5 ρr-ρ0r=0(ϕ>φ>0)ρc(θ>φ>ϕ)0(2π>φ>θ)R<r<R+te-ρc(ϕ>φ>0)0(2π>φ>ϕ)R-h<r<R0r<R-h

where ρ0r is the mean density for radius r, ρc is the density of the crust, and φ is the colatitude from the crater center. The degree-2 order-0 coefficient of the density perturbation is expressed as follows62:6 ρ20r=54∫0πρr,φ3cos2φsinφ-sinφdφ

We can obtain.7 ρ20r=54ρccosθsin2θ-cosϕsin2ϕR<r<R+te-54ρccosϕsin2ϕR-h<r<R0r<R-h

The degree-l order-m spherical harmonic coefficients of the gravitational potential at an arbitrary distance, s(s≥R), are given by62:8 glm(s)=4πM2l+1sl∫0Rrl+2ρlmrdr

where s is the distance from the center. Therefore, if this density perturbation is supported by an infinitely rigid lithosphere, the degree-2 order-0 gravitational potential coefficient at the object surface can be described as follows:9 g20=πR2ρc5Mtecosθsin2θ-h+tecosϕsin2ϕ

If we assume that h=te(cosϕ-cosθ)/(1-cosϕ)), this gravitational coefficient is always negative, g20<0, which gives Q<0.

The actual shapes of the eject blankets in Fig. 3 are not simple annulus shapes. However, this simple annulus provides a good approximation of the ejecta distribution because the wavelength of the degree-2 gravitational potential is sufficiently longer than that of the basin radius and mostly averages inhomogeneity or asymmetry in the ejecta distribution. As a simple example, if we define θ=1.3ϕ or θ=2.5ϕ instead of θ=2ϕ, the value of Q changes by only approximately 20–30%. Inhomogeneity in the tangential direction should have little effect on g20 if the total mass of ejecta as a function from the center is equal to the simple annulus shape because Eq. (6) does not depend on the tangential direction.

A mass excess/deficit that is compensated for by isostasy does not create a free-air gravity anomaly, whereas the ejecta mass supported by the bending and membrane stresses of the lithosphere creates a gravity anomaly. Therefore, the gravity anomaly decreases by a factor of (1-Cn), where the degree of isostatic compensation, Cn, is described as follows63:10 Cn=1-3ρm2n+1ρ¯σN3-4N2+τN-2+N-(1-ν)N-(1-ν)-3ρm2n+1ρ¯-1

where N=nn+1, n is the degree of spherical harmonics of the horizontal width of the load (n=2πR/w, where w is the horizontal width of the load),ρ¯ is the mean density of the satellite, and τ and σ are nondimensional parameters that are defined as follows:11 τ≡EtR2gρm-ρcandσ≡τ121-ν2tR2

where t is the elastic thickness;ρm is the density of the mantle underlying the lithospheric plate; E and ν are the Young’s modulus and Poisson’s ratio of the lithosphere, respectively; and g is the surface gravity. Here, Cn=0 indicates that the load is fully supported by the lithosphere, and Cn=1 indicates that the load is fully compensated by isostasy. We assumed that E = 9 GPa, ν =0.32 (ref. 64), ρm = 1000 kg m−3, and ρc = 930 kg m−3 (ref. 65). We assumed that w= 1000 km because the ejecta blanket spreads across several basin radii. Note that the value of w hardly affects Cn: Cn=0.871 (w= 3000 km),Cn=0.875 (w= 1000 km), and Cn=0.868 (w= 500 km) for t= 10 km, and Cn=0.992 (w= 3000 km),Cn=0.993 (w= 1000 km), and Cn=0.993 (w= 500 km) for t= 0.5 km. Consequently, we assume Cn=0.87 for t= 10 km and Cn=0.99 for t= 0.5 km.

If the value of Cn is homogeneous inside and outside the basin, the gravity anomaly is always negative (g20<0); thus, reorientation that orients the basin center to the tidal axis does not occur. At least Cn outside the basin must be smaller than Cn inside the basin to obtain g20>0. Nimmo et al.12 assumed that isostasy is archived inside the Sputnik Planitia basin. Similarly, in the case of Ganymede, the gravity anomaly becomes positive (g20>0) if (i) isostasy is archived inside the basin, which has the same attribute as h=0, and (ii) it is not archived outside the basin (i.e., the load of the ejecta is partially supported by the lithosphere). This assumption is natural for the following reasons. First, numerical simulations26 or impact experiments66 have proposed that the impactor can penetrate Europa’s ice shell and create conduits to the underlying ocean. Similarly, if Ganymede has an ocean and its ice shell is sufficiently thin, isostasy is easily achieved within the basin floor by the inflow of water. The furrow-forming impact alone may not create a large ocean because the size of a local melt pool created by the impact is up to 6 times the impactor radii67, which is less than 10% of the total volume of Ganymede. However, isostasy would be achieved within the basin floor even if Ganymede’s interior is warm convective ice without an ocean. Numerical simulations of warm convective ice without the ocean11,27 have shown that cold lithospheric material is replaced by warm subsurface material and that the postimpact lithospheric thickness becomes quite thin within the transient crater radius when the transient crater size is sufficiently larger than the lithospheric thickness. The lithospheric thickness recovered by cooling from the surface is roughly determined by the characteristic thermal diffusion depth, κtd, where κ is the thermal diffusivity and where td is the timeframe68. The relaxation time of a load with a horizontal width of 100 km is 25 years (see below). If we use κ= 2.9 × 10–5 for ice69,70 and td= 25 years, we obtain κtd= 150 m. An elastic thickness of t= 150 m does not support a load above the lithosphere at all (Cn=1). In other words, the mass anomaly within the basin floor quickly disappears before the lithospheric thickness recovers. Therefore, isostasy can be achieved within the basin floor of the furrow system, regardless of the two surface models used. Second, numerical simulations11,27 have shown that the postimpact lithospheric thickness and temperature do not become thin and warm outside the basin. Numerical simulations for the Sputnik Planitia71 also show that hot preimpact thermal structures or cold thermal structures with an ocean greater than 150 km thick produce basins that are broadly isostatically compensated inside the basin and that are not fully isostatically compensated outside the basin. Therefore, it is likely that isostasy is not fully archived outside the basin. Third, furrow formation is unlikely to cause a significant change in the rheology of the lithosphere overall before and after the furrow-forming impact because individual furrows occur in the brittle part of the lithosphere6,7, and the brittle part above the brittle‒ductile transition layer should already include many cracks and faults due to impacts or tectonics. Notably, if the lithosphere is fully damaged by impact and isostasy is fully and broadly archived even outside the basin, any gravity anomalies will not be created by ejecta. If so, we should conclude that a coincidence or another hypothesis, such as the impactor-origin silicate mascon beneath the basin or the variation in the lithospheric thickness, is responsible for the current position of the center of the furrow system coinciding with the tidal axis.

Based on the above assumption (h=0), the dimensionless parameter becomes.

Q=3πGρc(1-Cn)teΩ2RΔk2cosθsin2θ-cosϕsin2ϕ

or12 Q=3Gρc(1-Cn)Ω2RΔk2Vejecta2R2cos2ϕ+cosϕcosθ-sin2θ

The radius of the transient crater of the Gilgamesh basin is 135 km (ref. 46), which indicates a load of Q=0.027 (Δk2=0.6) or 0.54 (Δk2=0.03). The Gilgamesh basin does not cause significant reorientation of Ganymede, although a reorientation caused by the Gilgamesh basin has previously been proposed28. As another example, the mound at the sub-Jovian point, discovered by the Juno flyby and legacy data reanalysis, with a height of 3 km and an oval of 450 km × 750 km (refs. 49,50), has a value of Q=0.018 (Δk2=0.6) or 0.36 (Δk2=0.03), assuming h=-3 km, te= 0, Rϕ=300 km, and Cn=0.87 (t = 10 km). Following Eq. (41) in Matsuyama and Nimmo16, the mound leads to very small reorientation of less than 2.6° if its latitude is 45°. We obtain a value of Q=0.17 for the mound if we assume that h=-3 km, te= 0, Rϕ=300 km, Cn=0. 99 (t = 0.5 km), and Δk2=0.005.

The degree 2 Love number of Ganymede

The value of Δk2=k2T∗-k2T is determined by the internal structure of Ganymede and has a significant effect on the value of Q. Many studies have investigated the tidal response of Ganymede, assuming various interior structure models72–74. For example, Moore and Schubert72 reported that k2T is approximately 0.01 if the present-day Ganymede has no ocean and its mantle has a high viscosity of η> 1014 Pa s or approximately 0.4 to 0.6 if Ganymede has an ocean or its mantle has a low viscosity of η< 1012 Pa s. Although it is unclear whether Ganymede’s interior is differentiated or undifferentiated, we estimate the Love number on the basis of a four-layer Ganymede model. Here, we assume that Ganymede consists of (1) a 10 km-thick elastic icy shell, (2) a viscoelastic icy mantle, (3) a viscoelastic rocky mantle, and (4) a liquid metallic core, where the rigidities of the icy shell and rocky mantle are 10 and 100 GPa, respectively; the densities of the icy shell, icy mantle, rocky mantle, and metallic core are 1050, 1050, 1050, 3100, and 5150 kg/m3, respectively; the outer radii of the icy shell, icy mantle, rocky mantle, and metallic core are 2638, 2628, 1745, and 710 km, respectively; and the viscosity of the rocky mantle is 1020 Pa s. The value of k2T∗ is defined as the value when the elastic icy shell is assumed to be a liquid ocean layer. Those rigidities, densities, outer radii, and viscosities are from Kamata et al.73. We use a code released by Kamata75 to calculate the periodic spheroidal deformation of a planetary body to obtain the Love number. As a result, when the viscosity of the icy mantle is 1011, 1013, or 1015 Pa s, we obtain (k2T,k2T∗) = (0.579,0.604), (0.335,0.605), or (0.018,0.605), respectively. If the subsurface ocean layer (with a density of 1050 kg/m3) is present instead of the viscoelastic icy mantle, we can obtain (k2T,k2T∗) =(0.579,0.604). Considering the large uncertainties, we studied the two cases of Δk2=0.03 and Δk2=0.6 for t = 10 km. If we assume a 0.5 km-thick elastic icy shell instead of a 10 km-thick elastic icy shell, when the viscosity of the icy mantle is 1011, 1013, or 1015 Pa s, we obtain (k2T,k2T∗) =(0.599,0.604), (0.446, 0.605), or (0.018,0.605); therefore, we study the two cases of Δk2=0.005 and Δk2=0.6 for t = 0.5 km. This combination of Δk2=0.6 and t = 0.5 km may be unlikely because the very thin lithosphere and cold interior are not consistent.

Tidal decay of rotation

Following a large impact event, tidal friction quickly dampens the motions of the satellite’s orientation, such as librational oscillations and nonsynchronous rotation, and the satellite attains its new orientation on a short timescale19. The timescales of the tidal decay of librational rotation and nonsynchronous rotation have been estimated in previous studies19,76,77, and all three studies were equivalent. Following Noyelles et al.77, the timescale is expressed as follows:13 Tdamp=23Qdk2TGCno3R5,

where C is the maximum moment of inertia about the center of mass; no is the mean orbital motion; and Qd is the specific dissipation factor of the satellite. Assuming a homogenous interior (C=0.4MR2) and reasonable values for warm icy satellites,k2T=0.5 and Qd=100 (ref. 78), we can obtain Tdamp=870 years. If Ganymede was cold and rigid (Qd/k2T ~ 105), we can obtain Tdamp=4.4 million years.

Viscous relaxation of the load and bulge

The mass anomaly not supported by the lithosphere disappears owing to viscous relaxation on a timescale of η/ρmgw (ref. 68), where w is the horizontal width of the load. The icy materials underlying the elastic layer are convecting, and their reference viscosity is as low as η= 1012 to 1017 Pa s (ref. 79). For example, a 100 km-wide load not supported by the lithosphere would vanish in less than 25 years. Similarly, the density perturbations below the lithosphere created by heating and redistribution of the impact crater would vanish on a short timescale. If the relaxation time of the remnant tidal/rotational bulge is given by η/ρmgR (ref. 17), the bulge would vanish in 1 year. However, this estimate is suitable only for an isoviscous case, whereas the internal structure of Ganymede is much more complicated and, in particular, may have a rigid outer shell that can maintain a remnant bulge over geological time. In short, the actual bulge relaxation timescale is extremely uncertain. It is estimated that old craters on the dark terrain of Ganymede have relaxed viscously for 1 billion years55,80. The timescale of relaxation of the mass anomaly of the ejecta blankets and the tidal/rotational bulge may also be approximately 1 billion years.

The Valhalla basin

The original crater radius of the Valhalla basin was estimated to be 500 km on the basis of the mapping of ejecta and secondaries81. The thickness of the lithosphere was estimated to range from 15 to 20 km during the formation of the Valhalla basin6. Assuming Rϕ=500 km, t= 20 km,h= 0, w= 500 km, n= 15, Cn=0.66, and physical parameters for Callisto, we can obtain Q=199.8 (Δk2=0.03) or Q=10.0 (Δk2=0.6) for the Valhalla basin, which should lead to a significant reorientation of Callisto. Possible explanations for this include Callisto's rotational period, incomplete isostatic compensation within the basin floor, and other surface loads. If Callisto's rotational period was 50 h (currently 400 h), we could obtain Q=0.16 (Δk2=0.6). The Valhalla basin is one of the oldest surface features on Callisto21,82, and Callisto at the time of the formation of the Valhalla basin may not have yet reached a synchronous rotation state. However, this scenario may be unlikely because Callisto attains a synchronous rotation state in 2 million years if we use Qd/k2 ~ 20 (ref. 83). If Callisto was cold and rigid (Qd/k2 ~ 105), Callisto obtained a synchronous rotation state 1 billion years after its formation. Note that because Ganymede attains a synchronous rotation state at 0.66 million years for Qd=100 after its formation76, Ganymede’s rotational period at the time of furrow formation should not differ much from the present period. If isostasy is not achieved within the basin floor of the Valhalla basin (i.e., h>0), incomplete isostatic compensation of the negative topographic profile within the basin would decrease Q. Finally, as many areas of Callisto have not yet been imaged with sufficient resolution, we cannot rule out the possibility that large undiscovered loads exist on Callisto’s surface. Additionally, Callisto may have had a large remnant bulge, such as the equatorial bulge of Iapetus, because of Callisto’s undifferentiated and inert interior.

Note for the calculation of the ejecta blanket

Many studies have focused on the shape of the ejecta blanket, which shows that the mass distribution of ejecta is strongly sensitive to the impact incidence angle and highly asymmetric84–88. Therefore, most naturally, the impact incidence angle is responsible for the asymmetric load to explain the location of the center of the furrow system. However, there could be various alternative explanations, as we described above.

It is not unnatural that a thick ejecta blanket would be superimposed on where the furrows lie because this is also the case in the Valhalla basin. The Valhalla basin can be divided into three distinct zones outward from the center: a central smooth zone (rb < 360 km, where rb is the distance from the basin center), an inner ridge and trough zone (360 km < rb < 950 km), and an outer graben zone (950 km < rb < 1900 km). The crater density in the central smooth zone and inner ridge and trough zone is 3.5 times lower than that in the adjacent unmodified cratered terrain, whereas the crater density in the outer graben zone is intermediate between that in the inner zones and unmodified cratered terrain and apparently increases linearly outward; moreover, the deficiency of craters is interpreted to be due to obliteration by continuous ejecta blankets3,21. If troughs and grabens were created first, they would be damaged by ejecta particles. Therefore, troughs and grabens were likely created after the accumulation of ejecta. Similarly, no structures older than the furrows have been found on the Ganymede surface, and no craters cut by the furrows have been found3. Therefore, all craters on Ganymede were likely created after furrow formation (perhaps, all ancient craters created before furrow formation were removed by ejecta). Compared with ejecta particles flying in a ballistic trajectory, the timing of furrow formation formed by the asthenospheric flow would be delayed.

The thickness of the ejecta blankets at the rim is 100 km in our calculation; however, this is unrealistic. This is mainly because our ejecta model does not consider any subsequent movement after the landing of ejecta particles, which creates a point of discontinuity in the thickness along the basin rim. In reality, some dispersion induced by mass wasting or re-ejection should occur. The ejecta volume at a point of discontinuity is not large; thus, we expect that it hardly affects our results. Nevertheless, we recalculated the ejecta distribution in the vicinity of the basin rim, assuming a realistic surface slope. Specifically, we removed ejecta steeper than a given slope (λ) around the basin rim, calculated the moment of inertia tensor of the mass distribution of the ejecta blanket whose steep part was removed, and remapped the ejecta blanket after the most stable reorientation (Supplementary Fig. 1). Here, we assume that the removed ejecta do not create a mass anomaly. Although this calculation does not assume a particular geologic process, it is known that a portion of the ejecta rim moves into the basin floor due to lateral movement, such as mass wasting. Thus, the ejecta mass within the basin does not create a free-air gravity anomaly because we assume that isostasy is archived inside the basin. The 15 examples in this figure correspond to the 15 examples shown on the right side of Figs. 3, 4, and 5. As a result, among the 15 examples, only #3 for λ≤ 20 degrees is different from the original example. For the rest, we find that almost the same principal axis rotation state is obtained. It is unclear how much gradient is actually allowed on Ganymede’s surface; however, this figure shows that the discontinuities around the basin rim hardly affect our results. This is mainly because the unrealistically thick ejecta rim is created by the accumulation of ejecta in a very small area, and the removed ejecta mass was approximately less than 10% of the total mass when λ≥ 20 degrees and less than 20% when λ= 5 degrees.

We performed similar numerical simulations for Pluto. The simulations are almost the same as those for Ganymede, other than the trajectories of the ejecta particles. The trajectories were solved via the equation of motion for the circular restricted three-body motion (Pluto, Charon, and an ejecta particle)48. Supplementary Fig. 2 shows a similar probability as a function of the transient crater radius, similar to that in Fig. 7. The mean radius of the Sputnik Planitia is 550 km (ref. 13), which corresponds to a transient crater radius of 420 km. For this transient crater radius, the probability that the impact center is within 5° from the point that deviates poleward by 20° from the tidal axis is approximately 20–30%. Therefore, the center of the Sputnik Planitia of Pluto, which, together with the furrow system, deviates poleward by 20° from the tidal axis, can be explained by a similar model.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-69914-2.

Acknowledgements

We appreciate Shunichi Kamata and Ren Ikeya for helpful comments. This study was partly supported by the JSPS Grants-in-Aid for Scientific Research (Nos. 20K14538 and 20H04614) and Hyogo Science and Technology Association.

Author contributions

N.H. carried out everything.

Code availability

The software for calculating ejecta blankets is available at GitHub (https://github.com/naoyukihirata/ganymede-ejecta).

Competing interests

The author declares no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Smith BA The Jupiter system through the eyes of voyager 1 Science 1979 204 4396 951 972 10.1126/science.204.4396.951 17800430
Smith, B. A. et al. The Jupiter system through the eyes of voyager 1. Science 204(4396), 951–972 (1979).17800430 10.1126/science.204.4396.951
2. Smith BA The Galilean satellites and Jupiter: voyager 2 imaging science results Science 1979 206 4421 927 950 10.1126/science.206.4421.927 17733910
Smith, B. A. et al. The Galilean satellites and Jupiter: voyager 2 imaging science results. Science 206(4421), 927–950 (1979).17733910 10.1126/science.206.4421.927
3. Passey QR Shoemaker EM Morrison D Craters and basins on Ganymede and Callisto: Morphological indicators of crustal evolution Satellites of Jupiter 1982 University of Arizona Press 379 434
Passey, Q. R. & Shoemaker, E. M. Craters and basins on Ganymede and Callisto: Morphological indicators of crustal evolution. In Satellites of Jupiter (ed. Morrison, D.) 379–434 (University of Arizona Press, 1982).
4. Shoemaker EM Wolfe R Morrison D Cratering time scales for the Galilean satellites Satellites of Jupiter 1982 Tucson University of Arizona Press 277 339
Shoemaker, E. M. & Wolfe, R. Cratering time scales for the Galilean satellites. In Satellites of Jupiter (ed. Morrison, D.) 277–339 (University of Arizona Press, Tucson, 1982).
5. Pappalardo RT Collins GC Head J Helfenstein P McCord TB Moore JM Prockter LM Schenk PM Spencer JR Bagenal F Geology of Ganymede Jupiter: The Planet, Satellites and Magnetosphere 2004 Cambridge University Press 363 396
Pappalardo, R. T. et al. Geology of Ganymede. In Jupiter: The Planet, Satellites and Magnetosphere (eds Bagenal, F. et al.) 363–396 (Cambridge University Press, 2004).
6. McKinnon WB Melosh HJ Evolution of planetary lithospheres: Evidence from multiringed structures on Ganymede and Callisto Icarus 1980 44 454 471 10.1016/0019-1035(80)90037-8
McKinnon, W. B. & Melosh, H. J. Evolution of planetary lithospheres: Evidence from multiringed structures on Ganymede and Callisto. Icarus 44, 454–471 (1980).10.1016/0019-1035(80)90037-8
7. Melosh HJ A simple mechanical model of Valhalla Basin Callisto. J. Geophys. Res. Solid Earth 1982 87 B3 1880 1890 10.1029/JB087iB03p01880
Melosh, H. J. A simple mechanical model of Valhalla Basin. Callisto. J. Geophys. Res. Solid Earth 87(B3), 1880–1890 (1982).10.1029/JB087iB03p01880
8. Schenk PM McKinnon WB Ring geometry on Ganymede and Callisto Icarus 1987 72 209 234 10.1016/0019-1035(87)90126-6
Schenk, P. M. & McKinnon, W. B. Ring geometry on Ganymede and Callisto. Icarus 72, 209–234 (1987).10.1016/0019-1035(87)90126-6
9. Zuber MT Parmentier EM A geometric analysis of surface deformation: Implications for the tectonic evolution of Ganymede Icarus 1984 60 200 210 10.1016/0019-1035(84)90148-9
Zuber, M. T. & Parmentier, E. M. A geometric analysis of surface deformation: Implications for the tectonic evolution of Ganymede. Icarus 60, 200–210 (1984).10.1016/0019-1035(84)90148-9
10. Collins, G.C., Patterson, G.W., Head, J.W., Pappalardo, R.T., Prockter, L.M., Lucchitta, B.K. & Kay, J.P. Global geologic map of Ganymede: U.S. Geological Survey Scientific Investigations Map 3237, pamphlet 4 p., 1 sheet, scale 1:15,000,000, 10.3133/sim3237 (2013).
11. Hirata N Suetsugu R Ohtsuki K A global system of furrows on Ganymede indicative of their creation in a single impact event Icarus 2020 352 113941 10.1016/j.icarus.2020.113941
Hirata, N., Suetsugu, R. & Ohtsuki, K. A global system of furrows on Ganymede indicative of their creation in a single impact event. Icarus 352, 113941 (2020).10.1016/j.icarus.2020.113941
12. Nimmo F Reorientation of Sputnik Planitia implies a subsurface ocean on Pluto Nature 2016 540 94 96 10.1038/nature20148 27851735
Nimmo, F. et al. Reorientation of Sputnik Planitia implies a subsurface ocean on Pluto. Nature 540, 94–96 (2016).27851735 10.1038/nature20148
13. Keane JT Matsuyama I Kamata S Steckloff JK Reorientation and faulting of Pluto due to volatile loading within Sputnik Planitia Nature 2016 540 90 93 10.1038/nature20120 27851731
Keane, J. T., Matsuyama, I., Kamata, S. & Steckloff, J. K. Reorientation and faulting of Pluto due to volatile loading within Sputnik Planitia. Nature 540, 90–93 (2016).27851731 10.1038/nature20120
14. Johnson PE Keane JT Young LA Matsuyama I New constraints on Pluto's Sputnik Planitia ice sheet from a coupled reorientation-climate model Planet. Sci. J. 2021 2 5 194 10.3847/PSJ/ac1d42
Johnson, P. E., Keane, J. T., Young, L. A. & Matsuyama, I. New constraints on Pluto’s Sputnik Planitia ice sheet from a coupled reorientation-climate model. Planet. Sci. J. 2(5), 194 (2021).10.3847/PSJ/ac1d42
15. Gold T Instability of the Earth's axis of rotation Nature 1955 175 526 529 10.1038/175526a0
Gold, T. Instability of the Earth’s axis of rotation. Nature 175, 526–529 (1955).10.1038/175526a0
16. Matsuyama I Nimmo F Rotational stability of tidally deformed planetary bodies J. Geophys. Res. Planets 2007 112 E11003 10.1029/2007JE002942
Matsuyama, I. & Nimmo, F. Rotational stability of tidally deformed planetary bodies. J. Geophys. Res. Planets 112, E11003 (2007).10.1029/2007JE002942
17. Nimmo F Matsuyama I Reorientation of icy satellites by impact basins Geophys. Res. Lett. 2007 34 L19203 10.1029/2007GL030798
Nimmo, F. & Matsuyama, I. Reorientation of icy satellites by impact basins. Geophys. Res. Lett. 34, L19203 (2007).10.1029/2007GL030798
18. Matsuyama I Nimmo F Mitrovica JX Planetary reorientation Ann. Rev. Earth Planet. Sci. 2014 42 605 364 10.1146/annurev-earth-060313-054724
Matsuyama, I., Nimmo, F. & Mitrovica, J. X. Planetary reorientation. Ann. Rev. Earth Planet. Sci. 42, 605–364 (2014).10.1146/annurev-earth-060313-054724
19. Melosh HJ Large impact craters and the moon's orientation Earth Planet. Sci. Lett. 1975 26 353 360 10.1016/0012-821X(75)90011-4
Melosh, H. J. Large impact craters and the moon’s orientation. Earth Planet. Sci. Lett. 26, 353–360 (1975).10.1016/0012-821X(75)90011-4
20. Melosh HJ McKinnon WB The mechanics of ringed basin formation Geophys. Res. Lett. 1978 5 985 988 10.1029/GL005i011p00985
Melosh, H. J. & McKinnon, W. B. The mechanics of ringed basin formation. Geophys. Res. Lett. 5, 985–988 (1978).10.1029/GL005i011p00985
21. Schenk PM The geology of Callisto J. Geophys. Res. 1995 100 19023 19040 10.1029/95JE01855
Schenk, P. M. The geology of Callisto. J. Geophys. Res. 100, 19023–19040 (1995).10.1029/95JE01855
22. Kadel SD Chuang FC Greeley R Moore JM Geological history of the tyre region of Europa: A regional perspective on Europan surface features and ice thickness J. Geophys. Res. Planets 2000 105 22657 22669 10.1029/1999JE001203
Kadel, S. D., Chuang, F. C., Greeley, R. & Moore, J. M. Geological history of the tyre region of Europa: A regional perspective on Europan surface features and ice thickness. J. Geophys. Res. Planets 105, 22657–22669 (2000).10.1029/1999JE001203
23. Schenk PM Thickness constraints on the icy shells of the galilean satellites from a comparison of crater shapes Nature 2002 417 419 421 10.1038/417419a 12024207
Schenk, P. M. Thickness constraints on the icy shells of the galilean satellites from a comparison of crater shapes. Nature 417, 419–421 (2002).12024207 10.1038/417419a
24. Bray VJ Schenk PM Melosh HJ Morgan JV Collins GS Ganymede crater dimensions–Implications for central peak and central pit formation and development Icarus 2012 217 115 129 10.1016/j.icarus.2011.10.004
Bray, V. J., Schenk, P. M., Melosh, H. J., Morgan, J. V. & Collins, G. S. Ganymede crater dimensions–Implications for central peak and central pit formation and development. Icarus 217, 115–129 (2012).10.1016/j.icarus.2011.10.004
25. Bray VJ Collins GS Morgan JV Melosh HJ Schenk PM Hydrocode simulation of Ganymede and Europa cratering trends–How thick is Europa’s crust? Icarus 2014 231 394 406 10.1016/j.icarus.2013.12.009
Bray, V. J., Collins, G. S., Morgan, J. V., Melosh, H. J. & Schenk, P. M. Hydrocode simulation of Ganymede and Europa cratering trends–How thick is Europa’s crust?. Icarus 231, 394–406 (2014).10.1016/j.icarus.2013.12.009
26. Cox R Bauer AW Impact breaching of Europa’s ice: Constraints from numerical modeling J. Geophys. Res. Planets 2015 120 1708 1719 10.1002/2015JE004877
Cox, R. & Bauer, A. W. Impact breaching of Europa’s ice: Constraints from numerical modeling. J. Geophys. Res. Planets 120, 1708–1719 (2015).10.1002/2015JE004877
27. Silber EA Johnson BC Impact crater morphology and the structure of Europa’s ice shell J. Geophys. Res. Planets 2017 122 2685 2701 10.1002/2017JE005456
Silber, E. A. & Johnson, B. C. Impact crater morphology and the structure of Europa’s ice shell. J. Geophys. Res. Planets 122, 2685–2701 (2017).10.1002/2017JE005456
28. Murchie S Head J Global reorientation and its effect on tectonic patterns on Ganymede Geophys. Res. Lett. 1986 13 345 348 10.1029/GL013i004p00345
Murchie, S. & Head, J. Global reorientation and its effect on tectonic patterns on Ganymede. Geophys. Res. Lett. 13, 345–348 (1986).10.1029/GL013i004p00345
29. Ojakangas GW Stevenson DJ Polar wander of an ice shell on Europa Icarus 1989 81 242 270 10.1016/0019-1035(89)90053-5
Ojakangas, G. W. & Stevenson, D. J. Polar wander of an ice shell on Europa. Icarus 81, 242–270 (1989).10.1016/0019-1035(89)90053-5
30. Rubincam D Polar wander on Triton and Pluto due to volatile migration Icarus 2003 163 469 478 10.1016/S0019-1035(03)00080-0
Rubincam, D. Polar wander on Triton and Pluto due to volatile migration. Icarus 163, 469–478 (2003).10.1016/S0019-1035(03)00080-0
31. Schenk P Matsuyama I Nimmo F True polar wander on Europa from global-scale small-circle depressions Nature 2008 453 368 371 10.1038/nature06911 18480819
Schenk, P., Matsuyama, I. & Nimmo, F. True polar wander on Europa from global-scale small-circle depressions. Nature 453, 368–371 (2008).18480819 10.1038/nature06911
32. McKinnon, W.B., Mohit, P.S., Greenhagen, B.T., & Bland, M.T. Polar Wander on Ganymede and Callisto – A solution to the apex-antapex cratering conundrum. In 45th Lunar and Planetary Science Conference, Houston, abstract #1869 (2014).
33. Ballantyne, H.A., Asphaug, E., Denton, C.A., Emsenhuber, A., & Jutzi, M. Sputnik Planitia as an impactor remnant indicative of an ancient rocky mascon in an oceanless Pluto. Nature Astronomy in press (2024).
34. Willemann R Reorientation of planets with elastic lithospheres Icarus 1984 60 701 709 10.1016/0019-1035(84)90174-X
Willemann, R. Reorientation of planets with elastic lithospheres. Icarus 60, 701–709 (1984).10.1016/0019-1035(84)90174-X
35. Nimmo F Pappalardo RT Diapir-induced reorientation of Saturn's moon Enceladus Nature 2006 441 614 616 10.1038/nature04821 16738654
Nimmo, F. & Pappalardo, R. T. Diapir-induced reorientation of Saturn’s moon Enceladus. Nature 441, 614–616 (2006).16738654 10.1038/nature04821
36. Melosh HJ Planetary surface processes 2011 Cambridge University Press
Melosh, H. J. Planetary surface processes (Cambridge University Press, 2011).
37. Golombek MP Constraints on the expansion of Ganymede and the thickness of the lithosphere J. Geophys. Res. Planets 1982 87 A77 A83 10.1029/JB087iS01p00A77
Golombek, M. P. Constraints on the expansion of Ganymede and the thickness of the lithosphere. J. Geophys. Res. Planets 87, A77–A83 (1982).10.1029/JB087iS01p00A77
38. Golombek MP Banerdt WB Early thermal profiles and lithospheric strength of Ganymede from extensional tectonic features Icarus 1986 68 252 265 10.1016/0019-1035(86)90022-9
Golombek, M. P. & Banerdt, W. B. Early thermal profiles and lithospheric strength of Ganymede from extensional tectonic features. Icarus 68, 252–265 (1986).10.1016/0019-1035(86)90022-9
39. McKinnon WB Parmentier EM Burns JA Ganymede and Calisto Satellites 1986 Tucson Matthews University of Arizona Press 718 763
McKinnon, W. B. & Parmentier, E. M. Ganymede and Calisto. In Satellites (ed. Burns, J. A.) 718–763 (Matthews University of Arizona Press, Tucson, 1986).
40. Nimmo F Pappalardo RT Furrow flexure and ancient heat flux on Ganymede Geophys. Res. Lett. 2004 31 L19701 10.1029/2004GL020763
Nimmo, F. & Pappalardo, R. T. Furrow flexure and ancient heat flux on Ganymede. Geophys. Res. Lett. 31, L19701 (2004).10.1029/2004GL020763
41. Hirata N Rotational effect as the possible cause of the east-west asymmetric crater rims on Ryugu observed by LIDAR data Icarus 2021 354 114073 10.1016/j.icarus.2020.114073
Hirata, N. et al. Rotational effect as the possible cause of the east-west asymmetric crater rims on Ryugu observed by LIDAR data. Icarus 354, 114073 (2021).10.1016/j.icarus.2020.114073
42. Hirata N Secondary cratering from Rheasilvia as the possible origin of Vesta's equatorial troughs J. Geophys. Res. Planets 2023 128 E2022JE007473 10.1029/2022JE007473
Hirata, N. Secondary cratering from Rheasilvia as the possible origin of Vesta’s equatorial troughs. J. Geophys. Res. Planets 128, E2022JE007473 (2023).10.1029/2022JE007473
43. Housen KR Holsapple KA Ejecta from impact craters Icarus 2011 211 856 875 10.1016/j.icarus.2010.09.017
Housen, K. R. & Holsapple, K. A. Ejecta from impact craters. Icarus 211, 856–875 (2011).10.1016/j.icarus.2010.09.017
44. Raducan SD Davison TM Collins GS Ejecta distribution and momentum transfer from oblique impacts on asteroid surfaces Icarus 2022 374 114793 10.1016/j.icarus.2021.114793
Raducan, S. D., Davison, T. M. & Collins, G. S. Ejecta distribution and momentum transfer from oblique impacts on asteroid surfaces. Icarus 374, 114793 (2022).10.1016/j.icarus.2021.114793
45. Bierhaus EB Dones L Alvarellos JL Zahnle K The role of ejecta in the small crater populations on the mid-sized saturnian satellites Icarus 2012 218 602 621 10.1016/j.icarus.2011.12.011
Bierhaus, E. B., Dones, L., Alvarellos, J. L. & Zahnle, K. The role of ejecta in the small crater populations on the mid-sized saturnian satellites. Icarus 218, 602–621 (2012).10.1016/j.icarus.2011.12.011
46. Singer KN McKinnon WB Nowicki LT Secondary craters from large impacts on Europa and Ganymede: Ejecta size-velocity distributions on icy worlds, and the scaling of ejected blocks Icarus 2013 226 865 884 10.1016/j.icarus.2013.06.034
Singer, K. N., McKinnon, W. B. & Nowicki, L. T. Secondary craters from large impacts on Europa and Ganymede: Ejecta size-velocity distributions on icy worlds, and the scaling of ejected blocks. Icarus 226, 865–884 (2013).10.1016/j.icarus.2013.06.034
47. Arakawa M Okazaki M Nakamura M Jutzi M Yasui M Hasegawa S Dispersion and shattering strength of rocky and frozen planetesimals studied by laboratory experiments and numerical simulations Icarus 2022 373 114777 10.1016/j.icarus.2021.114777
Arakawa, M. et al. Dispersion and shattering strength of rocky and frozen planetesimals studied by laboratory experiments and numerical simulations. Icarus 373, 114777 (2022).10.1016/j.icarus.2021.114777
48. Murray CD Dermott SF Solar system dynamics 1999 New York Cambridge University Press 67
Murray, C. D. & Dermott, S. F. Solar system dynamics 67 (Cambridge University Press, New York, 1999).
49. Schenk, P., McKinnon, W.B., Moore, J. & Nimmo, F. The topography of Ganymede (and Callisto): Geology, global characteristics and future exploration. In 52nd Lunar and Planetary Science Conference, Houston, abstract #2228 (2021).
50. Ravine MA Ganymede observations by JunoCam on Juno Perijove 34 Geophys. Res. Lett. 2022 49 E2022GL099211 10.1029/2022GL099211 37034393
Ravine, M. A. et al. Ganymede observations by JunoCam on Juno Perijove 34. Geophys. Res. Lett. 49, E2022GL099211 (2022).37034393 10.1029/2022GL099211
51. Zahnle K Schenk P Sobieszczyk S Dones L Levison HF Differential cratering of synchronously rotating satellites by ecliptic comets Icarus 2001 153 111 129 10.1006/icar.2001.6668
Zahnle, K., Schenk, P., Sobieszczyk, S., Dones, L. & Levison, H. F. Differential cratering of synchronously rotating satellites by ecliptic comets. Icarus 153, 111–129 (2001).10.1006/icar.2001.6668
52. Schultz, P. H., Gault, D. E. & Crawford, D. Impacts of hemispherical granular targets: Implications for global impacts. In Reports of Planetary Geology and Geophysics Program, pp. 380–381 (NASA, 1986).
53. Li J-Y Ejecta from the DART-produced active asteroid Dimorphos Nature 2023 616 452 456 10.1038/s41586-023-05811-4 36858074
Li, J.-Y. et al. Ejecta from the DART-produced active asteroid Dimorphos. Nature 616, 452–456 (2023).36858074 10.1038/s41586-023-05811-4
54. Palguta J Anderson JD Schubert G Moore WB Mass anomalies on Ganymede Icarus 2006 180 428 441 10.1016/j.icarus.2005.08.020
Palguta, J., Anderson, J. D., Schubert, G. & Moore, W. B. Mass anomalies on Ganymede. Icarus 180, 428–441 (2006).10.1016/j.icarus.2005.08.020
55. Bland MT Singer KN McKinnon WB Schenk PM Viscous relaxation of Ganymede’s impact craters: Constraints on heat flux Icarus 2017 296 275 288 10.1016/j.icarus.2017.06.012
Bland, M. T., Singer, K. N., McKinnon, W. B. & Schenk, P. M. Viscous relaxation of Ganymede’s impact craters: Constraints on heat flux. Icarus 296, 275–288 (2017).10.1016/j.icarus.2017.06.012
56. Burkhard LML Thomas N Comment on ‘A global system of furrows on Ganymede indicative of their creation in a single impact event’ by Hirata, N., Suetsugu, R. and Ohtsuki, K. (Icarus 352, 2020): Investigating the influence of the hypothesized ancient impactor on Ganymede's orbital eccentricity Icarus 2024 416 116088 10.1016/j.icarus.2024.116088
Burkhard, L. M. L. & Thomas, N. Comment on ‘A global system of furrows on Ganymede indicative of their creation in a single impact event’ by Hirata, N., Suetsugu, R. and Ohtsuki, K. (Icarus 352, 2020): Investigating the influence of the hypothesized ancient impactor on Ganymede’s orbital eccentricity. Icarus 416, 116088 (2024).10.1016/j.icarus.2024.116088
57. Hirata N Differential impact cratering of Saturn's satellites by heliocentric impactors J. Geophys. Res. Planets 2016 121 111 117 10.1002/2015JE004940
Hirata, N. Differential impact cratering of Saturn’s satellites by heliocentric impactors. J. Geophys. Res. Planets 121, 111–117 (2016).10.1002/2015JE004940
58. Xu L Hirata N Miyamoto H Ray craters on Ganymede: Implications for cratering apex-antapex asymmetry and surface modification processes Icarus 2017 295 140 148 10.1016/j.icarus.2017.06.003
Xu, L., Hirata, N. & Miyamoto, H. Ray craters on Ganymede: Implications for cratering apex-antapex asymmetry and surface modification processes. Icarus 295, 140–148 (2017).10.1016/j.icarus.2017.06.003
59. Anderson JD Schubert G Jacobson RA Lau EL Moore WB Palguta JL Discovery of mass anomalies on Ganymede Science 2004 305 5686 989 991 10.1126/science.1099050 15310898
Anderson, J. D. et al. Discovery of mass anomalies on Ganymede. Science 305(5686), 989–991 (2004).15310898 10.1126/science.1099050
60. Zubarev A Nadezhdina I Oberst J Hussmann H Stark A New Ganymede control point network and global shape model Planet. Space Sci 2015 117 246 249 10.1016/j.pss.2015.06.022
Zubarev, A., Nadezhdina, I., Oberst, J., Hussmann, H. & Stark, A. New Ganymede control point network and global shape model. Planet. Space Sci 117, 246–249 (2015).10.1016/j.pss.2015.06.022
61. Casajus LG Ermakov AI Zannoni M Keane JT Stevenson D Buccino DR Durante D Parisi M Park RS Tortora P Bolton SJ Gravity field of ganymede after the juno extended mission Geophys. Res. Lettt. 2022 49 e2022GL099475 10.1029/2022GL099475
Casajus, L. G. et al. Gravity field of ganymede after the juno extended mission. Geophys. Res. Lettt. 49, e2022GL099475 (2022).10.1029/2022GL099475
62. Kaula WM An introduction to planetary physics: The terrestrial planets 1968 Hoboken Wiley
Kaula, W. M. An introduction to planetary physics: The terrestrial planets (Wiley, Hoboken, 1968).
63. Turcotte DL Willemann RJ Haxby WF Norberry J Role of membrane stresses in the support of planetary topography J. Geophys. Res. Solid Earth 1981 86 3951 3959 10.1029/JB086iB05p03951
Turcotte, D. L., Willemann, R. J., Haxby, W. F. & Norberry, J. Role of membrane stresses in the support of planetary topography. J. Geophys. Res. Solid Earth 86, 3951–3959 (1981).10.1029/JB086iB05p03951
64. Schulson EM Duval P Creep and fracture of ice 2009 Cambridge Cambridge University Press 401
Schulson, E. M. & Duval, P. Creep and fracture of ice 401 (Cambridge University Press, Cambridge, 2009).
65. Feistel R Wagner W A new equation of state for H2O ice Ih J. Phys. Chem. Ref. Data 2006 35 1021 1047 10.1063/1.2183324
Feistel, R. & Wagner, W. A new equation of state for H2O ice Ih. J. Phys. Chem. Ref. Data 35, 1021–1047 (2006).10.1063/1.2183324
66. Harriss KH Burchell MJ Hypervelocity impacts into ice-topped layered targets: Investigating the effects of ice crust thickness and subsurface density on crater morphology Meteor. Planet. Sci. 2017 52 1505 1522 10.1111/maps.12913
Harriss, K. H. & Burchell, M. J. Hypervelocity impacts into ice-topped layered targets: Investigating the effects of ice crust thickness and subsurface density on crater morphology. Meteor. Planet. Sci. 52, 1505–1522 (2017).10.1111/maps.12913
67. Barr AC Canup RM Origin of the Ganymede-Callisto dichotomy by impacts during the late heavy bombardment Nat. Geosci. 2010 3 164 10.1038/ngeo746
Barr, A. C. & Canup, R. M. Origin of the Ganymede-Callisto dichotomy by impacts during the late heavy bombardment. Nat. Geosci. 3, 164 (2010).10.1038/ngeo746
68. Turcotte DL Schubert G In Geodynamics 2002 2 Cambridge University Press
Turcotte, D. L. & Schubert, G. In Geodynamics 2nd edn. (Cambridge University Press, 2002).
69. Andersson O Inaba A Thermal conductivity of crystalline and amorphous ices and its implications on amorphization and glassy water Phys. Chem. Chem. Phys. 2005 7 1441 1449 10.1039/b500373c 19787966
Andersson, O. & Inaba, A. Thermal conductivity of crystalline and amorphous ices and its implications on amorphization and glassy water. Phys. Chem. Chem. Phys. 7, 1441–1449 (2005).19787966 10.1039/b500373c
70. Shulman LM The heat capacity of water ice in interstellar or interplanetary conditions Astron. Astrophys 2004 416 187 190 10.1051/0004-6361:20031746
Shulman, L. M. The heat capacity of water ice in interstellar or interplanetary conditions. Astron. Astrophys 416, 187–190 (2004).10.1051/0004-6361:20031746
71. Johnson BC Bowling TJ Trowbridge AJ Freed AM Formation of the Sputnik planum basin and the thickness of Pluto's subsurface ocean Geophys. Res. Lett. 2016 43 10068 10077 10.1002/2016GL070694
Johnson, B. C., Bowling, T. J., Trowbridge, A. J. & Freed, A. M. Formation of the Sputnik planum basin and the thickness of Pluto’s subsurface ocean. Geophys. Res. Lett. 43, 10068–10077 (2016).10.1002/2016GL070694
72. Moore WB Schubert G The tidal response of Ganymede and Callisto with and without liquid water oceans Icarus 2003 166 223 226 10.1016/j.icarus.2003.07.001
Moore, W. B. & Schubert, G. The tidal response of Ganymede and Callisto with and without liquid water oceans. Icarus 166, 223–226 (2003).10.1016/j.icarus.2003.07.001
73. Kamata S Matsuyama I Nimmo F Tidal resonance in icy satellites with subsurface oceans J. Geophys. Res. Planets 2015 120 1528 1542 10.1002/2015JE004821
Kamata, S., Matsuyama, I. & Nimmo, F. Tidal resonance in icy satellites with subsurface oceans. J. Geophys. Res. Planets 120, 1528–1542 (2015).10.1002/2015JE004821
74. Kamata S Kimura J Matsumoto K Nimmo F Kuramoto K Namiki N Tidal deformation of Ganymede: Sensitivity of Love numbers on the interior structure J. Geophys. Res. Planets 2016 121 1362 1375 10.1002/2016JE005071
Kamata, S. et al. Tidal deformation of Ganymede: Sensitivity of Love numbers on the interior structure. J. Geophys. Res. Planets 121, 1362–1375 (2016).10.1002/2016JE005071
75. Kamata S Poroviscoelastic gravitational dynamics J. Geophys. Res. Planets 2023 128 e2022JE007700 10.1029/2022JE007700
Kamata, S. Poroviscoelastic gravitational dynamics. J. Geophys. Res. Planets 128, e2022JE007700 (2023).10.1029/2022JE007700
76. Peale SJ Burns JA Rotation Histories of the Natural Satellites Planetary Satellites 1977 University of Arizona Press 87 112
Peale, S. J. Rotation Histories of the Natural Satellites. In Planetary Satellites (ed. Burns, J. A.) 87–112 (University of Arizona Press, 1977).
77. Noyelles B Karatekin Ö Rambaux N The rotation of Mimas Astron. Astrophys 2011 536 A61 10.1051/0004-6361/201117558
Noyelles, B., Karatekin, Ö. & Rambaux, N. The rotation of Mimas. Astron. Astrophys 536, A61 (2011).10.1051/0004-6361/201117558
78. Goldreich P Soter S Q in the solar system Icarus 1966 5 375 389 10.1016/0019-1035(66)90051-0
Goldreich, P. & Soter, S. Q in the solar system. Icarus 5, 375–389 (1966).10.1016/0019-1035(66)90051-0
79. Hussmann H Sotin C Lunine JI Schubert G Interiors and evolution of icy satellites Treatise on Geophysics 2015 2 Amsterdam Elsevier 605 635
Hussmann, H., Sotin, C. & Lunine, J. I. Interiors and evolution of icy satellites. In Treatise on Geophysics 2nd edn (ed. Schubert, G.) 605–635 (Elsevier, Amsterdam, 2015).
80. Singer KN Bland MT Schenk PM McKinnon WB Relaxed impact craters on Ganymede: Regional variation and high heat flows Icarus 2018 306 214 224 10.1016/j.icarus.2018.01.012
Singer, K. N., Bland, M. T., Schenk, P. M. & McKinnon, W. B. Relaxed impact craters on Ganymede: Regional variation and high heat flows. Icarus 306, 214–224 (2018).10.1016/j.icarus.2018.01.012
81. Schenk PM Ridolfi FJ Morphology and scaling of ejecta deposits on icy satellites Geophys. Res. Lett. 2002 29 1590 10.1029/2001GL013512
Schenk, P. M. & Ridolfi, F. J. Morphology and scaling of ejecta deposits on icy satellites. Geophys. Res. Lett. 29, 1590 (2002).10.1029/2001GL013512
82. Greeley R Klemaszewski JE Wagner R Galileo views of the geology of Callisto Planet. Space Sci. 2000 48 9 829 853 10.1016/S0032-0633(00)00050-7
Greeley, R., Klemaszewski, J. E. & Wagner, R. Galileo views of the geology of Callisto. Planet. Space Sci. 48(9), 829–853 (2000).10.1016/S0032-0633(00)00050-7
83. Downey BG Nimmo F Matsuyama I Inclination damping on Callisto Monthly Notices Royal Astronom. Soc. 2020 499 40 51 10.1093/mnras/staa2802
Downey, B. G., Nimmo, F. & Matsuyama, I. Inclination damping on Callisto. Monthly Notices Royal Astronom. Soc. 499, 40–51 (2020).10.1093/mnras/staa2802
84. Gault DE Wedekind JA Experimental studies of oblique impact Lunar Planet. Sci. Conf. Proc. 1978 9 3843 3875
Gault, D. E. & Wedekind, J. A. Experimental studies of oblique impact. Lunar Planet. Sci. Conf. Proc. 9, 3843–3875 (1978).
85. Herrick RR Forsberg-Taylor NK The shape and appearance of craters formed by oblique impact on the Moon and Venus Meteor. Planet. Sci. 2003 38 1551 1578 10.1111/j.1945-5100.2003.tb00001.x
Herrick, R. R. & Forsberg-Taylor, N. K. The shape and appearance of craters formed by oblique impact on the Moon and Venus. Meteor. Planet. Sci. 38, 1551–1578 (2003).10.1111/j.1945-5100.2003.tb00001.x
86. Herrick RR Hessen KK The planforms of low-angle impact craters in the northern hemisphere of Mars Meteor. Planet. Sci. 2006 41 1483 1495 10.1111/j.1945-5100.2006.tb00431.x
Herrick, R. R. & Hessen, K. K. The planforms of low-angle impact craters in the northern hemisphere of Mars. Meteor. Planet. Sci. 41, 1483–1495 (2006).10.1111/j.1945-5100.2006.tb00431.x
87. Krishna N Kumar PS Impact spallation processes on the Moon: A case study from the size and shape analysis of ejecta boulders and secondary craters of Censorinus crater Icarus 2016 264 274 299 10.1016/j.icarus.2015.09.033
Krishna, N. & Kumar, P. S. Impact spallation processes on the Moon: A case study from the size and shape analysis of ejecta boulders and secondary craters of Censorinus crater. Icarus 264, 274–299 (2016).10.1016/j.icarus.2015.09.033
88. Luo X-Z Zhu M-H Ding M Ejecta pattern of oblique impacts on the moon from numerical simulations J. Geophys. Res. Planets 2022 127 E2022JE007333 10.1029/2022JE007333
Luo, X.-Z., Zhu, M.-H. & Ding, M. Ejecta pattern of oblique impacts on the moon from numerical simulations. J. Geophys. Res. Planets 127, E2022JE007333 (2022).10.1029/2022JE007333
