Mantle flow throughout the lunar evolutionary history could be constrained by shear wave splitting if it is reliably observed in the moonquake waveforms. However, despite the long-standing availability of Apollo seismic data, it remains unclear whether the S-waves from moonquakes exhibit real splitting due to seismic anisotropy or whether apparent splitting arises instead from seismic wave propagation through isotropic heterogeneous media. To date, no systematic investigation has addressed this question. To fill this knowledge gap, we analyze direct S-wave splitting from deep moonquakes with high signal-to-noise ratios. Stacked waveforms from deep moonquakes show that the measurement results are inconsistent across different clusters at the same station and across different stations for the same cluster as well as across different frequency bands. Further analysis of unstacked individual events within the same cluster, differing in location by only a few kilometers, also exhibits large variability in the splitting parameters. Numerical simulations of seismic wave propagation through heterogeneous media demonstrate that the observed S-wave splitting of deep moonquakes can be explained by apparent splitting caused by scattering in addition to the potential intrinsic anisotropy. Therefore, even stacking waveforms of deep moonquakes do not improve the reliability of the splitting measurements. Our results suggest that there are no remnants of large-scale lateral mantle flow in the lunar mantle sampled by these deep moonquakes, or any existing anisotropy is obscured by the stronger effects of structural heterogeneity.
{"title":"Direct S-Wave Splitting From Deep Moonquakes: Real or Apparent?","authors":"Jing Shi, Jiaqi Li, Tao Wang","doi":"10.1029/2025JE009223","DOIUrl":"https://doi.org/10.1029/2025JE009223","url":null,"abstract":"<p>Mantle flow throughout the lunar evolutionary history could be constrained by shear wave splitting if it is reliably observed in the moonquake waveforms. However, despite the long-standing availability of Apollo seismic data, it remains unclear whether the <i>S</i>-waves from moonquakes exhibit real splitting due to seismic anisotropy or whether apparent splitting arises instead from seismic wave propagation through isotropic heterogeneous media. To date, no systematic investigation has addressed this question. To fill this knowledge gap, we analyze direct <i>S</i>-wave splitting from deep moonquakes with high signal-to-noise ratios. Stacked waveforms from deep moonquakes show that the measurement results are inconsistent across different clusters at the same station and across different stations for the same cluster as well as across different frequency bands. Further analysis of unstacked individual events within the same cluster, differing in location by only a few kilometers, also exhibits large variability in the splitting parameters. Numerical simulations of seismic wave propagation through heterogeneous media demonstrate that the observed <i>S</i>-wave splitting of deep moonquakes can be explained by apparent splitting caused by scattering in addition to the potential intrinsic anisotropy. Therefore, even stacking waveforms of deep moonquakes do not improve the reliability of the splitting measurements. Our results suggest that there are no remnants of large-scale lateral mantle flow in the lunar mantle sampled by these deep moonquakes, or any existing anisotropy is obscured by the stronger effects of structural heterogeneity.</p>","PeriodicalId":16101,"journal":{"name":"Journal of Geophysical Research: Planets","volume":"130 11","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145580752","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
S. Verkercke, L. Morrissey, J. Y. Chaufray, A. Georgiou, A. Ricketts, F. Leblanc
Mercury's exosphere is sustained by the continuous ejection of atoms from its surface, driven by solar wind, micro-meteoroid impacts, and surface heating. Due to its 3:2 spin-orbit resonance, some longitudes experience greater solar exposure, creating temperature variations from ∼90 to 700 K. This resonance also creates least exposed longitudes, called cold longitudes. These variations, combined with surface-solar interactions, lead to complex exospheric dynamics. Observations from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft revealed localized enhancements in sodium column density at cold longitudes near the aphelion, where these regions rotate into the day-side. Existing models do not explain this cold-pole enhancement but assume a smooth, impermeable surface, neglecting the highly porous nature of Mercury's regolith. This porosity allows for subsurface diffusion of volatiles through adsorption and desorption across regolith grains, influenced by temperature and species-specific surface binding energies. We couple a subsurface transport model with a 3D Monte Carlo-based Exosphere Global Model to investigate whether sodium accumulated in Mercury's regolith could explain the cold pole enhancement. Results suggest that at cold longitudes, low temperatures favor sodium retention near the surface. The gradual heating induces a release of sodium from the subsurface, producing the observed localized enhancements. This mechanism reconciles MESSENGER's findings with physical processes and highlights the significance of subsurface reservoirs in volatile dynamics. We highlight the need to consider regolith structure and subsurface processes in exosphere modeling. These results improve our understanding of Mercury's volatile cycle but also offer broader insights into the behavior of surface-bound species on airless planetary bodies.
{"title":"On the Origins of Mercury's Sodium Exospheric Cold Poles Enhancement: A Multiscale Exosphere Global Model","authors":"S. Verkercke, L. Morrissey, J. Y. Chaufray, A. Georgiou, A. Ricketts, F. Leblanc","doi":"10.1029/2025JE009309","DOIUrl":"https://doi.org/10.1029/2025JE009309","url":null,"abstract":"<p>Mercury's exosphere is sustained by the continuous ejection of atoms from its surface, driven by solar wind, micro-meteoroid impacts, and surface heating. Due to its 3:2 spin-orbit resonance, some longitudes experience greater solar exposure, creating temperature variations from ∼90 to 700 K. This resonance also creates least exposed longitudes, called cold longitudes. These variations, combined with surface-solar interactions, lead to complex exospheric dynamics. Observations from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft revealed localized enhancements in sodium column density at cold longitudes near the aphelion, where these regions rotate into the day-side. Existing models do not explain this cold-pole enhancement but assume a smooth, impermeable surface, neglecting the highly porous nature of Mercury's regolith. This porosity allows for subsurface diffusion of volatiles through adsorption and desorption across regolith grains, influenced by temperature and species-specific surface binding energies. We couple a subsurface transport model with a 3D Monte Carlo-based Exosphere Global Model to investigate whether sodium accumulated in Mercury's regolith could explain the cold pole enhancement. Results suggest that at cold longitudes, low temperatures favor sodium retention near the surface. The gradual heating induces a release of sodium from the subsurface, producing the observed localized enhancements. This mechanism reconciles MESSENGER's findings with physical processes and highlights the significance of subsurface reservoirs in volatile dynamics. We highlight the need to consider regolith structure and subsurface processes in exosphere modeling. These results improve our understanding of Mercury's volatile cycle but also offer broader insights into the behavior of surface-bound species on airless planetary bodies.</p>","PeriodicalId":16101,"journal":{"name":"Journal of Geophysical Research: Planets","volume":"130 11","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025JE009309","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145522255","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
G. Lopez-Reyes, J. A. Manrique, E. Clavé, A. Ollila, O. Beyssac, P. Pilleri, S. Bernard, E. Dehouck, M. Veneranda, S. K. Sharma, M. Nachon, J. Aramendia, O. Forni, F. Rull, T. Acosta-Maeda, A. Brown, S. M. Angel, K. Castro, E. Cloutis, L. Coloma, J. Comellas, D. Delapp, R. Jakubek, S. Julve-Gonzalez, E. Kelly, J. M. Madariaga, G. Montagnac, I. Poblacion, S. Schröder, S. A. Connell, I. Reyes-Rodriguez, Z. U. Wolf, S. Maurice, O. Gasnault, S. Clegg, A. Cousin, R. C. Wiens, the SuperCam Raman Working Group and the SuperCam Team
The Mars 2020 Perseverance rover introduced Raman spectroscopy to in situ planetary exploration for the first time when it landed in Jezero crater on Mars in February 2021. The SuperCam instrument onboard Perseverance is a multi-analytical tool capable of acquiring time-resolved Raman data from Martian targets at standoff distances of a few meters. This is a particularly challenging task due to the operational constraints, the harsh conditions on the Martian surface, and especially the very fine-grained nature of the Martian soil. To address these challenges, the SuperCam Raman team has invested significant effort into optimizing both the acquisition and post-processing of Raman data collected on Mars, as detailed in this work. Additionally, this paper reviews and discusses the detections made by SuperCam Raman during the first 1,000 sols (almost 3 Earth years) of the Mars 2020 mission. During this period, SuperCam Raman data provided key insights into the mineralogy of Jezero throughout the Crater, Delta, and Margin Campaigns. Key detections include olivine, carbonates, perchlorates, and sulfates (such as anhydrite), identified in both abraded patches and natural surfaces. The high specificity of Raman spectroscopy enables the unequivocal identification of these minerals, allowing for rapid and direct interpretation of Jezero's mineralogy, especially when combined with other techniques from SuperCam or others on the rover. Furthermore, this paper compiles the spectra acquired from the SuperCam Calibration Target samples on Mars, including studies on the degradation of the Ertalyte (PET), an organic polymer sample and analyses of diamond, apatite, and other reference materials.
{"title":"SuperCam Raman Activities at Jezero Crater, Mars: Observational Strategies, Data Processing, and Mineral Detections During the First 1000 Sols","authors":"G. Lopez-Reyes, J. A. Manrique, E. Clavé, A. Ollila, O. Beyssac, P. Pilleri, S. Bernard, E. Dehouck, M. Veneranda, S. K. Sharma, M. Nachon, J. Aramendia, O. Forni, F. Rull, T. Acosta-Maeda, A. Brown, S. M. Angel, K. Castro, E. Cloutis, L. Coloma, J. Comellas, D. Delapp, R. Jakubek, S. Julve-Gonzalez, E. Kelly, J. M. Madariaga, G. Montagnac, I. Poblacion, S. Schröder, S. A. Connell, I. Reyes-Rodriguez, Z. U. Wolf, S. Maurice, O. Gasnault, S. Clegg, A. Cousin, R. C. Wiens, the SuperCam Raman Working Group and the SuperCam Team","doi":"10.1029/2025JE008943","DOIUrl":"https://doi.org/10.1029/2025JE008943","url":null,"abstract":"<p>The Mars 2020 Perseverance rover introduced Raman spectroscopy to in situ planetary exploration for the first time when it landed in Jezero crater on Mars in February 2021. The SuperCam instrument onboard Perseverance is a multi-analytical tool capable of acquiring time-resolved Raman data from Martian targets at standoff distances of a few meters. This is a particularly challenging task due to the operational constraints, the harsh conditions on the Martian surface, and especially the very fine-grained nature of the Martian soil. To address these challenges, the SuperCam Raman team has invested significant effort into optimizing both the acquisition and post-processing of Raman data collected on Mars, as detailed in this work. Additionally, this paper reviews and discusses the detections made by SuperCam Raman during the first 1,000 sols (almost 3 Earth years) of the Mars 2020 mission. During this period, SuperCam Raman data provided key insights into the mineralogy of Jezero throughout the Crater, Delta, and Margin Campaigns. Key detections include olivine, carbonates, perchlorates, and sulfates (such as anhydrite), identified in both abraded patches and natural surfaces. The high specificity of Raman spectroscopy enables the unequivocal identification of these minerals, allowing for rapid and direct interpretation of Jezero's mineralogy, especially when combined with other techniques from SuperCam or others on the rover. Furthermore, this paper compiles the spectra acquired from the SuperCam Calibration Target samples on Mars, including studies on the degradation of the Ertalyte (PET), an organic polymer sample and analyses of diamond, apatite, and other reference materials.</p>","PeriodicalId":16101,"journal":{"name":"Journal of Geophysical Research: Planets","volume":"130 11","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2025JE008943","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145580837","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Philip Skemer, Andrew J. Cross, Bradford J. Foley, Keith D. Putirka
Earth's particular style of plate-tectonics—characterized by localized deformation along dynamic plate boundaries and long-lived stable plate interiors—appears to be unique among rocky objects in the solar system. However, it is entirely unknown how common plate tectonics and related lithospheric phenomena are among the vast population of exoplanets discovered astronomically or assumed to exist throughout the Universe. In this study, we explore the effect of planetary composition on mylonitization—a set of microphysical processes that is commonly associated with shear localization and plate boundary deformation on Earth. A model for planet compositions, based on stellar spectroscopy, is used to define a plausible range of theoretical mineral abundances in the mantles of rocky Earth-sized exoplanets. These mineral abundances, along with experimental rock rheology, are used to model microphysical evolution with two-phase mixing. The model is then used to determine the effect of composition on the time-scales for shear zone formation. We demonstrate that lithospheres composed of sub-equal proportions of two mineral phases will form shear zones over relatively short time-scales, a more favorable condition for forming Earth-like plate boundaries. In contrast, lithospheres that are nearly monomineralic may require unrealistically long time-scales to form plate boundary shear zones. Using this approach, we identify specific nearby stars with the optimal range of compositions to be targeted by future astronomical missions, including the Habitable Worlds Observatory.
{"title":"The Effect of Composition on Shear Localization in Planetary Lithospheres","authors":"Philip Skemer, Andrew J. Cross, Bradford J. Foley, Keith D. Putirka","doi":"10.1029/2025JE009106","DOIUrl":"https://doi.org/10.1029/2025JE009106","url":null,"abstract":"<p>Earth's particular style of plate-tectonics—characterized by localized deformation along dynamic plate boundaries and long-lived stable plate interiors—appears to be unique among rocky objects in the solar system. However, it is entirely unknown how common plate tectonics and related lithospheric phenomena are among the vast population of exoplanets discovered astronomically or assumed to exist throughout the Universe. In this study, we explore the effect of planetary composition on mylonitization—a set of microphysical processes that is commonly associated with shear localization and plate boundary deformation on Earth. A model for planet compositions, based on stellar spectroscopy, is used to define a plausible range of theoretical mineral abundances in the mantles of rocky Earth-sized exoplanets. These mineral abundances, along with experimental rock rheology, are used to model microphysical evolution with two-phase mixing. The model is then used to determine the effect of composition on the time-scales for shear zone formation. We demonstrate that lithospheres composed of sub-equal proportions of two mineral phases will form shear zones over relatively short time-scales, a more favorable condition for forming Earth-like plate boundaries. In contrast, lithospheres that are nearly monomineralic may require unrealistically long time-scales to form plate boundary shear zones. Using this approach, we identify specific nearby stars with the optimal range of compositions to be targeted by future astronomical missions, including the Habitable Worlds Observatory.</p>","PeriodicalId":16101,"journal":{"name":"Journal of Geophysical Research: Planets","volume":"130 11","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145522098","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sara Faggi, Shohei Aoki, Giuliano Liuzzi, Geronimo Villanueva, Hideo Sagawa, Michael J. Mumma, Ann Carine Vandaele, Frank Daerden, Sébastien Viscardy
Hydrogen chloride (HCl) was recently discovered in the Mars atmosphere using the ESA's ExoMars Trace Gas Orbiter (TGO) onboard ESA's ExoMars mission. Its discovery is the first confirmation of an active presence of any chlorine-bearing species in the modern Mars atmosphere. TGO permitted investigations of HCl altitude profiles with high precision and showed that water vapor and ice clouds play an important role in the production and temporary loss of HCl. TGO cannot always sample the Martian atmosphere near the surface, and when those measurements are possible, they are highly affected by the increase in dust opacity, nor can TGO observe at equatorial latitudes with high cadence, due to orbital constraints, so its measurements are not suitable to obtain instantaneous global coverage. In this work, we present a methodic investigation of the Martian atmosphere, in support of the ExoMars TGO mission, targeting HCl and water using iSHELL at NASA/InfraRed Telescope Facility. Our observations mapped the Martian atmosphere, exploring three seasons in Martian Year 36. We observed the beginning of an increase in the HCl abundances around LS = 249°–301°, followed by a drop in the abundances around LS = 319°. We confirmed a strong correlation between the spatial distribution of water vapor and HCl—both globally and locally—suggesting that water vapor plays an important role in the production of HCl, in agreement with previous studies. Our observations also suggest the presence of two competing processes involving the dust, one supporting HCl production and another one contributing to its destruction.
{"title":"Following the HCl Cycle Over Three Martian Seasons in Mars Year 36","authors":"Sara Faggi, Shohei Aoki, Giuliano Liuzzi, Geronimo Villanueva, Hideo Sagawa, Michael J. Mumma, Ann Carine Vandaele, Frank Daerden, Sébastien Viscardy","doi":"10.1029/2025JE009105","DOIUrl":"https://doi.org/10.1029/2025JE009105","url":null,"abstract":"<p>Hydrogen chloride (HCl) was recently discovered in the Mars atmosphere using the ESA's ExoMars Trace Gas Orbiter (TGO) onboard ESA's ExoMars mission. Its discovery is the first confirmation of an active presence of any chlorine-bearing species in the modern Mars atmosphere. TGO permitted investigations of HCl altitude profiles with high precision and showed that water vapor and ice clouds play an important role in the production and temporary loss of HCl. TGO cannot always sample the Martian atmosphere near the surface, and when those measurements are possible, they are highly affected by the increase in dust opacity, nor can TGO observe at equatorial latitudes with high cadence, due to orbital constraints, so its measurements are not suitable to obtain instantaneous global coverage. In this work, we present a methodic investigation of the Martian atmosphere, in support of the ExoMars TGO mission, targeting HCl and water using iSHELL at NASA/InfraRed Telescope Facility. Our observations mapped the Martian atmosphere, exploring three seasons in Martian Year 36. We observed the beginning of an increase in the HCl abundances around <i>L</i><sub><i>S</i></sub> = 249°–301°, followed by a drop in the abundances around <i>L</i><sub><i>S</i></sub> = 319°. We confirmed a strong correlation between the spatial distribution of water vapor and HCl—both globally and locally—suggesting that water vapor plays an important role in the production of HCl, in agreement with previous studies. Our observations also suggest the presence of two competing processes involving the dust, one supporting HCl production and another one contributing to its destruction.</p>","PeriodicalId":16101,"journal":{"name":"Journal of Geophysical Research: Planets","volume":"130 11","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145522096","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
E. R. Gehringer, G. Ferrari, A. M. Morgan, M. C. Palucis
The discovery of rounded, river-transported clasts near an alluvial fan on Mars by the Mars Science Laboratory (MSL) Curiosity rover provided the first ground-based evidence for sustained fluvial activity on Mars. However, questions remain about the paleoenvironment and provenance of these clasts. Clast shape and size are two quantitative metrics that can be used to reconstruct the history of sedimentary deposits. However, previous studies have not considered martian basaltic lithologies and thus may not be applicable to Mars. We use data from two alluvial fans in the Mojave Desert and a series of rotating drum experiments to quantify the difference between fluvial abrasion of basalt and granite to determine the effect of lithology on the relationship between clast roundness and transportation distance. We find that basalt rounds faster than granite both in the field and in the lab and that this effect is more pronounced in the field. Our rounding coefficient (pc) values show the agreement between granitic clasts in the lab (0.028 ± 0.0008) and the field (0.029 ± 0.003 and 0.034 ± 0.007) as well as the faster rounding rate we see from the field (0.048 ± 0.001) than experimental basalt (0.039 ± 0.001). This may be due to the greater susceptibility of basalt to chemical weathering. When applying our relationship to rounded clasts observed in the Gale crater, we find that our transport values agree with previous interpretations that point to a sediment provenance near the crater rim.
{"title":"Using Clast Morphometry to Reconstruct Basaltic Sediment Transport History","authors":"E. R. Gehringer, G. Ferrari, A. M. Morgan, M. C. Palucis","doi":"10.1029/2025JE009154","DOIUrl":"https://doi.org/10.1029/2025JE009154","url":null,"abstract":"<p>The discovery of rounded, river-transported clasts near an alluvial fan on Mars by the Mars Science Laboratory (MSL) <i>Curiosity</i> rover provided the first ground-based evidence for sustained fluvial activity on Mars. However, questions remain about the paleoenvironment and provenance of these clasts. Clast shape and size are two quantitative metrics that can be used to reconstruct the history of sedimentary deposits. However, previous studies have not considered martian basaltic lithologies and thus may not be applicable to Mars. We use data from two alluvial fans in the Mojave Desert and a series of rotating drum experiments to quantify the difference between fluvial abrasion of basalt and granite to determine the effect of lithology on the relationship between clast roundness and transportation distance. We find that basalt rounds faster than granite both in the field and in the lab and that this effect is more pronounced in the field. Our rounding coefficient (<i>p</i><sub><i>c</i></sub>) values show the agreement between granitic clasts in the lab (0.028 ± 0.0008) and the field (0.029 ± 0.003 and 0.034 ± 0.007) as well as the faster rounding rate we see from the field (0.048 ± 0.001) than experimental basalt (0.039 ± 0.001). This may be due to the greater susceptibility of basalt to chemical weathering. When applying our relationship to rounded clasts observed in the Gale crater, we find that our transport values agree with previous interpretations that point to a sediment provenance near the crater rim.</p>","PeriodicalId":16101,"journal":{"name":"Journal of Geophysical Research: Planets","volume":"130 11","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-11-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145522097","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kar Wai Cheng, Maxim D. Ballmer, Paul J. Tackley, Amir Khan
<p>A basal magma ocean (BMO) has been suggested to form as a consequence of cumulate overturn following crystallization of an initial global magma ocean. As a consequence of enrichment in iron and heat-producing elements (HPE), the BMO stabilizes gravitationally at the base of the planetary mantle, affecting the efficiency of mantle convection and crustal production. Seismic data collected during the recent Mars InSight mission provide compelling evidence for the presence of a molten silicate layer at the bottom of the Martian mantle. Here, we study the role of a BMO on the long-term evolution of Mars's interior using two-dimensional geodynamic models in spherical annulus geometry. We compare our model predictions with available observational constraints, including crustal HPE enrichment, seismically constrained present-day crustal and lithospheric thicknesses, and mantle temperatures. Among the parameters explored, we find that the interstitial porosity, which controls the initial HPE distribution between the BMO and the mantle, plays an important role on mantle thermal structure throughout Mars's evolution. Models with an intermediate interstitial porosity (<span></span><math>