Edwin S. Kite

Optical Images of an Exosolar Planet 25 Light Years from Earth

Fomalhaut, a bright star 7.7 parsecs (25 light-years) from Earth, harbors a belt of cold dust with a structure consistent with gravitational sculpting by an orbiting planet. Here, we present optical observations of an exoplanet candidate, Fomalhaut b. Fomalhaut b lies about 119 astronomical units (AU) from the star and 18 AU of the dust belt, matching predictions of its location. Hubble Space Telescope observations separated by 1.73 years reveal counterclockwise orbital motion. Dynamical models of the interaction between the planet and the belt indicate that the planets mass is at most three times that of Jupiter; a higher mass would lead to gravitational disruption of the belt, matching predictions of its location. The flux detected at 0.8 μm is also consistent with that of a planet with mass no greater than a few times that of Jupiter. The brightness at 0.6 μm and the lack of detection at longer wavelengths suggest that the detected flux may include starlight reflected off a circumplanetary disk, with dimension comparable to the orbits of the Galilean satellites. We also observe variability of unknown origin at 0.6 μm.

FOMALHAUT'S DEBRIS DISK AND PLANET: CONSTRAINING THE MASS OF FOMALHAUT B FROM DISK MORPHOLOGY

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhauts debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M pl 101.5 AU, and an orbital eccentricity e pl = 0.11-0.13. These conclusions are independent of Fom bs photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom bs orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a ≈ 133 AU lies at the periphery of Fom bs chaotic zone, and the mean disk eccentricity of e ≈ 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planets chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of ~ 100 Myr, and model them separately from their dust grain progeny; the latters orbits are strongly affected by radiation pressure and their lifetimes are limited to ~ 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut bs nominal space velocity does not bear this out, but the astrometric uncertainties may be large. If the apsidal misalignment proves real, our calculated upper mass limit of 3M J still holds. If the orbits are aligned, our model predicts M pl = 0.5M J, a pl = 115 AU, and e pl = 0.12. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M ⊕ of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Possible Disintegrating Short-period Super-Mercury Orbiting KIC 12557548

We report on the discovery of stellar occultations, observed with Kepler, which recur periodically at 15.685 hr intervals, but which vary in depth from a maximum of 1.3% to a minimum that can be less than 0.2%. The star that is apparently being occulted is KIC 12557548, a V = 16 mag K dwarf with T_eff, s ≃ 4400 K. The out-of-occultation behavior shows no evidence for ellipsoidal light variations, indicating that the mass of the orbiting object is less than ~3 M_J (for an orbital period of 15.7 hr). Because the eclipse depths are highly variable, they cannot be due solely to transits of a single planet with a fixed size. We discuss but dismiss a scenario involving a binary giant planet whose mutual orbit plane precesses, bringing one of the planets into and out of a grazing transit. This scenario seems ruled out by the dynamical instability that would result from such a configuration. We also briefly consider an eclipsing binary, possibly containing an accretion disk, that either orbits KIC 12557548 in a hierarchical triple configuration or is nearby on the sky, but we find such a scenario inadequate to reproduce the observations. The much more likely explanation—but one which still requires more quantitative development—involves macroscopic particles escaping the atmosphere of a slowly disintegrating planet not much larger than Mercury in size. The particles could take the form of micron-sized pyroxene or aluminum oxide dust grains. The planetary surface is hot enough to sublimate and create a high-Z atmosphere; this atmosphere may be loaded with dust via cloud condensation or explosive volcanism. Atmospheric gas escapes the planet via a Parker-type thermal wind, dragging dust grains with it. We infer a mass-loss rate from the observations of order 1 M_⊕ Gyr^(–1), with a dust-to-gas ratio possibly of order unity. For our fiducial 0.1 M_⊕ planet (twice the mass of Mercury), the evaporation timescale may be ~0.2 Gyr. Smaller mass planets are disfavored because they evaporate still more quickly, as are larger mass planets because they have surface gravities too strong to sustain outflows with the requisite mass-loss rates. The occultation profile evinces an ingress-egress asymmetry that could reflect a comet-like dust tail trailing the planet; we present simulations of such a tail.

GEODYNAMICS AND RATE OF VOLCANISM ON MASSIVE EARTH-LIKE PLANETS

We provide estimates of volcanism versus time for planets with Earth-like composition and masses 0.25–25 M⊕, 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. 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 10 Gyrold τ Ceti system, and tidal forcing can be shown to be weak, this would be evidence for plate tectonics.

Growth and form of the mound in Gale Crater, Mars: Slope wind enhanced erosion and transport

Ancient sediments provide archives of climate and habitability on Mars. Gale Crater, the landing site for the Mars Science Laboratory (MSL), hosts a 5-km-high sedimentary mound (Mount Sharp/Aeolis Mons). Hypotheses for mound formation include evaporitic, lacustrine, fluviodeltaic, and aeolian processes, but the origin and original extent of Gale’s mound is unknown. Here we show new measurements of sedimentary strata within the mound that indicate ∼3° outward dips oriented radially away from the mound center, inconsistent with the first three hypotheses. Moreover, although mounds are widely considered to be erosional remnants of a once crater-filling unit, we find that the Gale mound’s current form is close to its maximal extent. Instead we propose that the mound’s structure, stratigraphy, and current shape can be explained by growth in place near the center of the crater mediated by wind-topography feedbacks. Our model shows how sediment can initially accrete near the crater center far from crater-wall katabatic winds, until the increasing relief of the resulting mound generates mound-flank slope winds strong enough to erode the mound. The slope wind enhanced erosion and transport (SWEET) hypothesis indicates mound formation dominantly by aeolian deposition with limited organic carbon preservation potential, and a relatively limited role for lacustrine and fluvial activity. Morphodynamic feedbacks between wind and topography are widely applicable to a range of sedimentary and ice mounds across the Martian surface, and possibly other planets.

Seasonal melting and the formation of sedimentary rocks on Mars, with predictions for the Gale Crater mound

A model for the formation and distribution of sedimentary rocks on Mars is proposed. In this model (ISEE-Mars), 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 CO_2 atmosphere, dusty snow, and solar luminosity reduced by 23%. For these conditions snow melts only near the equator, when obliquity and eccentricity are high, and when perihelion occurs near equinox. These requirements for melting are satisfied by 0.01–20% of the probability distribution of Mars’ past spin–orbit parameters. This fraction is small, consistent with the geologic record of metastable surface liquid water acting as a “wet-pass filter” of Mars climate history, only recording orbital conditions that permitted surface liquid water. Total melt production is sufficient to account for observed aqueous alteration. 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 model Mars paleoclimate that only rarely produced enough meltwater to precipitate aqueous cements (sulfates, carbonates, phyllosilicates and silica) and indurate sediment. This is consistent with observations suggesting that surface aqueous alteration on Mars was brief and at low water/rock ratio. 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 Curiosity rover at Gale Crater’s mound (Mount Sharp, Aeolis Mons). Gale Crater’s mound is predicted to be a hemispheric maximum for snowmelt on Mars.

Low palaeopressure of the martian atmosphere estimated from the size distribution of ancient craters

Decay of the CO2-dominated atmosphere is an important component of long-term environmental change on Mars, but direct constraints on paleoatmospheric pressure P are few. Of particular interest is the climate that allowed rivers to flow early in Mars history, which was affected by P via direct and indirect greenhouse effects. The size of craters embedded within ancient layered sediments is a proxy for P: the smaller the minimum-sized craters that form, the thinner the past atmosphere. Here we use high-resolution orthophotos and Digital Terrain Models (DTMs) to identify ancient craters among the river deposits of Aeolis close to Gale crater, and compare their sizes to models of atmospheric filtering of impactors by thicker atmospheres. We obtain an upper limit of P <= 760+/-70 mbar, rising to P <= 1640+/-180 mbar if rimmed circular mesas are excluded. Our work assumes target properties appropriate for desert alluvium: if sediment developed bedrock-like rock-mass strength by early diagenesis, the upper limit increases by a factor of up to 2. If Mars did not have a stable multibar atmosphere at the time that the rivers were flowing, the warm-wet CO2 greenhouse of Pollack et al. (1987) is ruled out, and long-term average temperatures were most likely below freezing.The decay of the martian atmosphere—which is dominated by carbon dioxide—is a component of the long-term environmental change on Mars from a climate that once allowed rivers to flow to the cold and dry conditions of today. The minimum size of craters serves as a proxy for palaeopressure of planetary atmospheres, because thinner atmospheres permit smaller objects to reach the surface at high velocities and form craters. The Aeolis Dorsa region near Gale crater on Mars contains a high density of preserved ancient craters interbedded with river deposits and thus can provide constraints on atmospheric density at the time of fluvial activity. Here we use high-resolution images and digital terrain models from the Mars Reconnaissance Orbiter to identify ancient craters in deposits in Aeolis Dorsa that date to about 3.6 Gyr ago and compare their size distribution with models of atmospheric filtering of impactors. We obtain an upper limit of 0.9 ± 0.1 bar for the martian atmospheric palaeopressure, rising to 1.9 ± 0.2 bar if rimmed circular mesas—interpreted to be erosionally-resistant fills or floors of impact craters—are excluded. We assume target properties appropriate for desert alluvium: if sediment had rock-mass strength similar to bedrock at the time of impact, the paleopressure upper limit increases by a factor of up to two. If Mars did not have a stable multibar atmosphere at the time that the rivers were flowing—as suggested by our results—then a warm and wet CO_2/H_2O greenhouse is ruled out, and long-term average temperatures were most likely below freezing.

Geophysical and geochemical constraints on geoneutrino fluxes from Earth's mantle

Knowledge of the amount and distribution of radiogenic heating in the mantle is crucial for understanding the dynamics of the Earth, including its thermal evolution, the style and planform of mantle convection, and the energetics of the core. Although the flux of heat from the surface of the planet is robustly estimated, the contributions of radiogenic heating and secular cooling remain poorly defined. Constraining the amount of heat-producing elements in the Earth will provide clues to understanding nebula condensation and planetary formation processes in early Solar System. Mantle radioactivity supplies power for mantle convection and plate tectonics, but estimates of mantle radiogenic heat production vary by a factor of more than 20. Recent experimental results demonstrate the potential for direct assessment of mantle radioactivity through observations of geoneutrinos, which are emitted by naturally occurring radionuclides. Predictions of the geoneutrino signal from the mantle exist for several established estimates of mantle composition. Here we present novel analyses, illustrating surface variations of the mantle geoneutrino signal for models of the deep mantle structure, including those based on seismic tomography. These variations have measurable differences for some models, allowing new and meaningful constraints on the dynamics of the planet. An ocean based geoneutrino detector deployed at several strategic locations will be able to discriminate between competing compositional models of the bulk silicate Earth.

CLIMATE INSTABILITY ON TIDALLY LOCKED EXOPLANETS

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 planets maintain stable surface liquid water over geologic 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 habitable-zone rocky planets, are expected to be in synchronous rotation (tidally locked). In this paper, we investigate two hypothetical feedbacks that can destabilize climate on planets in synchronous rotation. (1) If small changes in pressure alter the temperature distribution across a planets surface such that the weathering rate goes up when the pressure goes down, a runaway positive feedback occurs involving increasing weathering rate near the substellar point, decreasing pressure, and increasing substellar surface temperature. We call this feedback enhanced substellar weathering instability (ESWI). (2) When decreases in pressure increase the fraction of 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 will occur. This substellar dissolution feedback can also cause a runaway climate shift. We use an idealized energy balance model to map out the conditions under which these instabilities may occur. In this simplified model, the weathering runaway can shrink the habitable zone and cause geologically rapid 103-fold atmospheric pressure shifts within the habitable zone. Mars may have undergone a weathering runaway in the past. Substellar dissolution is usually a negative feedback or weak positive feedback on changes in atmospheric pressure. It can only cause runaway changes for small, deep oceans and highly soluble atmospheric gases. Both instabilities are suppressed if the atmosphere has a high radiative efficiency. Our results are most relevant for atmospheres that are thin, have low greenhouse-gas radiative efficiency, and have a principal greenhouse gas that is also the main constituent of the atmosphere. ESWI also requires land near the substellar point, and tectonic resurfacing (volcanism, mountain-building) is needed for large jumps in pressure. These results identify a new pathway by which habitable-zone planets can undergo rapid climate shifts and become uninhabitable.

Pacing early Mars river activity: Embedded craters in the Aeolis Dorsa region imply river activity spanned ≳(1–20) Myr

The impactor flux early in Mars history was much higher than today, so sedimentary sequences include many buried craters. In combination with models for the impactor flux, observations of the number of buried craters can constrain sedimentation rates. Using the frequency of crater-river interactions, we find net sedimentation rate \lesssim 20-300 {\mu}m/yr at Aeolis Dorsa. This sets a lower bound of 1-15 Myr on the total interval spanned by fluvial activity around the Noachian-Hesperian transition. We predict that Gale Craters mound (Aeolis Mons) took at least 10-100 Myr to accumulate, which is testable by the Mars Science Laboratory.