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.

On January 1st 2019, the New Horizons spacecraft flew by the classical Kuiper belt object (486958) Arrokoth (provisionally designated 2014 MU69), possibly the most primitive object ever explored by a spacecraft. The I/F of Arrokoth is analyzed and fit with a photometric function that is a linear combination of the Lommel-Seeliger (lunar) and Lambert photometric functions. Arrokoth has a geometric albedo of p_V = 0.21_(-0.04)^(+0.05) at a wavelength of 550 nm and ~0.24 at 610 nm. Arrokoth's geometric albedo is greater than the median but consistent with a distribution of cold classical Kuiper belt objects whose geometric albedos were determined by fitting a thermal model to radiometric observations. Thus, Arrokoth's geometric albedo adds to the orbital and spectral evidence that it is a cold classical Kuiper belt object. Maps of the normal reflectance and hemispherical albedo of Arrokoth are presented. The normal reflectance of Arrokoth's surface varies with location, ranging from ~0.10-0.40 at 610 nm with an approximately Gaussian distribution. Both Arrokoth's extrema dark and extrema bright surfaces are correlated to topographic depressions. Arrokoth has a bilobate shape and the two lobes have similar normal reflectance distributions: both are approximately Gaussian, peak at ~0.25 at 610 nm, and range from ~0.10-0.40, which is consistent with co-formation and co-evolution of the two lobes. The hemispherical albedo of Arrokoth varies substantially with both incidence angle and location, the average hemispherical albedo at 610 nm is 0.063 +/- 0.015. The Bond albedo of Arrokoth at 610 nm is 0.062 +/- 0.015.

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.

Kuiper belt object (KBO) Eris is exceptionally bright with a greater visible geometric albedo than any other known KBO. Its infrared reflectance spectrum is dominated by methane, which should form tholins that darken the surface on timescales much shorter than the age of the Solar System. Thus one or more ongoing processes probably maintain its brightness. Eris is predicted to have a primarily nitrogen atmosphere that is in vapor pressure equilibrium with nitrogen-ice and is collisional (not ballistic). Eris's eccentric orbit is expected to result in two atmospheric regimes: (1) a period near perihelion when the atmosphere is global (analogous to the atmospheres of Mars, Triton, and Pluto) and (2) a period near aphelion when only a local atmosphere exists near the warmest region (analogous to the atmosphere of Io). A numerical model developed to simulate Eris's thermal and volatile evolution in the local atmosphere regime is presented. The model conserves energy, mass, and momentum while maintaining vapor pressure equilibrium. It is adaptable to other local, collisional, sublimation atmospheres, which in addition to Io and Eris, may occur on several volatile-bearing KBOs. The model was applied for a limiting case where Eris is fixed at aphelion and has an initial nitrogen-ice mass everywhere equal to the precipitable column of nitrogen in Pluto's atmosphere during the New Horizons encounter (the resultant mass if the Pluto atmosphere collapsed uniformly onto the surface). The model results indicate that (1) transport of nitrogen in the local, collisional, sublimation atmosphere regime is significant, (2) changes of Eris's albedo or color from nitrogen transport may be observable, and (3) uniform collapse of a global, nitrogen atmosphere likely cannot explain Eris's anomalous albedo in the present epoch. Seasonal volatile transport remains a plausible hypothesis to explain...

A search for temporal changes on Pluto and Charon was motivated by (1) the discovery of young surfaces in the Pluto system that imply ongoing or recent geologic activity, (2) the detection of active plumes on Triton during the Voyager 2 flyby, and (3) the abundant and detailed information that observing geologic processes in action provides about the processes. A thorough search for temporal changes using New Horizons images was completed. Images that covered the same region were blinked and manually inspected for any differences in appearance. The search included full-disk images such that all illuminated regions of both bodies were investigated and higher resolution images such that parts of the encounter hemispheres were investigated at finer spatial scales. Changes of appearance between different images were observed but in all cases were attributed to variability of the imaging parameters (especially geometry) or artifacts. No differences of appearance that are strongly indicative of a temporal change were found on the surface or in the atmosphere of either Pluto or Charon. Limits on temporal changes as a function of spatial scale and temporal interval during the New Horizons encounter are determined. The longest time interval constraint is one Pluto/Charon rotation period (~6.4 Earth days). Contrast reversal and high-phase bright features that change in appearance with solar phase angle are identified. The change of appearance of these features is most likely due to the change in phase angle rather than a temporal change. Had active plumes analogous to the plumes discovered on Triton been present on the encounter hemispheres of either Pluto or Charon, they would have been detected. The absence of active plumes may be due to temporal variability (i.e., plumes do occur but none were active on the encounter hemispheres during the epoch of the New Horizons encounter ...

Several trace constituents can reach saturation vapor pressure in Pluto's present day atmosphere. As such, we describe a search for discrete could features in Pluto's atmosphere using New Horizons data obtained on 14-15 July 2015, during the Pluto flyby closest approach. We report that Pluto's present day atmosphere is at least largely (>99% by surface area) free of discrete clouds. We also report a handful of features that may plausibly be clouds, all of which were detected obliquely and at high phase observing geometry. No cloud candidates were identified away from the terminators or in low phase (backscattering geometry) images.

The exploration of the Pluto-Charon system by the New Horizons spacecraft represents the first opportunity to understand the distribution of albedo and other photometric properties of the surfaces of objects in the Solar System's "Third Zone" of distant ice-rich bodies. Images of the entire illuminated surface of Pluto and Charon obtained by the Long Range Reconnaissance Imager (LORRI) camera provide a global map of Pluto that reveals surface albedo variegations larger than any other Solar System world except for Saturn's moon Iapetus. Normal reflectances on Pluto range from 0.08-1.0, and the low-albedo areas of Pluto are darker than any region of Charon. Charon exhibits a much blander surface with normal reflectances ranging from 0.20-0.73. Pluto's albedo features are well-correlated with geologic features, although some exogenous low-albedo dust may be responsible for features seen to the west of the area informally named Tombaugh Regio. The albedo patterns of both Pluto and Charon are latitudinally organized, with the exception of Tombaugh Regio, with darker regions concentrated at the Pluto's equator and Charon's northern pole The phase curve of Pluto is similar to that of Triton, the large moon of Neptune believed to be a captured Kuiper Belt Object (KBO), while Charon's is similar to that of the Moon. Preliminary Bond albedos are 0.25+/-0.03 for Charon and 0.72+/-0.07 for Pluto. Maps of an approximation to the Bond albedo for both Pluto and Charon are presented for the first time. Our work shows a connection between very high albedo (near unity) and planetary activity, a result that suggests the KBO Eris may be currently active.

The DEAP-1 low-background liquid argon detector was used to measure scintillation pulse shapes of electron and nuclear recoil events and to demonstrate the feasibility of pulse-shape discrimination (PSD) down to an electron-equivalent energy of 20 keV. In the surface dataset using a triple-coincidence tag we found the fraction of beta events that are misidentified as nuclear recoils to be $<1.4\times 10^{-7}$ (90% C.L.) for energies between 43-86 keVee and for a nuclear recoil acceptance of at least 90%, with 4% systematic uncertainty on the absolute energy scale. The discrimination measurement on surface was limited by nuclear recoils induced by cosmic-ray generated neutrons. This was improved by moving the detector to the SNOLAB underground laboratory, where the reduced background rate allowed the same measurement with only a double-coincidence tag. The combined data set contains $1.23\times10^8$ events. One of those, in the underground data set, is in the nuclear-recoil region of interest. Taking into account the expected background of 0.48 events coming from random pileup, the resulting upper limit on the electronic recoil contamination is $<2.7\times10^{-8}$ (90% C.L.) between 44-89 keVee and for a nuclear recoil acceptance of at least 90%, with 6% systematic uncertainty on the absolute energy scale. We developed a general mathematical framework to describe PSD parameter distributions and used it to build an analytical model of the distributions observed in DEAP-1. Using this model, we project a misidentification fraction of approx. $10^{-10}$ for an electron-equivalent energy threshold of 15 keV for a detector with 8 PE/keVee light yield. This reduction enables a search for spin-independent scattering of WIMPs from 1000 kg of liquid argon with a WIMP-nucleon cross-section sensitivity of $10^{-46}$ cm$^2$, assuming negligible contribution from nuclear recoil backgrounds.

We present a geomorphologic map of Titan's polar terrains. The map was generated from a combination of Cassini Synthetic Aperture Radar (SAR) and Imaging Science Subsystem imaging products, as well as altimetry, SARTopo and radargrammetry topographic datasets. In combining imagery with topographic data, our geomorphologic map reveals a stratigraphic sequence from which we infer process interactions between units. In mapping both polar regions with the same geomorphologic units, we conclude that processes that formed the terrains of the north polar region also acted to form the landscape we observe at the south. Uniform, SAR-dark plains are interpreted as sedimentary deposits, and are bounded by moderately dissected uplands. These plains contain the highest density of filled and empty lake depressions, and canyons. These units unconformably overlay a basement rock that outcrops as mountains and SAR-bright dissected terrains at various elevations across both poles. All these units are then superposed by surficial units that slope towards the seas, suggestive of subsequent overland transport of sediment. From estimates of the depths of the embedded empty depressions and canyons that drain into the seas, the SAR-dark plains must be >600 m thick in places, though the thickness may vary across the poles. At the lowest elevations of each polar region, there are large seas, which are currently liquid methane/ethane filled at the north and empty at the south. The large plains deposits and the surrounding hillslopes may represent remnant landforms that are a result of previously vast polar oceans, where larger liquid bodies may have allowed for a sustained accumulation of soluble and insoluble sediments, potentially forming layered sedimentary deposits. Coupled with vertical crustal movements, the resulting layers would be of varying solubilities and erosional resistances.

Saturn's moon Enceladus offers a unique opportunity in the search for life and habitable environments beyond Earth, a key theme of the National Research Council's 2013-2022 Decadal Survey. A plume of water vapor and ice spews from Enceladus's south polar region. Cassini data suggest that this plume, sourced by a liquid reservoir beneath the moon's icy crust, contain organics, salts, and water-rock interaction derivatives. Thus, the ingredients for life as we know it-- liquid water, chemistry, and energy sources-- are available in Enceladus's subsurface ocean. We have only to sample the plumes to investigate this hidden ocean environment. We present a New Frontiers class, solar-powered Enceladus orbiter that would take advantage of this opportunity, Testing the Habitability of Enceladus's Ocean (THEO). Developed by the 2015 Jet Propulsion Laboratory Planetary Science Summer School student participants under the guidance of TeamX, this mission concept includes remote sensing and in situ analyses with a mass spectrometer, a sub-mm radiometer-spectrometer, a camera, and two magnetometers. These instruments were selected to address four key questions for ascertaining the habitability of Enceladus's ocean within the context of the moon's geological activity: (1) How are the plumes and ocean connected? (2) Are the abiotic conditions of the ocean suitable for habitability? (3) How stable is the ocean environment? (4) Is there evidence of biological processes? By taking advantage of the opportunity Enceladus's plumes offer, THEO represents a viable, solar-powered option for exploring a potentially habitable ocean world of the outer solar system.

Ethane is expected to be the dominant photochemical product on Titan's surface and, in the absence of a process that sequesters it from exposed surface reservoirs, a major constituent of its lakes and seas. Absorption of Cassini's 2.2 cm radar by Ligeia Mare however suggests that this north polar sea is dominated by methane. In order to explain this apparent ethane deficiency, we explore the possibility that Ligeia Mare is the visible part of an alkanofer that interacted with an underlying clathrate layer and investigate the influence of this interaction on an assumed initial ethane-methane mixture in the liquid phase. We find that progressive liquid entrapment in clathrate allows the surface liquid reservoir to become methane-dominated for any initial ethane mole fraction below 0.75. If interactions between alkanofers and clathrates are common on Titan, this should lead to the emergence of many methane-dominated seas or lakes.

Many current and future dark matter and neutrino detectors are designed to measure scintillation light with a large array of photomultiplier tubes (PMTs). The energy resolution and particle identification capabilities of these detectors depend in part on the ability to accurately identify individual photoelectrons in PMT waveforms despite large variability in pulse amplitudes and pulse pileup. We describe a Bayesian technique that can identify the times of individual photoelectrons in a sampled PMT waveform without deconvolution, even when pileup is present. To demonstrate the technique, we apply it to the general problem of particle identification in single-phase liquid argon dark matter detectors. Using the output of the Bayesian photoelectron counting algorithm described in this paper, we construct several test statistics for rejection of backgrounds for dark matter searches in argon. Compared to simpler methods based on either observed charge or peak finding, the photoelectron counting technique improves both energy resolution and particle identification of low energy events in calibration data from the DEAP-1 detector and simulation of the larger MiniCLEAN dark matter detector.