Numerical modeling the formation of impact craters: Implications for the structure of Europa's ice shell

Craters produced by hypervelocity impacts are an invaluable tool for studying planetary surfaces. The observed impact crater depth-diameter (d-D) on the Galilean moon Europa exhibits three distinct transition regimes, two of which may correspond to the presence of warm convecting ice at depths of 7-8 km and a liquid ocean at 19-25 km, respectively [1]. In our study, we use iSALE2D to model formation of impact craters on Europa to investigate thickness and internal structure of its ice shell. This study is different from previous modeling studies [2,3] in that we consider the both fully conductive ice shell over ocean, as well as conductive lid overlying warm convecting ice, to discern the boundary conditions at the interface between the ice and the underlying ocean. Moreover, our model includes implementation of the full viscoelastic-plastic rheology for ice. Our results suggest that both conductive shell over ocean and conductive lid over warm convective ice are equally probable on Europa. We will discuss the implications and relevance of these results. The plausible scenarios are either a 6 - 7 km thick conductive ice lid overlying warm convecting ice at 265 K, or an 8 km completely conductive ice shell over ocean. Acknowledgements: We gratefully acknowledge the developers of iSALE-2D (www.isale-code.de), the simulation code used in our research, including G. Collins, K. Wünnermann, D. Elbeshausen, B. Ivanov and J. Melosh. References: [1] Schenk P. (2002) Nature, 417, 419-421. [2] Bray V.J. et al. (2014) Icarus, 231, 394-406. [3] Cox R. and Beuer A.W. (2015) JGR - Planets, 120(10), 1708-1719.