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 Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

We present a two-component Machine Learning (ML) based approach for classifying astronomical images by data-quality via an examination of sources detected in the images and image pixel values from representative sources within those images. The first component, which uses a clustering algorithm, creates a proper and small fraction of the image pixels to determine the quality of the observation. The representative images (and associated tables) are ~800 times smaller than the original images, significantly reducing the time required to train our algorithm. The useful information in the images is preserved, permitting them to be classified in different categories, but the required storage is reduced. The second component, which is a deep neural network model, classifies the representative images. Using ground-based telescope imaging data, we demonstrate that the method can be used to separate 'usable' images from those that present some problems for scientific projects -- such as images that were taken in sub-optimal conditions. This method uses two different data sets as input to a deep model and provides better performance than if we only used the images' pixel information. The method may be used in cases where large and complex data sets should be examined using deep models. Our automated classification approach achieves 97% agreement when compared to classification generated via manual image inspection. We compare our method with traditional results and show that the method improves the results by about 10%, and also presents more comprehensive outcomes.

The outer Solar System object (486958) Arrokoth (provisional designation 2014 MU$_{69}$) has been largely undisturbed since its formation. We study its surface composition using data collected by the New Horizons spacecraft. Methanol ice is present along with organic material, which may have formed through radiation of simple molecules. Water ice was not detected. This composition indicates hydrogenation of carbon monoxide-rich ice and/ or energetic processing of methane condensed on water ice grains in the cold, outer edge of the early Solar System. There are only small regional variations in color and spectra across the surface, suggesting Arrokoth formed from a homogeneous or well-mixed reservoir of solids. Microwave thermal emission from the winter night side is consistent with a mean brightness temperature of 29$\pm$5 K.

We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

The International Virtual Observatory Alliance (IVOA) held its bi-annual Interoperability Meetings in May 2019, and in October 2019 following the ADASS XXIX conference. We provide a brief report on the status of the IVOA and the activities of the Interoperability Meetings.

A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

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.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

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.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

The application of machine learning (ML) methods to the analysis of astrophysical datasets is on the rise, particularly as the computing power and complex algorithms become more powerful and accessible. As the field of ML enjoys a continuous stream of breakthroughs, its applications demonstrate the great potential of ML, ranging from achieving tens of millions of times increase in analysis speed (e.g., modeling of gravitational lenses or analysing spectroscopic surveys) to solutions of previously unsolved problems (e.g., foreground subtraction or efficient telescope operations). The number of astronomical publications that include ML has been steadily increasing since 2010. With the advent of extremely large datasets from a new generation of surveys in the 2020s, ML methods will become an indispensable tool in astrophysics. Canada is an unambiguous world leader in the development of the field of machine learning, attracting large investments and skilled researchers to its prestigious AI Research Institutions. This provides a unique opportunity for Canada to also be a world leader in the application of machine learning in the field of astrophysics, and foster the training of a new generation of highly skilled researchers.

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.

Looking at the orbits of small bodies with large semimajor axes, we are compelled to see patterns. Some of these patterns are noted as strong indicators of new or hidden processes in the outer Solar System, others are substantially generated by observational biases, and still others may be completely overlooked. We can gain insight into the current and past structure of the outer Solar System through a careful examination of these orbit patterns. In this chapter, we discuss the implications of the observed orbital distribution of distant trans-Neptunian objects (TNOs). We start with some cautions on how observational biases must affect the known set of TNO orbits. Some of these biases are intrinsic to the process of discovering TNOs, while others can be reduced or eliminated through careful observational survey design. We discuss some orbital element correlations that have received considerable attention in the recent literature. We examine the known TNOs in the context of the gravitational processes that the known Solar System induces in orbital distributions. We discuss proposed new elements of the outer Solar System, posited ancient processes, and the types of TNO orbital element distributions that they predict to exist. We conclude with speculation.

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.

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.

The Hubble Space Telescope (HST) measured the colors of eight Kuiper Belt Objects (KBOs) that will be observed by the New Horizons spacecraft including its 2019 close fly-by target the Cold Classical KBO (486958) 2014 MU69. We find that the photometric colors of all eight objects are red, typical of the Cold Classical dynamical population within which most reside. Because 2014 MU69 has a similar color to that of other KBOs in the Cold Classical region of the Kuiper Belt, it may be possible to use the upcoming high-resolution New Horizons observations of 2014 MU69 to draw conclusions about the greater Cold Classical population. Additionally, HST found none of these KBOs to be binary within separations of ~0.06 arcsec (~2000 km at 44 AU range) and {\Delta}m less than or equal to 0.5. This conclusion is consistent with the lower fraction of binaries found at relatively wide separations. A few objects appear to have significant photometric variability, but our observations are not of sufficient signal-to-noise or time duration for further interpretation.

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.

NASA's New Horizons spacecraft will conduct a close flyby of the cold classical Kuiper Belt Object (KBO) designated (486958) 2014 MU69 on January 1, 2019. At a heliocentric distance of 44 AU, "MU69" will be the most distant object ever visited by a spacecraft. To enable this flyby, we have developed an extremely high precision orbit fitting and uncertainty processing pipeline, making maximal use of the Hubble Space Telescope's Wide Field Camera 3 (WFC3) and pre-release versions of the ESA Gaia Data Release 2 (DR2) catalog. This pipeline also enabled successful predictions of a stellar occultation by MU69 in July 2017. We describe how we process the WFC3 images to match the Gaia DR2 catalog, extract positional uncertainties for this extremely faint target (typically 140 photons per WFC3 exposure), and translate those uncertainties into probability distribution functions for MU69 at any given time. We also describe how we use these uncertainties to guide New Horizons, plan stellar occultions of MU69, and derive MU69's orbital evolution and long-term stability.

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.

The recent discovery by Pan-STARRS1 of 1I/2017 U1 (`Oumuamua), on an unbound and hyperbolic orbit, offers a rare opportunity to explore the planetary formation processes of other stars, and the effect of the interstellar environment on a planetesimal surface. 1I/`Oumuamua's close encounter with the inner Solar System in 2017 October was a unique chance to make observations matching those used to characterize the small-body populations of our own Solar System. We present near-simultaneous g$^\prime$, r$^\prime$, and J photometry and colors of 1I/`Oumuamua from the 8.1-m Frederick C. Gillett Gemini North Telescope, and $gri$ photometry from the 4.2 m William Herschel Telescope. Our g$^\prime$r$^\prime$J observations are directly comparable to those from the high-precision Colours of the Outer Solar System Origins Survey (Col-OSSOS), which offer unique diagnostic information for distinguishing between outer Solar System surfaces. The J-band data also provide the highest signal-to-noise measurements made of 1I/`Oumuamua in the near-infrared. Substantial, correlated near-infrared and optical variability is present, with the same trend in both near-infrared and optical. Our observations are consistent with 1I/`Oumuamua rotating with a double-peaked period of $8.10 \pm 0.42$ hours and being a highly elongated body with an axial ratio of at least 5.3:1, implying that it has significant internal cohesion. The color of the first interstellar planetesimal is at the neutral end of the range of Solar System $g-r$ and $r-J$ solar-reflectance colors: it is like that of some dynamically excited objects in the Kuiper belt and the less-red Jupiter Trojans.

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.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

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.

This work explores the orbital distribution of minor bodies in the outer Solar System emplaced as a result of a Nice model migration from the simulations of Brasser & Morbidelli (2013). This planetary migration scatters a planetesimal disk from between 29-34 AU and emplaces a population of objects into the Kuiper belt region. From the 2:1 Neptune resonance and outward, the test particles analyzed populate the outer resonances with orbital distributions consistent with trans-Neptunian object (TNO) detections in semi-major axis, inclination, and eccentricity, while capture into the closest resonances is too efficient. The relative populations of the simulated scattering objects and resonant objects in the 3:1 and 4:1 resonances are also consistent with observed populations based on debiased TNO surveys, but the 5:1 resonance is severely underpopulated compared to population estimates from survey results. Scattering emplacement results in the expected orbital distribution for the majority of the TNO populations, however the origin of the large observed population in the 5:1 resonance remains unexplained.

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.

Photometry of moving sources typically suffers from reduced signal-to-noise (SNR) or flux measurements biased to incorrect low values through the use of circular apertures. To address this issue we present the software package, TRIPPy: TRailed Image Photometry in Python. TRIPPy introduces the pill aperture, which is the natural extension of the circular aperture appropriate for linearly trailed sources. The pill shape is a rectangle with two semicircular end-caps, and is described by three parameters, the trail length and angle, and the radius. The TRIPPy software package also includes a new technique to generate accurate model point-spread functions (PSF) and trailed point-spread functions (TSF) from stationary background sources in sidereally tracked images. The TSF is merely the convolution of the model PSF, which consists of a moffat profile, and super sampled lookup table. From the TSF, accurate pill aperture corrections can be estimated as a function of pill radius with a accuracy of 10 millimags for highly trailed sources. Analogous to the use of small circular apertures and associated aperture corrections, small radius pill apertures can be used to preserve signal-to-noise of low flux sources, with appropriate aperture correction applied to provide an accurate, unbiased flux measurement at all SNR.

Circumstellar disks have long been regarded as windows into planetary systems. The advent of high sensitivity, high resolution imaging in the submillimetre where both the solid and gas components of disks can be detected opens up new possibilities for understanding the dynamical histories of these systems and therefore, a better ability to place our own solar system, which hosts a highly evolved debris disk, in context. Comparisons of dust masses from protoplanetary and debris disks have revealed a stark downturn in mass in millimetre-sized grains around a stellar age of 10 Myr, ostensibly in the "transition disk" phase, suggesting a period of rapid accretion of such grains onto planetesimals. This rapid formation phase is in keeping with radionucleide studies of Kuiper Belt Objects in the solar system. Importantly, this suggests that any thermal gradients in the gas of disks of this era will be "frozen in" to the planetesimals as they rapidly accrete from the solids and ices in their vicinity. Measurements of radial gradients in thermal tracers such as DHO, DCN and other tracers can therefore provide insight into the nascent solar system's abudances. In studies of dynamical evolution of the solar system, it is tacitly assumed that such abundances can reveal the location of formation for bodies now found in the asteroid belt and Kuiper belt. Similarly, evidence of gas detected from collisional evolution in young debris disks could potentially reveal how rapidly objects have dynamically evolved in those systems, most of which will be significantly younger than the solar system.

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 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 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 Next Generation Virgo Cluster Survey is an optical imaging survey covering 104 deg^2 centered on the Virgo cluster. Currently, the complete survey area has been observed in the u*giz-bands and one third in the r-band. We present the photometric redshift estimation for the NGVS background sources. After a dedicated data reduction, we perform accurate photometry, with special attention to precise color measurements through point spread function-homogenization. We then estimate the photometric redshifts with the Le Phare and BPZ codes. We add a new prior which extends to iAB = 12.5 mag. When using the u*griz-bands, our photometric redshifts for 15.5 \le i \lesssim 23 mag or zphot \lesssim 1 galaxies have a bias |\Delta z| < 0.02, less than 5% outliers, and a scatter \sigma_{outl.rej.} and an individual error on zphot that increase with magnitude (from 0.02 to 0.05 and from 0.03 to 0.10, respectively). When using the u*giz-bands over the same magnitude and redshift range, the lack of the r-band increases the uncertainties in the 0.3 \lesssim zphot \lesssim 0.8 range (-0.05 < \Delta z < -0.02, \sigma_{outl.rej} ~ 0.06, 10-15% outliers, and zphot.err. ~ 0.15). We also present a joint analysis of the photometric redshift accuracy as a function of redshift and magnitude. We assess the quality of our photometric redshifts by comparison to spectroscopic samples and by verifying that the angular auto- and cross-correlation function w(\theta) of the entire NGVS photometric redshift sample across redshift bins is in agreement with the expectations.

Many nearby main-sequence stars have been searched for debris using the far-infrared Herschel satellite, within the DEBRIS, DUNES and Guaranteed-Time Key Projects. We discuss here 11 stars of spectral types A to M where the stellar inclination is known and can be compared to that of the spatially-resolved dust belts. The discs are found to be well aligned with the stellar equators, as in the case of the Sun's Kuiper belt, and unlike many close-in planets seen in transit surveys. The ensemble of stars here can be fitted with a star-disc tilt of ~<10 degrees. These results suggest that proposed mechanisms for tilting the star or disc in fact operate rarely. A few systems also host imaged planets, whose orbits at tens of AU are aligned with the debris discs, contrary to what might be expected in models where external perturbers induce tilts.

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.

We report the discovery of 2010 GB$_{174}$, a likely new member of the Inner Oort Cloud (IOC). 2010 GB$_{174}$ is one of 91 Trans Neptunian Objects (TNOs) and Centaurs discovered in a 76 deg$^2$ contiguous region imaged as part of the Next Generation Virgo Cluster Survey (NGVS) --- a moderate ecliptic latitude survey reaching a mean limiting magnitude of $g^\prime \simeq 25.5$ --- using MegaPrime on the 3.6m Canada France Hawaii Telescope. 2010 GB$_{174}$ is found to have an orbit with semi-major axis $a\simeq350.8$ AU, inclination $i \simeq 21.6^\circ$ and pericentre $q\sim48.5$ AU. This is the second largest perihelion distance among known solar system objects. Based on the sky coverage and depth of the NGVS, we estimate the number of IOC members with sizes larger than 300 km ($H_V \le 6.2$ mag) to be $\simeq 11\,000$. A comparison of the detection rate from the NGVS and the PDSSS (a characterized survey that `re-discovered' the IOC object Sedna) gives, for an assumed a power-law LF for IOC objects, a slope of $\alpha \simeq 0.7 \pm 0.2$, with only two detections in this region this slope estimate is highly uncertain.

Here we present high cadence photometry taken by the Acquisition Camera on Gemini South, of a close passage by the $\sim540$ km radius Kuiper Belt Object, (50000) Quaoar, of a r'=20.2 background star. Observations before and after the event show that the apparent impact parameter of the event was $0.019\pm0.004$", corresponding to a close approach of $580\pm120$ km to the centre of Quaoar. No signatures of occultation by either Quaoar's limb or its potential atmosphere are detectable in the relative photometry of Quaoar and the target star, which were unresolved during closest approach. From this photometry we are able to put constraints on any potential atmosphere Quaoar might have. Using a Markov chain Monte Carlo and likelihood approach, we place pressure upper limits on sublimation supported, isothermal atmospheres of pure N$_2$, CO, and CH$_4$. For N$_2$ and CO, the upper limit surface pressures are 1 and 0.7 $\mu{bar}$ respectively. The surface temperature required for such low sublimation pressures is $\sim33$ K, much lower than Quaoar's mean temperature of $\sim44$ K measured by others. We conclude that Quaoar cannot have an isothermal N$_2$ or CO atmosphere. We cannot eliminate the possibility of a CH$_4$ atmosphere, but place upper surface pressure and mean temperature limits of $\sim138$ nbar and $\sim44$ K respectively.

Color measurements and albedo distributions introduce non-intuitive observational biases in size-color relationships among Kuiper Belt Objects (KBOs) that cannot be disentangled without a well characterized sample population with systematic photometry. Peixinho et al. report that the form of the KBO color distribution varies with absolute magnitude, H. However, Tegler et al. find that KBO color distributions are a property of object classification. We construct synthetic models of observed KBO colors based on two B-R color distribution scenarios: color distribution dependent on H magnitude (H-Model) and color distribution based on object classification (Class-Model). These synthetic B-R color distributions were modified to account for observational flux biases. We compare our synthetic B-R distributions to the observed 'Hot' and 'Cold' detected objects from the Canada-France Ecliptic Plane Survey and the Meudon Multicolor Survey. For both surveys, the Hot population color distribution rejects the H-Model, but is well described by the Class-Model. The Cold objects reject the H-Model, but the Class-Model (while not statistically rejected) also does not provide a compelling match for data. Although we formally reject models where the structure of the color distribution is a strong function of H magnitude, we also do not find that a simple dependence of color distribution on orbit classification is sufficient to describe the color distribution of classical KBOs.

Here we present observations of 7 large Kuiper Belt Objects. From these observations, we extract a point source catalog with $\sim0.01"$ precision, and astrometry of our target Kuiper Belt Objects with $0.04-0.08"$ precision within that catalog. We have developed a new technique to predict the future occurrence of stellar occultations by Kuiper Belt Objects. The technique makes use of a maximum likelihood approach which determines the best-fit adjustment to cataloged orbital elements of an object. Using simulations of a theoretical object, we discuss the merits and weaknesses of this technique compared to the commonly adopted ephemeris offset approach. We demonstrate that both methods suffer from separate weaknesses, and thus, together provide a fair assessment of the true uncertainty in a particular prediction. We present occultation predictions made by both methods for the 7 tracked objects, with dates as late as 2015. Finally, we discuss observations of three separate close passages of Quaoar to field stars, which reveal the accuracy of the element adjustment approach, and which also demonstrate the necessity of considering the uncertainty in stellar position when assessing potential occultations.

We present the discovery of a long-term stable L5 (trailing) Neptune Trojan in data acquired to search for candidate Trans-Neptunian objects for the New Horizons spacecraft to fly by during an extended post-Pluto mission. This Neptune Trojan, 2011 HM102, has the highest inclination (29.4 degrees) of any known member of this population. It is intrinsically brighter than any single L5 Jupiter Trojan at H~8.18. We have determined its gri colors (a first for any L5 Neptune Trojan), which we find to be similar to the moderately red colors of the L4 Neptune Trojans, suggesting similar surface properties for members of both Trojan clouds. We also present colors derived from archival data for two L4 Neptune Trojans (2006 RJ103 and 2007 VL305), better refining the overall color distribution of the population. In this document we describe the discovery circumstances, our physical characterization of 2011 HM102, and this object's implications for the Neptune Trojan population overall. Finally, we discuss the prospects for detecting 2011 HM102 from the New Horizons spacecraft during their close approach in mid- to late-2013.

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.

We present observations of the debris disk around gamma Doradus, an F1V star, from the Herschel Key Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/Submillimetre). The disk is well-resolved at 70, 100 and 160 micron, resolved along its major axis at 250 micron, detected but not resolved at 350 micron, and confused with a background source at 500 micron. It is one of our best resolved targets and we find it to have a radially broad dust distribution. The modelling of the resolved images cannot distinguish between two configurations: an arrangement of a warm inner ring at several AU (best-fit 4 AU) and a cool outer belt extending from ~55 to 400 AU or an arrangement of two cool, narrow rings at ~70 AU and ~190 AU. This suggests that any configuration between these two is also possible. Both models have a total fractional luminosity of ~10^{-5} and are consistent with the disk being aligned with the stellar equator. The inner edge of either possible configuration suggests that the most likely region to find planets in this system would be within ~55 AU of the star. A transient event is not needed to explain the warm dust's fractional luminosity.

Debris disks have been found primarily around intermediate and solar mass stars (spectral types A-K) but rarely around low mass M-type stars. We have spatially resolved a debris disk around the remarkable M3-type star GJ581 hosting multiple planets using deep PACS images at 70, 100 and 160 microns as part of the DEBRIS Program on the Herschel Space Observatory. This is the second spatially resolved debris disk found around an M-type star, after the one surrounding the young star AU Mic (12 Myr). However, GJ 581 is much older (2-8 Gyr), and is X-ray quiet in the ROSAT data. We fit an axisymmetric model of the disk to the three PACS images and found that the best fit model is for a disk extending radially from 25+/-12 AU to more than 60 AU. Such a cold disk is reminiscent of the Kuiper Belt but it surrounds a low mass star (0.3 M_sol) and its fractional dust luminosity L_dust/L_* of \sim 10^-4 is much higher. The inclination limits of the disk found in our analysis make the masses of the planets small enough to ensure the long-term stability of the system according to some dynamical simulations. The disk is collisionally dominated down to submicron-sized grains and the dust cannot be expelled from the system by radiation or wind pressures because of the low luminosity and low X-ray luminosity of GJ581. We suggest that the correlation between low-mass planets and debris disks recently found for G-type stars also applies to M-type stars.Finally, the known planets, of low masses and orbiting within 0.3 AU from the star, cannot dynamically perturb the disk over the age of the star, suggesting that an additional planet exists at larger distance that is stirring the disk to replenish the dust.

The majority of debris discs discovered so far have only been detected through infrared excess emission above stellar photospheres. While disc properties can be inferred from unresolved photometry alone under various assumptions for the physical properties of dust grains, there is a degeneracy between disc radius and dust temperature that depends on the grain size distribution and optical properties. By resolving the disc we can measure the actual location of the dust. The launch of Herschel, with an angular resolution superior to previous far-infrared telescopes, allows us to spatially resolve more discs and locate the dust directly. Here we present the nine resolved discs around A stars between 20 and 40 pc observed by the DEBRIS survey. We use these data to investigate the disc radii by fitting narrow ring models to images at 70, 100 and 160 {\mu}m and by fitting blackbodies to full spectral energy distributions. We do this with the aim of finding an improved way of estimating disc radii for unresolved systems. The ratio between the resolved and blackbody radii varies between 1 and 2.5. This ratio is inversely correlated with luminosity and any remaining discrepancies are most likely explained by differences to the minimum size of grain in the size distribution or differences in composition. We find that three of the systems are well fit by a narrow ring, two systems are borderline cases and the other four likely require wider or multiple rings to fully explain the observations, reflecting the diversity of planetary systems.

This paper describes Herschel observations of the nearby (8.5pc) G5V multi-exoplanet host star 61 Vir at 70-500micron carried out by the DEBRIS survey. These reveal emission that is extended out to >15arcsec with a morphology that can be fitted by a nearly edge-on (77deg inclination) radially broad (from 30AU to >100AU) debris disk of fractional luminosity 2.7x10^-5, with two unrelated sources nearby that are more prominent at longer wavelengths. Chance alignment with a background object seen at 1.4GHz provides potential for confusion, but the star's 1.4"/yr proper motion allows Spitzer 70micron images to confirm that what we are interpreting as disk emission really is circumstellar. Although the exact shape of the disk's inner edge is not well constrained, the region inside 30AU must be significantly depleted in planetesimals. This is readily explained if there are additional planets in the 0.5-30AU region, but is also consistent with collisional erosion. We also find tentative evidence that the presence of detectable debris around nearby stars correlates with the presence of the lowest mass planets that are detectable in current radial velocity surveys. Out of an unbiased sample of the nearest 60 G stars, 11 are known to have planets, of which 6 (including 61 Vir) have planets that are all less massive than Saturn, and 4 of these have evidence for debris. The debris toward one of these planet-hosts (HD20794) is reported here for the first time. This fraction (4/6) is higher than that expected for nearby field stars (15%), and implies that systems that form low-mass planets are also able to retain bright debris disks. We suggest that this correlation could arise because such planetary systems are dynamically stable and include regions that are populated with planetesimals in the formation process where the planetesimals can remain unperturbed over Gyr timescales.

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.

It is very difficult to find archival images of solar system objects. While regular archive searches can find images at a fixed location, they cannot find images of moving targets. Archival images have become increasingly useful to galactic and stellar astronomers the last few years but, until now, solar system researchers have been at a disadvantage in this respect. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for images containing that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, several ESO instruments and the SDSS for a total of 6.5 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos

While regular archive searches can find images at a fixed location, they cannot find images of moving targets such as asteroids. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for observations of that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, and several ESO instruments for a total of 1.6 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos

The widely-separated, near-equal mass binaries hosted by the cold Classical Kuiper Belt are delicately bound and subject to disruption by many perturbing processes. We use analytical arguments and numerical simulations to determine their collisional lifetimes given various impactor size distributions, and include the effects of mass-loss and multiple impacts over the lifetime of each system. These collisional lifetimes constrain the population of small (R > ~1 km) objects currently residing in the Kuiper Belt, and confirm that the size distribution slope at small size cannot be excessively steep - likely q < ~3.5. We track mutual semi-major axis, inclination, and eccentricity evolution through our simulations, and show that it is unlikely that the wide binary population represents an evolved tail of the primordially-tight binary population. We find that if the wide binaries are a collisionally-eroded population, their primordial mutual orbit planes must have preferred to lie in the plane of the solar system. Finally, we find that current limits on the size distribution at small radii remain high enough that the prospect of detecting dust-producing collisions in real-time in the Kuiper Belt with future optical surveys is feasible.

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.

The D/H enrichment observed in Saturn's satellite Enceladus is remarkably similar to the values observed in the nearly-isotropic comets. Given the predicted strong variation of D/H with heliocentric distance in the solar nebula, this observation links the primordial source region of the nearly-isotropic comets with the formation location of Enceladus. That is, comets from the nearly-isotropic class were most likely fed into their current reservoir, the Oort cloud, from a source region near the formation location of Enceladus. Dynamical simulations of the formation of the Oort cloud indicate that Uranus and Neptune are, primarily, responsible for the delivery of material into the Oort cloud. In addition, Enceladus formed from material that condensed from the solar nebula near the location at which Saturn captured its gas envelope, most likely at or near Saturn's current location in the solar system. The coupling of these lines of evidence appears to require that Uranus and Neptune were, during the epoch of the formation of the Oort cloud, much closer to the current location of Saturn than they are currently. Such a configuration is consistent with the Nice model of the evolution of the outer solar system. Further measurements of the D/H enrichment in comets, particularly in ecliptic comets, will provide an excellent discriminator among various models of the formation of the outer solar system.

The existence of extremely wide binaries in the low-inclination component of the Kuiper Belt provides a unique handle on the dynamical history of this population. Some popular frameworks of the formation of the Kuiper Belt suggest that planetesimals were moved there from lower semi-major axis orbits by scattering encounters with Neptune. We test the effects such events would have on binary systems, and find that wide binaries are efficiently destroyed by the kinds of scattering events required to create the Kuiper Belt with this mechanism. This indicates that a binary-bearing component of the cold Kuiper Belt was emplaced through a gentler mechanism or was formed in situ.

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 present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

Two populations of minor bodies in the outer Solar System remain particularly elusive: Scattered Disk objects and Sedna-like objects. These populations are important dynamical tracers, and understanding the details of their spatial- and size-distributions will enhance our understanding of the formation and on-going evolution of the Solar System. By using newly-derived limits on the maximum heliocentric distances that recent pencil-beam surveys for Trans-Neptunian Objects were sensitive to, we determine new upper limits on the total numbers of distant SDOs and Sedna-like objects. While generally consistent with populations estimated from wide-area surveys, we show that for magnitude-distribution slopes of {\alpha} > 0.7-1.0, these pencil-beam surveys provide stronger upper limits than current estimates in literature.

Digital co-addition of astronomical images is a common technique for increasing signal-to-noise and image depth. A modification of this simple technique has been applied to the detection of minor bodies in the Solar System: first stationary objects are removed through the subtraction of a high-SN template image, then the sky motion of the Solar System bodies of interest is predicted and compensated for by shifting pixels in software prior to the co-addition step. This "shift-and-stack" approach has been applied with great success in directed surveys for minor Solar System bodies. In these surveys, the shifts have been parameterized in a variety of ways. However, these parameterizations have not been optimized and in most cases cannot be effectively applied to data sets with long observation arcs due to objects' real trajectories diverging from linear tracks on the sky. This paper presents two novel probabilistic approaches for determining a near-optimum set of shift-vectors to apply to any image set given a desired region of orbital space to search. The first method is designed for short observational arcs, and the second for observational arcs long enough to require non-linear shift-vectors. Using these techniques and other optimizations, we derive optimized grids for previous surveys that have used "shift-and-stack" approaches to illustrate the improvements that can be made with our method, and at the same time derive new limits on the range of orbital parameters these surveys searched. We conclude with a simulation of a future applications for this approach with LSST, and show that combining multiple nights of data from such next-generation facilities is within the realm of computational feasibility.

GRB 051103 is considered to be a candidate soft gamma repeater (SGR) extragalactic giant magnetar flare by virtue of its proximity on the sky to M81/M82, as well as its time history, localization, and energy spectrum. We have derived a refined interplanetary network localization for this burst which reduces the size of the error box by over a factor of two. We examine its time history for evidence of a periodic component, which would be one signature of an SGR giant flare, and conclude that this component is neither detected nor detectable under reasonable assumptions. We analyze the time-resolved energy spectra of this event with improved time- and energy resolution, and conclude that although the spectrum is very hard, its temporal evolution at late times cannot be determined, which further complicates the giant flare association. We also present new optical observations reaching limiting magnitudes of R > 24.5, about 4 magnitudes deeper than previously reported. In tandem with serendipitous observations of M81 taken immediately before and one month after the burst, these place strong constraints on any rapidly variable sources in the region of the refined error ellipse proximate to M81. We do not find any convincing afterglow candidates from either background galaxies or sources in M81, although within the refined error region we do locate two UV bright star forming regions which may host SGRs. A supernova remnant (SNR) within the error ellipse could provide further support for an SGR giant flare association, but we were unable to identify any SNR within the error ellipse. These data still do not allow strong constraints on the nature of the GRB 051103 progenitor, and suggest that candidate extragalactic SGR giant flares will be difficult, although not impossible, to confirm.

A novel method of generating artificial scintillation noise is developed and used to evaluate occultation rates and false positive rates for surveys probing the Kuiper Belt with the method of serendipitous stellar occultations. A thorough examination of survey design shows that: (1) diffraction-dominated occultations are critically (Nyquist) sampled at a rate of 2 Fsu^{-1}, corresponding to 40 s^{-1} for objects at 40 AU, (2) occultation detection rates are maximized when targets are observed at solar opposition, (3) Main Belt Asteroids will produce occultations lightcurves identical to those of Kuiper Belt Objects if target stars are observed at solar elongations of: 116 deg < epsilon < 125 deg, or 131 deg < epsilon < 141 deg, and (4) genuine KBO occultations are likely to be so rare that a detection threshold of >7-8 sigma should be adopted to ensure that viable candidate events can be disentangled from false positives.

We analyze photometric data in V and I for the globular cluster (GC) systems in five of the giant ellipticals in the Coma Cluster: NGC 4874, 4881, 4889, 4926, and IC 4051. We find that the GC luminosity functions are quite similar to one another, with a turnover derived from a composite sample of more than 9,000 GCs at V = 27.71 +- 0.07 (M_V = -7.3). Both a simple Gaussian curve and an evolved Schechter function fit the bright half of the GCLF equally well, though the Coma GCLF is broader and has a higher ``cutoff mass'' (M_c ~ 3 x 10^6 M_Sun) than in any of the Virgo giants. These five Coma members exhibit a huge range in GC specific frequency, from a low of S_N = 0.7 for NGC 4881 up to 12 for IC 4051 and NGC 4874. No single formation scenario appears able to account for these differences in otherwise-similar galaxies. The supergiant NGC 4874 has the richest globular cluster system known, probably holding more than 30,000 clusters; its true extent is not yet determined and may extend well out into the Coma potential well. For the three biggest GC systems (NGC 4874, 4889, IC 4051), all three populations are dominated by red, metal-rich clusters. Their metallicity distributions also may all have the normal bimodal form, with the two sequences at the expected mean colors <V-I>(blue) = 0.98 and <V-I>(red) = 1.15. However, the color distributions and relative numbers of metal-rich clusters show intriguing counterexamples to a trend established by Peng et al. 2008 (ApJ 681, 197) for the Virgo galaxies. At the very highest-density and most massive regimes represented by the Coma supergiants, formation of metal-rich clusters seems to have been especially favored.

We have performed a survey of the Kuiper belt covering ~ 1/3 a square degree of the sky using Suprime-cam on the Subaru telescope, to a limiting magnitude of m(R)~ 26.8 (50% threshold) and have found 36 new KBOs. We have confirmed that the luminosity function of the Kuiper belt must break as previously observed (Bernstein et al. 2004; Fuentes & Holman 2008). From the luminosity function, we have inferred the underlying size distribution and find that it is consistent with a large object power-law slope q1~4.8 that breaks to a slope q2~1.9 at object diameter Db~60 km assuming 6% albedos. We have found no conclusive evidence that the size distribution of KBOs with inclinations i<5 is different than that of those with i>5. We discuss implications of this measurement for early accretion in the outer solar system and Neptune migration scenarios.

We have derived a model of the Kuiper belt luminosity function exhibited by a broken power-law size distribution. This model allows direct comparison of the observed luminosity function to the underlying size distribution. We discuss the importance of the radial distribution model in determining the break diameter. We determine a best-fit break-diameter of the Kuiper belt size-distribution of 30<Db<90 km via a maximum-likelihood fit of our model to the observed luminosity function. We also confirm that the observed luminosity function for m(R) ~ 21-28 is consistent with a broken power-law size distribution, and exhibits a break at m(R)=26.0+0.7-1.8.

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.

The results of a search for sub-km Kuiper Belt Objects (KBOs) with the method of serendipitous stellar occultations are reported. Photometric time series were obtained on the 1.8m telescope at the Dominion Astrophysical Observatory (DAO) in Victoria, BC, and were analyzed for the presence of occultation events. Observations were performed at 40 Hz and included a total of 5.0 star-hours for target stars in the ecliptic open cluster M35 (beta=0.9deg), and 2.1 star-hours for control stars in the off-ecliptic open cluster M34 (beta=25.7deg). To evaluate the recovery fraction of the analysis method, and thereby determine the limiting detectable size, artificial occultation events were added to simulated time series (1/f scintillation-like power-spectra), and to the real data. No viable candidate occultation events were detected. This limits the cumulative surface density of KBOs to 3.5e10 deg^{-2} (95% confidence) for KBOs brighter than m_R=35.3 (larger than ~860m in diameter, assuming a geometric albedo of 0.04 and a distance of 40 AU). An evaluation of TNO occultations reported in the literature suggests that they are unlikely to be genuine, and an overall 95%-confidence upper limit on the surface density of 2.8e9 deg^{-2} is obtained for KBOs brighter than m_R=35 (larger than ~1 km in diameter, assuming a geometric albedo of 0.04 and a distance of 40 AU) when all existing surveys are combined.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

We present optical measurements of the faint end of the luminosity function in the core of the Coma cluster. Dwarf galaxies are detected down to a limiting magnitude of R approx. 25.75 in images taken with the Hubble Space Telescope. This represents the faintest determination of the Coma luminosity function to date. With the assumption that errors due to cosmic variance are small, evidence is found for a steep faint end slope with alpha < -2. Such a value is expected in theories in which reionization and other feedback processes are dependent on density.

We attempt to detect optical GRB afterglows on images taken by the Canada France Hawaii Telescope for the Very Wide survey, component of the Legacy Survey. To do so, a Real Time Analysis System called "Optically Selected GRB Afterglows" has been installed on a dedicated computer in Hawaii. This pipeline automatically and quickly analyzes Megacam images and extracts from them a list of variable objects which is displayed on a web page for validation by a member of the collaboration. The Very Wide survey covers 1200 square degrees down to i'=23.5. This paper briefly explain the RTAS process.

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 new (B,I) photometry for the globular cluster systems in eight Brightest Cluster Galaxies (BCGs), obtained with the ACS/WFC camera on the Hubble Space Telescope. In the very rich cluster systems that reside within these giant galaxies, we find that all have strongly bimodal color distributions All the BCGs show population gradients, with much higher relative numbers of red clusters within 5 kpc of their centers, consistent with their having formed at later times than the blue, metal-poor population. A striking new feature of the color distributions emerging from our data is that for the brightest clusters (M_I < -10.5) the color distribution becomes broad and less obviously bimodal. we suggest that it may be a characteristic of many BCGs. Furthermore, the blue (metal-poor) clusters become progressively redder with increasing luminosity, following a mass/metallicity scaling relation Z ~ M^0.55. We argue that these GCS characteristics are consistent with a hierarchical-merging formation picture in which the metal-poor clusters formed in protogalactic clouds or dense starburst complexes with gas masses in the range 10^7 - 10^10 M_Sun, but where the more massive clusters on average formed in bigger clouds with deeper potential wells where more pre-enrichment could occur.

We determine the period, p = 11.52 \pm 0.14 h, and a light curve peak-to-peak amplitude, a = 0.029 \pm 0.003 magnitudes, of the Neptunian irregular satellite Nereid. If the light curve variation is due to albedo variations across the surface, rather than solely to the shape of Nereid variations, the rotation period would be a factor of two shorter. In either case, such a rotation period and light curve amplitude, together with Nereid's orbital period, p=360.14 days, imply that Nereid is almost certainly in a regular rotation state, rather than the chaotic rotation state suggested by Schaefer and Schaefer (1988,2000) and Dobrovolskis (1995). Assuming that Nereid is perfectly spherical, the albedo variation is 3% across the observed surface. Assuming a uniform geometric albedo, the observed cross sectional area varies by 3%. We caution that the lightcurve found in this paper only sets limits on the combination of albedo and physical irregularity and that we cannot determine the orientation of Nereid's spin axis from our data.

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.

Deep HST/WFPC2 (V,I) photometry is used to investigate the globular cluster system (GCS) in NGC 4874, the central cD galaxy of the Coma cluster. The luminosity function of the clusters displays its normal Gaussian-like shape and turnover level. Other features of the system are surprising: the GCS is (a) spatially extended, with core radius r_c = 22 kpc, (b) entirely metal-poor (a narrow, unimodal metallicity distribution with mean [Fe/H] = -1.5), and (c) modestly populated, with specific frequency S_N = 3.7 +- 0.5. We suggest on the basis of some simple models that as much as half of this galaxy might have accreted from low-mass satellites, but no single one of the three classic modes of galaxy formation (accretion, disk mergers, in situ formation) can supply a fully satisfactory formation picture. Even when they are used in combination, strong challenges to these models remain. The principal anomaly in this GCS is essentially the complete lack of metal-rich clusters. If these were present in normal (M87-like) numbers in addition to the metal-poor ones that are already there, then the GCS in total would more closely resemble what we see in many other giant E galaxies.

We have used deep HST/WFPC2 images in V (F606W) and I (F814W) to measure the luminosity distribution of the globular clusters in NGC 4874, the central cD galaxy of the Coma cluster. We find the ``turnover'' point of the globular cluster luminosity function (GCLF) to lie at V = 27.88 +/- 0.12, while the overall GCLF shape matches the standard Gaussian-like form with dispersion SIGMA_V = 1.49 +/- 0.12. We use the GCLF as a standard candle by matching the turnover points in NGC 4874 and another Coma elliptical, IC 4051, with those of the giant ellipticals in the Virgo cluster (M87 and five others). The result is (m-M)_(Coma - Virgo) = 4.06 +/- 0.11 magnitudes, which converts to a Coma distance d = 102 Mpc if the Virgo distance modulus is (m-M)_0 = 30.99 +/- 0.04. The Hubble constant which emerges from our GCLF measurement is then H0 = (69 +/- 9) km/s/Mpc. We confirm this H0 value with a novel presentation of the ``Hubble diagram'' for GCLFs in giant E galaxies. Measurements of additional GCLFs in the Coma ellipticals, as well as calibrating galaxies in Virgo and Fornax, have excellent potential to refine this result in the near future.

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.

The universality of the globular cluster luminosity function (GCLF) contrasts the variation seen in the specific frequency ($S_N$). The variation in $S_N$ has been shown to follow a linear relation with $L_X$ for brightest cluster galaxies (Blakeslee, J. 1997, Harris, W.E. et al. 1998). Further, the variation of $S_N$ with galactic radius within individual giant ellipticals is seen to be a constant fraction of the gas density (McLaughlin 1998). There are now a number of galaxies for which direct mass estimates based on the radial velocities of the globular clusters are available. By comparing the mass of galaxies determined in this way with the number of clusters within these galaxies we show that the fraction of mass which is converted into globular clusters is constant independent of galaxy type or environment. This implies that the process of globular cluster formation is not influenced by the host galaxy and supports the notion of the universal GCLF.

The two brightest elliptical galaxies in the Eridanus A group, NGC 1400 and NGC 1407, have been observed in both the Washington T_1 and Kron-Cousins I filters to obtain photometry of their globular cluster systems (GCSs). This group of galaxies is of particular interest due to its exceptionally high M/L value, previously estimated at ~3000h, making this cluster highly dark-matter-dominated. NGC 1400's radial velocity (549 km/s) is extremely low compared to that of the central galaxy of Eridanus A (NGC 1407 with $v_\odot$ = 1766 km/s) and the other members of the system, suggesting that it is a foreground galaxy projected by chance onto the cluster. Using the shapes of the globular cluster luminosity functions, however, we derive distances of 17.6 +/- 3.1 Mpc to NGC 1407 and 25.4 +/- 7.0 Mpc to NGC 1400. These results support earlier conclusions that NGC 1400 is at the distance of Eridanus A and therefore has a large peculiar velocity. Specific frequencies are also derived for these galaxies, yielding values of S_N = 4.0 +/- 1.3 for NGC 1407 and S_N = 5.2 +/- 2.0 for NGC 1400. In this and other respects, these two galaxies have GCSs which are consistent with those observed in other galaxies.

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.

We have obtained spectra for globular cluster candidates in M104 with LDSS-2 on the William Herschel Telescope, confirming 34 objects as M104 globular clusters. We find a cluster velocity dispersion of $\sim$ 260 km/sec, and the Projected Mass Estimator gives a mass of 5.0 (3.5,6.7) $\times 10^{11}$ M$_\odot$ for M104 within a projected radius of $\sim 330^{''}$ (14 kpc for D=8.55 Mpc). Our best estimate for the mass-to-light ratio is M/L$_{V_T}$= 16$^{+5.5}_{-5.0}$ within the same radius. Considering all of the possible sources of uncertainty, we find a lower limit of M/L_V= 5.3, which is larger than the M/L$_V$ found from rotation curve analyses inside 180$^{''}$. We thus conclude that the mass-to-light ratio increases with radius, or in other words that M104 possesses a dark matter halo. There is a marginal detection of rotation in the M104 cluster system at the 92.5% confidence level; larger samples will be needed to investigate this possibility. Interestingly, the M104 globular cluster and planetary nebulae (PNe) kinematics are roughly consistent inside $\sim 100^{''}$. Finally, we find a mean cluster metallicity of [Fe/H] = $-$0.70 $\pm$ 0.3, which is more typical of clusters in gE/cD galaxies than it is of clusters in other spirals.

We show that the globular cluster luminosity function (GCLF) of the inner halo of the Milky Way is statistically different from the GCLF of the outer halo. We also find a similar difference between the inner and outer halo population of M31. We assert that this difference is evidence for some form of dynamical evolution of the cluster population and/or a dependence of GCLF shape on the environment in which the cluster population formed. We also find that the turnover luminosity of the GCLF is unaffected by these differences and further assert that this stability of the turnover luminosity affirms its usefulness as an indicator of cosmic distance.

We present the results of a $BVI$ photometric study of the globular cluster NGC 6517 using CCD images obtained at the prime focus of the CFHT. From the resultant color-magnitude diagram, we infer that the cluster is differentially reddened and develop a method for determining the amplitude of this effect, which amounts to \( \sim 0.4 \) mag in \( (V-I) \) across our $\sim 4$ arc-min field. From the corrected $V_0,$ $(V-I)_0$ color-magnitude diagram, we derive the metallicity and distance modulus of the cluster. Our values for the mean foreground reddening ($E(B-V) = 1.10 \pm 0.10$) and metallicity (\feh\ = -1.58 $\pm$ 0.05) agree to within uncertainties with previous determinations. Our apparent distance modulus of 18.3 $\pm$ 0.2, however, is larger than the sole previous determination. When the distance modulus is corrected for absorption, assuming that $R_V = 3.8 $ for this line of sight, we find that NGC 6517 is within 3 kpc of the Galactic center, making it a bulge member. This makes NGC 6517 only the fourth metal poor cluster in the bulge for which an accurate CMD is available.