JWST near-infrared spectroscopy of the Lucy Jupiter Trojan flyby targets: Evidence for OH absorption, aliphatic organics, and CO2
Abstract
We present observations obtained with the Near Infrared Spectrograph on JWST of the five Jupiter Trojans that will be visited by the Lucy spacecraft — the Patroclus–Menoetius binary, Eurybates, Orus, Leucus, and Polymele. The measured 1.7–5.3 m reflectance spectra, which provide increased wavelength coverage, spatial resolution, and signal-to-noise ratio over previous ground-based spectroscopy, reveal several distinct absorption features. We detect a broad OH band centered at 3 m that is most prominent on the less-red objects Eurybates, Patroclus–Menoetius, and Polymele. An additional absorption feature at 3.3–3.6 m, indicative of aliphatic organics, is systematically deeper on the red objects Orus and Leucus. The collisional fragment Eurybates is unique in displaying an absorption band at 4.25 m that we attribute to bound or trapped CO2. Comparisons with other solar system small bodies reveal broad similarities in the 2.7–3.6 m bands with analogous features on Centaurs, Kuiper belt objects (KBOs), and the active asteroid 238P. In the context of recent solar system evolution models, which posit that the Trojans initially formed in the outer solar system, the significant attenuation of the 2.7–3.6 m absorption features on Trojans relative to KBOs may be the result of secondary thermal processing of the Trojans’ surfaces at the higher temperatures of the Jupiter region. The CO2 band manifested on the surface of Eurybates suggests that CO2 may be a major constituent in the bulk composition of Trojans, but resides in the subsurface or deeper interior and is largely obscured by refractory material that formed from the thermophysical processes that were activated during their inward migration.
1 Introduction
Target | (km)aaDiameter and visible geometric albedo measurements are taken from the literature: Mueller et al. (2010) for Patroclus, Buie et al. (2018) for Polymele, Buie et al. (2021) for Leucus, and Mottola et al. (2023) for Eurybates and Orus. | aaDiameter and visible geometric albedo measurements are taken from the literature: Mueller et al. (2010) for Patroclus, Buie et al. (2018) for Polymele, Buie et al. (2021) for Leucus, and Mottola et al. (2023) for Eurybates and Orus. | Swarm | Color |
---|---|---|---|---|
617 PatroclusbbIn this paper, Patroclus refers to the combined Patroclus–Menoetius binary. The listed diameter is for the effective projected-area-equivalent sphere used for thermal modeling. | 146 | 0.045 | L5 | less-red |
3548 EurybatesccEurybates is a collisional fragment with a significantly more neutral visible color than the other less-red objects. | 69 | 0.044 | L4 | less-red |
21900 Orus | 60 | 0.040 | L4 | red |
11351 Leucus | 44 | 0.037 | L4 | red |
15094 Polymele | 22 | 0.074 | L4 | less-red |
Notes.
Target | UT Date | (au)aaHeliocentric distance , distance from JWST , phase angle , and -band apparent magnitude of the target at the time of the JWST observation. | (au)aaHeliocentric distance , distance from JWST , phase angle , and -band apparent magnitude of the target at the time of the JWST observation. | (deg)aaHeliocentric distance , distance from JWST , phase angle , and -band apparent magnitude of the target at the time of the JWST observation. | (mag)aaHeliocentric distance , distance from JWST , phase angle , and -band apparent magnitude of the target at the time of the JWST observation. | G235M (s)bbTotal exposure times across both dithers in the two grating settings. | G395M (s)bbTotal exposure times across both dithers in the two grating settings. |
---|---|---|---|---|---|---|---|
617 Patroclus | 2022 September 7 | 4.86 | 4.51 | 11.7 | 15.7 | 146 | 350 |
3548 Eurybates | 2022 November 9 | 4.76 | 4.14 | 10.1 | 17.0 | 408 | 1342 |
21900 Orus | 2022 September 5 | 4.94 | 4.84 | 11.9 | 17.7 | 379 | 1167 |
11351 Leucus | 2022 November 9 | 5.02 | 4.37 | 9.3 | 18.2 | 700 | 2889 |
15094 Polymele | 2022 November 21 | 4.80 | 4.08 | 9.1 | 18.8 | 2101 | 3326 |
Notes.
In 2021 October, NASA’s Lucy spacecraft set out on its decade-long mission to perform the first ever flyby exploration of the Jupiter Trojans. The Trojans are confined to two vast swarms along Jupiter’s orbit, centered on the L4 and L5 Lagrangian points that lead and trail the gas giant by roughly one-sixth of its orbit. These enigmatic small bodies are situated on the edge between the inner and outer solar system and hold important clues to the evolution of the entire planetary system. In the classical picture of solar system evolution, the Trojans formed either in situ or on the outer edge of the nearby main asteroid belt and should therefore bear a genetic connection to the inner solar system (Marzari & Scholl, 1998; Marzari et al., 2003). In contrast, more recent dynamical instability models posit a drastically different scenario in which the Trojans originated in the primordial disk of icy planetesimals located beyond the initial orbits of the ice giants and were later scattered inward and captured into resonance during an epoch of chaotic restructuring of the giant planets’ orbits (e.g., Gomes et al., 2005; Morbidelli et al., 2005; Tsiganis et al., 2005; Levison et al., 2008). Within this dynamical framework, the Trojans formed alongside the outer solar system small bodies, including the present day Kuiper belt objects (KBOs). Direct comparisons of the Trojans’ observable properties with those of other small body populations promise to place the formation and evolution of Trojans definitively within our broader understanding of the early dynamical history of the solar system.
Past ground- and space-based observations showed that the reflectance spectra of Jupiter Trojans appear featureless at wavelengths below 2.5 m (e.g., Dotto et al., 2006; Fornasier et al., 2007; Emery et al., 2011; Sharkey et al., 2019; Wong et al., 2019a), while spectroscopy at longer wavelengths revealed absorption around 3.0 m that may be indicative of fine-grained water ice or an N–H stretch feature, as well as weaker spectral features in the 3.3–3.6 m region that are usually associated with organics (Brown, 2016). Photometric surveys uncovered a color bimodality that divides the Trojans into the so-called “less-red” and “red” subpopulations. The typical spectra of these two subpopulations have different spectral slopes at visible wavelengths that roughly correspond to the P- and D-type asteroids, respectively (Szabó et al., 2007; Roig et al., 2008; Wong et al., 2014; Wong & Brown, 2015). Lucy will rendezvous with five Trojan targets that sample the full diversity of surface types (Levison et al., 2021). The 617 Patroclus–Menoetius binary and 15094 Polymele have colors consistent with the less-red group, while 21900 Orus and 11351 Leucus represent the red subpopulation. 3548 Eurybates is a member of a collisional family (e.g., Brož & Rozehnal, 2011) and has a significantly more neutral visible color that is similar to that of C-type asteroids (Fornasier et al., 2007). Table 1 lists the targets and their properties, in order of decreasing size.
Lucy is equipped with a low spectral resolution infrared imaging spectrometer that will obtain full-disk spectra of each Trojan target from 0.9 to 3.8 m at a spectral resolving power of roughly 250 (Olkin et al., 2021). At close approach, Lucy will carry out spatially resolved scans of portions of the surface at a spatial resolution of a few kilometers, as well as multicolor imaging at higher spatial resolution. These datasets will provide a unique opportunity to explore the detailed chemical composition and surface geology of these Trojans at fine spatial scales. In preparation for the encounters in 2027–2033, a multifaceted campaign of ground- and space-based observations is being undertaken to maximize our understanding of the Trojan targets prior to the flybys and compile complementary datasets that will be indispensable for the interpretation of the Lucy observations.
As part of this effort, near-infrared spectroscopic observations of the five Lucy Trojan targets were obtained at a spectral resolving power 1000 with the Near-Infrared Spectrograph (NIRSpec) instrument (Jakobsen et al., 2022; Böker et al., 2023) on JWST. These high-quality 1.7–5.3 m spectra are ideally suited to probe the wavelength regions where previous detections of spectral features have been reported, while also significantly expanding the wavelength coverage beyond 4 m into a region that is out of reach of both ground-based observatories and Lucy’s instruments. In this manner, the JWST observations will help paint a more complete picture of the Lucy flyby targets. By placing these JWST spectra within the geological context that Lucy will provide, we will obtain a powerful tool for understanding the surface composition, geology, and evolutionary history of these objects. More broadly, the insights gained from a joint analysis of the JWST and Lucy flyby data will permit significantly more sophisticated interpretations of the myriad future JWST observations of asteroids, irregular satellites, and KBOs, most of which will never benefit from spacecraft visits.
2 JWST Observations
Near-infrared observations of the five Trojan targets were obtained using the NIRSpec integral field unit (IFU) as part of JWST Cycle 1 General Observers Program #2574 (PI: M. Brown). All five visits utilized the same observing strategy: after slewing and instrument setup, a pair of dithered exposures was collected with the G395M grating, followed immediately by a pair of shorter-duration exposures with the G235M grating. These two gratings provide continuous, overlapping wavelength coverage from 1.7 to 5.3 m at a spectral resolving power of around 1000. The observations employed the IRS2 readout method, which reduces the level of correlated noise caused by electronic drifts during detector readout (so-called noise; Moseley et al. 2010; Rauscher et al. 2012). Details of the observations are summarized in Table 2.
The uncalibrated data products were processed through a custom spectral extraction pipeline. The conversion from the raw 2D detector readouts to fully calibrated 3D spectral data cubes was handled by v1.11.3 of the official JWST calibration pipeline jwst (Bushouse et al., 2022). The calibration reference files required for bias subtraction, flat-fielding, wavelength calibration, flux calibration, and other relevant data processing routines were automatically downloaded from context jwst_1100.pmap of the JWST Calibration Reference Data System. Several adjustments were made to the default data processing workflow. After the stacks of nondestructive detector readouts (uncal files) were converted to 2D countrate images (rate files), we manually corrected for residual low-level read noise by subtracting the 200 pixel wide moving median of the off-sky pixels from each detector column. This process removed the visible column striping that typically persists on NIRSpec full detector images even when the IRS2 readout method is selected. After the spectral data cubes were constructed for each dithered exposure, we did not proceed with the final stage of the jwst pipeline, which aligns the moving target between dither positions and generates a single combined data cube for each NIRSpec grating setting. Instead, we extracted the spectrum from each dithered exposure separately, which allowed us to inspect the individual spectra to examine minor spectral features and distinguish them from correlated noise or other low-level instrumental artifacts.


The output IFU data cubes (s3d files) contain a stack of 2D wavelength slices with a field of view and a spatial pixel scale of . Figure 1 shows the median-averaged slices of Patroclus–Menoetius and Eurybates from the first dithered exposure with the G235M grating; the slices have been resampled to the sky-projected coordinate frame by the JWST pipeline. The point-spread function (PSF) of the Patroclus–Menoetius binary has a discernibly extended shape. According to the binary orbit ephemerides in Grundy et al. (2018), Menoetius was situated roughly (i.e., 2.5 pixels) to the west of Patroclus at the time of observation (toward the right in Figure 1). Meanwhile, the additional source apparent on the image roughly 7 pixels to the right of the binary is a background object. We experimented with using the webbpsf package (Perrin et al., 2014) to simulate the combined PSF of the binary, but found that the generated PSFs were not an adequate match to the observed PSFs, resulting in significant residuals in the modeling and unreliable flux extraction for the two binary components. Instead, we carried out spectral extraction on the blended binary as a single source. Previous differential photometry of Patroclus and Menoetius during mutual events did not reveal any significant systematic differences in color between the binary components (Wong & Brown, 2019), suggesting large-scale similarity in surface composition. Hereafter, we refer to the binary as Patroclus for the sake of brevity. The remaining four targets show radially symmetric PSFs that are consistent with point sources. In the case of Orus, the target was blended with a background star during one of the two G235M exposures, and the corresponding data cube was not included in our analysis.
Spectral extraction was carried out using an empirical PSF fitting technique that has been applied to other NIRSpec IFU observations of minor bodies (see, for example, Brown & Fraser, 2023; Emery et al., 2024; Grundy et al., 2024, for more detailed descriptions of the methodology). The source PSF in each wavelength slice was modeled from the median frame that was constructed using a sliding window spanning . The width of this window was chosen to be sufficiently wide to average out interpixel jitter and uncorrected bad pixels, while simultaneously being narrow enough to reflect the local shape of the PSF accurately, which exhibits measurable shape oscillations across the wavelength range due to undersampling of the PSF at the NIRSpec IFU pixel scale, drifts in the centroid of the dispersed spectral trace across detector rows, and corresponding resampling artifacts in the cube building process. The background region was defined as all pixels outside of a pixel box centered on the target centroid position in the IFU field of view, and the median background flux was subtracted from the median frame to generate the PSF template. This template was set to zero in the background region and renormalized to a sum of unity, before being fit to the wavelength slice with a flux scaling factor and an additive background level.
We optimized the size of the extraction region within which the PSF template was defined to produce the cleanest spectrum. We found that including the extended diffraction wings of the source PSF (see Figure 1) contributes additional scatter to the resultant spectra, so we selected extraction regions that primarily sampled the peak of the source PSF, where the signal-to-noise ratio is the highest and the PSF templates can be reliably modeled. Pixels outside of the extraction region but inside the background region were masked during the fitting process. For Patroclus, we used a square extraction region of side length 11 pixels in order to encompass the PSFs of both binary companions. For Eurybates, Orus, and Leucus, a pixel box was chosen, while for the faintest target Polymele, a smaller pixel yielded reduced correlated noise in the extracted spectrum when compared to larger apertures. The pairs of dithered exposure spectra were passed through a 20 pixel wide moving median filter to remove outliers before being averaged together to produce the final combined irradiance spectra. The irradiance spectra of the five Trojan targets are plotted in the left panel of Figure 2. At short wavelengths, the irradiance closely mirrors the shape of the incident sunlight spectrum, while at m, the thermal emission from the targets’ surfaces dominates. At the spectral resolution of the G235M and G395M gratings, solar lines are also visible throughout the wavelength range.
To derive the reflectance spectra, we divided the targets’ irradiance spectra by the spectrum of the G2V-type solar analog star P330E, which was observed as part of Program #1538 (PI: K. Gordon). Each target spectrum was paired with a stellar spectrum that was produced with the same extraction aperture size. This process self-consistently corrected for both common-mode instrumental systematics present in the target and stellar spectra and wavelength-dependent flux biases due to the monotonically increasing size of the source PSF within our fixed extraction apertures. There are strong correlated noise artifacts at the long-wavelength end of the G235M spectra, and we removed the portions between 3.05 and 3.15 m prior to analyzing and plotting the reflectance spectra.
The spectra of all targets show significant contributions from thermal emission at wavelengths longer than 4 m. We modeled and removed this thermal contribution using the Near Earth Asteroid Thermal Model (NEATM; Harris 1998), following the methods described in McGraw et al. (2022). When implementing NEATM, we fixed the targets’ sizes and albedos to literature values (see Table 2). We assumed a linear reflectance continuum, nominally pinned to the spectrum at 2.6 m, and varied the beaming parameter in the model to determine the level of thermal emission that best fits the measured spectrum at longer wavelengths. This model thermal contribution was then subtracted to produce the final reflectance spectra. We also varied the slope and vertical position of the assumed reflected continuum to evaluate the range of reasonable spectral shapes. In all cases, fixing the continuum slope to that of a line between 2.6 and 3.7 m (wavelengths that are dominated by reflected flux) minimized the resultant residuals in the modeled spectra. After removal of the thermal component, some artifacts became apparent at the longest wavelengths, manifesting as unphysical breaks in the continuum shape. These artifacts were likely caused by instrumental systematics that were not corrected by dividing the solar analog spectrum, and we trimmed the spectra beyond 5.0 m.
3.0 m OH Band | 3.4 m Organic Band | 4.25 m CO2 Band | |||||
---|---|---|---|---|---|---|---|
Target | Center (m) | Depth (%) | 2.9 m Depth (%) | Center (m) | Depth (%) | Center (m) | Depth (%) |
617 Patroclus | |||||||
3548 Eurybates | |||||||
21900 Orus | |||||||
11351 Leucus | |||||||
15094 Polymele |
The wavelength grid of the IFU data cubes generated by the JWST calibration pipeline is oversampled by a factor 2.2 relative to the true spectral resolution of the G235M and G395M gratings.111Tabulated dispersion and resolution curves for all NIRSpec gratings can be found at https://jwst-docs.stsci.edu/jwst-near-infrared-spectrograph/nirspec-instrumentation/nirspec-dispersers-and-filters. We used median binning to downsample the spectra by a factor of 3. This binning factor is used for all figures displaying the JWST Trojan reflectance spectra. The final binned and thermal-corrected reflectance spectra are plotted in the right panel of Figure 2; here, the targets are arranged in order of increasing visible spectral slope (i.e., redness). The flux uncertainties computed in our PSF template fitting procedure are comparable to the inherent scatter in the resultant spectra. The spectra obtained with the G395M grating were scaled to match the reflectance level in the G235M spectra at 3.0 m. For most targets, the required scaling factor was small (), while for Orus a larger 2% downward adjustment of the G395M spectrum was needed, which can be attributed to low-level systematic flux offsets between the two dither positions and the aforementioned removal of one of the two G235M dithered exposures due to background star blending.
3 Surface Characterization
The JWST/NIRSpec observations represent a substantial leap in our knowledge of Trojan spectra. Previous observations of these five Trojan targets have largely been confined to wavelengths below 2.5 m. Sharkey et al. (2019) presented 0.7–2.5 m reflectance spectra obtained with the NASA Infrared Telescope Facility (IRTF). In the overlapping wavelength range (1.7–2.5 m), most of their spectra have continuum slopes that agree with our JWST spectra. For Polymele, the IRTF spectrum shows a somewhat steeper median slope than what we obtained from JWST; however, the former spectrum has a very low signal-to-noise ratio within the wavelength range of interest, allowing for a wide range of spectral slopes, including neutral values that are fully consistent with the continuum shape shown in Figure 2.
Our spectra offer improved precision relative to previously published results and expand our purview beyond 4 m into a hitherto unexplored region of wavelength space. The reflectance spectra display several distinct spectral features across the 2.5–4.5 m wavelength range. In this section, we characterize these features through comparisons with published results of other minor bodies, with the goal of identifying plausible chemical species on the surfaces and placing the Trojans within the broader context of inner and outer solar system small body populations.
3.1 3 m OH Band
The most prominent feature in the Trojan spectra is the broad absorption band that is centered around 3 m. This band is deepest on Eurybates and clearly discernible on the other less-red objects Patroclus and Polymele. Meanwhile, on the red objects Orus and Leucus, the presence of this band is considerably more marginal. A characteristic property of the absorption is its distinct short-wavelength edge at roughly 2.7 m. Previously, Brown (2016) detected a portion of this 3 m band in Keck/NIRSPEC spectra of 16 Trojans, with less-red objects displaying deeper band depths than red objects. The significantly higher signal-to-noise ratio JWST/NIRSpec spectra confirm the earlier detections and robustly corroborate the trend between surface color and 3 m band depth.
We determined the center of this broad absorption by dividing each spectrum by a linear continuum fit across the feature, fitting a third-order polynomial to the bottom of the band, and taking the location of the minimum of the best-fit polynomial as the band center. For Orus and Leucus, we excluded the region around the 3.4 m organic band from the polynomial fit (see Section 3.2). The band depth was calculated by averaging the continuum-divided spectrum within a 0.05 m wide window centered around the band center. This window width was selected so that 25–30 channels are averaged. This process was repeated 10,000 times, each time resampling the data points using normal distributions with widths equal to the respective uncertainties. The reported values and uncertainties are the mean and standard deviation of the resulting parameter distributions. We also computed the band depth at 2.9 m to enable more direct comparisons with ground-based spectra of near-Earth and main belt asteroids, whose band depths are often reported at this wavelength in the literature. Table 3 shows the results of our band analysis. The continuum-subtracted Trojan spectra in this region are presented in Figure 3.

The continuous wavelength coverage of the JWST spectra enables us to perform a more detailed characterization of the 3 m band than was previously possible. High-quality spectra of small bodies covering the 3 m region are scant in the literature. Ground-based observations in this wavelength range are challenging due to significant telluric absorption, which almost entirely suppresses the incoming light at 2.5–2.9 m. While a significant body of IRTF/SpeX spectra is available for main belt asteroids (Takir & Emery, 2012), including many that show superficially similar 3 m absorption features to the Trojans, the gap in wavelength coverage precludes definitive compositional inferences. The only available continuous reflectance spectra in this region were obtained from space telescopes or spacecraft flyby missions. In Figure 3, we plot a collection of representative spectra that show various types of absorption features that have been observed on main belt asteroids.
The spectrum of Pallas, derived from measurements made by the AKARI spacecraft (Usui et al., 2019), shows a distinctive sharp 2.7 m absorption that is diagnostic of hydrated silicates — specifically, the O–H stretch feature in phyllosilicates (e.g., Lebofsky, 1980; Jones et al., 1990; Rivkin et al., 2002). The spectra of Ceres (from globally averaged observations by the Dawn spacecraft; De Sanctis et al. 2018) and Themis (compiled from AKARI and ground-based IRTF spectra; Rivkin & Emery 2010; Usui et al. 2019) display the same 2.7 m hydrated silicate band. For Trojans generally, we can exclude hydrated silicates as a major component in the surface regolith due to the absence of a sharp absorption band centered at 2.7–2.8 m as seen in Ceres, Palas, and Themis.
Themis and Ceres show a second spectral feature around 3.0–3.2 m; at even longer wavelengths, there are other features from organics and/or carbonates, which will be discussed in the following subsection. Several interpretations have been forwarded for the 3.0–3.2 m band. In the case of Themis models incorporating very fine water frost ( m) coating underlying refractory grains provide a good match (Campins et al. 2010; Rivkin & Emery 2010; but see McKay et al. 2017; O’Rourke et al. 2020 for arguments against this interpretation). Meanwhile, N–H stretch absorption from ammoniated surface species has been hypothesized to be a primary contributor to the feature on Ceres (e.g., King et al., 1992; De Sanctis et al., 2015). The lack of absorption around 3.0–3.2 m in the Trojan spectra leads us to argue against the presence of fine-grained water frost or ammoniated species as a distinct surface element on Trojans.
We have also included the spectra of two active objects: the Rosetta spectrum of the Jupiter-family comet 67P/Churyumov–Gerasimenko (67P; Poch et al., 2020) and the recently published JWST/NIRSpec spectrum of the active outer main belt asteroid 238P/Read (238P; Kelley et al., 2023). 67P lacks the hydrated silicate feature at 2.7 m while having a complex of absorption bands spanning 2.8–3.6 m that is attributed to a combination of ammoniated salts and organics. In contrast, the spectrum of 238P (shown separately in the bottom panel of Figure 3) has a single broad 2.7–3.6 m band that is quite comparable to the feature seen in the less-red Trojan spectra, albeit with the reflectance minimum situated at somewhat longer wavelengths. Although 238P is much smaller (0.6 km) than the Trojans there were observed by JWST, the similarity in spectral shapes is noteworthy. 238P is more dynamically akin to typical outer main belt asteroids than to the Jupiter-family comets, indicating that it was not recently implanted into the inner solar system from the Kuiper belt (Haghighipour, 2009; Hsieh et al., 2016). Nevertheless, it is still possible that the origin of 238P lies in the outer Solar System, as dynamical instability models predict a nonnegligible contingent of primordial trans-Neptunian planetesimals among the outer main belt asteroid population (e.g., Vokrouhlický et al., 2016). Future spectroscopic surveys of more small outer main belt asteroids may uncover other objects with Trojan-like 2.7–3.6 m absorption bands.
Given the predictions from dynamical instability models that place the origin of the Trojans in the outer solar system, we now turn to KBOs for further insight into the 3 m feature. Prior to JWST, only a handful of KBOs had been spectroscopically studied at wavelengths longer than 3 m. During the first year of science operations, JWST observed more than 70 KBOs, providing a treasure trove of 1–5 m reflectance spectra spanning all dynamical classes and a wide range of sizes. For comparison with our Trojan spectra, we analyzed publicly available NIRSpec observations of a pair of KBOs that were obtained as part of General Observers Program #2418 (PI: N. Pinilla-Alonso; Pinilla-Alonso et al. 2024). Photometric surveys of KBOs have shown that objects smaller than 1000 km in diameter can be broadly divided into two groups based on their visible colors (e.g., Fraser & Brown, 2012; Wong & Brown, 2017; Fraser et al., 2023). The two KBO targets we selected—2004 NT33 (420 km; Vilenius et al. 2014) and 2004 TY364 (510 km; Lellouch et al. 2013)—are respectively from the so-called red and very-red subgroups within the hot classical population.
In addition, we included two Centaurs (i.e., inward-scattered KBOs that cross the giant planet region on unstable orbits) that were observed as part of Program #2418 and published in Licandro et al. (2024): Okyrhoe and 2003 WL7 (35 and 103 km, respectively; Duffard et al. 2014). These Centaurs have perihelion and aphelion distances of au and au, and au and au, respectively. All four observations were obtained with the low spectral resolution prism grating (0.6–5.3 m; 30–300), and we processed the observations following the same methodology that was outlined in Section 2. The prism-mode spectrum of the G-type solar analog star SNAP-2 (Program #1128; PI: N. Luetzgendorf; publicly available data) was used to produce the reflectance spectra, which we plot in Figure 4.

The spectra of 2004 NT33 and 2004 TY364 differ broadly in several ways. 2004 NT33 shows strong water ice absorption bands at 1.5, 2, and 3 m, including a distinct feature at 1.65 m and the Fresnel reflection peak at 3.1 m, both of which are indicative of the crystalline form of water ice (see Mastrapa et al. 2013 for a review of amorphous and crystalline water ice on solar system objects). Meanwhile, the spectrum of 2004 TY364 is considerably less water ice rich, with only marginal detections of the 1.5 and 2 m absorption features and no signs of the Fresnel peak. It shows a deep and broad 2.7–3.3 m absorption feature and an additional rounded spectral band at 3.3–3.6 m that corresponds to aliphatic organics (see Section 3.2). The pair of narrow features between 2.6 and 2.8 m are overtones of the band of CO2 ice (e.g., Quirico & Schmitt, 1997; Baratta & Palumbo, 1998), which have been detected in the JWST spectra of many KBOs (de Prá et al., 2024; Pinilla-Alonso et al., 2024). Moving inward to the Centaurs, we find that the spectrum of 2003 WL7 is largely consistent with 2004 NT33, showing clear water ice absorption bands, albeit with somewhat reduced depths. In contrast, the water ice features on Okyrhoe are strongly attenuated, and the shallow 3 m band has a more asymmetric profile that gently rises toward the continuum at 3.6 m.
The four KBO and Centaur spectra shown in Figure 4 are broadly representative of the range of 3 m band shapes seen on similarly sized bodies throughout the trans-Neptunian region (Licandro et al., 2024; Pinilla-Alonso et al., 2024). Placing these spectra alongside our results for the Trojans, we note the strong similarity in the general shape and width of the broad 3 m absorption features to the observed feature on Eurybates, Patroclus, and Polymele, with consistent short-wavelength edges located at roughly 2.7 m. This comparison suggests that the same molecular vibration — namely, the fundamental O–H stretch mode (e.g., Mastrapa et al., 2009) — is responsible for the 3 m absorption across all of these objects. In the case of KBOs and Centaurs, the primary contributor to the 3 m band is water ice (both amorphous and crystalline), which is known to be widely distributed across the populations based on detections of the 1.5 and 2 m bands (e.g., Barkume et al., 2008; Barucci et al., 2011).
The identity of the OH-bearing component on the Trojans’ surfaces is less clear. The absence of the 2 m water ice band precludes a significant water ice fraction across the surface, consistent with previous observations (Emery & Brown, 2003; Yang & Jewitt, 2007). This finding also aligns with the results from thermophysical modeling of Trojans, which demonstrate that water ice should be depleted down to 10 m depth on average, though the sublimation front at polar regions may be considerably shallower, up to a depth 10 cm (Guilbert-Lepoutre, 2014). While a coating of fine-grained water frost similar to the proposed component on Themis can effectively suppress the 1.5 and 2 m absorption features (Rivkin & Emery, 2010), it produces an OH band at 3.1 m that is too narrow to match the observed feature on Trojans. Instead, larger micron-sized grains of water ice are necessary to reproduce the broad rounded shape and the short-wavelength edge at 2.7 m that are seen on Trojans, KBOs, and Centaurs, as well as the comets Tempel 1 (Sunshine et al., 2006), Hartley 2 (Protopapa et al., 2014), and 238P (Kelley et al., 2023).
Exposure of more pristine subsurface material through impact gardening has been proposed as a possible avenue by which otherwise unexpected water ice features may emerge on small dark bodies in the middle solar system. Spatially resolved observations of the irregular Saturnian satellite Phoebe, which has a similar size to the largest Trojans and is thought to have originated from the outer solar system, have shown that while water ice absorption is manifested across the entire surface, they are most prominent within impact basins and are correlated with patches of higher optical geometric albedo (Clark et al., 2005; Buratti et al., 2008; Fraser & Brown, 2018).
Extending this argument to the Trojans, which have experienced a similar level of collisional evolution as the main belt asteroids (Marchi et al., 2023), surface turnover from cratering impacts may periodically expose small localized patches of water-rich interior material amid the predominant water-ice-poor dark regolith. However, it is unclear if the present day impactor flux is sufficient to maintain pristine water ice on the surface, given the significantly higher rate of sublimation loss at the higher temperatures of the Trojan region (Guilbert-Lepoutre, 2014; Lisse et al., 2022). We note that isolated patches of water ice have been detected in spatially resolved spectra of the comets Tempel 1 (Sunshine et al., 2006) and 67P (Barucci et al., 2016), which experience much higher average surface temperatures than the Trojans. In those cases, the water ice is concentrated in regions that are topographically shielded from insolation or freshly exposed due to physical processes such as mass wasting. It follows that analogous surface configurations on the Trojans may preserve water ice for extended periods. In the context of this proposed scenario, the difference between the 3 m spectra of less-red and red Trojans would indicate that the former have intrinsically more water ice in their subsurface than the latter and suggests distinct initial surface compositions. The future Lucy flybys will provide spatially resolved optical and near-infrared photometry and spectroscopy of the Trojan targets’ surfaces and allow us to assess whether or not localized regions of exposed water ice are present on these objects.
An alternative explanation for the observed 3 m OH band on the less-red Trojans is organic refractory material formed by the irradiation of carbon-bearing ices. Laboratory experiments simulating the effect of solar processing on H2O CH3OH NH3 ice mixtures (Urso et al., 2020) produced complex residues that display a broad, asymmetric absorption band extending from 2.7–2.8 m to 3.6–3.8 m, in good agreement with the absorption feature seen on Trojans. This interpretation has also been presented to explain the absorption feature observed on 238P (Kelley et al., 2023). In this scenario, the OH band around 3 m is produced through the dissociation of H2O and CH3OH. Additional bands that contribute to extending the overall absorption feature to longer wavelengths are due to C–H and N–H stretch modes.
The KBO–Centaur–Trojan comparison has important implications for our understanding of solar system dynamical evolution. Numerical integrations that simulate the effect of current dynamical instability models on the primordial trans-Neptunian disk invariably show that the inward-scattered planetesimals, some of which ultimately become captured into resonance by Jupiter, pass through an intermediate stage when they lie on highly eccentric giant-planet-crossing orbits (e.g., Morbidelli et al., 2005; Nesvorný et al., 2013; Roig & Nesvorný, 2015). Centaurs therefore serve as an ideal proxy for assessing how Trojans may have evolved during their inward migration.
The spectrum of Okyrhoe is unique among the body of published JWST spectra of KBOs, with its extremely muted () 3 m OH band and the characteristic asymmetric shape of its broad 2.7–3.6 m absorption. Its spectrum also shows the strongest resemblance to the spectra of less-red Trojans out of all the small body spectra we have presented. Okyrhoe has one of the closest current orbits among the Centaurs hitherto observed by JWST, and we hypothesize that the difference between Okyrhoe and the more distant Centaur 2003 WL7, which has a spectrum that is much more consistent with that of typical KBOs, is indicative of differential evolution of Centaurs as a function of their received insolation. Within the framework of dynamical instability models of solar system evolution, we can apply this argument to the entire Trojan population, which has been located at 5 au for 4 Gyr. It follows that the cumulative thermal processing of Trojan surfaces may have further attenuated the 3 m OH band beyond the level apparent on Okyrhoe, producing the final spectra we see today.
It is important to note that the KBOs observed by JWST so far are significantly larger than the Trojans and Centaurs. As such, we are unable to disentangle possible size–dependent trends in the initial bulk composition and subsequent surface evolution of these objects directly. Future JWST observations of additional low-perihelion Centaurs and KBOs smaller than 300 km in diameter will enable a more robust assessment of the proposed evolutionary trajectory that links the KBOs and Centaurs to the Trojans and provide a definitive empirical test of the purported outer solar system origin of Trojans.
3.2 The 3.2–4.0 m region

Moving to slightly longer wavelengths, we find a distinct absorption band centered around 3.4 m. In contrast to the aforementioned 3 m band, this feature is more prominent on the red objects Orus and Leucus than it is on the less-red targets. To examine the shape of these features better, we fit polynomial continuum functions to the spectra from 3.1 to 4.2 m, excluding the region between 3.3 and 3.7 m. Figure 5 shows the continuum-subtracted spectra, where the 3.4 m band is clearly discernible on all targets except for Polymele; it is possible that the relatively poor signal-to-noise ratio of the Polymele spectrum may be obscuring a weak absorption in this wavelength range that has a similar depth to the analogous features seen on the other less-red Trojans. To characterize these bands quantitatively, we first divided the spectra by the local linear continuum fits, anchored to the edges of the band (3.25 and 3.57 m). We then used a Gaussian to model the band and followed the same 10,000-iteration Monte Carlo procedure as with the 3.0 m band. The resulting band center and depth values are reported in Table 3.
Analogous spectral features have been observed across a wide range of solar system small body populations, including main belt asteroids (e.g., Campins et al., 2010; Rivkin & Emery, 2010; Licandro et al., 2011; Takir & Emery, 2012; De Sanctis et al., 2018), the Jupiter-family comet 67P (Raponi et al., 2020), regular and irregular Jovian and Saturnian satellites (e.g., McCord et al., 1997; Cruikshank et al., 2008, 2014; Clark et al., 2012; Brown & Rhoden, 2014), and KBOs (e.g., Grundy et al., 2016; Stern et al., 2019; Brown & Fraser, 2023). In Figure 5, we plot the continuum-subtracted spectra of Ceres, Themis, 67P, Okyrhoe, and 2004 TY364. All of these spectra show a distinct 3.3–3.6 m absorption feature that is characteristic of C–H stretch modes from aliphatic organics. Several individual component bands in this region are resolved in the high-precision Rosetta spectrum of 67P (Raponi et al., 2020), including the asymmetric stretching modes of the aliphatic CH3 (3.38 m) and CH2 (3.42 m) groups and the blended symmetric stretching band of CH2 and CH3 (3.47 m). The presence of a distinct absorption feature at 3.3 m in the spectrum of 67P and the relative broadening of Ceres’ organic absorption band toward shorter wavelengths have been attributed to the presence of aromatic hydrocarbons in addition to the aliphatic organics (e.g., De Sanctis et al., 2019; Raponi et al., 2020). No evidence for aromatic hydrocarbons can be found in the Trojan, KBO, and Centaur spectra.
Laboratory experiments have demonstrated that aliphatic organic compounds are readily produced through ion bombardment and ultraviolet irradiation of ice mixtures containing various combinations of H2O, CH4, CH3OH, CO, NH3, and N2 (e.g., Strazzulla et al., 2001; Moore & Hudson, 2003; Palumbo et al., 2004; Cruikshank et al., 2005; Materese et al., 2014; Urso et al., 2020). An alternative pathway for forming aliphatic organics is through global aqueous alteration, a scenario that has been proposed for even small primitive objects, such as the parent bodies of the carbonaceous chondrite meteorites (e.g., Schulte & Shock, 2004; Vinogradoff et al., 2018). In addition to yielding the 3.3–3.6 m absorption bands observed throughout the middle and outer solar system, these processes also readily reproduce the reddened and darkened surfaces that are characteristic of Trojans, KBOs, Centaurs, and irregular satellites.
Within this interpretation, the systematically stronger aliphatic organic bands on Orus and Leucus suggest that the primordial surface composition of red Trojans was significantly richer in CH-bearing species than the less-red Trojans. The complex interplay between relative ice abundance, grain size, and radiation dose that is evident from these past experiments, as well as the multiplicity of possible chemical pathways that produce the same aliphatic organic compounds, makes it impossible to link the observed organic absorption on Trojans to a specific initial ice mixture. While initial surface ice inventories incorporating hypervolatile species such as N2 and CH4 necessitate a cold outer solar system origin, the ability of mixtures containing CH3OH to reproduce the broad 2.7–3.6 m spectral feature means that a middle solar system origin for Trojans cannot be fully excluded, as methanol ice can persist on airless surfaces at 5–20 au on kiloyear to megayear timescales (e.g., Wong & Brown, 2016; Lisse et al., 2022).
In the spectrum of Ceres, the 3.4 m band is accompanied by an additional major absorption feature at 3.7–4.0 m. Past modeling of Ceres’ surface has revealed that this pair of bands is indicative of organic compounds mixed with a significant amount of carbonate-bearing species, which absorb strongly at both 3.4 and 3.9–4.0 m (De Sanctis et al., 2015). The longer-wavelength band is particularly prominent within the bright regions in Occator Crater, where the absorption takes on a distinctive triangular shape with a reflectance minimum around 3.95 m (De Sanctis et al., 2016; Carrozzo et al., 2018). The intrinsic scatter in the Trojan spectra throughout this wavelength range is poorer than at shorter wavelengths, and there are signs of increased correlated noise features. We do not detect any hints of a broad dip in reflectance and place upper limits of a few percent on the presence of a carbonate absorption at 3.7–4.0 m.
3.3 4.25 m CO2 band

The reflectance spectrum of Eurybates contains a distinct absorption band at 4.25 m that is not seen on any other Trojan. This wavelength corresponds to the vibrational mode of CO2 ice (e.g., Sandford & Allamandola, 1990; Hansen, 1997). Analogous spectral features are ubiquitous across the middle and outer solar system, from giant planet satellites (e.g., McCord et al., 1998; Cruikshank et al., 2010; Pinilla-Alonso et al., 2011) to Centaurs and KBOs (e.g., Cruikshank et al., 1993; Harrington Pinto et al., 2022; de Prá et al., 2024; Pinilla-Alonso et al., 2024). Figure 6 shows a sample of reflectance spectra displaying 4.25 m absorption, plotted alongside the continuum-subtracted spectrum of Eurybates; the high spectral resolution JWST/NIRSpec spectra of Europa and Ganymede were taken from Trumbo & Brown (2023) and Trumbo et al. (2023), respectively. We analyzed the CO2 band shape on Eurybates using a similar methodology to our analysis of the 3.4 m band (Section 3.2). After subtracting the best-fit linear continuum through the edges of the band at 4.15 and 4.33 m, we carried out a Monte Carlo analysis using a Gaussian band profile. The band center and depth are reported in Table 3.
On the cold surfaces of KBOs and Centaurs, this band, centered around 4.268 m, is attributed to pure CO2 ice (e.g., de Prá et al., 2024). However, at the typical temperatures of the Jupiter Trojans, CO2 ice is extremely unstable, sublimating on day-long timescales (e.g., Lisse et al., 2022). Any exposed CO2 ice on the surface would be ephemeral, such as in the case of 67P, where isolated patches of CO2 ice in the southern hemisphere were found to persist during its cold season, when the local region was in shadow, before sublimating away on week-long timescales upon exposure to sunlight (Filacchione et al., 2016). Another example illustrating the instability of CO2 ice in the Trojan region is the recent JWST observation of the active Centaur 39P/Oterna, which revealed CO2 gas emission while the target was at 5.82 au (Harrington Pinto et al., 2023). It follows that the CO2 on the surface Eurybates must be bound to some refractory component. Type II CO2 clathrates (CO2 bound in a crystalline water ice structure) are a salient possibility and have been forwarded as an explanation for the 4.27 m band seen on Phoebe (Cruikshank et al., 2010). Laboratory studies have shown that the CO2 band center for type II clathrates drifts to shorter wavelengths with increasing temperature, reaching 4.265 m at the typical surface temperatures of Trojans (Blake et al., 1991). Meanwhile, shorter-wavelength absorption features may be due to mineral-bound CO2 (Hibbitts et al., 2000, 2003) or irradiation-formed carbonic acid (Jones et al., 2014), which have band centers near 4.25 m. Our measured band center for Eurybates ( m) is broadly consistent with both of these possibilities. CO2 can also be produced as a component of the organic-rich residue that forms from the irradiation of carbon-rich ice mixtures — the same process that yields the 3.4 m aliphatic organic features described in Section 3.2 (Urso et al., 2020).
The absence of a discernible CO2 feature on the other Trojan targets — both less-red and red — presents an interesting conundrum. Eurybates is the only member of a collisional family among the five targets observed with JWST, so we posit that its visible CO2 absorption is related to its impact history. It is important to note that the estimated age of the Eurybates family (1–4 Gyr; Marschall et al. 2022) is significantly longer than typical space weathering timescales for Trojans provided in the literature, which are on the order of yr (e.g., Melita et al., 2015). Therefore, with respect to irradiation-related alteration of surface materials, Eurybates is considered mature. It follows that the differences between Eurybates and the other Trojans must be due to fundamental distinctions in the physical and/or chemical properties of the surfaces.
Here, we once again use the unique spectrum of Okyrhoe to develop a hypothesis for this phenomenon. While Okyrhoe shows a muted 3 m feature with a depth that is intermediate between that of Trojans and KBOs, its spectrum lacks a discernible CO2 absorption band. A subset of Centaurs has been observed to display cometary activity, with the activity concentrated among low-perihelion objects (e.g., Jewitt, 2009, 2015; Wong et al., 2019b). The exact trigger for the onset of Centaur activity is still not fully understood, but a leading hypothesis points toward the transition between amorphous to crystalline water ice and the concomitant release of trapped volatiles (e.g., Capria et al., 2000; Jewitt, 2009; Wierzchos et al., 2017). While Okyrhoe has never been observed to be active, its relatively close orbit places it in a region where the amorphous-to-crystalline transition can take place (Guilbert-Lepoutre, 2012). Physical modeling of cometary activity has demonstrated that fine-grained refractory material is entrained in the outgassed volatiles and can be redeposited across the surface (e.g., Jewitt, 2002). The obscuration of Okyrhoe’s surface due to past activity may explain the distinctly muted 3 m OH band and the absence of 4.25 m CO2 absorption in its spectrum.
If Trojans initially formed as icy planetesimals in the primordial trans-Neptunian region, as current dynamical instability models suggest, they likely became active during their inward migration from the outer solar system. The increased temperatures and space weathering would have led to thermal processing of the initially volatile-rich surfaces. Due to the observed prevalence of CO2 ice across the KBO population (de Prá et al., 2024; Pinilla-Alonso et al., 2024), it is possible that all Trojans developed surfaces rich in bound CO2 (e.g., through the aforementioned irradiation process). Meanwhile, the activated outgassing at these higher temperatures would have caused simultaneous physical obscuration of the spectral features across the population. In the special case of Eurybates, the later family forming impact destroyed the deposited lag material on the parent body and exposed material from deeper layers, thereby providing an unobscured view of the 4.25 m bound CO2 band. Analogous JWST observations of other members of the Eurybates family will enable us to determine whether the CO2 feature is indeed a distinctive feature of the collisional fragments. Recent photometric observations have confirmed a second collisional family within the Trojans — the Ennomos family (Wong & Brown, 2023). Comparisons of the two families could assess whether there are systematic differences in the way CO2 is bound to the surfaces of Trojans.
4 Conclusions
We have presented new 1.7–5.3 m JWST/NIRSpec reflectance spectra of the five Jupiter Trojans that will be visited by the Lucy spacecraft in 2027–2033. These observations provide enhanced sensitivity and wavelength coverage relative to previously published results. Across the five spectra, we have identified and characterized three major absorption bands, which we summarize below.
-
1.
A broad OH absorption feature centered around 3 m is discernible across the five targets. The feature is systematically deeper on the less-red Trojans, particularly Eurybates, and has a sharp short-wavelength edge at roughly 2.7 m. This band is inconsistent with the sharper phyllosilicate feature seen on many large main belt asteroids (e.g., Ceres, Pallas, and Themis), instead being more comparable to the rounded absorption features that are characteristic of water-ice-rich Centaurs and KBOs, as well as the active outer main belt asteroid 238P. While small patches of exposed water ice may account for the presence of this OH band on Trojans, an alternative explanation is that the OH-bearing components within an irradiation residue were formed from the solar processing of volatile ice mixtures rich in carbon and OH-bearing molecules such as methanol.
-
2.
Between 3.3 and 3.6 m, we find a distinct absorption band that is indicative of aliphatic organic compounds. The depth of this feature is stronger on the red Trojans Orus and Leucus than on the less-red object, suggesting that the red Trojans initially formed with a higher relative abundance of CH-bearing species on their surfaces. We do not find evidence for carbonates on the Trojans, which lack the additional absorption band at 3.7–4.0 m.
-
3.
A 4.25 m absorption band is present only in the spectrum of Eurybates. This feature indicates the presence of CO2 on the surface; as pure CO2 ice is highly unstable to sublimation loss at 5.2 au, the CO2 on Eurybates is likely bound to or trapped within some less volatile material, e.g., water ice clathrates. The observation that the only target to display this CO2 feature is a collisional fragment leads us to hypothesize that CO2 may be ubiquitous across the Trojan population, with the spectral signature being obscured on the uncollided objects due to some overlaying lag material.
Although these JWST spectra have significantly enriched our understanding of the Trojans’ surface composition, we do not yet have definitive proof that exclusively links these objects to the outer solar system, as predicted by recent dynamical instability models. While all of the three spectral features outlined above are commonplace throughout the KBO and Centaur population, the existence of other nearby small bodies with analogous features (e.g., CO2 on Jovian and Saturnian satellites) demonstrates that there may be active chemical processes within the giant planet region that are capable of producing the observed spectral features on Trojans. Nevertheless, the tentative trend between decreasing OH band strength and decreasing heliocentric distance across the KBOs, Centaurs, and Trojans suggests a unified evolutionary trajectory. In particular, the strongly attenuated OH absorption on the close-in Centaur Okyrhoe, with a band depth that is intermediate between those seen on KBOs and Trojans, points toward secondary thermal processing as a potential explanation for the observed variation in spectral shape. Possible activation of outgassing on the Trojans during their purported inward migration, analogous to the process on present day active Centaurs, may also account for the absence of a discernible CO2 band on all but the collisional fragments within the Trojan population.
Looking ahead to the Lucy flybys, we expect that these JWST/NIRSpec spectra will continue to provide crucial insights. For example, because the onboard instruments do not cover the 4–5 m region, our knowledge of the unique CO2 band on Eurybates will be indispensable when assessing any systematic differences that are observed at shorter wavelengths in comparison to the other flyby targets. Conversely, the fine spatial resolution provided by Lucy across the Trojans’ surfaces will refine our interpretations of the JWST spectra. The spacecraft will be able to probe for isolated patches of exposed water ice and yield compositional maps that link the OH and aliphatic organic bands to specific geological features. Likewise, using thermal infrared measurements, we will unveil the detailed physical properties of the regolith and search for the telltale signs of past cometary activity on Trojans. Ultimately, the findings from the Lucy mission must be synthesized with spectroscopic results for Trojans, Centaurs, KBOs, and other small bodies in order to place the Trojans within the broader picture of solar system formation and evolution properly.
References
- Astropy Collaboration (2022) Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167
- Astropy Collaboration (2013) Astropy Collaboration, Price-Whelan, A. M., Sipőcz, B. M., et al. 2018, AJ, 156, 123
- Astropy Collaboration (2018) Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33
- Baratta & Palumbo (1998) Baratta, G., & Palumbo, M. 1998, JOSAA, 15, 3076
- Barkume et al. (2008) Barkume, K. M., Brown, M. E., & Schaller, E. L. 2008, AJ, 135, 55
- Barucci et al. (2011) Barucci, M. A., Alvarez-Candal, A., Merlin, F., et al. 2011, Icar, 214, 297
- Barucci et al. (2016) Barucci, M. A., Filacchione, G., Fornasier, S., et al. 2016, A&A, 595, A102
- Blake et al. (1991) Blake, D., Allamandola, L., Sandford, S., Hudgins, D., & Freund, F. 1991, Sci, 254, 548
- Böker et al. (2023) Böker, T., Beck, T. L., Birkmann, S. M., et al. 2023, PASP, 135, 038001
- Brown (2016) Brown, M. E. 2016, AJ, 152, 159
- Brown & Fraser (2023) Brown, M. E., & Fraser, W. C. 2023, PSJ, 4, 130
- Brown & Rhoden (2014) Brown, M. E., & Rhoden, A. R. 2014, ApJ, 793, L44
- Brož & Rozehnal (2011) Brož, M., & Rozehnal, J. 2011, MNRAS, 414, 565
- Buie et al. (2021) Buie, M. W., Keeney, B. A., Strauss, R. H., et al. 2021, PSJ, 2, 202
- Buie et al. (2018) Buie, M. W., Zangari, A. M., Marchi, S., Levison, H. F., & Mottola, S. 2018, AJ, 155, 245
- Buratti et al. (2008) Buratti, B. J., Soderlund, K., Bauer, J., et al. 2008, Icar, 193, 309
- Bushouse et al. (2022) Bushouse, H., Eisenhamer, J., Dencheva, N., et al. 2023, JWST Calibration Pipeline, v1.11.3, Zenodo, doi:10.5281/zenodo.8157276
- Campins et al. (2010) Campins, H., Hargrove, K., Pinilla-Alonso, N., et al. 2010, Natur, 464, 1320
- Capria et al. (2000) Capria, M. T., Coradini, A., De Sanctis, M. C., & Orosei, R. 2000, AJ, 119, 3112
- Carrozzo et al. (2018) Carrozzo, F. G., De Sanctis, M. C., Raponi, A., et al. 2018, SciA, 4, e1701645
- Clark et al. (2005) Clark, R. N., Brown, R. H., Jaumann, R. et al. 2005, Natur, 435, 66
- Clark et al. (2012) Clark, R. N., Cruikshank, D. P., Jaumann, R., et al. 2012, Icar, 218, 831
- Cruikshank et al. (2014) Cruikshank, D. P., Dalle Ore, C. M., Clark, R. N., & Pendelton, Y. J. 2014, Icar, 233, 306
- Cruikshank et al. (2005) Cruikshank, D. P., Imanaka, H., & Dalle Ore, C. M. 2005, AdSpR, 36, 178
- Cruikshank et al. (2010) Cruikshank, D. P., Meyer, A. W., Brown, R. H., et al. 2010, Icar, 206, 561
- Cruikshank et al. (1993) Cruikshank, D. P., Roush, T. L., Owen, T. C., et al. 1993, Sci, 261, 742
- Cruikshank et al. (2008) Cruikshank, D. P., Wegryn, E., Dalle Ore, C. M., et al. 2008, Icar, 193, 334
- de Prá et al. (2024) de Prá, M. N., Hénault, E., Pinilla-Alonso, N., et al. 2024, NatAs, in press
- De Sanctis et al. (2018) De Sanctis, M. C., Ammannito, E., Carrozzo, F. G., et al. 2018, M&PS, 53, 1844
- De Sanctis et al. (2015) De Sanctis, M. C., Ammannito, E., Raponi, A., et al. 2015, Natur, 528, 241
- De Sanctis et al. (2016) De Sanctis, M. C., Raponi, A., Ammannito, E., et al. 2016, Natur, 536, 54
- De Sanctis et al. (2019) De Sanctis, M. C., Vinogradoff, V., Raponi, A., et al. 2019, MNRAS, 482, 2407
- Dotto et al. (2006) Dotto, E., Fornasier, S., Barucci, M. A., et al. 2006, Icar, 183, 420
- Duffard et al. (2014) Duffard, R., Pinilla-Alonso, N., Santos-Sanz, P., et al. 2014, A&A, 564, A92
- Emery & Brown (2003) Emery, J. P., & Brown, R. H. 2003, Icar, 164, 104
- Emery et al. (2011) Emery, J. P., Burr, D. M., & Cruikshank, D. P. 2011, AJ, 141, 25
- Emery et al. (2024) Emery, J. P., Wong, I., Brunetto, R., et al. 2024, Icar, in press (arXiv:2309.15230)
- Filacchione et al. (2016) Filacchione, G., Raponi, A., Capaccioni, R., et al. 2016, Sci, 354, 1563
- Fornasier et al. (2007) Fornasier, S., Dotto, E., Hainaut, O., et al. 2007, Icar, 190, 622
- Fraser & Brown (2012) Fraser, W. C., & Brown, M. E. 2012, ApJ, 749, 33
- Fraser & Brown (2018) Fraser, W. C., & Brown, M. E. 2018, AJ, 156, 23
- Fraser et al. (2023) Fraser, W. C., Pike, R. E., Marsset, M., et al. 2023, PSJ, 4, 80
- Gomes et al. (2005) Gomes, R., Levison, H. F., Tsiganis, K., & Morbidelli, A. 2005, Natur, 435, 466
- Grundy et al. (2016) Grundy, W. M., Binzel, R. P., Buratti, B. J., et al. 2016, Sci, 351, 9189
- Grundy et al. (2018) Grundy, W. M., Noll, K. S., Buie, M. W., & Levison, H. F. 2018, Icar, 305, 198
- Grundy et al. (2024) Grundy, W. M., Wong, I., Glein, C. R., et al. 2024, Icar, 411, 115923
- Guilbert-Lepoutre (2012) Guilbert-Lepoutre, A. 2012, AJ, 144, 97
- Guilbert-Lepoutre (2014) Guilbert-Lepoutre, A. 2014, Icar, 231, 232
- Haghighipour (2009) Haghighipour, N. 2009, M&PS, 44, 1863
- Hansen (1997) Hansen, G. B. 1997, JGR, 102, 21569
- Harrington Pinto et al. (2023) Harrington Pinto, O., Kelley, M. S. P., Villanueva, G. L., et al. 2023, PSJ, 4, 208
- Harrington Pinto et al. (2022) Harrington Pinto, O., Womack, M., Fernandez, Y., & Bauer, J. 2022, PSJ, 3, 247
- Harris (1998) Harris, A. W. 1998, Icar, 131, 291
- Harris et al. (2020) Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Natur, 585, 357
- Hibbitts et al. (2000) Hibbitts, C. A., McCord, T. B., & Hansen, G. B. 2000, JGR, 105, 22541
- Hibbitts et al. (2003) Hibbitts, C. A., Pappalardo, R. T., Hansen, G. B., & McCord, T. B. 2003 JGRE, 108, 5036
- Hsieh et al. (2016) Hsieh, H. H., & Haghighipour, N. 2016, Icar, 277, 19
- Hunter (2007) Hunter, J. D. 2007, CSE, 9, 90
- Jakobsen et al. (2022) Jakobsen, P., Ferruit, P., Alves de Oliveria, C., et al. 2022, A&A, 661, A80
- Jewitt (2009) Jewitt, D. 2009, AJ, 137, 4296
- Jewitt (2015) Jewitt, D. 2015, AJ, 150, 201
- Jewitt (2002) Jewitt, D., C., 2002, AJ, 123, 1039
- Jones et al. (2014) Jones, B. M., Kaiser, R. I., & Strazzulla, G. 2014, ApJ, 788, 170
- Jones et al. (1990) Jones, T. D., Lebofsky, L. A., Lewis, J. S., & Marley, M. S. 1990, Icar, 88, 172
- Kelley et al. (2023) Kelley, M. S. P., Hsieh, H. H., Bodewits, D., et al. 2023, Natur, 619, 720
- King et al. (1992) King, T. V. V., Clark, R. N., Calvin, W. M., Sherman, D. M., & Brown, R. H. 1992, Sci, 255, 1551
- Lebofsky (1980) Lebofsky, L. A. 1980, AJ, 85, 573
- Lellouch et al. (2013) Lellouch, E., Santos-Sanz, P., Lacerda, P., et al. 2013, A&A, 557, A60
- Levison et al. (2008) Levison, H. F., Morbidelli, A., Van Laerhoven, C., Gomes, R., & Tsiganis, K. 2008, Icar, 196, 258
- Levison et al. (2021) Levison, H. F., Olkin, C. B., Noll, K. S., et al. 2021, PSJ, 2, 171
- Licandro et al. (2011) Licandro, J., Campins, H., Kelley, M., et al. 2011, A&A, 525, A34
- Licandro et al. (2024) Licandro, J., Pinilla-Alonso, N., Holler, B. J., et al. 2024, NatAs, in revision
- Lisse et al. (2022) Lisse, C. M., Gladstone, G. R., Young, L. A., et al. 2022, PSJ, 3, 112
- Marchi et al. (2023) Marchi, S., Nesvorný, D., Vokrouhlický, D., Bottke, W. F., & Levison, H. F. 2023, AJ, 166, 221
- Marschall et al. (2022) Marschall, R., Nesvorný, D., Deienno, R., et al. 2022, AJ, 164, 167
- Marzari & Scholl (1998) Marzari, F., & Scholl, H. 1998, A&A, 339, 278
- Marzari et al. (2003) Marzari, F., Tricarico, P., & Scholl, H. 2003, Icar, 162, 453
- Mastrapa et al. (2009) Mastrapa, R. M., Sandford, S. A., Roush, T. L., Cruikshank, D. P., & Dalle Ore, C. M. 2009, ApJ, 701, 1347
- Mastrapa et al. (2013) Mastrapa, R. M. E., Grundy, W. M., & Gudipati, M. S. 2013, in The Science of Solar System Ices, ed. M. S. Gudipati & J. Castillo-Rogez (Berlin: Springer), 371
- Materese et al. (2014) Materese, C. K., Cruikshank, D. P., Sandford, S. A., et al. 2014, ApJ, 788, 111
- McCord et al. (1997) McCord, T. B., Carlson, R., Smythe, W., et al. 1997, Sci, 278, 271
- McCord et al. (1998) McCord, T. B., Hansen, G. B., Clark, R. N., et al. 1998, JGR, 103, 8603
- McGraw et al. (2022) McGraw, L. E., Emery, J. P., Thomas, C. A., et al. 2022, PSJ, 3, 243
- McKay et al. (2017) McKay, A. J., Bodewits, D., & Li, J.-Y. 2017, Icar, 286, 308
- Melita et al. (2015) Melita, M. D., Kaňuchová, Z., Brunetto, R., & Strazzulla, G. 2015, Icar, 248, 222
- Moore & Hudson (2003) Moore, M. H., & Hudson, R. L. 2003, Icar, 161, 486
- Morbidelli et al. (2005) Morbidelli, A., Levison, H. F., Tsiganis, K., & Gomes, R. 2005, Natur, 435, 462
- Moseley et al. (2010) Moseley, S., Arendt, R. G., Fixsen, D. J., et al. 2010, SPIE, 7742, 77421B
- Mottola et al. (2023) Mottola, S., Hellmich, S., Buie, M. W., et al. 2023, PSJ, 4, 18
- Mueller et al. (2010) Mueller, M., Marchis, F., Emery, J. P., et al. 2010, Icar, 205, 505
- Nesvorný et al. (2013) Nesvorný, D., Vokrouhlick]’y, D., & Morbidelli, A. 2013, ApJ, 768, 45
- Olkin et al. (2021) Olkin, C. B., Levison, H. F., Vincent, M., et al. 2021, PSJ, 2, 172
- O’Rourke et al. (2020) O’Rourke, L., Müller, T. G., Biver, N., et al. 2020, ApJ, 898, L45
- Palumbo et al. (2004) Palumbo, M. E., Ferini, G., & Baratta, G. A. 2004, AdSpR, 33, 49
- Perrin et al. (2014) Perrin, M. D., Sivaramakrishnan, A., Lajoie, C.-P., et al. 2014, Proc. SPIE, 9143, 91433X
- Pinilla-Alonso et al. (2024) Pinilla-Alonso, N., Brunetto, R., de Prá, M. N., et al. 2024, NatAs, in revision
- Pinilla-Alonso et al. (2011) Pinilla-Alonso, N., Roush, T. L., Marzo, G. A., et al. 2011, Icar, 215, 75
- Poch et al. (2020) Poch, O., Istiqomah, I., Quirico, E., et al. 2020, Sci, 367, 1212
- Protopapa et al. (2014) Protopapa, S., Sunshine, J. M., Feaga, L. M., et al. 2014, Icar, 238, 191
- Raponi et al. (2020) Raponi, A., Ciarniello, M., Capaccioni, F., et al. 2020, NatAs, 4, 500
- Quirico & Schmitt (1997) Quirico, E., & Schmitt, B. 1997, Icar, 127, 354
- Rauscher et al. (2012) Rauscher, B. J., Arendt, R. G., Fixsen, D. J., et al. 2012, Proc. SPIE, 8453, 84531F
- Rivkin & Emery (2010) Rivkin, A. S., & Emery, J. P. 2010, Natur, 464, 1322
- Rivkin et al. (2002) Rivkin, A. S., Howell, E. S., Vilas, F., & Lebofsky, L. A. 2002, in Asteroids III, ed. W. F. Bottke Jr., et al. (Univ. of Arizona), 235
- Roig & Nesvorný (2015) Roig, F., & Nesvorný, D. 2015, AJ, 150 186
- Roig et al. (2008) Roig, F., Ribeiro, A. O., & Gil-Hutton, R. 2008, A&A, 483, 911
- Sandford & Allamandola (1990) Sandford, S. A., & Allamandola, L. J. 1990, ApJ, 355, 357
- Schulte & Shock (2004) Schulte, M., & Shock, E. 2004, M&PS, 39, 1577
- Sharkey et al. (2019) Sharkey, B. N. L., Reddy, V., Sanchez, J. A., Izawa, M. R. M., & Emery, J. P. 2019, AJ, 158, 204
- Stern et al. (2019) Stern, S. A., Weaver, H. A., Spencer, J. R., et al. 2019, Sci, 364, 9771
- Strazzulla et al. (2001) Strazzulla, G., Baratta, G. A., & Palumbo, M. E. 2001, AcSpA, 57, 825
- Sunshine et al. (2006) Sunshine, J. M., A’Hearn, M. F., Groussin, O., et al. 2006, Sci, 311, 1453
- Szabó et al. (2007) Szabó, Gy. M., Ivezić, Ž., Jurić, M., & Lupton, R. 2007, MNRAS, 377, 1393
- Takir & Emery (2012) Takir, D., & Emery, J. P. 2012, Icar, 219, 641
- Trumbo & Brown (2023) Trumbo, S. K., & Brown, M. E. 2023, Sci, 381, 1308
- Trumbo et al. (2023) Trumbo, S. K., Brown, M. E., Bockelée-Morvan, D., et al. 2023, SciA, 9, 3724
- Tsiganis et al. (2005) Tsiganis, K., Gomes, R., Morbidelli, A., & Levison, H. F. 2005, Natur, 435, 459
- Urso et al. (2020) Urso, R. G., Vuitton, V., Danger, G., et al. 2020, A&A, 644, A115
- Usui et al. (2019) Usui, F., Hasegawa, S., Ootsubo, T., & Onaka, T. 2019, PASJ, 71, 1
- Vilenius et al. (2014) Vilenius, E., Kiss, C., Müller, T., et al. 2014, A&A, 564, A35
- Vinogradoff et al. (2018) Vinogradoff, V., Bernard, S., Le Guillou, C., & Remusat, L. 2018, Icar, 305, 358
- Virtanen et al. (2020) Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, NatMe, 17, 261
- Vokrouhlický et al. (2016) Vokrouhlický, D., Bottke, W. F., & Nesvorný, D. 2016, AJ, 152, 39
- Wierzchos et al. (2017) Wierzchos, K., Womack, M., & Sarid, G. 2017, AJ, 153, 230
- Wong & Brown (2015) Wong, I., & Brown, M. E. 2015, AJ, 150, 174
- Wong & Brown (2016) Wong, I., & Brown, M. E. 2016, AJ, 152, 90
- Wong & Brown (2017) Wong, I., & Brown, M. E. 2017, AJ, 153, 145
- Wong & Brown (2019) Wong, I., & Brown, M. E. 2019, AJ, 157, 203
- Wong et al. (2019a) Wong, I., Brown, M. E., Blacksberg, J., et al. 2019, AJ, 157, 161
- Wong et al. (2019b) Wong, I., Mishra, A., & Brown, M. E., 2019, AJ, 157, 225
- Wong & Brown (2023) Wong, I., & Brown, M. E. 2023, AJ, 165, 15
- Wong et al. (2014) Wong, I., Brown, M. E., & Emery, J. P. 2014, AJ, 148, 112
- Yang & Jewitt (2007) Yang, B., & Jewitt, D. 2007, AJ, 134, 223