Generation of equatorial jets by large-scale latent heating on the giant planets

Three-dimensional numerical simulations show that large-scale latent heating resulting from condensation of water vapor can produce multiple zonal jets similar to those on the gas giants (Jupiter and Saturn) and ice giants (Uranus and Neptune). For plausible water abundances (3-5 times solar on Jupiter/Saturn and 30 times solar on Uranus/Neptune), our simulations produce ∼20 zonal jets for Jupiter and Saturn and 3 zonal jets on Uranus and Neptune, similar to the number of jets observed on these planets. Moreover, these Jupiter/Saturn cases produce equatorial superrotation whereas the Uranus/Neptune cases produce equatorial subrotation, consistent with the observed equatorial-jet direction on these planets. Sensitivity tests show that water abundance, planetary rotation rate, and planetary radius are all controlling factors, with water playing the most important role; modest water abundances, large planetary radii, and fast rotation rates favor equatorial superrotation, whereas large water abundances favor equatorial subrotation regardless of the planetary radius and rotation rate. Given the larger radii, faster rotation rates, and probable lower water abundances of Jupiter and Saturn relative to Uranus and Neptune, our simulations therefore provide a possible mechanism for the existence of equatorial superrotation on Jupiter and Saturn and the lack of superrotation on Uranus and Neptune. Nevertheless, Saturn poses a possible difficulty, as our simulations were unable to explain the unusually high speed (∼) of that planet’s superrotating jet. The zonal jets in our simulations exhibit modest violations of the barotropic and Charney-Stern stability criteria. Overall, our simulations, while idealized, support the idea that latent heating plays an important role in generating the jets on the giant planets.     

Capture of planetesimals by the giant planets

We investigate the possibility of gravitational capture of planetesimals as temporary or permanent satellites of Uranus and Neptune during the process of planetary growth. The capture mechanism is based in the enhancement of the Hill's sphere of action not only due to the mass acquired by the planet, but also by the variation of the planet-Sun distance as a consequence of the scattering of planetesimals by the planets of the outer solar system. Our calculations indicate that satellite capture was very important, specially during the first stages of the accretion process, contributing in a significant way to the planetary growth.     

From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets

We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk-planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1-5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1-2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.     

Study of dynamical processes at the final stages of planetary system formation: Resonance motion of giant planets

According to current observational data, planets of many exoplanetary systems have resonant motion. The formation of resonance configurations is studied within a unified model of planetary migration. Planets in the observed systems 24 Sex, HD 37124, HD 73526, HD 82943, HD 128311, HD 160691, Kepler 9, NN Ser, which are moving in the 2: 1 resonance, could have been captured into this resonance due to both the Type I and II migration with a wide range of parameters. The migration conditions are defined for the formation of HD 45364 and HD 200964 that are in the 3: 2 and 4: 3 first-order resonances, correspondingly. The results obtained for HD 200964 show that planets can be captured in the first-order resonances, when the outer-to-inner orbital period ratios for the planets are less than 3: 2, only if Type I migration rates are large, and the mass of at least one planet is substantially less than the modern masses of the observed giant planets. The formation of the HD 102272, HD 108874, HD 181433 and HD 202206 systems with planets in high-order resonances is considered. The capture into these resonances can be realized with very slow Type II migration. Possible bounds for migration parameters are considered. In particular, it has been found that the capture of HD 108874 into the 4: 1 resonance is possible only if the angle between the plane of planetary orbits and the plane of sky is appreciably less than 90°, i.e., the planetary masses are a few times larger than the minimum values. The capture of HD 202206 into the 5: 1 resonance is possible at low migration rates; however, another mechanism is required to explain the high observed eccentricity of the inner planet (for example, strong gravitational interaction between the planets). Resonant configurations can be disrupted due to the interaction between planets and remaining fragments of the planetesimal disk as, for example, may occur in the three-planet system 47 UMa. The specific orbital features observed for this system are explained.     

A seasonal climate model of the atmospheres of the giant planets at the voyager encounter time: I. Saturn's stratosphere

A radiative seasonal model which incorporates a multilayer radiative transfer treatment at wave-lengths longward of 7 μm is presented and applied to Saturn's stratosphere. Opacities due to H₂-He, CH₄, C₂H₂, and C₂H₆ are included. Season-dependent insolation is shown to produce a strong hemispheric asymmetry decreasing with depth at the Voyager encounter times, and seasonal amplitudes of 30°K at the poles are predicted in the high stratosphere. The ring-modulated dependence of the insolation and the orbital eccentricity are shown to have a significant effect. Calculations agree closely with the Voyager 1 and 2 radio occultation ingress profiles recorded at 76°S and 36.5°S for CH₄/H₂ = 3.5 + 1.4/? 1.0 × 10^?3;the estimated errors include modeling systematic errors and uncertainties in the occultations profiles. The possible role of aerosols in the stratospheric heating is analyzed. The Voyager 2 egress profile recorded at 31°S cannot be reproduced by calculations. Some constraints on the C₂H₂ and C₂H₆ abundances are derived. The upper portion of the occultation profiles (p < 3mbar) can be matched for C₂H₂/H₂ = 1.0 + 1.3/?0.6 × 10^?7, C₂H₆/H₂ = 1.5 + 1.8/?0.9 × 10^?6 at 76°S and C₂H₂/H₂ = 4 + 6/?4 × 10^?8, C₂H₆/H₂ = 6 + 9/?6 × 10^?7 at 36.5°N. At the northern occultation latitude, the discrepancy with the concentrations derived from analysis of IRIS spectra by R. Courtin, D. Gautier, A. Marten, B. Bézard, and R. Hanel (1984, Astrophys. J.287) can be explained by a sharp variation of the mixing ratios of these gases with altitude in the upper stratosphere. Other interpretations are discussed.     

Magnetospheric radio emission from extrasolar giant planets: the role of the host stars

We present a new analysis of the expected magnetospheric radio emission from extrasolar giant planets (EGPs) for a distance limited sample of the nearest known extrasolar planets. Using recent results on the correlation between stellar X-ray flux and mass-loss rates from nearby stars, we estimate the expected mass-loss rates of the host stars of extrasolar planets that lie within 20 pc of the Earth. We find that some of the host stars have mass-loss rates that are more than 100 times that of the Sun and, given the expected dependence of the planetary magnetospheric radio flux on stellar wind properties, this has a very substantial effect. Using these results and extrapolations of the likely magnetic properties of the extrasolar planets, we infer their likely radio properties.
We compile a list of the most promising radio targets and conclude that the planets orbiting Tau Bootes, Gliese 86, Upsilon Andromeda and HD 1237 (as well as HD 179949) are the most promising candidates, with expected flux levels that should be detectable in the near future with upcoming telescope arrays. The expected emission peak from these candidate radio emitting planets is typically ∼40–50 MHz. We also discuss a range of observational considerations for detecting EGPs.     

Young Jupiters are faint: new models of the early evolution of giant planets

Here we show preliminary calculations of the cooling and contraction of a 2 M_J planet. These calculations, which are being extended to 1–10 M_J, differ from other published “cooling tracks” in that they include a core accretion‐gas capture formation scenario, the leading theory for the formation of gas giant planets.We find that the initial post‐accretionary intrinsic luminosity of the planet is ∼3 times less than previously published models which use arbitrary initial conditions. These differences last a few tens of millions of years. Young giant planets are intrinsically fainter than has been previously appreciated. We also discuss how uncertainties in atmospheric chemistry and the duration of the formation time of giant planets lead to challenges in deriving planetary physical properties from comparison with tabulated model values. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Upgraded technique to analyze the vertical structure of the aerosol component of the atmospheres of giant planets

To analyze the behavior of the optical thickness of aerosols or the ratio of the optical thicknesses of the aerosol and gas components in the spectral absorption bands of atmospheric gases with depth, we developed a software package. The package structure includes the units for the following operations: (1) to calculate the Legendre expansion coefficients x _i of the phase function and the volume scattering coefficient σ₀ of the polydisperse medium with the specified refractive index and the size distribution function N(r) with the use of the code developed by M.I. Mishchenko; (2) to generate the array containing the pairs of the single scattering albedo ω and the geometric albedo A _g for the models of a semi-infinite homogeneous medium with the parameters determined in the previous step; (3) to determine the single scattering albedo values from the comparison of the calculated and measured values of the geometric albedo for each of the measured points in the examined absorption band of the atmospheric gas (accounting for the change of the phase function due to Rayleigh scattering); (4) to calculate the spectral values of the effective optical depths τ_eff of the levels, where the intensity field of light diffusely reflected by the investigated atmosphere is formed; (5) to derive the scattering and absorption components of the effective optical depth (τ _eff ^s and τ _eff ^v ) from the values of ω and τ_eff; (6) from the values of τ _eff ^v to determine the amount of the absorbing gas NL (in km-amg) along the line of sight and, from these values, the atmospheric pressure p(NL) and the gas component of the scattering portion of the optical depth τ _g (λ₀) at the wavelength λ₀ = 887.2 nm; (7) from the values of τ _eff ^s (λ) and τ _g (λ, NL), to find the aerosol component τ _a (λ, NL); (8) to build the plots of τ _a (λ) or the ratio τ _a (λ)/τ _g (λ) reduced to λ₀ = 887.2 nm versus the pressure. The software package was validated in the analysis of the spectrophotometric data obtained in the measurements of the whole disk of Jupiter in the profiles of the strong absorption bands of methane centered at λλ 841.6, 864, and 887.2 nm under the assumption on two versions of the size distribution of particles (the modified gamma distribution and the log-normal one). It was found that the model with the gamma distribution is analyzed several times more quickly than the log-normal one that yields the close results in computations performed for the same medium.     

The accretion of brown dwarfs and planets by giant stars – II. Solar-mass stars on the red giant branch

This paper extends our previous study of planet/brown dwarf accretion by giant stars to solar-mass stars located on the red giant branch. The model assumes that the planet is dissipated at the bottom of the convective envelope of the giant star. The evolution of the giant is then followed in detail. We analyse the effects of different accretion rates and different initial conditions. The computations indicate that the accretion process is accompanied by a substantial expansion of the star, and, in the case of high accretion rates, hot bottom burning can be activated. The possible observational signatures that accompany the engulfing of a planet are also extensively investigated. They include the ejection of a shell and a subsequent phase of IR emission, an increase in the ⁷Li surface abundance and a potential stellar metallicity enrichment, spin-up of the star because of the deposition of orbital angular momentum, the possible generation of magnetic fields and the related X-ray activity caused by the development of shear at the base of the convective envelope, and the effects on the morphology of the horizontal branch in globular clusters. We propose that the IR excess and high Li abundance observed in 4–8 per cent of the G and K giants originate from the accretion of a giant planet, a brown dwarf or a very low-mass star.     

A three-dimensional model of moist convection for the giant planets II: Saturn's water and ammonia moist convective storms

Moist convective storms constitute a key aspect in the global energy budget of the atmospheres of the giant planets. Among them, Saturn is known to develop the largest scale convective storms in the Solar System, the Great White Spots (GWS) which occur rarely and have been detected once every 30 years approximately. On the average, Saturn seems to show much less convective storms than Jupiter with smaller size and reduced frequency and intensity. Here we present detailed simulations of the onset and development of storms at the Equator and mid-latitudes of Saturn. These are the regions where most of the recent convective activity of the planet has been observed. We use a 3D anelastic model with parameterized microphysics (Hueso and Sánchez-Lavega, 2001, Icarus 151, 257) studying the onset and evolution of water and ammonia moist convective storms up to sizes of a few hundred km. Water storms, while more difficult to initiate than in Jupiter, can be very energetic, arriving to the 150 mbar level and developing vertical velocities on the order of 150 m s⁻¹. Ammonia storms develop easier but with a much smaller intensity unless very large abundances of ammonia (10 times solar) are present in Saturn's atmosphere. The Coriolis forces play a major role in the morphology and properties of water based storms.     

Calculations of the accretion and evolution of giant planets: The effects of solid cores

The evolution of the giant planets is calculated under the general hypothesis that the solid cores formed first, by accretion of small particles, and that these cores later gravitationally attracted their gaseous envelopes from the solar nebula. The evolution passes through the following phases. (1) Planetesimals accrete to form a core of rocky and icy material. (2) When the core mass has grown to a few tenths of an Earth mass, a gaseous envelope in hydrostatic equilibrium begins to form around the core. (3) The core and envelope continue to grow until the “critical” core mass is reached, beyond which point the envelope increases in mass much more rapidly than the core. (4) The envelope mass increases quickly to its present value and prodices a relatively high luminosity, derived from gravitational contraction. (5) Accretion of both core and envelope terminates, and the planet contracts and cools to its present state on a time scale of 5 × 10⁹ years. Evolutionary calculations of phases (2) through (5) are presented, based on solutions of the time-dependent stellar structure equations in spherical symmetry. The physical considerations that determine the critical core mass are discussed; its value is found to depend strongly on the core accretion rate but only weakly on surface boundary conditions. Evolutionary tracks up to the present state are presented for objects of Uranus and Saturn mass.     

The deep wind structure of the giant planets: Results from an anelastic general circulation model

The giant gas planets have hot convective interiors, and therefore a common assumption is that these deep atmospheres are close to a barotropic state. Here we show using a new anelastic general circulation model that baroclinic vorticity contributions are not negligible, and drive the system away from an isentropic and therefore barotropic state. The motion is still aligned with the direction of the axis of rotation as in a barotropic rotating fluid, but the wind structure has a vertical shear with stronger winds in the atmosphere than in the interior. This shear is associated with baroclinic compressibility effects. Most previous convection models of giant planets have used the Boussinesq approximation, which assumes the density is constant in depth; however, Jupiter's actual density varies by four orders of magnitude through its deep molecular envelope. We therefore developed a new general circulation model (based on the MITgcm) that is anelastic and thereby incorporates this density variation. The model's geometry is a full 3D sphere down to a small inner core. It is nonhydrostatic, uses an equation of state suitable for hydrogen-helium mixtures (SCVH), and is driven by an internal heating profile. We demonstrate the effect of compressibility by comparing anelastic and Boussinesq cases. The simulations develop a mean state that is geostrophic and hydrostatic including the often neglected, but significant, vertical Coriolis contribution. This leads to modification of the standard thermal wind relation for a deep compressible atmosphere. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, and provide a mechanism for the transport of heat poleward, which can cause the observed flat meridional emission. We address the issue of over-forcing which is common in such convection models and analyze the dependence of our results on this; showing that the vertical wind structure is not very sensitive to the Rayleigh number. We also study the effect of rotation, showing how the transition from a rapidly to a slowly rotating system affects the dynamics.     

A tale of planet formation: from dust to planets

The characterization of exoplanets and their birth protoplanetary disks has enormously advanced in the last decade. Benefitting from that, our global understanding of the planet formation processes has been substantially improved. In this review, we first summarize the cutting-edge states of the exoplanet and disk observations. We further present a comprehensive panoptic view of modern core accretion planet formation scenarios, including dust growth and radial drift, planetesimal formation by the streaming instability, core growth by planetesimal accretion and pebble accretion. We discuss the key concepts and physical processes in each growth stage and elaborate on the connections between theoretical studies and observational revelations. Finally, we point out the critical questions and future directions of planet formation studies.     

A 2D study of the effects of the size of a solid core on the equatorial flow in giant planets

Zonal winds simulated in two-dimensional computer models of turbulent convection in the equatorial plane of giant planets have greater surface amplitudes for cases with smaller solid cores, and therefore larger buoyancy driving, all other properties being equal. This differential rotation in radius is maintained by the convergence of angular momentum flux, which occurs because of the convective flow that develops due to the effects of planetary rotation and density stratification. The superposition of the convective flow and the stronger zonal flow produces wave-like, instead of cellular convection.     

Formation of the regular satellites of giant planets in an extended gaseous nebula II: satellite migration and survival

For a satellite to survive in the disk the time scale of satellite migration must be longer than the time scale for gas dissipation. For large satellites (∼1000 km) migration is dominated by the gas tidal torque. We consider the possibility that the redistribution of gas in the disk due to the tidal torque of a satellite with mass larger than the inviscid critical mass causes the satellite to stall and open a gap (W.R. Ward, 1997, Icarus 26, 261-281). We adapt the inviscid critical mass criterion to include gas drag, and m-dependent nonlocal deposition of angular momentum. We find that such a model holds promise of explaining the survival of satellites in the subnebula, the mass versus distance relationship apparent in the saturnian and uranian satellite systems, the concentration of mass in Titan, and the observation that the satellites of Jupiter get rockier closer to the planet whereas those of Saturn become increasingly icy. It is also possible that either weak turbulence (close to the planet) or gap-opening satellite tidal torque removes gas on a similar time scale (10⁴-10⁵ years) as the orbital decay time of midsized (200-700 km) regular satellites forming in the inner disk (inside the centrifugal radius (I. Mosqueira and P.R. Estrada, 2003, Icarus, this issue)). We argue that Saturn’s satellite system bridges the gap between those of Jupiter and Uranus by combining the formation of a Galilean-sized satellite in a gas optically thick subnebula with a strong temperature gradient, and the formation of smaller satellites, closer to the planet, in a disk with gas optical depth ?1, and a weak temperature gradient.Using an optically thick inner disk (given gaseous opacity), and an extended, quiescent, optically thin outer disk, we show that there are regions of the disk of small net tidal torque (even zero) where satellites (Iapetus-sized or larger) may stall far from the planet. For our model these outer regions of small net tidal torque correspond roughly to the locations of Callisto and Iapetus. Though the precise location depends on the (unknown) size of the transition region between the inner and outer disks, the result that Saturn’s is found much farther out (at ∼3rcS, where rcS is Saturn’s centrifugal radius) than Jupiter’s (at ∼ 2rcJ, where rcJ is Jupiter’s centrifugal radius) is mostly due to Saturn’s less massive outer disk and larger Hill radius. However, despite the large separation between Ganymede and Callisto and Titan and Iapetus, the long formation and migration time scales for Callisto and Iapetus (I. Mosqueira and P.R. Estrada, 2003, Icarus, this issue) makes it possible (depending on the details of the damping of acoustic waves) that the tidal torque of Ganymede and Titan clears the gas disk out to their location, thus stranding Callisto and Iapetus far from the planet. Either way, our model provides an explanation for the presence of regular satellites outside the centrifugal radii of Jupiter and Saturn, and the absence of such a satellite for Uranus.     

Validation of Sagnier's theory by its application to the major planets

The understanding of quantum mechanical phenomena has come to rely heavily on theory framed in terms of operators and their eigenvalue equations. This paper investigates the utility of that technique as related to the reciprocity principle in diffuse reflection. The reciprocity operator is shown to be unitary and Hermitian; hence, its eigenvectors form a complete orthonormal basis. The relevant eigenvalue is found to be infinitely degenerate. A superposition of the eigenfunctions found from solution by separation of variables is inadequate to form a general solution that can be fitted to a one-dimensional boundary condition, because the difficulty of resolving the reciprocity operator into a superposition of independent one-dimensional operators has yet to be overcome. A particular lunar application in the form of a failed prediction of limb-darkening of the full Moon from brightness versus phase illustrates this problem. A general solution is derived which fully exploits the determinative powers of the reciprocity operator as an unresolved two-dimensional operator. However, a solution based on a sum of one-dimensional operators, if possible, would be much more powerful. A close association is found between the reciprocity operator and the particle-exchange operator of quantum mechanics, which may indicate the direction for further successful exploitation of the approach based on the operational calculus.