Source author record

Kathryn Volk

Kathryn Volk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

Col-OSSOS: The Two Types of Kuiper Belt Surfaces

The Colours of the Outer Solar System Origins Survey (Col-OSSOS) has gathered high quality, near-simultaneous (g-r) and (r-J) colours of 92 Kuiper Belt Objects (KBOs) with (u-g) and (r-z) gathered for some. We present the current state of the survey and data analysis. Recognizing that the optical colours of most icy bodies broadly follow the reddening curve, we present a new projection of the optical-NIR colours, which rectifies the main non-linear features in the optical-NIR along the ordinates. We find evidence for a bifurcation in the projected colours which presents itself as a diagonal empty region in the optical-NIR. A reanalysis of past colour surveys reveals the same bifurcation. We interpret this as evidence for two separate surface classes: the BrightIR class spans the full range of optical colours and broadly follows the reddening curve, while the FaintIR objects are limited in optical colour, and are less bright in the NIR than the BrightIR objects. We present a two class model. Objects in each class consist of a mix of separate blue and red materials, and span a broad range in colour. Spectra are modelled as linear optical and NIR spectra with different slopes, that intersect at some transition wavelength. The underlying spectral properties of the two classes fully reproduce the observed structures in the UV-optical-NIR colour space ($0.4\lesssimλ\lesssim1.4 \mbox{ $μ$m}$), including the bifurcation observed in the Col-OSSOS and H/WTSOSS datasets, the tendency for cold classical KBOs to have lower (r-z) colours than excited objects, and the well known bimodal optical colour distribution.

preprint2022arXiv

Dynamical Population of Comet Reservoirs

The Oort cloud and the scattered disk are the two primary reservoirs for long-period and short-period comets, respectively. In this review, we assess the known observational constraints on these reservoirs' properties and their formation. In addition, we discuss how the early orbital evolution of the giant planets generated the modern scattered disk from the early, massive planetesimal disk and how $\sim$5\% of this material was captured into the Oort cloud. We review how the Sun's birth environment and dynamical history within the Milky Way alters the formation and modern structure of the Oort cloud. Finally, we assess how the coming decade's anticipated observing campaigns may provide new insights into the formation and properties of the Oort cloud and scattered disk.

preprint2022arXiv

Free Inclinations for Transneptunian Objects in the Main Kuiper Belt

There is a complex inclination structure present in the transneptunian object (TNO) orbital distribution in the main classical belt region (between orbital semimajor axes of 39 and 48 au). The long-term gravitational effects of the giant planets make TNO orbits precess, but non-resonant objects maintain a nearly constant 'free' inclination ($I_\text{free}$) with respect to a local forced precession pole. Because of the likely cosmogonic importance of the distribution of this quantity, we tabulate free inclinations for all main-belt TNOs, each individually computed using barycentric orbital elements with respect to each object's local forcing pole. We show that the simplest method, based on the Laplace-Lagrange secular theory, is unable to give correct forcing poles for objects near the $ν_{18}$ secular resonance, resulting in poorly conserved $I_\text{free}$ values in much of the main belt. We thus instead implemented an averaged Hamiltonian to obtain the expected nodal precession for each TNO, yielding significantly more accurate free inclinations for non-resonant objects. For the vast majority (96\%) of classical belt TNOs, these $I_\text{free}$ values are conserved to $<1^\circ$ over 4 Gyr numerical simulations, demonstrating the advantage of using this well-conserved quantity in studies of the TNO population and its primordial inclination profile; our computed distributions only reinforce the idea of a very co-planar surviving 'cold' primordial population, overlain by a large $I$-width implanted 'hot' population.

preprint2021arXiv

Contemporaneous Multi-Wavelength and Precovery Observations of Active Centaur P/2019 LD2 (ATLAS

Gateway Centaur and Jupiter co-orbital P/2019 LD2 (ATLAS) (Sarid et al. 2019) provides the first opportunity to observe the migration of a Solar System small body from a Centaur orbit to a Jupiter Family Comet (JFC) four decades from now (Kareta et al., 2020; Hsieh et al. 2020). The Gateway transition region is beyond where water ice can power cometary activity, and coma production there is as poorly understood as in all Centaurs. We present contemporaneous multi-wavelength observations of LD2 from 2020 July 2-4: Gemini-North visible imaging, NASA IRTF near-infrared spectroscopy, and ARO SMT millimeter-wavelength spectroscopy. Precovery DECam images limit the nucleu's effective radius to <=1.2 km and no large outbursts were seen in archival Catalina Sky Survey observations. LD2's coma has g'-r'=0.70+/-0.07, r'-i'=0.26+/-0.07, a dust production rate of ~10-20 kg/s, and an outflow velocity between v~0.6-3.3 m/s. We did not detect CO towards LD2 on 2020 July 2-3, with a 3-sigma upper limit of Q(CO) < 4.4 * 10^27 mol/s (<200 kg/s). Near-infrared spectra show evidence for water ice at the 1-10% level depending on grain size. Spatial profiles and archival data are consistent with sustained activity. The evidence supports the hypothesis that LD2 is a typical small Centaur that will become a typical JFC, and thus it is critical to understanding the transition between these two populations. Finally, we discuss potential strategies for a community-wide, long baseline monitoring effort.

preprint2021arXiv

OSSOS: The Eccentricity and Inclination Distributions of the Stable Neptunian Trojans

The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $σ_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $σ_e$ of 0.045 and a typical sin(i) x Gaussian inclination distribution with a width $σ_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $σ_{e-cold}$ ~ 0.02 / $σ_{i-cold}$ ~ 6$^\circ$ and $σ_{e-hot}$~ 0.05 / $σ_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 members and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range.

preprint2020arXiv

Col-OSSOS: Compositional homogeneity of three Kuiper belt binaries

The surface characterization of Trans-Neptunian Binaries (TNBs) is key to understanding the properties of the disk of planetesimals from which these objects formed. In the optical wavelengths, it has been demonstrated that most equal-sized component systems share similar colors, suggesting they have a similar composition. The color homogeneity of binary pairs contrasts with the overall diversity of colors in the Kuiper belt, which was interpreted as evidence that Trans-Neptunian Objects (TNOs) formed from a locally homogeneous and globally heterogeneous protoplanetary disk. In this paradigm, binary pairs must have formed early, before the dynamically hot TNOs were scattered out from their formation location. The latter inferences, however, relied on the assumption that the matching colors of the binary components imply matching composition. Here, we test this assumption by examining the component-resolved photometry of three TNBs found in the Outer Solar System Origins Survey: 505447 (2013 SQ99), 511551 (2014 UD225) and 506121 (2016 BP81), across the visible and J-band near-infrared wavelength range. We report similar colors within 2 sigma for the binary pairs suggestive of similar reflectance spectra and hence surface composition. This advocates for gravitational collapse of pebble clouds as a possible TNO formation route. We however stress that several similarly small TNOs, including at least one binary, have been shown to exhibit substantial spectral variability in the near-infrared, implying color equality of binary pairs is likely to be violated in some cases.

preprint2020arXiv

Dynamical instabilities in systems of multiple short-period planets are likely driven by secular chaos: a case study of Kepler-102

We investigated the dynamical stability of high-multiplicity Kepler and K2 planetary systems. Our numerical simulations find instabilities in $\sim20\%$ of the cases on a wide range of timescales (up to $5\times10^9$ orbits) and over an unexpectedly wide range of initial dynamical spacings. To identify the triggers of long-term instability in multi-planet systems, we investigated in detail the five-planet Kepler-102 system. Despite having several near-resonant period ratios, we find that mean motion resonances are unlikely to directly cause instability for plausible planet masses in this system. Instead, we find strong evidence that slow inward transfer of angular momentum deficit (AMD) via secular chaos excites the eccentricity of the innermost planet, Kepler-102 b, eventually leading to planet-planet collisions in $\sim80\%$ of Kepler-102 simulations. Kepler-102 b likely has a mass $>\sim0.1M_{\oplus}$, hence a bulk density exceeding about half Earth's, in order to avoid dynamical instability. To investigate the role of secular chaos in our wider set of simulations, we characterize each planetary system's AMD evolution with a "spectral fraction" calculated from the power spectrum of short integrations ($\sim5\times10^6$ orbits). We find that small spectral fractions ($\lesssim0.01$) are strongly associated with dynamical stability on long timescales ($5\times10^9$ orbits) and that the median time to instability decreases with increasing spectral fraction. Our results support the hypothesis that secular chaos is the driver of instabilities in many non-resonant multi-planet systems, and also demonstrate that the spectral analysis method is an efficient numerical tool to diagnose long term (in)stability of multi-planet systems from short simulations.

preprint2020arXiv

Machine Learning Classification of Kuiper Belt Populations

In the outer solar system, the Kuiper Belt contains dynamical sub-populations sculpted by a combination of planet formation and migration and gravitational perturbations from the present-day giant planet configuration. The subdivision of observed Kuiper Belt objects (KBOs) into different dynamical classes is based on their current orbital evolution in numerical integrations of their orbits. Here we demonstrate that machine learning algorithms are a promising tool for reducing both the computational time and human effort required for this classification. Using a Gradient Boosting Classifier, a type of machine learning regression tree classifier trained on features derived from short numerical simulations, we sort observed KBOs into four broad, dynamically distinct populations - classical, resonant, detached, and scattering - with a >97 per cent accuracy for the testing set of 542 securely classified KBOs. Over 80 per cent of these objects have a $>3σ$ probability of class membership, indicating that the machine learning method is classifying based on the fundamental dynamical features of each population. We also demonstrate how, by using computational savings over traditional methods, we can quickly derive a distribution of class membership by examining an ensemble of object clones drawn from the observational errors. We find two major reasons for misclassification: inherent ambiguity in the orbit of the object - for instance, an object that is on the edge of resonance - and a lack of representative examples in the training set. This work provides a promising avenue to explore for fast and accurate classification of the thousands of new KBOs expected to be found by surveys in the coming decade.

preprint2020arXiv

OSSOS XX: The Meaning of Kuiper Belt Colors

Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.

preprint2019arXiv

Carbon Chain Depletion of 2I/Borisov

The composition of comets in the Solar System come in multiple groups thought to encode information about their formation in different regions of the outer protosolar disk. The recent discovery of the second interstellar object, 2I/Borisov, allows for spectroscopic investigations into its gas content and a preliminary classification of it within the Solar System comet taxonomies to test the applicability of planetesimal formation models to other stellar systems. We present spectroscopic and imaging observations from 2019 September 20th to October 26th at the Bok, MMT, and LBT telescopes. We identify CN in the comet's spectrum and set precise upper limits on the abundance of C2 on all dates. We use a Haser model to convert our integrated fluxes to production rates and find Q(CN) = 5.0 +/- 2.0 * 10^24 mol/s on September 20th and Q(CN) = 1.1 - 1.9 * 10^24 mol/s on later dates, both consistent with contemporaneous observations. We set our lowest upper limit on a C2 production rate, Q(C2) < 1.6 * 10^23 mol/s, on October 10th. The measured ratio upper limit for that date, Q(C2)/Q(CN) < 0.095 indicates that 2I/Borisov is strongly in the (carbon chain) 'depleted' taxonomic group. The only comparable Solar System comets have detected ratios near this limit, making 2I/Borisov statistically likely to be more depleted than any known comet. Most 'depleted' comets are Jupiter Family Comets, perhaps indicating a similiarity in formation conditions between the most depleted of the JFCs and 2I/Borisov. More work is needed to understand the applicability of our knowledge of Solar System comet taxonomies onto interstellar objects, and we discuss future work that could help clarify the usefulness of the approach.

preprint2019arXiv

OSSOS XVIII: Constraining migration models with the 2:1 resonance using the Outer Solar System Origins Survey

Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $σ_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $α=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history.

preprint2019arXiv

Physical Characterization of the December 2017 Outburst of the Centaur 174P/Echeclus

The Centaurs are the small solar system bodies intermediate between the active inner solar system Jupiter Family Comets and their inactive progenitors in the trans-Neptunian region. Among the fraction of Centaurs which show comet-like activity, 174P/Echeclus is best known for its massive 2005 outburst in which a large apparently active fragment was ejected above the escape velocity from the primary nucleus. We present visible imaging and near-infrared spectroscopy of Echeclus during the first week after its December 2017 outburst taken at the Faulkes North & South Telescopes and the NASA IRTF, the largest outburst since 2005. The coma was seen to be highly asymmetric. A secondary peak was seen in the near-infrared 2D spectra, which is strongly hinted at in the visible images, moving hyperbolically with respect to the nucleus. The retrieved reflectance spectrum of Echelcus is consistent with the unobscured nucleus but becomes bluer when a wider extraction aperture is used. We find that Echeclus's coma is best explained as dominated by large blue dust grains, which agrees with previous work. We also conducted a high-resolution orbital integration of Echeclus's recent evolution and found no large orbital changes that could drive its modern evolution. We interpret the second peak in the visible and near-infrared datasets as a large cloud of larger-than-dust debris ejected at the time of outburst. If Echeclus is typical of the Centaurs, there may be several debris ejection or fragmentation events per year on other Centaurs that are going unnoticed.

preprint2016arXiv

Corralling a distant planet with extreme resonant Kuiper belt objects

The four longest period Kuiper belt objects have orbital periods close to integer ratios with each other. A hypothetical planet with orbital period $\sim$17,117 years, semimajor axis $\sim$665 AU, would have N/1 and N/2 period ratios with these four objects. The orbital geometries and dynamics of resonant orbits constrain the orbital plane, the orbital eccentricity and the mass of such a planet, as well as its current location in its orbital path.

preprint2016arXiv

OSSOS III - Resonant Trans-Neptunian Populations: Constraints from the first quarter of the Outer Solar System Origins Survey

The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic Plane Survey (CFEPS; Gladman et al. 2012), with which we find broad agreement. We confirm that the 5:2 resonance is more populated than models of the outer Solar System's dynamical history predict; our minimum population estimate shows that the high eccentricity (e>0.35) portion of the resonance is at least as populous as the 2:1 and possibly as populated as the 3:2 resonance. One OSSOS block was well-suited to detecting objects trapped at low libration amplitudes in Neptune's 3:2 resonance, a population of interest in testing the origins of resonant TNOs. We detected three 3:2 objects with libration amplitudes below the cutoff modeled by CFEPS; OSSOS thus offers new constraints on this distribution. The OSSOS detections confirm that the 2:1 resonance has a dynamically colder inclination distribution than either the 3:2 or 5:2 resonances. Using the combined OSSOS and CFEPS 2:1 detections, we constrain the fraction of 2:1 objects in the symmetric mode of libration to be 0.2-0.85; we also constrain the fraction of leading vs. trailing asymmetric librators, which has been theoretically predicted to vary depending on Neptune's migration history, to be 0.05-0.8. Future OSSOS blocks will improve these constraints.

preprint2016arXiv

OSSOS: IV. Discovery of a dwarf planet candidate in the 9:2 resonance with Neptune

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

preprint2016arXiv

The Outer Solar System Origins Survey: I. Design and First-Quarter Discoveries

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

preprint2015arXiv

Consolidating and Crushing Exoplanets: Did it happen here?

The Kepler mission results indicate that systems of tightly-packed inner planets (STIPs) are present around of order 5% of FGK field stars (whose median age is ~5 Gyr). We propose that STIPs initially surrounded nearly all such stars and those observed are the final survivors of a process in which long-term metastability eventually ceases and the systems proceed to collisional consolidation or destruction, losing roughly equal fractions of systems every decade in time. In this context, we also propose that our Solar System initially contained additional large planets interior to the current orbit of Venus, which survived in a metastable dynamical configuration for 1-10% of the Solar System's age. Long-term gravitational perturbations caused the system to orbit cross, leading to a cataclysmic event which left Mercury as the sole surviving relic.

preprint2013arXiv

Do Centaurs preserve their source inclinations?

The Centaurs are a population of small, planet-crossing objects in the outer solar system. They are dynamically short-lived and represent the transition population between the Kuiper belt and the Jupiter family short-period comets. Dynamical models and observations of the physical properties of the Centaurs indicate that they may have multiple source populations in the trans-Neptunian region. It has been suggested that the inclination distribution of the Centaurs may be useful in distinguishing amongst these source regions. The Centaurs, however, undergo many close encounters with the giant planets during their orbital evolution; here we show that these encounters can substantially determine the inclination distribution of the Centaurs. Almost any plausible initial inclination distribution of a Kuiper belt source results in Centaurs having inclinations peaked near 10-20 degrees. Our studies also find that the Kuiper belt is an extremely unlikely source of the retrograde Centaur that has been observed.

preprint2012arXiv

The effect of orbital evolution on the Haumea (2003 EL61) collisional family

The Haumea family is currently the only identified collisional family in the Kuiper belt. We numerically simulate the long-term dynamical evolution of the family to estimate a lower limit of the family's age and to assess how the population of the family and its dynamical clustering are preserved over Gyr timescales. We find that the family is not younger than 100 Myr, and its age is at least 1 Gyr with 95% confidence. We find that for initial velocity dispersions of 50-400 m/s, approximately 20-45% of the family members are lost to close encounters with Neptune after 3.5 Gyr of orbital evolution. We apply these loss rates to two proposed models for the formation of the Haumea family, a graze-and-merge type collision between two similarly sized, differentiated KBOs or the collisional disruption of a satellite orbiting Haumea. For the graze-and-merge collision model, we calculate that >85% of the expected mass in surviving family members within 150 m/s of the collision has been identified, but that one to two times the mass of the known family members remains to be identified at larger velocities. For the satellite-break-up model, we estimate that the currently identified family members account for ~50% of the expected mass of the family. Taking observational incompleteness into account, the observed number of Haumea family members is consistent with either formation scenario at the 1 sigma level, however both models predict more objects at larger relative velocities (>150 m/s) than have been identified.

preprint2011arXiv

Inclination Mixing in the Classical Kuiper Belt

We investigate the long-term evolution of the inclinations of the known classical and resonant Kuiper belt objects (KBOs). This is partially motivated by the observed bimodal inclination distribution and by the putative physical differences between the low- and high-inclination populations. We find that some classical KBOs undergo large changes in inclination over gigayear timescales, which means that a current member of the low-inclination population may have been in the high-inclination population in the past, and vice versa. The dynamical mechanisms responsible for the time-variability of inclinations are predominantly distant encounters with Neptune and chaotic diffusion near the boundaries of mean motion resonances. We reassess the correlations between inclination and physical properties including inclination time-variability. We find that the size-inclination and color-inclination correlations are less statistically significant than previously reported (mostly due to the increased size of the data set since previous works with some contribution from inclination variability). The time-variability of inclinations does not change the previous finding that binary classical KBOs have lower inclinations than non-binary objects. Our study of resonant objects in the classical Kuiper belt region includes objects in the 3:2, 7:4, 2:1, and eight higher-order mean motion resonances. We find that these objects (some of which were previously classified as non-resonant) undergo larger changes in inclination compared to the non-resonant population, indicating that their current inclinations are not generally representative of their original inclinations. They are also less stable on gigayear timescales.