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

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

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.

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.

New Horizons observations suggest that CH4 on Pluto has a complex history, involving reservoirs of different composition, thickness and stability controlled by volatile processes occurring on different timescales. In order to interpret these observations, we use a Pluto volatile transport model able to simulate the cycles of N2 and CH4 ices over millions of years. By assuming fixed solid mixing ratios, we explore how changes in surface albedos, emissivities and thermal inertias impact volatile transport. This work is therefore a direct and natural continuation of the work by Bertrand et al. (2018), which only explored the N2 cycles. Results show that bright CH4 deposits can create cold traps for N2 ice outside Sputnik Planitia, leading to a strong coupling between the N2 and CH4 cycles. Depending on the assumed albedo for CH4 ice, the model predicts CH4 ice accumulation (1) at the same equatorial latitudes where the Bladed Terrain Deposits are observed, supporting the idea that these CH4-rich deposits are massive and perennial, or (2) at mid-latitudes (25{\deg}N-70{\deg}N), forming a thick mantle which is consistent with New Horizons observations. In our simulations, both CH4 ice reservoirs are not in an equilibrium state and either one can dominate the other over long timescales, depending on the assumptions made for the CH4 albedo. This suggests that long-term volatile transport exists between the observed reservoirs. The model also reproduces the formation of N2 deposits at mid-latitudes and in the equatorial depressions surrounding the Bladed Terrain, as observed by New Horizons. At the poles, only seasonal CH4 and N2 deposits are obtained in Pluto's current orbital configuration. Finally, we show that Pluto's atmosphere always contained, over the last astronomical cycles, enough gaseous CH4 to absorb most of the incoming Lyman-flux.

Observations of Pluto's surface made by the New Horizons spacecraft indicates present-day nitrogen ice glaciation in and around the basin known as Sputnik Planum. Motivated by these observations, we have developed an evolutionary glacial flow model of solid nitrogen ice taking into account its published thermophysical and rheologies properties. This model assumes that glacial ice layers flow laminarly and have low aspect ratios which permits a vertically integrated mathematical formulation. We assess the conditions for the validity of laminar nitrogen ice motion by revisiting the problem of the onset of solid-state buoyant convection of nitrogen ice for a variety of bottom thermal boundary conditions. Subject to uncertainties in nitrogen ice rheology, nitrogen ice layers are estimated to flow laminarly for thicknesses less than 400-1000 meters. The resulting mass-flux formulation for when the nitrogen ice flows as a laminar dry glacier is characterized by an Arrhenius-Glen functional form. The flow model developed is used here to qualitatively answer some questions motivated by observed glacial flow features found on Sputnik Planum. We find that the wavy transverse dark features found along the northern shoreline of Sputnik Planum may be a transitory imprint of shallow topography just beneath the ice surface suggesting the possibility that a major shoreward flow event happened relatively recently within the last few hundred years. Model results also support the interpretation that the prominent darkened features resembling flow lobes observed along the eastern shoreline of the Sputnik Planum basin may be a result of wet nitrogen glacial ice flowing into the basin from the pitted highlands of eastern Tombaugh Regio.

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.

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.