Active eruptions from the south polar region of Saturn's small (~500 km diameter) moon Enceladus are concentrated along a series of lineaments known as the `tiger stripes', thought to be partially open fissures that connect to the liquid water ocean beneath the ice shell. Whereas aspects of the tiger stripes have been addressed in previous work, no study to date simultaneously explains why they should be located only at the south pole, why there are multiple approximately parallel and regularly spaced fractures, and what accounts for their spacing of ~35 km. Here we propose that secular cooling and the resulting ice shell thickening and global tensile stresses cause the first fracture to form at one of the poles, where the ice shell is thinnest due to tidal heating. The tensile stresses are thereby partially relieved, preventing a similar failure at the opposite pole. We propose that subsequent activity then concentrates in the vicinity of the first fracture as the steadily erupted water ice loads the flanks of the open fissure, causing bending in the surrounding elastic plate and further tensile failure in bands parallel to the first fracture, leading to a cascading sequence of parallel fissures until the conditions no longer permit through-going fractures.

Tidal interactions between planets or stars and the bodies that orbit them dissipate energy in their interiors. The energy dissipated drives internal heating and a fraction of that energy will be released as seismic energy. Here we formalize a model to describe the tidally-driven seismic activity on planetary bodies based on tidal dissipation.

A model for the formation and distribution of sedimentary rocks on Mars is proposed. The rate-limiting step is supply of liquid water from seasonal melting of snow or ice. The model is run for a O(10^2) mbar pure CO2 atmosphere, dusty snow, and solar luminosity reduced by 23%. For these conditions snow only melts near the equator, and only when obliquity >40 degrees, eccentricity >0.12, and perihelion occurs near equinox. These requirements for melting are satisfied by 0.01-20% of the probability distribution of Mars' past spin-orbit parameters. Total melt production is sufficient to account for aqueous alteration of the sedimentary rocks. The pattern of seasonal snowmelt is integrated over all spin-orbit parameters and compared to the observed distribution of sedimentary rocks. The global distribution of snowmelt has maxima in Valles Marineris, Meridiani Planum and Gale Crater. These correspond to maxima in the sedimentary-rock distribution. Higher pressures and especially higher temperatures lead to melting over a broader range of spin-orbit parameters. The pattern of sedimentary rocks on Mars is most consistent with a Mars paleoclimate that only rarely produced enough meltwater to precipitate aqueous cements and indurate sediment. The results suggest intermittency of snowmelt and long globally-dry intervals, unfavorable for past life on Mars. This model makes testable predictions for the Mars Science Laboratory rover at Gale Crater. Gale Crater is predicted to be a hemispheric maximum for snowmelt on Mars.

Feedbacks that can destabilize the climates of synchronously-rotating rocky planets may arise on planets with strong day-night surface temperature contrasts. Earth-like habitable-zone (HZ) planets maintain stable surface liquid water over geological time. This requires equilibrium between the temperature-dependent rate of greenhouse-gas consumption by weathering,and greenhouse-gas resupply by other processes. Detected small-radius exoplanets, and anticipated M-dwarf HZ rocky planets, are expected to be tidally locked. We investigate two feedbacks that can destabilize climate on tidally-locked planets. (1) If small changes in pressure alter the temperature distribution across a planet's surface such that the weathering rate increases when the pressure decreases, a positive feedback occurs involving increasing weathering rate near the substellar point, decreasing pressure, and increasing substellar surface temperature. (2) When decreases in pressure increase the surface area above the melting point (through reduced advective cooling of the substellar point), and the corresponding increase in volume of liquid causes net dissolution of the atmosphere, a further decrease in pressure occurs. We use an idealized energy balance model to map out the conditions under which these instabilities may occur. The weathering runaway can shrink the habitable zone, and cause geologically rapid 10^3-fold pressure shifts within the HZ. Mars may have undergone a weathering runaway in the past. Substellar dissolution is usually a negative feedback or weak positive feedback on changes in pressure. Both instabilities are suppressed if the atmosphere has a high radiative efficiency. Our results are most relevant for atmospheres that are thin and have low greenhouse-gas radiative efficiency. These results identify a new pathway by which HZ planets can undergo rapid climate shifts and become uninhabitable.

Channel networks on the plateau adjacent to Juventae Chasma have the highest drainage densities reported on Mars.We model frozen precipitation on the Juventae plateau,finding that the trigger for forming these channel networks could have been ephemeral lakeshore precipitation,and that they do not require past temperatures higher than today.If short-lived and localized events explain some dendritic channel networks on Mars, this would weaken the link between dendritic valley networks and surface climate conditions that could sustain life. Our analysis uses MRAMS simulations and HiRISE DTMs.We model localized weather systems driven by water vapor release from ephemeral lakes during outflow channel formation.At Juventae Chasma,mean snowfall reaches a maximum of 0.9mm/hr water equivalent on the SW rim of the chasm.Radiative effects of the thick cloud cover raise maximum (minimum, mean) plateau surface temperatures by up to 24K(9K, 17K)locally.The key result is that the area of maximum modeled precipitation shows a striking correspondence to the mapped Juventae plateau channel networks.Three independent methods show this fit is unlikely to be due to chance.We use a snowpack energy balance model to show that if the snow has the albedo of dust(0.28), and for a solar luminosity of 0.8($\equiv$3.0Gya), then if the atmospheric greenhouse effect is unchanged from(6K warmer than)today only 0.4%(21%)of lake-induced precipitation events produce snowpack that undergoes melting.However, warming from associated dense cloud cover would allow melting over a wider range of conditions.In these localized precipitation scenarios, global temperatures need not be higher than today, and the rest of the planet remains dry.

We use the Mars Regional Atmospheric Modeling System (MRAMS) to simulate lake storms on Mars, finding that intense localized precipitation will occur for lake size >=10^3 km^2. Mars has a low-density atmosphere, so deep convection can be triggered by small amounts of latent heat release. In our reference simulation, the buoyant plume lifts vapor above condensation level, forming a 20km-high optically-thick cloud. Ice grains grow to 200 microns radius and fall near (or in) the lake at mean rates up to 1.5 mm/hr water equivalent (maximum rates up to 6 mm/hr water equivalent). Because atmospheric temperatures outside the surface layer are always well below 273K, supersaturation and condensation begin at low altitudes above lakes on Mars. In contrast to Earth lake-effect storms, lake storms on Mars involve continuous precipitation, and their vertical velocities and plume heights exceed those of tropical thunderstorms on Earth. Convection does not reach above the planetary boundary layer for lakes <<10^3 km^2 or for atmospheric pressure >O(10^2) mbar. Instead, vapor is advected downwind with little cloud formation. Precipitation occurs as snow, and the daytime radiative forcing at the land surface due to plume vapor and storm clouds is too small to melt snow directly (<+10 W/m^2). However, if orbital conditions are favorable, then the snow may be seasonally unstable to melting and produce runoff to form channels. We calculate the probability of melting by running thermal models over all possible orbital conditions and weighting their outcomes by probabilities given by Laskar et al., 2004. We determine that for an equatorial vapor source, sunlight 15% fainter than at present, and snowpack with albedo 0.28 (0.35), melting may occur with 4%(0.1%) probability. This rises to 56%(12%) if the ancient greenhouse effect was modestly (6K) greater than today.

To ascertain whether magnetic dynamos operate in rocky exoplanets more massive or hotter than the Earth, we developed a parametric model of a differentiated rocky planet and its thermal evolution. Our model reproduces the established properties of Earth's interior and magnetic field at the present time. When applied to Venus, assuming that planet lacks plate tectonics and has a dehydrated mantle with an elevated viscosity, the model shows that the dynamo shuts down or never operated. Our model predicts that at a fixed planet mass, dynamo history is sensitive to core size, but not to the initial inventory of long-lived, heat-producing radionuclides. It predicts that rocky planets larger than 2.5 Earth masses will not develop inner cores because the temperature-pressure slope of the iron solidus becomes flatter than that of the core adiabat. Instead, iron "snow" will condense near or at the top of these cores, and the net transfer of latent heat upwards will suppress convection and a dynamo. More massive planets can have anemic dynamos due to core cooling, but only if they have mobile lids (plate tectonics). The lifetime of these dynamos is shorter with increasing planet mass but longer with higher surface temperature. Massive Venus-like planets with stagnant lids and more viscous mantles will lack dynamos altogether. We identify two alternative sources of magnetic fields on rocky planets: eddy currents induced in the hot or molten upper layers of planets on very short period orbits, and dynamos in the ionic conducting layers of "ocean" planets with ~10% mass in an upper mantle of water (ice).

We provide estimates of volcanism versus time for planets with Earth-like composition and masses from 0.25 to 25 times Earth, as a step toward predicting atmospheric mass on extrasolar rocky planets. Volcanism requires melting of the silicate mantle. We use a thermal evolution model, calibrated against Earth, in combination with standard melting models, to explore the dependence of convection-driven decompression mantle melting on planet mass. Here we show that (1) volcanism is likely to proceed on massive planets with plate tectonics over the main-sequence lifetime of the parent star; (2) crustal thickness (and melting rate normalized to planet mass) is weakly dependent on planet mass; (3) stagnant lid planets live fast (they have higher rates of melting than their plate tectonic counterparts early in their thermal evolution) but die young (melting shuts down after a few Gyr); (4) plate tectonics may not operate on high mass planets because of the production of buoyant crust which is difficult to subduct; and (5) melting is necessary but insufficient for efficient volcanic degassing - volatiles partition into the earliest, deepest melts, which may be denser than the residue and sink to the base of the mantle on young, massive planets. Magma must also crystallize at or near the surface, and the pressure of overlying volatiles must be fairly low, if volatiles are to reach the surface. If volcanism is detected in the Tau Ceti system, and tidal forcing can be shown to be weak, this would be evidence for plate tectonics.

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.