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.

We use data from Cassini's Composite Infrared Spectrometer CIRS to map Mimas' surface temperatures and its thermophysical properties. This provides a dramatic improvement on the work in Howett et al. (2011), where the values were determined at only two regions on Mimas (one inside, and another outside of the anomalous region). We use all spatially-resolved scans made by CIRS' focal plane 3 (FP3, 600 to 1100 cm-1) of Mimas' surface, which are largely daytime observations but do include one nighttime one. The resulting temperature maps confirm the presence and location of Mimas' previously discovered thermally anomalous region. No other thermally anomalous regions were discovered, although we note that the surface coverage is incomplete on Mimas' leading and anti-Saturn hemisphere. The thermal inertia map confirms that the anomalous region has a notably higher thermal inertia than its surroundings: 98 +/- 42 J m-2 K-1 s-1/2 inside of the anomaly, compared to 34 +/- 32 J m-2 K-1 s-1/2 outside. The albedo inside and outside of the anomalous region agrees within their uncertainty: 0.45+/-0.08 inside compared to 0.41 +/- 0.07 outside the anomaly. Interestingly the albedo appears brighter inside the anomaly region, which may not be surprising given this region does appear brighter at some UV wavelengths (0.338 microns, see Schenk et al., 2011). However, this result should be treated with caution because, as previously stated, statistically the albedo of these two regions is the same when their uncertainties are considered. These thermal inertia and albedo values determined here are consistent with those found by Howett et al. (2011), who determined the thermal inertia inside the anomaly to be 66 +/- 23 J m-2 K-1 s-1/2 and less than 16 J m-2 K-1 s-1/2 outside, with albedos that varied from 0.49 to 0.70.

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.

Motivated by the Kepler K2 time series of Titan, we present an aperture optimization technique for extracting photometry of saturated moving targets with high temporally- and spatially-varying backgrounds. Our approach uses $k$-means clustering to identify interleaved families of images with similar Point-Spread Function and saturation properties, optimizes apertures for each family independently, then merges the time series through a normalization procedure. By applying $k$-means aperture optimization to the K2 Titan data, we achieve $\leq$0.33% photometric scatter in spite of background levels varying from 15% to 60% of the target's flux. We find no compelling evidence for signals attributable to atmospheric variation on the timescales sampled by these observations. We explore other potential applications of the $k$-means aperture optimization technique, including testing its performance on a saturated K2 eclipsing binary star. We conclude with a discussion of the potential for future continuous high-precision photometry campaigns for revealing the dynamical properties of Titan's atmosphere.

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].

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.

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.

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.

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.

On July 14th 2015, NASA's New Horizons mission gave us an unprecedented detailed view of the Pluto system. The complex compositional diversity of Pluto's encounter hemisphere was revealed by the Ralph/LEISA infrared spectrometer on board of New Horizons. We present compositional maps of Pluto defining the spatial distribution of the abundance and textural properties of the volatiles methane and nitrogen ices and non-volatiles water ice and tholin. These results are obtained by applying a pixel-by-pixel Hapke radiative transfer model to the LEISA scans. Our analysis focuses mainly on the large scale latitudinal variations of methane and nitrogen ices and aims at setting observational constraints to volatile transport models. Specifically, we find three latitudinal bands: the first, enriched in methane, extends from the pole to 55deg N, the second dominated by nitrogen, continues south to 35deg N, and the third, composed again mainly of methane, reaches 20deg N. We demonstrate that the distribution of volatiles across these surface units can be explained by differences in insolation over the past few decades. The latitudinal pattern is broken by Sputnik Planitia, a large reservoir of volatiles, with nitrogen playing the most important role. The physical properties of methane and nitrogen in this region are suggestive of the presence of a cold trap or possible volatile stratification. Furthermore our modeling results point to a possible sublimation transport of nitrogen from the northwest edge of Sputnik Planitia toward the south.

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.

Cassini's Composite Infrared Spectrometer (CIRS) observed both of Rhea's polar regions during two flybys on 2013/03/09 and 2015/02/10. The results show Rhea's southern winter pole is one of the coldest places directly observed in our solar system: temperatures of 25.4+/-7.4 K and 24.7+/-6.8 K are inferred. The surface temperature of the northern summer pole is warmer: 66.6+/-0.6 K. Assuming the surface thermophysical properties of both polar regions are comparable then these temperatures can be considered a summer and winter seasonal temperature constraint for the polar region. These observations provide solar longitude coverage at 133 deg and 313 deg for the summer and winter poles respectively, with additional winter temperature constraint at 337 deg. Seasonal models with bolometric albedos of 0.70-0.74 and thermal inertias of 1-46 MKS can provide adequate fits to these temperature constraints. Both these albedo and thermal inertia values agree (within error) with those previously observed on both Rhea's leading and trailing hemispheres. Investigating the seasonal temperature change of Rhea's surface is particularly important, as the seasonal wave is sensitive to deeper surface temperatures (~10cm to m) than the more commonly reported diurnal wave (<1cm). The low thermal inertia derived here implies that Rhea's polar surfaces are highly porous even at great depths. Analysis of a CIRS 10 to 600 cm-1 stare observation, taken between 16:22:33 and 16:23:26 UT on 2013/03/09 centered on 71.7 W, 58.7 S provides the first analysis of a thermal emissivity spectrum on Rhea. The results show a flat emissivity spectrum with negligible emissivity features. A few possible explanations exist for this flat emissivity spectrum, but the most likely for Rhea is that the surface is both highly porous and composed of small particles (less than approximately 50 um).

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.

We present the distance measurement to z = 0.32 using the 11th data release of the Sloan Digital Sky Survey-III Baryon Acoustic Oscillation Survey (BOSS). We use 313,780 galaxies of the low-redshift (LOWZ) sample over 7,341 square-degrees to compute $D_V = (1264 \pm 25)(r_d/r_{d,fid})$ - a sub 2% measurement - using the baryon acoustic feature measured in the galaxy two-point correlation function and power-spectrum. We compare our results to those obtained in DR10. We study observational systematics in the LOWZ sample and quantify potential effects due to photometric offsets between the northern and southern Galactic caps. We find the sample to be robust to all systematic effects found to impact on the targeting of higher-redshift BOSS galaxies, and that the observed north-south tensions can be explained by either limitations in photometric calibration or by sample variance, and have no impact on our final result. Our measurement, combined with the baryonic acoustic scale at z = 0.57, is used in Anderson et al. (2013a) to constrain cosmological parameters.