A youthful Titan implied by improved impact simulations

The small number of impact craters found on Titan suggests that its surface is relatively young. Previous work estimated its surface age to be between 200 and 1000 Myr. This estimate, however, is based on crater scaling laws for water and sand, which are not representative of the composition of Titan’s icy surface. Titan’s surface is likely composed of water ice, methane clathrates, or a combination of both. Here, we perform impact simulations for impactors of various sizes that strike an icy target with a 0-15 km thick methane clathrate cap layer. We derive new crater scaling laws based on our numerical results, and find that Titan’s surface age is 300-340 Myr, assuming heliocentric impactors and surface clathrates. This age, which represents the crater retention age, indicates a relatively youthful surface, suggesting that active endogenic and/or exogenic processes have recently reshaped Titan’s surface.

💡 Research Summary

Titan, Saturn’s largest moon, exhibits an unusually low density of impact craters—only about 90 identifiable features despite its size and age. Previous estimates of the surface (crater‑retention) age, ranging from 200 to 1000 Myr, relied on crater scaling laws derived for water or sand targets, which do not accurately represent Titan’s icy crust that likely contains a mixture of water ice and methane‑clathrate (a solid that traps methane molecules). Recognizing this mismatch, the authors performed a series of high‑resolution impact simulations using the iSALE‑2D shock physics code to quantify how a methane‑clathrate cap layer modifies crater formation.

Key methodological points:

• Impact velocity was set to 10.5 km s⁻¹, the mean heliocentric impact speed on Titan, with vertical impacts simulated and later corrected for the most probable 45° incidence angle (multiplying by sin 0.38).
• Impactor diameters spanned 1.5 km (the smallest size that survives atmospheric entry) to 25 km (the size inferred for the largest observed crater, Menrva).
• Four target configurations were examined: a pure water‑ice crust (0 km clathrate) and three cases with methane‑clathrate caps of 5, 10, and 15 km thickness. Thermal profiles for each cap thickness were taken from Kalousová & Sotin (2020), assuming 1 mm ice grains and conductive heat transport.
• Material strength models incorporated temperature‑dependent behavior for both water ice and methane‑clathrate, based on laboratory data indicating that clathrate is 20–30 × stronger than ice at 250–287 K but weakens sharply near its dissociation temperature.
• Because no dedicated equation of state (EOS) exists for methane‑clathrate, the ANEOS for water‑ice was used for both materials, a common practice in prior Titan impact studies.
• For impactors larger than 5 km, a subsurface ocean 100 km deep was added to assess its influence; sensitivity tests showed that EOS choice and ocean presence altered crater diameters by ≤10 %, contributing less than 2 % uncertainty to the final age estimate.

Results revealed a clear “jump” in crater size when the transient crater breaches the clathrate layer. Small impactors (< 2.5 km) produce craters confined within the clathrate cap, so the strength contrast between clathrate and ice has little effect. For larger impactors (≈ 2.5–4 km), the insulating nature of the clathrate raises temperatures at depth, weakening the underlying ice and yielding final crater diameters about 1.7 × larger than in a pure‑ice target. Very large impactors (> 20 km) excavate to depths where both scenarios encounter similarly weak ice, so crater size becomes insensitive to the cap thickness.

From the simulated data, the authors derived four empirical scaling relationships (log D_crater = a log d_imp + b) specific to each clathrate thickness. These new laws sit between the classic water‑ice and sand scaling curves (Artemieva & Lunine 2005; Korycansky & Zahnle 2005) but predict noticeably larger diameters for craters ≤ 50 km, reflecting the thermal‑strength effect of the clathrate layer.

To translate crater size into surface age, the study adopted the heliocentric impactor flux model of Zahnle et al. (2003) unchanged, but substituted the new scaling laws when converting observed crater diameters (≥ 20 km) to impactor sizes. Matching the observed 90‑crater inventory yields a crater‑retention age of 300–340 Myr, provided a methane‑clathrate cap is present. This age is at the younger end of earlier estimates and suggests that Titan’s surface has been actively reshaped within the last few hundred million years.

Uncertainties stem primarily from the assumed clathrate thickness and its strength model; variations in these parameters could shift the age by roughly ±10 %. Nonetheless, sensitivity analyses indicate that EOS selection, ocean inclusion, and neglect of porosity each contribute < 2 % to the total error budget.

The implications are significant: a relatively youthful surface points to ongoing endogenic processes (e.g., internal heating, cryovolcanism, subsurface ocean dynamics) and/or vigorous exogenic activity (continuous methane precipitation, clathrate formation, and erosion). Moreover, the findings underscore the necessity of incorporating realistic target material properties when interpreting crater statistics on icy worlds. Future work should aim to obtain laboratory measurements of methane‑clathrate behavior under Titan‑relevant pressures and temperatures, and to acquire higher‑resolution radar or lidar data that could directly detect clathrate‑controlled geomorphology.

In summary, by integrating methane‑clathrate physics into impact simulations, the authors provide a refined crater scaling framework and a more constrained, younger surface age for Titan, reinforcing the view that this moon remains geologically active on a timescale of a few hundred million years.