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 $\sigma_{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 $\sigma_e$ of 0.045 and a typical sin(i) x gaussian inclination distribution with a width $\sigma_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 $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 member 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.

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.

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.

The flyby of Pluto and Charon by the New Horizons spacecraft provided high-resolution images of cratered surfaces embedded in the Kuiper belt, an extensive region of bodies orbiting beyond Neptune. Impact craters on Pluto and Charon were formed by collisions with other Kuiper belt objects (KBOs) with diameters from ~40 kilometers to ~300 meters, smaller than most KBOs observed directly by telescopes. We find a relative paucity of small craters less than approximately 13 kilometers in diameter, which cannot be explained solely by geological resurfacing. This implies a deficit of small KBOs (less than 1 to 2 kilometers in diameter). Some surfaces on Pluto and Charon are likely greater than 4 billion years old, thus their crater records provide information on the size-frequency distribution of KBOs in the early Solar System.

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 $\sigma_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 $\alpha=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.

There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Our top recommendation is to join LSST and gain access to the full data releases by hosting a data centre, which could be done by adding to the CADC, which is already highly involved in hosting planetary data and supporting computational modelling for orbital studies. We also support MSE, which can provide spectroscopy and thus compositional information for thousands of small bodies. We support a Canadian-led microsatellite, POEP, which will provide small body sizes by measuring occultations. We support the idea of piggybacking space weather instruments on other astronomical space probes to provide data for the space weather community. Many Canadian planetary scientists are involved in space probe missions, but through haphazard and temporary arrangements like co-appointments at US institutions, so we would like the community to support Canadian researchers to participate in these large, international missions.

In this white paper (WP), we highlight several examples of small and moderate aperture telescopes that are being used for education and/or research. We further discuss potential costs for establishing new, small observatories, as well as joining existing international consortia. The WP includes a brief overview of select observing sites, with a discussion on how small telescopes at exceptional observing locations can be competitive, under certain circumstances, with larger and more expensive facilities located at poorer sites. In addition to research, these facilities enable many different types of educational experiences for wide range of people, from high school students to undergraduates to graduate students to postdocs. Canada should remain committed to partnering with large, international observatories such as CFHT, Gemini, and TMT, but it should also negotiate international agreements and commit funding to expand the use of small and moderate research observatories at domestic and international sites through coordination with the NRC, the Tri-Council, and the Canadian Foundation for Innovation. Both capital and operational costs (with site rental costs allowed) need to be included in support possibilities. CASCA should establish and maintain a small to moderate telescope expression of interest database that would help to facilitate Canadian institutions in organizing consortia, particularly for smaller institutions. The astronomical community should work with the NRC to make existing facilities more accessible to the astronomical community for research. This could involve, for example, automating the Plaskett and/or providing travel funds for supporting classical observing modes. Finally, a small to moderate aperture facility in the Arctic would be a world-class observatory and should be advanced over the next decade.

[Highly abridged, from executive summary] As much as NewSpace presents opportunities, there are significant challenges that must be overcome, requiring engagement with policy makers to influence domestic and international space governance. Failure to do so could result in a range of long-lasting negative outcomes for science and space stewardship. How will the Canadian astronomical community engage with NewSpace? What are the implications for NewSpace on the astro-environment, including Earth orbits, lunar and cis-lunar orbits, and surfaces of celestial bodies? This white paper analyzes the rapid changes in space use and what those changes could mean for Canadian astronomers. Our recommendations are as follows: Greater cooperation between the astronomical and the Space Situational Awareness communities is needed. Build closer ties between the astronomical community and Global Affairs Canada (GAC). Establish a committee for evaluating the astro-environmental impacts of human space use, including on and around the Moon and other bodies. CASCA and the Tri-Council should coordinate to identify programs that would enable Canadian astronomers to participate in pay-for-use services at appropriate funding levels. CASCA should continue to foster a relationship with CSA, but also build close ties to the private space industry. Canadian-led deep space missions are within Canada's capabilities, and should be pursued.

While collisional families are common in the asteroid belt, only one is known in the Kuiper belt, linked to the dwarf planet Haumea. The characterization of Haumea's family helps to constrain its origin and, more generally, the collisional history of the Kuiper belt. However, the size distribution of the Haumea family is difficult to constrain from the known sample, which is affected by discovery biases. Here, we use the Outer Solar System Origins Survey (OSSOS) Ensemble to look for Haumea family members. In this OSSOS XVI study we report the detection of three candidates with small ejection velocities relative to the family formation centre. The largest discovery, 2013 UQ15, is conclusively a Haumea family member, with a low ejection velocity and neutral surface colours. Although the OSSOS Ensemble is sensitive to Haumea family members to a limiting absolute magnitude (Hr) of 9.5 (inferred diameter of ~90 km), the smallest candidate is significantly larger, Hr = 7.9. The Haumea family members larger than approximately 20 km in diameter must be characterized by a shallow H-distribution slope in order to produce only these three large detections. This shallow size distribution suggests that the family formed in a graze-and-merge scenario, not a catastrophic collision.

We use published models of the early Solar System evolution with a slow, long-range and grainy migration of Neptune to predict the orbital element distributions and the number of modern-day Centaurs. The model distributions are biased by the OSSOS survey simulator and compared with the OSSOS Centaur detections. We find an excellent match to the observed orbital distribution, including the wide range of orbital inclinations which was the most troublesome characteristic to fit in previous models. A dynamical model, in which the original population of outer disk planetesimals was calibrated from Jupiter Trojans, is used to predict that OSSOS should detect 11+/-4 Centaurs with semimajor axis a<30 au, perihelion distance q>7.5 au and diameter D>10 km (absolute magnitude H_r<13.7 for a 6% albedo). This is consistent with 15 actual OSSOS Centaur detections with H_r<13.7. The population of Centaurs is estimated to be 21,000+/-8,000 for D>10 km. The inner scattered disk at 50<a<200 au should contain (2.0+/-0.8)x10^7 D>10 km bodies and the Oort cloud should contain (5.0+/-1.9)x10^8 D>10 km comets. Population estimates for different diameter cutoffs can be obtained from the size distribution of Jupiter Trojans (N(>D) proportional to D^(-2.1) for 5<D<100 km). We discuss model predictions for the Large Synoptic Survey Telescope observations of Centaurs.

Most known trans-Neptunian objects (TNOs) gravitationally scattering off the giant planets have orbital inclinations consistent with an origin from the classical Kuiper belt, but a small fraction of these "scattering TNOs" have inclinations that are far too large (i > 45 deg) for this origin. These scattering outliers have previously been proposed to be interlopers from the Oort cloud or evidence of an undiscovered planet. Here we test these hypotheses using N-body simulations and the 69 centaurs and scattering TNOs detected in the Outer Solar Systems Origins Survey and its predecessors. We confirm that observed scattering objects cannot solely originate from the classical Kuiper belt, and we show that both the Oort cloud and a distant planet generate observable highly inclined scatterers. Although the number of highly inclined scatterers from the Oort Cloud is ~3 times less than observed, Oort cloud enrichment from the Sun's galactic migration or birth cluster could resolve this. Meanwhile, a distant, low-eccentricity 5 Earth-mass planet replicates the observed fraction of highly inclined scatterers, but the overall inclination distribution is more excited than observed. Furthermore, the distant planet generates a longitudinal asymmetry among detached TNOs that is less extreme than often presumed, and its direction reverses across the perihelion range spanned by known TNOs. More complete models that explore the dynamical origins of the planet are necessary to further study these features. With observational biases well-characterized, our work shows that the orbital distribution of detected scattering bodies is a powerful constraint on the unobserved distant solar system.

The migration of Neptune's resonances through the proto-Kuiper belt has been imprinted in the distribution of small bodies in the outer Solar System. Here we analyze five published Neptune migration models in detail, focusing on the high pericenter distance (high-q) trans-Neptunian Objects (TNOs) near Neptune's 5:2 and 3:1 mean-motion resonances, because they have large resonant populations, are outside the main classical belt, and are relatively isolated from other strong resonances. We compare the observationally biased output from these dynamical models with the detected TNOs from the Outer Solar System Origins Survey, via its Survey Simulator. All of the four new OSSOS detections of high-q non-resonant TNOs are on the Sunward side of the 5:2 and 3:1 resonances. We show that even after accounting for observation biases, this asymmetric distribution cannot be drawn from a uniform distribution of TNOs at 2sigma confidence. As shown by previous work, our analysis here tentatively confirms that the dynamical model that uses grainy slow Neptune migration provides the best match to the real high-q TNO orbital data. However, due to extreme observational biases, we have very few high-q TNO discoveries with which to statistically constrain the models. Thus, this analysis provides a framework for future comparison between the output from detailed, dynamically classified Neptune migration simulations and the TNO discoveries from future well-characterized surveys. We show that a deeper survey (to a limiting r-magnitude of 26.0) with a similar survey area to OSSOS could statistically distinguish between these five Neptune migration models.

The Colours of the Outer Solar System Origins Survey (Col-OSSOS) is acquiring near-simultaneous $g$, $r$, and $J$ photometry of unprecedented precision with the Gemini North Telescope, targeting nearly a hundred trans-Neptunian objects (TNOs) brighter than $m_r=23.6$ mag discovered in the Outer Solar System Origins Survey. Combining the optical and near-infrared photometry with the well-characterized detection efficiency of the Col-OSSOS target sample will provide the first flux-limited compositional dynamical map of the outer Solar System. In this paper, we describe our observing strategy and detail the data reduction processes we employ, including techniques to mitigate the impact of rotational variability. We present optical and near-infrared colors for 35 TNOs. We find two taxonomic groups for the dynamically excited TNOs, the neutral and red classes, which divide at $g-r \simeq 0.75$. Based on simple albedo and orbital distribution assumptions, we find that the neutral class outnumbers the red class, with a ratio of 4:1 and potentially as high as 11:1. Including in our analysis constraints from the cold classical objects, which are known to exhibit unique albedos and $r-z$ colors, we find that within our measurement uncertainty, our observations are consistent with the primordial Solar System protoplanetesimal disk being neutral-class-dominated, with two major compositional divisions in $grJ$ color space.

In preparation for the Jan 1/2019 encounter between the New Horizons spacecraft and the Kuiper Belt object 2014 MU69, we provide estimates of the expected impact crater surface density on the Kuiper Belt object. Using the observed crater fields on Charon and Pluto down to the resolution limit of the 2015 New Horizons flyby of those bodies and estimates of the orbital distribution of the crater-forming projectiles, we calculate the number of craters per unit area formed as a function of the time a surface on 2014 MU69 has been exposed to bombardment. We find that if the shallow crater-size distribution from 1-15 km exhibited on Pluto and Charon is indeed due to the sizes of Kuiper Belt projectiles, 2014 MU69 should exhibit a surface that is only lightly cratered below 1 km scale, despite being bombarded for $\sim4$~billion years. Its surface should therefore be more clearly indicative of its accretionary environment. In addition, this object may the first observed for which the majority of the bombardment is from exogenic projectiles moving less than or near the speed of sound in the target materials, implying morphologies more akin to secondary craters elsewhere in the Solar System. Lastly, if the shallow Kuiper Belt size distribution implied from the Pluto and Charon imaging is confirmed at MU69, then we conclude that this size distribution is a preserved relic of its state $\simeq$4.5~Gyr ago and provides a direct constraint on the planetesimal formation process itself.

We present variability measurements and partial light curves of Trans-Neptunian Objects (TNOs) from a two-night pilot study using Hyper Suprime-Cam (HSC) on the Subaru Telescope (Maunakea, Hawai'i, USA). Subaru's large aperture (8-m) and HSC's large field of view (1.77 square degrees) allow us to obtain measurements of multiple objects with a range of magnitudes in each telescope pointing. We observed 65 objects with m_r = 22.6--25.5 mag in just six pointings, allowing 20--24 visits of each pointing over the two nights. Our sample, all discovered in the recent Outer Solar System Origins Survey (OSSOS), span absolute magnitudes H_r = 6.2--10.8 mag and thus investigates smaller objects than previous light curve projects have typically studied. Our data supports the existence of a correlation between light curve amplitude and absolute magnitude seen in other works, but does not support a correlation between amplitude and orbital inclination. Our sample includes a number of objects from different dynamical populations within the trans-Neptunian region, but we do not find any relationship between variability and dynamical class. We were only able to estimate periods for 12 objects in the sample and found that a longer baseline of observations is required for reliable period analysis. We find that 31 objects (just under half of our sample) have variability greater than 0.4 magnitudes during all of the observations; in smaller 1.25 hr, 1.85 hr and 2.45 hr windows, the median variability is 0.13, 0.16 and 0.19 mags, respectively. The fact that variability on this scale is common for small TNOs has important implications for discovery surveys (such as OSSOS or the Large Synoptic Survey Telescope) and color measurements.

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 (Hr>5) dynamically excited Trans-Neptunian Objects (TNOs) and centaurs acquired as part of the Colours of the Outer Solar System Origins Survey (Col-OSSOS). We combine our dataset 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 (Tegler et al. 2016). 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.

Context. Centaurs are icy objects in transition between the transneptunian region and the inner solar system, orbiting the Sun in the giant planet region. Some Centaurs display cometary activity, which cannot be sustained by the sublimation of water ice in this part of the solar system, and has been hypothesized to be due to the crystallization of amorphous water ice. Aims. In this work, we look at Centaurs discovered by the Outer Solar System Origins Survey (OSSOS) and search for cometary activity. Tentative detections would improve understanding of the origins of activity among these objects. Methods. We search for comae and structures by fitting and subtracting both Point Spread Functions (PSF) and Trailed point-Spread Functions (TSF) from the OSSOS images of each Centaur. When available, Col-OSSOS images were used to search for comae too. Results. No cometary activity is detected in the OSSOS sample. We track the recent orbital evolution of each new Centaur to confirm that none would actually be predicted to be active, and we provide size estimates for the objects. Conclusions. The addition of 20 OSSOS objects to the population of 250 known Centaurs is consistent with the currently understood scenario, in which drastic drops in perihelion distance induce changes in the thermal balance prone to trigger cometary activity in the giant planet region.

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System.

We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100~au. The detection of these objects implies a 9:1 resonance population of $1.1\times10^4$ objects with $H_r<8.66$ ($D~\gtrsim~100$~km) on similar orbits (95\% confidence range of $\sim0.4-3\times10^4$). Integrations over 4~Gyr of an ensemble of clones spanning these objects' orbit fit uncertainties reveal that they both have median resonance occupation timescales of $\sim1$~Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last $\sim1$~Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000-4500 objects with $H_r<8.66$; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper belt's dynamical history.

All surveys include observational biases, which makes it impossible to directly compare properties of discovered trans-Neptunian Objects (TNOs) with dynamical models. However, by carefully keeping track of survey pointings on the sky, detection limits, tracking fractions, and rate cuts, the biases from a survey can be modelled in Survey Simulator software. A Survey Simulator takes an intrinsic orbital model (from, for example, the output of a dynamical Kuiper belt emplacement simulation) and applies the survey biases, so that the biased simulated objects can be directly compared with real discoveries. This methodology has been used with great success in the Outer Solar System Origins Survey (OSSOS) and its predecessor surveys. In this chapter, we give four examples of ways to use the OSSOS Survey Simulator to gain knowledge about the true structure of the Kuiper Belt. We demonstrate how to statistically compare different dynamical model outputs with real TNO discoveries, how to quantify detection biases within a TNO population, how to measure intrinsic population sizes, and how to use upper limits from non-detections. We hope this will provide a framework for dynamical modellers to statistically test the validity of their models.

The scattering trans-Neptunian Objects (TNOs) can be measured to smaller sizes than any other distant small-body population. We use the largest sample yet obtained, 68 discoveries, primarily by the Outer Solar System Origins Survey (OSSOS), to constrain the slope of its luminosity distribution, with sensitivity to much fainter absolute $H$ magnitudes than previous work. Using the analysis technique in Shankman et al. (2016), we confirm that a single slope for the $H$-distribution is not an accurate representation of the scattering TNOs and Centaurs, and that a break in the distribution is required, in support of previous conclusions. A bright-end slope of $\alpha_b=0.9$ transitioning to a faint-end slope $\alpha_f$ of 0.4-0.5 with a differential number contrast $c$ from 1 (a knee) to 10 (a divot) provides an acceptable match to our data. We find that break magnitudes $H_b$ of 7.7 and 8.3, values both previously suggested for dynamically hot Kuiper belt populations, are equally non-rejectable for a range of $\alpha_f$ and $c$ in our statistical analysis. Our preferred divot $H$-distribution transitions to $\alpha_f=0.5$ with a divot of contrast $c=3$ at $H_b=8.3$, while our preferred knee $H$-distribution transitions to $\alpha_f=0.4$ at $H_b=7.7$. The intrinsic population of scattering TNOs required to match the OSSOS detections is $3\times10^6$ for $H_r<12$, and $9\times10^4$ for $H_r<8.66$ ($D\gtrsim100$~km), with Centaurs having an intrinsic population two orders of magnitude smaller.

Several different classes of trans-Neptunian objects (TNOs) have been identified based on their optical and near-infrared colors. As part of the Colours of the Outer Solar System Origins Survey, we have obtained $g$, $r$, and $z$ band photometry of 26 TNOs using Subaru and Gemini Observatories. Previous color surveys have not utilized $z$ band reflectance, and the inclusion of this band reveals significant surface reflectance variations between sub-populations. The colors of TNOs in $g-r$ and $r-z$ show obvious structure, and appear consistent with the previously measured bi-modality in $g-r$. The distribution of colors of the two dynamically excited surface types can be modeled using the two-component mixing models from Fraser \& Brown (2012). With the combination of $g-r$ and $r-z$, the dynamically excited classes can be separated cleanly into red and neutral surface classes. In $g - r$ and $r - z$, the two dynamically excited surface groups are also clearly distinct from the cold classical TNO surfaces, which are red, with $g-r\gtrsim$0.85 and $r-z\lesssim$0.6, while all dynamically excited objects with similar $g-r$ colors exhibit redder $r-z$ colors. The $z$ band photometry makes it possible for the first time to differentiate the red excited TNO surfaces from the red cold classical TNO surfaces. The discovery of different $r-z$ colors for these cold classical TNOs makes it possible to search for cold classical surfaces in other regions of the Kuiper belt and to completely separate cold classical TNOs from the dynamically excited population, which overlaps in orbital parameter space.

The accumulating, but small, set of large semi-major axis trans-Neptunian objects (TNOs) shows an apparent clustering in the orientations of their orbits. This clustering must either be representative of the intrinsic distribution of these TNOs, or else arise as a result of observation biases and/or statistically expected variations for such a small set of detected objects. The clustered TNOs were detected across different and independent surveys, which has led to claims that the detections are therefore free of observational bias. This apparent clustering has led to the so-called "Planet 9" hypothesis that a super-Earth currently resides in the distant solar system and causes this clustering. The Outer Solar System Origins Survey (OSSOS) is a large program that ran on the Canada-France-Hawaii Telescope from 2013--2017, discovering more than 800 new TNOs. One of the primary design goals of OSSOS was the careful determination of observational biases that would manifest within the detected sample. We demonstrate the striking and non-intuitive biases that exist for the detection of TNOs with large semi-major axes. The eight large semi-major axis OSSOS detections are an independent dataset, of comparable size to the conglomerate samples used in previous studies. We conclude that the orbital distribution of the OSSOS sample is consistent with being detected from a uniform underlying angular distribution.

The cold classical Kuiper belt objects have low inclinations and eccentricities and are the only Kuiper belt population suspected to have formed in situ. Compared with the dynamically excited populations, which exhibit a broad range of colours and a low binary fraction of ~10% cold classical Kuiper belt objects typically have red optical colours with ~30% of the population found in binary pairs; the origin of these differences remains unclear. We report the detection of a population of blue-coloured, tenuously bound binaries residing among the cold classical Kuiper belt objects. Here we show that widely separated binaries could have survived push-out into the cold classical region during the early phases of Neptune's migration. The blue binaries may be contaminants, originating at ~38 au, and could provide a unique probe of the formative conditions in a region now nearly devoid of objects. The idea that the blue objects, which are predominantly binary, are the products of push-out requires that the planetesimals formed entirely as multiples. Plausible formation routes include planetesimal formation via pebble accretion and subsequent binary production through dynamic friction and binary formation during the collapse of a cloud of solids.

We report the discovery of the minor planet 2013 SY$_{99}$, on an exceptionally distant, highly eccentric orbit. With a perihelion of 50.0 au, 2013 SY$_{99}$'s orbit has a semi-major axis of $730 \pm 40$ au, the largest known for a high-perihelion trans-Neptunian object (TNO), well beyond those of (90377) Sedna and 2012 VP$_{113}$. Yet, with an aphelion of $1420 \pm 90$ au, 2013 SY$_{99}$'s orbit is interior to the region influenced by Galactic tides. Such TNOs are not thought to be produced in the current known planetary architecture of the Solar System, and they have informed the recent debate on the existence of a distant giant planet. Photometry from the Canada-France-Hawaii Telescope, Gemini North and Subaru indicate 2013 SY$_{99}$ is $\sim 250$ km in diameter and moderately red in colour, similar to other dynamically excited TNOs. Our dynamical simulations show that Neptune's weak influence during 2013 SY$_{99}$'s perihelia encounters drives diffusion in its semi-major axis of hundreds of astronomical units over 4 Gyr. The overall symmetry of random walks in semi-major axis allow diffusion to populate 2013 SY$_{99}$'s orbital parameter space from the 1000-2000 au inner fringe of the Oort cloud. Diffusion affects other known TNOs on orbits with perihelia of 45 to 49 au and semi-major axes beyond 250 au, providing a formation mechanism that implies an extended population, gently cycling into and returning from the inner fringe of the Oort cloud.

We present deep imaging observations, orbital dynamics, and dust tail model analyses of the double-component asteroid P/2016 J1 (J1-A and J1-B). The observations were acquired at the Gran Telescopio Canarias (GTC) and the Canada-France-Hawaii Telescope (CFHT) from mid March to late July, 2016. A statistical analysis of backward-in-time integrations of the orbits of a large sample of clone objects of P/2016 J1-A and J1-B shows that the minimum separation between them occurred most likely $\sim$2300 days prior to the current perihelion passage, i.e., during the previous orbit near perihelion. This closest approach was probably linked to a fragmentation event of their parent body. Monte Carlo dust tail models show that those two components became active simultaneously $\sim$250 days before the current perihelion, with comparable maximum loss rates of $\sim$0.7 kg s$^{-1}$ and $\sim$0.5 kg s$^{-1}$, and total ejected masses of 8$\times$10$^{6}$ kg and 6$\times$10$^{6}$ kg for fragments J1-A and J1-B, respectively. In consequence, the fragmentation event and the present dust activity are unrelated. The simultaneous activation times of the two components and the fact that the activity lasted 6 to 9 months or longer, strongly indicate ice sublimation as the most likely mechanism involved in the dust emission process.

The orbital element distribution of trans-Neptunian objects (TNOs) with large pericenters has been suggested to be influenced by the presence of an undetected, large planet at >200 AU from the Sun. To find additional observables caused by this scenario, we here present the first detailed emplacement simulation in the presence of a massive ninth planet on the distant Kuiper Belt. We perform 4 Gyr N-body simulations with the currently known Solar System planetary architecture, plus a 10 Earth mass planet with similar orbital parameters to those suggested by Trujillo & Sheppard (2014) or Batygin & Brown (2016), and 10^5 test particles in an initial planetesimal disk. We find that including a distant superearth-mass planet produces a substantially different orbital distribution for the scattering and detached TNOs, raising the pericenters and inclinations of moderate semimajor axis (50<a<500 AU) objects. We test whether this signature is detectable via a simulator with the observational characteristics of four precisely characterized TNO surveys. We find that the qualitatively very distinct Solar System models that include a ninth planet are essentially observationally indistinguishable from an outer Solar System produced solely by the four giant planets. We also find that the mass of the Kuiper Belt's current scattering and detached populations is required to be 3-10 times larger in the presence of an additional planet. We do not find any evidence for clustering of orbital angles in our simulated TNO population. Wide-field, deep surveys targeting inclined high-pericenter objects will be required to distinguish between these different scenarios.

We explore the distant giant planet hypothesis by integrating the large semi-major axis, large pericenter Trans-Neptunian Objects (TNOs) in the presence of the giant planets and an external perturber whose orbit is consistent with the proposed distant, eccentric, and inclined giant planet, so called planet 9. We find that TNOs with semi-major axes greater than 250 au experience some longitude of perihelion shepherding, but that a generic outcome of such evolutions is that the TNOs evolve to larger pericenter orbits, and commonly get raised to retrograde inclinations. This pericenter and inclination evolution requires a massive disk of TNOs (tens of M$_\Earth$) in order to explain the detection of the known sample today. Some of the highly inclined orbits produced by the examined perturbers will be inside of the orbital parameter space probed by prior surveys, implying a missing signature of the 9th planet scenario. The distant giant planet scenarios explored in this work do not reproduce the observed signal of simultaneous clustering in argument of pericenter, longitude of the ascending node, and longitude of perihelion in the region of the known TNOs.

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.

We report the orbital distribution of the Trans-Neptunian objects (TNOs) discovered during the High Ecliptic Latitude (HiLat) extension of the Canada-France Ecliptic Plane Survey (CFEPS), conducted from June 2006 to July 2009. The HiLat component was designed to address one of the shortcomings of ecliptic surveys (like CFEPS), their lack of sensitivity to high-inclination objects. We searched 701~deg$^2$ of sky ranging from 12$^\circ$ to 85$^\circ$ ecliptic latitude and discovered \lKBO TNOs, with inclinations between 15$^\circ$ to 104$^\circ$. This survey places a very strong constraint on the inclination distribution of the hot component of the classical Kuiper Belt, ruling out any possibility of a large intrinsic fraction of highly inclined orbits. Using the parameterization of \citet{2001AJ....121.2804B}, the HiLat sample combined with CFEPS imposes a width $14^\circ \le \sigma \le 15.5^\circ$, with a best match for $\sigma = 14.5^\circ$. HiLat discovered the first retrograde TNO, 2008~KV$_{42}$, with an almost polar orbit with inclination 104$^\circ$, and (418993), a scattering object with perihelion in the region of Saturn's influence, with $a \sim 400$~AU and $i = 68^\circ$.

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.

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.

We propose that two of the most surprising results so far among exoplanet discoveries are related: the existences of both hot Jupiters and the high frequency of multi-planet systems with periods $P\lesssim200$~days. In this paradigm, the vast majority of stars rapidly form along with multiple close-in planets in the mass range of Mars to super-Earths/mini-Neptunes. Such systems of tightly-packed inner planets (STIPs) are metastable, with the time scale of the dynamical instability having a major influence on final planet types. In most cases, the planets consolidate into a system of fewer, more massive planets, but long after the circumstellar gas disk has dissipated. This can yield planets with masses above the traditional critical core of $\sim$10 $M_\oplus$, yielding short-period giants that lack abundant gas. A rich variety of physical states are also possible given the range of collisional outcomes and formation time of the close-in planets. However, when dynamical consolidation occurs before gas dispersal, a critical core can form that then grows via gas capture into a short-period gas giant. In this picture the majority of Hot and Warm Jupiters formed locally, rather than migrating down from larger distances.

We measure the absolute magnitude, $H$, distribution, $dN(H) \propto 10^{\alpha H}$ of the scattering Trans-Neptunian Objects (TNOs) as a proxy for their size-frequency distribution. We show that the H-distribution of the scattering TNOs is not consistent with a single-slope distribution, but must transition around $H_g \sim 9$ to either a knee with a shallow slope or to a divot, which is a differential drop followed by second exponential distribution. Our analysis is based on a sample of 22 scattering TNOs drawn from three different TNO surveys, the Canada-France Ecliptic Plane Survey (CFEPS, Petit et al. 2011), Alexandersen et al. (2014), and the Outer Solar System Origins Survey (OSSOS, Bannister et al. 2016), all of which provide well characterized detection thresholds, combined with a cosmogonic model for the formation of the scattering TNO population. Our measured absolute magnitude distribution result is independent of the choice of cosmogonic model. Based on our analysis, we estimate that number of scattering TNOs is (2.4-8.3)$\times 10^5$ for $H_r < 12$. A divot $H$-distribution is seen in a variety of formation scenarios and may explain several puzzles in Kuiper Belt science. We find that a divot $H$-distribution simultaneously explains the observed scattering TNO, Neptune Trojan, Plutino, and Centaur $H$-distributions while simultaneously predicting a large enough scattering TNO population to act as the sole supply of the Jupiter-Family Comets.

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.

The Canada-France Ecliptic Plane Survey discovered four trans-Neptunian objects with semi-major axes near the 5:1 resonance, revealing a large and previously undetected intrinsic population. Three of these objects are currently resonant with Neptune, and the fourth is consistent with being an object that escaped the resonance at some point in the past. The non-resonant object may be representative of a detached population that is stable at slightly lower semi-major axes than the 5:1 resonance. We generated clones of these objects by resampling the astrometric uncertainty and examined their behavior over a 4.5 Gyr numerical simulation. The majority of the clones of the three resonant objects (>90%) spend a total of 10^7 years in resonance during their 4.5 Gyr integrations; most clones experience multiple periods of resonance capture. Our dynamical integrations reveal an exchange between the 5:1 resonance, the scattering objects, and other large semi-major axis resonances, especially the 4:1, 6:1, and 7:1. The multiple capture events and relatively short resonance lifetimes after capture suggest that these objects are captured scattering objects that stick in the 5:1 resonance. These 5:1 resonators may be representative of a temporary population, requiring regular contributions from a source population. We examined the dynamical characteristics (inclination, eccentricity, resonant island, libration amplitude) of the detected objects and their clones in order to provide an empirical model of the orbit structure of the 5:1 resonance. This resonance is dynamically hot and includes primarily symmetric librators. Given our orbit model, the intrinsic population necessary for the detection of these three objects in the 5:1 resonance is 1900(+3300 -1400, 95% confidence) with H_g< 8 and e > 0.5.

The planet candidate Fomalhaut b is bright in optical light but undetected in longer wavelengths, requiring a large, reflective dust cloud. The most recent observations find an extremely eccentric orbit (e ~ 0.8), indicating that Fomalhaut b cannot be the planet that is constraining the system's eccentric debris ring. An irregular satellite swarm around a super-Earth has been proposed, however, explaining the well-constrained debris ring requires an additional planet on an orbit that crosses that of the putative super-Earth. This paper expands upon a second theory: Fomalhaut b is a transient dust cloud produced by a catastrophic collision between planetesimals in the disk. We perform collisional probability simulations of the Fomalhaut debris disk based on the structure of our Kuiper belt, finding that the catastrophic disruption rate of d ~ 100 km bodies in the high-eccentricity scattering component is several per decade. This model paints a picture of the Fomalhaut system as having recently (within ~10-100 Myr) experienced a dynamical instability within its planetary system, which scattered a massive number of planetesimals onto large, high eccentricity orbits similar to that of Fom b. If Fomalhaut b is indeed a dust cloud produced by such a collision, we should soon see another appear, while Fomalhaut b will expand until it is either resolved or becomes too faint to be seen.

The Trans-Neptunian Objects (TNOs) may preserve evidence of planet building in their orbital and size-distributions. While all populations show steep size-distributions for large objects, recently relative deficit of Neptunian Trojans and scattering objects with diameters $D<100\,\mathrm{km}$ were detected. We have investigated this deficit with a 32 square degree survey, detecting 77 TNOs to a limiting $r$-band magnitude of 24.6. Our Plutinos sample (18 objects in 3:2 mean motion resonance with Neptune) also shows a deficit of $D<100\,\mathrm{km}$ objects. We reject a single power-law size-distribution and find that the Plutinos favour a divot. The Plutinos are thus added the list of populations with a deficit of $D<100\,\mathrm{km}$ objects. The fact that three independent samples of three different populations show this trend suggests that it is a real feature, possibly shared by all hot TNO populations as a remnant of "born big" planetesimal formation processes. We surmise the existence of $9000\pm3000$ Plutinos with $H_r\leq8.66$ and $37000^{+12000}_{-10000}$ Plutinos with $H_r\leq10.0$. Our survey also discovered one temporary Uranian Trojan, one temporary Neptunian Trojans and one stable Neptunian Trojan, from which we derive populations of $110^{+500}_{-100}$, $210^{+900}_{-200}$ and $150^{+600}_{-140}$ , respectively, with $H_r\leq10.0$. The Neptunian Trojans are thus less numerous than the main belt asteroids, which has over 700 asteroids with $H_r\leq10.0$. With such numbers, the temporary Neptunian Trojans cannot be previously stable Trojans that happen to be escaping the resonance now; they must be captured from another reservoir. With three 3:1 and one 4:1 resonators, we add to the growing evidence that the high-order resonances are more populated than expected.

The Hungaria asteroids are interior to the main asteroid belt, with semimajor axes between 1.8 and 2 AU, low eccentricities and inclinations of 16-35 degrees. Small asteroids in the Hungaria region are dominated by a collisional family associated with (434) Hungaria. The dominant spectral type of the Hungaria group is the E or X-type (Warner et al, 2009), mostly due to the E-type composition of Hungaria and its genetic family. It is widely believed the E-type asteroids are related to the aubrite meteorites, also known as enstatite achondrites (Gaffey et al, 1992). Here we explore the hypothesis that aubrites originate in the Hungaria family. In order to test this connection, we compare model Cosmic Ray Exposure ages from orbital integrations of model meteoroids with those of aubrites. We show that long CRE ages of aubrites (longest among stony meteorite groups) reflect the delivery route of meteoroids from Hungarias to Earth being different than those from main-belt asteroids. We find that the meteoroids from Hungarias predominantly reach Earth by Yarkovsky-drifting across the orbit of Mars, with no assistance from orbital resonances. We conclude that the CRE ages of aubrites are fully consistent with a dominant source at the inner boundary of the Hungaria family at 1.7 AU. From here, meteoroids reach Earth through the Mars-crossing region, with relatively quick delivery times favored due to collisions (with Hungarias and the inner main-belt objects). We find that, after Vesta, (434) Hungaria is the best candidate for an asteroidal source of an achondrite group.

Trojan objects share a planet's orbit, never straying far from the triangular Lagrangian points, 60 degrees ahead of (L4) or behind (L5) the planet. We report the detection of a Uranian Trojan; in our numerical integrations, 2011 QF99 oscillates around the Uranian L4 Lagrange point for >70 kyr and remains co-orbital for ~1 Myr before becoming a Centaur. We constructed a Centaur model, supplied from the transneptunian region, to estimate temporary co-orbital capture frequency and duration (to factor of two accuracy), finding that at any time 0.4% and 2.8% of the population will be Uranian and Neptunian co-orbitals, respectively. The co-orbital fraction (~2.4%) among Centaurs in the IAU Minor Planet Centre database is thus as expected under transneptunian supply.

Because of their proximity within the transneptunian region, the plutinos (objects in the 3:2 mean-motion resonance with Neptune) are numerous in flux-limited catalogs, and well-studied theoretically. We perform detailed modelling of the on-sky detection biases for plutinos, with special attention to those that are simultaneously in the Kozai resonance. In addition to the normal 3:2 resonant argument libration, Kozai plutinos also show periodic oscillations in eccentricity and inclination, coupled to the argument of perihelion (omega) oscillation. Due to the mean-motion resonance, plutinos avoid coming to pericenter near Neptune's current position in the ecliptic plane. Because Kozai plutinos are restricted to certain values of omega, perihelion always occurs out of the ecliptic plane, biasing ecliptic surveys against finding these objects. The observed Kozai plutino fraction (fkoz) has been measured by several surveys, finding values between 8% and 25%, while the true fkoz has been predicted to be between 10% and 30% by different giant planet migration simulations. We show that the observed fkoz varies widely depending on the ecliptic latitude and longitude of the survey, so debiasing to find the true ratio is complex. Even a survey that covers most or all of the sky will detect an apparent Kozai fraction that is different from the true fkoz. We present a map of the on-sky plutino Kozai fraction that would be detected by all-sky flux-limited surveys. This will be especially important for the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) and Large Synoptic Survey Telescope (LSST) projects, which may detect large numbers of plutinos as they sweep the sky. The Kozai fraction and the distribution of the orbital elements of Kozai plutinos may be a diagnostic of giant planet migration; future migration simulations should provide details on their resonant Kozai populations.

Via joint analysis of a calibrated telescopic survey, which found scattering Kuiper Belt objects, and models of their expected orbital distribution, we measure the form of the scattering object's size distribution. Ruling out a single power-law at greater than 99% confidence, we constrain the form of the size distribution and find that, surprisingly, our analysis favours a very sudden decrease (a divot) in the number distribution as diameters decrease below 100 km, with the number of smaller objects then rising again as expected via collisional equilibrium. Extrapolating at this collisional equilibrium slope produced enough kilometer-scale scattering objects to supply the nearby Jupiter-Family comets. Our interpretation is that this divot feature is a preserved relic of the size distribution made by planetesimal formation, now "frozen in" to portions of the Kuiper Belt sharing a "hot" orbital inclination distribution, explaining several puzzles in Kuiper Belt science. Additionally, we show that to match today's scattering-object inclination distribution, the supply source that was scattered outward must have already been vertically heated to of order 10 degrees.

The GJ 581 system has been amply studied since its discovery in 2005: the number of known planets in the system has increased and their orbital parameters are among the most precisely determined for radial velocity detected exoplanets. We have acquired MOST space-based photometry during 2007 and 2009, with the aims of measuring the stellar variability and searching for transits of GJ 581e, respectively. We quantify our sensitivity to shallow transit signals using Monte Carlo simulations, and perform a transit search within the 3$\sigma$ transit windows corresponding to both the circular and Keplerian orbit ephemerides. Our analysis rules out transits for a planet with an orbital period of 3.15 days (GJ 581 e) having a radius larger than 1.62 $R_{\oplus}$ (or a density lower than 2.39 g cm$^{-3}$ for an orbital inclination of 90$^{\circ}$) to 2$\sigma$ confidence. Thus, if the planet transits, we can exclude hydrogen, helium and water theoretical model compositions. The MOST photometry also allows us to rule out transits of GJ 581b within the Keplerian orbit-derived transit window for impact parameter values smaller than $\sim$0.4 and confirm previous results which exclude transits for this planet within the circular orbit-derived transit window, for all plausible interior compositions. We find that the stellar brightness of GJ 581 is stable to within 1%, a characteristic which is favourable to the development of life in the habitable zone of the system. In the 2009 photometry, we detect a stellar signal with a period of 5.586 $\pm$ 0.051 days, which is close to the orbital period of GJ 581b ($P=$5.37 days). However, further monitoring of the system is necessary to verify the nature of this variation.

The transneptunian objects (TNOs) trapped in mean-motion resonances with Neptune were likely emplaced there during planet migration late in the giant-planet formation process. We perform detailed modelling of the resonant objects detected in the Canada-France Ecliptic Plane Survey (CFEPS) in order to provide population estimates and, for some resonances, constrain the complex internal orbital element distribution. Detection biases play a critical role because phase relationships with Neptune make object discovery more likely at certain longitudes. This paper discusses the 3:2, 5:2, 2:1, 3:1, 5:1, 4:3, 5:3, 7:3, 5:4, and 7:4 mean-motion resonances, all of which had CFEPS detections, along with our upper limit on 1:1 Neptune Trojans (which is consistent with their small population estimated elsewhere). For the plutinos (TNOs in the 3:2 resonance) we refine the orbital element distribution given in Kavelaars et al. (2009) and show that steep H-magnitude distributions (N(H) proportional to 10aH, with a=0.8-0.9) are favoured in the range Hg=8-9, and confirm that this resonance does not share the inclination distribution of the classical Kuiper Belt. We give the first population estimate for the 5:2 resonance and find that, to within the uncertainties, the population is equal to that of the 3:2 (13,000 TNOs with Hg < 9.16), whereas the 2:1 population is smaller by a factor of 3-4 compared to the other two resonances. We also measure significant populations inhabiting the 4:3, 5:3, 7:3, 5:4, 7:4, 3:1, and 5:1 resonances, with Hg < 9.16 (D >100 km) populations in the thousands. We compare our intrinsic population and orbital-element distributions with several published models of resonant-TNO production; the most striking discrepancy is that resonances beyond the 2:1 are in reality more heavily populated than in published models.

The Kepler Mission recently identified 997 systems hosting candidate extrasolar planets, many of which are super-Earths. Realizing these planetary systems are candidates to host extrasolar asteroid belts, we use mid-infrared data from the Wide-field Infrared Survey Explorer (WISE) to search for emission from dust in these systems. We find excesses around eight stars, indicating the presence of warm to hot dust (~100-500 K), corresponding to orbital distances of 0.1-10 AU for these solar-type stars. The strongest detection, KOI 1099, demands ~500 K dust interior to the orbit of its exoplanet candidate. One star, KOI 904, may host very hot dust (~1200 K, corresponding to 0.02 AU). Although the fraction of these exoplanet-bearing stars with detectable warm excesses (~3%) is similar to that found by Spitzer surveys of solar-type field stars, the excesses detectable in the WISE data have much higher fractional luminosities (Ldust/L*) than most known debris disks, implying that the fraction with debris disks of comparable luminosity may actually be significantly higher. It is difficult to explain the presence of dust so close to the host stars, generally corresponding to dust rings at radii <0.3 AU; both the collisional and Poynting-Robertson drag timescales to remove dust from the system are hundreds of years or less at these distances. Assuming a steady-state for these systems implies large mass consumption rates with these short removal timescales, meaning that the dust production mechanism in these systems must almost certainly be episodic in nature.

The low-inclination component of the Classical Kuiper Belt is host to a population of extremely widely-separated binaries. These systems are similar to other Trans-Neptunian binaries (TNBs) in that the primary and secondary components of each system are of roughly equal size. We have performed an astrometric monitoring campaign of a sample of seven wide-separation, long-period TNBs and present the first-ever well-characterized mutual orbits for each system. The sample contains the most eccentric (2006 CH69, e=0.9) and the most widely-separated, weakly bound (2001 QW322, a/Rh~0.22) binary minor planets known, and also contains the system with lowest-measured mass of any TNB (2000 CF105, M~1.85E17 kg). Four systems orbit in a prograde sense, and three in a retrograde sense. They have a different mutual inclination distribution compared to all other TNBs, preferring low mutual-inclination orbits. These systems have geometric r-band albedos in the range of 0.09-0.3, consistent with radiometric albedo estimates for larger solitary low-inclination Classical Kuiper Belt objects, and we limit the plausible distribution of albedos in this region of the Kuiper Belt. We find that gravitational collapse binary formation models produce a similar orbital distribution to that currently observed, which along with a confluence of other factors supports formation of the cold Classical Kuiper Belt in situ through relatively rapid gravitational collapse rather than slow hierarchical accretion. We show that these binary systems are sensitive to disruption via collisions, and their existence suggests that the size distribution of TNOs at small sizes remains relatively shallow.

We report the orbital distribution of the trans-neptunian objects (TNOs) discovered during the Canada-France Ecliptic Plane Survey, whose discovery phase ran from early 2003 until early 2007. The follow-up observations started just after the first discoveries and extended until late 2009. We obtained characterized observations of 321 sq.deg. of sky to depths in the range g ~ 23.5--24.4 AB mag. We provide a database of 169 TNOs with high-precision dynamical classification and known discovery efficiency. Using this database, we find that the classical belt is a complex region with sub-structures that go beyond the usual splitting of inner (interior to 3:2 mean-motion resonance [MMR]), outer (exterior to 2:1 MMR), and main (in between). The main classical belt (a=40--47 AU) needs to be modeled with at least three components: the `hot' component with a wide inclination distribution and two `cold' components (stirred and kernel) with much narrower inclination distributions. The hot component must have a significantly shallower absolute magnitude (Hg) distribution than the other two components. With 95% confidence, there are 8000+1800-1600 objects in the main belt with Hg <= 8.0, of which 50% are from the hot component, 40% from the stirred component and 10% from the kernel; the hot component's fraction drops rapidly with increasing Hg. Because of this, the apparent population fractions depend on the depth and ecliptic latitude of a trans-neptunian survey. The stirred and kernel components are limited to only a portion of the main belt, while we find that the hot component is consistent with a smooth extension throughout the inner, main and outer regions of the classical belt; the inner and outer belts are consistent with containing only hot-component objects. The Hg <= 8.0 TNO population estimates are 400 for the inner belt and 10,000 for the outer belt within a factor of two.

We conducted a satellite search around the dwarf planet 1 Ceres using Hubble Space Telescope and ground-based Palomar data. No candidate objects were found orbiting Ceres in its entire stability region down to ~500km from the surface of Ceres. Assuming a satellite would have the same albedo as Ceres, which has a visual geometric albedo of 0.07-0.10, our detection limit is sensitive to satellites larger than 1-2 km in diameter.

We report here the discovery of an new near-equal mass Trans-Neptunian Binaries (TNBs) L5c02 and the and the putative detection of a second TNB (L4k12) among the year two and three detections of the Canada-France-Eclipic Plane Survey (CFEPS). These new binaries (internal designation L4k12 and L5c02) have moderate separations of 0.4" and 0.6" respectively. The follow-up observation confirmed the binarity of L5c02, but L4k12 are still lack of more followup observations. L4k12 has a heliocentric orbital inclination of ~ 35?degree, marking this system as having the highest heliocentric orbital inclination among known near-equal mass binaries. Both systems are members of the classical main Kuiper belt population. Based on the sample of objects searched we determine that the fraction of near-equal mass wide binaries with separations > 0.4" is 1.5% to 20% in the cold main classical Kuiper belt and, if our detection the binarity L4k12 holds, 3% to 43% in the hot main classical objects are binary. In this manuscript we describe our detection process, the sample of objects surveyed, our confirmation observations.

We have searched Microvariability and Oscillations of STars (MOST) satellite photometry obtained in 2004, 2005, and 2007 of the solar-type star HD 209458 for Trojan asteroid swarms dynamically coupled with the system's transiting "hot Jupiter" HD 209458b. Observations of the presence and nature of asteroids around other stars would provide unique constraints on migration models of exoplanetary systems. Our results set an upper limit on the optical depth of Trojans in the HD 209458 system that can be used to guide current and future searches of similar systems by upcoming missions. Using cross-correlation methods with artificial signals implanted in the data, we find that our detection limit corresponds to a relative Trojan transit depth of 1\times10-4, equivalent to ~1 lunar mass of asteroids, assuming power-law Trojan size distributions similar to Jupiter's Trojans in our solar system. We confirm with dynamical interpretations that some asteroids could have migrated inward with the planet to its current orbit at 0.045 AU, and that the Yarkovsky effect is ineffective at eliminating objects of > 1 m in size. However, using numerical models of collisional evolution we find that, due to high relative speeds in this confined Trojan environment, collisions destroy the vast majority of the asteroids in <10 Myr. Our modeling indicates that the best candidates to search for exoTrojan swarms in 1:1 mean resonance orbits with "hot Jupiters" are young systems (ages of about 1 Myr or less). Years of Kepler satellite monitoring of such a system could detect an asteroid swarm with a predicted transit depth of 3\times10-7.

Multiple impact basins formed on the Moon about 3.8 Gyr ago in what is known as the lunar cataclysm or late heavy bombardment. Many workers currently interpret the lunar cataclysm as an impact spike primarily caused by main-belt asteroids destabilized by delayed planetary migration. We show that morphologically fresh (class 1) craters on the lunar highlands were mostly formed during the brief tail of the cataclysm, as they have absolute crater number density similar to that of the Orientale basin and ejecta blanket. The connection between class 1 craters and the cataclysm is supported by the similarity of their size-frequency distribution to that of stratigraphically-identified Imbrian craters. Majority of lunar craters younger than the Imbrium basin (including class 1 craters) thus record the size-frequency distribution of the lunar cataclysm impactors. This distribution is much steeper than that of main-belt asteroids. We argue that the projectiles bombarding the Moon at the time of the cataclysm could not have been main-belt asteroids ejected by purely gravitational means.

We have performed an ecliptic imaging survey of the Kuiper belt with our deepest and widest field achieving a limiting flux of m(g') = 26.4, with a sky coverage of 3.0 square-degrees. This is the largest coverage of any other Kuiper belt survey to this depth. We detect 72 objects, two of which have been previously observed. We have improved the Bayesian maximum likelihood fitting technique presented in Gladman et al. (1998) to account for calibration and sky density variations and have used this to determine the luminosity function of the Kuiper belt. Combining our detections with previous surveys, we find the luminosity function is well represented by a single power-law with slope alpha = 0.65 +/- 0.05 and an on ecliptic sky density of 1 object per square-degree brighter than m(R)=23.42 +/- 0.13. Assuming constant albedos, this slope suggests a differential size-distribution slope of 4.25 +/- 0.25, which is steeper than the Dohnanyi slope of 3.5 expected if the belt is in a state of collisional equilibrium. We find no evidence for a roll-over or knee in the luminosity function and reject such models brightward of m(R) ~ 24.6.

We have examined the fate of impact ejecta liberated from the surface of Mercury due to impacts by comets or asteroids, in order to study (1) meteorite transfer to Earth, and (2) re-accumulation of an expelled mantle in giant-impact scenarios seeking to explain Mercury's large core. In the context of meteorite transfer, we note that Mercury's impact ejecta leave the planet's surface much faster (on average) than other planet's in the Solar System because it is the only planet where impact speeds routinely range from 5-20 times the planet's escape speed. Thus, a large fraction of mercurian ejecta may reach heliocentric orbit with speeds sufficiently high for Earth-crossing orbits to exist immediately after impact, resulting in larger fractions of the ejecta reaching Earth as meteorites. We calculate the delivery rate to Earth on a time scale of 30 Myr and show that several percent of the high-speed ejecta reach Earth (a factor of -3 less than typical launches from Mars); this is one to two orders of magnitude more efficient than previous estimates. Similar quantities of material reach Venus. These calculations also yield measurements of the re-accretion time scale of material ejected from Mercury in a putative giant impact (assuming gravity is dominant). For mercurian ejecta escaping the gravitational reach of the planet with excess speeds equal to Mercury's escape speed, about one third of ejecta re-accretes in as little as 2 Myr. Thus collisional stripping of a silicate proto-mercurian mantle can only work effectively if the liberated mantle material remains in small enough particles that radiation forces can drag them into the Sun on time scale of a few million years, or Mercury would simply re-accrete the material.

We have made a comprehensive transit search for exoplanets down to ~1.5 - 2 Earth radii in the HD 189733 system, based on 21-days of nearly uninterrupted broadband optical photometry obtained with the MOST (Microvariability & Oscillations of STars) satellite in 2006. We have searched these data for realistic limb-darkened transits from exoplanets other than the known hot Jupiter, HD 189733b, with periods ranging from about 0.4 days to one week. Monte Carlo statistical tests of the data with synthetic transits inserted into the data-set allow us to rule out additional close-in exoplanets with sizes ranging from about 0.15 - 0.31 RJ (Jupiter radii), or 1.7 - 3.5 RE (Earth radii) on orbits whose planes are near that of HD 189733b. These null results constrain theories that invoke lower-mass hot Super-Earth and hot Neptune planets in orbits similar to HD 189733b due to the inward migration of this hot Jupiter. This work also illustrates the feasibility of discovering smaller transiting planets around chromospherically active stars.

We have made a comprehensive transit search for exoplanets down to about 2 Earth radii in the HD 209458 system, based on nearly uninterrupted broadband optical photometry obtained with the MOST (Microvariability and Oscillations of Stars) satellite, spanning 14 days in 2004 and 44 days in 2005. We have searched these data for limb-darkened transits at periods other than that of the known giant planet, from about 0.5 days to 2 weeks. Monte Carlo statistical tests of the data with synthetic transits inserted allow us to rule out additional close-in exoplanets with sizes ranging from about 0.20-0.36 RJ (Jupiter radii), or 2.2-4.0 RE (Earth radii) on orbits whose planes are near that of HD 209458b. These null results constrain theories that invoke lower mass planets in orbits similar to HD 209458b to explain its anomalously large radius, and those that predict "hot Earths" due to the inward migration of HD 209458b.

We propose that the presence of additional planets in extrasolar planetary systems can be detected by long-term transit timing studies. If a transiting planet is on an eccentric orbit then the presence of another planet causes a secular advance of the transiting planet's pericenter over and above the effect of general relativity. Although this secular effect is impractical to detect over a small number of orbits, it causes long-term differences in when future transits occur, much like the long-term decay observed in pulsars. Measuring this transit-timing delay would thus allow the detection of either one or more additional planets in the system or the first measurements of non-zero oblateness ($J_2$) of the central stars.

We calculate the current spatial distribution of projectile delivery to the Earth and Moon using numerical orbital dynamics simulations of candidate impactors drawn from a debiased Near-Earth-Object (NEO) model. Surprisingly, we find that the average lunar impact velocity is 20 km/s, which has ramifications for converting observed crater densities to impactor size distributions. We determine that current crater production on the leading hemisphere of the Moon is 1.29 +/- 0.01 that of the trailing when considering the ratio of craters within 30 degrees of the apex to those within 30 degrees of the antapex and that there is virtually no nearside-farside asymmetry. As expected, the degree of leading-trailing asymmetry increases when the Moon's orbital distance is decreased. We examine the latitude distribution of impactor sites and find that for both the Earth and Moon there is a small deficiency of time-averaged impact rates at the poles. The ratio between deliveries within 30 degrees of the pole to that of a 30 degree band centered on the equator is nearly unity for Earth (<1%)(0.992 +/- 0.001) but detectably non-uniform for the Moon (~10%)(0.912 +/- 0.004). The terrestrial arrival results are examined to determine the degree of AM/PM asymmetry to compare with meteorite fall times (of which there seems to be a PM excess). Our results show that the impact flux of objects derived from the NEOs in the AM hours is ~2 times that of the PM hemisphere, further supporting the assertion that meteorite-dropping objects are recent ejections from the main asteroid belt rather than young fragments of NEOs.

We report the discovery of the first trans-neptunian object, designated 2004 XR190, with a nearly-cirular orbit beyond the 2:1 mean-motion resonance. Fitting an orbit to 23 astrometric observations spread out over 12 months yields an orbit of a=57.2\pm0.4, e=0.08\pm0.04, and i=46.6 deg. All viable orbits have perihelia distances q>49 AU. The very high orbital inclination of this extended scattered disk object might be explained by several models, but its existence again points to a large as-yet undiscovered population of transneptunian objects with large orbital perihelia and inclination.

We present the data acquisition strategy and characterization procedures for the Canada-France Ecliptic Plane Survey (CFEPS), a sub-component of the Canada-France-Hawaii Telescope Legacy Survey. The survey began in early 2003 and as of summer 2005 has covered 430 square degrees of sky within a few degrees of the ecliptic. Moving objects beyond the orbit of Uranus are detected to a magnitude limit of $m_R$=23 -- 24 (depending on the image quality). To track as large a sample as possible and avoid introducing followup bias, we have developed a multi-epoch observing strategy that is spread over several years. We present the evolution of the uncertainties in ephemeris position and orbital elements as the objects progress through the epochs. We then present a small 10-object sample that was tracked in this manner as part of a preliminary survey starting a year before the main CFEPS project. We describe the CFEPS survey simulator, to be released in 2006, which allows theoretical models of the Kuiper Belt to be compared with the survey discoveries since CFEPS has a well-documented pointing history with characterized detection efficiencies as a function of magnitude and rate of motion on the sky. Using the pre-survey objects we illustrate the usage of the simulator in modeling the classical Kuiper Belt.

We present simulations of Triton's post-capture orbit that confirm the importance of Kozai-type oscillations in its orbital elements. In the context of the tidal orbital evolution model, these variations require average pericenter distances much higher than previously published, and the timescale for the tidal orbital evolution of Triton becomes longer than the age of the Solar System. Recently-discovered irregular satellites present a new constraint on Triton's orbital history. Our numerical integrations of test particles indicate a timescale for Triton's orbital evolution to be less than $10^5$ yrs for a reasonable number of distant satellites to survive Triton's passage. This timescale is inconsistent with the exclusively tidal evolution (time scale of $>10^8$ yrs), but consistent with the interestion with the debris from satellite-satellite collisions. Any major regular satellites will quickly collide among themselves after being perturbed by Triton, and the resulting debris disk would eventually be swept up by Triton; given that the total mass of the Uranian satellite system is 40% of that of Triton, large scale evolution is possible. This scenario could have followed either collisional or the recently-discussed three-body-interaction-based capture.

We use a serendipitously discovered overdensity of extremely red objects (EROs) to study the morphologies and cumulative surface number density of EROs in a dense environment. Our extremely deep imaging allows us to select very faint EROs, reaching Ks=21, or ~2 magnitudes fainter than the L* of passively evolving ellipticals at z=1.5. We find that the shape of the ERO cumulative surface number density in our overdense field mimics that of the field ERO population over all magnitudes down to Ks=21 but with a factor of 3-4 higher normalization. The excellent seeing in our images (0.4" in Ks and 0.6" in R) allows for morphological classification of the brighter (Ks<19) EROs and we find a mix of morphologies including interacting systems and disks; the fraction of pure bulges (at most 38%), galaxies with disks (at least 46%), and interacting systems (at least 21%) is consistent with morphological fractions in field ERO studies. The similarity in the shape of the cumulative surface density and morphological mix between our overdense field and the field ERO population suggests that ERO galaxies in overdense regions at z~1-2 may not have had an appreciably different history than those in the field.

We present BVRI colors of 13 Jovian and 8 Saturnian irregular satellites obtained with the 2.56m Nordic Optical Telescope on La Palma, the 6.5m Magellan Baade Telescope on La Campanas, and the 6m MMT on Mt. Hopkins. The observations were performed between December 2001 to March 2002. Nearly all of the known irregular satellites can be divided into two distinct classes based on their colors. One, the grey color class, has the similar colors to the C-type asteroid, and the other, the light red color class, has colors similar to P/D-type asteroids. We also find at least one object, the Jovian irregular J XXIII Kalyke, that has colors similar to the red colored Centaurs/TNOs, although its classification is unsecure. We also find that there is a correlation between the physical properties and dynamical properties of the irregular satellites. Most of the dynamical clusters have homogeneous colors, which points to single homogeneous progenitors being cratered or fragmented as the source of each individual cluster. The heterogeneous colored clusters are most easily explained by assuming that there are several dynamical clusters in the area, rather than just one.

By telescopic tracking, we have established that the orbit of the trans-neptunian object (2000 CR$_{105}$) has a perihelion of $\simeq$44 AU, and is thus outside the domain controlled by strong gravitational close encounters with Neptune. Because this object is on a very large, eccentric orbit (with semimajor axis $a\simeq$216 AU and eccentricity $e\simeq$0.8) this object must have been placed on this orbit by a gravitational perturbation which is {\it not} direct gravitational scattering off of any of the giant planets (on their current orbits). The existence of this object may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced it, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of trans-neptunian objects in an `extended scattered disk' with perihelia above the previously-discussed 38 AU boundary.

We have conducted pencil-beam searches for outer solar system objects to a limiting magnitude of R ~ 26. Five new trans-neptunian objects were detected in these searches. Our combined data set provides an estimate of ~90 trans-neptunian objects per square degree brighter than ~ 25.9. This estimate is a factor of 3 above the expected number of objects based on an extrapolation of previous surveys with brighter limits, and appears consistent with the hypothesis of a single power-law luminosity function for the entire trans-neptunian region. Maximum likelihood fits to all self-consistent published surveys with published efficiency functions predicts a cumulative sky density Sigma(<R) obeying log10(Sigma) = 0.76(R-23.4) objects per square degree brighter than a given magnitude R.

Motivated by a desire to understand the size distributon of the Kuiper Belt, an observing program was conducted at the Mount Palomar 5-m telescope from 1994-1996. The observations consisted of follow-up observations of known objects (in order to improve their very indeterminate orbits), and deep exposures on a single field to search for small objects below the limiting magnitudes of other surveys. Eighteen object recoveries were successfully obtained over the course of the follow-up program. Data reduction of the deep fields consisted of a software recombination of many fields shifted at different angular rates in order to detect objects at differing heliocentric distances. We set an upper limit of <1 object per 0.05 square degrees in the ecliptic brighter than magnitude R ~ 25. The lack of detected objects in this work serves to help constrain the number density of Kuiper Belt comets.