# Col-OSSOS: Color and Inclination Are Correlated throughout the Kuiper Belt

Marsset M., Fraser W. C., Pike R. E., Bannister M. T., Schwamb M. E., Volk K., Kavelaars J. J., Alexandersen M., Chen Y.-T., Gladman B. J., et al.

► To cite this version:

Marsset M., Fraser W. C., Pike R. E., Bannister M. T., Schwamb M. E., et al. Col-OSSOS: Color and Inclination Are Correlated throughout the Kuiper Belt. The Astronomical Journal, 2019, 157 (3), pp.94. 10.3847/1538-3881/aaf72e hal-02372817

HAL Id: hal-02372817

https://hal.science/hal-02372817v1

Submitted on 20 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Col-OSSOS: Color and Inclination Are Correlated throughout the Kuiper Belt

Michaël Marsset<fnref n="1" />,<fnref n="2" /> ID, Wesley C. Fraser<fnref n="1" />,<fnref n="3" /> ID, Rosemary E. Pike<fnref n="4" /> ID, Michele T. Bannister<fnref n="1" />,<fnref n="3" />,<fnref n="5" /> ID, Megan E. Schwamb<fnref n="6" /> ID, Kathryn Volk<fnref n="7" /> ID, J. J. Kavelaars<fnref n="3" />,<fnref n="5" /> ID, Mike Alexandersen<fnref n="4" /> ID, Ying-Tung Chen<fnref n="4" /> ID, Brett J. Gladman<fnref n="8" />, Stephen D. J. Gwyn<fnref n="3" />, Matthew J. Lehner<fnref n="4" />,<fnref n="9" />,<fnref n="10" /> ID, Nuno Peixinho<fnref n="11" /> ID, Jean-Marc Petit<fnref n="12" /> ID, and Shiang-Yu Wang<fnref n="4" />

## Abstract

Both physical and dynamical properties must be considered to constrain the origins of the dynamically excited distant solar system populations. We present high-precision (g-r) colors for 25 small ($H_r > 5$) dynamically excited trans-Neptunian objects (TNOs) and centaurs acquired as part of the Colours of the Outer Solar System Origins Survey. We combine our data set with previously published measurements and consider a set of 229 colors of outer solar system objects on dynamically excited orbits. The overall color distribution is bimodal and can be decomposed into two distinct classes, termed gray and red, that each has a normal color distribution. The two color classes have different inclination distributions: red objects have lower inclinations than the gray ones. This trend holds for all dynamically excited TNO populations. Even in the worst-case scenario, biases in the discovery surveys cannot account for this trend; it is intrinsic to the TNO population. Considering that TNOs are the precursors of centaurs, and that their inclinations are roughly preserved as they become centaurs, our finding solves the conundrum of centaurs being the only outer solar system population identified so far to exhibit this property. The different orbital distributions of the gray and red dynamically excited TNOs provide strong evidence that their colors are due to different formation locations in a disk of planetesimals with a compositional gradient.

**Key words:** Kuiper belt: general – minor planets, asteroids: general – surveys

**Supporting material:** machine-readable table

## 1. Introduction

Trans-Neptunian objects (TNOs) represent some of the most unaltered remnants of the planetary formation process. The surface of these objects can be used as tracers of their dynamical evolution which, in turn, helps us understand the current structure of the Kuiper Belt (e.g., Gladman 2005; Petit et al. 2011). Numerical models have been used to assess the evolution of TNOs from their formation locations onto their current orbits (e.g., Malhotra 1995; Gomes 2003; Levison et al. 2008; Nesvorny 2015a, 2015b). Different formation locations might correspond to distinct surface compositions, though this remains disputed (Gil-Hutton 2002; Stern 2002; Trujillo & Brown 2002; Delsanti et al. 2004; Santos-Sanz et al. 2009). As such, both physical and dynamical properties must be considered to constrain the origins of TNOs. Unlike the majority of the known dwarf-planet-sized bodies, the smaller TNOs are too faint to be studied through optical and infrared spectroscopy.

For these objects, one must instead rely on what broadband colors reveal.

The TNOs can be divided into two broad dynamical groups of objects: the dynamically quiescent cold classicals, and the dynamically excited hot population (Tegler & Romanishin 2000; Brown 2001). The latter comprises several subpopulations: (1) the hot classicals located at 40–48 au between Neptune’s 3:2 and 2:1 mean-motion orbital resonances (MMRs; Levison & Stern 2001; Doressoundiram et al. 2002; Gladman et al. 2008), (2) the resonant objects locked in MMRs with Neptune (Malhotra 1995; Levison et al. 2008; Gladman et al. 2012), (3) the scattering population on dynamically unstable and highly eccentric orbits (Lykawka & Mukai 2007; Gladman et al. 2008; Gomes et al. 2008), and (4) the detached TNOs that experienced gravitational excitation in the past and are now on stable orbits unperturbed by Neptune (Gladman et al. 2002; Emel'yanenko et al. 2003; Gomes et al. 2005; Brasser & Schwamb 2015). Centaurs are thought to be former TNOs that evolved into inner orbits with semimajor axes and perihelia between the orbits of Jupiter and Neptune (e.g., Levison & Duncan 1997). We consider them here as part of the larger group of dynamically excited TNOs.

There are indications that colors and orbital properties of the dynamically excited TNOs are related. The largest dynamically

<figure><img src="image_3.png" /></figure>

Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

<footnote n="1">¹ Astrophysics Research Centre, Queen's University Belfast, Belfast BT7 1NN, UK; mmarsset@mit.edu</footnote>

<footnote n="2">² Department of Earth, Atmospheric and Planetary Sciences, MIT, 77 Massachusetts Avenue, Cambridge, MA 02139, USA</footnote>

<footnote n="3">³ NRC-Herzberg Astronomy and Astrophysics, National Research Council of Canada, 5071 West Saanich Road, Victoria, BC V9E 2E7, Canada</footnote>

<footnote n="4">⁴ Institute of Astronomy and Astrophysics, Academia Sinica; 11F of AS/NTUAstronomy-Mathematics Building, No.1, Section 4, Roosevelt Road, Taipei 10617, Taiwan, Republic of China</footnote>

<footnote n="5">⁵ Department of Physics and Astronomy, University of Victoria, Elliott Building, 3800 Finnerty Road, Victoria, BC V8P 5C2, Canada</footnote>

<footnote n="6">⁶ Gemini Observatory, Northern Operations Center, 670 North A'ohoku Place, Hilo, HI 96720, USA</footnote>

<footnote n="7">⁷ Lunar and Planetary Laboratory, The University of Arizona, 1629 E. University Boulevard, Tucson, AZ 85721, USA</footnote>

<footnote n="8">⁸ Department of Physics and Astronomy, University of British Columbia, Vancouver, BC V6T 1Z1, Canada</footnote>

<footnote n="9">⁹ Department of Physics and Astronomy, University of Pennsylvania, 209 S. 33rd Street, Philadelphia, PA 19104, USA</footnote>

<footnote n="10">¹⁰ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA</footnote>

<footnote n="11">¹¹ CITEUC—Centre for Earth and Space Science Research of the University of Coimbra, Geophysical and Astronomical Observatory of the University of Coimbra, 3040-004 Coimbra, Portugal</footnote>

<footnote n="12">¹² Institut UTINAM UMR6213, CNRS, Univ. Bourgogne Franche-Comté, OSU Theta F-25000 Besançon, France</footnote>

<footnote>Received 2018 September 7; revised 2018 December 4; accepted 2018 December 5; published 2019 February 4</footnote>
excited TNOs uniformly exhibit neutral, ice-dominated surfaces (Fornasier et al. 2004; Barkume et al. 2008; Barucci et al. 2008, 2011; Brown 2012), possibly because they have sufficient gravity to retain their volatiles (Schaller & Brown 2007). In contrast, smaller dynamically excited TNOs have a bimodal color distribution, with classes termed gray and red that are divided around optical color (g-r) = 0.78 (Tegler & Romanishin 1998, 2000, 2003; Peixinho et al. 2003, 2012, 2015; Tegler et al. 2003, 2008, 2016; Fraser & Brown 2012; Lacerda et al. 2014; Wong & Brown 2017); this is fully quantified in Section 3. A correlation between colors and orbital inclinations was found for the population of classical TNOs (Tegler & Romanishin 2000; Hainaut & Delsanti 2002; Trujillo & Brown 2002; Doressoundiram et al. 2005; Peixinho et al. 2008). However, the data sets in these works were contaminated by the presence in their sample of peculiar compositional groups, including neutral-color Haumea family members (Brown et al. 2007; Ragozzine & Brown 2007; Schaller & Brown 2008; Snodgrass et al. 2010) and red cold classical objects (Tegler & Romanishin 2000; Brown 2001). A careful removal of these objects diminishes the strength of the correlation between color and inclination, but the correlation remains statistically significant (Peixinho et al. 2015). More recently, Tegler et al. (2016) divided the observed population of centaurs into two groups by color and analyzed their orbital inclination distributions. They found that red centaurs have lower inclinations than the gray ones. As minor planets on trans-neptunian orbits are perturbed into the centaur region, they roughly preserve the inclinations they had before becoming centaurs (Volk & Malhotra 2013). As such, the different orbital distributions of the two color classes of centaurs should be reflective of a similar trend in the more distant TNOs. Such a trend might be the origin of the observed color versus inclination correlation of the hot classical population. Two color groups with different inclination distributions could indeed manifest as such a correlation. In short, different inclinations for the gray and red dynamically hot objects are anticipated but not yet observed in the Kuiper Belt.

We re-examine here the orbital inclination distribution of the two color classes of the dynamically excited TNOs. We use a large data set of high-quality color measurements that includes optical (g-r) colors obtained for 44 objects as part of the Colours of the Outer Solar System Origin Survey (Col-OSSOS; Schwamb et al. 2018), 25 of which are new measurements (Section 2). We further use robust color measurements extracted from several works reporting optical colors of TNOs, bringing our sample to 229 dynamically excited TNOs. We show in Section 3 that red dynamically excited TNOs have smaller orbital inclinations than their gray counterparts and that this trend appears to be present in every dynamical population. The centaurs are no longer the only population to exhibit this peculiarity. This has strong implications for the origins of the dynamically excited TNO populations (Section 4).

## 2. Sample Selection, Observations, and Color Measurements

### 2.1. Sample Selection

Our analysis uses both newly obtained colors and color measurements from the literature for TNOs confirmed to be on dynamically excited orbits. We consider the subset of 44 dynamically excited objects from Outer Solar System Origin Survey (OSSOS; Bannister et al. 2016, 2018) for which our

team acquired optical and near-infrared colors as part of the Col-OSSOS survey (Schwamb et al. 2018). Measurements for 19 of these objects are previously published colors from Schwamb et al. (2018), and 25 are new color measurements reported here. Those include measurements for four objects that we previously published in Pike et al. (2017) and reprocessed here using the latest version of our photometric pipeline. We complement this Col-OSSOS data set with available high-quality color measurements from the literature. These include optical colors from Peixinho et al. (2015), which compiles data from many authors, and Tegler et al. (2016), as well as *Hubble Space Telescope* (HST) data from Fraser & Brown (2012) and Fraser et al. (2015). Only colors observed with filters in the ~500–850 nm range were considered, i.e., Johnson-Cousins (V-R) and (R-I) and HST (F606w–F814w). One object from Peixinho et al. (2015), 2002 GO9 (Crantor), has (R-I) color inconsistent with its (V-R) color and measured spectra (e.g., Alvarez-Candal et al. 2007). We consider only the (V-R) measurement reported in Peixinho et al. (2015) for that target.

Only small objects are considered for this analysis. The size transition from ice-rich surfaces on large TNOs to potentially volatile depleted surfaces on small TNOs is a matter of debate, but the transition is generally proposed to be at absolute magnitude $H_r \sim 6$ (Fraser & Brown 2012; Peixinho et al. 2012). Here, we find that including objects in the $5 < H_r < 6$ mag range does not affect the conclusions of our analysis (see Section 3). We therefore adopt a magnitude cut of $H_r > 5.0$, corresponding to diameters of $\lesssim 440$ km assuming a visible geometrical albedo of 0.09 (Lacerda et al. 2014).

We classified all of the TNOs in our sample using the procedure outlined in Gladman et al. (2008), which we briefly review here. Each object's orbit was fit using the Bernstein & Khushalani (2000) orbit-fitting routine. Starting from that best-fit orbit, a Monte-Carlo search was performed for the minimum and maximum semimajor axis orbits that do not produce orbit-fit residuals more than 1.5 times worse than those of the best-fit orbit. All three orbits (the best-fit, minimum-a, and maximum-a) were integrated forward in time for 10 Myr and classified according to their evolution. A secure classification is one for which all three orbits behave the same way, and insecure classification is one where at least one orbit behaves differently. We have 202 secure classifications and 27 insecure classifications. Most of the insecure classifications are a result of an object being in or near an MMR rather than having a particularly large uncertainty in semimajor axis because all of the TNOs in our data set have reasonably long (>3 oppositions) observational arcs.

Objects belonging to compositionally distinct TNO populations, which we term as interlopers, were rejected from our sample. Specifically, we excluded all known Haumea family members, as well as gray-colored OSSOS/Col-OSSOS objects orbiting within the Haumea family cluster (confirmed members have $41.6 < a < 44.2$ au, $0.08 < e < 0.19$, $24.2 < i < 29^\circ 1$; Brown et al. 2007; Ragozzine & Brown 2007; Schaller & Brown 2008; Snodgrass et al. 2010). We further rejected cold classical objects, which are known for exhibiting different colors (Tegler & Romanishin 2000; Pike et al. 2017) and albedos (Brucker et al. 2009; Vilenius et al. 2014) from those of the dynamically excited TNOs. Gulbis et al. (2010) and Petit et al. (2011) found that the orbital inclination distribution of cold classicals can be best reproduced by a Gaussian of width $2.0^{+0.6}_{-0.5}$ or $\sin(i)$ times a Gaussian of width $2.6$, respectively. Adopting an inclination cut of $i > 5^\circ$ ensures that contamination is minimal.
The same inclination cut was applied to all of the objects in our data set, independently of their dynamical classification, in order to perform a consistent analysis throughout all of the populations. We further verified that changing the inclination cut does not affect the conclusions of our analysis (Section 3). Centaurs with Tisserand parameter $T_J < 3$ were also excluded from our sample as their orbits are strongly coupled to Jupiter (e.g., Gladman et al. 2008) and thus have likely experienced larger orbital perturbations and changes in inclination (e.g., Volk & Malhotra 2013). The smallest ($H > 11$) centaurs and scattering disk objects appear to have altered colors compared to the rest of the centaur and scattering population (Jewitt 2015), so they were also excluded. No other populations in our data set contain objects with $H > 11$. Finally, we also rejected the two known objects with retrograde orbits because their origin population remains elusive (Gladman et al. 2009; Chen et al. 2016). Orbital elements and dynamical classifications for the resulting data set of 229 objects are reported in Appendix A.

## 2.2. Observations, Color Measurements, and Spectral Slopes

All of the new color measurements presented in this paper were acquired through the Col-OSSOS large program (PI: W. Fraser) on Gemini North between 2014 August and 2016 November. Col-OSSOS collects near-simultaneous *g*-, *r*-, and *J*-band photometry of a magnitude-limited ($r < 23.6$) subset of the OSSOS sample. Optical measurements were acquired with the Gemini Multi-Object Spectrograph (GMOS; Hook et al. 2004), and the *J*-band sequence was obtained with the Near InfraRed Imager and Spectrometer (NIRI; Hodapp et al. 2003). We observed using a *rgJgr* sequence to account for any brightness variation due to light-curve effects.

A complete description of the data reduction and photometric extraction for Col-OSSOS is provided in Schwamb et al. (2018). We summarize it briefly here. Processing of the Col-OSSOS data was performed with the Gemini-IRAF package. We first used the bias images acquired as part of the GMOS calibration plan to remove the bias offset from the science frames. Master sky flats were produced from a set of twilight flats and used to flatten the science frames. Photometric measurements were achieved using the TRIPPpy software package, which makes use of a pill-shaped aperture, built from the convolution product of a circular aperture with a line describing the direction and rate of motion of the target (Fraser et al. 2016). This approach minimizes the background contribution to the total flux in the photometric aperture.

Next, the Col-OSSOS data were calibrated to the Sloan Digital Sky Survey (SDSS) magnitude catalog (SDSS Collaboration et al. 2017) in cases where the field of view of the target was covered by SDSS; in other cases, we calibrated to the Pan-STARRS magnitude catalog (Magnier et al. 2013, 2016). We used in-frame catalog stars to determine the color conversion between the SDSS or Pan-STARRS and GMOS filter sets, using the relationships provided in Schwamb et al. (2018). We then fit the calibrated *g* and *r* magnitudes with a linear light-curve model assuming a constant (*g*-*r*) color across the observing sequence. The fit was weighted by the uncertainty of the individual photometric measurements calculated as the quadratic sum of the photometric uncertainty, the zero-point uncertainty, and the aperture correction uncertainty. The (*g*-*r*) color was then derived from our best fit, and the uncertainty on the color computed as the uncertainty on the fit. New colors obtained from our observations are reported in Table 1.

All of the color measurements from Col-OSSOS and the literature were converted to spectral slope (s), defined as the percentage increase in reflectance per $10^3$ Å change in wavelength normalized to 550 nm, using the Synphot tool in the STSDAS software package.<fnref n="13" /> Reported spectral slopes are the mean values from all the measurements of a target weighted by the inverse of the squared uncertainties. We consider only well-quantified measurements with a spectral slope uncertainty of $\Delta s < 12\%/(10^3$ Å), because larger error bars do not allow clear differentiation between the two color classes. The full list of measurements we consider is provided in Table 4 (Appendix A).

# 3. Red Dynamically Excited TNOs have Lower Orbital Inclinations

The whole data set is shown in Figure 1. The data suggest a relationship between color and inclination. In particular, there is a paucity of red TNOs with high orbital inclinations in our data set. In the following sections, we assess whether this trend is statistically robust. To do so, we divide our data set into two groups by color (Section 3.1) and test whether the two color groups have different inclination distributions (Section 3.2). An alternative approach would be to divide our data set into two inclination groups and compare the color distribution of these two groups. We take the first approach as our data set is clearly bimodal in color, as discussed in the next section. Effects of survey biases on our analysis are considered in Section 3.3.

## 3.1. Color Classification

Visually, the overall color distribution of all dynamically excited TNOs in our sample appears to be bimodal (Figure 2). The test of normality rejects a single Gaussian distribution at the 99.8% confidence level (CL). We therefore test if it can be decomposed into multiple normal distributions. We calculated the Bayesian information criterion (BIC) for a set of *N*-component Gaussian mixture models (GMM), and confirmed that the best BIC is found for a two-component model. The GMM decomposes our sample into two normal distributions, with respective means $s = 11.5\%/(10^3$ Å) and $s = 29.0\%/(10^3$ Å), and standard deviations $\sigma_s = 5.8\%/(10^3$ Å) and $\sigma_s = 8.4\%/(10^3$ Å). This confirms a previously known property of the Kuiper Belt, which is that the dynamically excited red TNOs exhibit a broader range of colors than the gray ones. The origin of this feature remains unclear (e.g., Brown et al. 2011), and we do not investigate it further, as this goes beyond the scope of this paper. The two normal distributions intersect at $s = 20.6\%/(10^3$ Å), corresponding to optical colors ($g-r) = 0.78$, ($V-R) = 0.56$, and (F606W-F814W) = -0.11. Therefore, we classify all objects with $s + \delta s < 20.6\%/(10^3$ Å) as gray and all objects with $s - \delta s > 20.6\%/(10^3$ Å) as red, where $\delta s$ is the uncertainty on the spectral slope. Objects at the boundary between the two color groups were not classified as they are ambiguous. Our classification is consistent with previous works that placed the bifurcation between gray and red TNOs at comparable spectral slope values (Table 2).

<footnote n="13">¹³ www.stsci.edu/institute/software_hardware/stsdas</footnote>
<table><caption>**Table 1**
New Col-OSSOS Color Measurements and Observing Circumstances</caption><thead><tr><th>Target<br>Designation</th><th>MPC<br>Designation</th><th>OSSOS<br>r' mag<sup>a</sup></th><th>H<sub>r</sub><sup>a</sup></th><th>Col-OSSOS<br>(g'-r')</th><th>Δ (au)</th><th>r (au)</th><th>α</th><th>Mean<br>MJD</th></tr></thead><tbody><tr><td>o3e01</td><td>2002GG166</td><td>21.50 ± 0.09</td><td>7.73</td><td>0.63 ± 0.01</td><td>21.69</td><td>20.70</td><td>0.52</td><td>57133.47513</td></tr><tr><td>o3e02</td><td>2013GH137</td><td>23.34 ± 0.14</td><td>8.32</td><td>0.77 ± 0.03</td><td>31.45</td><td>30.46</td><td>0.29</td><td>57519.33308</td></tr><tr><td>o3e05</td><td>2013GW136</td><td>22.69 ± 0.07</td><td>7.42</td><td>0.78 ± 0.01</td><td>32.62</td><td>31.65</td><td>0.46</td><td>57130.40010</td></tr><tr><td>o3e09</td><td>2013GY136</td><td>22.94 ± 0.05</td><td>7.32</td><td>0.55 ± 0.01</td><td>35.22</td><td>34.24</td><td>0.38</td><td>57132.42388</td></tr><tr><td>o3e23PD</td><td>2001FO185</td><td>23.37 ± 0.08</td><td>7.09</td><td>0.87 ± 0.06</td><td>42.03</td><td>41.03</td><td>0.20</td><td>57516.79552</td></tr><tr><td>o3e29</td><td>2013GO137</td><td>23.46 ± 0.08</td><td>7.09</td><td>0.78 ± 0.04</td><td>42.53</td><td>41.53</td><td>0.05</td><td>57139.41226</td></tr><tr><td>o3e39</td><td>2013GP136</td><td>23.07 ± 0.07</td><td>6.42</td><td>0.74 ± 0.02</td><td>44.75</td><td>43.75</td><td>0.13</td><td>57133.96175</td></tr><tr><td>o3e44</td><td>2013GG138</td><td>23.26 ± 0.09</td><td>6.34</td><td>1.18 ± 0.03</td><td>48.39</td><td>47.39</td><td>0.16</td><td>57135.40627</td></tr><tr><td>o3o09</td><td>2013JB65</td><td>23.23 ± 0.06</td><td>8.13</td><td>0.76 ± 0.02</td><td>31.98</td><td>31.00</td><td>0.47</td><td>57518.49578</td></tr><tr><td>o3o11(1)</td><td>2013JK64</td><td>22.94 ± 0.05</td><td>7.69</td><td>1.11 ± 0.01</td><td>33.43</td><td>32.48</td><td>0.61</td><td>57514.57625</td></tr><tr><td>o3o11(2)</td><td>2013JK64</td><td>22.94 ± 0.05</td><td>7.69</td><td>0.98 ± 0.01</td><td>33.28</td><td>32.47</td><td>1.03</td><td>57130.63127</td></tr><tr><td>o3o14</td><td>2013J064</td><td>23.54 ± 0.08</td><td>8.00</td><td>0.59 ± 0.02</td><td>35.74</td><td>34.92</td><td>0.93</td><td>57130.50551</td></tr><tr><td>o3o15</td><td>2013JD65</td><td>23.48 ± 0.07</td><td>7.90</td><td>0.83 ± 0.03</td><td>35.70</td><td>34.85</td><td>0.88</td><td>57132.55296</td></tr><tr><td>o3o18</td><td>2013JE64</td><td>23.56 ± 0.15</td><td>7.94</td><td>0.41 ± 0.13</td><td>36.65</td><td>35.66</td><td>0.36</td><td>57519.50134</td></tr><tr><td>o3o20PD</td><td>2007JF43</td><td>21.15 ± 0.02</td><td>5.27</td><td>1.01 ± 0.01</td><td>37.96</td><td>37.08</td><td>0.74</td><td>57135.50416</td></tr><tr><td>o3o21</td><td>2013JR65</td><td>23.51 ± 0.12</td><td>7.53</td><td>0.46 ± 0.03</td><td>38.83</td><td>37.86</td><td>0.47</td><td>57550.43376</td></tr><tr><td>o3o27</td><td>2013JJ65</td><td>23.40 ± 0.08</td><td>7.22</td><td>1.10 ± 0.03</td><td>41.43</td><td>40.62</td><td>0.82</td><td>57131.54637</td></tr><tr><td>o3o28</td><td>2013JN65</td><td>23.42 ± 0.21</td><td>7.23</td><td>0.57 ± 0.01</td><td>41.17</td><td>40.20</td><td>0.38</td><td>57514.47993</td></tr><tr><td>o3o29</td><td>2013JL64</td><td>23.26 ± 0.12</td><td>7.03</td><td>0.69 ± 0.03</td><td>42.36</td><td>41.39</td><td>0.40</td><td>57515.53610</td></tr><tr><td>o3o34</td><td>2013JH64</td><td>22.70 ± 0.04</td><td>5.60</td><td>0.70 ± 0.01</td><td>51.42</td><td>50.55</td><td>0.56</td><td>57133.53409</td></tr><tr><td>o5s32</td><td>2015RJ277</td><td>23.21 ± 0.04</td><td>7.12</td><td>0.63 ± 0.01</td><td>39.62</td><td>38.78</td><td>0.76</td><td>57696.25289</td></tr><tr><td>o5s45</td><td>2015RG277</td><td>23.15 ± 0.03</td><td>6.79</td><td>0.97 ± 0.02</td><td>42.51</td><td>41.62</td><td>0.60</td><td>57687.33189</td></tr><tr><td>o5s05<sup>b</sup></td><td>2015RV245</td><td>23.21 ± 0.04</td><td>10.10</td><td>0.61 ± 0.04</td><td>20.67</td><td>19.81</td><td>1.42</td><td>57696.36785</td></tr><tr><td>o5s06<sup>b</sup></td><td>2015RW245</td><td>22.90 ± 0.03</td><td>8.53</td><td>0.74 ± 0.02</td><td>26.51</td><td>25.57</td><td>0.74</td><td>57688.30839</td></tr><tr><td>o5s16PD<sup>b</sup></td><td>2004PB112</td><td>22.99 ± 0.03</td><td>7.39</td><td>0.77 ± 0.04</td><td>35.62</td><td>34.76</td><td>0.82</td><td>57695.35870</td></tr><tr><td>o5t04<sup>b</sup></td><td>2015RU245</td><td>22.99 ± 0.04</td><td>9.32</td><td>0.81 ± 0.01</td><td>22.53</td><td>21.65</td><td>1.17</td><td>57690.36639</td></tr></tbody></table>

**Notes.**

a OSSOS *r'* and $H_r$ mag are at discovery, at an earlier observing geometry than the Col-OSSOS color measurement.

b Object with previously published colors in Pike et al. (2017) and reprocessed here using the latest version of TRIPPy (Fraser et al. 2016).

Δ = heliocentric range; *r* = distance at discovery; α = phase angle.

## 3.2. Statistical Discrepancy between the Inclination Distributions of Gray and Red TNOs

### 3.2.1. Statistical Tests

Four main tests were used to assess the orbital inclination versus spectral slope distribution of our data set. The one-dimensional two-sample (2-s) Kolmogorov–Smirnov (KS; Chakravarti et al. 1967) and Anderson–Darling (AD; Stephens 1974) tests were used to determine whether the inclination distributions of the two TNO color groups defined in Section 3.1 come from a single parent distribution. These two tests are broadly similar and also complementary. The KS test is more sensitive to the center of the distributions, while the AD test gives more weight to the tails. The AD test is generally considered more powerful than the KS test. These two tests are nonparametric with no reference to a Gaussian. We therefore report 1 minus the *p*-value of the test (in percentage) as the level of confidence that the two distributions are different, rather than a number of *σ* deviation from a mean. We also used the student's T-test of means (Snedecor & Cochran 1989) to determine whether the mean inclinations of the gray and red TNOs are equal. Finally, we derived the Spearman coefficient (Zwillinger & Kokoska 2000) of our data set to search for a correlation between orbital inclination and spectral slope. Unlike the other three tests, the Spearman test does not require dividing our data set into two color groups. All of the results discussed in the following sections are summarized in Table 3. *p*-values were determined by computing each statistic for sets

of inclination/spectral slope pairs constructed by independently randomly sampling the observed distribution of each parameter (Efron & Tibshirani 1993). We sampled with replacement so the same value of inclination or spectral slope could be sampled several times. The reported CL then indicates the probability that a random sample produces a larger test statistic than the observed population.

### 3.2.2. Full Data Set

We first consider the full data set as a whole. All statistical tests indicate a strong (CL > 99.7%) discrepancy between the inclination distributions of the gray and red TNOs, as well as between their mean inclinations. The Spearman test indicates that the spectral slope anti-correlates with inclination at the same significance level (Table 3). The color versus inclination trend is particularly obvious when comparing the ratios of gray-to-red objects between low and high inclinations. From Figure 1, it appears that the inclination transition from where gray and red TNOs are roughly equally numbered to where red TNOs are almost nonexistent is located around *i* = 21°. We therefore chose to compare the ratios of red-to-gray objects above and below this value. Only five TNOs—i.e., 11% of the objects—can be unambiguously classified as red at *i* > 21° (Figure 1). In contrast, red objects account for 42% of the low-inclination (5° < *i* < 21°) population. To assess the likelihood of producing such an observation from a random population of objects, we bootstrap with replacement from the observed TNO samples of $N_g$ and $N_r$ objects, where $N_g$ and $N_r$ are the observed
<table><caption>**Table 3**

Confidence Levels for Gray and Red TNOs Having Different Inclination Distributions</caption><thead><tr><th></th><th>hc</th><th>res</th><th>cen+sca+det</th><th>all</th><th>all-cen</th></tr></thead><tbody><tr><td>2-sample KS test</td><td>98.1%</td><td>91.2%</td><td>97.8%</td><td>&gt;<strong>99.9%</strong></td><td>&gt;<strong>99.9%</strong></td></tr><tr><td>2-sample AD test</td><td>93.6%</td><td>98.5%</td><td>96.7%</td><td>&gt;<strong>99.9%</strong></td><td>&gt;<strong>99.9%</strong></td></tr><tr><td>T-test of means</td><td>99.5%</td><td>99.2%</td><td>97.9%</td><td>&gt;<strong>99.9%</strong></td><td>&gt;<strong>99.9%</strong></td></tr><tr><td>Spearman correlation (i &gt; 5°)</td><td>93.2%</td><td>95.0%</td><td>85.4%</td><td><strong>99.7%</strong></td><td>99.6%</td></tr><tr><td>Spearman correlation (5° &lt; i &lt; 21°)</td><td>60.3%</td><td>65.5%</td><td>82.2%</td><td>82.1%</td><td>82.5%</td></tr></tbody></table>

**Note.** hc = hot classical; res = resonant objects; cen = centaurs; sca = scattering disk objects; det = detached objects; all = full data set; all-cen = full data set without centaurs. The bold indicates a ≥99.7% confidence level.

agreement with previous works that reported an inclination versus color correlation for classical TNOs (Tegler & Romanishin 2000; Hainaut & Delsanti 2002; Trujillo & Brown 2002; Doressoudiram et al. 2005; Peixinho et al. 2008). These works however could have been contaminated by the presence in the sample of objects belonging to peculiar compositional groups, including cold classical objects. Here, we show that the observed discrepancy between gray and red hot classicals holds when interlopers are removed from this population.

A similar trend is found in the resonant group and in the group of centaurs/scattering/detached objects at a comparable significance level (Table 3). We also observe that out of the 102 resonant TNOs in our data set, where 49 are gray and 31 are red (the rest being unclassified), the 14 highest-inclination objects all belong to the gray class. When bootstrapping two random samples of 49 and 31 objects from the resonant population, we find that this occurs in less than 0.1% of the simulations, thus strongly supporting the conclusion that the gray and red resonant TNOs have different inclination distributions. Finding that resonant objects follow the same trend as other dynamical populations of TNOs is interesting. It implies that the initial orbital inclinations these objects had before being captured into orbital resonances with Neptune were not completely randomized by the capture process and the dynamical mechanisms, such as Kozai oscillation (e.g., Gomes et al. 2005), that take place in those resonances.

We stress that the group of centaurs/scattering/detached TNOs is strongly biased by the presence of centaurs in the sample. Excluding those objects completely removes the trend of gray and red TNOs having different orbital inclinations in that group. We note however that the scattering disk population does not contain any red objects at i > 21°, similar to what is observed in the resonant population and, to a lesser extent, the hot classical (1 interloper) and the centaur (2 interlopers) populations (Figure 3). Scattering disk objects might therefore follow the same trend as other TNO populations, but the current small number of scattering TNOs with well-measured colors prevents us from drawing any firm conclusion. Detached objects seem to constitute the only population that does not follow the trend, but their sample is also too small for conclusive analysis.

### 3.2.4. Summary of Statistical Tests

To summarize, we divided our data set into two color classes at 20.6%/(10³ Å) and found a highly significant difference between the orbital inclination distributions of the gray and red TNOs. This discrepancy cannot be detected at the CL > 99.7% significance level in the individual dynamical population of TNOs, but it does show up in each of these populations at lower statistical significance. This indicates that (1) this trend is common to every population of dynamically excited TNOs and (2) that our overall data set is not biased by one particular dynamical population. We

emphasize that the results presented here slightly vary depending on the adopted magnitude and inclination cuts used to define the analyzed sample, as well as the spectral slope value of 20.6%/(10³ Å) we adopt to differentiate between gray and red TNOs. We find however that a reasonable modification of those values does not affect the conclusions of our study in any significant way (see Appendix B).

### 3.3. Is the Observed Color-inclination Distribution an Artifact of the Discovery Surveys?

#### 3.3.1. Analytical Argument

The data set built for this work uses photometric measurements from a large variety of discovery surveys, each with its own detection biases and discovery efficiency. This raises the question as to whether these observing biases could be responsible for the observed inclination discrepancy between gray and red TNOs reported in this work. A discovery bias in favor of red TNOs comes from the combination of their higher optical albedos (Fraser et al. 2014; Lacerda et al. 2014) and redder colors with respect to gray objects. The observed ratio of red-to-gray objects depends on the filter used for the observations. Observations performed in broadband filters at longer wavelengths will discover a higher fraction of red TNOs than observations made at shorter wavelengths. This can introduce a color-inclination correlation if low-latitude discovery surveys (which are biased toward lower inclination TNOs) were performed at longer wavelengths than high-latitude surveys. To test whether this effect can be responsible for the observed distribution, let us consider the worst-case scenario, which would be if off-ecliptic TNO surveys were made in the V-band (550 nm), while on-ecliptic surveys were made in the R-band (660 nm).

The number of objects with albedos between a and a + da, between radii R and R + dR, and heliocentric distances between r and r + dr detected by the survey is

$$n(R, r, a) = A f(R) g(r) h(a) da dr dR, \quad (1)$$

where f(R) is the TNO size distribution, g(r) is the Kuiper Belt radial distribution, h(a) is the TNO albedo distribution, and A is a constant. Assuming a power-law distribution in R, f(R) ∝ R⁻q, and recasting the expression for R in terms of apparent magnitude m, it can be shown (Schwamb et al. 2018) that the number of objects with albedo a₀ ≤ a ≤ a₁ and magnitude m₀ ≤ m ≤ m₁ observed within the boundaries r₀ and r₁ of the Kuiper Belt is given by

$$N(m_0 < m < m_1) = C \int_{a_0}^{a_1} h(a) a^{5\alpha/2} da \\ \times \int_{r_0}^{r_1} g(r) r^{-5\alpha} \Delta^{-5\alpha} dr 10^{\alpha m}, \quad (2)$$
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Same as Figure 1 for the individual dynamical populations of the dynamically excited TNOs. Top left: hot classicals; middle: resonant objects; right: combined group of centaurs, scattering disk, and detached objects. Bottom: the combined group is divided into centaurs (left), scattering disk objects (middle), and detached objects (right). The absence of red TNOs on highly inclined orbits appears to be a common property of all dynamical populations. Centaurs are the only population for which different orbital inclinations of the two color classes have been shown so far (Tegler et al. 2016). Our data set reveals a similar trend to that of Tegler et al. (2016) for hot classicals, resonant objects, and centaurs, but with a lower statistical significance (see Table 3). The small number of scattering disk and detached objects in our sample does not permit us to statistically confirm the existence of a matching trend in those populations. However, we do observe a lack of red scattering disk objects on highly inclined orbits ($i > 21°$).</figcaption></figure>

where C is a constant, $\Delta$ is the geocentric distance to the object, and $\alpha = (q - 1)/5 \approx 1.0$ (Fraser et al. 2014) is the power-law slope of the apparent magnitude distribution. Here, we have assumed that the albedo, radial, and size distributions are independent. For simplicity, this expression does not reflect the apparent latitudinal and longitudinal variation of TNO density of

the resonant populations. We are only interested in the ratio of gray-to-red objects observed within a given survey. Assuming that within each orbital class, separate color populations share the same spatial distribution, the latitudinal and longitudinal variations will be reflected in both color populations equally. Thus, for our derivation such variations can be ignored. In the case of a
population composed of two color groups, where each group has its own albedo distribution, but similar size and radial distributions, and where the intrinsic number of objects in the two groups is given by $A_2 = \gamma A_1$, where $\gamma$ is a constant, the observed ratio of populations 1 and 2 is

$$ \frac{N_1}{N_2}(m_0 < m < m_1) = \gamma \frac{\int_{a_0}^{a_1} h_1(a) a^{\frac{5\alpha}{2}} da}{\int_{a_0}^{a_1} h_2(a) a^{\frac{5\alpha}{2}} da} \quad (3) $$

Here, we have assumed that the fraction of objects in the two groups is not a function of the discovery survey's observing depth. This is true as long as the size distribution can be modeled as a single power law ($f(R) \propto R^{-q}$), i.e., as long as the limiting magnitude of the surveys is not sensitive to the break (or divot) magnitude (see, e.g., Shankman et al. 2013, 2016; Fraser et al. 2014; Lawler et al. 2018b). A fraction of objects in our data set, mainly centaurs and scattering disk objects, are fainter than the magnitude break measured for these objects ($H = 7.7$ or 8.3; Lawler et al. 2018b). The presence of a break in the size distribution is explored numerically below in Section 3.3.2.

For simplicity, we consider the simple scenario where the two color groups each have a unique albedo, $a_1$ and $a_2$. Then, the albedo distribution is given by the Dirac delta function $h_n(a) = \delta(a_n - a)$, and Equation (3) becomes

$$ \frac{N_1}{N_2}(m_0 < m < m_1) = \gamma a_1^{\frac{5\alpha}{2}} a_2^{-\frac{5\alpha}{2}} \quad (4) $$

In the V-band (550 nm), gray and red TNOs have mean albedos of 6% and 12%, respectively (Fraser et al. 2014; Lacerda et al. 2014), implying a red-to-gray ratio of $N_{\text{red}}/N_{\text{gray}} = 5.66\gamma$. In this work, we derived mean spectral slopes of $11.4\%/(10^3 \AA)$ and $29.1\%/(10^3 \AA)$ for the gray and red color classes of TNOs, respectively. This implies albedos of 6.7% and 15.7% in the R-band (660 nm), and $N_{\text{red}}/N_{\text{gray}} = 8.41\gamma$. Therefore, in the V spectral band the observed red-to-gray ratio is 67% that in the R spectral band.

In our data set, we measure a red-to-gray ratio of 0.821 for $5^\circ < i < 10^\circ$ (23 red objects, 28 gray, and 17 unclassified). Assuming the worst-case scenario, where all on-ecliptic discovery surveys were performed in R and all off-ecliptic surveys were performed in V, the expected red-to-gray ratio at high inclination would be 0.550. Yet, the observed ratio for $i > 21^\circ$ is 0.119 (5 red objects, 42 gray, and 6 unclassified). If considering g- (480 nm) and r-band (620 nm) filters instead of V and R, then the observed ratio in g is 56% that in r, and the expected ratio in our data set would be 0.450 at $i > 21^\circ$. Choosing different inclination cuts for the two subsamples slightly varies these numbers but does not affect the conclusion. Therefore, it appears that even the observationally heavily biased scenario considered here cannot account for the observed relationship between color and inclination.

### 3.3.2. Survey Simulation

We further test the effects of discovery biases on the color versus inclination distribution of TNOs using the OSSOS survey simulator (Bannister et al. 2016, 2018; Lawler et al. 2018a). This simulator replicates the detection characteristics (pointings, detection and tracking efficiencies) of a given survey to determine which input model objects would be

detected by the survey. We use this simulator to model a heavily biased survey composed of a single on-ecliptic block performed in the r-band and two high-latitude blocks, with mean inclinations of $12^\circ$ and $22^\circ$, performed in g. All of the blocks are equally sensitive to the input objects in our simulation. Considering TNOs' mean ($g-r) = 0.78$ color, the observational depth was set to $m_r = 24.0$ for the on-ecliptic block, and $m_g = 24.75$ for the off-ecliptic ones. We use as input populations the Canadian-France Ecliptic Plane Survey (CFEPS) L7 model of the Kuiper Belt (Petit et al. 2011) and the improved verison of the model of Kaib et al. (2011, hereafter K11) by Shankman et al. (2016) of the scattering disk and centaur populations.

Each input object was cloned several times, allowing its orbital parameters a, e, and i to vary by 10% with respect to their nominal values, while randomizing its position angles ($\omega$, $\Omega$, and M). This cloning ensures that a sufficient number of simulated TNOs are randomly produced. We ran the simulator until a sample of 10,000 objects sampling the full parameter space was detected. Simulations were run independently for the on-ecliptic block and the high-latitude ones to ensure the same number of objects (5000) were detected by the blocks. Optical ($g-r$) colors were drawn from the double normal distribution measured in Section 3.1 and randomly attributed to the simulated objects. The input color ratio is irrelevant here; we are interested only in the relative color ratios between low-i and high-i objects in the output survey population. In order to determine the detectability of the simulated objects, an absolute magnitude distribution was assumed in the magnitude range of interest ($5 < H_r < 11$). We used two different distributions: a single-slope distribution with an exponent of $\alpha = 1.0$ and a broken-slope distribution with slope exponents of $\alpha_b = 1.0$ and $\alpha_f = 0.2$ for the bright and faint end of the distribution, respectively (Fraser et al. 2014). The break magnitude of the broken-slope distribution was set to $H_r = 7.7$ for the gray TNOs and $H_r = 6.9$ for the red ones in order to account for their respective mean albedos of 6% and 12%.

The population of detected objects in our simulated survey reveals a trend of increasing color with decreasing inclination (Figure 4), in agreement with our survey setup. We first consider the population drawn from the single-slope distribution. The red-to-gray ratio of the objects found in the on-ecliptic block is 54% that in the off-ecliptic blocks for the L7 model and 62% for the K11 model, in good agreement with the predicted value of 56% derived in Section 3.3.1. However, a fraction of the high-i TNOs in our simulation are detected by the low-latitude block. Because of this, the color discrepancy between the low-i and high-i objects in our simulation is weaker than predicted analytically (Section 3.3.1). The red-to-gray ratio of the high-inclination ($>21^\circ$) objects is 0.689 that of the low-inclination ($i < 10^\circ$) ones for the L7 model, and 0.854 for the K11 model. This strengthens our proposal that survey biases cannot account for the observed color versus inclination distribution. Finally, we find that the presence of a break in the magnitude distribution increases the red-to-gray ratio of the population detected in the high-latitude blocks relative to the on-ecliptic one; the on-ecliptic block ratio is 63% that in the off-ecliptic blocks for the L7 model, and 69% for the K11 model. This reinforces our conclusion that the observed color versus inclination correlation cannot be a result of the survey detection biases.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Inclination-dependent color ratio observed in our TNO sample (blue) compared to models of a heavily color-biased survey (Section 3.3.2; black). The survey model cannot reproduce the observed color ratio, implying it is intrinsic to the TNO population. Error bars on the TNO sample are derived from the number of color-unclassified TNOs (Section 3.1) in each inclination bin of 5°. The simulated survey observed a uniform color-ratio model on the ecliptic in *r* and at high latitudes in *g*. The indents in model color ratios at *i* ≈ 12°, 22° correspond to the mean inclinations of our simulated survey high-latitude sky blocks. The dotted vertical lines indicate theoretical extreme color-ratio cases observable for our simulation (Section 3.3.1).</figcaption></figure>

### 3.3.3. Real-life Surveys

We assess a subset of past surveys to further show that the correlation we report of color versus inclination distribution of the dynamically excited TNO populations is not the result of a discovery bias. Indeed, real surveys are far from being as biased as in the worst-case scenario we explored in Sections 3.3.1 and 3.3.2. We select the surveys with the most discoveries and an obtainable pointing history, including Chiang & Brown (1999), Gladman et al. (1998, 2001), Trujillo et al. (2001), Allen et al. (2002), Bernstein et al. (2004), Petit et al. (2006), Fuentes & Holman (2008), Fuentes et al. (2009), Fraser & Kavelaars (2009), Fraser et al. (2008, 2010), Alexandersen et al. (2016), the Deep Ecliptic Survey (DES; Millis et al. 2002; Elliot et al. 2005), CFEPS (Petit et al. 2011, 2017), and OSSOS (Bannister et al. 2016, 2018). This list is by no means exhaustive, but we ensure it covers well the range of ecliptic latitudes and of filters sampled by past surveys, and it includes the surveys that have reported the most detections of TNOs with $H_r > 5.0$ (the DES, CFEPS, and OSSOS). We omit the two surveys of Schwamb et al. (2010) and Rabinowitz et al. (2012), which were up to large ecliptic latitudes, in red-band and in broadband red+blue filters, respectively, as we could not retrieve their precise pointing history.

Figure 5 highlights that these surveys did not favor the detection of red TNOs on highly inclined orbits. On the contrary, most of the surveys performed far from the ecliptic plane used red-band filters, whereas the majority of surveys using blue-band and broadband blue+red filters were

conducted close to the ecliptic. We should therefore expect objects detected on highly inclined orbits to be preferentially red compared to the low-inclination ones, which is contrary to the correlation we report. We also note that only a small fraction of the known TNOs were detected in the high-latitude fields, in agreement with their narrow orbital inclination distribution. This implies that even if a discovery bias exists in our data set, it probably only has a marginal effect on its color-inclination distribution.

## 4. Discussion

Considering the commonly accepted idea that centaurs are dynamically evolved TNOs, it would appear surprising if the inclination dichotomy between gray and red centaurs found by Tegler et al. (2016) did not occur in the Kuiper Belt. Indeed, how would such a signal appear after the scattering of centaurs from the Kuiper Belt? Here, we have shown that this dichotomy also exists in the Kuiper Belt, and that it appears to be a common feature to all dynamical classes of the dynamically hot TNOs. This removes the conundrum of centaurs being the only population observed so far to exhibit that property.

Finding different orbital properties for the gray and red dynamically hot TNOs opens up interesting prospects for understanding the origin of the color diversity of small TNOs. So far, two main interpretations have been put forward to explain this diversity. The first hypothesis suggests that all TNOs were originally similar, and evolutionary—collisional and resurfacing—processes altered them differently. The second hypothesis is that the two TNO color classes are composed of intrinsically different objects, with distinct compositions, that probably formed at different locations in the protoplanetary disk.

Orbital inclination is directly linked to the strength of the collisional environment. As such, previous reports of a color versus inclination correlation in the classical belt (Tegler & Romanishin 2000; Hainaut & Delsanti 2002; Trujillo & Brown 2002; Doressoundiram et al. 2005; Peixinho et al. 2008) were interpreted to be the result of collisional resurfacing processes (Gil-Hutton 2002; Stern 2002; Trujillo & Brown 2002; Moroz et al. 2003; Delsanti et al. 2004; Santos-Sanz et al. 2009). In that scenario, we would expect a direct correlation between color and inclination because the collisional kinetic energy directly correlates with inclination. Yet, we find no such correlation in our data set for orbital inclinations $5° < i < 21°$, where both color classes are well sampled by surveys (Section 3.2.2). As such, the strong Spearman correlation found for the complete set of TNOs most likely comes from the absence of red TNOs above the $21°$ inclination. The lack of TNO color variation with rotation (Jewitt & Luu 2001) and the lack of a correlation between color and orbital eccentricity (Thébault & Doressoundiram 2003; Thébault 2003) also appear to contradict a collisional origin for the inclination versus color relationship. Numerical simulations reveal that the kinetic energy TNOs receive from collisions always correlates better with eccentricity than inclination. Applying the same statistical tests as presented in Section 3.2 to the eccentricity versus color distribution in our data set reveals no correlation between these parameters.

It therefore appears that the dynamically excited TNOs formed in two intrinsically different populations, each with its
<figure><img src="image_6.png" /><figcaption>**Figure 5.** On-sky pointing coordinates of the TNO discovery surveys listed in Section 3.3.3. We show the observing filter in which these surveys were performed, and schematically indicate their number of discoveries. Only the survey fields with positive detections are displayed for clarity. The dotted black curve and the blue curve indicate the location of the ecliptic and the galactic planes, respectively. Almost all of the surveys performed with blue-band and broadband blue+red filters were conducted close to the ecliptic plane. Most high-latitude observations come from the high-latitude component of CFEPS (Petit et al. 2017), which was performed with a red-band filter. These past observations imply that trans-Neptunian objects detected on highly inclined orbits should be preferentially red compared to the low-inclination ones. This is contrary to the correlation we report.</figcaption></figure>

own orbital inclination and color distributions. Dynamical simulations have shown that the high-inclination TNOs could have formed closer to the Sun and been emplaced at their current location following a phase of planetary migration (Gomes 2003). Our study therefore suggests that gray TNOs accreted closer to the Sun compared to the red objects. During migration, their orbits were more excited by Neptune, which widened their orbital inclination distribution. The presence of gray class interlopers in the cold classical region, the so-called blue binaries (Fraser et al. 2017), may however complicate this picture. The fact that these interlopers are all gray suggests that the gray class originally bordered the inner edge of the cold classical region (Schwamb et al. 2018). Future detailed modeling of planetary migration in a color-varying planetesimal disk will provide predictions which can be tested against our observed color-inclination distribution.

Different formation locations and emplacement efficiency as the origin of the two TNO color groups also appears compatible with the observation that gray TNOs are proportionally more numerous in the dynamically excited TNO populations. In our data set, gray TNOs represent 61%, 64%, 83%, and 87% of the total number of objects in the resonant, centaur, scattering, and detached populations, respectively, whereas they represent only 56% of the hot classicals. We however stress that this observation must be considered with extreme caution as it relies on a heavily biased data set built from various surveys with different detection biases. The Col-OSSOS survey will, upon completion, provide a complete fully debiased sample of TNO colors from which the intrinsic ratio of gray-to-red objects in each dynamical population will be derived (Schwamb et al. 2018).

**5. Conclusion**

We report new (*g*-*r*) colors for 25 small (*H*<sub>mag</sub> > 5) dynamically excited TNOs observed in the Col-OSSOS survey. Combined with previously published color measurements from Col-OSSOS and other surveys, we consider the colors of a set of 229 dynamically excited TNOs and centaurs. We show that, when dividing our data set into two color classes at a spectral

slope of *s* = 20.6%/((10³ Å)), red objects have a lower orbital inclination distribution than their gray counterparts. This is true whether or not centaurs are included in our sample. This trend appears to be common to every dynamical population in the Kuiper Belt. The scattering and detached populations are too sparsely sampled to confirm a matching trend in those populations. Even in the worst-case scenario, observing biases in the discovery surveys cannot account for this correlation: it is intrinsic to the underlying TNO and centaur populations. Considering that TNOs are the precursors of centaurs, and that their inclinations are roughly preserved as they become centaurs (Volk & Malhotra 2013), our finding solves the conundrum of centaurs being the only outer solar system population identified so far to exhibit this property (Tegler et al. 2016). Finally, the observed inclination difference between gray and red TNOs tends to favor the hypothesis of different formation locations for these objects, and a compositional gradient in the original protodisk of planetesimals, rather than different evolutionary pathways.

The authors acknowledge the sacred nature of Maunakea, and appreciate the opportunity to observe from the mountain. This work is based on observations from the Large and Long Program GN-LP-1 (2014B through 2016B), obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). We thank the staff at Gemini North for their dedicated support of the Col-OSSOS program. Data was processed using the Gemini-IRAF package. STSDAS and PyRAF are products of the Space Telescope Science Institute, which is operated by AURA for NASA. This research used the facilities of the Canadian Astronomy Data Centre operated by the National Research Council of Canada with the support of the Canadian Space Agency. The authors thank C. Shankman for sharing his
Python implementation of the OSSOS survey simulator, and S. Lawler and N. Kaib for providing their model of the scattering disk and centaur populations. M.T.B. appreciates support from UK STFC grant ST/L000709/1. M.E.S. was supported by Gemini Observatory. K.V. acknowledges support from NASA grants NNX14AG93G and NNX15AH59G. N.P. acknowledges funding from the Portuguese FCT—Foundation for Science and Technology (ref: SFRH/BGCT/113686/2015). CITEUC is funded by National Funds through FCT—Foundation for Science and Technology (project: UID/ Multi/00611/2013) and FEDER—European Regional Development Fund through COMPETE 2020—Operational Programme Competitiveness and Internationalisation (project: POCI-01-0145-FEDER-006922). W.C.F. acknowledges support from UK Science and Technology Facilities Council grant ST/P0003094/1.

Facilities: Gemini:Gillett (GMOS-N, NIRI).

Software: Astropy (Astropy Collaboration et al. 2013), IRAF (Tody 1986), Matplotlib (Hunter 2007), NumPy (van der Walt et al. 2011), OSSOS survey simulator (Lawler et al. 2018a), pyraf, Scikit-learn (Pedregosa et al. 2011), SciPy (van der Walt et al. 2011), Sextractor (Bertin & Arnouts 1996), STSDAS, TRIPPy (Fraser et al. 2016).

## Appendix A
### Full Data Set

The complete set of 229 outer solar system objects on dynamically-excited orbits considered in this work is provided in Table 4, together with the orbital and physical properties of these objects.

<table><caption>**Table 4**
Orbit, Magnitude, and Measured Spectral Slope of the Outer Solar System Objects Studied in This Work</caption><thead><tr><th>MPC Number</th><th>Target Designation</th><th>Orbit Classification</th><th>H<sub>r</sub></th><th>a (au)</th><th>e</th><th>inc (°)</th><th>Spectral slope (%/10<sup>3</sup> Å)</th><th>References</th></tr></thead><tbody><tr><td colspan="9" style="text-align:center; font-weight:bold;">Hot Classical</td></tr><tr><td rowspan="5">181708</td><td>1993FW</td><td>cla (i)</td><td>6.85</td><td>43.72</td><td>0.05</td><td>7.75</td><td>22.23 ± 6.16</td><td>2</td></tr><tr><td>1996RQ20</td><td>cla</td><td>7.17</td><td>43.86</td><td>0.10</td><td>31.67</td><td>20.64 ± 9.66</td><td>2</td></tr><tr><td>1996TS66</td><td>cla</td><td>6.23</td><td>43.89</td><td>0.13</td><td>7.35</td><td>31.62 ± 10.39</td><td>2</td></tr><tr><td>1997CV29</td><td>cla (i)</td><td>7.45</td><td>42.18</td><td>0.05</td><td>8.03</td><td>31.50 ± 2.50</td><td>2</td></tr><tr><td>1997QH4</td><td>cla</td><td>7.28</td><td>42.70</td><td>0.08</td><td>13.19</td><td>29.57 ± 9.68</td><td>2</td></tr><tr><td rowspan="2">385191</td><td>1997RT5</td><td>cla</td><td>7.21</td><td>41.31</td><td>0.03</td><td>12.73</td><td>14.20 ± 8.10</td><td>2</td></tr><tr><td>1998FS144</td><td>cla</td><td>6.99</td><td>41.77</td><td>0.02</td><td>9.87</td><td>20.68 ± 5.84</td><td>2</td></tr><tr><td rowspan="4">181855</td><td>1998WT31</td><td>cla</td><td>7.68</td><td>45.98</td><td>0.18</td><td>28.68</td><td>9.38 ± 10.16</td><td>2</td></tr><tr><td>1999CB119</td><td>cla</td><td>7.04</td><td>46.91</td><td>0.13</td><td>8.75</td><td>34.40 ± 7.22</td><td>2</td></tr><tr><td>1999CL119</td><td>cla</td><td>6.00</td><td>46.87</td><td>0.01</td><td>23.29</td><td>16.71 ± 1.82</td><td>2,3</td></tr><tr><td>1999CL158</td><td>cla (i)</td><td>6.88</td><td>41.44</td><td>0.21</td><td>10.03</td><td>5.16 ± 3.29</td><td>2</td></tr><tr><td>79983</td><td>1999DF9</td><td>cla</td><td>6.10</td><td>46.38</td><td>0.14</td><td>9.82</td><td>36.41 ± 5.46</td><td>2</td></tr><tr><td>118379</td><td>1999HC12</td><td>cla</td><td>8.08</td><td>45.30</td><td>0.23</td><td>15.35</td><td>3.16 ± 7.90</td><td>2</td></tr><tr><td rowspan="2">38083</td><td>1999HX11 (Rhadamanthus)</td><td>cla (i)</td><td>7.67</td><td>38.98</td><td>0.15</td><td>12.75</td><td>9.34 ± 4.75</td><td>2</td></tr><tr><td>1999RY214</td><td>cla</td><td>7.26</td><td>45.13</td><td>0.18</td><td>13.69</td><td>19.47 ± 6.64</td><td>2</td></tr><tr><td rowspan="10">86177</td><td>1999RY215</td><td>cla</td><td>7.02</td><td>45.26</td><td>0.24</td><td>22.18</td><td>2.65 ± 1.38</td><td>2,3</td></tr><tr><td>1999XY143</td><td>cla</td><td>6.20</td><td>43.00</td><td>0.08</td><td>7.17</td><td>27.53 ± 6.43</td><td>2</td></tr><tr><td>2000CJ105</td><td>cla</td><td>5.99</td><td>44.16</td><td>0.11</td><td>11.59</td><td>36.53 ± 7.18</td><td>2</td></tr><tr><td>2000CO105</td><td>cla</td><td>5.89</td><td>47.05</td><td>0.15</td><td>19.27</td><td>21.27 ± 8.56</td><td>2</td></tr><tr><td>2000CP104</td><td>cla</td><td>6.86</td><td>44.25</td><td>0.10</td><td>9.48</td><td>21.43 ± 5.53</td><td>2</td></tr><tr><td>2000KK4</td><td>cla</td><td>6.27</td><td>41.26</td><td>0.08</td><td>19.13</td><td>26.81 ± 7.47</td><td>2</td></tr><tr><td>2000OH67</td><td>cla</td><td>7.49</td><td>44.16</td><td>0.02</td><td>5.63</td><td>24.60 ± 5.30</td><td>3</td></tr><tr><td>2001CZ31</td><td>cla</td><td>5.88</td><td>45.19</td><td>0.12</td><td>10.22</td><td>4.00 ± 0.70</td><td>3</td></tr><tr><td>2001FO185 (o3e23PD)</td><td>cla</td><td>7.09</td><td>46.45</td><td>0.12</td><td>10.64</td><td>26.50 ± 3.70</td><td>1</td></tr><tr><td rowspan="10">469361</td><td>2001HY65</td><td>cla</td><td>6.35</td><td>43.16</td><td>0.12</td><td>17.14</td><td>20.10 ± 1.55</td><td>2,3</td></tr><tr><td>2001KA77</td><td>cla</td><td>5.41</td><td>47.41</td><td>0.09</td><td>11.92</td><td>35.09 ± 1.29</td><td>2,3</td></tr><tr><td>2001PK47</td><td>cla</td><td>7.49</td><td>39.76</td><td>0.07</td><td>8.73</td><td>8.70 ± 3.10</td><td>3</td></tr><tr><td>2001QC298</td><td>cla</td><td>6.49</td><td>46.27</td><td>0.12</td><td>30.57</td><td>-1.25 ± 1.45</td><td>2,3</td></tr><tr><td>2001QR297</td><td>cla</td><td>6.77</td><td>44.40</td><td>0.03</td><td>5.14</td><td>24.60 ± 4.10</td><td>3</td></tr><tr><td>2001RY143 (o4h48)</td><td>cla</td><td>6.80</td><td>42.08</td><td>0.15</td><td>6.91</td><td>29.50 ± 2.10</td><td>1</td></tr><tr><td>2002PD155</td><td>cla</td><td>7.42</td><td>43.23</td><td>0.01</td><td>5.77</td><td>18.00 ± 5.90</td><td>3</td></tr><tr><td>2003LD9</td><td>cla</td><td>6.92</td><td>47.26</td><td>0.17</td><td>6.97</td><td>31.30 ± 3.30</td><td>3</td></tr><tr><td>2003QQ91</td><td>cla</td><td>7.87</td><td>38.77</td><td>0.07</td><td>5.41</td><td>15.30 ± 9.00</td><td>2</td></tr><tr><td rowspan="3">307616</td><td>2003QW90</td><td>cla</td><td>5.02</td><td>43.77</td><td>0.08</td><td>10.36</td><td>32.36 ± 5.55</td><td>2</td></tr><tr><td>2003UY413</td><td>cla</td><td>7.01</td><td>47.25</td><td>0.21</td><td>20.74</td><td>25.50 ± 9.70</td><td>2</td></tr><tr><td>2004PA112</td><td>cla</td><td>7.54</td><td>38.92</td><td>0.11</td><td>32.93</td><td>6.00 ± 1.90</td><td>3</td></tr><tr><td rowspan="3">504847</td><td>2010RE188 (o3l18)</td><td>cla</td><td>6.18</td><td>46.01</td><td>0.15</td><td>6.75</td><td>8.50 ± 0.60</td><td>1</td></tr><tr><td>2013GG138 (o3e44)</td><td>cla</td><td>6.34</td><td>47.46</td><td>0.03</td><td>24.61</td><td>44.60 ± 1.70</td><td>1</td></tr><tr><td>2013GO137 (o3e29)</td><td>cla</td><td>7.09</td><td>41.42</td><td>0.09</td><td>29.25</td><td>20.70 ± 2.30</td><td>1</td></tr><tr><td rowspan="2">500886</td><td>2013JN65 (o3o28)</td><td>cla (i)</td><td>7.23</td><td>40.67</td><td>0.01</td><td>19.64</td><td>7.80 ± 0.70</td><td>1</td></tr><tr><td>2013JR65 (o3o21)</td><td>cla</td><td>7.53</td><td>46.20</td><td>0.19</td><td>11.71</td><td>0.40 ± 2.10</td><td>1</td></tr><tr><td rowspan="2">505448</td><td>2013SA100 (o3l79)</td><td>cla</td><td>5.75</td><td>46.30</td><td>0.17</td><td>8.48</td><td>13.10 ± 0.48</td><td>1</td></tr><tr><td>2013SZ99 (o3l15)</td><td>cla</td><td>7.52</td><td>38.28</td><td>0.02</td><td>19.84</td><td>14.70 ± 1.20</td><td>1</td></tr><tr><td></td><td>2014UH225 (o4h29)</td><td>cla</td><td>7.30</td><td>38.64</td><td>0.04</td><td>29.53</td><td>7.80 ± 1.30</td><td>1</td></tr></tbody></table>
<table><caption>**Table 4**
(Continued)</caption><thead><tr><th>MPC<br>Number</th><th>Target<br>Designation</th><th>Orbit<br>Classification</th><th>H<sub>r</sub></th><th>a (au)</th><th>e</th><th>inc (°)</th><th>Spectral slope<br>(%/10<sup>3</sup> Å)</th><th>References</th></tr></thead><tbody><tr><td>511553</td><td>2014UK225 (o4h19)</td><td>cla (i)</td><td>7.43</td><td>43.53</td><td>0.13</td><td>10.69</td><td>31.40 ± 1.30</td><td>1</td></tr><tr><td>511554</td><td>2014UL225 (o4h20)</td><td>cla</td><td>7.24</td><td>46.34</td><td>0.20</td><td>7.95</td><td>6.40 ± 2.10</td><td>1</td></tr><tr><td></td><td>2015RG277 (o5s45)</td><td>cla</td><td>6.79</td><td>42.96</td><td>0.01</td><td>12.09</td><td>32.80 ± 1.20</td><td>1</td></tr><tr><td></td><td>2015RJ277 (o5s32)</td><td>cla</td><td>7.12</td><td>46.69</td><td>0.20</td><td>5.52</td><td>11.70 ± 0.70</td><td>1</td></tr><tr><td colspan="9">Resonant Objects</td></tr><tr><td></td><td>2006RJ103</td><td>1:1</td><td>7.67</td><td>30.04</td><td>0.03</td><td>8.16</td><td>6.16 ± 2.69</td><td>2</td></tr><tr><td></td><td>2007VL305</td><td>1:1</td><td>8.42</td><td>30.03</td><td>0.06</td><td>28.12</td><td>15.77 ± 6.10</td><td>2</td></tr><tr><td></td><td>1998UU43</td><td>4:3</td><td>7.24</td><td>36.46</td><td>0.13</td><td>9.58</td><td>16.40 ± 7.23</td><td>2</td></tr><tr><td>143685</td><td>2003SS317</td><td>4:3</td><td>7.88</td><td>36.46</td><td>0.24</td><td>5.91</td><td>31.31 ± 3.89</td><td>2</td></tr><tr><td></td><td>2004TX357</td><td>4:3</td><td>8.78</td><td>36.58</td><td>0.21</td><td>16.26</td><td>11.65 ± 2.77</td><td>2</td></tr><tr><td></td><td>2005ER318</td><td>4:3</td><td>7.92</td><td>36.52</td><td>0.16</td><td>10.42</td><td>11.09 ± 2.55</td><td>2</td></tr><tr><td></td><td>2014UD229 (o4h13)</td><td>4:3</td><td>8.18</td><td>36.40</td><td>0.15</td><td>6.85</td><td>16.40 ± 0.80</td><td>1</td></tr><tr><td></td><td>2014UX228 (o4h18)</td><td>4:3</td><td>7.35</td><td>36.36</td><td>0.17</td><td>20.66</td><td>7.00 ± 1.50</td><td>1</td></tr><tr><td>15789</td><td>1993SC</td><td>3:2</td><td>7.02</td><td>39.52</td><td>0.19</td><td>5.16</td><td>39.84 ± 8.77</td><td>2</td></tr><tr><td>20108</td><td>1995QZ9</td><td>3:2</td><td>8.14</td><td>39.45</td><td>0.15</td><td>19.56</td><td>16.40 ± 5.60</td><td>2</td></tr><tr><td></td><td>1996RR20</td><td>3:2</td><td>6.92</td><td>39.53</td><td>0.17</td><td>5.31</td><td>37.28 ± 11.41</td><td>2</td></tr><tr><td>15875</td><td>1996TP66</td><td>3:2</td><td>7.23</td><td>39.38</td><td>0.33</td><td>5.69</td><td>32.58 ± 8.05</td><td>2</td></tr><tr><td>118228</td><td>1996TQ66</td><td>3:2</td><td>7.54</td><td>39.41</td><td>0.12</td><td>14.66</td><td>36.47 ± 11.89</td><td>2</td></tr><tr><td>24952</td><td>1997QJ4</td><td>3:2</td><td>7.61</td><td>39.36</td><td>0.23</td><td>16.57</td><td>6.35 ± 6.26</td><td>2</td></tr><tr><td>91133</td><td>1998HK151</td><td>3:2</td><td>7.17</td><td>39.36</td><td>0.23</td><td>5.94</td><td>7.81 ± 3.90</td><td>2</td></tr><tr><td>91205</td><td>1998US43</td><td>3:2</td><td>8.06</td><td>39.26</td><td>0.13</td><td>10.61</td><td>1.32 ± 4.33</td><td>2</td></tr><tr><td></td><td>1998WS31</td><td>3:2</td><td>8.21</td><td>39.36</td><td>0.20</td><td>6.74</td><td>18.76 ± 9.68</td><td>2</td></tr><tr><td>69990</td><td>1998WU31</td><td>3:2</td><td>8.22</td><td>39.23</td><td>0.19</td><td>6.58</td><td>11.54 ± 5.32</td><td>2</td></tr><tr><td></td><td>1998WV31</td><td>3:2</td><td>7.88</td><td>39.25</td><td>0.27</td><td>5.72</td><td>17.63 ± 8.15</td><td>2</td></tr><tr><td>69986</td><td>1998WW24</td><td>3:2</td><td>8.30</td><td>39.43</td><td>0.23</td><td>13.93</td><td>17.08 ± 3.99</td><td>2,3</td></tr><tr><td></td><td>1998WZ31</td><td>3:2</td><td>8.26</td><td>39.50</td><td>0.17</td><td>14.60</td><td>3.16 ± 6.60</td><td>2</td></tr><tr><td>47171</td><td>1999TC36 (Lempo)</td><td>3:2</td><td>5.15</td><td>39.41</td><td>0.23</td><td>8.41</td><td>34.21 ± 3.43</td><td>2</td></tr><tr><td></td><td>1999TR11</td><td>3:2</td><td>8.37</td><td>39.40</td><td>0.24</td><td>17.14</td><td>40.30 ± 11.82</td><td>2</td></tr><tr><td></td><td>2000FV53</td><td>3:2</td><td>7.81</td><td>39.18</td><td>0.16</td><td>17.34</td><td>2.17 ± 2.05</td><td>2,3</td></tr><tr><td>47932</td><td>2000GN171</td><td>3:2</td><td>5.92</td><td>39.35</td><td>0.28</td><td>10.81</td><td>25.16 ± 3.04</td><td>2,4</td></tr><tr><td></td><td>2000YH2</td><td>3:2 (i)</td><td>8.05</td><td>39.23</td><td>0.30</td><td>12.91</td><td>6.00 ± 1.50</td><td>3</td></tr><tr><td>131318</td><td>2001FL194</td><td>3:2</td><td>7.88</td><td>39.31</td><td>0.17</td><td>13.70</td><td>6.00 ± 1.90</td><td>3</td></tr><tr><td>469362</td><td>2001KB77</td><td>3:2</td><td>7.52</td><td>39.57</td><td>0.28</td><td>17.52</td><td>14.03 ± 5.02</td><td>2,4</td></tr><tr><td>469372</td><td>2001QF298</td><td>3:2</td><td>5.25</td><td>39.38</td><td>0.11</td><td>22.35</td><td>3.82 ± 3.30</td><td>2,4</td></tr><tr><td>139775</td><td>2001QG298</td><td>3:2</td><td>7.31</td><td>39.35</td><td>0.19</td><td>6.49</td><td>24.60 ± 1.10</td><td>3</td></tr><tr><td>126155</td><td>2001YJ140</td><td>3:2</td><td>7.51</td><td>39.43</td><td>0.29</td><td>5.98</td><td>10.40 ± 0.80</td><td>3</td></tr><tr><td>55638</td><td>2002VE95</td><td>3:2</td><td>5.41</td><td>39.29</td><td>0.29</td><td>16.33</td><td>39.36 ± 2.41</td><td>2,4</td></tr><tr><td>84719</td><td>2002VR128</td><td>3:2</td><td>5.19</td><td>39.42</td><td>0.26</td><td>14.01</td><td>26.03 ± 1.66</td><td>2,4</td></tr><tr><td></td><td>2003SR317 (o3l13PD)</td><td>3:2</td><td>7.98</td><td>39.43</td><td>0.17</td><td>8.35</td><td>12.00 ± 1.10</td><td>1</td></tr><tr><td></td><td>2003TH58</td><td>3:2</td><td>7.17</td><td>39.36</td><td>0.09</td><td>27.93</td><td>6.33 ± 4.21</td><td>2</td></tr><tr><td>120216</td><td>2004EW95</td><td>3:2</td><td>6.46</td><td>39.38</td><td>0.32</td><td>29.28</td><td>3.07 ± 0.62</td><td>2,3,4</td></tr><tr><td>469708</td><td>2005GE187</td><td>3:2</td><td>7.51</td><td>39.34</td><td>0.33</td><td>18.25</td><td>31.37 ± 1.97</td><td>2,3</td></tr><tr><td></td><td>2005TV189</td><td>3:2</td><td>7.72</td><td>39.41</td><td>0.19</td><td>34.39</td><td>9.80 ± 1.60</td><td>3</td></tr><tr><td>444745</td><td>2007JF43 (o3o20PD)</td><td>3:2</td><td>5.27</td><td>39.36</td><td>0.18</td><td>15.08</td><td>35.10 ± 0.60</td><td>1</td></tr><tr><td></td><td>2010TJ182 (o4h07)</td><td>3:2</td><td>7.68</td><td>39.65</td><td>0.28</td><td>9.50</td><td>9.40 ± 0.90</td><td>1</td></tr><tr><td></td><td>2013GH137 (o3e02)</td><td>3:2</td><td>8.32</td><td>39.44</td><td>0.23</td><td>13.47</td><td>20.00 ± 2.00</td><td>1</td></tr><tr><td></td><td>2013JB65 (o3o09)</td><td>3:2</td><td>8.13</td><td>39.40</td><td>0.19</td><td>24.90</td><td>19.90 ± 1.30</td><td>1</td></tr><tr><td></td><td>2013JD65 (o3o15)</td><td>3:2</td><td>7.90</td><td>39.37</td><td>0.09</td><td>13.02</td><td>23.90 ± 2.00</td><td>1</td></tr><tr><td>500883</td><td>2013JJ65 (o3o27)</td><td>3:2</td><td>7.22</td><td>39.37</td><td>0.26</td><td>19.82</td><td>40.30 ± 1.80</td><td>1</td></tr><tr><td></td><td>2014UO229 (o4h11)</td><td>3:2</td><td>8.25</td><td>39.45</td><td>0.16</td><td>10.09</td><td>16.06 ± 0.57</td><td>1</td></tr><tr><td></td><td>2014UV228 (o4h09)</td><td>3:2</td><td>8.49</td><td>39.49</td><td>0.23</td><td>10.13</td><td>13.00 ± 1.40</td><td>1</td></tr><tr><td></td><td>2014UX229 (o4h05)</td><td>3:2</td><td>8.04</td><td>39.64</td><td>0.34</td><td>15.97</td><td>12.20 ± 1.00</td><td>1</td></tr><tr><td></td><td>1999CX131</td><td>5:3</td><td>7.17</td><td>42.24</td><td>0.23</td><td>9.76</td><td>14.40 ± 6.08</td><td>2</td></tr><tr><td>126154</td><td>2001YH140</td><td>5:3</td><td>5.63</td><td>42.33</td><td>0.14</td><td>11.08</td><td>14.10 ± 1.94</td><td>2</td></tr><tr><td>149349</td><td>2002VA131</td><td>5:3</td><td>6.72</td><td>42.27</td><td>0.24</td><td>7.07</td><td>38.16 ± 4.07</td><td>2</td></tr><tr><td></td><td>2005SE278</td><td>5:3</td><td>7.06</td><td>42.31</td><td>0.11</td><td>6.89</td><td>36.22 ± 5.72</td><td>2</td></tr><tr><td>434709</td><td>2006CJ69</td><td>5:3</td><td>7.62</td><td>42.18</td><td>0.23</td><td>17.92</td><td>37.17 ± 2.85</td><td>2</td></tr><tr><td>469306</td><td>1999CD158</td><td>7:4</td><td>5.10</td><td>43.66</td><td>0.14</td><td>25.48</td><td>18.33 ± 8.51</td><td>2</td></tr><tr><td>118378</td><td>1999HT11</td><td>7:4 (i)</td><td>7.14</td><td>43.74</td><td>0.11</td><td>5.05</td><td>32.37 ± 4.44</td><td>2</td></tr><tr><td>60620</td><td>2000FD8</td><td>7:4</td><td>6.71</td><td>43.75</td><td>0.22</td><td>19.51</td><td>28.97 ± 2.29</td><td>2,3</td></tr><tr><td>182222</td><td>2000YU1</td><td>7:4 (i)</td><td>6.68</td><td>43.59</td><td>0.10</td><td>6.11</td><td>18.98 ± 3.75</td><td>2</td></tr><tr><td></td><td>2003QX91</td><td>7:4</td><td>8.83</td><td>43.66</td><td>0.25</td><td>27.68</td><td>11.50 ± 7.30</td><td>3</td></tr></tbody></table>
<table><caption>**Table 4**
(Continued)</caption><thead><tr><th>MPC<br>Number</th><th>Target<br>Designation</th><th>Orbit<br>Classification</th><th>H<sub>r</sub></th><th>a (au)</th><th>e</th><th>inc (°)</th><th>Spectral slope<br>(%/10<sup>3</sup> Å)</th><th>References</th></tr></thead><tbody><tr><td></td><td>2004OQ15</td><td>7:4</td><td>7.22</td><td>43.65</td><td>0.12</td><td>9.73</td><td>6.71 ± 3.62</td><td>2</td></tr><tr><td></td><td>2005SF278</td><td>7:4</td><td>6.56</td><td>43.73</td><td>0.19</td><td>13.34</td><td>33.59 ± 4.89</td><td>2</td></tr><tr><td>437915</td><td>2002GD32</td><td>9:5</td><td>6.16</td><td>44.48</td><td>0.14</td><td>6.59</td><td>34.38 ± 5.19</td><td>2</td></tr><tr><td>511555</td><td>2014UM225 (o4h31)</td><td>9:5</td><td>7.21</td><td>44.48</td><td>0.10</td><td>18.30</td><td>22.30 ± 1.90</td><td>1</td></tr><tr><td></td><td>1997SZ10</td><td>2:1</td><td>8.37</td><td>47.89</td><td>0.36</td><td>11.80</td><td>28.36 ± 9.51</td><td>2</td></tr><tr><td>26308</td><td>1998SM165</td><td>2:1</td><td>5.96</td><td>47.66</td><td>0.37</td><td>13.50</td><td>25.15 ± 2.21</td><td>2,3</td></tr><tr><td>137295</td><td>1999RB216</td><td>2:1</td><td>7.12</td><td>47.55</td><td>0.29</td><td>12.69</td><td>17.31 ± 1.43</td><td>2,3</td></tr><tr><td>130391</td><td>2000JG81</td><td>2:1</td><td>8.02</td><td>47.42</td><td>0.28</td><td>23.46</td><td>11.16 ± 3.14</td><td>2,3</td></tr><tr><td></td><td>2004TV357</td><td>2:1</td><td>6.76</td><td>47.72</td><td>0.28</td><td>9.76</td><td>-0.62 ± 0.80</td><td>2,3</td></tr><tr><td></td><td>2005CA79</td><td>2:1</td><td>5.24</td><td>47.78</td><td>0.22</td><td>11.67</td><td>15.60 ± 0.77</td><td>2,3</td></tr><tr><td>308379</td><td>2005RS43</td><td>2:1</td><td>5.18</td><td>47.80</td><td>0.20</td><td>10.01</td><td>7.89 ± 0.67</td><td>2,3</td></tr><tr><td>470083</td><td>2006SG369</td><td>2:1</td><td>7.79</td><td>47.78</td><td>0.37</td><td>13.59</td><td>28.62 ± 3.73</td><td>2</td></tr><tr><td></td><td>2013GW136 (o3e05)</td><td>2:1</td><td>7.42</td><td>47.74</td><td>0.34</td><td>6.66</td><td>21.20 ± 0.60</td><td>1</td></tr><tr><td>500877</td><td>2013JE64 (o3o18)</td><td>2:1</td><td>7.94</td><td>47.76</td><td>0.28</td><td>8.34</td><td>-2.60 ± 8.50</td><td>1</td></tr><tr><td>182397</td><td>2001QW297</td><td>9:4</td><td>6.94</td><td>51.67</td><td>0.23</td><td>17.04</td><td>27.25 ± 6.53</td><td>2</td></tr><tr><td>181867</td><td>1999CV118</td><td>7:3</td><td>7.35</td><td>52.77</td><td>0.29</td><td>5.48</td><td>25.97 ± 7.00</td><td>2</td></tr><tr><td>95625</td><td>2002GX32</td><td>7:3 (i)</td><td>7.77</td><td>53.07</td><td>0.38</td><td>13.94</td><td>37.59 ± 10.41</td><td>2</td></tr><tr><td>79978</td><td>1999CC158</td><td>12:5 (i)</td><td>5.70</td><td>53.88</td><td>0.28</td><td>18.75</td><td>24.74 ± 5.47</td><td>2</td></tr><tr><td>119878</td><td>2002CY224</td><td>12:5 (i)</td><td>6.18</td><td>53.85</td><td>0.35</td><td>15.73</td><td>32.65 ± 3.96</td><td>2</td></tr><tr><td>69988</td><td>1998WA31</td><td>5:2</td><td>7.56</td><td>55.18</td><td>0.43</td><td>9.46</td><td>17.97 ± 11.42</td><td>2</td></tr><tr><td>26375</td><td>1999DE9</td><td>5:2</td><td>5.03</td><td>55.41</td><td>0.42</td><td>7.62</td><td>20.45 ± 3.32</td><td>2</td></tr><tr><td>60621</td><td>2000FE8</td><td>5:2</td><td>6.76</td><td>55.29</td><td>0.40</td><td>5.87</td><td>12.19 ± 1.94</td><td>2</td></tr><tr><td>119068</td><td>2001KC77</td><td>5:2</td><td>7.10</td><td>55.09</td><td>0.36</td><td>12.90</td><td>20.80 ± 1.20</td><td>2</td></tr><tr><td></td><td>2001XQ254</td><td>5:2</td><td>7.92</td><td>55.37</td><td>0.44</td><td>7.11</td><td>14.71 ± 2.90</td><td>2</td></tr><tr><td>135571</td><td>2002GG32</td><td>5:2</td><td>7.43</td><td>55.29</td><td>0.35</td><td>14.68</td><td>38.31 ± 4.25</td><td>2</td></tr><tr><td>143707</td><td>2003UY117</td><td>5:2</td><td>5.71</td><td>55.55</td><td>0.41</td><td>7.54</td><td>21.97 ± 0.80</td><td>2,4</td></tr><tr><td></td><td>2004EG96</td><td>5:2</td><td>8.15</td><td>55.53</td><td>0.42</td><td>16.21</td><td>10.52 ± 2.64</td><td>2</td></tr><tr><td></td><td>2004HO79</td><td>5:2</td><td>7.38</td><td>55.21</td><td>0.41</td><td>5.62</td><td>26.40 ± 9.80</td><td>2</td></tr><tr><td></td><td>2004TT357</td><td>5:2</td><td>7.86</td><td>55.34</td><td>0.43</td><td>8.98</td><td>12.95 ± 3.40</td><td>2</td></tr><tr><td></td><td>2005SD278</td><td>5:2</td><td>6.17</td><td>55.50</td><td>0.28</td><td>17.85</td><td>18.02 ± 1.55</td><td>2</td></tr><tr><td></td><td>2009YG19</td><td>5:2</td><td>6.02</td><td>55.67</td><td>0.41</td><td>5.15</td><td>26.60 ± 4.90</td><td>4</td></tr><tr><td></td><td>2013GY136 (o3e09)</td><td>5:2</td><td>7.32</td><td>55.53</td><td>0.41</td><td>10.88</td><td>6.40 ± 0.60</td><td>1</td></tr><tr><td></td><td>2013JK64 (o3o11)</td><td>5:2</td><td>7.69</td><td>55.25</td><td>0.41</td><td>11.08</td><td>35.97 ± 0.48</td><td>1</td></tr><tr><td></td><td>2000CQ105</td><td>13:5 (i)</td><td>6.27</td><td>56.99</td><td>0.39</td><td>19.68</td><td>3.73 ± 1.69</td><td>2,3</td></tr><tr><td></td><td>1999RJ215</td><td>11:4 (i)</td><td>7.86</td><td>59.27</td><td>0.41</td><td>19.71</td><td>13.85 ± 1.82</td><td>2,3</td></tr><tr><td>500879</td><td>2013JH64 (o3o34)</td><td>11:4</td><td>5.60</td><td>59.23</td><td>0.38</td><td>13.73</td><td>15.60 ± 0.80</td><td>1</td></tr><tr><td>136120</td><td>2003LG7</td><td>3:1</td><td>8.58</td><td>62.16</td><td>0.48</td><td>20.10</td><td>12.99 ± 4.66</td><td>2</td></tr><tr><td>385607</td><td>2005E0297</td><td>3:1</td><td>7.47</td><td>62.70</td><td>0.34</td><td>25.04</td><td>13.77 ± 4.80</td><td>2</td></tr><tr><td></td><td>2006QJ181</td><td>3:1</td><td>7.15</td><td>62.58</td><td>0.49</td><td>20.06</td><td>11.99 ± 3.59</td><td>2</td></tr><tr><td>126619</td><td>2002CX154</td><td>11:3 (i)</td><td>7.51</td><td>71.27</td><td>0.47</td><td>15.98</td><td>21.27 ± 6.51</td><td>2</td></tr><tr><td></td><td>2007TA418</td><td>15:4 (i)</td><td>7.39</td><td>72.42</td><td>0.50</td><td>21.96</td><td>14.00 ± 2.30</td><td>3</td></tr><tr><td>184212</td><td>2004PB112 (o5s16PD)</td><td>27:4 (i)</td><td>7.39</td><td>107.52</td><td>0.67</td><td>15.43</td><td>20.30 ± 2.50</td><td>1</td></tr><tr><td></td><td>2007TC434 (o4h39)</td><td>9:1</td><td>7.13</td><td>129.63</td><td>0.69</td><td>26.47</td><td>11.90 ± 1.30</td><td>1</td></tr><tr><td>148209</td><td>2000CR105</td><td>20:1 (i)</td><td>6.44</td><td>221.58</td><td>0.80</td><td>22.76</td><td>13.22 ± 8.40</td><td>2</td></tr><tr><td>474640</td><td>2004VN112</td><td>36:1 (i)</td><td>6.22</td><td>327.38</td><td>0.85</td><td>25.55</td><td>11.50 ± 4.25</td><td>2</td></tr></tbody></table><table><thead><tr><th colspan="9">Centaurs</th></tr></thead><tbody><tr><td>2060</td><td>1977UB (Chiron)</td><td>cen</td><td>6.53</td><td>13.85</td><td>0.39</td><td>6.91</td><td>2.62 ± 0.88</td><td>2</td></tr><tr><td>5145</td><td>1992AD (Pholus)</td><td>cen</td><td>7.55</td><td>20.31</td><td>0.57</td><td>24.70</td><td>52.82 ± 8.23</td><td>2</td></tr><tr><td>7066</td><td>1993HA2 (Nessus)</td><td>cen</td><td>9.42</td><td>24.54</td><td>0.52</td><td>15.65</td><td>43.35 ± 2.99</td><td>2,3</td></tr><tr><td>8405</td><td>1995GO (Asbolus)</td><td>cen</td><td>9.51</td><td>18.06</td><td>0.62</td><td>17.62</td><td>13.61 ± 2.57</td><td>2</td></tr><tr><td>10199</td><td>1997CU26 (Chariklo)</td><td>cen</td><td>6.83</td><td>15.77</td><td>0.17</td><td>23.39</td><td>14.43 ± 0.91</td><td>2</td></tr><tr><td>49036</td><td>1998QM107 (Pelion)</td><td>cen</td><td>10.42</td><td>20.02</td><td>0.14</td><td>9.37</td><td>20.84 ± 10.25</td><td>2</td></tr><tr><td>52975</td><td>1998TF35 (Cyllarus)</td><td>cen</td><td>8.71</td><td>26.16</td><td>0.38</td><td>12.64</td><td>34.99 ± 1.29</td><td>2,3</td></tr><tr><td>31824</td><td>1999UG5 (Elatus)</td><td>cen</td><td>10.67</td><td>12.59</td><td>0.41</td><td>5.60</td><td>24.71 ± 1.08</td><td>2,3</td></tr><tr><td>121725</td><td>1999XX143 (Aphidas)</td><td>cen</td><td>8.83</td><td>17.94</td><td>0.46</td><td>6.78</td><td>18.70 ± 2.50</td><td>3</td></tr><tr><td>54598</td><td>2000QC243 (Bienor)</td><td>cen</td><td>7.96</td><td>16.50</td><td>0.20</td><td>20.74</td><td>8.93 ± 4.18</td><td>2</td></tr><tr><td>32532</td><td>2001PT13 (Thereus)</td><td>cen</td><td>9.58</td><td>10.64</td><td>0.20</td><td>20.35</td><td>10.31 ± 2.43</td><td>2</td></tr><tr><td>119315</td><td>2001SQ73</td><td>cen</td><td>9.24</td><td>17.45</td><td>0.18</td><td>17.43</td><td>8.83 ± 0.77</td><td>2,3</td></tr><tr><td>427507</td><td>2002DH5</td><td>cen</td><td>10.34</td><td>22.03</td><td>0.37</td><td>22.47</td><td>2.61 ± 3.63</td><td>2</td></tr><tr><td>55576</td><td>2002GB10 (Amycus)</td><td>cen</td><td>7.90</td><td>25.03</td><td>0.39</td><td>13.34</td><td>37.88 ± 2.16</td><td>2,4</td></tr><tr><td>83982</td><td>2002GO9 (Crantor)</td><td>cen</td><td>8.96</td><td>19.40</td><td>0.28</td><td>12.77</td><td>24.61 ± 1.08</td><td>2,3</td></tr><tr><td>95626</td><td>2002GZ32</td><td>cen</td><td>6.88</td><td>23.06</td><td>0.22</td><td>15.02</td><td>23.40 ± 5.41</td><td>2</td></tr></tbody></table>
<table><caption>**Table 4**
(Continued)</caption><thead><tr><th>MPC<br>Number</th><th>Target<br>Designation</th><th>Orbit<br>Classification</th><th>H<sub>r</sub></th><th>a (au)</th><th>e</th><th>inc (°)</th><th>Spectral slope<br>(%/10<sup>3</sup> Å)</th><th>References</th></tr></thead><tbody><tr><td rowspan="2">250112</td><td>2002KY14</td><td>cen</td><td>10.19</td><td>12.55</td><td>0.31</td><td>19.46</td><td>36.28 ± 2.23</td><td>2,4</td></tr><tr><td>2002PQ152</td><td>cen</td><td>9.50</td><td>25.79</td><td>0.19</td><td>9.35</td><td>40.40 ± 6.70</td><td>4</td></tr><tr><td></td><td>2002QX47</td><td>cen</td><td>8.69</td><td>25.43</td><td>0.37</td><td>7.28</td><td>2.00 ± 3.90</td><td>4</td></tr><tr><td>120061</td><td>2003CO1</td><td>cen</td><td>9.38</td><td>20.72</td><td>0.47</td><td>19.74</td><td>11.03 ± 5.08</td><td>2</td></tr><tr><td rowspan="2">136204</td><td>2003WL7</td><td>cen</td><td>8.76</td><td>20.14</td><td>0.26</td><td>11.17</td><td>13.10 ± 1.56</td><td>2,4</td></tr><tr><td>2004QQ26</td><td>cen</td><td>9.91</td><td>22.93</td><td>0.15</td><td>21.45</td><td>5.50 ± 1.90</td><td>3</td></tr><tr><td>160427</td><td>2005RL43</td><td>cen</td><td>7.99</td><td>24.51</td><td>0.04</td><td>12.26</td><td>31.30 ± 1.30</td><td>3</td></tr><tr><td>447178</td><td>2005RO43</td><td>cen</td><td>7.16</td><td>28.85</td><td>0.52</td><td>35.46</td><td>8.84 ± 0.78</td><td>3,4</td></tr><tr><td>248835</td><td>2006SX368</td><td>cen</td><td>9.27</td><td>22.08</td><td>0.46</td><td>36.30</td><td>11.00 ± 2.10</td><td>4</td></tr><tr><td>309139</td><td>2006XQ51</td><td>cen</td><td>10.15</td><td>15.85</td><td>0.38</td><td>31.61</td><td>4.90 ± 3.00</td><td>4</td></tr><tr><td>187661</td><td>2007JG43</td><td>cen</td><td>9.20</td><td>23.95</td><td>0.40</td><td>33.13</td><td>7.50 ± 0.70</td><td>3</td></tr><tr><td rowspan="2">341275</td><td>2007RG283</td><td>cen</td><td>8.54</td><td>19.91</td><td>0.23</td><td>28.76</td><td>11.00 ± 3.20</td><td>4</td></tr><tr><td>2007RH283</td><td>cen</td><td>8.49</td><td>15.92</td><td>0.34</td><td>21.35</td><td>7.35 ± 0.66</td><td>3,4</td></tr><tr><td rowspan="2">281371</td><td>2008FC76</td><td>cen</td><td>9.41</td><td>14.69</td><td>0.31</td><td>27.12</td><td>29.82 ± 2.34</td><td>2,4</td></tr><tr><td>2010TH</td><td>cen</td><td>9.05</td><td>18.55</td><td>0.32</td><td>26.70</td><td>10.00 ± 2.10</td><td>4</td></tr><tr><td>449097</td><td>2012UT68</td><td>cen</td><td>9.44</td><td>20.25</td><td>0.38</td><td>15.41</td><td>32.70 ± 2.50</td><td>4</td></tr><tr><td>463368</td><td>2012VU85</td><td>cen</td><td>8.39</td><td>29.23</td><td>0.31</td><td>15.07</td><td>29.00 ± 5.00</td><td>4</td></tr><tr><td rowspan="3">459865</td><td>2013XZ8</td><td>cen</td><td>9.56</td><td>13.43</td><td>0.37</td><td>22.53</td><td>8.90 ± 2.10</td><td>4</td></tr><tr><td>2014UJ225 (o4h01)</td><td>cen</td><td>10.29</td><td>23.19</td><td>0.38</td><td>21.32</td><td>12.50 ± 1.10</td><td>1</td></tr><tr><td>2015RV245 (o5s05)</td><td>cen</td><td>10.10</td><td>21.98</td><td>0.48</td><td>15.39</td><td>9.90 ± 2.30</td><td>1</td></tr><tr><td colspan="9" style="text-align: center; font-weight: bold;">Scattering Disk Objects</td></tr><tr><td>15874</td><td>1996TL66</td><td>sca</td><td>5.34</td><td>83.42</td><td>0.58</td><td>23.98</td><td>1.56 ± 3.89</td><td>2</td></tr><tr><td>33128</td><td>1998BU48</td><td>sca</td><td>7.18</td><td>33.27</td><td>0.38</td><td>14.25</td><td>29.03 ± 10.54</td><td>2</td></tr><tr><td>181902</td><td>1999RD215</td><td>sca</td><td>8.11</td><td>122.40</td><td>0.69</td><td>25.93</td><td>12.70 ± 6.30</td><td>2</td></tr><tr><td>91554</td><td>1999RZ215</td><td>sca</td><td>8.07</td><td>101.73</td><td>0.70</td><td>25.50</td><td>12.71 ± 2.20</td><td>2,3</td></tr><tr><td rowspan="2">29981</td><td>1999TD10</td><td>sca</td><td>8.85</td><td>95.78</td><td>0.87</td><td>5.96</td><td>10.51 ± 0.74</td><td>2,3</td></tr><tr><td>2000CQ104</td><td>sca (i)</td><td>8.08</td><td>36.64</td><td>0.23</td><td>13.51</td><td>8.38 ± 3.08</td><td>2</td></tr><tr><td>60608</td><td>2000EE173</td><td>sca</td><td>8.11</td><td>49.20</td><td>0.54</td><td>5.95</td><td>13.52 ± 2.83</td><td>2,3</td></tr><tr><td>87269</td><td>2000OO67</td><td>sca</td><td>9.33</td><td>562.10</td><td>0.96</td><td>20.07</td><td>25.98 ± 6.34</td><td>2</td></tr><tr><td>87555</td><td>2000QB243</td><td>sca</td><td>8.77</td><td>34.74</td><td>0.56</td><td>6.79</td><td>12.57 ± 2.56</td><td>2,3</td></tr><tr><td>82158</td><td>2001FP185</td><td>sca</td><td>6.20</td><td>215.84</td><td>0.84</td><td>30.80</td><td>19.62 ± 3.13</td><td>2</td></tr><tr><td rowspan="3">82155</td><td>2001FZ173</td><td>sca (i)</td><td>6.08</td><td>85.33</td><td>0.62</td><td>12.72</td><td>19.43 ± 2.10</td><td>2</td></tr><tr><td>2001KG77</td><td>sca (i)</td><td>8.57</td><td>61.67</td><td>0.45</td><td>15.48</td><td>7.90 ± 6.30</td><td>2</td></tr><tr><td>2002GB32</td><td>sca</td><td>7.90</td><td>206.81</td><td>0.83</td><td>14.19</td><td>19.17 ± 1.71</td><td>2</td></tr><tr><td>469442</td><td>2002GG166 (o3e01)</td><td>sca</td><td>7.73</td><td>34.42</td><td>0.59</td><td>7.71</td><td>11.10 ± 0.60</td><td>1</td></tr><tr><td>73480</td><td>2002PN34</td><td>sca</td><td>8.73</td><td>30.90</td><td>0.57</td><td>16.64</td><td>14.20 ± 5.50</td><td>2</td></tr><tr><td>65489</td><td>2003FX128 (Ceto)</td><td>sca</td><td>6.55</td><td>100.62</td><td>0.82</td><td>22.30</td><td>13.86 ± 0.87</td><td>2,3</td></tr><tr><td>506479</td><td>2003HB57</td><td>sca</td><td>7.64</td><td>159.63</td><td>0.76</td><td>15.50</td><td>14.38 ± 2.30</td><td>2</td></tr><tr><td>149560</td><td>2003QZ91</td><td>sca</td><td>8.03</td><td>41.51</td><td>0.48</td><td>34.82</td><td>12.80 ± 2.60</td><td>3</td></tr><tr><td>469750</td><td>2005PU21</td><td>sca</td><td>6.38</td><td>174.64</td><td>0.83</td><td>6.17</td><td>32.80 ± 2.01</td><td>2</td></tr><tr><td>145474</td><td>2005SA278</td><td>sca</td><td>6.53</td><td>92.28</td><td>0.64</td><td>16.27</td><td>8.10 ± 1.60</td><td>3</td></tr><tr><td rowspan="7">308933</td><td>2006SQ372</td><td>sca</td><td>7.89</td><td>799.41</td><td>0.97</td><td>19.46</td><td>25.44 ± 1.04</td><td>2,3</td></tr><tr><td>2007TG422</td><td>sca</td><td>6.44</td><td>503.03</td><td>0.93</td><td>18.59</td><td>14.44 ± 3.00</td><td>2</td></tr><tr><td>2007VJ305</td><td>sca</td><td>6.96</td><td>192.12</td><td>0.82</td><td>11.98</td><td>15.50 ± 2.26</td><td>2</td></tr><tr><td>2013JO64 (o3o14)</td><td>sca</td><td>8.00</td><td>143.31</td><td>0.75</td><td>8.58</td><td>8.90 ± 1.40</td><td>1</td></tr><tr><td>2013UR15 (o3l01)</td><td>sca</td><td>10.89</td><td>55.82</td><td>0.72</td><td>22.25</td><td>13.70 ± 2.90</td><td>1</td></tr><tr><td>2014UQ229 (o4h03)</td><td>sca</td><td>9.55</td><td>49.87</td><td>0.78</td><td>5.68</td><td>35.90 ± 1.20</td><td>1</td></tr><tr><td>2015RU245 (o5t04)</td><td>sca</td><td>9.32</td><td>30.99</td><td>0.29</td><td>13.75</td><td>22.90 ± 0.90</td><td>1</td></tr><tr><td>2015RW245 (o5s06)</td><td>sca</td><td>8.53</td><td>56.47</td><td>0.53</td><td>13.30</td><td>18.30 ± 1.10</td><td>1</td></tr><tr><td colspan="9" style="text-align: center; font-weight: bold;">Detached Objects</td></tr><tr><td>181874</td><td>1999CF119</td><td>det</td><td>7.23</td><td>88.73</td><td>0.56</td><td>19.71</td><td>13.19 ± 6.01</td><td>2</td></tr><tr><td rowspan="2">40314</td><td>1999HW11</td><td>det (i)</td><td>6.96</td><td>52.68</td><td>0.26</td><td>17.20</td><td>13.38 ± 6.14</td><td>2</td></tr><tr><td>1999KR16</td><td>det</td><td>5.82</td><td>48.76</td><td>0.30</td><td>24.82</td><td>38.31 ± 1.37</td><td>2,3</td></tr><tr><td rowspan="2">60458</td><td>2000AF255</td><td>det</td><td>5.96</td><td>48.59</td><td>0.25</td><td>30.88</td><td>41.96 ± 3.08</td><td>2,3</td></tr><tr><td>2000CM114</td><td>det</td><td>7.20</td><td>59.49</td><td>0.40</td><td>19.70</td><td>14.54 ± 2.18</td><td>2</td></tr><tr><td>118702</td><td>2000OM67</td><td>det</td><td>7.33</td><td>98.33</td><td>0.60</td><td>23.36</td><td>15.67 ± 3.31</td><td>2</td></tr><tr><td>469333</td><td>2000PE30</td><td>det</td><td>6.09</td><td>54.31</td><td>0.34</td><td>18.42</td><td>4.07 ± 3.64</td><td>2</td></tr><tr><td rowspan="3">182223</td><td>2000YC2</td><td>det (i)</td><td>7.33</td><td>58.47</td><td>0.38</td><td>19.88</td><td>0.86 ± 6.01</td><td>2</td></tr><tr><td>2001FM194</td><td>det</td><td>7.71</td><td>54.05</td><td>0.36</td><td>28.69</td><td>12.15 ± 3.79</td><td>2</td></tr><tr><td>2001QX322</td><td>det (i)</td><td>6.48</td><td>58.10</td><td>0.39</td><td>28.53</td><td>14.41 ± 1.83</td><td>2,3</td></tr><tr><td>2003FZ129</td><td></td><td>det</td><td>7.22</td><td>61.79</td><td>0.39</td><td>5.79</td><td>9.88 ± 2.29</td><td>2</td></tr></tbody></table>
<table><caption>**Table 4**
(Continued)</caption><thead><tr><th>MPC<br>Number</th><th>Target<br>Designation</th><th>Orbit<br>Classification</th><th>H<sub>r</sub></th><th>a (au)</th><th>e</th><th>inc (°)</th><th>Spectral slope<br>(%/10<sup>3</sup> Å)</th><th>References</th></tr></thead><tbody><tr><td></td><td>2004OJ14</td><td>det</td><td>7.25</td><td>55.30</td><td>0.29</td><td>22.46</td><td>16.60 ± 2.33</td><td>2</td></tr><tr><td>470309</td><td>2007JK43</td><td>det*</td><td>7.20</td><td>46.16</td><td>0.49</td><td>44.89</td><td>14.30 ± 0.80</td><td>3</td></tr><tr><td>496315</td><td>2013GP136 (o3e39)</td><td>det (i)</td><td>6.42</td><td>150.15</td><td>0.73</td><td>33.54</td><td>18.50 ± 1.00</td><td>1</td></tr><tr><td></td><td>2013JL64 (o3o29)</td><td>det (i)</td><td>7.03</td><td>56.77</td><td>0.37</td><td>27.67</td><td>15.20 ± 1.90</td><td>1</td></tr></tbody></table>

**Note.** Dynamical classification: *cla* = classical belt; *sca* = scattering disk; *cen* = centaur; *x:y* = resonators, where *x* and *y* indicate the specific mean-motion resonance. (i) Insecure classification. *2007JK43 is nearly on a Uranus crossing orbit and has a semimajor axis that falls within the classical belt. It is classified as detached because no close encounter happens during the 10 My integration of its orbit. Orbital elements are in the barycentric reference frame and were computed using the Bernstein & Khushalani (2000) orbit-fitting procedure. Dynamical classification uses the Gladman et al. (2008) classification scheme. See Bannister et al. (2016) for additional information. *H*<sub>r</sub> are the absolute magnitudes provided in “Ref.” and converted to *r*-band using the Synphot/STSDAS tool. Peixinho et al. (2015) do not report the absolute magnitude for several objects in their data set. We retrieved these magnitudes from the database of Minor Bodies in the Outer Solar System (MBOSS) (Hainaut et al. 2012) when available, or from the original papers (Boehnhardt et al. 2002; Benecchi et al. 2009, 2011; Sheppard 2010) otherwise.

**References.** 1 = Col-OSSOS (Pike et al. 2017, Schwamb et al. 2018 and this work); 2 = Peixinho et al. (2015) and references herein; 3: Fraser & Brown (2012), Fraser et al. (2015); 4: Tegler et al. (2016).

(This table is available in its entirety in machine-readable form.)

## Appendix B

### Effects of Sample Selection and Color Classification on Our Analysis

We explore here the effects of sample selection and color classification on our statistical analysis. More specifically, we test whether our statistical results are sensitive to the magnitude and inclination cuts used to define our data set. We further test whether varying the spectral slope value used to divide our data set into the gray and red color classes make any difference.

Large TNOs are characterized by different surfaces from those of small objects (see introduction, Section 1). The transition size between large and small TNO surfaces is usually proposed at *H*<sub>r</sub> ~ 6 (Fraser & Brown 2012; Peixinho et al. 2012), i.e., at a diameter of 280 km assuming a visible albedo of 0.09 (Lacerda et al. 2014). In this work, our statistical analysis was performed including all objects with *H*<sub>r</sub> > 5 (Section 2.1). We test here whether our analysis was biased by the presence in our data set of large 5 < *H*<sub>r</sub> < 6 TNOs by repeating our analysis considering only objects with *H*<sub>r</sub> > 6. Our results are summarized at Table 5. They are statistically similar to those presented in Section 3.

We further test the effects of changing the inclination cut used for the selection of our data set. Our inclination cut of *i* = 5° should ensure minimal contamination of our sample by cold classicals (Section 2.1). Let us assume, however, that our sample comprises a significant fraction of cold classical objects up to *i* = 6°. Repeating our statistical analysis only for objects located at orbital inclinations of *i* > 6° reveals similar results to that of Section 3 (Table 5).

Finally, we test if our statistics are sensitive to the spectral slope value *s*<sub>lim</sub> used to divide our data set into the gray and red TNO classes. By adopting the two most extreme spectral slope values from Table 2 (17.4% and 23.8%/(10<sup>3</sup> Å)), we find that our results hold for the overall data set (Table 6). Choosing the lowest spectral slope value weakens our results for the hot classical population as one previously unclassified high-*i* object falls into the red class. Choosing the highest value weakens our results for the centaur/scattering/detached objects as a significantly higher fraction of low-*i* objects fall into the gray class.

<table><caption>**Table 5**
Effects of Magnitude and Inclination Cuts</caption><thead><tr><th rowspan="2"></th><th colspan="5">cen</th></tr><tr><th>hc</th><th>res</th><th>+sca+det</th><th>all</th><th>all-cen</th></tr></thead><tbody><tr><td></td><td colspan="5" style="text-align: center;"><i>H</i><sub>r</sub> > 6</td></tr><tr><td>2-sample KS test</td><td>98.0%</td><td>94.1%</td><td>99.3%</td><td>&gt;<b>99.9%</b></td><td><b>99.9%</b></td></tr><tr><td>2-sample AD test</td><td>95.8%</td><td>98.9%</td><td>99.1%</td><td>&gt;<b>99.9%</b></td><td>&gt;<b>99.9%</b></td></tr><tr><td>T-test of means</td><td>99.3%</td><td>99.3%</td><td>99.4%</td><td>&gt;<b>99.9%</b></td><td>&gt;<b>99.9%</b></td></tr><tr><td>Spearman<br>correlation</td><td>96.0%</td><td>96.8%</td><td>89.7%</td><td><b>99.8%</b></td><td>99.3%</td></tr><tr><td></td><td colspan="5" style="text-align: center;"><i>i</i> > 6°</td></tr><tr><td>2-sample KS test</td><td>98.4%</td><td>87.3%</td><td>96.5%</td><td>&gt;<b>99.9%</b></td><td><b>99.7%</b></td></tr><tr><td>2-sample AD test</td><td>96.4%</td><td>89.2%</td><td>95.1%</td><td>&gt;<b>99.9%</b></td><td><b>99.9%</b></td></tr><tr><td>T-test of means</td><td>99.4%</td><td>97.4%</td><td>97.2%</td><td>&gt;<b>99.9%</b></td><td>&gt;<b>99.9%</b></td></tr><tr><td>Spearman<br>correlation</td><td>92.6%</td><td>89.2%</td><td>77.9%</td><td>98.8%</td><td>97.7%</td></tr></tbody></table>

**Note.** Same acronyms as those in Table 3.

<table><caption>**Table 6**
Effects of Color Classification</caption><thead><tr><th></th><th>hc</th><th>res</th><th>cen+sca+det</th><th>all</th><th>all-cen</th></tr></thead><tbody><tr><td></td><td colspan="5" style="text-align: center;"><i>s</i><sub>lim</sub> = 17.4 %/(10<sup>3</sup> Å)</td></tr><tr><td>2-sample KS test</td><td>65.8%</td><td>92.8%</td><td>97.8%</td><td>&gt;<b>99.9%</b></td><td>99.5%</td></tr><tr><td>2-sample AD test</td><td>70.1%</td><td>97.2%</td><td>97.1%</td><td>&gt;<b>99.9%</b></td><td>99.5%</td></tr><tr><td>T-test of means</td><td>93.2%</td><td>98.5%</td><td>98.1%</td><td>&gt;<b>99.9%</b></td><td><b>99.8%</b></td></tr><tr><td></td><td colspan="5" style="text-align: center;"><i>s</i><sub>lim</sub> = 23.8 %/(10<sup>3</sup> Å)</td></tr><tr><td>2-sample KS test</td><td>99.3%</td><td>87.5%</td><td>78.5%</td><td><b>99.9%</b></td><td>98.6%</td></tr><tr><td>2-sample AD test</td><td>94.5%</td><td>96.7%</td><td>72.7%</td><td><b>99.8%</b></td><td><b>99.8%</b></td></tr><tr><td>T-test of means</td><td>99.0%</td><td>98.2%</td><td>90.6%</td><td>&gt;<b>99.9%</b></td><td><b>99.9%</b></td></tr></tbody></table>

**Note.** Same acronyms as those in Table 3.

### ORCID iDs

Michaël Marsset https://orcid.org/0000-0001-8617-2425

Wesley C. Fraser https://orcid.org/0000-0001-6680-6558

Rosemary E. Pike https://orcid.org/0000-0003-4797-5262

Michele T. Bannister https://orcid.org/0000-0003-3257-4490
Megan E. Schwamb https://orcid.org/0000-0003-4365-1455

Kathryn Volk https://orcid.org/0000-0001-8736-236X

J. J. Kavelaars https://orcid.org/0000-0001-7032-5255

Mike Alexandersen https://orcid.org/0000-0003-4143-8589

Ying-Tung Chen https://orcid.org/0000-0001-7244-6069

Matthew J. Lehner https://orcid.org/0000-0003-4077-0985

Nuno Peixinho https://orcid.org/0000-0002-6830-476X

Jean-Marc Petit https://orcid.org/0000-0003-0407-2266

## References

Alexandersen, M., Gladman, B., Kavelaars, J. J., et al. 2016, AJ, 152, 111

Allen, R. L., Bernstein, G. M., & Malhotra, R. 2002, AJ, 124, 2949

Alvarez-Candal, A., Barucci, M. A., Merlin, F., Guilbert, A., & de Bergh, C. 2007, A&A, 475, 369

Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33

Bannister, M. T., Gladman, B. J., Kavelaars, J. J., et al. 2018, ApJS, 236, 18

Bannister, M. T., Kavelaars, J. J., Petit, J.-M., et al. 2016, AJ, 152, 70

Barkume, K. M., Brown, M. E., & Schaller, E. L. 2008, AJ, 135, 55

Barucci, M. A., Alvarez-Candal, A., Merlin, F., et al. 2011, Icar, 214, 297

Barucci, M. A., Brown, M. E., Emery, J. P., & Merlin, F. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizaon Press), 143

Benecchi, S. D., Noll, K. S., Grundy, W. M., et al. 2009, Icar, 200, 292

Benecchi, S. D., Noll, K. S., Stephens, D. C., Grundy, W. M., & Rawlins, J. 2011, Icar, 213, 693

Bernstein, G., & Khushalani, B. 2000, AJ, 120, 3323

Bernstein, G. M., Trilling, D. E., Allen, R. L., et al. 2004, AJ, 128, 1364

Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393

Boehnhardt, H., Delsanti, A., Barucci, A., et al. 2002, A&A, 395, 297

Brasser, R., & Schwamb, M. E. 2015, MNRAS, 446, 3788

Brown, M. E. 2001, AJ, 121, 2804

Brown, M. E. 2012, AREPS, 40, 467

Brown, M. E., Barkume, K. M., Ragozzine, D., & Schaller, E. L. 2007, Natur, 446, 294

Brown, M. E., Schaller, E. L., & Fraser, W. C. 2011, ApJL, 739, L60

Brucker, M. J., Grundy, W. M., Stansberry, J. A., et al. 2009, Icar, 201, 284

Chakravarti, I. M., Laha, R. G., & J. R. 1967, Handbook of Methods of Applied Statistics, Vol. 1 (New York: Wiley)

Chen, Y.-T., Lin, H. W., Holman, M. J., et al. 2016, ApJL, 827, L24

Chiang, E. I., & Brown, M. E. 1999, AJ, 118, 1411

Delsanti, A., Hainaut, O., Jourdeuil, E., et al. 2004, A&A, 417, 1145

Doressoundiram, A., Peixinho, N., de Bergh, C., et al. 2002, AJ, 124, 2279

Doressoundiram, A., Peixinho, N., Doucet, C., et al. 2005, Icar, 174, 90

Duncan, M. J., & Levison, H. F. 1997, Sci, 276, 1670

Efron, B., & Tibshirani, R. J. 1993, An Introduction to the Bootstrap (London/ Boca Raton, FL: Chapman and Hall/CRC Press)

Elliott, J. L., Kern, S. D., Clancy, K. B., et al. 2005, AJ, 129, 1117

Emel'yanenko, V. V., Asher, D. J., & Bailey, M. E. 2003, MNRAS, 338, 443

Fornasier, S., Dotto, E., Barucci, M. A., & Barbieri, C. 2004, A&A, 422, L43

Fraser, W., Alexandersen, M., Schwamb, M. E., et al. 2016, AJ, 151, 158

Fraser, W. C., Bannister, M. T., Pike, R. E., et al. 2017, NatAs, 1, 0088

Fraser, W. C., & Brown, M. E. 2012, ApJ, 749, 33

Fraser, W. C., Brown, M. E., & Glass, F. 2015, ApJ, 804, 31

Fraser, W. C., Brown, M. E., Morbidelli, A., Parker, A., & Batygin, K. 2014, ApJ, 782, 100

Fraser, W. C., Brown, M. E., & Schwamb, M. E. 2010, Icar, 210, 944

Fraser, W. C., & Kavelaars, J. J. 2009, AJ, 137, 72

Fraser, W. C., Kavelaars, J. J., Holman, M. J., et al. 2008, Icar, 195, 827

Fuentes, C. I., George, M. R., & Holman, M. J. 2009, ApJ, 696, 91

Fuentes, C. I., & Holman, M. J. 2008, AJ, 136, 83

Gil-Hutton, R. 2002, P&SS, 50, 57

Gladman, B. 2005, Sci, 307, 71

Gladman, B., Holman, M., Grav, T., et al. 2002, Icar, 157, 269

Gladman, B., Kavelaars, J., Petit, J. M., et al. 2009, ApJL, 697, L91

Gladman, B., Kavelaars, J. J., Nicholson, P. D., Loredo, T. J., & Burns, J. A. 1998, AJ, 116, 2042

Gladman, B., Kavelaars, J. J., Petit, J.-M., et al. 2001, AJ, 122, 1051

Gladman, B., Lawler, S. M., Petit, J. M., et al. 2012, AJ, 144, 23

Gladman, B., Marsden, B. G., & Vanlaerhoven, C. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizaon Press), 43

Gomes, R. S. 2003, Icar, 161, 404

Gomes, R. S., Fern Ndez, J. A., Gallardo, T., & Brunini, A. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizaon Press), 259

Gomes, R. S., Gallardo, T., Fernández, J. A., & Brunini, A. 2005, CeMDA, 91, 109

Gulbis, A. A. S., Elliot, J. L., Adams, E. R., et al. 2010, AJ, 140, 350

Hainaut, O. R., Boehnhardt, H., & Protopapa, S. 2012, A&A, 546, 115

Hainaut, O. R., & Delsanti, A. C. 2002, A&A, 389, 641

Hodapp, K. W., Jensen, J. B., Irwin, E. M., et al. 2003, PASP, 115, 1388

Hook, I. M., Jørgensen, I., Allington-Smith, J. R., et al. 2004, PASP, 116, 425

Hunter, J. D. 2007, CSE, 9, 90

Jewitt, D. 2015, AJ, 150, 201

Jewitt, D. C., & Luu, J. X. 2001, AJ, 122, 2099

Kaib, N. A., Roskar, R., & Quinn, T. 2011, Icar, 215, 491

Lacerda, P., Fornasier, S., Lellouch, E., et al. 2014, ApJL, 793, L2

Lawler, S. M., Kavelaars, J. J., Alexandersen, M., et al. 2018a, FrASS, 5, 14

Lawler, S. M., Shankman, C., Kavelaars, J. J., et al. 2018b, AJ, 155, 197

Levison, H., Morbidelli, A., Vanlaerhoven, C., Gomes, R., & Tsiganis, K. 2008, Icar, 196, 258

Levison, H. F., & Duncan, M. J. 1997, Icar, 127, 13

Levison, H. F., & Stern, S. A. 2001, AJ, 121, 1730

Lykawka, P. S., & Mukai, T. 2007, Icar, 189, 213

Magnier, E. A., Schlafly, E., Finkbeiner, D., et al. 2013, ApJS, 205, 20

Magnier, E. A., Schlafly, E. F., Finkbeiner, D. P., et al. 2016, arXiv:1612.05242

Malhotra, R. 1995, AJ, 110, 420

Millis, R. L., Buie, M. W., Wasserman, L. H., et al. 2002, AJ, 123, 2083

Moroz, L. V., Baratta, G., Distefano, E., et al. 2003, EM&P, 92, 279

Nesvorny, D. 2015a, AJ, 150, 73

Nesvorny, D. 2015b, AJ, 150, 68

Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, Journal of Machine Learning Research, 12, 2825

Peixinho, N., Delsanti, A., & Doressoundiram, A. 2015, A&A, 577, A35

Peixinho, N., Delsanti, A., Guilbert-Lepoutre, A., Gafeira, R., & Lacerda, P. 2012, A&A, 546, A86

Peixinho, N., Doressoundiram, A., Delsanti, A., et al. 2003, A&A, 410, L29

Peixinho, N., Lacerda, P., & Jewitt, D. 2008, AJ, 136, 1837

Petit, J.-M., Holman, M. J., Gladman, B. J., et al. 2006, MNRAS, 365, 429

Petit, J. M., Kavelaars, J. J., Gladman, B. J., et al. 2011, AJ, 142, 131

Petit, J.-M., Kavelaars, J. J., Gladman, B. J., et al. 2017, AJ, 153, 236

Pike, R. E., Fraser, W. C., Schwamb, M. E., et al. 2017, AJ, 154, 101

Rabinowitz, D., Schwamb, M. E., Hadjiyska, E., & Tourtellotte, S. 2012, AJ, 144, 140

Ragozzine, D., & Brown, M. E. 2007, AJ, 134, 2160

Santos-Sanz, P., Ortiz, J. L., Barrera, L., & Boehnhardt, H. 2009, A&A, 494, 693

Schaller, E. L., & Brown, M. E. 2007, ApJL, 659, L61

Schaller, E. L., & Brown, M. E. 2008, ApJL, 684, L107

Schwamb, M. E., Brown, M. E., Rabinowitz, D. L., & Ragozzine, D. 2010, ApJ, 720, 1691

Schwamb, M. E., Fraser, W. C., Bannister, M. T., et al. 2018, ApJS, submitted (arXiv:1809.08501)

SDSS Collaboration, Albareti, F. D., Allende Prieto, C., et al. 2017, ApJS, 233, 25

Shankman, C., Gladman, B. J., Kaib, N., Kavelaars, J. J., & Petit, J. M. 2013, ApJL, 764, L2

Shankman, C., Kavelaars, J. J., Gladman, B. J., et al. 2016, AJ, 151, 31

Sheppard, S. S. 2010, AJ, 139, 1394

Snedecor, G. W., & Cochran, W. G. 1989, Statistical Methods (8th ed.; Iowa City, IA: Univ. Iowa Press)

Snodgrass, C., Carry, B., Dumas, C., & Hainaut, O. 2010, A&A, 511, A72

Stephens, M. A. 1974, J. Am. Stat. Assoc., 69, 730

Stern, S. A. 2002, AJ, 124, 2297

Tegler, S. C., Bauer, J. M., Romanishin, W., & Peixinho, N. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizaon Press), 105

Tegler, S. C., & Romanishin, W. 1998, Natur, 392, 49

Tegler, S. C., & Romanishin, W. 2000, Natur, 407, 979

Tegler, S. C., & Romanishin, W. 2003, Icar, 161, 181

Tegler, S. C., Romanishin, W., Consolmagno, G. J., & J. S. 2016, AJ, 152, 210

Tegler, S. C., Romanishin, W., & Consolmagno, S. J. G. J. 2003, ApJL, 599, L49

Thébault, P. 2003, EM&P, 92, 233

Thébault, P., & Doressoundiram, A. 2003, Icar, 162, 27

Tody, D. 1986, Proc. SPIE, 627, 733
Trujillo, C. A., & Brown, M. E. 2002, ApJ, 566, L125

Trujillo, C. A., Jewitt, D. C., & Luu, J. X. 2001, AJ, 122, 457

van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, CSE, 13, 22

Vilenius, E., Kiss, C., Müller, T., et al. 2014, A&A, 564, A35

Volk, K., & Malhotra, R. 2013, Icar, 224, 66

Wong, I., & Brown, M. E. 2017, AJ, 153, 145

Zwillinger, D., & Kokoska, S. 2000, CRC Standard Probability and Statistics Tables and Formulae (New York: Chapman and Hall)