Rosetta observations of 67P/Churyumov-Gerasimenko (67P) reveal that most changes occur in the fallback-generated smooth terrains, vast deposits of granular material blanketing the comet's northern hemisphere. These changes express themselves both morphologically and spectrally across the nucleus, yet we lack a model that describes their formation and evolution. Here we present a self-consistent model that thoroughly explains the activity and mass loss from Hapi's smooth terrains. Our model predicts the removal of dust via re-radiated solar insolation localized within depression scarps that are substantially more ice-rich than previously expected. We couple our model with numerous Rosetta observations to thoroughly capture the seasonal erosion of Hapi's smooth terrains, where local scarp retreat gradually removes the uppermost dusty mantle. As sublimation-regolith interactions occur on rocky planets, comets, icy moons and KBOs, our coupled model and observations provide a foundation for future understanding of the myriad of sublimation-carved worlds.

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.

The James Webb Space Telescope (JWST), scheduled for launch in 2018, is the successor to the Hubble Space Telescope (HST) but with a significantly larger aperture (6.5 m) and advanced instrumentation focusing on infrared science (0.6-28.0 $\mu$m ). In this paper we examine the potential for scientific investigation of Titan using JWST, primarily with three of the four instruments: NIRSpec, NIRCam and MIRI, noting that science with NIRISS will be complementary. Five core scientific themes are identified: (i) surface (ii) tropospheric clouds (iii) tropospheric gases (iv) stratospheric composition and (v) stratospheric hazes. We discuss each theme in depth, including the scientific purpose, capabilities and limitations of the instrument suite, and suggested observing schemes. We pay particular attention to saturation, which is a problem for all three instruments, but may be alleviated for NIRCam through use of selecting small sub-arrays of the detectors - sufficient to encompass Titan, but with significantly faster read-out times. We find that JWST has very significant potential for advancing Titan science, with a spectral resolution exceeding the Cassini instrument suite at near-infrared wavelengths, and a spatial resolution exceeding HST at the same wavelengths. In particular, JWST will be valuable for time-domain monitoring of Titan, given a five to ten year expected lifetime for the observatory, for example monitoring the seasonal appearance of clouds. JWST observations in the post-Cassini period will complement those of other large facilities such as HST, ALMA, SOFIA and next-generation ground-based telescopes (TMT, GMT, EELT).

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.

The spatial distribution of the tropospheric methane on Titan was measured using near-infrared spectroscopy. Ground-based observations at 1.5$\mu{\rm m}$ (H-band) were performed during the same night using instruments with adaptive optics at both the W. M. Keck Observatory and at the Paranal Observatory on 17 July 2014 UT. The integral field observations with SINFONI on the VLT covered the entire H-band at moderate resolving power, $R=\lambda/\Delta\lambda\approx1,500$, while the Keck observations were performed with NIRSPAO near 1.55254$\mu{\rm m}$ at higher resolution, $R\approx25,000$. The moderate resolution observations are used for flux calibration and for the determination of model parameters that can be degenerate in the interpretation of high resolution spectra. Line-by-line calculations of CH$_4$ and CH$_3$D correlated $k$ distributions from the HITRAN 2012 database were used, which incorporate revised line assignments near 1.5$\mu{\rm m}$. We fit the surface albedo and aerosol distributions in the VLT SINFONI observations that cover the entire H-band window and used these quantities to constrain the models of the high-resolution Keck NIRSPAO spectra when retrieving the methane abundances. Cassini VIMS images of the polar regions, acquired on 20 July 2014 UT, are used to validate the assumption that the opacity of tropospheric aerosol is relatively uniform below 10 km. We retrieved methane abundances at latitudes between 42$^{\circ}$ S and 80$^{\circ}$ N. The tropospheric methane in the Southern mid-latitudes was enhanced by a factor of $\sim$10--40% over the nominal profile that was measured using the GCMS on Huygens. The Northern hemisphere had $\sim$90% of the nominal methane abundance up to polar latitudes (80$^{\circ}$N). These measurements suggest that a source of saturated polar air is equilibrating with dryer conditions at lower latitudes.

We investigate the effects of varying Saturn's orbit on the atmospheric circulation and surface methane distribution of Titan. Using a new general circulation model of Titan's atmosphere, we simulate its climate under four characteristic configurations of orbital parameters that correspond to snapshots over the past 42 kyr, capturing the amplitude range of long-period cyclic variations in eccentricity and longitude of perihelion. The model, which covers pressures from the surface to 0.5 mbar, reproduces the present-day temperature profile and tropospheric superrotation. In all four simulations, the atmosphere efficiently transports methane poleward, drying out the low- and mid-latitudes, indicating that these regions have been desert-like for at least tens of thousands of years. Though circulation patterns are not significantly different, the amount of surface methane that builds up over either pole strongly depends on the insolation distribution; in the present-day, methane builds up preferentially in the north, in agreement with observations, where summer is milder but longer. The same is true, to a lesser extent, for the configuration 14 kyr ago, while the south pole gains more methane in the case for 28 kyr ago, and the system is almost symmetric 42 kyr ago. This confirms the hypothesis that orbital forcing influences the distribution of surface liquids, and that the current observed asymmetry could have been partially or fully reversed in the past. The evolution of the orbital forcing implies that the surface reservoir is transported on timescales of $\sim$30 kyr, in which case the asymmetry reverses with a period of $\sim$125 kyr. Otherwise, the orbital forcing does not produce a net asymmetry over longer timescales, and is not a likely mechanism for generating the observed dichotomy.

In this first paper of a series on the formation and abundance of substellar mass dwarf black holes (DBHs), we present a heuristic for deducing the stability of non-rotating matter embedded in a medium against collapse and the formation of a black hole. We demonstrate the heuristic's accuracy for a family of spherical mass distributions whose stability is known through other means. We also present the applications of this heuristic that could be applied to data sets of various types of simulations, including the possible formation of DBHs in the expanding gases of extreme astrophysical phenomena including Type Ia and Type II supernovae, hypernovae, and in the collision of two compact objects. These papers will also explore the observational and cosmological implications of DBHs, including estimates of the total masses of these objects bound in galaxies and ejected into the intergalactic medium. Applying our formalism to a Type II supernova simulation, we have found regions in one data set that are within a factor of three to four of both the density and mass necessary to create a DBH.

Fresh insights and powerful numerical tools are revitalizing the theoretical exploration of the supernova mechanism. The realization that the protoneutron star is Rayleigh-Taylor unstable at various times and radii and, hence, that a multi-dimensional perspective is required is one agent of this revolution. However, a new physical understanding of the nature of explosions (even spherical explosions) that are driven by neutrino heating and that escape from deep within a gravitational potential well is also emerging. This, together with the new multi-dimensional approach, promises to establish a paradigm within which supernova explosions and their consequences can be studied in the future.