Search

Sparse and serendipitous asteroid photometry obtained by wide field surveys such as the Asteroid Terrestrial-impact Last Alert System (\ATLAS) is a valuable resource for studying the properties of large numbers of small Solar System bodies. We have gathered a large database of \ATLAS photometry in wideband optical cyan and orange filters, consisting of 9.6\e{7} observations of 4.5\e{5} main belt asteroids and Jupiter Trojans. We conduct a phase curve analysis of these asteroids considering each apparition separately, allowing us to accurately reject outlying observations and to remove apparitions and asteroids not suitable for phase curve determination. We obtain a dataset of absolute magnitudes and phase parameters for over 100,000 selected asteroids observed by \ATLAS, $\sim66,000$ of which had sufficient measurements to derive colours in the \ATLAS filters. To demonstrate the power of our dataset we consider the properties of the Nysa-Polana complex, for which the \ATLAS colours and phase parameters trace the S-like and C-like compositions amongst family members. We also compare the properties of the leading and trailing groups of Jupiter Trojans, finding no significant differences in their phase parameters or colours as measured by \ATLAS, supporting the consensus that these groups were captured from a common source population during planetary migration. Furthermore, we identify $\sim9000$ asteroids that exhibit large shifts in derived absolute magnitude between apparitions, indicating that these objects have both elongated shapes and spin axes with obliquity $\sim 90$ degrees., Comment: 23 pages, 19 figures, published in MNRAS. Minor copyedit changes and updated links to data repositories

The interstellar object population of the Milky Way is a product of its stars. However, what is in fact a complex structure in the Solar neighbourhood has traditionally in ISO studies been described as smoothly distributed. Using a debiased stellar population derived from the Gaia DR3 stellar sample, we infer that the velocity distribution of ISOs is far more textured than a smooth Gaussian. The moving groups caused by Galactic resonances dominate the distribution. 1I/`Oumuamua and 2I/Borisov have entirely normal places within these distributions; 1I is within the non-coeval moving group that includes the Matariki (Pleiades) cluster, and 2I within the Coma Berenices moving group. We show that for the composition of planetesimals formed beyond the ice line, these velocity structures also have a chemodynamic component. This variation will be visible on the sky. We predict that this richly textured distribution will be differentiable from smooth Gaussians in samples that are within the expected discovery capacity of the Vera C. Rubin Observatory. Solar neighbourhood ISOs will be of all ages and come from a dynamic mix of many different populations of stars, reflecting their origins from all around the Galactic disk., Comment: Submitted to AJ

The impact of the DART spacecraft into Dimorphos, moon of the asteroid Didymos, changed Dimorphos' orbit substantially, largely from the ejection of material. We present results from twelve Earth-based facilities involved in a world-wide campaign to monitor the brightness and morphology of the ejecta in the first 35 days after impact. After an initial brightening of ~1.4 magnitudes, we find consistent dimming rates of 0.11-0.12 magnitudes/day in the first week, and 0.08-0.09 magnitudes/day over the entire study period. The system returned to its pre-impact brightness 24.3-25.3 days after impact through the primary ejecta tail remained. The dimming paused briefly eight days after impact, near in time to the appearance of the second tail. This was likely due to a secondary release of material after re-impact of a boulder released in the initial impact, through movement of the primary ejecta through the aperture likely played a role., Comment: 16 pages, 5 Figures, accepted in the Astrophysical Journal Letters (ApJL) on October 16, 2023

Using the absolute detection calibration and abundant detections of the OSSOS (Outer Solar System Origin Survey) project, we provide population measurements for the main Kuiper Belt. For absolute magnitude $H_r<8.3$, there are 30,000 non-resonant main-belt objects, with twice as many hot-component objects than cold, and with total mass of 0.014 $M_\Earth$, only 1/7 of which is in the cold belt (assuming a cold-object albedo about half that of hot component objects). We show that transneptunian objects with $5.5 < H_r < 8.3$ (rough diameters 400--100~km) have indistinguishable absolute magnitude (size) distributions, regardless of being in the cold classical Kuiper belt (thought to be primordial) or the `hot' population (believed to be implanted after having been formed elsewhere). We discuss how this result was not apparent in previous examinations of the size distribution due to the complications of fitting assumed power-law functional forms to the detections at differing depths. This shared size distribution is surprising in light of the common paradigm that the hot population planetesimals formed in a higher density environment much closer to the Sun, in an environment that also (probably later) formed larger (dwarf planet and bigger) objects. If this paradigm is correct, % << BG added clause our result implies that planetesimal formation was relatively insensitive to the local disk conditions and that the subsequent planet-building process in the hot population did not modify the shape of the planetesimal size distribution in this 50--300~km range., Comment: 27 pages, 3 figures

The Milky Way is thought to host a huge population of interstellar objects (ISOs), numbering approximately $10^{15}\mathrm{pc}^{-3}$ around the Sun, which are formed and shaped by a diverse set of processes ranging from planet formation to galactic dynamics. We define a novel framework: firstly to predict the properties of this Galactic ISO population by combining models of processes across planetary and galactic scales, and secondly to make inferences about the processes modelled, by comparing the predicted population to what is observed. We predict the spatial and compositional distribution of the Galaxy's population of ISOs by modelling the Galactic stellar population with data from the APOGEE survey and combining this with a protoplanetary disk chemistry model. Selecting ISO water mass fraction as an example observable quantity, we evaluate its distribution both at the position of the Sun and averaged over the Galactic disk; our prediction for the Solar neighbourhood is compatible with the inferred water mass fraction of 2I/Borisov. We show that the well-studied Galactic stellar metallicity gradient has a corresponding ISO compositional gradient. We also demonstrate the inference part of the framework by using the current observed ISO composition distribution to constrain the parent star metallicity dependence of the ISO production rate. This constraint, and other inferences made with this framework, will improve dramatically as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) progresses and more ISOs are observed. Finally, we explore generalisations of this framework to other Galactic populations, such as that of exoplanets., Comment: Accepted to AJ

The detached transneptunian objects (TNOs) are those with semimajor axes beyond the 2:1 resonance with Neptune, which are neither resonant nor scattering. Using the detached sample from the OSSOS telescopic survey, we produce the first studies of their orbital distribution based on matching the orbits and numbers of the known TNOs after accounting for survey biases. We show that the detached TNO perihelion ($q$) distribution cannot be uniform, but is instead better matched by two uniform components with a break near $q\approx40$ au. We produce parametric two-component models that are not rejectable by the OSSOS data set, and estimate that there are $36,\!000^{+12,000}_{-9,000}$ detached TNOs with absolute magnitudes $H_r < 8.66$ ($D \gtrsim 100$ km) and semimajor axes $48 < a < 250$ au (95% confidence limits). Although we believe these heuristic two-parameter models yield a correct population estimate, we then use the same methods to show that the perihelion distribution of a detached disk created by a simulated rogue planet matches the $q$ distribution even better, suggesting that the temporary presence of other planets in the early Solar System is a promising model to create today's large semimajor axis TNO population. This numerical model results in a detached TNO population estimate of $48,\!000^{+15,000}_{-12,000}$. Because this illustrates how difficult-to-detect $q>50$ au objects are likely present, we conclude that there are $(5 \pm 2)\times10^4$ dynamically detached TNOs, which are thus roughly twice as numerous as the entire transneptunian hot main belt., Comment: Accepted for publication in The Planetary Science Journal. 16 pages, 8 figures

We calculate the upper bounds of the population of theoretically stable Centaur orbits between Uranus and Neptune. These small bodies are on low-eccentricity, low-inclination orbits in two specific bands of semi-major axis, centred at $\sim$24.6 au and $\sim$25.6 au. They exhibit unusually long Gyr-stable lifetimes in previously published numerical integrations, orders of magnitude longer than that of a typical Centaur. Despite the increased breadth and depth of recent solar system surveys, no such objects have been found. Using the Outer Solar System Origins Survey (OSSOS) survey simulator to calculate the detection efficiency for these objects in an ensemble of fully characterised surveys, we determine that a population of 72 stable Centaurs with absolute magnitude $H_{r}\leq10$ ($95\%$ confidence upper limit) could remain undetected. The upcoming Legacy Survey of Space and Time (LSST) will be able to detect this entire intrinsic population due to its complete coverage of the ecliptic plane. If detected, these objects will be interesting dynamically-accessible mission targets -- especially as comparison of the stable Centaur orbital phase space to the outcomes of several modern planetary migration simulations suggests that these objects could be close to primordial in nature., Comment: Accepted to PSJ. 8 pages, 3 figures. Comments welcome

The Vera C. Rubin Observatory is expected to start the Legacy Survey of Space and Time (LSST) in early to mid-2025. This multi-band wide-field synoptic survey will transform our view of the solar system, with the discovery and monitoring of over 5 million small bodies.The final survey strategy chosen for LSST has direct implications on the discoverability and characterization of solar system minor planets and passing interstellar objects. Creating an inventory of the solar system is one of the four main LSST science drivers. The LSST observing cadence is a complex optimization problem that must balance the priorities and needs of all the key LSST science areas. To design the best LSST survey strategy, a series of operation simulations using the Rubin Observatory scheduler have been generated to explore the various options for tuning observing parameters and prioritizations. We explore the impact of the various simulated LSST observing strategies on studying the solar system's small body reservoirs. We examine what are the best observing scenarios and review what are the important considerations for maximizing LSST solar system science. In general, most of the LSST cadence simulations produce +/-5% or less variations in our chosen key metrics, but a subset of the simulations significantly hinder science returns with much larger losses in the discovery and light curve metrics., Comment: Accepted to ApJS, 103 pages (including references), 43 figures, 9 Tables. Videos will be available in the online journal formatted and published version of the paper [v2.0 submission corrects the author list metadata from the arxiv initial submission and updates the abstract]

Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstrated the activation process of an asteroid from an impact under precisely known impact conditions. Here we report the observations of the DART impact ejecta with the Hubble Space Telescope (HST) from impact time T+15 minutes to T+18.5 days at spatial resolutions of ~2.1 km per pixel. Our observations reveal a complex evolution of ejecta, which is first dominated by the gravitational interaction between the Didymos binary system and the ejected dust and later by solar radiation pressure. The lowest-speed ejecta dispersed via a sustained tail that displayed a consistent morphology with previously observed asteroid tails thought to be produced by impact. The ejecta evolution following DART's controlled impact experiment thus provides a framework for understanding the fundamental mechanisms acting on asteroids disrupted by natural impact., Comment: accepted by Nature

The LCO Outbursting Objects Key (LOOK) Project uses the telescopes of the Las Cumbres Observatory (LCO) Network to: (1) to systematically monitor a sample of Dynamically New Comets over the whole sky, and (2) use alerts from existing sky surveys to rapidly respond to and characterize detected outburst activity in all small bodies. The data gathered on outbursts helps to characterize each outburst's evolution with time, assess the frequency and magnitude distribution of outbursts in general, and contributes to the understanding of outburst processes and volatile distribution in the Solar System. The LOOK Project exploits the synergy between current and future wide-field surveys such as ZTF, PanSTARRS, and LSST as well as rapid-response telescope networks such as LCO, and serves as an excellent testbed for what will be needed the much larger number of objects coming from Rubin Observatory. We will describe the LOOK Project goals, the planning and target selection (including the use of NEOexchange as a Target and Observation Manager or "TOM"), and results from the first phase of observations, including the detection of activity and outbursts on the giant comet C/2014 UN271 (Bernardinelli-Bernstein) and the discovery and follow-up of outbursts on comets. Within these outburst discoveries, we present a high cadence of 7P/Pons-Winnecke with days, a large outburst on 57P/duToit-Neujmin-Delporte, and evidence that comet P/2020 X1 (ATLAS) was in outburst when discovered., Comment: 35 pages, 11 figures, accepted for publication in PSJ

In the present-day Kuiper Belt, the number of compositional classes and the orbital distributions of these classes hold important cosmogonic implications for the Solar System. In a companion paper by Fraser et al., we demonstrate that the observed color distribution of small (H>6) Trans-Neptunian Objects (TNOs) can be accounted for by the existence of only two composition classes, named brightIR and faintIR, where the range of colors in each class is governed by a mixture of two material end members. Here, we investigate the orbital distribution of the two color classes identified by Fraser et al. and find that the orbital inclinations of the brightIR class objects are correlated with their optical colors. Using the output of numerical simulations investigating the orbital evolution of TNOs during their scattering phase with Neptune, we show that this correlation could reflect a composition gradient in the early protoplanetary disk, in the range of heliocentric distances over which TNOs from the brightIR class accreted. However, tensions between this interpretation and the existence of blue contaminants among cold classical TNOs, and possible alternative origins for the detected correlation, currently bear uncertainty on our proposed interpretation., Comment: 19 pages, 4 figures, 5 tables. Submitted, comments welcome!

The distribution of surface classes of resonant trans-Neptunian objects (TNOs) provides constraints on the protoplanetesimal disk and giant planet migration. To better understand the surfaces of TNOs, the Colours of the Outer Solar System Origins Survey (Col-OSSOS) acquired multi-band photometry of 102 TNOs, and found that the surfaces of TNOs can be well described by two surface classifications, BrightIR and FaintIR. These classifications both include optically red members and are differentiated predominantly based on whether their near-infrared spectral slope is similar to their optical spectral slope. The vast majority of cold classical TNOs, with dynamically quiescent orbits, have the FaintIR surface classification, and we infer that TNOs in other dynamical classifications with FaintIR surfaces share a common origin with the cold classical TNOs. Comparison between the resonant populations and the possible parent populations of cold classical and dynamically excited TNOs reveal that the 3:2 has minimal contributions from the FaintIR class, which could be explained by the $\nu_8$ secular resonance clearing the region near the 3:2 before any sweeping capture occurred. Conversely, the fraction of FaintIR objects in the 4:3 resonance, 2:1 resonance, and the resonances within the cold classical belt, suggest that the FaintIR surface formed in the protoplanetary disk between 34.6 and 47 au, though the outer bound depends on the degree of resonance sweeping during migration. The presence and absence of the FaintIR surfaces in Neptune's resonances provides critical constraints for the history of Neptune's migration, the evolution of the $\nu_8$, and the surface class distribution in the initial planetesimal disk, Comment: 19 pages, 8 figures. in Press at PSJ

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

Dynamically excited objects within the Kuiper belt show a bimodal distribution in their surface colors, and these differing surface colors may be a tracer of where these objects formed. In this work we explore radial color distributions in the primordial planetesimal disk and implications for the positions of ice line/color transitions within the Kuiper belt's progenitor populations. We combine a full dynamical model of the Kuiper belt's evolution due to Neptune's migration with precise surface colors measured by the Colours of the Outer Solar System Origins Survey in order to examine the true color ratios within the Kuiper belt and the ice lines within the primordial disk. We investigate the position of a dominant, surface color changing ice-line, with two possible surface color layouts within the initial disk; (1) inner neutral surfaces and outer red, and (2) inner red surfaces and outer neutral. We performed simulations with a primordial disk that truncates at 30 au. By radially stepping the color transition out through 0.5 au intervals we show that both disk configurations are consistent with the observed color fraction. For an inner neutral, outer red primordial disk we find that the color transition can be at $28^{+2}_{-3}$ au at a 95% confidence level. For an inner red, outer neutral primordial disk the color transition can be at $27^{+3}_{-3}$ au at a 95% confidence level., Comment: 18 pages, 8 figures, accepted for publication in PSJ

Planetesimals inevitably bear the signatures of their natal environment, preserving in their composition a record of the metallicity of their system's original gas and dust, albeit one altered by the formation process. When planetesimals are dispersed from their system of origin, this record is carried with them. As each star is likely to contribute at least $10^{12}$ interstellar objects, the Galaxy's drifting population of interstellar objects (ISOs) provides an overview of the properties of its stellar population through time. Using the EAGLE cosmological simulation and models of protoplanetary formation, our modelling predicts an ISO population with a bimodal distribution in their water mass fraction. Objects formed in low-metallicity, typically older, systems have a higher water fraction than their counterparts formed in high-metallicity protoplanetary disks, and these water-rich objects comprise the majority of the population. Both detected ISOs seem to belong to the lower water fraction population; these results suggest they come from recently formed systems. We show that the population of ISOs in galaxies with different star formation histories will have different proportions of objects with high and low water fractions. This work suggests that it is possible that the upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time will detect a large enough population of ISOs to place useful constraints on models of protoplanetary disks, as well as galactic structure and evolution., Comment: Accepted by ApJ Letters, 9 pages

We present the discovery of 2013 VZ70, the first known horseshoe coorbital companion of Saturn. Observed by the Outer Solar System Origins Survey (OSSOS) for 4.5 years, the orbit of 2013 VZ70 is determined to high precision, revealing that it currently is in `horseshoe' libration with the planet. This coorbital motion will last at least thousands of years but ends ~10 kyr from now; 2013 VZ70 is thus another example of the already-known `transient coorbital' populations of the giant planets, with this being the first known prograde example for Saturn (temporary retrograde coorbitals are known for Jupiter and Saturn). We present a theoretical steady state model of the scattering population of trans-Neptunian origin in the giant planet region (2--34 au), including the temporary coorbital populations of the four giant planets. We expose this model to observational biases using survey simulations in order to compare the model to the real detections made by a set of well-characterized outer Solar System surveys. While the observed number of coorbitals relative to the scattering population is higher than predicted, we show that the number of observed transient coorbitals of each giant planet relative to each other is consistent with a transneptunian source., Comment: 27 pages (double-line spaced manuscript format), 4 figures, 3 tables. Published in The Planetary Science Journal, 2:212 (11pp), 2021 October. Open Access journal article available at https://iopscience.iop.org/article/10.3847/PSJ/ac1c6b

The cold main classical Kuiper Belt consists of those non-resonant small solar system bodies with low orbital inclinations and orbital semi-major axes between 42.4 and 47.7 au. These objects likely formed \textit{in situ} and the population has experienced minimal collisional modification since formation. Using the Outer Solar System Origins Survey (OSSOS) ensemble sample and characterization, combined with constraints from deeper surveys and supported by evidence from the Minor Planet Center catalog and the Deep Ecliptic Survey, we determine the absolute magnitude $H_r$ distribution of the cold classical belt from $H_r\simeq5$ to 12 (roughly diameters of 400 km to 20 km). We conclude that the cold population's $H_r$ distribution exhibits an exponential cutoff at large sizes. Exponential cutoffs at large sizes are not a natural outcome of pair-wise particle accretion but exponentially tapered power-law size distributions are a feature of numerical simulations of planetesimal formation via a streaming instability. Our observation of an exponential cutoff agrees with previous observational inferences that no large objects ($D \gtrsim 400$ km) exist in the cold population. We note that the asymptotic slope of the $H_r$ distribution is consistent with $\alpha \sim 0.4$ and this asymptotic slope is also found in streaming instability modelling of planetesimal formation and is thus not necessarily associated with achieving collisional equilibrium. Studies of the transneptunian region are providing the parameters that will enable future streaming-instability studies to determine the initial conditions of planetesimal formation in the $\approx$45 au region of the Sun's protoplanetary disk., Comment: 20 pages, 2 figures, accepted in AAS Journals

Interstellar objects (ISOs), the parent population of 1I/Oumuamua and 2I/Borisov, are abundant in the interstellar medium of the Milky Way. This means that the interstellar medium, including molecular cloud regions, has three components: gas, dust, and ISOs. From the observational constraints for the field density of ISOs drifting in the solar neighbourhood, we infer a typical molecular cloud of 10 pc diameter contains some 10$^{18}$ ISOs. At typical sizes ranging from hundreds of metres to tens of km, ISOs are entirely decoupled from the gas dynamics in these molecular clouds. Here we address the question of whether ISOs can follow the collapse of molecular clouds. We perform low-resolution simulations of the collapse of molecular clouds containing initially static ISO populations toward the point where stars form. In this proof-of-principle study, we find that the interstellar objects definitely follow the collapse of the gas -- and many become bound to the new-forming numerical approximations to future stars (sinks). At minimum, 40\% of all sinks have one or more ISO test particles gravitationally bound to them for the initial ISO distributions tested here. This value corresponds to at least $10^{10}$ actual interstellar objects being bound after three initial free-fall times. Thus, ISOs are a relevant component of star formation. We find that more massive sinks bind disproportionately large fractions of the initial ISO population, implying competitive capture of ISOs. Sinks can also be solitary, as their ISOs can become unbound again -- particularly if sinks are ejected from the system. Emerging planetary systems will thus develop in remarkably varied environments, ranging from solitary to richly populated with bound ISOs., Comment: 15 pages, 7 figures, accepted for ApJ

Here, we present results on the intrinsic collision probabilities, $ P_I$, and range of collision speeds, $V_I$, as a function of the heliocentric distance, $r$, in the trans-Neptunian region. The collision speed is one of the parameters, that serves as a proxy to a collisional outcome e.g., complete disruption and scattering of fragments, or formation of crater, where both processes are directly related to the impact energy. We utilize an improved and de-biased model of the trans-Neptunian object (TNO) region from the "Outer Solar System Origins Survey" (OSSOS). It provides a well-defined orbital distribution model of TNOs, based on multiple opposition observations of more than 1000 bodies. In this work we compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using the \"{O}pik's approach, as revised and modified by Wetherill for non-circular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in ($r, V_I$) phase space. Although collisions are calculated to happen within $r\sim 20 - 200$~AU and $V_I \sim 0.1$~km/s to as high as $V_I\sim9$~km/s, most of the collisions are likely to happen at low relative velocities $V_I<1$~km/s and are dominated by the main classical belt., Comment: 13 pages, 6 figures

The cold classical Kuiper Belt Objects (KBOs) possess a high, $\gtrsim30\%$ binary fraction. Widely separated and dynamically fragile, these binary systems have been useful in tracing the origins of KBOs. A new class of binaries was recently identified by their colours. The so-called blue binaries are unanimously members of the less red compositional class, and exhibit a 100% binary fraction. They appear to be push-out survivors, emplaced in the classical region during Neptune's phases of outward migration. The presence of these binary systems implies that the majority of objects that formed near the cold classical region formed as binaries. Here we present new optical colour measurements of cold classical KBOs from the Colours of the Outer Solar System Origins Survey, including colours of a blue binary discovered by the Solar System Origins Legacy Survey -- 2015 RJ277. The increased size of the colours sample has resulted in order-of-magnitude decrease in the probability that the binaries and singles sample share the same colour distribution. From the Anderson-Darling statistic, this probability is only a 0.3%, while it is only 0.002% when utilizing the difference of means statistic. We find a hint that the blue binaries have inflated free inclinations compared to their red counterparts, consistent with the push-out origin for these bodies., Comment: Accepted for publication in the PSJ. 4 Figures

The detections of 1I/Oumuamua and 2I/Borisov within just two years demonstrate impressively that interstellar objects (ISOs) must be common in the Milky Way. Once released from their parent system, these ISOs travel for Gyr through interstellar space. While often imagined as empty, interstellar space contains gas and dust most prominent in the form of molecular clouds. Performing numerical simulations, we test how often ISOs cross such molecular clouds. We find that the ISOs pass amazingly often through molecular clouds. In the solar neighbourhood, ISOs typically spend 0.1-0.2% of their journey inside molecular clouds, for relative slow ISOs ($<$ 5 km/s) this can increase to 1-2%, equivalent to 10 - 20 Myr per Gyr. Thus the dynamically youngest ISOs spend the longest time in molecular clouds. In other words, molecular clouds must mainly contain relatively young ISOs ($<$ 1-2 Gyr). Thus the half-life of the seeding process by ISOs is substantially shorter than a stellar lifetime. The actual amount of time spent in MCs decreases with distance to the Galactic Centre. We find that ISOs pass so often through MCs that backtracing their path to find their parent star beyond 250 Myr seems beyond the point. Besides, we give a first estimate of the ISO density depending on the galactic centre distance based on the stellar distribution., Comment: 12 pages, 7 figures, accepted by ApJ

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 members and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range., Comment: 15 pages, 7 figures, v3, accepted for publication in Icarus

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 42r*, with 30

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., Comment: 10 pages, 4 figures, accepted to The Planetary Science Journal

We report the clear detection of C$_2$ and of abundant NH$_2$ in the first prominently active interstellar comet, 2I/Borisov. We observed 2I on three nights in November 2019 at optical wavelengths 4800--9300 \AA with the Multi-Unit Spectroscopic Explorer (MUSE) integral-field spectrograph on the ESO/Very Large Telescope. These data, together with observations close in time from both 0.6-m TRAPPIST telescopes, provide constraints on the production rates of species of gas in 2I's coma. From the MUSE detection on all epochs of several bands of the optical emission of the C$_2$ Swan system, a rich emission spectrum of NH$_2$ with many highly visible bands, and the red (1-0) bandhead of CN, together with violet CN detections by TRAPPIST, we infer production rates of $Q$(C$_2$) = $1.1\times10^{24}$ mol s$^{-1}$, $Q$(NH$_2$) = $4.8\times10^{24}$ mol s$^{-1}$ and $Q$(CN) = $(1.8\pm0.2)\times 10^{24}$ mol s$^{-1}$. In late November at 2.03~au, 2I had a production ratio of C$_2$/CN$=0.61$, only barely carbon-chain depleted, in contrast to earlier reports measured further from the Sun of strong carbon-chain depletion. Thus, 2I has shown evolution in its C$_2$ production rate: a parent molecule reservoir has started sublimating. At $Q$(NH$_2$)/$Q$(CN) = 2.7, this second interstellar object is enriched in NH$_2$, relative to the known Solar System sample., Comment: 5 figures. Submitted to AAS Journals

We observed the interstellar comet 2I/Borisov using the Neil Gehrels-Swift Observatory's Ultraviolet/Optical Telescope. We obtained images of the OH gas and dust surrounding the nucleus at six epochs spaced before and after perihelion (-2.56 AU to 2.54 AU). Water production rates increased steadily before perihelion from $(7.0\pm1.5)\times10^{26}$ molecules s$^{-1}$ on Nov. 1, 2019 to $(10.7\pm1.2)\times10^{26}$ molecules s$^{-1}$ on Dec. 1. This rate of increase in water production rate is quicker than that of most dynamically new comets and at the slower end of the wide range of Jupiter-family comets. After perihelion, the water production rate decreased to $(4.9\pm0.9)\times10^{26}$ molecules s$^{-1}$ on Dec. 21, which is much more rapidly than that of all previously observed comets. Our sublimation model constrains the minimum radius of the nucleus to 0.37 km, and indicates an active fraction of at least 55% of the surface. $A(0)f\rho$ calculations show a variation between 57.5 and 105.6 cm with a slight trend peaking before the perihelion, lower than previous and concurrent published values. The observations confirm that 2I/Borisov is carbon-chain depleted and enriched in NH$_2$ relative to water., Comment: 12 pages, 3 figures, 2 tables, submitted to ApJL

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., Comment: 23 pages, 12 figures, 2 tables, accepted for publication in the Astronomical Journal

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., Comment: 15 pages, 3 figures, to appear in The Trans-Neptunian Solar System (Chapter 3)

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., Comment: 11 pages, 3 figures, 3 tables

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 5010 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

Uranus and Neptune, and their diverse satellite and ring systems, represent the least explored environments of our Solar System, and yet may provide the archetype for the most common outcome of planetary formation throughout our galaxy. Ice Giants will be the last remaining class of Solar System planet to have a dedicated orbital explorer, and international efforts are under way to realise such an ambitious mission in the coming decades. In 2019, the European Space Agency released a call for scientific themes for its strategic science planning process for the 2030s and 2040s, known as Voyage 2050. We used this opportunity to review our present-day knowledge of the Uranus and Neptune systems, producing a revised and updated set of scientific questions and motivations for their exploration. This review article describes how such a mission could explore their origins, ice-rich interiors, dynamic atmospheres, unique magnetospheres, and myriad icy satellites, to address questions at the heart of modern planetary science. These two worlds are superb examples of how planets with shared origins can exhibit remarkably different evolutionary paths: Neptune as the archetype for Ice Giants, whereas Uranus may be atypical. Exploring Uranus' natural satellites and Neptune's captured moon Triton could reveal how Ocean Worlds form and remain active, redefining the extent of the habitable zone in our Solar System. For these reasons and more, we advocate that an Ice Giant System explorer should become a strategic cornerstone mission within ESA's Voyage 2050 programme, in partnership with international collaborators, and targeting launch opportunities in the early 2030s., Comment: 34 pages, 9 figures, accepted to Planetary and Space Science

The Large Synoptic Survey Telescope (LSST) is expected to increase known small solar system object populations by an order of magnitude or more over the next decade, enabling a broad array of transformative solar system science investigations to be performed. In this white paper, we discuss software tools and infrastructure that we anticipate will be needed to conduct these investigations and outline possible approaches for implementing them. Feedback from the community or contributions to future updates of this work are welcome. Our aim is for this white paper to encourage further consideration of the software development needs of the LSST solar system science community, and also to be a call to action for working to meet those needs in advance of the expected start of the survey in late 2022., Comment: 82 pages

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

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more., Comment: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

The discovery of 1I/`Oumuamua confirmed that planetesimals must exist in great numbers in interstellar space. Originally generated during planet formation, they are scattered from their original systems and subsequently drift through interstellar space. As a consequence they should seed molecular clouds with at least hundred-metre-scale objects. We consider how the galactic background density of planetesimals, enriched from successive generations of star and system formation, can be incorporated into forming stellar systems. We find that at minimum of the order of 10$^{7}$ `Oumuamua-sized and larger objects, plausibly including hundred-kilometre-scale objects, should be present in protoplanetary disks. At such initial sizes, the growth process of these seed planetesimals in the initial gas- and dust-rich protoplanetary disks is likely to be substantially accelerated. This could resolve the tension between accretionary timescales and the observed youth of fully-fledged planetary systems. Our results strongly advocate that the population of interstellar planetesimals should be taken into account in future studies of planet formation. As not only the Galaxy's stellar metallicity increased over time but also the density of interstellar objects, we hypothesize that this enriched seeding accelerates and enhances planetary formation after the first couple of generations of planetary systems., Comment: Accepted to ApJL. Conversation welcome

A foundational goal of the Large Synoptic Survey Telescope (LSST) is to map the Solar System small body populations that provide key windows into understanding of its formation and evolution. This is especially true of the populations of the Outer Solar System -- objects at the orbit of Neptune $r > 30$AU and beyond. In this whitepaper, we propose a minimal change to the LSST cadence that can greatly enhance LSST's ability to discover faint distant Solar System objects across the entire wide-fast-deep (WFD) survey area. Specifically, we propose that the WFD cadence be constrained so as to deliver least one sequence of $\gtrsim 10$ visits per year taken in a $\sim 10$ day period in any combination of $g, r$, and $i$ bands. Combined with advanced shift-and-stack algorithms (Whidden et al. 2019) this modification would enable a nearly complete census of the outer Solar System to $\sim 25.5$ magnitude, yielding $4-8$x more KBO discoveries than with single-epoch baseline, and enabling rapid identification and follow-up of unusual distant Solar System objects in $\gtrsim 5$x greater volume of space. These increases would enhance the science cases discussed in Schwamb et al. (2018) whitepaper, including probing Neptune's past migration history as well as discovering hypothesized planet(s) beyond the orbit of Neptune (or at least placing significant constraints on their existence)., Comment: White Paper submitted in response to the Call for LSST Cadence Optimization White Papers

We propose an ecliptic Deep Drilling Field that will discover some 10,000~small and faint Kuiper Belt Objects (KBOs) --- primitive rocky/icy bodies that orbit at the outside of our Solar System and uniquely record the processes of planetary system formation and evolution. The primary goals are to measure the KBO size and shape distributions down to 25~km, a size that probes both the early and ongoing evolution of this population. These goals can be met with around 10~hours total of on-sky time (five separate fields that are observed for 2.1~hours each). Additional science will result from downstream observations that provide colors and orbit refinement, for a total time request of 40~hours over the ten year LSST main survey., Comment: White Paper submitted in response to the Call for LSST Cadence Optimization White Papers

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., Comment: 49 pages, 14 figures, 7 tables

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., Comment: 24 pages, 5 figures, accepted for publication in AJ

Making an inventory of the Solar System is one of the four fundamental science requirements for the Large Synoptic Survey Telescope (LSST). The current baseline footprint for LSST's main Wide-Fast-Deep (WFD) Survey observes the sky below 0$^\circ$ declination, which includes only half of the ecliptic plane. Critically, key Solar System populations are asymmetrically distributed on the sky: they will be entirely missed, or only partially mapped, if only the WFD occurs. We propose a Northern Ecliptic Spur (NES) mini survey, observing the northern sky up to +10$^\circ$ ecliptic latitude, to maximize Solar System science with LSST. The mini survey comprises a total area of $\sim$5800 deg$^2$/604 fields, with 255 observations/field over the decade, split between g,r, and z bands. Our proposed survey will 1) obtain a census of main-belt comets; 2) probe Neptune's past migration history, by exploring the resonant structure of the Kuiper belt and the Neptune Trojan population; 3) explore the origin of Inner Oort cloud objects and place significant constraints on the existence of a hypothesized planet beyond Neptune; and 4) enable precise predictions of KBO stellar occultations. These high-ranked science goals of the Solar System Science Collaboration are only achievable with this proposed northern survey., Comment: White Paper submitted in response to the Call for LSST Cadence Optimization White Papers

Making an inventory of the Solar System is one of the four pillars that the requirements for the Large Synoptic Survey Telescope (LSST) are built upon. The choice between same-filter nightly pairs or different-filter nightly pairs in the Wide-Fast-Deep (WFD) Survey will have a dramatic effect on the ability of the Moving Object Pipeline System (MOPS) to detect certain classes of Solar System objects; many of the possible filter pairings would result in significant ($\sim50\%$ or more) loss of Solar System object detections. In particular, outer Solar System populations can be significantly redder than those in the inner Solar System, and nightly pairs in $r$-band will result in the deepest survey for the outer Solar System. To maximize the potential for outer Solar System science, we thus advocate for ensuring that the WFD survey contains a sufficient number of $r$-$r$ nightly pairs for each field during a discovery season to ensure detection and linking using MOPS. We also advocate for adding additional spectral energy distributions (SEDs) that more accurately model outer Solar System populations to the pipeline for evaluating the outputs of the LSST operations simulator. This will enable a better estimate of how many Solar System population detections are lost or gained for different filter choices in the WFD survey., Comment: LSST Cadence Optimization White Paper

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., Comment: Accepted to ApJS; on-line supplemental files will be available with the AJS published version of the paper

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., Comment: Invited paper, special issue Data: Insights and Challenges in a Time of Abundance. Data tables and example survey simulator are in the supplementary materials (see arXiv source under Downloads > Other formats)