Lunar Trailblazer Spacecraft Tracking and Mission Recovery Attempt: Characterization of Status and Behavior of a Non‐Cooperative Object in Cis‐Lunar Space

Earth and Space Science American Geophysical Union (AGU) 13:8 (2026)

Authors:

BL Ehlmann, J Bellerose, G Lantoine, E Furlan, E Scire, S Fajardo‐Acosta, L Bennett, M Sanchez Net, TJW Lazio, M Brozović, J Masiero, AV Steckel, PA Burke, M Kimura, S Foxman, MPM Zaw, LM Lee, F Clarke, M Hauge, D McDonald, J Adler, R Strauss, D Trilling, CS Edwards, MC Nolan, D Lyster, L Robinson, AT Klesh, CC Seybold

Abstract:

Abstract Unexpectedly following launch, the Lunar Trailblazer mission experienced software anomalies that led it to orient solar panels away from the sun and lose communication with Earth. This paper describes efforts to determine the spacecraft state and attempt recovery of the mission's science at the Moon. First, ground observatories at optical and radar wavelengths were engaged to maintain custody of the spacecraft and knowledge of its trajectory. Second, viability of recovery of the mission science objectives was established via testbed work to understand system behavior in fault conditions and determination of trajectories and maneuvers that would enable lunar orbit insertion. Third, optical photometry and radar doppler broadening were employed to determine Lunar Trailblazer's spin and orientation, using approaches similar to those in asteroid studies, to establish when solar panels might again receive sufficient power to boot the spacecraft and initialize its radio. Fourth, X‐band‐capable groundstations in addition to the NASA Deep Space Network were engaged to monitor for the spacecraft's radio carrier signal nearly continually, including crowd‐sourced monitoring and tip‐and‐cue style commanding. Lunar Trailblazer left the Earth‐Moon system and is in a 14‐year Earth return, heliocentric orbit. As it moved further away from Earth prospects for recovery became formidable; ultimately, the ability of the telecom system to return telemetry to Earth would have been insufficient to enable actions to recover the spacecraft, and the recovery attempt ended 6 July 2025. Lunar Trailblazer's mission recovery efforts illuminate capabilities in characterizing a 1–3.5 m 3 size non‐cooperative object in cis‐lunar space. Plain Language Summary Lunar Trailblazer is a NASA small satellite mission to study the Moon. Software errors after launch caused the spacecraft to orient solar panels away from the sun and the radio to reinitialize in a state that did not allow commanding or telemetry receipt; having depleted its batteries, it is now in a low power, “browned out” state. This paper describes the Lunar Trailblazer team recovery efforts to reestablish Earth‐spacecraft communication and return to the Moon to achieve the science objectives. The team used ground observatories at optical and radar wavelengths to keep track of the trajectory of the spacecraft, determine its orientation, and compute trajectories for Earth‐Moon system return. Along with power model and testbed data on power use during the boot/initialization sequence, this was used to determine when the spacecraft solar panels could next receive power sufficient to initialize the radio. Ultimately, while enough propellant is in the tanks that could have enabled a return to Earth through calendar year 2025, the inability to receive telemetry due to distance ended recovery attempts. The procedures employed to characterize Lunar Trailblazer's state can be used on other cis‐lunar spacecraft and near‐Earth asteroids. Key Points Over a 130‐day mission recovery attempt, Lunar Trailblazer was a non‐cooperative cis‐lunar object, characterized using Earth ground assets Radar and optical observatory data enabled reconstruction of spacecraft orientation (spin rate, spin pole) and precise long‐term trajectory An ad‐hoc global X‐band ground station network that included crowd‐sourced signal monitoring and tip‐and‐cue style commanding was created

Direct Imaging Discovery of Giant Exoplanet β Pictoris d: A Decade-long Game of Hide-and-seek

The Astrophysical Journal Letters American Astronomical Society 1006:1 (2026) L10-L10

Authors:

Ben J Sutlieff, Markus J Bonse, Valentin Christiaens, Clémence Fontanive, Elisabeth C Matthews, Luke T Parker, Tim D Pearce, Jayne L Birkby, Beth A Biller, Trent J Dupuy, Emily O Garvin, Leyla Iskandarli, Jens Kammerer, Yifan Zhou, Robert J De Rosa, Aarynn L Carter, Sasha Hinkley, Matthew A Kenworthy, William O Balmer, Iain Hammond, James Mang, Caroline V Morley, Mark J Neeser, Olivier Absil, Anthony Boccaletti

Abstract:

We report the direct imaging discovery of a third exoplanet in the β Pictoris (β Pic) system. We detect β Pictoris d (β Pic d) in noncoronagraphic observations obtained with the Very Large Telescope (VLT) Enhanced Resolution Imager and Spectrograph (ERIS), as well as multi-epoch archival datasets from the JWST Near Infrared Camera (NIRCam) and VLT/SPHERE. Astrometric measurements over an 11 yr baseline demonstrate that it is consistent with a gravitationally bound source with orbital motion. Joint multi-planet orbit fits of all three planets in the system yield a semimajor axis of 26.0−6.1+2.2 au and inclination 89.0−0.6+0.7 deg for planet d. β Pic d has a larger orbital semimajor axis than the other known planets in the system, but is coplanar with the inner two planets, and its orbit is consistent with sculpting the inner edge of the debris disk. β Pic d has a contrast of ΔL′=12.11±0.15 mag, with colors and luminosity that closely match those of 51 Eri b, another exoplanet in the β Pic moving group. Its VLT/ERIS and JWST/NIRCam colors are distinct from those of free-floating planetary-mass objects of a similar age and temperature. Its red F410M − F444W color indicates strong CO2 absorption in its atmosphere and suggests significant enhancement in metals compared to free-floating objects. From the ATMO hot-start evolutionary models, we estimate an effective temperature of 600−60+45 K and mass of 2.4 ± 0.6 MJup, which also closely matches similar estimates for 51 Eri b. β Pic d is among the lowest-mass exoplanets imaged from the ground. This discovery highlights the deep sensitivity achievable with ground-based imaging in the mid-infrared and the discovery potential of future high-contrast observations with the Extremely Large Telescope.

Degradation of plume-deposited organics at Enceladus and implications for future surface missions

Copernicus Publications (2026)

Authors:

Tom Nordheim, Robert Grayson, Leonardo Regoli, Edith Fayolle, Morgan Cable, Shannon MacKenzie, Kevin Hand, Jasmina Wiemann, Carly Howett, Mathieu Choukroun, Christopher Paranicas

Abstract:

Saturn’s ocean moon Enceladus is among the most compelling targets for future astrobiology-focused missions. Material erupted from the moon’s south polar “Tiger Stripe” fractures is thought to originate from a global subsurface ocean, providing a natural pathway for ocean-derived compounds to be sampled without penetrating the ice shell. Cassini observations demonstrated that Enceladus’ plume contains water vapor, ice grains, salts, organics, phosphates, and other compounds of astrobiological interest. Consequently, future landed or plume-sampling missions to Enceladus will seek to determine whether ocean-sourced organic molecules, and potentially biosignatures, can be preserved and detected in plume material.However, once plume grains are erupted and deposited onto the surface, they are exposed to multiple space weathering agents. Solar ultraviolet radiation, Saturnian magnetospheric charged particles, and galactic cosmic rays can alter or destroy organic molecules contained within surface material. The degree of processing depends not only on radiation environment and molecular susceptibility, but also on the rate at which fresh plume fallout buries previously deposited material. Identifying locations and depths where organics may remain minimally processed is therefore key for mission planning and interpretation of future measurements.Here we model the degradation and preservation of plume-deposited organics on Enceladus, using the amino acids glycine and phenylalanine as representative organic molecules and potential biosignature proxies. We combine solar UV photodestruction calculations, Geant4-based charged-particle irradiation simulations, and estimates of local plume deposition rates to produce a coupled model of irradiation and burial across the surface. Solar UV photolysis rates are calculated using experimental rate coefficients adjusted for seasonal variation in photon flux at selected locations on Enceladus, including equatorial and southern-hemisphere sites at leading and trailing longitudes. Because even small abundances of non-water-ice contaminants can strongly influence UV attenuation, we also examine the effect of trace tholin-like material mixed into the ice. Charged-particle processing is estimated using Cassini-derived magnetospheric electron spectra and a numerical transport model for the galactic cosmic ray flux at Saturn.Our results show that preservation of plume-deposited organics is highly sensitive to burial rate and season. At low to moderate surface deposition rates, less than approximately 0.002 mm yr⁻¹, newly deposited surface material is strongly processed by solar UV radiation before burial can shield it from further exposure. In such regions, there may be no practical “safe sampling depth” at which pristine or minimally processed organics can be accessed, because material has already been altered before reaching depth. By contrast, in regions receiving the highest predicted plume deposition rates, especially during local winter when solar UV exposure is reduced, burial can occur rapidly enough to protect a significant fraction of the organic inventory. Under these conditions, the near-surface stratigraphy may consist of millimeter-scale layers that alternate between relatively well-preserved and heavily photoprocessed material, reflecting seasonal variations in irradiation during emplacement.Our findings have implications for the mission design and sampling strategy of future Enceladus missions. A landed mission seeking minimally processed ocean-derived organics should prioritize regions of high plume fallout, where burial is most effective, and should consider the seasonal context of deposition and sampling. The most favorable strategy may be to collect plume grains directly before surface deposition, e.g., from a funnel on the surface or plume fly-through, thereby avoiding prolonged radiolytic and photolytic alteration. If surface sampling is required, targeting the uppermost approximately 0.1 mm of material in high-deposition regions during local winter may maximize the probability of detecting less-processed organics. More broadly, our results demonstrate that preservation potential at Enceladus is governed by the competition between surface irradiation and plume-driven burial, and that this balance should be incorporated into future landing site selection, sampling-depth requirements, and biosignature-detection strategies.

Fine Layering Effects on Thermal Infrared Emissivity of CI Simulant Materials 

(2026)

Authors:

Emma-Catherine Belhadfa, Neil Bowles, Katherine Shirley

Abstract:

Introduction: Thermal infrared emissivity measurements of asteroid regolith analogs are challenging owing to atmospheric water vapor absorption, sample heating requirements, and the need for controlled atmospheric conditions [1], yet they provide fundamental constraints on surface thermal properties that cannot be obtained from reflectance spectroscopy alone [1]. While diffuse reflectance measurements have demonstrated that minimal fine dust coverage can dominate spectral signatures [2], spacecraft-based thermal emission instruments like the OSIRIS-REx Thermal Emission Spectrometer (OTES) observe different physical processes related to thermal emission rather than scattered light [3]. The disconnect between laboratory studies and spacecraft observations has thus limited our ability to interpret thermal infrared spectra of asteroid surfaces. Previous work using Space Resource Technology's CI simulant showed that 7-10 wt% fine dust coverage could impose fine-dominated reflectance features on coarse substrates [2], but the corresponding thermal emission properties remained uncharacterized. To bridge this gap, we conducted systematic thermal emissivity measurements of layered CI simulant materials using Oxford’s PASCALE instrument [4] under nitrogen atmosphere, constraining how dust deposition mechanisms affect the thermal emission processes observed by spacecraft instruments at airless bodies like asteroid (101955) Bennu. Methods: We measured thermal emission of layered CI simulant [5] samples using PASCALE under nitrogen atmosphere across 2000-400 cm⁻¹ (5-25 µm), eliminating atmospheric water vapor interference. Six layering configurations were tested, using 10 wt% fines (5% emissivity variations from unity), while the fluffy group shows more subdued but consistent spectral signatures. All method-dependent variations exceed the 2% measurement precision, demonstrating that dust deposition mechanism leaves diagnostic thermal emission signatures that can distinguish (and potentially identify) natural surface processes on airless body surfaces. Discussion: The separation between fluffy and compact layering methods demonstrates that thermal emission spectroscopy can distinguish surface formation processes on airless bodies. These results provide constraints missing from reflectance-only studies, by characterizing thermal emission properties relevant to spacecraft observations like OTES. The ability to spectrally distinguish between natural deposition processes offers new frameworks for understanding regolith evolution and thermophysical properties on asteroid surfaces. Summary: This study establishes thermal emissivity as a diagnostic tool for identifying dust deposition mechanisms on asteroid surfaces, demonstrating that layering processes leave distinct spectral signatures. References: [1] Salisbury et al. (1991) Icarus 92, 280-297. [2] Belhadfa et al. (2026) MaPs, In Prep. [3] Christensen P. R. et al. (2018) Space Science Reviews (Vol. 214, Issue 5). [4] Donaldson Hanna et al. (2019) Icarus 319, 701-723. [5] Landsman Z. et al. (2020) EPSC.  

Galileo PPR Thermal Inertia & Albedo Measurements of Europa, Ganymede and Callisto

Copernicus Publications (2026)

Authors:

Sarah Howes, Carly Howett, Duncan Lyster

Abstract:

IntroductionThe presence of endogenic hotspots provides a measure of the level of geologic activity of icy moons, since they are indicative of ongoing resurfacing processes. However, to avoid misinterpreting thermal abnormalities, it is first necessary to understand passive thermal emission that is governed by the physical structure of materials. Two important thermophysical properties in such analysis are bolometric Bond albedo and thermal inertia: if diurnal temperature variations can be accurately modeled by adjusting these two parameters, a passive rather than endogenic origin is possible. In this work, we aim to constrain the thermal inertia and Bond albedo across the surfaces of Europa, Ganymede, and Callisto using brightness temperature observations recorded by the Galileo Photopolarimeter-Radiometer (PPR) instrument [1]. By presenting estimated values and their uncertainties for these thermophysical properties, we prepare for future thermal measurements carried out by both Europa Clipper and JUICE.MethodsWe first take average diurnal temperatures of the surfaces of Europa, Ganymede, and Callisto using brightness temperatures recorded by PPR. In this analysis, we use 29 datasets for Europa, notably increased beyond previous efforts [2], comparable to [3].  With less observations available, only 7 datasets each were used for both Ganymede and Callisto. By translating the observed radiance into brightness temperatures, the variation in temperature with local time is determined for latitude and longitude bins across each of the three moons. These diurnal curves are compared to those predicted by a 1-D thermal model [4] to determine what thermal inertias and Bond albedos can fit the data within a reduced chi-squared cut-off of χ2red ≤1.0. This analysis is used to extensively quantify the uncertainty of the two thermophysical parameters derived from PPR data.ResultsEuropa: Ensuring closed upper and lower limits within our χ2red cut-off, we map albedo and thermal inertia for 33% and 24% of Europa's surface area into 6°x6° latitude/longitude bins (Fig. 1, left). We find a range of 0.375-0.75 for albedo and 20-110 J m-2 K-1 s-1/2 (MKS) for thermal inertia, agreeing with [2] and [3]. Our uncertainty analysis indicates well-constrained estimates for albedo, with the average higher and lower ranges overlapping within their uncertainties: Ahigh = 0.11 ± 0.06 and Alow = 0.18 ± 0.12. Uncertainties for thermal inertia remain poorly constrained, with average higher and lower ranges being ­Γhigh = 103(+153/-103) MKS and Γ­low = 17 ± 9 MKS.Ganymede: Due to less available surface coverage, Ganymede’s surface is divided into 18°x6° longitude/latitude bins in order to meet the diurnal fitting routine requirements (Fig. 1, right). Preliminary results indicate the Bond albedo remains nearly uniform, with an average of 0.42 ± 0.07 across the surface and agreeing well with previous work [6]. Specifically, our χ2red-deduced higher and lower uncertainty ranges of albedo are: Ahigh = 0.10 ± 0.06 and Alow = 0.14 ± 0.08. No apparent distinct surface variations in thermal inertia are as of yet discerned, with global values distributed across a range of 20-70 MKS. Upper limit thermal inertia estimates align with [6] within χ2red-deduced uncertainty ranges of ­Γhigh = 57 ± 35 MKS and Γ­low = 25 ± 11 MKS.Callisto: With less PPR coverage available, Callisto fits are performed across grouped hemispheric regions, where the Jovian and anti-Jovian hemispheres are analyzed separately. Each hemisphere is divided into 10° latitude strips. Preliminary results for thermal inertia and Bond albedo indicate an overall agreement with lower-bound estimations from previous literature [5,6]. Further constraints are expected to be obtained at the time of the conference.ConclusionThese results give an indication of the albedo and thermal inertia variation across Europa, Ganymede and Callisto. They aid in preparing for the arrival of Europa Clipper and JUICE to the Jupiter system by improving estimates for passive surface thermal properties and providing uncertainties of their values. This will enable future work to help discern what temperatures may lie above those expected from passive emission alone, providing a critical first step in the search for endogenic heating anomalies. Future work aims to refine the Ganymede and Callisto calculations using a two-component (ice and non-ice) analysis, and characterize the causes behind thermal inertia variations by modeling microphysical ice states for the three Galilean moons.Fig. 1: Thermal Inertia (top) and Bond albedo (bottom) map of Europa (left) and preliminary results for Ganymede (right).[1] Russell, E.E. et al., 1992. Space Sci Rev 60, 531-563.; [2] Rathbun, J.A. et al., 2010. Icarus 210, 763-769.; [3] Lange, L. et al., 2026. arXiv:2604.14374; [4] Lyster, D. et al., 2025. pp. EPSC-DPS2025-1479; [5] Meyer, C. et al., 2026. Planet. Sci. J. 7, 10; [6] Spencer, J.R. et al., 1989. Icarus 78(2), 337-354.