Dr Carly Howett

Abstract:

Introduction: Understanding surface temperatures on airless planetary bodies is crucial for interpreting thermal observations and constraining surface properties. We present TEMPEST (Thermal Evolution Model for Planetary Environment Surface Temperatures), a modular, open-source Python model that simulates diurnal and non-periodic thermal evolution on irregular bodies. Unlike traditional 1D periodic solvers, TEMPEST handles transient heating events such as eclipses, non-synchronous rotations such as tumbling asteroids, and seasonal variations. Key capabilities include surface roughness modelling via hemispherical craters, multiple thermal conduction schemes, and modular scattering using lookup tables (LUTs). TEMPEST has been used to analyse data from the Lucy mission and has been validated against the well-established Spencer 1D thermal model, thermprojrs [1].Figure 1: TEMPEST allows the user to select a facet to view any of its time varying properties including insolation, temperature and radiance. The diurnal temperature curves (right) are those of the corresponding outlined facets selected by the user in the interactive pane (left).Methods: TEMPEST calculates surface temperatures by solving a surface energy balance that includes solar flux, thermal emission, vertical heat conduction, and (optionally) radiative self-heating. Figure 1 shows the user interface once the model has completed a run. Key components include:Thermal solvers: Includes a standard 1D periodic conduction scheme influenced by the widely used thermprojrs [1] and a non-equilibrium solver, designed for better performance and stability in non-periodic cases. Scattering treatments: Utilises precomputed LUTs for various scattering laws (e.g., Lambertian, Lommel-Seeliger). This structure allows users to incorporate empirical bi-directional reflectance function (BRDF) data (e.g., from goniometer measurements of lunar regolith) or test the impact of different scattering assumptions, which can be particularly important for investigating the temperature of shadowed regions, as shown in Figure 2. The modularity also facilitates user modification for specific research needs. Surface roughness: Implemented via hemispherical sub-facet craters with adjustable rim angle to match roughness with a specified RMS slope angle. Non-periodic and time-dependent conditions: Supports time-dependent boundary conditions, including periodic scenarios such as eclipses and seasonal variations due to orbital eccentricity, as well as non-periodic cases including tumbling rotation, endogenic heating, and, or other user-defined transient heating scenarios. Designed for efficient parallel execution, the model runs effectively on multi-core personal computers and can efficiently simulate shape models with tens of thousands of facets. It has also been deployed on high-performance computing clusters for larger-scale models on the order of 1 million facets. Input configuration files are simple and flexible, allowing integration into larger analysis pipelines.Figure 2: An example insolation curve from a 1666 facet model of the bilobate comet 67P. The effects of scattered light can be seen either side of the main peak, this is particularly important for permanently shadowed regions. The selected facet is shown with a blue outline; sunlight direction is shown with a yellow arrow.Results: We validated TEMPEST by comparing temperature time series with Spencer’s 1D model thermprojrs [1] under idealised conditions, showing consistent results – see Figure 3. Applied to high-resolution shape models of 67P/Churyumov-Gerasimenko and 101955 Bennu, the model produces detailed temperature maps that reflect the significant influence of self-shadowing and local geometry, quantifying, for example, the temperature reduction in shadowed craters. Non-periodic simulations have been run to explore rotational transitions and eclipse effects, enabling new modes of comparison with observational datasets. The modular scattering and roughness components offer a powerful way to assess how sub-resolution scale parameters impact apparent thermal inertia and surface radiative behaviour. TEMPEST is already being used to interpret thermal data from recent missions, including Lucy, and can be adapted for upcoming datasets from targets like those of Comet Interceptor and Europa Clipper.Figure 3: TEMPEST shows good agreement with ‘industry standard’ thermophysical models in 1 dimension.TEMPEST is open-source and available at:github.com/duncanLyster/TEMPEST/Acknowledgement: This work was made possible by support from the UK Science and Technology Facilities Council. References:[1] Spencer, J.R., Lebofsky, L.A., and Sykes, M.V., 1989. Systematic biases in radiometric diameter determinations. Icarus, 78(2), pp.337-354.[2] Lyster, D., Howett, C., & Penn, J. (2024). Predicting surface temperatures on airless bodies: An open-source Python tool. EPSC Abstracts, 18, EPSC2024-1121.[3] Lyster, D.G., Howett, C.J.A., Spencer, J.R., Emery, J.P., Byron, B., et al. (2025). Thermal Modelling of the Flyby of Binary Main Belt Asteroid (152830) Dinkinesh by NASA’s Lucy Mission. Submitted to EPSC Abstracts, 2025.

Abstract:

Introduction: The Lucy mission's first asteroid flyby provided a unique and unexpected opportunity to study a binary asteroid system up close. Originally expected to encounter a single target, Dinkinesh, the discovery of its small, tidally locked moon, Selam, introduced additional opportunity and complexity to the interpreting flyby observations [1]. We present thermal modelling of the binary system, quantifying how the presence of Selam influenced radiance measurements and indicating its possible impact on thermal inertia estimates. Thermal inertia (TI) offers insight into surface properties such as grain size and regolith structure. Determining the TI of Dinkinesh adds to our understanding of small S-type asteroids and enables comparison within a binary, potentially revealing differences driven by tidal effects or surface evolution.Methods: We modelled the flyby geometry and instrument measurements using the new TESBY (Thermal Emissions Spectrometer flyBY) module of TEMPEST (the Thermophysical Equilibrium Model for Planetary Environment Surface Temperatures) [2] to simulate the thermal radiance of both bodies and assess their combined effect on interpretation of data from the Lucy Thermal Emission Spectrometer (L’TES) instrument [3].The Thermal Model: Dinkinesh and its satellite, Selam, were modelled in TEMPEST. A stereo-photogrammetric shape model is available for the primary target – Dinkinesh [4], with ~2 m lateral and ~0.5 m vertical resolution, covering ~60% of the surface. This shape model was downsampled to a dimensionally accurate model with 1266 facets with a resolution of ~35 m. A sphere of representative diameter (230 m [1]) was used for the satellite Selam.Figure 1: TESBY visualization of flyby. Global view of the flyby trajectory (left), and the FOV of the instrument (centre), with corresponding L’LORRI image for comparison [1] taken 0.54 minutes before closest approach (right). Input is the TEMPEST [5] result for the shape model of Dinkinesh, and representative diameter sphere for Selam. Parameters used: solar distance = 2.19 AU, rotation periods = 3.74 hours (Dinkinesh) and 52.7 hours (Selam) [1] thermal inertia (provisional) = 40 J m-2 s-1/2 K-1, geometric albedo = 0.27Flyby geometry: Building on the TEMPEST framework, the TESBY module is given the geometry information for the flyby and the thermal data from TEMPEST. Based on the 7.3 mrad Field-of-View (FOV) of the L’TES instrument [3] TESBY produces simulated radiance measurements by computing a weighted sum of blackbody curves from each visible facet, based on its temperature, projected area, and emission angle. Matching these modelled radiances to the instrument data allows us to fit for the thermal inertia of the asteroid. A complicating factor in this study is that the sensitivity of L’TES is not uniform across its FOV, including this effect in the model is the subject of ongoing work.Figure 2: Preliminary modelled radiance results (blue line) compared to L’TES observation (red) using the same model settings as Fig. 1. Scaled radiances (dotted line) are also provided (see main text for more information).Results: An example of the currently predicted model radiance is given by Figure 2. As it shows, there is a notable offset between the predicted and observed radiances. Accounting for the position of the targets in the L’TES FOV is expected to resolve the observed discrepancy in absolute radiance levels. However, as the scaled model shows, the predicted radiances are able to capture the shape of the L’TES radiance.We find that due to the slower rotation rate of Selam, the maximum surface temperatures on the satellite can be as much as 25 K higher than those on Dinkinesh (Fig. 1), meaning that despite the small size (lobe diameter of only 230 m, compared with 790 m for Dinkinesh [1]), the contribution to measured radiance is significant. This effect is highlighted by investigation of the integrated radiances of the targets throughout the flyby (Fig. 3), where the entry and exit of Selam within the FOV is visible, as well as the dip in integrated radiance while Selam is partially eclipsed by Dinkinesh. Our results demonstrate the importance of considering the full system in flyby analysis, informing techniques for similar encounters in the future. This work highlights how the thermal signature of even a small secondary body can significantly impact observations, shaping our understanding of asteroid surface properties and thermal environments.Continued analysis will focus on the use of TEMPEST/TESBY to constrain the thermal inertia of this binary asteroid from L’TES flyby observations.  Figure 3: Variation in integrated wavelength for Dinkinesh (target, blue), Selam (satellite, red) and combined effect (green). Radiances were integrated over the 200–1500 cm⁻¹ spectral range. The results show that despite its small size, Selam makes a significant difference to the spectral radiance, particularly at shorter wavelengths. The dip in combined spectral radiance at observations 3315-3320 is due to Selam being eclipsed by Dinkinesh.The thermal model code is open source and available at: github.com/duncanLyster/TEMPEST/Acknowledgement: This work was made possible by support from the UK Science and Technology Facilities Council.  References: [1] Levison, H.F., Marchi, S., Noll, K.S. et al. A contact binary satellite of the asteroid (152830) Dinkinesh. Nature 629, 1015–1020 (2024).[2] Lyster, D., Howett, C., & Penn, J. (2024). Predicting surface temperatures on airless bodies: An open-source Python tool. EPSC Abstracts, 18, EPSC2024-1121.[3] Christensen, P. R., et al. The Lucy Thermal Emission Spectrometer (L’TES) Instrument, Space Sci. Rev. (2023)[4] Preusker, F. et al. (2024). Shape Model of Asteroid (152830) Dinkinesh from Photogrammetric Analysis of Lucy’s Frame Camera L’LORRI. 55th Lunar and Planetary Science Conference, Abstract #1903.[5] Lyster, D., Howett, C., & Penn, J. (2025). TEMPEST: A Modular Thermophysical Model for Airless Bodies with Support for Surface Roughness and Non-Periodic Heating. Submitted to EPSC Abstracts, 2025

Abstract:

Introduction: Enceladus Thermal Mapper (ETM) is an Oxford-built high-heritage instrument that is being developed for outer solar system operations. ETM is based upon the design of Lunar Thermal Mapper (LTM, launched on Lunar Trailblazer, Fig. 1). It has a strong heritage story, including MIRMIS (on Comet Interceptor), Compact Modular Sounder (on TechDemoSat-1) and filters shared with Lunar Diviner (on Lunar Reconnaissance Orbiter). ETM is a miniaturized thermal infrared multispectral imager, with space for 15 spectral channels (bandpasses) that can be tailored to the mission requirements. It consists of a five-mirror telescope and optical system and an uncooled microbolometer detector array. Real-time calibration is achieved using a motorized mirror to point to an onboard blackbody target and empty space. ETM has an IFOV of 35 mm, so assuming a 100 30 km orbit it will have a spatial resolution of 40 to 70 m/pixel and a swath width of 14 - 27 km. ETM Updates: Through UK Space Agency funding we have developed three areas of ETM: its filter profile, radiation tolerance and sensitivity to Enceladus-like surfaces. Filters: ETM is a push broom thermal mapper, which works by the detector being swept over a surface. Each of the detector’s 15 channels is made up 16 rows, which are coadded to increase the signal to noise. A recently completed preliminary study has updated ETM’s bandpasses to include filters between 6.25 mm and 200 mm to enable it to detect Enceladus’ polar winter (170 K). Depending on the mission goals not all channels need to be utilised to achieve this, making some available for additional studies (e.g. searching for salt). Radiation: The radiation environments of Enceladus are vastly different to those of the Moon. Recent radiation testing and analysis showed that the majority of ETM’s existing design is already highly radiation tolerant. With some additional shielding and one component change all parts can reach the radiation hardness required to operate in the Saturn-system. The additional shielding may be provided by the spacecraft structure, depending on the adopted design. Sensitivity: ETM’s sensitivity to cryogenic surfaces is currently predicted through a well-characterised model. However, as part of the LTM calibration campaign we plan to directly measure its sensitivity to

Abstract:

NASA’s Cassini mission revealed endogenic activity at the south pole of Saturn’s moon Enceladus. The activity is concentrated along four fractures in Enceladus’ ice shell, which are much warmer than their surroundings and the source of Enceladus’ plumes. This work provides a review of the current state of knowledge of the energy and mass lost by Enceladus through this activity. Specifically, we discuss the composition of the plumes, along with their spatial and temporal variation. The mass flux loss predicted for the three plume constituents (gas, dust and charged particles) is reviewed and a total mass flux of ejected material that subsequently escapes Enceladus is estimated to be 2.1×1011 kg over a Saturn year. Given that Enceladus’ ocean is predicted to be 1019 kg this loss is sustainable in the very long term (∼1.5 billion Earth years). However, unless a resupply mechanism (such as serpentinization) exists molecular hydrogen is expected to be depleted within ∼1 million Earth years. The difficulty in determining Enceladus’ current heat flow is outlined, along with the advantages and disadvantages of the various techniques used to derive it. We find a robust lower limit for Enceladus’ exogenic production is 7.3 GW. Tidal heating models show endogenic emission of this level is sustainable, and Enceladus may have long-term near-surface heating (a result supported by studies of Enceladus’ geology). Finally, we offer suggestions for future observations, instrumentation, and missions. Enceladus remains a high-priority target for NASA, and as such it is highly likely that we will return to study this enigmatic world. Hopefully these missions will answer some of the questions that remain.