Search

The Vera C. Rubin Observatory is due to commence the 10-year Legacy Survey of Space and Time (LSST) at the end of 2025. To detect transient/variable sources and identify solar system objects (SSOs), the processing pipelines require templates of the static sky to perform difference imaging. During the first year of the LSST, templates must be generated as the survey progresses, otherwise SSOs cannot be discovered nightly. The incremental template generation strategy has not been finalized; therefore, we use the Metric Analysis Framework (MAF) and a simulation of the survey cadence (one_snap_v4.0_10yrs}) to explore template generation in Year 1. We have assessed the effects of generating templates over timescales of days-weeks, when at least four images of sufficient quality are available for $\geq90\%$ of the visit. We predict that SSO discoveries will begin $\sim$2-3 months after the start of the survey. We find that the ability of the LSST to discover SSOs in real-time is reduced in Year 1. This is especially true for detections in areas of the sky that receive fewer visits, such as the North Ecliptic Spur (NES), and in less commonly used filters, such as the $u$ and $g$-bands. The lack of templates in the NES dominates the loss of real-time SSO discoveries; across the whole sky the MAF Main-Belt asteroid (MBA) discovery metric decreases by up to $63\%$ compared to the baseline observing strategy, whereas the metric decreases by up to $79\%$ for MBAs in the NES alone., Comment: Submitted for consideration in the ApJS Rubin LSST Survey Strategy Optimization Focus Issue

The Legacy Survey of Space and Time (LSST) at Vera C. Rubin Observatory is planned to begin in the Fall of 2025. The LSST survey cadence has been designed via a community-driven process regulated by the Survey Cadence Optimization Committee (SCOC), which recommended up to 3% of the observing time to carry out Target of Opportunity (ToO) observations. Experts from the scientific community, Rubin Observatory personnel, and members of the SCOC were brought together to deliver a recommendation for the implementation of the ToO program during a workshop held in March 2024. Four main science cases were identified: gravitational wave multi-messenger astronomy, high energy neutrinos, Galactic supernovae, and small potentially hazardous asteroids possible impactors. Additional science cases were identified and briefly addressed in the documents, including lensed or poorly localized gamma-ray bursts and twilight discoveries. Trigger prioritization, automated response, and detailed strategies were discussed for each science case. This document represents the outcome of the Rubin ToO 2024 workshop, with additional contributions from members of the Rubin Science Collaborations. The implementation of the selection criteria and strategies presented in this document has been endorsed in the SCOC Phase 3 Recommendations document (PSTN-056). Although the ToO program is still to be finalized, this document serves as a baseline plan for ToO observations with the Rubin Observatory.

We utilize serendipitous observations from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the Zwicky Transient Facility (ZTF) in addition to targeted follow-up observations from the Las Cumbres Observatory (LCO) and Liverpool Telescope to analyze the first observed instance of cometary activity by the newly-discovered Jupiter-family comet C/2023 RN3 (ATLAS), whose orbital dynamics place it close to residing on a Centaur-like orbit. Across our 7-month baseline, we observe an epoch of cometary activity commencing in August 2023 with an increase in brightness of >5.4 mag. The lightcurve of 2023 RN3 indicates the presence of continuous cometary activity across our observations, suggesting the onset of a new period of sustained activity. We find no evidence of any outbursts on top of the observed brightening, nor do we find any significant color evolution across our observations. 2023 RN3 is visibly extended in LCO and Liverpool Telescope observations, indicating the presence of a spatially-extended coma. Numerical integration of 2023 RN3's orbit reveals the comet to have recently undergone a slight increase in semimajor axis due to a planetary encounter with Jupiter, however whether this orbital change could trigger 2023 RN3's cometary activity is unclear. Our estimate for the maximum dust production metric of Afrho ~400 cm is consistent with previous measurements for the Jupiter-family comet and Centaur populations., Comment: 20 pages, 9 figures

Centaurs are small Solar System objects on chaotic orbits in the giant planet region, forming an evolutionary continuum with the Kuiper belt objects and Jupiter-family comets. Some Centaurs are known to exhibit cometary activity, though unlike comets this activity tends not to correlate with heliocentric distance and the mechanism behind it is currently poorly understood. We utilize serendipitous observations from the Asteroid Terrestrial-impact Last Alert System (ATLAS), Zwicky Transient Facility (ZTF), Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), Dark Energy Survey (DES), and Gaia in addition to targeted follow-up observations from the Las Cumbres Observatory, TRAnsiting Planets and PlanetesImals Small Telescope South (TRAPPIST-South), and Gemini North telescope to analyze an unexpected brightening exhibited by the known active Centaur (2060) Chiron in 2021. This is highly indicative of a cometary outburst. As of 2023 February, Chiron has still not returned to its pre-brightening magnitude. We find Chiron's rotational lightcurve, phase curve effects, and possible high-albedo surface features to be unlikely causes of this observed brightening. We consider the most likely cause to be an epoch of either new or increased cometary activity, though we cannot rule out a possible contribution from Chiron's reported ring system, such as a collision of as-yet unseen satellites shepherding the rings. We find no evidence for coma in our Gemini or TRAPPIST-South observations, though this does not preclude the possibility that Chiron is exhibiting a coma that is too faint for observation or constrained to the immediate vicinity of the nucleus., Comment: 39 pages, 14 figures, 14 tables. Has been accepted for publication in PSJ

On 26 September 2022, the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, the satellite of binary near-Earth asteroid (65803) Didymos. This demonstrated the efficacy of a kinetic impactor for planetary defense by changing the orbital period of Dimorphos by 33 minutes (Thomas et al. 2023). Measuring the period change relied heavily on a coordinated campaign of lightcurve photometry designed to detect mutual events (occultations and eclipses) as a direct probe of the satellite's orbital period. A total of 28 telescopes contributed 224 individual lightcurves during the impact apparition from July 2022 to February 2023. We focus here on decomposable lightcurves, i.e. those from which mutual events could be extracted. We describe our process of lightcurve decomposition and use that to release the full data set for future analysis. We leverage these data to place constraints on the post-impact evolution of ejecta. The measured depths of mutual events relative to models showed that the ejecta became optically thin within the first ~1 day after impact, and then faded with a decay time of about 25 days. The bulk magnitude of the system showed that ejecta no longer contributed measurable brightness enhancement after about 20 days post-impact. This bulk photometric behavior was not well represented by an HG photometric model. An HG1G2 model did fit the data well across a wide range of phase angles. Lastly, we note the presence of an ejecta tail through at least March 2023. Its persistence implied ongoing escape of ejecta from the system many months after DART impact., Comment: 52 pages, 5 tables, 9 figures, accepted to PSJ

Jovian co-orbitals share Jupiter's orbit in 1:1 mean motion resonance. This includes $>$10,000 so-called Trojan asteroids surrounding the leading (L4) and trailing (L5) Lagrange points, viewed as stable groups dating back to planet formation. Via a massive numerical study we identify for the first time some Trojans which are certainly only `metastable'; instead of being primordial, they are recent captures from heliocentric orbits into moderately long-lived (10 kyr - 100 Myr) metastable states that will escape back to the scattering regime. We have also identified (1) the first two jovian horseshoe co-orbitals that exist for many resonant libration periods, and (2) eight jovian quasi-satellites with metastable lifetimes of 4-130 kyr. Our perspective on the Trojan population is thus now more complex as Jupiter joins the other giant planets in having known metastable co-orbitals which are in steady-state equilibrium with the planet-crossing Centaur and asteroid populations, in agreement with theoretical estimates.

Long-period comet C/2018 F4 (PANSTARRS) was observed to show duplicity of its inner region in 2020 September, suggestive of a splitting event. We here present analyses of our observations of the comet taken from the LOOK project and the University of Hawaii 2.2 m telescope after the discovery of the splitting. The two fragments Components A and B, estimated to be $\sim\!60$ m to 4 km in radius, remained highly similar to each other in terms of brightness, colour, and dust morphology throughout our observing campaign from 2020 September to 2021 December. Our fragmentation model yielded that the two components split at a relative speed of $3.00 \pm 0.18$ m s$^{-1}$ in 2020 late April, implying a specific energy change of $\left(5.3 \pm 2.8 \right) \times 10^3$ J kg$^{-1}$, and that Component B was subjected to a stronger nongravitational acceleration than Component A in both the radial and normal directions of the orbit. The obtained splitting time is broadly consistent with the result from the dust morphology analysis, which further suggested that the dominant dust grains were millimeter-sized and ejected at speed $\sim\!2$ m s$^{-1}$. We postulate that the pre-split nucleus of the comet consisted of two lobes resembling the one of 67P, or that the comet used to be a binary system like main-belt comet 288P. Regardless, we highlight the possibility of using observations of split comets as a feasible manner to study the bilobate shape or binarity fraction of cometary nuclei., Comment: Accepted to AJ for publication

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]

The Double Asteroid Redirection Test (DART) spacecraft successfully performed the first test of a kinetic impactor for asteroid deflection by impacting Dimorphos, the secondary of near-Earth binary asteroid (65803) Didymos, and changing the orbital period of Dimorphos. A change in orbital period of approximately 7 minutes was expected if the incident momentum from the DART spacecraft was directly transferred to the asteroid target in a perfectly inelastic collision, but studies of the probable impact conditions and asteroid properties indicated that a considerable momentum enhancement ($\beta$) was possible. In the years prior to impact, we used lightcurve observations to accurately determine the pre-impact orbit parameters of Dimorphos with respect to Didymos. Here we report the change in the orbital period of Dimorphos as a result of the DART kinetic impact to be -33.0 +/- 1.0 (3$\sigma$) minutes. Using new Earth-based lightcurve and radar observations, two independent approaches determined identical values for the change in the orbital period. This large orbit period change suggests that ejecta contributed a significant amount of momentum to the asteroid beyond what the DART spacecraft carried., Comment: Accepted by Nature

Cometary activity may be driven by ices with very low sublimation temperatures, such as carbon monoxide ice, which can sublimate at distances well beyond 20 au. This point is emphasized by the discovery of Oort cloud comet C/2014 UN$_{271}$ (Bernardinelli-Bernstein), and its observed activity out to $\sim$26 au. Through observations of this comet's optical brightness and behavior, we can potentially discern the drivers of activity in the outer solar system. We present a study of the activity of comet Bernardinelli-Bernstein with broad-band optical photometry taken at 19-20 au from the Sun (2021 June to 2022 February) as part of the LCO Outbursting Objects Key (LOOK) Project. Our analysis shows that the comet's optical brightness during this period was initially dominated by cometary outbursts, stochastic events that ejected $\sim10^7$ to $\sim10^8$ kg of material on short (< 1 day) timescales. We present evidence for three such outbursts occurring in 2021 June and September. The nominal nuclear volumes excavated by these events are similar to the 10-100 m pit-shaped voids on the surfaces of short-period comet nuclei, as imaged by spacecraft. Two out of three Oort cloud comets observed at large pre-perihelion distances exhibit outburst behavior near 20 au, suggesting such events may be common in this population. In addition, quiescent CO-driven activity may account for the brightness of the comet in 2022 January to February, but that variations in the cometary active area (i.e., the amount of sublimating ice) with heliocentric distance are also possible., Comment: Accepted for publication in The Astrophysical Journal Letters 13 pages, 6 figures, 1 table

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

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

We report results from new and archival observations of the newly discovered active asteroid (248370) 2005 QN_137, which has been determined to be a likely main-belt comet based on a subsequent discovery that it is recurrently active near perihelion. From archival data analysis, we estimate g'-, r'-, i'-, and z'-band absolute magnitudes for the nucleus of H_g=16.62+/-0.13, H_r=16.12+/-0.10, H_i=16.05+/-0.11, and H_z=15.93+/-0.08, corresponding to nucleus colors of g'-r'=0.50+/-0.16, r'-i'=0.07+/-0.15, and i'-z'=0.12+/-0.14, an equivalent V-band absolute magnitude of H_V=16.32+/-0.08, and a nucleus radius of r_n=1.6+/-0.2 km (using a V-band albedo of p_V=0.054+/-0.012). Meanwhile, we find mean near-nucleus coma colors when 248370 was active of g'-r'=0.47+/-0.03, r'-i'=0.10+/-0.04, and i'-z'=0.05+/-0.05, and similar mean dust tail colors, suggesting that no significant gas coma is present. We find approximate ratios between the scattering cross-sections of near-nucleus dust (within 5000 km of the nucleus) and the nucleus of A_d/A_n=0.7+/-0.3 on 2016 July 22, and 1.8

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

Vera C. Rubin Observatory will be a key facility for small body science in planetary astronomy over the next decade. It will carry out the Legacy Survey of Space and Time (LSST), observing the sky repeatedly in u, g, r, i, z, and y over the course of ten years using a 6.5 m effective diameter telescope with a 9.6 square degree field of view, reaching approximately r = 24.5 mag (5-{\sigma} depth) per visit. The resulting dataset will provide extraordinary opportunities for both discovery and characterization of large numbers (10--100 times more than currently known) of small solar system bodies, furthering studies of planetary formation and evolution. This white paper summarizes some of the expected science from the ten years of LSST, and emphasizes that the planetary astronomy community should remain invested in the path of Rubin Observatory once the LSST is complete., Comment: White paper submitted to the 2020 Planetary Astronomy Decadal Survey (7 pages, 1 figure)

Pluto's terrains display a diversity of crater retention ages ranging from areas with no identifiable craters to heavily cratered terrains. This variation in crater densities is consistent with geologic activity occurring throughout Pluto's history and also a variety of resurfacing styles, including both exogenic and endogenic processes. Using estimates of impact flux and cratering rates over time, Pluto's heavily cratered terrains appear to be relatively ancient, 4 Ga or older. Charon's smooth plains, informally named Vulcan Planitia, did experience early resurfacing, but there is a relatively high spatial density of craters on Vulcan Planitia and almost all overprint the other types of volcanic or tectonic features. Both Vulcan Planitia and the northern terrains on Charon are also estimated to be ancient, 4 Ga or older. The craters on Pluto and Charon also show a distinct break in their size-frequency distributions (SFDs), where craters smaller than approximately 10-15 km in diameter have a shallower SFD power-law slope than those larger than this break diameter. This SFD shape on Pluto and Charon is different than what is observed on the Earth's Moon, and gives the Kuiper belt impactor SFD a different shape than that of the asteroid belt., Comment: 44 pages, 14 figures, 2 tables

We demonstrate dynamical pathways from main-belt asteroid and Centaur orbits to those in co-orbital motion with Jupiter, including the retrograde (inclination $i>90^o$) state. We estimate that at any given time, there should be $\sim1$ kilometer-scale or larger escaped asteroid in a transient direct (prograde) orbit with semimajor axis near that of Jupiter's ($a\simeq a_J$), with proportionally more smaller objects as determined by their size distribution. Most of these objects would be in the horseshoe dynamical state, which are hard to detect due to their moderate eccentricities (spending most of their time beyond 5 AU) and longitudes relative to Jupiter being spread nearly all over the sky. We also show that $\approx$1% of the transient asteroid co-orbital population is on retrograde orbits with Jupiter. This population, like the recently identified asteroid (514107) 2015 BZ$_{509}$, can spend millions of years with $a\simeq a_J$ including tens or hundreds of thousands of years formally in the retrograde 1:-1 co-orbital resonance. Escaping near-Earth asteroids (NEAs) are thus likely the precursors to the handful of known high-inclination objects with $a\simeq a_J$. We compare the production of jovian co-orbitals from escaping NEAs with those from incoming Centaurs. We find that temporary direct co-orbitals are likely dominated by Centaur capture, but we only find production of (temporary) retrograde jovian co-orbitals (including very long-lived ones) from the NEA source. We postulate that the primordial elimination of the inner Solar System's planetesimal population could provide a supply route for a metastable outer Solar System reservoir for the high-inclination Centaurs., Comment: Accepted at AJ

Vatira-class near-Earth objects (NEOs) have orbits entirely interior to the orbit of Venus with aphelia $0.307

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

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., Comment: V2. Incorrect Fig 1 altered. Submitted to ApJ Letters Dec 21/2018. 4 figures, 1 Table

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)

Given rapid photodissociation and photodegradation, the recently discovered organics in the Martian subsurface and atmosphere were probably delivered in geologically recent times. Possible parent bodies are C-type asteroids, comets, and interplanetary dust particles (IDPs). The dust infall rate was estimated, using different methods, to be between $0.71$ and $2.96 \times 10^6$ kg/yr (Nesvorny et al., 2011, Borin et al., 2017, Crismani et al., 2017); assuming a carbon content of 10% (Flynn, 1996), this implies an IDP carbon flux of $0.07 - 0.3 \times 10^6$ kg/yr. We calculate for the first time the carbon flux from impacts of asteroids and comets. To this end, we perform dynamical simulations of impact rates on Mars. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We separately add comets and asteroids to the simulations as massless test particles, based on their current orbital elements, yielding Mars impact rates of $4.34\times10^{-3}$ comets/Myr and 3.3 asteroids/Myr. We estimate the global carbon flux on Mars from cometary impacts to be $\sim 0.013 \times 10^{6}$~kg/yr within an order of magnitude, while asteroids deliver $\sim 0.05 \times 10^6$~kg/yr. These values correspond to $\sim 4-19 \%$ and $\sim 17-71 \%$, respectively, of the IDP-borne carbon flux estimated by Nesvorny et al. 2011, Borin et al. 2017 and Crismani et al. 2017. Unlike the spatially homogeneous IDP infall, impact ejecta are distributed locally, concentrated around the impact site. We find organics from asteroids and comets to dominate over IDP-borne organics at distances up to 150~km from the crater center. Our results may be important for the interpretation of in situ detections of organics on Mars., Comment: 30 pages, 5 figures, accepted for publication in Icarus

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

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., Comment: 6 pages main text, 4 pages references, 3 pages of figures for main text. 16 pages of supplementary material including 4 figures