Thermal state of the Martian interior at present day as constrained by elastic lithosphere thickness estimates and recent volcanic activity

Previous studies have constrained the lithosphere at the north and south poles of Mars to be thick and cold, with elastic thicknesses of 330 to 450km [1], and >150km [2], respectively. The elastic thickness characterizes the stiffness of the lithosphere in response to loading and is directly linked to the thermal state of the lithosphere and the surface heat flow. Thus, elastic thickness estimates at the north and south poles provide crucial constraints on the present-day surface heat flow on Mars. Additional information on the present-day planetary thermal state comes from evidence of ongoing melting in the mantle, as indicated by the presence of both young lava flows in Tharsis and Elysium provinces and an active mantle plume beneath Elysium Planitia [3,4,5]. In this study we explore the thermal evolution of Mars using global 3D geodynamic models. These models improve upon our previous work [6] by including updated interior structure information from the InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission [7,8] and by considering constraints on the present-day thermal state of the planet as noted above. Thermal evolution models using the most recent crustal thickness estimates [8,9], require that the crust contains more than half of the total amount of heat producing elements (HPEs) to explain localized recent volcanic activity on Mars [8]. We find that the crustal thickness variations control the surface heat flow and the elastic thickness pattern, as well as the location of melting zones in the present-day Martian mantle. The strongest constraint for the thermal history and present-day state of the interior is given by the elastic thickness at the north pole. While at the south pole, all models show values >150km, compatible with the latest estimate [2], only a few models present an elastic thickness >300km at the north pole, with values still lower than the recent estimate of [1].  A larger elastic thickness at the north pole could indicate: 1) a northern crust less enriched in HPEs, 2) a colder lithosphere due to a weaker blanketing effect caused by a thinner or higher-conductivity crust on the northern hemisphere, 3) ongoing viscoelastic relaxation, suggesting that the observed surface deflection beneath the north polar cap is not the final one [1], or a combination thereof. In contrast to the cold lithosphere inferred for the Martian polar regions, recent volcanic activity suggests a warmer interior beneath Tharsis and Elysium provinces [3,4]. This reveals an important spatial variability in the thermal state and thickness of the Martian lithosphere. Our work shows that only a narrow range of models can match elastic thickness estimates at the polar caps and explain Mars’ recent volcanic activity, thereby providing important insights into the structure and thermal evolution of the interior.References:[1] Broquet et al., 2020. [2] Broquet et al., 2021. [3] Voigt et al., 2023. [4] Hauber et al., 2011. [5] Broquet & Andrews-Hanna, 2023. [6] Plesa et al., 2018. [7] Stähler et al., 2021. [8] Knapmeyer-Endrun et al., 2021. [9] Wieczorek et al., 2023.