Abstract Characterizing Martian surface temperatures and near-surface thermal behavior is essential for interpreting regolith thermophysical properties, surface–atmosphere interactions, and supporting mission operations. Traditional studies primarily utilize one-dimensional (1D) models, such as the widely adopted KRC framework, which assume strictly vertical heat transfer and laterally uniform terrain. While these models capture average diurnal trends, they neglect lateral heat transport, slope-dependent solar forcing, and realistic terrain-induced thermal variability. In this study, we present a comprehensive three-dimensional (3D) thermophysical model developed using the finite element method. By integrating high-resolution Digital Terrain Models (DTMs), the model explicitly accounts for surface slopes, roughness, and terrain-controlled solar illumination. We validated the model against KRC simulations under identical conditions and in-situ ground truth data from the InSight lander. To demonstrate its utility, the model was applied to a topographically complex sub-region within Jezero Crater at the Perseverance landing site. Our results show that while 3D simulations align with 1D mean diurnal trends, they resolve significant spatial temperature heterogeneities that 1D approaches fail to capture. This framework provides a more accurate representation of local thermal processes, offering critical applications for rover-scale data interpretation, landing site selection, future mission planning and science.
Raju et al. (Fri,) studied this question.