The TRAPPIST-1 system is sufficiently closely packed that tides raised by one planet on another are significant. We investigate whether this source of tidal heating is comparable to eccentricity tides raised by the star. Assuming a homogeneous body with a Maxwell rheology, we find that energy dissipation from stellar tides always dominates over that from planet-planet tides across a range of viscosities. TRAPPIST-1 g may experience the greatest proportion of planet-planet tidal heating, where it can account for between 2% and 20% of the total amount of tidal heating, for high ($10^{21}$ Pa s) and low viscosity ($10^{14}$ Pa s) regimes, respectively. If planet-planet tidal heating is to exceed that from stellar eccentricity tides, orbital eccentricities must be no more than e=$10^{-3}$ to $10^{-4}$ for most of the TRAPPIST-1 planets.

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

Isostatic equilibrium is commonly defined as the state achieved when there are no lateral gradients in hydrostatic pressure, and thus no lateral flow, at depth within the lower viscosity mantle that underlies a planetary body's outer crust. In a constant-gravity Cartesian framework, this definition is equivalent to the requirement that columns of equal width contain equal masses. Here we show, however, that this equivalence breaks down when the spherical geometry of the problem is taken into account. Imposing the "equal masses" requirement in a spherical geometry, as is commonly done in the literature, leads to significant lateral pressure gradients along internal equipotential surfaces, and thus corresponds to a state of disequilibrium. Compared with the "equal pressures" model we present here, the "equal masses" model always overestimates the compensation depth--by ~27% in the case of the lunar highlands and by nearly a factor of two in the case of Enceladus.

As first noted 200 years ago by Laplace, the Moon's rotational and tidal bulges are significantly larger than expected, given the Moon's present orbital and rotational state. This excess deformation has been ascribed to a fossil figure, frozen in when the Moon was closer to the Earth. However, the observed figure is only consistent with an eccentric and non-synchronous orbit, contrary to our understanding of the Moon's formation and evolution. Here, we show that lunar mascons and impact basins have a significant contribution to the observed lunar figure. Removing their contribution reveals a misaligned fossil figure consistent with an early epoch of true polar wander (driven by the formation of the South Pole-Aitken Basin) and an early low-eccentricity, synchronous lunar orbit. This new self-consistent model solves a long-standing problem in planetary science, and will inform future studies of the Moon's dynamical evolution and early dynamo.

We present a model for the dispersal of protoplanetary disks by winds from either the central star or the inner disk. These winds obliquely strike the flaring disk surface and strip away disk material by entraining it in an outward radial-moving flow at the wind-disk interface which lies several disk scale heights above the mid-plane. The disk dispersal time depends on the entrainment velocity at which disk material flows into this turbulent shear layer interface. If the entrainment efficiency is ~10% of the local sound speed, a likely upper limit, the dispersal time at 1 AU is ~6 Myr for a disk with a surface density of 10^3 g cm^{-2}, a solar mass central star, and a wind with an outflow rate 10^{-8} Msun/yr and terminal velocity 200 km/s. When compared to photoevaporation and viscous evolution, wind stripping can be a dominant mechanism only for the combination of low accretion rates (< 10^{-8} Msun/yr) and wind outflow rates approaching these accretion rates. This case is unusual since generally outflow rates are < 0.1 of of accretion rates.

The areal centroids of the youngest polar deposits on Mars are offset from those of adjacent paleopolar deposits by 5-10 degrees. We test the hypothesis that the offset is the result of true polar wander (TPW), the motion of the solid surface with respect to the spin axis, caused by a mass redistribution within or on the surface of Mars. In particular, we consider TPW driven by late-stage volcanism during the late Hesperian to Amazonian. There is observational and qualitative support for this hypothesis: in both North and South, observed offsets lie close to a great circle 90 degrees from Tharsis, as expected for polar wander after Tharsis formed. We calculate the magnitude and direction of TPW produced by mapped late-stage lavas for a range of lithospheric thicknesses, lava thicknesses, eruption histories, and prior polar wander events. If Tharsis formed close to the equator, the stabilizing effect of a fossil rotational bulge located close to the equator leads to predicted TPW of <2 degrees, too small to account for observed offsets. If, however, Tharsis formed far from the equator, late-stage TPW driven by low-latitude, late-stage volcanism would be 6-33 degrees, similar to that inferred from the location of paleopolar deposits. 4.4+/-1.3x10^19 kg of young erupted lava can account for the offset of the Dorsa Argentea Formation from the present-day south rotation pole. This mass is consistent with prior mapping-based estimates and would imply a mass release of CO2 by volcanic degassing similar to that in the atmosphere at the present time. The South Polar Layered Deposits are offset from the spin axis in the opposite sense to the other paleopolar deposits. This can be explained by an additional contribution from a plume beneath Elysium. We conclude with a list of observational tests of the TPW hypothesis.

Photoevaporation may provide an explanation for the short lifetimes of disks around young stars. With the exception of neutral oxygen lines, the observed low-velocity forbidden line emission from T Tauri stars can be reproduced by photoevaporating models. The natural formation of a gap in the disk at several AU due to photoevaporation and viscous spreading provides a possible halting mechanism for migrating planets and an explanation for the abundance of observed planets at these radii.

We report detection of continuum emission at 850 and 450 micron from disks around four Classical T Tauri stars in the MBM 12 (L1457) young association. Using a simple model we infer masses of 0.0014-0.012 M_sun for the disk of LkHa 263 ABC, 0.005-0.021 M_sun for S18 ABab, 0.03-0.18 M_sun for LkHa 264 A, and 0.023-0.23 M_sun for LkHa 262. The disk mass found for LkHa 263 ABC is consistent with the 0.0018 M_sun inferred from the scattered light image of the edge-on disk around component C. Comparison to earlier 13CO line observations indicates CO depletion by up to a factor 300 with respect to dark-cloud values. The spectral energy distributions (SED) suggest grain growth, possibly to sizes of a few hundred micron, but our spatially unresolved data cannot rule out opacity as an explanation for the SED shape. Our observations show that these T Tauri stars are still surrounded by significant reservoirs of cold material at an age of 1-5 Myr. We conclude that the observed differences in disk mass are likely explained by binary separation affecting the initial value. With available accretion rate estimates we find that our data are consistent with theoretical expectations for viscously evolving disks having decreased their masses by ~30%.

The recent discovery of Jupiter-mass planets orbiting at a few AU from their stars compliments earlier detections of massive planets on very small orbits. The short period orbits strongly suggest that planet migration has occurred, with the likely mechanism being tidal interactions between the planets and the gas disks out of which they formed. The newly discovered long period planets, together with the gas giant planets in our solar system, show that migration is either absent or rapidly halted in at least some systems. We propose a mechanism for halting type-II migration at several AU in a gas disk. Photoevaporation of the disk by irradiation from the central star can produce a gap in the disk at a few AU, preventing planets outside the gap from migrating down to the star. This would result in an excess of systems with planets at or just outside the photoevaporation radius.

The evolution of a stellar disk under the influence of viscous evolution, photoevaporation from the central source, and photoevaporation by external stars is studied. We take the typical parameters of TTSs and the Trapezium Cluster conditions. The photoionizing flux from the central source is assumed to arise both from the quiescent star and accretion shocks at the base of stellar magnetospheric columns, along which material from the disk accretes. The accretion flux is calculated self-consistently from the accretion mass loss rate. We find that the disk cannot be entirely removed using only viscous evolution and photoionization from the disk-star accretion shock. However, when FUV photoevaporation by external massive stars is included the disk is removed in 10^6 -10^7yr; and when EUV photoevaporation by external massive stars is included the disk is removed in 10^5 - 10^6yr. An intriguing feature of photoevaporation by the central star is the formation of a gap in the disk at late stages of the disk evolution. As the gap starts forming, viscous spreading and photoevaporation work in resonance. There is no gap formation for disks nearby external massive stars because the outer annuli are quickly removed by the dominant EUV flux. On the other hand, at larger, more typical distances (d>>0.03pc) from the external stars the flux is FUV dominated. As a consequence, the disk is efficiently evaporated at two different locations; forming a gap during the last stages of the disk evolution.