Breaching the Ice: The Role of Impact Cratering in Facilitating Surface-Ocean Exchange on Europa

Europa is one of the prime targets for astrobiological exploration in upcoming years. Beneath its ice shell, Europa harbors a global subsurface ocean. Habilibility on the icy moons depends critically on the exchange of material between the oxidant-rich surface and the subsurface ocean. Impact cratering represents a primary mechanism for facilitating this ocean-surface exchange, either by directly breaching the ice shell to create transient water pathways or by fracturing and weakening the crust to facilitate surface material transport.

Full impact penetration events and subsequent impact-induced weak zones are strongly dependent on the ice shell structure and thickness. Europa’s ice shell thickness remains controversial, with estimates ranging from 1 to 47 km [1-2]. Recent numerical modeling of Callanish and Tyre multiring impact basins [3] has suggested an ice shell thickness exceeding 20 km, with a conductive lid of 6 to 8 km.

We employ the iSALE-3D shock physics codes [4-6] to simulate crater formation and determine the criteria for full ice shell penetration. We use the tabular 5-phase Equation of State (EOS) for water-ice [7] to describe physical changes under expanding temperature and pressure,  and compare these results with the derived analytical EOS which was adopted by most of previous research. We cover a wide parameter space, ranging from vertical (90°) to highly oblique (15°) impacts, with undifferentiated stony and icy impactors (1000 - 2700kg/m³) ranging from 50 m to 6 km in diameter. The mean impactor velocity at Europa is 26 km/s, and we include a broader impactor velocity range (5 - 40km/s). Furthermore, we also vary the ice shell thickness (up to 50km) to cover various geophysical scenarios. To accurately quantify material fragmentation and to trace impactor material, we utilize a newly developed disruption tracking module by [8].

Our results demonstrate that the threshold for full penetration, and thus the direct exchange,  is heavily dependent on the impactor-to-target size ratio and impact velocity. We find that impact angle governs total amount of  melt generation and highly oblique impacts cause shear heating and frictional melting. Additionally, we quantify the amount of subsurface ocean material transported to the surface and map the distribution of resulting thermal anomalies. 

In future work, we will couple geophysical modeling with our impact-induced thermal anomaly. We will also characterize the shear failure zones and fracture networks generated by impacts. These results can provide constraints for geophysical investigations and identify local gravity anomalies which are potentially detectable by the ESA’s JUICE and NASA’s Europa Clipper missions.

References:

[1] Bray V. J. et al. (2014) Icarus, 231, 394–406.

[2] Howell S. M. (2021) Planet. Sci. J., 2, 129.

[3] Wakita S. et al. (2024) Sci. Adv., 10, eadj8455.

[4] Amsden A. A. et al. (1980) LANL Report LA-8095, 101 pp.

[5] Collins G. S. et al. (2004) Meteorit. Planet. Sci., 39, 217–231.

[6] Wünnemann K. et al. (2006) Icarus, 180, 514–527.

[7] Senft L. E. and Stewart S. T. (2008) Meteorit. Planet. Sci., 43, 1993–2013.

[8] Röhlen R. et al. (2025) Icarus, 431, 116464.