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 Mars Reconnaissance Orbiter (MRO) spacecraft has been in orbit around Mars since March 2006. The Context Camera (CTX) on MRO has returned over 100,000 images of the planet at approximately 5-6 meters per pixel, providing nearly global coverage. During that time, Mars has gone through nearly 7 of its own years, changing solar distance from 1.38 to 1.67 AU and the corresponding solar flux by 45% due to its orbital eccentricity. Seasonal effects and transient phenomena affect atmospheric transparency. Combined with an aging detector, CTX images are difficult to mosaic seamlessly, for all of these changes prevent equalizing images to create a visually smoothly illuminated product. We have developed a method, based on previous work by other researchers for other datasets, to mitigate almost all photometric variations between images in order to create the appearance of an evenly illuminated, practically seamless mosaic. We describe how the process works, which uses a reference source to tie brightness values, and demonstrate its effects across Mars' surface. While the workflow developed for this product is applicable to other planetary bodies, it requires a reference source, which may not yet exist.

The flyby of Pluto and Charon by the New Horizons spacecraft provided high-resolution images of cratered surfaces embedded in the Kuiper belt, an extensive region of bodies orbiting beyond Neptune. Impact craters on Pluto and Charon were formed by collisions with other Kuiper belt objects (KBOs) with diameters from ~40 kilometers to ~300 meters, smaller than most KBOs observed directly by telescopes. We find a relative paucity of small craters less than approximately 13 kilometers in diameter, which cannot be explained solely by geological resurfacing. This implies a deficit of small KBOs (less than 1 to 2 kilometers in diameter). Some surfaces on Pluto and Charon are likely greater than 4 billion years old, thus their crater records provide information on the size-frequency distribution of KBOs in the early Solar System.

The dynamics of Mars' obliquity are believed to be chaotic, and the historical ~3.5 Gyr (late-Hesperian onward) obliquity probability density function (PDF) is high uncertain and cannot be inferred from direct simulation alone. Obliquity is also a strong control on post-Noachian Martian climate, enhancing the potential for equatorial ice/snow melting and runoff at high obliquities (> 40{\deg}) and enhancing the potential for desiccation of deep aquifers at low obliquities (< 25{\deg}). We developed a new technique using the orientations of elliptical craters to constrain the true late-Hesperian-onward obliquity PDF. To do so, we developed a forward model of the effect of obliquity on elliptic crater orientations using ensembles of simulated Mars impactors and ~3.5 Gyr-long Mars obliquity simulations. In our model, the inclinations and speeds of Mars crossing objects bias the preferred orientation of elliptic craters which are formed by low-angle impacts. Comparison of our simulation predictions with a validated database of elliptic crater orientations allowed us to invert for best-fitting obliquity history. We found that since the onset of the late-Hesperian, Mars' mean obliquity was likely low, between ~10{\deg} and ~30{\deg}, and the fraction of time spent at high obliquities > 40{\deg} was likely < 20%.

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.

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 exploration of the binary Pluto-Charon and its small satellites during the New Horizons flyby in 2015 revealed not only widespread geologic and compositional diversity across Pluto, but surprising complexity, a wide range of surface unit ages, evidence for widespread activity stretching across billion of years to the near-present, as well as numerous atmospheric puzzles, and strong atmospheric coupling with its surface. New Horizons also found an unexpected diversity of landforms on its binary companion, Charon. Pluto's four small satellites yielded surprises as well, including their unexpected rapid and high obliquity rotation states, high albedos, and diverse densities. Here we briefly review the findings made by New Horizons and the case for a follow up mission to investigate the Pluto system in more detail. As the next step in the exploration of this spectacular planet-satellite system, we recommend an orbiter to study it in considerably more detail, with new types of instrumentation, and to observe its changes with time. We further call for the in-depth study of Pluto orbiter missions as a precursor to the 2023 Planetary Science Decadal Survey.

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.

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

The identification of impact craters on planetary surfaces provides important information about their geological history. Most studies have relied on individual analysts who map and identify craters and interpret crater statistics. However, little work has been done to determine how the counts vary as a function of technique, terrain, or between researchers. Furthermore, several novel internet-based projects ask volunteers with little to no training to identify craters, and it was unclear how their results compare against the typical professional researcher. To better understand the variation among experts and to compare with volunteers, eight professional researchers have identified impact features in two separate regions of the moon. Small craters (diameters ranging from 10 m to 500 m) were measured on a lunar mare region and larger craters (100s m to a few km in diameter) were measured on both lunar highlands and maria. Volunteer data were collected for the small craters on the mare. Our comparison shows that the level of agreement among experts depends on crater diameter, number of craters per diameter bin, and terrain type, with differences of up to $\sim\pm45%$. We also found artifacts near the minimum crater diameter that was studied. These results indicate that caution must be used in most cases when interpreting small variations in crater size-frequency distributions and for craters $\le10$ pixels across. Because of the natural variability found, projects that emphasize many people identifying craters on the same area and using a consensus result are likely to yield the most consistent and robust information.

An empirical model for forecasting solar wind speed related geomagnetic events is presented here. The model is based on the estimated location and size of solar coronal holes. This method differs from models that are based on photospheric magnetograms (e.g., Wang-Sheeley model) to estimate the open field line configuration. Rather than requiring the use of a full magnetic synoptic map, the method presented here can be used to forecast solar wind velocities and magnetic polarity from a single coronal hole image, along with a single magnetic full-disk image. The coronal hole parameters used in this study are estimated with Kitt Peak Vacuum Telescope He I 1083 nm spectrograms and photospheric magnetograms. Solar wind and coronal hole data for the period between May 1992 and September 2003 are investigated. The new model is found to be accurate to within 10% of observed solar wind measurements for its best one-month periods, and it has a linear correlation coefficient of ~0.38 for the full 11 years studied. Using a single estimated coronal hole map, the model can forecast the Earth directed solar wind velocity up to 8.5 days in advance. In addition, this method can be used with any source of coronal hole area and location data.

An evaluation of the principal uncertainties in the computation of neutrino fluxes produced in cosmic ray showers in the atmosphere is presented. The neutrino flux predictions are needed for comparison with experiment to perform neutrino oscillation studies. The paper concentrates on the main limitations which are due to hadron production uncertainties. It also treats primary cosmic ray flux uncertainties, which are at a lower level. The absolute neutrino fluxes are found to have errors of around 15% in the neutrino energy region important for contained events underground. Large cancellations of these errors occur when ratios of fluxes are considered, in particular, the $\nu_\mu/\bar{\nu}_\mu$ ratio below $E_\nu=1$ GeV, the $(\nu_\mu+\bar{\nu}_\mu)/(\nu_e+\bar{\nu}_e)$ ratio below $E_\nu=10$ GeV and the up/down ratios above $E_\nu=1$ GeV are at the 1% level. A detailed breakdown of the origin of these errors and cancellations is presented.

A Monte-Carlo calculation of the atmospheric neutrino fluxes [1,2] has been extended to take account of the three-dimensional (3D) nature of the problem, including the bending of secondary particles in the geomagnetic field. Emphasis has been placed on minimizing the approximations when introducing the 3D considerations. In this paper, we describe the techniques used and quantify the effects of the small approximations which remain. We compare 3D and 1D calculations using the same physics input in order to evaluate the conditions under which the 3D calculation is required and when the considerably simpler 1D calculation is adequate. We find that the 1D and 3D results are essentially identical for neutrino energy greater than 5 GeV except for small effects in the azimuthal distributions due to bending of the secondary muon by the geomagnetic field during their propagation in the atmosphere.