Solar system

Research group

Research theme

Sub department

About

Related research groups

Spectral Variability and Compositional Insights from Asteroid (101955) Bennu’s Sampling Sites Using OTES Data 

(2025)

Authors:

Emma Belhadfa, Katherine Shirley, Neil Bowles

Abstract:

Introduction: During the Reconnaissance phase of NASA’s OSIRIS-REx mission, the Thermal Emission Spectrometer (OTES) acquired high–spatial resolution emissivity spectra over Bennu’s four prospective sampling sites [1, 2]. We analyse the calibrated OTES dataset archived in the Planetary Data System [3] to quantify compositional and mineralogical diversity across the original four candidate sample sites (Nightingale, Kingfisher, Osprey, and Sandpiper) and to explore possible drivers of Bennu’s surface heterogeneity, including implications for Bennu’s mineralogy and space-weathering history.  Figure 1: Site-Averaged Emissivity Spectra with Annotated Band Parameters Methods: Calibrated emissivity spectra (5.7-100 µm) were linked to corresponding OCAMS imagery [5] to place the thermal infrared measurements in geological context, by cross-referencing observation times. For every spectrum we derived four diagnostic band parameters: Christiansen Feature (CF), silicate stretching band, silicate bending band and spectral slope, following the methods outlined in [6]. Each site contains thousands of spectral observations (site-averaged for visualization in Figure 1). The corresponding band parameters were compared using three statistical models: Principal Component Analysis (PCA) [5], k-Nearest Neighbors (KNN) [7], and Analysis of Variance (ANOVA) [8]. The three methods compare the mean and variance of each individual observation per site, considering how the in-group variance (i.e. the spread within all observations of a single site) compares to the out-group variance (i.e. the spread from other sites).  Results: Significant differences in emissivity spectra emerged among the four sites. PCA indicated that the first three components explain 85.5% of spectral variance, distinguishing Kingfisher as notably unique, with Sandpiper and Osprey exhibiting the greatest similarity. The KNN analysis corroborated PCA findings, reaching optimal classification accuracy (47%) at k = 21. ANOVA highlighted significant variability among the sites, especially in the spectral slope parameter (F = 762.8), suggesting differences in particle size distribution and space weathering could be driving factors in the detected heterogeneity [9]. Band ratio analyses provided additional insight into site-specific mineralogical distinctions, notably the relationship between silicate features and aqueous alteration indicators [10].  Figure 2: Distributions of Band Parameters by Site Discussion: Variability in spectral parameters aligns with documented particle size frequency distributions and known space weathering spectral types across Bennu’s surface [9]. Nightingale, the mission’s selected sample site, captures representative global characteristics, contrasting with Kingfisher’s distinct compositional and physical attributes, potentially related to differences in Fe/Mg content and degree of aqueous alteration [10].  Conclusion: Integrative use of multiple statistical approaches confirms the compositional and physical diversity of Bennu's surface, as seen through the four prospective sites. These analyses provide a framework for interpreting returned sample data and offer insights into the connections between mineralogy, particle size, and space weathering processes on small airless body surfaces.  References: [1] Lauretta D. S. et al (2021) Sample Return Missions. [2] Hamilton V. et al. (2021) A&A (Vol. 650). [3] Christensen, P. R. et al. (2019) NASA Planetary Data System [4] Christensen P. R. et al. (2018) Space Science Reviews (Vol. 214, Issue 5). [5] Rizk B. et al (2018) Space Science Reviews (Vol. 214, Issue 1). [6] Xie B. et al (2022) Minerals (Vol. 508, Issue 12). [7] Kramer O. (2013) Intelligent Systems Reference Library (13-23). [8] Sawyer S. (2009) Journal of Manual & Manipulative Therapy. [9] Clark B. E. et al (2023) Icarus (Vol. 400). [10] Bates H. et al (2020) MaPS (Vol. 55, Issue 1). 
More details from the publisher

Super-Earth lava planet from birth to observation: photochemistry, tidal heating, and volatile-rich formation

(2025)

Authors:

Harrison Nicholls, Tim Lichtenberg, Richard D Chatterjee, Claire Marie Guimond, Emma Postolec, Raymond T Pierrehumbert

Abstract:

Larger-than-Earth exoplanets are sculpted by strong stellar irradiation, but it is unknown whence they originate. Two propositions are that they formed with rocky interiors and hydrogen-rich envelopes (‘gas-dwarfs’), or with bulk compositions rich in water-ices (‘water-worlds’) . Multiple observations of super-Earth L 98-59 d have revealed its low bulk-density, consistent with substantial volatile content alongside a rocky/metallic interior, and recent JWST spectroscopy evidences a high mean molecular weight atmosphere. Its density and composition make it a waymarker for disentangling the processes which separate super-Earths and sub-Neptunes across geological timescales. We simulate the possible pathways for L 98-59 d from birth up to the present day using a comprehensive evolutionary modelling framework. Emerging from our calculations is a novel self-limiting mechanism between radiative cooling, tidal heating, and mantle rheology, which we term the 'radiation-tide-rheology feedback'. Coupled numerical modelling yields self-limiting tidal heating estimates that are up to two orders of magnitude lower than previous calculations, and yet are still large enough to enable the extension of primordial magma oceans to Gyr timescales. Our analysis indicates that the planet formed with a large amount (>1.8 mass%) of sulfur and hydrogen, and a chemically-reducing mantle; inconsistent with both the canonical gas-dwarf and water-world scenarios. A thick atmosphere and tidal heating sustain a permanent deep magma ocean, allowing the dissolution and retention of volatiles within its mantle. Transmission features can be explained by in-situ photochemical production of SO2 in a high-molecular weight H2-H2S background. These results subvert the emerging gas-dwarf vs. water-world dichotomy of small planet categorisation, inviting a more nuanced classification framework. We show that interactions between planetary interiors and atmospheres shape their observable characteristics over billions of years.
More details from the publisher

TEMPEST: A Modular Thermophysical Model for Airless Bodies with Support for Surface Roughness and Non-Periodic Heating

(2025)

Authors:

Duncan Lyster, Carly Howett, Joseph Penn

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.
More details from the publisher

Temperature, Composition, and Cloud structure in Atmosphere of Neptune from MIRI-MRS and NIRSpec-IFU Observations

(2025)

Authors:

Michael Roman, Leigh Fletcher, Heidi Hammel, Oliver King, Glenn Orton, Naomi Rowe-Gurney, Patrick Irwin, Julianne Moses, Imke de Pater, Henrik Melin, Jake Harkett, Simon Toogood, Stefanie Milam

Abstract:

We present observations and analysis of Neptune’s atmosphere from JWST, providing new constraints on hydrocarbon abundances, cloud properties, and temperature structure across the planet’s disk.  JWST observed Neptune in June 2023 (program1249) as part of the Solar System Guaranteed Time Observations (GTO). Integral field spectroscopy (IFS) with the Near-Infrared Spectrograph (NIRSpec) and the Mid-Infrared Instrument/Medium Resolution Spectrometer (MIRI/MRS) were combined to provide nearly simultaneous and continuous spatial and spectral data between 1.66 and 28.70 microns.We show how wavelengths sensitive to the atmospheric temperatures reveal a structure consistent with Voyager [1] and ground-based imaging [2,3], with a sharply defined warm polar vortex. In contrast, wavelengths sensitive to stratospheric hydrocarbons (namely acetylene and ethane) show a marked enhancement in the northern winter hemisphere.Finally, we examine the distribution and vertical structure of clouds in context of the temperature and chemical structure. Scattered light in NIRSpec observations indicate variable discrete clouds extend to pressures of roughly 50 mbar at the northernmost latitudes and south pole. [1] Conrath, B. J., F. M. Flasar, and P. J. Gierasch. "Thermal structure and dynamics of Neptune's atmosphere from Voyager measurements." Journal of Geophysical Research: Space Physics 96, no. S01 (1991): 18931-18939.[2] Fletcher, Leigh N., Imke de Pater, Glenn S. Orton, Heidi B. Hammel, Michael L. Sitko, and Patrick GJ Irwin. "Neptune at summer solstice: zonal mean temperatures from ground-based observations, 2003–2007." Icarus 231 (2014): 146-167.[3] Roman, Michael T., Leigh N. Fletcher, Glenn S. Orton, Thomas K. Greathouse, Julianne I. Moses, Naomi Rowe-Gurney, Patrick GJ Irwin et al. "Subseasonal variation in Neptune’s mid-infrared emission." The Planetary Science Journal 3, no. 4 (2022): 78.
More details from the publisher

Temperature, Composition, and Cloud structure in Atmosphere of Uranus from MIRI-MRS and NIRSpec-IFU Spectra

(2025)

Authors:

Michael Roman, Leigh Fletcher, Heidi Hammel, Patrick Irwin, Oliver King, Naomi Rowe-Gurney, Julianne Moses, Glenn Orton, Imke de Pater, Henrik Melin, Jake Harkett, Matthew Hedman, Simon Toogood, Stefanie Milam

Abstract:

Introduction: Due to Uranus’ weak thermal radiance, the thermal and compositional structures of its atmosphere have remained poorly characterised. Here, using the unprecedented sensitivity of JWST's MIRI and NIRSpec instruments, we present an analysis of Uranus' spatially resolved spectrum spanning the near- and mid-infrared, revealing how temperatures, composition, and clouds vary across the planet's northern hemisphere.Observations: JWST observed Uranus on 8--9 January 2023 (program1248) as part of the Solar System Guaranteed Time Observations (GTO). Integral field spectroscopy (IFS) with the Near-Infrared Spectrograph (NIRSpec) and the Mid-Infrared Instrument/Medium Resolution Spectrometer (MIRI/MRS) were combined to provide nearly simultaneous and continuous spatial and spectral data between 1.66 and 28.70 microns.Temperatures: The nearly continuous spectral coverage offered by the combination of NIRSpec and MIRI provide constraints on the temperature structure from the stratosphere down to several bars. The average temperature-pressure vertical profile is largely consistent with that determined from Spitzer [1], but the spatially resolved JWST reveal how these temperatures vary with latitude in the stratosphere and cloud layer for the first time [2]. They also suggest the possibility of a sub-adiabatic cloud layer.Chemistry: Our radiative transfer analysis of MIRI-MRS spectra 1) provide new constraints on minor species in Uranus’ stratosphere and 2) reveals how various hydrocarbons vary as a function of latitude. The observed distributions are indicative of a combination of seasonal photochemistry [3] and dynamical processes, as we will briefly discuss.Clouds and hazes: Finally, we briefly examine the vertical cloud structure and its latitudinal variation as sensed by NIRSpec data. The data reveal the opacity of Uranus clouds and hazes spanning the transition from scattered sunlight to thermal emission for the first time. The overall vertical structure suggested by these new data largely agrees with that of prior work [3,4,5], but a comparison between observed and model spectra reveal interesting discrepancies and possibly a need for additional sources of opacity. [1] Orton, G.S., Fletcher, L.N., Moses, J.I., Mainzer, A.K., Hines, D., Hammel, H.B., Martin-Torres, F.J., Burgdorf, M., Merlet, C., Line, M.R.: Mid-infrared spectroscopy of uranus from the spitzer infrared spectrometer: 1. determination of the mean temperature structure of the upper troposphere and stratosphere. Icarus 243, 494–513 (2014)[2] Roman, M.T., Fletcher, L.N., Orton, G.S., Rowe-Gurney, N., Irwin, P.G.: Uranus in northern midspring: persistent atmospheric temperatures and circulations inferred from thermal imaging. The Astronomical Journal 159(2), 45 (2020)[3] Moses, J.I., Fletcher, L.N., Greathouse, T.K., Orton, G.S., Hue, V.: Seasonal stratospheric photochemistry on uranus and neptune. Icarus 307, 124–145 (2018)[4] Sromovsky, L.A., Karkoschka, E., Fry, P.M., Pater, I., Hammel, H.B.: The methane distribution and polar brightening on uranus based on hst/stis, keck-nirc2, and irtf/spex observations through 2015. Icarus 317, 266–306 (2019)189[5] Irwin, P.G., Teanby, N.A., Fletcher, L.N., Toledo, D., Orton, G.S., Wong, M.H.,Roman, M.T., Perez-Hoyos, S., James, A., Dobinson, J.: Hazy blue worlds:A holistic aerosol model for uranus and neptune, including dark spots[6] Roman, M.T., Banfield, D., Gierasch, P.J.: Aerosols and methane in the ice giant atmospheres inferred from spatially resolved, near-infrared spectra: I. uranus, 2001–2007. Icarus 310, 54–76 (2018)
More details from the publisher