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

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

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

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

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

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7].

Pluto's atmospheric haze settles out rapidly compared with geological timescales. It needs to be accounted for as a surface material, distinct from Pluto's icy bedrock and from the volatile ices that migrate via sublimation and condensation on seasonal timescales. This paper explores how a steady supply of atmospheric haze might affect three distinct provinces on Pluto. We pose the question of why they each look so different from one another if the same haze material is settling out onto all of them. Cthulhu is a more ancient region with comparatively little present-day geological activity, where the haze appears to simply accumulate over time. Sputnik Planitia is a very active region where glacial convection, as well as sublimation and condensation rapidly refresh the surface, hiding recently deposited haze from view. Lowell Regio is a region of intermediate age featuring very distinct coloration from the rest of Pluto. Using a simple model haze particle as a colorant, we are not able to match the colors in both Lowell Regio and Cthulhu. To account for their distinct colors, we propose that after arrival at Pluto's surface, haze particles may be less inert than might be supposed from the low surface temperatures. They must either interact with local materials and environments to produce distinct products in different regions, or else the supply of haze must be non-uniform in time and/or location, such that different products are delivered to different places.

We present new Hubble Space Telescope and ground-based Keck observations and new Keplerian orbit solutions for the mutual orbit of binary Jupiter Trojan asteroid (617) Patroclus and Menoetius, targets of NASA's Lucy mission. We predict event times for the upcoming mutual event season, which is anticipated to run from late 2017 through mid 2019.

A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location of this material on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons. The escape of Pluto's atmosphere provides a potential feed stock for production of complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon's winter pole was proposed as an explanation on the basis of an image of Charon's northern hemisphere, but not modelled quantitatively. Here we report images of the southern hemisphere illuminated by Pluto-shine and also images taken during the approach phase showing the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon, the supply and temporary cold-trapping of material escaping from Pluto, and, while cold-trapped, its photolytic processing into more complex and less volatile molecules. The model results are consistent with the proposed mechanism producing the observed colour pattern on Charon.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

We report HST lightcurve observations of the New Horizons (NH) spacecraft encounter KBO (486958) 2014 MU69 acquired near opposition in July 2017. In order to plan the optimum flyby sequence the NH mission planners needed to learn as much as possible about the target in advance of the encounter. Specifically, from lightcurve data, encounter timing could be adjusted to accommodate a highly elongated, binary, or rapidly rotating target. HST astrometric (Porter et al. 2018) and stellar occultation (Buie et al. 2018) observations constrained MU69's orbit and diameter (21-41 km for an albedo of 0.15-0.04), respectively. Photometry from the astrometric dataset suggested a variability of $\ge$0.3 mags, but they did not determine the period or provide shape information. We strategically spaced 24 HST orbits over 9 days to investigate rotation periods from approximately 2-100 hours and to better constrain the lightcurve amplitude. Until NH detected MU69 in its optical navigation images beginning in August 2018, this HST campaign provided the most accurate photometry to date. The mean variation in our data is 0.15 mags which suggests that MU69 is either nearly spherical (a:b axis ratio of 1:1.15), or its pole vector is pointed near the line of sight to Earth; this interpretation does not preclude a near-contact binary or bi-lobed object. However, image stacks do conclude that MU69 does not have a binary companion $\ge$2000 km with a sensitivity to 29th magnitude. We report with confidence that MU69 is not both rapidly rotating and highly elongated. We note that our results are consistent with the fly-by imagery and orientation of MU69 (Stern et al. 2019). The combined dataset also suggests that within the KBO lightcurve literature there are likely other objects which share a geometric configuration like MU69 resulting in an underestimate of the contact binary fraction for the CC Kuiper Belt.

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

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.

The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was ejected. The closest encounter, at 0.60pc (0.53-0.67pc, 90% confidence interval), was with the M2.5 dwarf HIP 3757 at a relative velocity of 24.7km/s, 1Myr ago. A more distant encounter (1.6pc) but with a lower encounter (ejection) velocity of 10.7km/s was with the G5 dwarf HD 292249, 3.8Myr ago. Two more stars have encounter distances and velocities intermediate to these. The encounter parameters are similar across six different non-gravitational trajectories for 'Oumuamua. Ejection of 'Oumuamua by scattering from a giant planet in one of the systems is plausible, but requires a rather unlikely configuration to achieve the high velocities found. A binary star system is more likely to produce the observed velocities. None of the four home candidates have published exoplanets or are known to be binaries. Given that the 7 million stars in Gaia DR2 with 6D phase space information is just a small fraction of all stars for which we can eventually reconstruct orbits, it is a priori unlikely that our current search would find 'Oumuamua's home star system. As 'Oumuamua is expected to pass within 1pc of about 20 stars and brown dwarfs every Myr, the plausibility of a home system depends also on an appropriate (low) encounter velocity.

The Trojan clouds contain a surprisingly diverse population of objects. It is only by understanding this diversity can we unravel their history. The more objects Lucy visits the better we can constrain the formation and evolution of the giant planets. While Lucy will visit four Trojans in the L$_4$ cloud, its trajectory only passes one known Trojan system, the PM binary, in L$_5$. There is ample time in the Lucy timeline to possibly add additional L$_5$ targets, thereby maximizing the mission's science return. In the next decade, the 8.4-m Large Synoptic Survey Telescope (LSST), is expected to detect nearly 300,000 Jupiter Trojans between approximately 16 and 24.5 r-band magnitudes in the $\sim$18,000 deg$^2$ main survey and the expected $\sim$4,000 deg$^2$ northern ecliptic extension. We explore potential opportunities for LSST to find new Lucy L$_5$ fly-by targets.

The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

We present a community-led assessment of the solar system investigations achievable with NASA's next-generation space telescope, the Wide Field InfraRed Survey Telescope (WFIRST). WFIRST will provide imaging, spectroscopic, and coronagraphic capabilities from 0.43-2.0 $\mu$m and will be a potential contemporary and eventual successor to JWST. Surveys of irregular satellites and minor bodies are where WFIRST will excel with its 0.28 deg$^2$ field of view Wide Field Instrument (WFI). Potential ground-breaking discoveries from WFIRST could include detection of the first minor bodies orbiting in the Inner Oort Cloud, identification of additional Earth Trojan asteroids, and the discovery and characterization of asteroid binary systems similar to Ida/Dactyl. Additional investigations into asteroids, giant planet satellites, Trojan asteroids, Centaurs, Kuiper Belt Objects, and comets are presented. Previous use of astrophysics assets for solar system science and synergies between WFIRST, LSST, JWST, and the proposed NEOCam mission are discussed. We also present the case for implementation of moving target tracking, a feature that will benefit from the heritage of JWST and enable a broader range of solar system observations.

We observed the July 17, 2017 stellar occultation of HD 168233 by the Kuiper Belt Object (486958) 2014 MU$_{69}$, the close flyby target of the extended New Horizons mission. Rather than capture a solid body occultation by the KBO itself, our program aimed to constrain the opacity of rings, moons, or other debris in the nearby environment. We used the Hubble Space Telescope Fine Guidance Sensors (HST FGS) instrument in TRANS F583W mode to collect 40 Hz time resolution photometry of the stellar occultation star for two HST orbits during this observation. We present the results of reduction and calibration of the HST FGS photometry, and set upper limits on rings or other dust opacity within the Hill sphere of (486958) 2014 MU$_{69}$ at distances ranging from $\sim$20,000 km to $\sim$75,000 km from the main body.

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

We report results from analyzing the B612 asteroid observation data taken by the sCMOS cameras on board of Planet SkySat-3 using the synthetic tracking technique. The analysis demonstrates the expected sensitivity improvement in the signal-to-noise ratio of the asteroids from properly stacking up the the short exposure images in post-processing.

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum.

We searched for dust or debris rings in the Pluto-Charon system before, during, and after the New Horizons encounter. Methodologies included searching for back-scattered light during the approach to Pluto (phase $\sim15^\circ$), in situ detection of impacting particles, a search for stellar occultations near the time of closest approach, and by forward-scattered light during departure (phase $\sim165^\circ$). A search using HST prior to the encounter also contributed to the results. No rings, debris, or dust features were observed, but our detection limits provide an improved picture of the environment throughout the Pluto-Charon system. Searches for rings in back-scattered light covered 35,000-250,000 km from the system barycenter, a zone that starts interior to the orbit of Styx, and extends to four times the orbital radius of Hydra. We obtained our firmest limits using the NH LORRI camera in the inner half of this region. Our limits on the normal $I/F$ of an unseen ring depends on the radial scale of the rings: $2\times10^{-8}$ ($3\sigma$) for 1500 km wide rings, $1\times10^{-8}$ for 6000 km rings, and $7\times10^{-9}$ for 12,000 km rings. Beyond $\sim100,000$ km from Pluto, HST observations limit normal $I/F$ to $\sim8\times10^{-8}$. Searches for dust from forward-scattered light extended from the surface of Pluto to the Pluto-Charon Hill sphere ($r_{\rm Hill}=6.4\times10^6$ km). No evidence for rings or dust was detected to normal $I/F$ limits of $\sim8.9\times10^{-7}$ on $\sim10^4$ km scales. Four occulation observations also probed the space interior to Hydra, but again no dust or debris was detected. Elsewhere in the solar system, small moons commonly share their orbits with faint dust rings. Our results suggest that small grains are quickly lost from the system due to solar radiation pressure, whereas larger particles are unstable due to perturbations by the known moons.

In this work we investigated the longitudinal (zonal) variability of H$_2$O and ammonia (NH$_3$) hydrate ices on the surface of Charon through analysis of the 1.65 $\mu$m and 2.21 $\mu$m absorption features, respectively. Near-infrared spectra presented here were obtained between 2015-07-14 and 2015-08-30 UT with the OSIRIS integral field spectrograph on Keck I. Spectra centered on six different sub-observer longitudes were obtained through the Hbb (1.473-1.803 $\mu$m) and Kbb (1.965-2.381 $\mu$m) filters. Gaussian functions were fit to the aforementioned bands to obtain information on band center, band depth, full width at half maximum, and band area. The shift in the band center of the temperature-dependent 1.65 $\mu$m feature was used to calculate the H$_2$O ice temperature. The mean temperature of the ice on the observable portion of Charon's surface is 45$\pm$14 K and we report no statistically significant variations in temperature across the surface. We hypothesize that the crystalline and amorphous phases of water ice reached equilibrium over 3.5 Gyr ago, with thermal recrystallization balancing the effects of irradiation amorphization. We do not believe that cryovolcanism is necessary to explain the presence of crystalline water ice on the surface of Charon. Absorption from ammonia species is detected between 12$^{\circ}$ and 290$^{\circ}$, in agreement with results from New Horizons. Ongoing diffusion of ammonia through the rocky mantle and upper layer of water ice is one possible mechanism for maintaining its presence in Charon's surface ice. Reduced Charon spectra corrected for telluric and solar absorption are available as supplementary online material.

NASA's New Horizons spacecraft observed (15810) 1994 JR$_1$, a 3:2 resonant Kupier Belt Object (KBO), using the LOng Range Reconnaissance Imager (LORRI) on November 2, 2015 from a distance of 1.85 AU, and again on April 7, 2016 from a distance of 0.71 AU. These were the first close observations of any KBO other than Pluto. Combining ground-based and Hubble Space Telecope (HST) observations at small phase angles and the LORRI observations at higher phase angles, we produced the first disk-integrated solar phase curve of a typical KBO from $\alpha$=0.6-58$^\circ$. Observations at these geometries, attainable only from a spacecraft in the outer Solar System, constrain surface properties such as macroscopic roughness and the single particle phase function. 1994 JR$_1$ has a rough surface with a 37$\pm$5$^\circ$ mean topographic slope angle and has a relatively rapid rotation period of 5.47$\pm$0.33 hours. 1994 JR$_1$ is currently 2.7 AU from Pluto; our astrometric points enable high-precision orbit determination and integrations which show that it comes this close to Pluto every 2.4 million years (10$^4$ heliocentric orbits), causing Pluto to perturb 1994 JR$_1$. During the November spacecraft observation, the KBO was simultaneously observed using HST in two colors, confirming its very red spectral slope. These observations have laid the groundwork for numerous potential future distant KBO observations in the New Horizons-Kuiper Belt Extended Mission.

We present results derived from the first multi-chord stellar occultation by the trans-Neptunian object (229762) 2007 UK$_{126}$, observed on 2014 November 15. The event was observed by the Research and Education Collaborative Occultation Network (RECON) project and International Occultation Timing Association (IOTA) collaborators throughout the United States. Use of two different data analysis methods obtain a satisfactory fit to seven chords, yelding an elliptical fit to the chords with an equatorial radius of $R=338_{-10} ^{+15}$ km and equivalent radius of $R_{eq}=319_{-7} ^{+14}$ km. A circular fit also gives a radius of $R=324_{-23} ^{+30}$ km. Assuming that the object is a Maclaurin spheroid with indeterminate aspect angle, and using two published absolute magnitudes for the body, we derive possible ranges for geometric albedo between $p_{V}=0.159_{-0.013} ^{+0.007}$ and $p_{R}=0.189_{-0.015}^{+0.009}$, and for the body oblateness between $\epsilon=0.105_{-0.040} ^{+0.050}$ and $\epsilon=0.118_{-0.048} ^{+0.055}$. For a nominal rotational period of 11.05 h, an upper limit for density of $\rho=1740$ kg~m$^{-3}$ is estimated for the body.

We present an analysis of surveying the inner Solar System for objects that may pose some threat to the Earth. Most of the analysis is based on understanding the capability provided by Sentinel, a concept for an infrared space-based telescope placed in a heliocentric orbit near the distance of Venus. From this analysis, we show 1) the size range being targeted can affect the survey design, 2) the orbit distribution of the target sample can affect the survey design, 3) minimum observational arc length during the survey is an important metric of survey performance, and 4) surveys must consider objects as small as D=15-30 m to meet the goal of identifying objects that have the potential to cause damage on Earth in the next 100 years. Sentinel will be able to find 50% of all impactors larger than 40 meters in a 6.5 year survey. The Sentinel mission concept is shown to be as effective as any survey in finding objects bigger than D=140 m but is more effective when applied to finding smaller objects on Earth-impacting orbits. Sentinel is also more effective at finding objects of interest for human exploration that benefit from lower propulsion requirements. To explore the interaction between space and ground search programs, we also study a case where Sentinel is combined with the Large Synoptic Survey Telescope and show the benefit of placing a space-based observatory in an orbit that reduces the overlap in search regions with a ground-based telescope. In this case, Sentinel+LSST can find more than 70% of the impactors larger than 40 meters assuming a 6.5 year lifetime for Sentinel and 10 years for LSST.

Approach images taken by the LORRI imaging system during the New Horizons spacecraft encounter have been used to determine the mean radii and shapes of Pluto and Charon. The primary observations are limb locations derived using three independent approaches. The resulting mean radii of Pluto and Charon are 1188.3 +/- 1.6 km and 606.0 +/- 1.0 km, respectively (2-sigma). The corresponding densities are 1854 +/- 11 kg m-3 and 1701 +/- 33 kg m-3 (2-sigma). The Charon radius value is consistent with previous Earth-based occultation estimates. The Pluto radius estimate is consistent with solar occultation measurements performed by the ALICE and Fine Sun Sensor instruments on New Horizons. Neither Pluto nor Charon show any evidence for tidal/rotational distortions; upper bounds on the oblateness are <0.6% and <0.5%, respectively.

We describe the discovery of a satellite in orbit about the dwarf planet (136472) Makemake. This satellite, provisionally designated S/2015 (136472) 1, was detected in imaging data collected with the Hubble Space Telescope's Wide Field Camera 3 on UTC April 27, 2015 at 7.80$\pm$0.04 magnitudes fainter than Makemake. It likely evaded detection in previous satellite searches due to a nearly edge-on orbital configuration, placing it deep within the glare of Makemake during a substantial fraction of its orbital period. This configuration would place Makemake and its satellite near a mutual event season. Insufficient orbital motion was detected to make a detailed characterization of its orbital properties, prohibiting a measurement of the system mass with the discovery data alone. Preliminary analysis indicates that if the orbit is circular, its orbital period must be longer than 12.4 days, and must have a semi-major axis $\gtrsim$21,000 km. We find that the properties of Makemake's moon suggest that the majority of the dark material detected in the system by thermal observations may not reside on the surface of Makemake, but may instead be attributable to S/2015 (136472) 1 having a uniform dark surface. This "dark moon hypothesis" can be directly tested with future JWST observations. We discuss the implications of this discovery for the spin state, figure, and thermal properties of Makemake and the apparent ubiquity of trans-Neptunian dwarf planet satellites.

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

The New Horizons mission has provided resolved measurements of Pluto's moons Styx, Nix, Kerberos, and Hydra. All four are small, with equivalent spherical diameters of $\approx$40 km for Nix and Hydra and ~10 km for Styx and Kerberos. They are also highly elongated, with maximum to minimum axis ratios of $\approx$2. All four moons have high albedos ( $\approx$50-90 %) suggestive of a water-ice surface composition. Crater densities on Nix and Hydra imply surface ages $\gtrsim$ 4 Ga. The small moons rotate much faster than synchronous, with rotational poles clustered nearly orthogonal to the common pole directions of Pluto and Charon. These results reinforce the hypothesis that the small moons formed in the aftermath of a collision that produced the Pluto-Charon binary.

The New Horizons spacecraft mapped colors and infrared spectra across the encounter hemispheres of Pluto and Charon. The volatile ices CH$_4$, CO, and N$_2$, that dominate Pluto's surface, have complicated spatial distributions resulting from sublimation, condensation, and glacial flow acting over seasonal and geological timescales. Pluto's H$_2$O ice "bedrock" is also mapped, with isolated outcrops occurring in a variety of settings. Pluto's surface exhibits complex regional color diversity associated with its distinct provinces. Charon's color pattern is simpler, dominated by neutral low latitudes and a reddish northern polar region. Charon near infrared spectra reveal highly localized areas with strong NH$_3$ absorption tied to small craters with relatively fresh-appearing impact ejecta.

We discuss two semi-independent calibration techniques used to determine the in-flight radiometric calibration for the New Horizons' Multi-spectral Visible Imaging Camera (MVIC). The first calibration technique compares the observed stellar flux to modeled values. The difference between the two provides a calibration factor that allows the observed flux to be adjusted to the expected levels for all observations, for each detector. The second calibration technique is a channel-wise relative radiometric calibration for MVIC's blue, near-infrared and methane color channels using observations of Charon and scaling from the red channel stellar calibration. Both calibration techniques produce very similar results (better than 7% agreement), providing strong validation for the techniques used. Since the stellar calibration can be performed without a color target in the field of view and covers all of MVIC's detectors, this calibration was used to provide the radiometric keywords delivered by the New Horizons project to the Planetary Data System (PDS). These keywords allow each observation to be converted from counts to physical units; a description of how these keywords were generated is included. Finally, mitigation techniques adopted for the gain drift observed in the near-infrared detector and one of the panchromatic framing cameras is also discussed.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

(3200) Phaethon exhibits both comet- and asteroid-like properties, suggesting it could be a rare transitional object such as a dormant comet or previously volatile-rich asteroid. This justifies detailed study of (3200) Phaethon's physical properties, as a better understanding of asteroid-comet transition objects can provide insight into minor body evolution. We therefore acquired time-series photometry of (3200) Phaethon over 15 nights from 1994 to 2013, primarily using the Tektronix 2048x2048 pixel CCD on the University of Hawaii 2.2-m telescope. We utilized light curve inversion to: (1) refine (3200) Phaethon's rotational period to P=3.6032+/-0.0008 h; (2) estimate a rotational pole orientation of lambda=+85+/-13 degrees and beta=-20+/-10 degrees; and (3) derive a shape model. We also used our extensive light curve dataset to estimate the slope parameter of (3200) Phaethon's phase curve as G~0.06, consistent with C-type asteroids. We discuss how this highly oblique pole orientation with a negative ecliptic latitude supports previous evidence for (3200) Phaethon's origin in the inner main asteroid belt as well as the potential for deeply buried volatiles fueling impulsive yet rare cometary outbursts.

The Deep Ecliptic Survey (DES) discovered hundreds of Kuiper Belt objects from 1998-2005. Follow-up observations yielded 304 objects with good dynamical classifications (Classical, Scattered, Centaur, or 16 mean-motion resonances with Neptune). The DES search fields are well documented, enabling us to calculate the probability of detecting objects with particular orbital parameters and absolute magnitudes at a randomized point in each orbit. Grouping objects together by dynamical class leads, we estimate the orbital element distributions (a, e, i) for the largest three classes (Classical, 3:2, and Scattered) using maximum likelihood. Using H-magnitude as a proxy for the object size, we fit a power law to the number of objects for 8 classes with at least 5 detected members (246 objects). The best Classical slope is alpha=1.02+/-0.01 (observed from 5<=H<=7.2). Six dynamical classes (Scattered plus 5 resonances) are consistent in slope with the Classicals, though the absolute number of objects is scaled. The exception to the power law relation are the Centaurs (non-resonant with perihelia closer than Neptune, and thus detectable at smaller sizes), with alpha=0.42+/-0.02 (7.5<H<11). This is consistent with a knee in the H-distribution around H=7.2 as reported elsewhere (Bernstein et al. 2004, Fraser et al. 2014). Based on the Classical-derived magnitude distribution, the total number of objects (H<=7) in each class are: Classical (2100+/-300 objects), Scattered (2800+/-400), 3:2 (570+/-80), 2:1 (400+/-50), 5:2 (270+/-40), 7:4 (69+/-9), 5:3 (60+/-8). The independent estimate for the number of Centaurs in the same H range is 13+/-5. If instead all objects are divided by inclination into "Hot" and "Cold" populations, following Fraser et al. (2014), we find that alphaHot=0.90+/-0.02, while alphaCold=1.32+/-0.02, in good agreement with that work.

We report on Hubble Space Telescope (HST) observations of three Kuiper Belt Objects (KBOs), discovered in our dedicated ground-based search campaign, that are candidates for long-range observations from the New Horizons spacecraft: 2011 JY31, 2011 HZ102, and 2013 LU35. Astrometry with HST enables both current and future critical accuracy improvements for orbit precision, required for possible New Horizons observations, beyond what can be obtained from the ground. Photometric colors of all three objects are red, typical of the Cold Classical dynamical population within which they reside; they are also the faintest KBOs to have had their colors measured. None are observed to be binary with HST above separations of ~0.02 arcsec (~700 km at 44 AU) and {\Delta}m less than or equal to 0.5.

Combining stellar occultation observations probing Pluto's atmosphere from 1988 to 2013 and models of energy balance between Pluto's surface and atmosphere, we conclude that Pluto's atmosphere does not collapse at any point in its 248-year orbit. The occultation results show an increasing atmospheric pressure with time in the current epoch, a trend present only in models with a high thermal inertia and a permanent N2 ice cap at Pluto's north rotational pole.

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

We report results from monitoring Pluto's 0.8 to 2.4 {\mu}m reflectance spectrum with IRTF/SpeX on 65 nights over the dozen years from 2001 to 2012. The spectra show vibrational absorption features of simple molecules CH4, CO, and N2 condensed as ices on Pluto's surface. These absorptions are modulated by the planet's 6.39 day rotation period, enabling us to constrain the longitudinal distributions of the three ices. Absorptions of CO and N2 are concentrated on Pluto's anti-Charon hemisphere, unlike absorptions of less volatile CH4 ice that are offset by roughly 90{\deg} from the longitude of maximum CO and N2 absorption. In addition to the diurnal variations, the spectra show longer term trends. On decadal timescales, Pluto's stronger CH4 absorption bands have been getting deeper, while the amplitude of their diurnal variation is diminishing, consistent with additional CH4 absorption at high northern latitudes rotating into view as the sub-Earth latitude moves north (as defined by the system's angular momentum vector). Unlike the CH4 absorptions, Pluto's CO and N2 absorptions appear to be declining over time, suggesting more equatorial or southerly distributions of those species. Comparisons of geometrically-matched pairs of observations favor geometric explanations for the observed secular changes in CO and N2 absorption, although seasonal volatile transport could be at least partly responsible. The case for a volatile transport contribution to the secular evolution looks strongest for CH4 ice, despite it being the least volatile of the three ices.

Hubble Space Telescope observations between 2001 and 2010 resolved the binary components of the Cold Classical transneptunian object (79360) Sila-Nunam (provisionally designated 1997 CS29). From these observations we have determined the circular, retrograde mutual orbit of Nunam relative to Sila with a period of 12.50995 \pm 0.00036 days and a semimajor axis of 2777 \pm 19 km. A multi-year season of mutual events, in which the two near-equal brightness bodies alternate in passing in front of one another as seen from Earth, is in progress right now, and on 2011 Feb. 1 UT, one such event was observed from two different telescopes. The mutual event season offers a rich opportunity to learn much more about this barely-resolvable binary system, potentially including component sizes, colors, shapes, and albedo patterns. The low eccentricity of the orbit and a photometric lightcurve that appears to coincide with the orbital period are consistent with a system that is tidally locked and synchronized, like the Pluto-Charon system. The orbital period and semimajor axis imply a system mass of (10.84 \pm 0.22) \times 10^18 kg, which can be combined with a size estimate based on Spitzer and Herschel thermal infrared observations to infer an average bulk density of 0.72 +0.37 -0.23 g cm^-3, comparable to the very low bulk densities estimated for small transneptunian binaries of other dynamical classes.

Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own protosolar nebula using the properties of the small body reservoirs in the solar system. In standard planet formation theory, after Jupiter and Saturn formed they scattered a significant number of remnant planetesimals into highly eccentric orbits. In this paper, we show that if there had been a massive, extended protoplanetary disk at that time, then the disk would have excited Kozai oscillations in some of the scattered objects, driving them into high-inclination (i > 50 deg), low-eccentricity orbits (q > 30 AU). The dissipation of the gaseous disk would strand a subset of objects in these high-inclination orbits; orbits that are stable on Gyr time scales. To date, surveys have not detected any Kuiper Belt Objects with orbits consistent with this dynamical mechanism. Using these non-detections by the Deep Ecliptic Survey (DES) and the Palomar Distant Solar System Survey we are able to rule out an extended gaseous protoplanetary disk (R_D > 80 AU) in our solar system at the time of Jupiter's formation. Future deep all sky surveys such as the Large Synoptic Survey Telescope (LSST) will all us to further constrain the size of the protoplanetary disk.

We present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK183, (149780) Altjira, 2001 QY297, 2003 QW111, and 2003 QY90 based on Hubble Space Telescope and Keck 2 laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits yields a sample of 22 mutual orbits for which the period P, semimajor axis a, and eccentricity e have been determined. These orbits have mutual periods ranging from 5 to over 800 days, semimajor axes ranging from 1,600 to 37,000 km, eccentricities ranging from 0 to 0.8, and system masses ranging from 2 x 10^17 to 2 x 10^22 kg. Based on the relative brightnesses of primaries and secondaries, most of these systems consist of near equal-sized pairs, although a few of the most massive systems are more lopsided. The observed distribution of orbital properties suggests that the most loosely-bound transneptunian binary systems are only found on dynamically cold heliocentric orbits. Of the 22 known binary mutual or-bits, orientation ambiguities are now resolved for 9, of which 7 are prograde and 2 are retro-grade, consistent with a random distribution of orbital orientations, but not with models predicting a strong preference for retrograde orbits. To the extent that other perturbations are not dominant, the binary systems undergo Kozai oscillations of their eccentricities and inclinations with periods of the order of tens of thousands to millions of years, some with strikingly high amplitudes.

The PHOT (Portable High-Speed Occultation Telescope) systems were developed for the specific purpose of observing stellar occultations by solar system objects. Stellar occultations have unique observing constraints: they may only be observable from certain parts of the globe; they often require a rapid observing cadence; and they require accurate timestamp information for each exposure. The PHOT systems consist of 14" telescopes, CCD cameras, camera mounting plates, GPS-based time standards, and data acquisition computers. The PHOT systems are similar in principle to the POETS systems (Portable Occultation, Eclipse and Transit Systems, described by Souza et al. 2006 and reported on by Gulbis et al. 2008), with the main differences being (a) different CCD/Cameras with slightly different specifications and (b) a stand-alone custom-built time standard used by PHOT, whereas POETS uses a commercial time-standard that is controlled from a computer. Since 2005, PHOT systems have been deployed on over two dozen occasions to sites in the US, Mexico, Chile, Namibia, South Africa, France, Austria, Switzerland, Australia and New Zealand, mounted on portable 14" telescopes or on larger stationary telescopes. Occultation light curves acquired from the 3.9-m AAT (Anglo-Australian Telescope) have produced photometric signal-to-noise ratios (SNR) of 333 per scale height for a stellar occultation by Pluto (Young et al. 2008). In this paper we describe the seven PHOT subsystems in detail (telescopes, cameras, timers and data stations) and present SNR estimates for actual and predicted occultations as functions of star brightness, telescope aperture and frame rate.

Using data from the Deep Ecliptic Survey (DES), we investigate the inclination distributions of objects in the Kuiper Belt. We present a derivation for observational bias removal and use this procedure to generate unbiased inclination distributions for Kuiper Belt objects (KBOs) of different DES dynamical classes, with respect to the Kuiper Belt Plane. Consistent with previous results, we find that the inclination distribution for all DES KBOs is well fit by the sum of two Gaussians, or a Gaussian plus a generalized Lorentzian, multiplied by sin i. Approximately 80% of KBOs are in the high-inclination grouping. We find that Classical object inclinations are well fit by sin i multiplied by the sum of two Gaussians, with roughly even distribution between Gaussians of widths 2.0 -0.5/+0.6 degrees and 8.1 -2.1/+2.6 degrees. Objects in different resonances exhibit different inclination distributions. The inclinations of Scattered objects are best matched by sin i multiplied by a single Gaussian that is centered at 19.1 -3.6/+3.9 degrees with a width of 6.9 -2.7/+4.1 degrees. Centaur inclinations peak just below 20 degrees, with one exceptionally high-inclination object near 80 degrees. The currently observed inclination distribution of the Centaurs is not dissimilar to that of the Scattered Extended KBOs and Jupiter-family comets, but is significantly different from the Classical and Resonant KBOs. While the sample sizes of some dynamical classes are still small, these results should begin to serve as a critical diagnostic for models of Solar System evolution.

This report arises from an ongoing program to monitor Neptune's largest moon Triton spectroscopically in the 0.8 to 2.4 micron range using IRTF/SpeX. Our objective is to search for changes on Triton's surface as witnessed by changes in the infrared absorption bands of its surface ices N2, CH4, H2O, CO, and CO2. We have recorded infrared spectra of Triton on 53 nights over the ten apparitions from 2000 through 2009. The data generally confirm our previously reported diurnal spectral variations of the ice absorption bands (Grundy & Young 2004). Nitrogen ice shows a large amplitude variation, with much stronger absorption on Triton's Neptune-facing hemisphere. We present evidence for seasonal evolution of Triton's N2 ice: the 2.15 micron absorption band appears to be diminishing, especially on the Neptune-facing hemisphere. Although it is mostly dissolved in N2 ice, Triton's CH4 ice shows a very different longitudinal variation from the N2 ice, challenging assumptions of how the two ices behave. Unlike Triton's CH4 ice, the CO ice does exhibit longitudinal variation very similar to the N2 ice, implying that CO and N2 condense and sublimate together, maintaining a consistent mixing ratio. Absorptions by H2O and CO2 ices show negligible variation as Triton rotates, implying very uniform and/or high latitude spatial distributions for those two non-volatile ices.

We explore the variability and cross-frequency correlation of the flux density and polarization of the blazar OJ287, using imaging at 43 GHz with the Very Long Baseline Array, as well as optical and near-infrared polarimetry. The polarization and flux density in both the optical waveband and the 43 GHz compact core increased by a small amount in late 2005, and increased significantly along with the near-IR polarization and flux density over the course of 10 days in early 2006. Furthermore, the values of the electric vector position angle (EVPA) at the three wavebands are similar. At 43 GHz, the EVPA of the blazar core is perpendicular to the flow of the jet, while the EVPAs of emerging superluminal knots are aligned parallel to the jet axis. The core polarization is that expected if shear aligns the magnetic field at the boundary between flows of disparate velocities within the jet. Using variations in flux density, percentage polarization, and EVPA, we model the inner jet as a spine-sheath system. The model jet contains a turbulent spine of half-width 1.2 degrees and maximum Lorentz factor of 16.5, a turbulent sheath with Lorentz factor of 5, and a boundary region of sheared field between the spine and sheath. Transverse shocks propagating along the fast, turbulent spine can explain the superluminal knots. The observed flux density and polarization variations are then compatible with changes in the direction of the inner jet caused by a temporary change in the position of the core if the spine contains wiggles owing to an instability. In addition, we can explain a stable offset of optical and near-IR percentage polarization by a steepening of spectral index with frequency, as supported by the data.

We present Hubble Space Telescope observations of six binary transneptunian systems: 2000 QL251, 2003 TJ58, 2001 XR254, 1999 OJ4, (134860) 2000 OJ67, and 2004 PB108. The mutual orbits of these systems are found to have periods ranging from 22 to 137 days, semimajor axes ranging from 2360 to 10500 km, and eccentricities ranging from 0.09 to 0.55. These orbital parameters enable estimation of system masses ranging from 0.2 to 9.7 x 10+18 kg. For reasonable assumptions of bulk density (0.5 to 2.0 g cm-3), the masses can be combined with visible photometry to constrain sizes and albedos. The resulting albedos are consistent with an emerging picture of the dynamically "Cold" Classical subpopulation having relatively high albedos, compared with comparablysized objects on more dynamically excited orbits.

We present observations of thermal emission from fifteen transneptunian objects (TNOs) made using the Spitzer Space Telescope. Thirteen of the targets are members of the Classical population: six dynamically hot Classicals, five dynamically cold Classicals, and two dynamically cold inner Classical Kuiper Belt Objects (KBOs). We fit our observations using thermal models to determine the sizes and albedos of our targets finding that the cold Classical TNOs have distinctly higher visual albedos than the hot Classicals and other TNO dynamical classes. The cold Classicals are known to be distinct from other TNOs in terms of their color distribution, size distribution, and binarity fraction. The Classical objects in our sample all have red colors yet they show a diversity of albedos which suggests that there is not a simple relationship between albedo and color. As a consequence of high albedos, the mass estimate of the cold Classical Kuiper Belt is reduced from approximately 0.01 Earth masses to approximately 0.001 Earth masses. Our results also increase significantly the sample of small Classical KBOs with known albedos and sizes from 21 to 32 such objects.

We report resolved photometry of the primary and secondary components of 23 transneptunian binaries obtained with the Hubble Space Telescope. V-I colors of the components range from 0.7 to 1.5 with a median uncertainty of 0.06 magnitudes. The colors of the primaries and secondaries are correlated with a Spearman rank correlation probability of 99.99991%, 5 sigma for a normal distribution. Fits to the primary vs. secondary colors are identical to within measurement uncertainties. The color range of binaries as a group is indistinguishable from that of the larger population of apparently single transneptunian objects. Whatever mechanism produced the colors of apparently single TNOs acted equally on binary systems. The most likely explanation is that the colors of transneptunian objects and binaries alike are primordial and indicative of their origin in a locally homogeneous, globally heterogeneous protoplanetary disk.

We are reporting on a new accurate photographic light curve of Pluto for 1933-1934 when the heliocentric distance was 40 AU. We used 43 B-band and V-band images of Pluto on 32 plates taken on 15 nights from 19 March 1933 to 10 March 1934. Most of these plates were taken with the Mount Wilson 60" and 100" telescopes, but 7 of the plates (now at the Harvard College Observatory) were taken with the 12" and 16" Metcalf doublets at Oak Ridge. The plates were measured with an iris diaphragm photometer, which has an average one-sigma photometric error on these plates of 0.08 mag as measured by the repeatability of constant comparison stars. The modern B and V magnitudes for the comparison stars were measured with the Lowell Observatory Hall 1.1-m telescope. The magnitudes in the plate's photographic system were converted to the Johnson B- and V-system after correction with color terms, even though they are small in size. We find that the average B-band mean opposition magnitude of Pluto in 1933-1934 was 15.73 +- 0.01, and we see a roughly sinusoidal modulation on the rotational period (6.38 days) with a peak-to-peak amplitude of 0.11 +- 0.03 mag. With this, we show that Pluto darkened by 5% from 1933-1934 to 1953-1955. This darkening from 1933-1934 to 1953-1955 cannot be due to changing viewing geometry (as both epochs had identical sub-Earth latitudes), so our observations must record a real albedo change over the southern hemisphere. The later darkening trend from 1954 to the 1980s has been explained by changing viewing geometry (as more of the darker northern hemisphere comes into view). Thus, we now have strong evidence for albedo changes on the surface of Pluto, and these are most easily explained by the systematic sublimation of frosts from the sunward pole that led to a drop in the mean surface albedo.

A four-body orbit solution for the Pluto system yields GM values of 870.3 +/- 3.7, 101.4 +/- 2.8, 0.039 +/- 0.034, and 0.021 +/- 0.042 km3 sec-2 for Pluto, Charon, Nix, and Hydra, respectively. Assuming a Charon-like density of 1.63 gm cm-3, the implied diameters for Nix and Hydra are 88 and 72 km, leading to visual geometric albedos of 0.08 and 0.18, respectively, though with considerable uncertainty. The eccentricity of Charon's orbit has a significant nonzero value; however, the 0.030 +/- 0.009 deg yr-1 rate at which the line of apsides precesses is insufficient to explain the difference in the longitude of periapsis seen in the orbits fitted to the 1992-1993 and 2002-2003 data sets. The mean orbital periods for Hydra, Nix, and Charon are in the ratios of 6.064 +/- 0.006 : 3.991 +/- 0.007 : 1, but we have not identified any resonant arguments that would indicate the existence of a mean motion resonance between any pairs of satellites.

We present 0.8 to 2.4 micron spectral observations of uranian satellites, obtained at IRTF/SpeX on 17 nights during 2001-2005. The spectra reveal for the first time the presence of CO2 ice on the surfaces of Umbriel and Titania, by means of 3 narrow absorption bands near 2 microns. Several additional, weaker CO2 ice absorptions have also been detected. No CO2 absorption is seen in Oberon spectra, and the strengths of the CO2 ice bands decline with planetocentric distance from Ariel through Titania. We use the CO2 absorptions to map the longitudinal distribution of CO2 ice on Ariel, Umbriel, and Titania, showing that it is most abundant on their trailing hemispheres. We also examine H2O ice absorptions in the spectra, finding deeper H2O bands on the leading hemispheres of Ariel, Umbriel, and Titania, but the opposite pattern on Oberon. Potential mechanisms to produce the observed longitudinal and planetocentric distributions of the two ices are considered.

We describe IRAC 3.6-8 mu_m observations and ground-based near-IR JHKs photometry from Mimir and 2MASS of the massive double cluster h & chi Persei complete to J=15.5 (M ~ 1.3 Msun. Within 25' of the cluster centers we detect \~11,000 sources with J < 15.5, ~7,000 sources with [4.5] < 15, and ~ 5000 sources with [8] < 14.5. In both clusters, the surface density profiles derived from the 2MASS data decline with distance from the cluster centers as expected for a bound cluster. Within 15' of the cluster centers, ~ 50% of the stars lie on a reddened ~ 13 Myr isochrone; at 15'-25' from the cluster centers, ~ 40% lie on this isochrone. Thus, the optical/2MASS color-magnitude diagrams indicate that h & chi Per are accompanied by a halo population with roughly the same age and distance as the two dense clusters. The double cluster lacks any clear IR excess sources for J < 13.5 (~ 2.7 Msun). Therefore, disks around high-mass stars disperse prior to ~ 10^{7} yr. At least 2-3% of the fainter cluster stars have strong IR excess at both [5.8] and [8]. About 4-8% of sources slightly more massive than the Sun (~ 1.4 Msun) have IR excesses at [8]. Combined with the lack of detectable excesses for brighter stars, this result suggests that disks around lower-mass stars have longer lifetimes. The IR excess population also appears to be larger at longer IRAC bands ([5.8], [8]) than at shorter IRAC/2MASS bands ($K_{s}$, [4.5]), a result consistent with an inside-out clearing of disks.

The Edgeworth-Kuiper belt encodes the dynamical history of the outer solar system. Kuiper belt objects (KBOs) bear witness to coagulation physics, the evolution of planetary orbits, and external perturbations from the solar neighborhood. We critically review the present-day belt's observed properties and the theories designed to explain them. Theories are organized according to a possible time-line of events. In chronological order, epochs described include (1) coagulation of KBOs in a dynamically cold disk, (2) formation of binary KBOs by fragmentary collisions and gravitational captures, (3) stirring of KBOs by Neptune-mass planets (``oligarchs''), (4) eviction of excess oligarchs, (5) continued stirring of KBOs by remaining planets whose orbits circularize by dynamical friction, (6) planetary migration and capture of Resonant KBOs, (7) creation of the inner Oort cloud by passing stars in an open stellar cluster, and (8) collisional comminution of the smallest KBOs. Recent work underscores how small, collisional, primordial planetesimals having low velocity dispersion permit the rapid assembly of ~5 Neptune-mass oligarchs at distances of 15-25 AU. We explore the consequences of such a picture. We propose that Neptune-mass planets whose orbits cross into the Kuiper belt for up to ~20 Myr help generate the high-perihelion members of the hot Classical disk and Scattered belt. By contrast, raising perihelia by sweeping secular resonances during Neptune's migration might fill these reservoirs too inefficiently when account is made of how little primordial mass might reside in bodies having sizes of order 100 km. These and other frontier issues in trans-Neptunian space are discussed quantitatively.

Observations of Pluto and its solar-tidal stability zone were made using the Advanced Camera for Surveys' (ACS) Wide Field Channel (WFC) on the Hubble Space Telescope on UT 2005 May 15 and UT 2005 May 18. Two small satellites of Pluto, provisionally designated S/2005 P 1 and S/2005 P 2, were discovered, as discussed by Weaver et al. (2006) and Stern et al. (2006a). Confirming observations of the newly discovered moons were obtained using the ACS in the High Resolution Channel (HRC) mode on 2006 Feb 15 (Mutchler et al. 2006). Both sets of observations provide strong constraints on the existence of any additional satellites in the Pluto system. Based on the May 2005 observations using the ACS/WFC, we place a 90%-confidence lower limit of m_V = 26.8 (m_V = 27.4 for a 50%-confidence lower limit) on the magnitude of undiscovered satellites greater than 5" (1.1x10^5 km) from Pluto. Using the 2005 Feb 15 ACS/HRC observations we place 90%-confidence lower limits on the apparent magnitude of any additional satellites of m_V = 26.4 between 3"-5" (6.9x10^4-1.1x10^5 km) from Pluto, m_V = 25.7 between 1"-3" (2.3x10^4-6.9x10^4 km) from Pluto, and m_V = 24. between 0.3"-1" (6.9x10^3-2.3x10^4 km) from Pluto. The 90%-confidence magnitude limits translate into upper limits on the diameters of undiscovered satellites of 29 km outside of 5" from Pluto, 36 km between 3"-5" from Pluto, 49 km between 1"-3" from Pluto, and 115 km between 0.3"-1" for a comet-like albedo of p_V = 0.04. If potential satellites are assumed to have a Charon-like albedo of p_V = 0.38, the diameter limits are 9 km, 12 km, 16 km, and 37 km, respectively.

We present new astrometry of Pluto's three satellites from images taken of the Pluto system during 2002-3 with the High Resolution Camera (HRC) mode of the Advanced Camera for Surveys (ACS) instrument on the Hubble Space Telescope. The observations were designed to produce an albedo map of Pluto but they also contain images of Charon and the two recently discovered satellites, S/2005 P1 and S/2005 P2. Orbits fitted to all three satellites are nearly co-planar and, for Charon and P2, have eccentricities consistent with zero. The orbit of the outermost satellite, P1, has a significant eccentricity of 0.0052 +/- 0.0011. Orbital periods of P1, P2, and Charon are 38.2065 +/- 0.0014, 24.8562 +/- 00013, and 6.3872304 +/- 0.0000011 days, respectively. The total system mass based on Charon's orbit is 1.4570 +/- 0.0009 x 10^22 kg. We confirm previous results that orbital periods are close to the ratio of 6:4:1 (P1:P2:Charon) indicative of mean-motion resonances, but our results formally preclude precise integer period ratios. The orbits of P1 and P2, being about the barycenter rather than Pluto, enable us to measure the Charon/Pluto mass ratio as 0.1165 +/- 0.0055. This new mass ratio implies a density of 1.65 +/- 0.06 g cm^-3 for Charon (603.6 km radius) and 2.03 +/- 0.06 g cm^-3 for Pluto (1153 km radius) thus adding confirmation that Charon is significantly less dense than Pluto. Finally, by stacking all images, we can extract globally averaged photometry. P1 has a mean opposition magnitude of V=24.39 +/- 0.09 and color of (B-V) = 0.64 +/- 0.12. P2 has a mean opposition magnitude of V=24.55 +/- 0.10 and color of (B-V) = 0.91 +/- 0.15.

Pluto's first known moon, Charon, was discovered in 1978 (Christy 1978) and has a diameter about half that of Pluto (Buie 1992,Young 1994, Sicardy 2005), which makes it larger relative to its primary than any other moon in the Solar System. Previous searches for other satellites around Pluto have been unsuccessful (Stern 1991, Stern 1994, Stern 2003), but they were not sensitive to objects <=150 km in diameter and there are no fundamental reasons why Pluto should not have more satellites (Stern 1994). Here we report the discovery of two additional moons around Pluto, provisionally designated S/2005 P1 (hereafter P1) and S/2005 P2 (hereafter P2), which makes Pluto the first Kuiper belt object (KBO) known to have multiple satellites. These new satellites are much smaller than Charon (diameter~1200 km), with P1 ranging in diameter from 60-165 km depending on the surface reflectivity, and P2 about 20% smaller than P1. Although definitive orbits cannot be derived, both new satellites appear to be moving in circular orbits in the same orbital plane as Charon, with orbital periods of ~38 days (P1) and ~25 days (P2). The implications of the discovery of P1 and P2 for the origin and evolution of the Pluto system, and for the satellite formation process in the Kuiper belt, are discussed in a companion paper (Stern 2006).

Our discovery of two new satellites of Pluto, designated S/2005 P 1 and S/2005 P 2 (henceforth, P1 and P2), combined with the constraints on the absence of more distant satellites of Pluto, reveal that Pluto and its moons comprise an unusual, highly compact, quadruple system. The two newly discovered satellites of Pluto have masses that are very small compared to both Pluto and Charon, creating a striking planet-satellite system architecture. These facts naturally raise the question of how this puzzling satellite system came to be. Here we show that P1 and P2's proximity to Pluto and Charon, along with their apparent locations in high-order mean-motion resonances, likely result from their being constructed from Plutonian collisional ejecta. We argue that variable optical depth dust-ice rings form sporadically in the Pluto system, and that rich satellite systems may be found, perhaps frequently, around other large Kuiper Belt objects.

We review ongoing efforts to identify occupants of mean-motion resonances (MMRs) and collisional families in the Edgeworth-Kuiper belt. Direct integrations of trajectories of Kuiper belt objects (KBOs) reveal the 1:1 (Trojan), 5:4, 4:3, 3:2 (Plutino), 5:3, 7:4, 9:5, 2:1 (Twotino), and 5:2 MMRs to be inhabited. Apart from the Trojan, resonant KBOs typically have large orbital eccentricities and inclinations. The observed pattern of resonance occupation is consistent with resonant capture and adiabatic excitation by a migratory Neptune; however, the dynamically cold initial conditions prior to resonance sweeping that are typically assumed by migration simulations are probably inadequate. Given the dynamically hot residents of the 5:2 MMR and the substantial inclinations observed in all exterior MMRs, a fraction of the primordial belt was likely dynamically pre-heated prior to resonance sweeping. A pre-heated population may have arisen as Neptune gravitationally scattered objects into trans-Neptunian space. The spatial distribution of Twotinos offers a unique diagnostic of Neptune's migration history. The Neptunian Trojan population may rival the Jovian Trojan population, and the former's existence is argued to rule out violent orbital histories for Neptune. Finally, lowest-order secular theory is applied to several hundred non-resonant KBOs with well-measured orbits to update proposals of collisional families. No convincing family is detected.

The Deep Ecliptic Survey is a project whose goal is to survey a large area of the near-ecliptic region to a faint limiting magnitude (R ~ 24) in search of objects in the outer solar system. We are collecting a large homogeneous data sample from the Kitt Peak Mayall 4-m and Cerro Tololo Blanco 4-m telescopes with the Mosaic prime-focus CCD cameras. Our goal is to collect a sample of 500 objects with good orbits to further our understanding of the dynamical structure of the outer solar system. This survey has been in progress since 1998 and is responsible for 272 designated discoveries as of March 2003. We summarize our techniques, highlight recent results, and describe publically available resources.

As part of our ongoing Deep Ecliptic Survey (DES) of the Kuiper belt, we report on the occupation of the 1:1 (Trojan), 4:3, 3:2, 7:4, 2:1, and 5:2 Neptunian mean-motion resonances (MMRs). The occupation of the 1:1 and 5:2 MMRs is not easily understood within the standard model of resonance sweeping by a migratory Neptune over an initially dynamically cold belt. Our dynamically hot, 5:2 resonant objects can librate with modest amplitudes of 90 deg within the resonance for at least 1 Gyr. Their trajectories cannot be explained by close encounters with Neptune alone, given the latter's current orbit. The hot orbits of such 5:2 resonant KBOs, unlike hot orbits of previously known resonant KBOs, may imply that these objects were pre-heated to large inclination and large eccentricity prior to resonance capture by a migratory Neptune. Our first discovered Neptunian Trojan, 2001QR322, may not owe its existence to Neptune's migration at all. The trajectory of 2001QR322 is remarkably stable; the object can undergo tadpole-type libration about Neptune's leading Lagrange (L4) point for at least 1 Gyr with a libration amplitude of 24 deg. Trojan capture probably occurred while Neptune accreted the bulk of its mass. For an assumed albedo of 12--4%, our Trojan is 130--230 km in diameter. Model-dependent estimates place the total number of Neptune Trojans resembling 2001QR322 at 20--60. Their existence might rule out violent orbital histories for Neptune.

Images of the trans-Neptunian objects 1997 CQ29 and 2000 CF105 obtained with the Hubble Space Telescope's WFPC2 camera show them to be binary. The two components of 1997 CQ29 were separated in our images by 0.20 arcsec in November 2001 and by 0.33 arcsec in June/July 2002. The corresponding minimum physical distances are 6100 km and 10,200 km. The companion to 2000 CF105 was 0.78 arcsec from the primary, at least 23,400 km. Six other objects in the trans-Neptunian region, including Pluto and its moon Charon, are known to be binaries; 1997 CQ29 and 2000 CF105 are the seventh and eighth known pair. Binarity appears to be a not-uncommon characteristic in this region of the solar system, with detectable companions present in 4% of the objects we have examined.