P076-0010
Rheology and Heat Transport in Europa's Ice Shell

Wednesday, 16 December 2020
Poster
Ionawr Munhoz-Boillot1,2, Falko Schulz1 and Ana-Catalina Plesa1, (1)German Aerospace Center (DLR), Berlin, Germany, (2)Technische Universität Berlin, Berlin, Germany
Abstract:
The efficiency of heat transfer in the outer shell of icy satellites is important to determine the evolution and thermal state of their interior with major implications for the cooling behavior of an internal ocean.

In this study, we systematically investigate thermal convection in the ice shell of Europa using an Arrhenius viscosity that takes into account the rheology of ice I material that is dependent on both grain size and strain rate. To this end, we employ the geodynamical code GAIA [1] with a mixed rheology approach similar to [2], and perform calculations in a 2D Cartesian box and spherical annulus geometry for two values of the ice shell thickness (i.e., 30 and 70 km). In our simulations, we test various constant grain size values (Fig. 1).

In a first series of simulations, we tested the importance of the dislocation creep mechanism for modeling convection in Europa’s ice shell. Our results show that, in a mixed diffusion-dislocation creep rheology, diffusion creep is the dominant heat transfer mechanism similar to the study of [3]. A pure dislocation creep rheology leads to a conductive ice shell. Dislocation creep may become dominant if its rheological prefactor increases by about 5 orders of magnitude, which, even taking into account the uncertainty associated with rheological measurements, is considered unrealistic.

Additional simulations that use a mixed diffusion-basal slip rheology show that, for ice shells, basal slip may be a relevant deformation mechanism in addition to diffusion creep (Fig. 1).

Another important aspect is that the Nusselt number and hence the efficiency of heat transfer is larger for a thick ice shell (70 km, Fig. 1a) compared to a thinner one (i.e., 30 km). However, the dimensional surface heat flow obtained for a thin ice shell is larger than for a thicker one, as shown in Fig. 1b. This is caused by the rescaling of the Nusselt number to a dimensional heat flow.

References: [1] Hüttig et al., PEPI 2013; [2] Schulz et al., GJI 2019; [3] Harel et al., Icarus 2020.