Heat flow is an important constraint on planetary formation and evolution. It has been suggested that Martian obliquity cycles might cause periodic collapses in atmospheric pressure, leading to corresponding decreases in regolith thermal conductivity (which is controlled by gas in the pore spaces). Geothermal heat would then build up in the subsurface, potentially affecting present-day heat flow - and thus the measurements made by a heat-flow probe such as the InSight HP$^{3}$ instrument. To gauge the order of magnitude of this effect, we model the diffusion of a putative heat pulse caused by thermal conductivity changes with a simple numerical scheme and compare it to the heat-flow perturbations caused by other effects. We find that an atmospheric collapse to 300 Pa in the last 40 kyr would lead to a present-day heat flow that is up to $2-8\%$ larger than the average geothermal background. Considering the InSight mission with expected $5-15\%$ error bars on the HP$^{3}$ measurement, this perturbation would only be significant in the best-case scenario of full instrument deployment, completed measurement campaign, and a well-modelled surface configuration. The prospects for detecting long-term climate perturbations via spacecraft heat-flow experiments remain challenging.

We used infrared data from the Lunar Reconnaissance Orbiter (LRO) Diviner Lunar Radiometer Experiment to globally map thermophysical properties of the Moon's regolith fines layer. Thermal conductivity varies from 7.4$\times$10$^{-4}$ W m$^{-1}$ K$^{-1}$ at the surface, to 3.4$\times$10$^{-3}$ W m$^{-1}$ K$^{-1}$ at depths of $\sim$1 m, given density values of 1100 kg m$^{-3}$ at the surface, to 1800 kg m$^{-3}$ at 1-m depth. On average, the scale height of these profiles is $\sim$7 cm, corresponding to a thermal inertia of 55 $\pm$2 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$ at 273 K, relevant to the diurnally active near-surface layer, $\sim$4-7 cm. The temperature-dependence of thermal conductivity and heat capacity leads to a $\sim$2$\times$ diurnal variation in thermal inertia at the equator. On global scales, the regolith fines are remarkably uniform, implying rapid homogenization by impact gardening of this layer on timescales $<$ 1 Gyr. Regional and local scale variations show prominent impact features $<$ 1 Gyr old, including higher thermal inertia ($>$ 100 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) in the interiors and ejecta of Copernican-aged impact craters, and lower thermal inertia ($<$ 50 J m$^{-2}$ K$^{-1}$ s$^{-1/2}$) within the lunar cold spots identified by Bandfield et al. (2014). Observed trends in ejecta thermal inertia provide a potential tool for age-dating craters of previously unknown age, complementary to the approach suggested by Ghent et al. (2014). Several anomalous regions are identified in the global 128 pixels-per-degree maps presented here, including a high-thermal inertia deposit near the antipode of Tycho crater.

We assess local-time variations of epithermal-neutron count rates measured by the Lunar Prospector Neutron Spectrometer. We investigate the nature of these variations and find no evidence to support the idea that such variations are caused by diurnal variations of hydrogen concentration across the lunar surface. Rather we find an anticorrelation between instrumental temperature and epithermal-neutron count rate. We have also found that the measured counts are dependent on the temperatures of the top decimeters of the lunar subsurface as constrained by the Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment temperature measurements. Finally, we have made the first measurement of the effective leakage depth for epithermal-neutrons of ~20 cm.