Source author record

Jonathan I. Lunine

Jonathan I. Lunine appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Jupiter's Temperature Structure: A Reassessment of the Voyager Radio Occultation Measurements

The thermal structure of planetary atmospheres is an essential input for predicting and retrieving the distribution of gases and aerosols, as well as the bulk chemical abundances. In the case of Jupiter, the temperature at a reference level - generally taken at 1 bar - serves as the anchor in models used to derive the planet's interior structure and composition. Most models assume the temperature measured by the Galileo probe (Seiff et al. 1998). However, those data correspond to a single location, an unusually clear, dry region, affected by local atmospheric dynamics. On the other hand, the Voyager radio occultation observations cover a wider range of latitudes, longitudes, and times (Lindal et al. 1981). The Voyager retrievals were based on atmospheric composition and radio refractivity data that require updating and were never properly tabulated: the few existing tabulations are incomplete and ambiguous. Here, we present a systematic electronic digitization of all available temperature profiles from Voyager, followed by their reanalysis, employing currently accepted values of the abundances and radio refractivities of atmospheric species. We find the corrected temperature at the 1 bar level to be up to 4 K greater than previously published values, i.e., 170.3{\pm}3.8 K at 12°S (Voyager 1 ingress) and 167.3{\pm}3.8 K at 0°N (Voyager 1 egress). This is to be compared with the Galileo probe value of 166.1{\pm}0.8 K at the edge of an unusual feature at 6.57°N. Altogether, this suggests that Jupiter's tropospheric temperatures may vary spatially by up to 7 K between 7°N and 12°S.

preprint2022arXiv

Revelations on Jupiter's Formation, Evolution and Interior: Challenges from Juno Results

The Juno mission has revolutionized and challenged our understanding of Jupiter. As Juno transitioned to its extended mission, we review the major findings of Jupiter's internal structure relevant to understanding Jupiter's formation and evolution. Results from Juno's investigation of Jupiter's interior structure imply that the planet has compositional gradients and is accordingly non-adiabatic, with a complex internal structure. These new results imply that current models of Jupiter's formation and evolution require a revision. In this paper, we discuss potential formation and evolution paths that can lead to an internal structure model consistent with Juno data, and the constraints they provide. We note that standard core accretion formation models, including the heavy-element enrichment during planetary growth is consistent with an interior that is inhomogeneous with composition gradients in its deep interior. However, such formation models typically predict that this region, which could be interpreted as a primordial dilute core, is confined to about 10% of Jupiter's total mass. In contrast, structure models that fit Juno data imply that this region contains 30% of the mass or more. One way to explain the origin of this extended region is by invoking a relatively long (about 2 Myrs) formation phase where the growing planet accretes gas and planetesimals delaying the runaway gas accretion. Alternatively, Jupiter's fuzzy core could be a result of a giant impact or convection post-formation. These novel scenarios require somewhat special and specific conditions. Clarity on the plausibility of such conditions could come from future high-resolution observations of planet-forming regions around other stars, from the observed and modeled architectures of extrasolar systems with giant planets, and future Juno data obtained during its extended mission.

preprint2021arXiv

Cold traps of hypervolatiles in the protosolar nebula at the origin of comet C/2016 R2 (PanSTARRS)'s peculiar composition

Recent observations of the long period comet C/2016 R2 (PanSTARRS) indicate an unusually high N2/CO abundance ratio, typically larger than 0.05, and at least 2-3 times higher than the one measured in 67P/Churyumov-Gerasimenko. Another striking compositional feature of this comet is its heavy depletion in H2O, compared to other comets. Here, we investigate the formation circumstances of a generic comet whose composition reproduces these two key features. We first envisage the possibility that this comet agglomerated from clathrates, but we find that such a scenario does not explain the observed low water abundance. We then alternatively investigate the possibility that the building blocks of the comet C/2016 R2 (PanSTARRS) agglomerated from grains and pebbles made of pure condensates via the use of a disk model describing the radial transport of volatiles. We show that N2/CO ratios reproducing the value estimated in this comet can be found in grains condensed in the vicinity of the CO and N2 icelines. Moreover, high CO/H2O ratios (>100 times the initial gas phase value) can be found in grains condensed in the vicinity of the CO iceline. If the building blocks of a comet assembled from such grains, they should present N2/CO and CO/H2O ratios consistent with the measurements made in comet C/2016 R2 (PanSTARRS)'s coma. Our scenario indicates that comet C/2016 R2 (PanSTARRS) formed in a colder environment than the other comets that share more usual compositions. Our model also explains the unusual composition of the interstellar comet 2l/Borisov.

preprint2021arXiv

Evidence for multiple Ferrel-like cells on Jupiter

Jupiter's atmosphere is dominated by multiple jet streams which are strongly tied to its 3D atmospheric circulation. Lacking a rigid bottom boundary, several models exist for how the meridional circulation extends into the planetary interior. Here we show, collecting evidence from multiple instruments of the Juno mission, the existence of mid-latitudinal meridional circulation cells which are driven by turbulence, similar to the Ferrel cells on Earth. Different than Earth, which contains only one such cell in each hemisphere, the larger, faster rotating Jupiter can incorporate multiple cells. The cells form regions of upwelling and downwelling, which we show are clearly evident in Juno's microwave data between latitude 60S and 60N. The existence of these cells is confirmed by reproducing the ammonia observations using a simplistic model. This study solves a long-standing puzzle regarding the nature of Jupiter's sub-cloud dynamics and provides evidence for 8 cells in each Jovian hemisphere.

preprint2020arXiv

The Saturn Ring Skimmer Mission Concept: The next step to explore Saturn's rings, atmosphere, interior, and inner magnetosphere

The innovative Saturn Ring Skimmer mission concept enables a wide range of investigations that address fundamental questions about Saturn and its rings, as well as giant planets and astrophysical disk systems in general. This mission would provide new insights into the dynamical processes that operate in astrophysical disk systems by observing individual particles in Saturn's rings for the first time. The Ring Skimmer would also constrain the origin, history, and fate of Saturn's rings by determining their compositional evolution and material transport rates. In addition, the Ring Skimmer would reveal how the rings, magnetosphere, and planet operate as an inter-connected system by making direct measurements of the ring's atmosphere, Saturn's inner magnetosphere and the material owing from the rings into the planet. At the same time, this mission would clarify the dynamical processes operating in the planet's visible atmosphere and deep interior by making extensive high-resolution observations of cloud features and repeated measurements of the planet's extremely dynamic gravitational field. Given the scientific potential of this basic mission concept, we advocate that it be studied in depth as a potential option for the New Frontiers program.

preprint2015arXiv

Nebular dead zone effects on the D/H ratio in chondrites and comets

Comets and chondrites show non-monotonic behaviour of their Deuterium to Hydrogen (D/H) ratio as a function of their formation location from the Sun. This is difficult to explain with a classical protoplanetary disk model that has a decreasing temperature structure with radius from the Sun. We want to understand if a protoplanetary disc with a dead zone, a region of zero or low turbulence, can explain the measured D/H values in comets and chondrites. We use time snapshots of a vertically layered disk model with turbulent surface layers and a dead zone at the midplane. The disc has a non-monotonic temperature structure due to increased heating from self-gravity in the outer parts of the dead zone. We couple this to a D/H ratio evolution model in order to quantify the effect of such thermal profiles on D/H enrichment in the nebula. We find that the local temperature peak in the disk can explain the diversity in the D/H ratios of different chondritic families. This disk temperature profile leads to a non-monotonic D/H enrichment evolution, allowing these families to acquire their different D/H values while forming in close proximity. The formation order we infer for these families is compatible with that inferred from their water abundances. However, we find that even for very young disks, the thermal profile reversal is too close to the Sun to be relevant for comets.

preprint2015arXiv

Orbital instability of close-in exomoons in non-coplanar systems

This work shows the dynamical instability that can happen to close-in satellites when planet oblateness is not accounted for in non-coplanar multiplanet systems. Simulations include two secularly interacting Jupiter-mass planets mutually inclined by 10 degrees, with the host planet either oblate or spherical. With a spherical host planet, moons within a critical planetocentric distance experience high inclinations and in some cases high eccentricities, while more distant moons orbit stably with low inclinations and eccentricities, as expected. These counter-intuitive dynamical phenomena disappear with an oblate host planet, in which case the moons' Laplace plane transitions from the host planet's equatorial plane to the host planet's precessing orbital plane as their semi-major axes increase, and all moons are dynamically stable with very mild changes in orbits. Direct perturbation from the perturbing planet has been investigated and ruled out as an explanation for the behavior of the innermost satellites, therefore leaving the central star's perturbation as the cause. Instability occurs while the nodal precession of the satellite and the central star (as seen from the host planet's frame) approaches the 1:1 secular resonance. In non-coplanar systems, around a non-oblate planet, the nodal precession of the moon becomes slow and comparable to that of the planet, giving rise to resonant configurations. The above effect needs to be taken into account in setting up numerical simulations.

preprint2015arXiv

The fate of ethane in Titan's hydrocarbon lakes and seas

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.

preprint2014arXiv

Carbon-rich planet formation in a solar composition disk

The C--to--O ratio is a crucial determinant of the chemical properties of planets. The recent observation of WASP 12b, a giant planet with a C/O value larger than that estimated for its host star, poses a conundrum for understanding the origin of this elemental ratio in any given planetary system. In this paper, we propose a mechanism for enhancing the value of C/O in the disk through the transport and distribution of volatiles. We construct a model that computes the abundances of major C and O bearing volatiles under the influence of gas drag, sublimation, vapor diffusion, condensation and coagulation in a multi--iceline 1+1D protoplanetary disk. We find a gradual depletion in water and carbon monoxide vapors inside the water's iceline with carbon monoxide depleting slower than water. This effect increases the gaseous C/O and decreases the C/H ratio in this region to values similar to those found in WASP 12b's day side atmosphere. Giant planets whose envelopes were accreted inside the water's iceline should then display C/O values larger than those of their parent stars, making them members of the class of so-called ``carbon-rich planets''.

preprint2014arXiv

Equilibrium composition between liquid and clathrate reservoirs on Titan

Hundreds of lakes and a few seas of liquid hydrocarbons have been observed by the Cassini spacecraft to cover the polar regions of Titan. A significant fraction of these lakes or seas could possibly be interconnected with subsurface liquid reservoirs of alkanes. In this paper, we investigate the interplay that would happen between a reservoir of liquid hydrocarbons located in Titan's subsurface and a hypothetical clathrate reservoir that progressively forms if the liquid mixture diffuses throughout a preexisting porous icy layer. To do so, we use a statistical-thermodynamic model in order to compute the composition of the clathrate reservoir that forms as a result of the progressive entrapping of the liquid mixture. This study shows that clathrate formation strongly fractionates the molecules between the liquid and the solid phases. Depending on whether the structure I or structure II clathrate forms, the present model predicts that the liquid reservoirs would be mainly composed of either propane or ethane, respectively. The other molecules present in the liquid are trapped in clathrates. Any river or lake emanating from subsurface liquid reservoirs that significantly interacted with clathrate reservoirs should present such composition. On the other hand, lakes and rivers sourced by precipitation should contain higher fractions of methane and nitrogen, as well as minor traces of argon and carbon monoxide.

preprint2014arXiv

GCM simulations of Titan's middle and lower atmosphere and comparison to observations

Simulation results are presented from a new general circulation model (GCM) of Titan, the Titan Atmospheric Model (TAM), which couples the Flexible Modeling System (FMS) spectral dynamical core to a suite of external/sub-grid-scale physics. These include a new non-gray radiative transfer module that takes advantage of recent data from Cassini-Huygens, large-scale condensation and quasi-equilibrium moist convection schemes, a surface model with "bucket" hydrology, and boundary layer turbulent diffusion. The model produces a realistic temperature structure from the surface to the lower mesosphere, including a stratopause, as well as satisfactory superrotation. The latter is shown to depend on the dynamical core's ability to build up angular momentum from surface torques. Simulated latitudinal temperature contrasts are adequate, compared to observations, and polar temperature anomalies agree with observations. In the lower atmosphere, the insolation distribution is shown to strongly impact turbulent fluxes, and surface heating is maximum at mid-latitudes. Surface liquids are unstable at mid- and low-latitudes, and quickly migrate poleward. The simulated humidity profile and distribution of surface temperatures, compared to observations, corroborate the prevalence of dry conditions at low latitudes. Polar cloud activity is well represented, though the observed mid-latitude clouds remain somewhat puzzling, and some formation alternatives are suggested.

preprint2014arXiv

New insights on Saturn's formation from its nitrogen isotopic composition

The recent derivation of a lower limit for the $^{14}$N/$^{15}$N ratio in Saturn's ammonia, which is found to be consistent with the Jovian value, prompted us to revise models of Saturn's formation using as constraints the supersolar abundances of heavy elements measured in its atmosphere. Here we find that it is possible to account for both Saturn's chemical and isotopic compositions if one assumes the formation of its building blocks at $\sim$45 K in the protosolar nebula, provided that the O abundance was $\sim$2.6 times protosolar in its feeding zone. To do so, we used a statistical thermodynamic model to investigate the composition of the clathrate phase that formed during the cooling of the protosolar nebula and from which the building blocks of Saturn were agglomerated. We find that Saturn's O/H is at least $\sim$34.9 times protosolar and that the corresponding mass of heavy elements ($\sim$43.1 \Mearth) is within the range predicted by semi-convective interior models.

preprint2014arXiv

Simulations of Titan's paleoclimate

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.

preprint2014arXiv

Water: from clouds to planets

Results from recent space missions, in particular Spitzer and Herschel, have lead to significant progress in our understanding of the formation and transport of water from clouds to disks, planetesimals, and planets. In this review, we provide the underpinnings for the basic molecular physics and chemistry of water and outline these advances in the context of water formation in space, its transport to a forming disk, its evolution in the disk, and finally the delivery to forming terrestrial worlds and accretion by gas giants. Throughout, we pay close attention to the disposition of water as vapor or solid and whether it might be subject to processing at any stage. The context of the water in the solar system and the isotopic ratios (D/H) in various bodies are discussed as grounding data point for this evolution. Additional advances include growing knowledge of the composition of atmospheres of extra-solar gas giants, which may be influenced by the variable phases of water in the protoplanetary disk. Further, the architecture of extra-solar systems leaves strong hints of dynamical interactions, which are important for the delivery of water and subsequent evolution of planetary systems. We conclude with an exploration of water on Earth and note that all of the processes and key parameters identified here should also hold for exoplanetary systems.

preprint2013arXiv

Targeting Young Stars with Kepler: Planet Formation, Migration Mechanisms and the Early History of Planetary Systems

This white paper discusses a repurposed mission for the Kepler spacecraft that focusses on solving outstanding problems in planet formation and evolution by targeting the study of the hot Jupiter population of young stars. This mission can solve the question of the mode of migration of hot Jupiters, address the problem of whether Jupiters form by hot-start (gravitational instability) or cold-start (core accretion) mechanisms, and provide a wealth of data on the early stages of planetary system evolution during the active phases of stars which impact planetary habitability. In one year of observations of three weeks dwell time per field, Kepler would increase by more than an order of magnitude the number of known hot Jupiters, which can be followed up with fast cadence observations to to search for transit timing variations and to perform asteroseismological characterization of the host stars. This mission scenario continues to operate Kepler in the photometric monitoring mode for which it was designed, and is generally flexible with regards to field selection enabling prioritization of fuel usage and attitude control constraints.

preprint2012arXiv

Building Terrestrial Planets

This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth.

preprint2012arXiv

Chemical Constraints on the Oxygen Abundances in Jupiter and Saturn

We perform a comparative analysis of the chemical kinetics of CO and $\rm PH_3$ in Jupiter and Saturn to assess the full set of constraints available on the troposphere water abundance in the two giant planets. For carbon monoxide we employ both a widely used CO kinetic scheme from Yung et al, and a newly identified CO chemical scheme from Visscher and Moses. For $\rm PH_3$ chemical scheme, we use the same chemical scheme as in Visscher and Fegley. Yung's chemical scheme for CO yields a water enrichment of 0.95 - 23.0 times solar abundance on Jupiter, and an upper limit of 14.0 for Saturn. Visscher's chemical scheme in contrast produces a water enrichment of 0.24 - 2.6 times solar abundance in Jupiter, and for Saturn an upper limit for water enrichment of 8.0. From this scheme, which takes advantage of the most up-to-date kinetics, we preclude high water enrichments on Jupiter and Saturn, and show that the kinetics approach yields Jovian bulk abundance in which values of C/O elevated relative to solar are admissible. Our result is consistent with recent reinterpretation of Galileo Probe data in which Jupiter formed in a water-depleted portion of the protoplanetary disk (Mousis et al).

preprint2012arXiv

Nebular water depletion as the cause of Jupiter's low oxygen abundance

Motivated by recent spectroscopic observations suggesting that atmospheres of some extrasolar giant-planets are carbon-rich, i.e. carbon/oxygen ratio (C/O) $\ge$ 1, we find that the whole set of compositional data for Jupiter is consistent with the hypothesis that it be a carbon-rich giant planet. We show that the formation of Jupiter in the cold outer part of an oxygen-depleted disk (C/O $\sim$1) reproduces the measured Jovian elemental abundances at least as well as the hitherto canonical model of Jupiter formed in a disk of solar composition (C/O = 0.54). The resulting O abundance in Jupiter's envelope is then moderately enriched by a factor of $\sim$2 $\times$ solar (instead of $\sim$7 $\times$ solar) and is found to be consistent with values predicted by thermochemical models of the atmosphere. That Jupiter formed in a disk with C/O $\sim$1 implies that water ice was heterogeneously distributed over several AU beyond the snow line in the primordial nebula and that the fraction of water contained in icy planetesimals was a strong function of their formation location and time. The Jovian oxygen abundance to be measured by NASA's Juno mission en route to Jupiter will provide a direct and strict test of our predictions.

preprint2012arXiv

Planetesimal Compositions in Exoplanet Systems

We have used recent surveys of the composition of exoplanet host stars to investigate the expected composition of condensed material in planetesimals formed beyond the snow line in the circumstellar nebulae of these systems. Of the major solid forming elements, we find that, as for the Sun, the C and O abundances (and particularly the C/O abundance ratio) have the most significant effect on the composition of icy planetesimals formed in these systems. The calculations use a self-consistent model for the condensation sequence of volatile ices from the nebula gas after refractory (silicate and metal) phases have condensed. The resultant proportions of refractory phases and ices were calculated for a range of nebular temperature structure and redox conditions. Planetesimals in systems with sub-solar C/O should be water ice-rich, with lower than solar mass fractions of refractory materials, while in super-solar C/O systems planetesimals should have significantly higher fractions of refractories, in some cases having little or no water ice. C-bearing volatile ices and clathrates also become increasingly important with increasing C/O depending on the assumed nebular temperatures. These compositional variations in early condensates in the outer portions of the nebula will be significant for the equivalent of the Kuiper Belt in these systems, icy satellites of giant planets and the enrichment (over stellar values) of volatiles and heavy elements in giant planet atmospheres.

preprint2012arXiv

The dual origin of the nitrogen deficiency in comets: selective volatile trapping in the nebula and postaccretion radiogenic heating

We propose a scenario that explains the apparent nitrogen deficiency in comets in a way consistent with the fact that the surfaces of Pluto and Triton are dominated by nitrogen-rich ice. We use a statistical thermodynamic model to investigate the composition of the successive multiple guest clathrates that may have formed during the cooling of the primordial nebula from the most abundant volatiles present in the gas phase. These clathrates agglomerated with the other ices (pure condensates or stoichiometric hydrates) and formed the building blocks of comets. We report that molecular nitrogen is a poor clathrate former, when we consider a plausible gas phase composition of the primordial nebula. This implies that its trapping into cometesimals requires a low disk temperature ($\sim$20 K) in order to allow the formation of its pure condensate. We find that it is possible to explain the lack of molecular nitrogen in comets as a consequence of their postformation internal heating engendered by the decay of short-lived radiogenic nuclides. This scenario is found consistent with the presence of nitrogen-rich ice covers on Pluto and Triton. Our model predicts that comets should present xenon-to-water and krypton-to-water ratios close to solar xenon-to-oxygen and krypton-to-oxygen ratios, respectively. In contrast, the argon-to-water ratio is predicted to be depleted by a factor of $\sim$300 in comets compared to solar argon-to-oxygen, as a consequence of poor trapping efficiency and radiogenic heating.

preprint2011arXiv

Carbon-rich Giant Planets: Atmospheric Chemistry, Thermal Inversions, Spectra, and Formation conditions

The recent inference of a carbon-rich atmosphere, with C/O >= 1, in the hot Jupiter WASP-12b motivates the exotic new class of carbon-rich planets (CRPs). We report a detailed study of the atmospheric chemistry and spectroscopic signatures of carbon-rich giant planets (CRGs) and the compositions of icy planetesimals required for their formation, and the apportionment of ices, rock, and volatiles in their envelopes. For C/O >= 1, most of the atmospheric oxygen is occupied by CO for T > 1400 K and pressure (P) < 1 bar, causing a substantial depletion in H2O, and an overabundance of CH4 compared to those obtained by assuming solar abundances (C/O = 0.54) in chemical equilibrium. These differences in chemistry cause distinctly observable signatures in spectra. We also find that a C/O >= 1 strongly depletes the abundances of TiO and VO available to form thermal inversions, which is adequate to rule out thermal inversions due to TiO/VO even in the most highly irradiated hot Jupiters, such as WASP-12b. Adopting stellar abundances (C/O = 0.44) for the primordial disk composition and low-temperature formation conditions (T <= 30 K) for WASP-12b leads to a C/O ratio of 0.27 in accreted planetesimals, and, consequently, in the planet's envelope. In contrast, a C/O ratio of 1 in the envelope of WASP-12b requires a substantial depletion of oxygen in the disk, by a factor of ~0.41 for the same formation conditions. This scenario also satisfies the constraints on the C/H and O/H ratios reported for WASP-12b. If, alternatively, hotter conditions prevailed in a stellar composition disk such that only H2O is condensed, the remaining gas can potentially have a C/O ~ 1. However, a high C/O in WASP-12b caused predominantly by gas accretion would preclude super-stellar C/H ratios which also fit the data.

preprint2011arXiv

Removal of Titan's Atmospheric Noble Gases by their Sequestration in Surface Clathrates

A striking feature of the atmosphere of Titan is that no heavy noble gases other than argon were detected by the Gas Chromatograph Mass Spectrometer (GCMS) aboard the Huygens probe during its descent to Titan's surface in January 2005. Here we provide an explanation of the mysterious absence or rarity of these noble gases in Titan's atmosphere: the thermodynamic conditions prevailing at the surface-atmosphere interface of the satellite allow the formation of multiple guest clathrates that preferentially store some species, including all heavy noble gases, over others. The clean water ice needed for formation of these clathrates could be delivered by successive episodes of cryovolcanic lavas that have been hypothesized to regularly cover the surface of Titan. The formation of clathrates in the porous lavas and their propensity for trapping Ar, Kr and Xe would progressively remove these species from the atmosphere of Titan over its history. In some circumstances, a global clathrate crust with an average thickness not exceeding a few meters could be sufficient on Titan for a complete removal of the heavy noble gases from the atmosphere.

preprint2010arXiv

Impact regimes and post-formation sequestration processes: implications for the origin of heavy noble gases in terrestrial planets

The difference between the measured atmospheric abundances of neon, argon, krypton and xenon for Venus, the Earth and Mars is striking. Because these abundances drop by at least two orders of magnitude as one moves outward from Venus to Mars, the study of the origin of this discrepancy is a key issue that must be explained if we are to fully understand the different delivery mechanisms of the volatiles accreted by the terrestrial planets. In this work, we aim to investigate whether it is possible to quantitatively explain the variation of the heavy noble gas abundances measured on Venus, the Earth and Mars, assuming that cometary bombardment was the main delivery mechanism of these noble gases to the terrestrial planets. To do so, we use recent dynamical simulations that allow the study of the impact fluxes of comets upon the terrestrial planets during the course of their formation and evolution. Assuming that the mass of noble gases delivered by comets is proportional to rate at which they collide with the terrestrial planets, we show that the krypton and xenon abundances in Venus and the Earth can be explained in a manner consistent with the hypothesis of cometary bombardment. In order to explain the krypton and xenon abundance differences between the Earth and Mars, we need to invoke the presence of large amounts of CO2-dominated clathrates in the Martian soil that would have efficiently sequestered these noble gases.

preprint2009arXiv

A primordial origin for the atmospheric methane of Saturn's moon Titan

The origin of Titan's atmospheric methane is a key issue for understanding the origin of the Saturnian satellite system. It has been proposed that serpentinization reactions in Titan's interior could lead to the formation of the observed methane. Meanwhile, alternative scenarios suggest that methane was incorporated in Titan's planetesimals before its formation. Here, we point out that serpentinization reactions in Titan's interior are not able to reproduce the deuterium over hydrogen (D/H) ratio observed at present in methane in its atmosphere, and would require a maximum D/H ratio in Titan's water ice 30% lower than the value likely acquired by the satellite during its formation, based on Cassini observations at Enceladus. Alternatively, production of methane in Titan's interior via radiolytic reactions with water can be envisaged but the associated production rates remain uncertain. On the other hand, a mechanism that easily explains the presence of large amounts of methane trapped in Titan in a way consistent with its measured atmospheric D/H ratio is its direct capture in the satellite's planetesimals at the time of their formation in the solar nebula. In this case, the mass of methane trapped in Titan's interior can be up to 1,300 times the current mass of atmospheric methane.

preprint2009arXiv

Elemental abundances and minimum mass of heavy elements in the envelope of HD 189733b

Oxygen (O) and carbon (C) have been inferred recently to be subsolar in abundance from spectra of the atmosphere of the transiting hot Jupiter HD 189733b. Yet, the mass and radius of the planet coupled with structure models indicate a strongly supersolar abundance of heavy elements in the interior of this object. Here we explore the discrepancy between the large amount of heavy elements suspected in the planet's interior and the paucity of volatiles measured in its atmosphere. We describe the formation sequence of the icy planetesimals formed beyond the snow line of the protoplanetary disk and calculate the composition of ices ultimately accreted in the envelope of HD 189733b on its migration pathway. This allows us to reproduce the observed volatile abundances by adjusting the mass of ices vaporized in the envelope. The predicted elemental mixing ratios should be 0.15--0.3 times solar in the envelope of HD 189733b if they are fitted to the recent O and C determinations. However, our fit to the minimum mass of heavy elements predicted by internal structure models gives elemental abundances that are 1.2--2.4 times oversolar in the envelope of HD189733b. We propose that the most likely cause of this discrepancy is irradiation from the central star leading to development of a radiative zone in the planet's outer envelope which would induce gravitational settling of elements. Hence, all strongly irradiated extrasolar planets should present subsolar abundances of volatiles. We finally predict that the abundances of nitrogen (N), sulfur (S) and phosphorus (P) are of $\sim$ $2.8 \times 10^{-5}$, $5.3 \times 10^{-6}$ and $1.8 \times 10^{-7}$ relative to H$_2$, respectively in the atmosphere of HD 189733b.

preprint2008arXiv

Clathration of Volatiles in the Solar Nebula and Implications for the Origin of Titan's atmosphere

We describe a scenario of Titan's formation matching the constraints imposed by its current atmospheric composition. Assuming that the abundances of all elements, including oxygen, are solar in the outer nebula, we show that the icy planetesimals were agglomerated in the feeding zone of Saturn from a mixture of clathrates with multiple guest species, so-called stochiometric hydrates such as ammonia hydrate, and pure condensates. We also use a statistical thermodynamic approach to constrain the composition of multiple guest clathrates formed in the solar nebula. We then infer that krypton and xenon, that are expected to condense in the 20-30 K temperature range in the solar nebula, are trapped in clathrates at higher temperatures than 50 K. Once formed, these ices either were accreted by Saturn or remained embedded in its surrounding subnebula until they found their way into the regular satellites growing around Saturn. In order to explain the carbon monoxide and primordial argon deficiencies of Titan's atmosphere, we suggest that the satellite was formed from icy planetesimals initially produced in the solar nebula and that were partially devolatilized at a temperature not exceeding 50 K during their migration within Saturn's subnebula. The observed deficiencies of Titan's atmosphere in krypton and xenon could result from other processes that may have occurred both prior or after the completion of Titan. Thus, krypton and xenon may have been sequestrated in the form of XH3+ complexes in the solar nebula gas phase, causing the formation of noble gas-poor planetesimals ultimately accreted by Titan. Alternatively, krypton and xenon may have also been trapped efficiently in clathrates located on the satellite's surface or in its atmospheric haze.

preprint2008arXiv

Determination of the minimum masses of heavy elements in the envelopes of Jupiter and Saturn

We calculate the minimum mass of heavy elements required in the envelopes of Jupiter and Saturn to match the observed oversolar abundances of volatiles. Because the clathration efficiency remains unknown in the solar nebula, we have considered a set of sequences of ice formation in which the fraction of water available for clathration is varied between 0 and 100 %. In all the cases considered, we assume that the water abundance remains homogeneous whatever the heliocentric distance in the nebula and directly derives from a gas phase of solar composition. Planetesimals then form in the feeding zones of Jupiter and Saturn from the agglomeration of clathrates and pure condensates in proportions fixed by the clathration efficiency. A fraction of Kr and Xe may have been sequestrated by the H3+ ion in the form of stable XeH3+ and KrH3+ complexes in the solar nebula gas phase, thus implying the formation of at least partly Xe- and Kr-impoverished planetesimals in the feeding zones of Jupiter and Saturn. These planetesimals were subsequently accreted and vaporized into the hydrogen envelopes of Jupiter and Saturn, thus engendering volatiles enrichments in their atmospheres, with respect to hydrogen. Taking into account both refractory and volatile components, and assuming plausible molecular mixing ratios in the gas phase of the outer solar nebula, we show that it is possible to match the observed enrichments in Jupiter and Saturn, whatever the clathration efficiency. Our calculations predict that the O/H enrichment decreases from 6.7 to 5.6 times solar (O/H) in the envelope of Jupiter and from 18.1 to 15.4 times solar (O/H) in the envelope of Saturn with the growing clathration efficiency in the solar nebula.

preprint2004arXiv

The formation and habitability of terrestrial planets in the presence of hot jupiters

`Hot jupiters,' giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form (Armitage 2003). We present results of simulations of terrestrial planet formation in the presence of hot jupiters, broadly defined as having orbital radii <= 0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g. 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with hot jupiters.

preprint2003arXiv

Making other Earths: Dynamical Simulations of Terrestrial Planet Formation and Water Delivery

We present results from 42 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary "embryos" (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 42 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 "habitable" planets between 0.9 and 1.1 AU. These planets range from dry worlds to "water worlds" with 100+ oceans of water (1 ocean = 1.5x10^24 g), and vary in mass between 0.23 and 3.85 Earth masses. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call "super embryos," 1-2 Earth mass protoplanets which can serve as the accretion seeds for 2+ Earth mass planets with large water contents.