Stagnant lid convection with temperature-dependent thermal conductivity and the thermal evolution of icy worlds

SUMMARY Convection is an efficient process to release heat from planetary interiors, but its efficiency depends on the detailed properties of planetary mantles and materials. A property whose impact has not yet been studied extensively is the temperature dependence of thermal conductivity. Because thermal conductivity controls heat fluxes, its variations with temperature may alter heat transfer. Here, I assess qualitatively and quantitatively the influence of temperature-dependent thermal conductivity on stagnant lid convection. Assuming that thermal conductivity varies as the inverse of temperature $(k \propto 1/T)$, which is the case for ice Ih, the main component of outer shells of solar System large icy bodies, I performed numerical simulations of convection in 3-D-Cartesian geometry with top-to-bottom viscosity and conductivity ratios in the ranges 105 ≤ Δη ≤ 108 and 1 ≤ Rk ≤ 10, respectively. These simulations indicate that with increasing Rk, and for given values of the Rayleigh number and Δη, heat flux is reduced by a factor Rk0.82, while the stagnant lid is thickening. These results have implications for the structures and thermal evolutions of large icy bodies, the impact of temperature-dependent conductivity being more important with decreasing surface temperature, Tsurf. The heat fluxes and thermal evolutions obtained with temperature-dependent conductivity are comparable to those obtained with constant conductivity, provided that the conductivity is fixed to its value at the bottom or in the interior of the ice shell, that is, around 2.0–3.0 W m−1 K−1, depending on the body. By contrast, temperature-dependent conductivity leads to thicker stagnant lids, by about a factor 1.6–1.8 at Pluto (Tsurf = 40 K) and a factor 1.2–1.4 at Europa (Tsurf = 100 K), and smaller interior temperatures. Overall, temperature-dependent thermal conductivity therefore provides more accurate descriptions of the thermal evolutions of icy bodies.