(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The observation that the orbits of long-period Kuiper Belt objects are anomalously clustered in physical space has recently prompted the Planet Nine hypothesis - the proposed existence of a distant and eccentric planetary member of our solar system. Within the framework of this model, a Neptune-like perturber sculpts the orbital distribution of distant Kuiper Belt objects through a complex interplay of resonant and secular effects, such that in addition to perihelion-circulating objects, the surviving orbits get organized into apsidally aligned and anti-aligned configurations with respect to Planet Nine's orbit. In this work, we investigate the role of Kuiper Belt initial conditions on the evolution of the outer solar system using numerical simulations. Intriguingly, we find that the final perihelion distance distribution depends strongly on the primordial state of the system, and demonstrate that a bimodal structure corresponding to the existence of both aligned and anti-aligned clusters is only reproduced if the initial perihelion distribution is assumed to extend well beyond $\sim 36$ AU. The bimodality in the final perihelion distance distribution is due to the existence of permanently stable objects, with the lower perihelion peak corresponding to the anti-aligned orbits and the higher perihelion peak corresponding to the aligned orbits. We identify the mechanisms which enable the persistent stability of these objects and locate the regions of phase space in which they reside. The obtained results contextualize the Planet Nine hypothesis within the broader narrative of solar system formation, and offer further insight into the observational search for Planet Nine.

The Saturnian irregular satellite, Phoebe, can be broadly described as a water-rich rock. This object, which presumably originated from the same primordial population shared by the dynamically excited Kuiper Belt Objects, has received high resolution spectral imaging during the Cassini flyby. We present a new analysis of the Visual Infrared Mapping Spectrometer observations of Phoebe, which critically, includes a geometry correction routine that enables pixel-by-pixel mapping of visible and infrared spectral cubes directly onto the Phoebe shape model, even when an image exhibits significant trailing errors. The result of our re-analysis is a successful match of 40 images, producing spectral maps covering the majority of Phoebe's surface, roughly a 3rd of which is imaged by high resolution observations (<22 km per pixel resolution). There is no spot on Phoebe's surface that is absent of water absorption. The regions richest in water are clearly associated with the Jason and South Pole impact basins. We find Phoebe exhibits only three spectral types, and a water-ice concentration that correlates with physical depth and visible albedo. The water-rich and water-poor regions exhibit significantly different crater size frequency distributions, and different large crater morphologies. We propose that Phoebe once had a water-poor surface whose water-ice concentration was enhanced by basin forming impacts which exposed richer subsurface layers. Finally, we demonstrate that the range of Phoebe's water-ice absorption spans the same range exhibited by dynamically excited Kuiper Belt Objects. The common water-ice absorption depths and primordial origins, and the association of Phoebe's water-rich regions with its impact basins, suggests the plausible idea that Kuiper Belt Objects also originated with water-poor surfaces that were enhanced through stochastic collisional modification.

Models of the Solar System's dynamical evolution predict the dispersal of primitive planetesimals from their formative regions amongst the gas-giant planets due to the early phases of planetary migration. Consequently, carbonaceous objects were scattered both into the outer asteroid belt and out to the Kuiper Belt. These models predict that the Kuiper Belt should contain a small fraction of objects with carbonaceous surfaces, though to date, all reported visible reflectance spectra of small Kuiper Belt Objects (KBOs) are linear and featureless. We report the unusual reflectance spectrum of a small KBO, (120216) 2004 EW95, exhibiting a large drop in its near-UV reflectance and a broad shallow optical absorption feature centered at ~700 nm which is detected at greater than 4-sigma significance. These features, confirmed through multiple epochs of spectral photometry and spectroscopy, have respectively been associated with ferric oxides and phyllosilicates. The spectrum bears striking resemblance to those of some C-type asteroids, suggesting that 2004 EW95 may share a common origin with those objects. 2004 EW95 orbits the Sun in a stable mean motion resonance with Neptune, at relatively high eccentricity and inclination, suggesting it may have been emplaced there by some past dynamical instability. These results appear consistent with the aforementioned model predictions and are the first to show a reliably confirmed detection of silicate material on a small KBO.

The Large Synoptic Survey Telescope (LSST) is uniquely equipped to search for Solar System bodies due to its unprecedented combination of depth and wide field coverage. Over a ten-year period starting in 2022, LSST will generate the largest catalog of Solar System objects to date. The main goal of the LSST Solar System Science Collaboration (SSSC) is to facilitate the efforts of the planetary community to study the planets and small body populations residing within our Solar System using LSST data. To prepare for future survey cadence decisions and ensure that interesting and novel Solar System science is achievable with LSST, the SSSC has identified and prioritized key Solar System research areas for investigation with LSST in this roadmap. The ranked science priorities highlighted in this living document will inform LSST survey cadence decisions and aid in identifying software tools and pipelines needed to be developed by the planetary community as added value products and resources before the planned start of LSST science operations.

While satellites of mid- to small-Kuiper belt objects tend to be similar in size and brightness to their primaries, the largest Kuiper belt objects preferentially have satellites with small fractional brightness. In the two cases where the sizes and albedos of the small faint satellites have been measured, these satellites are seen to be small icy fragments consistent with collisional formation. Here we examine Dysnomia and Vanth, the satellites of Eris and Orcus, respectively. Using the Atacama Large Millimeter Array, we obtain the first spatially resolved observations of these systems at thermal wavelengths. We find a diameter for Dysnomia of 700$\pm$115 km and for Vanth of 475+/-75 km, with albedos of 0.04_+0.02_-0.01 and 0.08+/-0.02 respectively. Both Dysnomia and Vanth are indistinguishable from typical Kuiper belt objects of their size. Potential implications for the formation of these types of satellites are discussed.

We present a daytime thermal image of Europa taken with the Atacama Large Millimeter Array. The imaged region includes the area northwest of Pwyll Crater, which is associated with a nighttime thermal excess seen by the Galileo Photopolarimeter Radiometer and with two potential plume detections. We develop a global thermal model of Europa and simulate both the daytime and nighttime thermal emission to determine if the nighttime thermal anomaly is caused by excess endogenic heat flow, as might be expected from a plume source region. We find that the nighttime and daytime brightness temperatures near Pwyll Crater cannot be matched by including excess heat flow at that location. Rather, we can successfully model both measurements by increasing the local thermal inertia of the surface.

The Hilda asteroids are primitive bodies in resonance with Jupiter whose origin and physical properties are not well understood. Current models posit that these asteroids formed in the outer Solar System and were scattered along with the Jupiter Trojans into their present-day positions during a chaotic episode of dynamical restructuring. In order to explore the surface composition of these enigmatic objects in comparison with an analogous study of Trojans (Emery et al. 2011), we present new near-infrared spectra (0.7-2.5 $\mu$m) of 25 Hilda asteroids. No discernible absorption features are apparent in the data. Synthesizing the bimodalities in optical color and infrared reflectivity reported in previous studies, we classify 26 of the 28 Hildas in our spectral sample into the so-called less-red and red sub-populations and find that the two sub-populations have distinct average spectral shapes. Combining our results with visible spectra, we find that Trojans and Hildas possess similar overall spectral shapes, suggesting that the two minor body populations share a common progenitor population. A more detailed examination reveals that while the red Trojans and Hildas have nearly identical spectra, less-red Hildas are systematically bluer in the visible and redder in the near-infrared than less-red Trojans, indicating a putative broad, shallow absorption feature between 0.5 and 1.0 $\mu$m. We argue that the less-red and red objects found in both Hildas and Trojans represent two distinct surface chemistries and attribute the small discrepancy between less-red Hildas and Trojans to the difference in surface temperatures between the two regions.

The hypothesis that a massive Planet Nine exists in the outer solar system on a distant eccentric orbit was inspired by observations showing that the objects with the most distant eccentric orbits in the Kuiper belt have orbits which are physically aligned, that is, they are clustered in longitude of perihelion and have similar orbital planes. Questions have remained, however, about the effects of observational bias on these observations, particularly on the longitudes of perihelion. Specifically, distant eccentric Kuiper belt objects tend to be faint and only observable near their perihelia, suggesting that the longitudes of perihelion of the known distant objects could be strongly biased by the limited number of locations in the sky where deep surveys have been carried out. We have developed a method to rigorously estimate the longitude of perihelion bias for Kuiper belt observations. We find that the probability that the 10 known Kuiper belt objects with semimajor axis beyond 230 AU are drawn from a population with uniform longitude of perihelion is 1.2%. Combined with the observation that the orbital poles of these object are also clustered, the overall probability of detecting these two independent clusterings in a randomly distributed sample is 0.025%. While observational bias is clearly present in these observations, it is unlikely to explain the observed alignment of the distant eccentric Kuiper belt objects.

The densities of mid-sized Kuiper belt objects are a key constraint into understanding the assembly of objects in the outer solar system. These objects are critical for understanding the currently unexplained transition from the smallest Kuiper belt objects with densities lower than that of water to the largest objects with significant rock content. Mapping this transition is made difficult by the uncertainties in the diameters of these objects, which maps into an even larger uncertainty in volume and thus density. The substantial collecting area of the Atacama Large Millimeter Array allows significantly more precise measurements of thermal emission from outer solar system objects and could potentially greatly improve the density measurements. Here we use new thermal observations of four objects with satellites to explore the improvements possible with millimeter data. We find that effects due to effective emissivity at millimeter wavelengths make it difficult to use the millimeter data directly to find diameters and thus volumes for these bodies. In addition, we find that when including the effects of model uncertainty, the true uncertainties on the sizes of outer solar system objects measured with radiometry are likely larger than those previously published. Substantial improvement in object sizes will likely require precise occultation measurements.

We conducted a two-night photometric survey of small Kuiper Belt objects (KBOs) near opposition using the wide-field Hyper Suprime-Cam instrument on the 8.2 m Subaru Telescope. The survey covered about 90 deg^2 of sky, with each field imaged in the g and i bands. We detected 356 KBOs, ranging in absolute magnitude from 6.5 to 10.4. Filtering for high-inclination objects within the hot KBO population, we show that the g-i color distribution is strongly bimodal, indicative of two color classes - the red and very red subpopulations. After categorizing objects into the two subpopulations by color, we present the first dedicated analysis of the magnitude distributions of the individual color subpopulations and demonstrate that the two distributions are roughly identical in shape throughout the entire size range covered by our survey. Comparing the color distribution of small hot KBOs with that of Centaurs, we find that they have similar bimodal shapes, thereby providing strong confirmation of previous explanations for the attested bimodality of Centaurs. We also show that the magnitude distributions of the two KBO color subpopulations and the two color subpopulations observed in the Jupiter Trojans are statistically indistinguishable. Finally, we discuss a hypothesis describing the origin of the KBO color bimodality based on our survey results.

Current models of Solar System evolution posit that the asteroid populations in resonance with Jupiter are comprised of objects scattered inward from the outer Solar System during a period of dynamical instability. In this paper, we present a new analysis of the absolute magnitude and optical color distribution of Hilda asteroids, which lie in the 3:2 mean motion resonance with Jupiter, with the goal of comparing the bulk properties with previously published results from an analogous study of Jupiter Trojans. We report an updated power law fit of the Hilda magnitude distribution through H=14. Using photometric data listed in the Sloan Moving Object Catalog, we confirm the previously-reported strong bimodality in visible spectral slope distribution, indicative of two sub-populations with differing surface compositions. When considering collisional families separately, we find that collisional fragments follow a unimodal color distribution with spectral slope values consistent with the bluer of the two sub-populations. The color distributions of Hildas and Trojans are comparable and consistent with a scenario in which the color bimodality in both populations developed prior to emplacement into their present-day locations. We propose that the shallower magnitude distribution of the Hildas is a result of an initially much larger Hilda population, which was subsequently depleted as smaller bodies were preferentially ejected from the narrow 3:2 resonance via collisions. Altogether, these observations provide a strong case supporting a common origin for Hildas and Trojans as predicted by current dynamical instability theories of Solar System evolution.

The trans-Neptunian region of the solar system exhibits an intricate dynamical structure, much of which can be explained by an instability-driven orbital history of the giant planets. However, the origins of a highly inclined, and in certain cases retrograde, population of trans-Neptunian objects remain elusive within the framework of this evolutionary picture. In this work, we show that the existence of a distant, Neptune-like planet that resides on an eccentric and mildly inclined orbit fully accounts for the anomalous component the trans-Neptunian orbital distribution. Adopting the same parameters for Planet Nine as those previously invoked to explain the clustering of distant Kuiper belt orbits in physical space, we carry out a series of numerical experiments which elucidate the physical process though which highly inclined Kuiper belt objects with semi-major axes smaller than 100 AU are generated. The identified dynamical pathway demonstrates that enigmatic members of the Kuiper belt such as Drac and Niku are derived from the extended scattered disk of the solar system.

Hi'iaka is the larger outer satellite of the dwarf planet Haumea. Using relative photometry from the Hubble Space Telescope and Magellan and a phase dispersion minimization analysis, we have identified the rotation period of Hi'iaka to be ~9.8 hrs (double-peaked). This is ~120 times faster than its orbital period, creating new questions about the formation of this system and possible tidal evolution. The rapid rotation suggests that Hi'iaka could have a significant obliquity and spin precession that could be visible in light curves within a few years. We then turn to an investigation of what we learn about the (presently unclear) formation of the Haumea system and family based on this unexpectedly rapid rotation rate. We explore the importance of the initial semi-major axis and rotation period in tidal evolution theory and find they strongly influence the time required to despin to synchronous rotation, relevant to understanding a wide variety of satellite and binary systems. We find that despinning tides do not necessarily lead to synchronous spin periods for Hi'iaka, even if it formed near the Roche limit. Therefore the short rotation period of Hi'iaka does not rule out significant tidal evolution. Hi'iaka's spin period is also consistent with formation near its current location and spin up due to Haumea-centric impactors.

The six-degree obliquity of the sun suggests that either an asymmetry was present in the solar system's formation environment, or an external torque has misaligned the angular momentum vectors of the sun and the planets. However, the exact origin of this obliquity remains an open question. Batygin & Brown (2016) have recently shown that the physical alignment of distant Kuiper Belt orbits can be explained by a 5-20 Earth-mass planet on a distant, eccentric, and inclined orbit, with an approximate perihelion distance of ~250 AU. Using an analytic model for secular interactions between Planet Nine and the remaining giant planets, here we show that a planet with similar parameters can naturally generate the observed obliquity as well as the specific pole position of the sun's spin axis, from a nearly aligned initial state. Thus, Planet Nine offers a testable explanation for the otherwise mysterious spin-orbit misalignment of the solar system.

One of the most enigmatic and hitherto unexplained properties of Jupiter Trojans is their bimodal color distribution. This bimodality is indicative of two sub-populations within the Trojans, which have distinct size distributions. In this paper, we present a simple, plausible hypothesis for the origin and evolution of the two Trojan color sub-populations. In the framework of dynamical instability models of early Solar System evolution, which suggest a common primordial progenitor population for both Trojans and Kuiper belt objects, we use observational constraints to assert that the color bimodalities evident in both minor body populations developed within the primordial population prior to the onset of instability. We show that, beginning with an initial composition of rock and ices, location-dependent volatile loss through sublimation in this primordial population could have led to sharp changes in the surface composition with heliocentric distance. We propose that the depletion or retention of H$_{2}$S ice on the surface of these objects was the key factor in creating an initial color bimodality. Objects that retained H$_{2}$S on their surfaces developed characteristically redder colors upon irradiation than those that did not. After the bodies from the primordial population were scattered and emplaced into their current positions, they preserved this primordial color bimodality to the present day. We explore predictions of the volatile loss model - in particular, the effect of collisions within the Trojan population on the size distributions of the two sub-populations - and propose further experimental and observational tests of our hypothesis

To date, reflectance spectra of Jupiter Trojan asteroids have revealed no distinctive absorption features. For this reason, the surface composition of these objects remains a subject of speculation. Spectra have revealed, however, that the Jupiter Trojan asteroids consist of two distinct sub-populations which differ in the optical to near-infrared colors. The origins and compositional differences between the two sub-populations remain unclear. Here we report the results from a 2.2-3.8 $\mu$m spectral survey of a collection of 16 Jupiter Trojan asteroids, divided equally between the two sub-populations. We find clear spectral absorption features centered around 3.1 $\mu$m in the less red population. Additional absorption consistent with expected from organic materials might also be present. No such features are see in the red population. A strong correlation exists between the strength of the 3.1 $\mu$m absorption feature and the optical to near-infrared color of the objects. While traditionally absorptions such as these in dark asteroids are modeled as being due to fine-grain water frost, we find it physically implausible that the special circumstances required to create such fine-grained frost would exist on a substantial fraction of the Jupiter Trojan asteroids. We suggest, instead, that the 3.1 $\mu$m absorption on Trojans and other dark asteroids could be due to N-H stretch features. Additionally, we point out that reflectivities derived from WISE observations show a strong absorption beyond 4$\mu$m for both populations. The continuum of 3.1 $\mu$m features and the common absorption beyond 4 $\mu$m might suggest that both sub-populations of Jupiter Trojan asteroids formed in the same general region of the early solar system.

We use an extensive suite of numerical simulations to constrain the mass and orbit of Planet Nine, the recently proposed perturber in a distant eccentric orbit in the outer solar system. We compare our simulations to the observed population of aligned eccentric high semimajor axis Kuiper belt objects and determine which simulation parameters are statistically compatible with the observations. We find that only a narrow range of orbital elements can reproduce the observations. In particular, the combination of semimajor axis, eccentricity, and mass of Planet Nine strongly dictates the semimajor axis range of the orbital confinement of the distant eccentric Kuiper belt objects. Allowed orbits, which confine Kuiper belt objects with semimajor axis beyond 380~AU, have perihelia roughly between 150 and 350~AU, semimajor axes between 380 and 980~AU, and masses between 5 and 20 Earth masses. Orbitally confined objects also generally have orbital planes similar to that of the planet, suggesting that the planet is inclined approximately 30~degrees to the ecliptic. We compare the allowed orbital positions and estimated brightness of Planet Nine to previous and ongoing surveys which would be sensitive to the planet's detection and use these surveys to rule out approximately two-thirds of the planet's orbit. Planet Nine is likely near aphelion with an approximate brightness of $22<V<25$. At opposition, its motion, mainly due to parallax, can easily be detected within 24 hours.

Haumea is a dwarf planet with two known satellites, an unusually high spin rate, and a large collisional family, making it one of the most interesting objects in the outer solar system. A fully self-consistent formation scenario responsible for the satellite and family formation is still elusive, but some processes predict the initial formation of many small moons, similar to the small moons recently discovered around Pluto. Deep searches for regular satellites around KBOs are difficult due to observational limitations, but Haumea is one of the few for which sufficient data exist. We analyze Hubble Space Telescope (HST) observations, focusing on a ten-consecutive-orbit sequence obtained in July 2010, to search for new very small satellites. To maximize the search depth, we implement and validate a non-linear shift-and-stack method. No additional satellites of Haumea are found, but by implanting and recovering artificial sources, we characterize our sensitivity. At distances between $\sim$10,000 km and $\sim$350,000 km from Haumea, satellites with radii as small as $\sim$10 km are ruled out, assuming an albedo ($p \simeq 0.7$) similar to Haumea. We also rule out satellites larger than $\gtrsim$40 km in most of the Hill sphere using other HST data. This search method rules out objects similar in size to the small moons of Pluto. By developing clear criteria for determining the number of non-linear rates to use, we find that far fewer shift rates are required ($\sim$35) than might be expected. The non-linear shift-and-stack method to discover satellites (and other moving transients) is tractable, particularly in the regime where non-linear motion begins to manifest itself.

Recent analyses have shown that distant orbits within the scattered disk population of the Kuiper belt exhibit an unexpected clustering in their respective arguments of perihelion. While several hypotheses have been put forward to explain this alignment, to date, a theoretical model that can successfully account for the observations remains elusive. In this work we show that the orbits of distant Kuiper belt objects cluster not only in argument of perihelion, but also in physical space. We demonstrate that the perihelion positions and orbital planes of the objects are tightly confined and that such a clustering has only a probability of 0.007% to be due to chance, thus requiring a dynamical origin. We find that the observed orbital alignment can be maintained by a distant eccentric planet with mass greater than ~10 Earth masses, whose orbit lies in approximately the same plane as those of the distant Kuiper belt objects, but whose perihelion is 180 degrees away from the perihelia of the minor bodies. In addition to accounting for the observed orbital alignment, the existence of such a planet naturally explains the presence of high perihelion Sedna-like objects, as well as the known collection of high semimajor axis objects with inclinations between 60 and 150 degrees whose origin was previously unclear. Continued analysis of both distant and highly inclined outer solar system objects provides the opportunity for testing our hypothesis as well as further constraining the orbital elements and mass of the distant planet.

We present a comprehensive analysis of spatially resolved moderate spectral resolution near infrared spectra obtained with the adaptive optics system at the Keck Observatory. We identify three compositionally distinct end member regions: the trailing hemisphere bullseye, the leading hemisphere upper latitudes, and a third component associated with leading hemisphere chaos units. We interpret the composition of the three end member regions to be dominated by irradiation products, water ice, and evaporite deposits or salt brines, respectively. The third component is associated with geological features and distinct from the geography of irradiation, suggesting an endogenous identity. Identifying the endogenous composition is of particular interest for revealing the subsurface composition. However, its spectrum is not consistent with linear mixtures of the salt minerals previously considered relevant to Europa. The spectrum of this component is distinguished by distorted hydration features rather than distinct spectral features, indicating hydrated minerals but making unique identification difficult. In particular, it lacks features common to hydrated sulfate minerals, challenging the traditional view of an endogenous salty component dominated by Mg-sulfates. Chloride evaporite deposits are one possible alternative.

We present an analysis of survey observations targeting the leading L4 Jupiter Trojan cloud near opposition using the wide-field Suprime-Cam CCD camera on the 8.2 m Subaru Telescope. The survey covered about 38 deg$^2$ of sky and imaged 147 fields spread across a wide region of the L4 cloud. Each field was imaged in both the $g'$ and the $i'$ band, allowing for the measurement of $g-i$ color. We detected 557 Trojans in the observed fields, ranging in absolute magnitude from $H=10.0$ to $H = 20.3$. We fit the total magnitude distribution to a broken power law and show that the power-law slope rolls over from $0.45\pm 0.05$ to $0.36^{+0.05}_{-0.09}$ at a break magnitude of $H_{b}=14.93^{+0.73}_{-0.88}$. Combining the best-fit magnitude distribution of faint objects from our survey with an analysis of the magnitude distribution of bright objects listed in the Minor Planet Center catalog, we obtain the absolute magnitude distribution of Trojans over the entire range from $H=7.2$ to $H=16.4$. We show that the $g-i$ color of Trojans decreases with increasing magnitude. In the context of the less-red and red color populations, as classified in Wong et al. 2014 using photometric and spectroscopic data, we demonstrate that the observed trend in color for the faint Trojans is consistent with the expected trend derived from extrapolation of the best-fit color population magnitude distributions for bright catalogued Trojans. This indicates a steady increase in the relative number of less-red objects with decreasing size. Finally, we interpret our results using collisional modeling and propose several hypotheses for the color evolution of the Jupiter Trojan population.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of GSC 6214-210 A and B, a solar-mass member of the 5-10 Myr Upper Scorpius association with a 15 $\pm$ 2 Mjup companion orbiting at $\approx$330 AU (2.2"). Previous photometry and spectroscopy spanning 0.3-5 $\mu$m revealed optical and thermal excess as well as strong H$\alpha$ and Pa~$\beta$ emission originating from a circum-substellar accretion disk around GSC 6214-210 B, making it the lowest mass companion with unambiguous evidence of a subdisk. Despite ALMA's unprecedented sensitivity and angular resolution, neither component was detected in our 880 $\mu$m (341 GHz) continuum observations down to a 3-$\sigma$ limit of 0.22 mJy/beam. The corresponding constraints on the dust mass and total mass are <0.15 Mearth and <0.05 Mjup, respectively, or <0.003% and <0.3% of the mass of GSC 6214-210 B itself assuming a 100:1 gas-to-dust ratio and characteristic dust temperature of 10-20 K. If the host star possesses a putative circum-stellar disk then at most it is a meager 0.0015% of the primary mass, implying that giant planet formation has certainly ceased in this system. Considering these limits and its current accretion rate, GSC 6214-210 B appears to be at the end stages of assembly and is not expected to gain any appreciable mass over the next few Myr.

Here we present additional photometry of targets observed as part of the Hubble Wide Field Camera 3 Test of Surfaces in the Outer Solar System. 12 targets were re-observed with the Wide Field Camera 3 in optical and NIR wavebands designed to compliment those used during the first visit. Additionally, all observations originally presented by Fraser and Brown (2012) were reanalyzed through the same updated photometry pipeline. A reanalysis of the optical and NIR colour distribution reveals a bifurcated optical colour distribution and only two identifiable spectral classes, each of which occupies a broad range of colours and have correlated optical and NIR colours, in agreement with our previous findings. We report the detection of significant spectral variations on 5 targets which cannot be attributed to photometry errors, cosmic rays, point spread function or sensitivity variations, or other image artifacts capable of explaining the magnitude of the variation. The spectrally variable objects are found to have a broad range of dynamical classes, absolute magnitudes, exhibit a broad range of apparent magnitude variations, and are found in both compositional classes. The spectrally variable objects with sufficiently accurate colours for spectral classification maintain their membership, belonging to the same class at both epochs. 2005 TV189 exhibits a sufficiently broad difference in colour at the two epochs that span the full range of colours of the neutral class. This strongly argues that the neutral class is one single class with a broad range of colours, rather than the combination of multiple overlapping classes.

The Jupiter Trojans are a significant population of minor bodies in the middle Solar System that have garnered substantial interest in recent years. Several spectroscopic studies of these objects have revealed notable bimodalities with respect to near-infrared spectra, infrared albedo, and color, which suggest the existence of two distinct groups among the Trojan population. In this paper, we analyze the magnitude distributions of these two groups, which we refer to as the red and less red color populations. By compiling spectral and photometric data from several previous works, we show that the observed bimodalities are self-consistent and categorize 221 of the 842 Trojans with absolute magnitudes in the range H<12.3 into the two color populations. We demonstrate that the magnitude distributions of the two color populations are distinct to a high confidence level (>95%) and fit them individually to a broken power law, with special attention given to evaluating and correcting for incompleteness in the Trojan catalog as well as incompleteness in our categorization of objects. A comparison of the best-fit curves shows that the faint-end power-law slopes are markedly different for the two color populations, which indicates that the red and less red Trojans likely formed in different locations. We propose a few hypotheses for the origin and evolution of the Trojan population based on the analyzed data.