We examine 141 N-body simulations of terrestrial planet late-stage accretion that use the Grand Tack scenario, coupling the collisional results with a hafnium-tungsten (Hf-W) isotopic evolution model. Accretion in the Grand Tack scenario results in faster planet formation than classical accretion models because of higher planetesimal surface density induced by a migrating Jupiter. Planetary embryos that grow rapidly experience radiogenic ingrowth of mantle $^{182}$W that is inconsistent with the measured terrestrial composition, unless much of the tungsten is removed by an impactor core that mixes thoroughly with the target mantle. For physically Earth-like surviving planets, we find that the fraction of equilibrating impactor core $k_\text{core} \geq 0.6$ is required to produce results agreeing with observed terrestrial tungsten anomalies (assuming equilibration with relatively large volumes of target mantle material; smaller equilibrating mantle volumes would require even larger $k_\text{core}$). This requirement of substantial core re-equilibration may be difficult to reconcile with fluid dynamical predictions and hydrocode simulations of mixing during large impacts, and hence this result does not favor the rapid planet building that results from Grand Tack accretion.

Silicon and Mg in differentiated rocky bodies exhibit heavy isotope enrichments that have been attributed to evaporation of partially or entirely molten planetesimals. We evaluate the mechanisms of planetesimal evaporation in the early solar system and the conditions that controlled attendant isotope fractionations. Energy balance at the surface of a body accreted within ~1 Myr of CAI formation and heated from within by 26Al decay results in internal temperatures exceeding the silicate solidus, producing a transient magma ocean with a thin surface boundary layer of order < 1 meter that would be subject to foundering. Bodies that are massive enough to form magma oceans by radioisotope decay (ge 0.1%) can retain hot rock vapor even in the absence of ambient nebular gas. We find that a steady-state rock vapor forms within minutes to hours and results from a balance between rates of magma evaporation and atmospheric escape. Vapor pressure buildup adjacent to the surfaces of the evaporating magmas would have inevitably led to an approach to equilibrium isotope partitioning between the vapor phase and the silicate melt. Numerical simulations of this near-equilibrium evaporation process for a body with a radius of ~ 700 km yield a steady-state far-field vapor pressure corresponding to 95% saturation. Approaches to equilibrium isotope fractionation between vapor and melt should have been the norm during planet formation due to the formation of steady-state rock vapor atmospheres and/or the presence of protostellar gas. We model the Si and Mg isotopic composition of bulk Earth and show that the best fit is for a carbonaceous chondrite-like source material with about 12% loss of Mg and 15% loss of Si resulting from near-equilibrium evaporation into the solar protostellar disk of hydrogen gas on timescales of 10,000 to 100,000 years.

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.

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.

As a long-term energy source, tidal heating in subsurface oceans of icy satellites can influence their thermal, rotational, and orbital evolution, and the sustainability of oceans. We present a new theoretical treatment for tidal heating in thin subsurface oceans with overlying incompressible elastic shells of arbitrary thickness. The stabilizing effect of an overlying shell damps ocean tides, reducing tidal heating. This effect is more pronounced on Enceladus than on Europa because the effective rigidity on a small body like Enceladus is larger. For the range of likely shell and ocean thicknesses of Enceladus and Europa, the thin shell approximation of Beuthe (2016) is generally accurate to less than about 4%.The time-averaged surface distribution of ocean tidal heating is distinct from that due to dissipation in the solid shell, with higher dissipation near the equator and poles for eccentricity and obliquity forcing respectively. This can lead to unique horizontal shell thickness variations if the shell is conductive. The surface displacement driven by eccentricity and obliquity forcing can have a phase lag relative to the forcing tidal potential due to the delayed ocean response. For Europa and Enceladus, eccentricity forcing generally produces greater tidal amplitudes due to the large eccentricity values relative to the obliquity values. Despite the small obliquity values, obliquity forcing generally produces larger phase lags due to the generation of Rossby-Haurwitz waves. If Europa's shell and ocean are respectively 10 and 100 km thick, the tide amplitude and phase lag are 26.5 m and $<1$ degree for eccentricity forcing, and $<2.5$ m and $<18$ degrees for obliquity forcing. Measurement of the obliquity phase lag (e.g. by Europa Clipper) would provide a probe of ocean thickness

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.

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.

Because of the high energies involved, giant impacts that occur during planetary accretion cause large degrees of melting. The depth of melting in the target body after each collision determines the pressure and temperature conditions of metal-silicate equilibration and thus geochemical fractionation that results from core-mantle differentiation. The accretional collisions involved in forming the terrestrial planets of the inner Solar System have been calculated by previous studies using N-body accretion simulations. Here we use the output from such simulations to determine the volumes of melt produced and thus the pressure and temperature conditions of metal-silicate equilibration, after each impact, as Earth-like planets accrete. For these calculations a parametrised melting model is used that takes impact velocity, impact angle and the respective masses of the impacting bodies into account. The evolution of metal-silicate equilibration pressures (as defined by evolving magma ocean depths) during Earth's accretion depends strongly on the lifetime of impact-generated magma oceans compared to the time interval between large impacts. In addition, such results depend on starting parameters in the N-body simulations, such as the number and initial mass of embryos. Thus, there is the potential for combining the results, such as those presented here, with multistage core formation models to better constrain the accretional history of the Earth.

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 constructed a 6-degrees of freedom rotational model of Titan as a 3-layer body consisting of a rigid core, a fluid global ocean, and a floating ice shell. The ice shell exhibits partially-compensated lateral thickness variations in order to simultaneously match the observed degree-two gravity and shape coefficients. The rotational dynamics are affected by the gravitational torque of Saturn, the gravitational coupling between the inner core and the shell, and the pressure coupling at the fluid-solid boundaries. Between 10 and 13% of our model Titans have an obliquity (due to a resonance with the 29.5-year periodic annual forcing) that is consistent with the observed value. The shells of the successful models have a mean thickness of 130 to 140 km, and an ocean of ~250 km thickness. Our simulations of the obliquity evolution show that the Cassini obliquity measurement is an instantaneous one, and does not represent a mean value. Future measurements of the time derivative of the obliquity would help to refine the interior models. We expect in particular a variation of roughly 7 arcmin over the duration of the Cassini mission.

We report results of a Hubble Space Telescope (HST) campaign with the Advanced Camera for Surveys to observe Europa at eastern elongation, i.e. Europa's leading side, on 2008 June 29. With five consecutive HST orbits, we constrain Europa's atmospheric \ion{O}{1} 1304 \A and \ion{O}{1} 1356 \A emissions using the prism PR130L. The total emissions of both oxygen multiplets range between 132 $\pm$ 14 and 226 $\pm$ 14 Rayleigh. An additional systematic error with values on the same order as the statistical errors may be due to uncertainties in modelling the reflected light from Europa's surface. The total emission also shows a clear dependence of Europa's position with respect to Jupiter's magnetospheric plasma sheet. We derive a lower limit for the O$_2$ column density of 6 $\times$ 10$^{18}$ m$^{-2}$. Previous observations of Europa's atmosphere with STIS in 1999 of Europa's trailing side show an enigmatic surplus of radiation on the anti-Jovian side within the disk of Europa. With emission from a radially symmetric atmosphere as a reference, we searched for an anti-Jovian vs sub-Jovian asymmetry with respect to the central meridian on the leading side, and found none. Likewise, we searched for departures from a radially symmetric atmospheric emission and found an emission surplus centered around 90 degree west longitude, for which plausible mechanisms exist. Previous work about the possibility of plumes on Europa due to tidally-driven shear heating found longitudes with strongest local strain rates which might be consistent with the longitudes of maximum UV emissions. Alternatively, asymmetries in Europa's UV emission can also be caused by inhomogeneous surface properties, inhomogeneous solar illuminations, and/or by Europa's complex plasma interaction with Jupiter's magnetosphere.

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.