Thermal convection in a spherical shell with melting/freezing at either or both of its boundaries

Abstract

In a number of geophysical or planetological settings, including Earth’s inner core, a silicate mantle crystallizing from a magma ocean, or an ice shell surrounding a deep water ocean—a situation possibly encountered in a number of Jupiter and Saturn’s icy satellites—a convecting crystalline layer is in contact with a layer of its melt. Allowing for melting/freezing at one or both of the boundaries of the solid layer is likely to affect the pattern of convection in the layer. We study here the onset of thermal convection in a viscous spherical shell with dynamically induced melting/freezing at either or both of its boundaries. It is shown that the behavior of each interface—permeable or impermeable—depends on the value of a dimensionless number P (one for each boundary), which is the ratio of a melting/freezing timescale over a viscous relaxation timescale. A small value of P corresponds to permeable boundary conditions, while a large value of P corresponds to impermeable boundary conditions. Linear stability analysis predicts a significant effect of semi-permeable boundaries when the number P characterizing either of the boundary is small enough: allowing for melting/freezing at either of the boundary allows the emergence of larger scale convective modes. The effect is particularly drastic when the outer boundary is permeable, since the degree 1 mode remains the most unstable even in the case of thin spherical shells. In the case of a spherical shell with permeable inner and outer boundaries, the most unstable mode consists in a global translation of the solid shell, with no deformation. In the limit of a full sphere with permeable outer boundary, this corresponds to the “convective translation” mode recently proposed for Earth’s inner core. As another example of possible application, we discuss the case of thermal convection in Enceladus’ ice shell assuming the presence of a global subsurface ocean, and found that melting/freezing could have an important effect on the pattern of convection in the ice shell.

Keywords

planetary mantle / thermal convection / phase change / linear stability analysis

References

[1] Abramovich M, Stegun I A. Handbook of Mathematical Functions, 1965 Washington DC: US Government Printing Office

[2] Alboussière T, Deguen R, Melzani M. Melting Induced Stratification above the Earth’s Inner Core due to Convective Translation. Nature, 2010, 466: 744-747.

[3] Chandrasekhar S. Hydrodynamic and Hydromagnetic Stability. International Series of Monographs on Physics, 1961

[4] Deguen R. Structure and Dynamics of Earth’s Inner Core. Earth Planet. Sci. Lett., 2012, 333–334: 211-225.

[5] Deguen R, Alboussière T, Cardin P. Thermal Convection in Earth’s Inner Core with Phase Change at Its Boundary. Geophys. J. Int., 2013

[6] Grott M, Sohl F, Hussmann H. Degree-One Convection and the Origin of Enceladus’ Dichotomy. Icarus, 2007, 191(1): 203-210.

[7] Irving J C E, Deuss A. Hemispherical Structure in Inner Core Velocity Anisotropy. Journal of Geophysical Research, 2011, 116 B4 B04307

[8] Kivelson M G, Khurana K K, Russell C T, . Galileo Magnetometer Measurements: A Stronger Case for a Subsurface Ocean at Europa. Science, 2000, 289(5483): 1340-1343.

[9] Labrosse S, Hernlund J W, Coltice N. A Crystallizing Dense Magma Ocean at the Base of the Earth’s Mantle. Nature, 2007, 450(7171): 866-869.

[10] McNamara A K, Zhong S. Degree-One Mantle Convection: Dependence on Internal Heating and Temperature-Dependent Rheology. Geophysical Research Letters, 2005, 32 1 L01301

[11] Mizzon H, Monnereau M. Implication of the Lopsided Growth for the Viscosity of Earth’s Inner Core. Earth Planet. Sci. Lett., 2013, 361: 391-401.

[12] Monnereau M, Calvet M, Margerin L, . Lopsided Growth of Earth’s Inner Core. Science, 2010, 328: 1014-1017.

[13] Monnereau M, Dubuffet F. Is Io’s Mantle Really Mol ten?. Icarus, 2002, 158(2): 450-459.

[14] Nimmo F, Pappalardo R T. Diapir-Induced Reorientation of Saturn’s Moon Enceladus. Nature, 2006, 441(7093): 614-616.

[15] Niu F L, Wen L X. Hemispherical Variations in Seismic Velocity at the Top of the Earth’s Inner Core. Nature, 2001, 410: 1081-1084.

[16] Porco C C, Helfenstein P, Thomas P C, . Cassini Observes the Active South Pole of Enceladus. Science, 2006, 311(5766): 1393-1401.

[17] Ribe N M. Analytical Approaches to Mantle Dynamics. Treatise on Geophysics, 2007, 7: 167-226.

[18] Schubert G, Anderson J D, Travis B J, . Enceladus: Present Internal Structure and Differentiation by Early and Long-Term Radiogenic Heating. Icarus, 2007, 188(2): 345-355.

[19] Solomatov V S. Fluid Dynamics of a Terrestrial Magma Ocean. Origin of the Earth and Moon, 2000, 1: 323-338.

[20] Spohn T, Schubert G. Oceans in the Icy Galilean Satellites of Jupiter?. Icarus, 2003, 161(2): 456-467.

[21] Stegman D R, Freeman J, May D A. Origin of Ice Diapirism, True Polar Wander, Subsurface Ocean, and Tiger Stripes of Enceladus Driven by Compositional Convection. Icarus, 2009, 202(2): 669-680.

[22] Tanaka S, Hamaguchi H. Degree One Heterogeneity and Hemispherical Variation of Anisotropy in the Inner Core from PKP(BC)-PKP(DF) Times. Journal of Geophysical Research, 1997, 102: 2925-2938.

[23] Tyler R H. Ocean Tides Heat Enceladus. Geophysical Research Letters, 2009, 36 15 L15205

[24] Tyler R H. Strong Ocean Tidal Flow and Heating on Moons of the Outer Planets. Nature, 2008, 456(7223): 770-772.

[25] Ulvrová M, Labrosse S, Coltice N, . Numerical Modelling of Convection Interacting with a Melting and Solidification Front: Application to the Thermal Evolution of the Basal Magma Ocean. Physics of the Earth and Planetary Interiors, 2012, 206(207): 51-66.

[26] Waite J H Jr., Lewis W S, Magee B A, . Liquid Water on Enceladus from Observations of Ammonia and 40ar in the Plume. Nature, 2009, 460(7254): 487-490.

[27] Zhong S, Zuber M T. Degree-1 Mantle Convection and the Crustal Dichotomy on Mars. Earth Planet. Sci. Lett., 2001, 189(1): 75-84.