Comparative thermal evolution of Uranus and Neptune

We extend a Jovian convective-cooling model to Uranus and Neptune. The model assumes that efficient interior convection prevails, so that escape of interior heat is governed by the radiative properties of the atmosphere. A comparison of the thermal evolution of Uranus and Neptune indicates that the larger amount of solar radiation absorbed in Uranus' atmosphere tends to differentially suppress the escape of interior heat. The model is shown to be consistent with recent infrared observations of the thermal balance of Uranus and Neptune, and with the presumed age of these planets.     

Uranus and Neptune interior models

This paper is concerned with the interior structure of Uranus and Neptune. Our approach is three-fold. First, a set of three-layer models for both Uranus and Neptune are constructed using a method similar to that used in the study of the terrestrial planets. The variations of the mass density (s) and flattening e(s) with fractional mean radius s for two representative models of Uranus and Neptune are calculated. The results are tabulated. A comparison of these models shows that these two planets are probably very similar to each other in their basic dynamical features. Such similarity is very seldom seen in our solar system. Secondly, we check the conformance between the theoretical results and observational data for the two planets. And thirdly, the 6th degree Stokes zonal parameters for Uranus and for Neptune are predicted, based on the interior models put forward in this paper.     

Giant planet formation

The accumulation of giant planets involves processes typical for terrestrial planet formation as well as gasdynamic processes that were previously known only in stars. The condensible element cores of the gas-giants grow by solid body accretion while envelope formation is governed by stellar-like equilibria and the dynamic departures thereof. Two hypotheses for forming Uranus/Neptune-type planets — at sufficiently large heliocentric distances while allowing accretion of massive gaseous envelopes, i.e. Jupiter-type planets at intermediate distances — have been worked out in detailed numerical calculations: (1) Hydrostatic gas-accretion models with time-dependent solid body accretion-rates show a slow-down of core-accretion at the appropriate masses of Uranus and Neptune. As a consequence, gas-accretion also stagnates and a window is opened for removing the solar nebula during a time of roughly constant envelope mass. (2) Gasdynamic calculations of envelope accretion for constant planetesimal accretion-rates show a dynamic transition to new envelope equilibria at the so called critical mass. For a wide range of solar nebula conditions the new envelopes have respective masses similar to those of Uranus and Neptune and are more tightly bound to the cores. The transitions occur under lower density conditions typical for the outer parts of the solar nebula, whereas for higher densities, i.e. closer to the Sun, gasdynamic envelope accretion sets in and is able to proceed to Jupiter-masses.     

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.     

Can magnetic fields be generated in the icy mantles of Uranus and Neptune?

It has been shown by us previously that a hydromagnetic dynamo can operate in the core of Uranus but probably not on Neptune. A similar analysis is made for the “icy” liquid mantles of both planets. It is concluded that pressure ionization and the associated increased conductivity of water is probably not enough to satisfy the necessary conditions for a dynamo on Uranus and that it is marginal for Neptune. On the other hand the expected presence of metallic water in a thick layer around the core of Neptune makes the operation of a dynamo on this planet plausible. A similar layer on Uranus might be too thin to play the same role. It appears that if a magnetic field is indeed present on Uranus it is probably generated in the core of the planet, while on Neptune it is more likely operating in the icy mantle.     

Models of the giant planets

Models of the giant planets were constructed based on the assumption that the hydrogen to helium ratio is solar in these planets. This assumption, together with arguments about the condensation sequence in the primitive solar nebula, yields models with a central core of rock and possibly ice surrounded by an envelope of hydrogen, helium, methane, ammonia, and water. These last three volatiles may be individually enhanced due to condensation at the period of core formation. Jupiter was found to have a core of about 40 earth masses and a water enhancement in the atmosphere of about 7.5 times the solar value. Saturn was found to have a core of 20 earth masses and a water enhancement in the atmosphere of about 25 times the solar value. Rock plus ice constitute 75–85% of the mass of Uranus and Neptune. Temperatures in the interiors of these planets are probably above the melting points, if there is an adiabatic relation throughout the interiors. Some aspects of the sensitivities of these results to uncertainties in rotational flattening are discussed.     

The phase diagram of water and the magnetic fields of Uranus and Neptune

The interior of giant planets can give valuable information on formation and evolution processes of planetary systems. However, the interior and evolution of Uranus and Neptune is still largely unknown. In this paper, we compare water-rich three-layer structure models of these planets with predictions of shell structures derived from magnetic field models. Uranus and Neptune have unusual non-dipolar magnetic fields contrary to that of the Earth. Extensive three-dimensional simulations of Stanley and Bloxham (Stanley, S., Bloxham, J. [2004]. Nature 428, 151-153) have indicated that such a magnetic field is generated in a rather thin shell of at most 0.3 planetary radii located below the H/He rich outer envelope and a conducting core that is fluid but stably stratified. Interior models rely on equation of state data for the planetary materials which have usually considerable uncertainties in the high-pressure domain. We present interior models for Uranus and Neptune that are based on ab initio equation of state data for hydrogen, helium, and water as the representative of all heavier elements or ices. Based on a detailed high-pressure phase diagram of water we can specify the region where superionic water should occur in the inner envelope. This superionic region correlates well with the location of the stably-stratified region as found in the dynamo models. Hence we suggest a significant impact of the phase diagram of water on the generation of the magnetic fields in Uranus and Neptune.     

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.     

Dynamical evolution of a cometary swarm in the outer planetary region

The dynamical evolution of bodies under the gravitational influence of the accreting proto-Uranus and proto-Neptune is investigated. The main aim of this study is to analyze the interrelations between the accretion of Uranus and Neptune with other processes of cosmological importance as, for example, the formation of a cometary reservoir from bodies placed into near-parabolic orbits by planetary perturbations and the scattering of bodies to the region of the terrestrial planets. Starting with a mass ratio (initial mass/present mass) of 0.1, Uranus and Neptune acquire masses close to their present ones in a time scale of 10⁸ years. Neptune is found to be the most important contributor of comets to the cometary reservoir. The time scale of bodies scattered by Neptune to reach near-parabolic orbits (semimajor axes a > 10⁴ AU)is about 10⁹ years. The contribution of Uranus was partially inhibited because a large part of the residual bodies of its accretion zone fell under the strong gravitational influence of Jupiter and Saturn. A significant fraction of the bodies dispersed by Uranus and Neptune reached the region of the terrestrial planets in a time scale of some 10⁸ years.     

Submillimeter and millimeter observations of Uranus and Neptune

We have made narrowband photometric measurements of Uranus and Neptune covering the wavelength range from 0.35 to 3.3 mm. The observations provide accurate comparative radiometry of these planets. Absolute calibration was referenced to Mars, and to Jupiter as a secondary standard. The results establish Uranus and Neptune as reliable secondary calibrators in their own right. We have combined our observations with other measurements made in the period 1978 through 1984 in the spectral range of 17 μm through 3 mm to form models for atmospheric temperature structure in the vertical range from 100 mbar to 8 bar. The simplest models imply that the tropospheres of both planets are consistent with “frozen” equilibrium H₂ and a mixing ratio of CH₄ of about 2% by volume in the deep atmosphere. There is some evidence in the Uranus data which implies the presence of discrete spectral lines. These could be due to CH₄ pure rotational or dimer transitions or to minor constituents such as CO, which remain uncondensed even at the cold temperatures in the atmosphere of Uranus.     

Some dynamical aspects of the accretion of Uranus and Neptune: The exchange of orbital angular momentum with planetesimals

The final stage of the accretion of Uranus and Neptune is numerically investigated. The four Jovian planets are considered with Jupiter and Saturn assumed to have reached their present sizes, whereas Uranus and Neptune are taken with initial masses 0.2 of their present ones. Allowance is made for the orbital variation of the Jovian planets due to the exchange of angular momentum with interacting bodies (“planetesimals”). Two possible effects that may have contributed to the accretion of Uranus and Neptune are incorporated in our model: (1) an enlarged cross section for accretion of incoming planetesimals due to the presence of extended gaseous envelopes and/or circumplanetary swarms of bodies; and (2) intermediate protoplanets in mid-range orbits between the Jovian planets. Significant radial displacements are found for Uranus and Neptune during their accretion and scattering of planetesimals. The orbital angular momentum budgets of Neptune, Uranus, and Saturn turn out to be positive; i.e., they on average gain orbital angular momentum in their interactions with planetesimals and hence they are displaced outwardly. Instead, Jupiter as the main ejector of bodies loses orbital angular momentum so it moves sunward. The gravitational stirring of planetesimals caused by the introduction of intermediate protoplanets has the effect that additional solid matter is injected into the accretion zones of Uranus and Neptune. For moderate enlargements of the radius of the accretion cross section (2–4 times), the accretion time scale of Uranus and Neptune are found to be a few 10⁸ years and the initial amount of solid material required to form them of a few times their present masses. Given the crucial role played by the size of the accretion cross section, questions as to when Uranus and Neptune acquired their gaseous envelopes, when the envelopes collapsed onto the solid cores, and how massive they were are essential in defining the efficiency and time scale of accretion of the two outer Jovian planets.     

Theoretical predictions of deuterium abundances in the Jovian planets

Using current concepts for the origin of the Jovian planets and current constraints on their interior structure, we argue that the presence of large amounts of “ice” (H₂O, CH₄, and NH₃) in Uranus and Neptune indicates temperatures low enough to condense these species at the time Uranus and Neptune formed. Yet such low temperatures imply orders-of-magnetude fractionation effects for deuterium into the “ice” component if isotopic equilibration can occur. Our models thus imply that Uranus and Neptune should have a D/H ratio at least four times primordial, contrary to observation for Uranus. We find that the Jovian and Saturnian D/H should be close to primordial regardless of formation scenario. The Uranus anomaly could indicate that there was a strong initial radial gradient in D/H in the primordial solar nebula, or that Uranus is so inactive that no significant mixing of its interior has occurred over the age of the solar system. Observation of Neptune's atmospheric D/H may help to resolve the problem.     

The internal structures and the relative rotation rates of Uranus and Neptune

The difference between the interior structures of Uranus and Neptune is presented, based on models which fit the observed mass, radius, and gravitational moments for the assumed rotation periods of these planets. If Uranus and Neptune are assumed to be as similar in internal structure as they are in mass and radius, the rotation period for Neptune must be shorter than that for Uranus. It is suggested that the true rotation period is given by Neptune's oblateness, while the photometric period corresponds to the motion of Rossby waves in the upper atmosphere.     

The formation of the outer planets

One of the outstanding problems in planetary cosmogony is to account for the depletion of hydrogen in the outer planets, Neptune and Uranus. It is suggested that these planets were originally similar to the major planets but that the settling towards the centre of grains, enriched by substances such as methane, ammonia and water because of the low temperatures, released enough energy to cause the evaporation of most of the hydrogen.     

Deuterium enrichment in giant planets

Recently published laboratory measurements of the isotopic exchange rate constant k⁻(T) between CD₄ and H₂ are used to calculate f(z)—the isotopic enrichment factor between CH₄ and H₂—at every level in the outer atmosphere of the giant planets. The variation of f(z) with local vertical velocity, temperature and pressure has been calculated under the assumption that atmospheres are convective and uncertainties have been calculated by error propagation. Considering only the random errors—mainly the uncertainty on k⁻(T)—the f values in the observable upper atmospheres of giant planets (i.e. at z = 0, P = 1 bar) are: f(0) = 1.25 ± 0.05, 1.38 ± 0.06, 1.68 ± 0.09, and 1.61 ± 0.08 for Jupiter, Saturn, Uranus, and Neptune, respectively. Additional systematic errors due to the uncertainty in calculating the vertical velocity in the framework of the mixing length Prandtl theory lead to an overall uncertainty on f(0) of ±0.12, ±0.15, ±0.23, and ±0.21 for each planet, respectively. The D/H ratios in H₂ derived from the measured CH₃D/CH₄ ratios in the upper atmosphere of the four giant planets are then recalculated. Uranus and Neptune seem to be enriched in deuterium with respect to the protosolar nebula but depleted relative to the Standard Mean Oceanic Water on the Earth (SMOW). However calculations based on current interior models of Neptune suggest that ices which formed the core of the planet had a D/H ratio of the order of the SMOW. The deuterium abundance in proto-Uranian ices remains uncertain. The case where water is a major constituent of the fluid envelope of Neptune is discussed. It is shown that the D/H ratio of the planet would then be higher than the value measured in hydrogen. Even in this case, the D/H ratio in proto-Neptunian ices is less than the recently revised value in P/Halley and less than the value measured in water of the Semarkona meteorite. These results suggest that the ices which formed the core of Neptune did not have an interstellar origin. Similarly, the comparison of the most recent determination of the D/H ratio in the atmosphere of Titan with the value of D/H in P/Halley suggests that this atmosphere was not formed by infalling comets but more likely from grains embedded in the sub-nebula of Saturn.     

试论行星磁场强度分布的成因

行星磁场强度的分布具有规律性,中间巨行星的强度值大,两侧类地行星和远日行星的强度值小。在类地行星和远日行星中,也是中间行星磁场强度值大,如地球和海王星。两侧行星磁场强度值小,如金星、火星、天王星和冥王星。水星磁场强度值比金星大是因为水星环较早的从A组星云环中分裂出来所致。行星磁场强度大小分布的规律性与星云环级式分裂、环体物质中聚即铁磁性物质较快中聚的规律性一致。所以,铁磁性物质多少,可能是行星磁场强度大小分布的一个主要原因。     

The geometric albedos of Uranus and Neptune between 2100 and 3350 Å

Absolutely calibrated spectra of Uranus, Neptune, and the solar analog stars 16 Cyg A and B between 2100 and 3350 Å are reported. The geometric albedos of both planets are close to the curve expected for a semi-infinite Rayleigh-Raman scattering atmosphere between 2100 and 2800 Å. Longward of 2800 Å the albedos fall below the Rayleigh-Raman values and connect smoothly to the ground-based photometry of J.S. Neff, D.C. Humm, J.T. Bergstralh, A.L. Cochran, W.D. Cochran, E.S. Barker, and R.G. Tull (1984, Icarus60, 221–235).Neptune is about 5.5% brighter than Uranus and shows slightly stronger Raman scattering signatures in the MgII lines at 2800 Å in accordance with the results of Neff et al. for the CaII H and K lines. This means that the stratospheric haze on Neptune is thinner than on Uranus. The fact that the Neptunian geometric albedo between 2100 and 2800 Å is so close to the ideal semi-infinite Rayleigh-Raman scattering atmosphere could be exploited for future absolute calibrations of other Solar System objects in this wavelength region.     

Evidence for Planetary Forcing of the Cosmic Ray Intensity and Solar Activity Throughout the Past 9400 Years

Paleo-cosmic-ray (PCR) records based on cosmogenic ¹⁰Be and ¹⁴C data are used to study the variations in cosmic-ray intensity and solar activity over the past 9400 years. There are four strong correlations with the motion of the Jovian planets; the probability of occurring by chance being <?10^?5. They are i) the PCR periodicities at 87, 350, 510, and 710 years, which closely approximate integer multiples of half the Uranus–Neptune synodic period; ii) eight periodicities in the torques calculated to be exerted by the planets on an asymmetric tachocline that approximate the periods observed in the PCR; iii) the maxima of the long-term PCR variations are coincident with syzygy (alignment) of the four Jovian planets in 5272 and 644 BP; and iv) in the time domain, the PCR intensity decreases during the first 60 years of the ≈?172 year Jose cycle (Jose, Astron. J. 70, 193, 1965) and increases in the remaining ≈?112 years in association with barycentric anomalies in the distance between the Sun and the center of mass of the solar system. Furthermore, sunspot and neutron-monitor data show that three anomalous sunspot cycles (4th, 7th, and 20th) and the long sunspot minimum of 2006 – 2009 CE coincided with the first and second barycentric anomalies of the 58th and 59th Jose cycles. Phase lags between the planetary and heliospheric effects are ≤?five years. The 20 largest Grand Minima during the past 9400 years coincided with the latter half of the Jose cycle in which they occurred. These correlations are not of terrestrial origin, nor are they due to the planets’ contributing directly to the cosmic-ray modulation process in the heliosphere. Low cosmic-ray intensity (higher solar activity) occurred when Uranus and Neptune were in superior conjunction (mutual cancellation), while high intensities occurred when Uranus–Neptune were in inferior conjunction (additive effects). Many of the prominent peaks in the PCR Fourier spectrum can be explained in terms of the Jose cycle, and the occurrence of barycentric anomalies.     

Uranus and Neptune: Shape and rotation

Both Uranus and Neptune are thought to have strong zonal winds with velocities of several 100 m s⁻¹. These wind velocities, however, assume solid-body rotation periods based on Voyager 2 measurements of periodic variations in the planets’ radio signals and of fits to the planets’ magnetic fields; 17.24 h and 16.11 h for Uranus and Neptune, respectively. The realization that the radio period of Saturn does not represent the planet’s deep interior rotation and the complexity of the magnetic fields of Uranus and Neptune raise the possibility that the Voyager 2 radio and magnetic periods might not represent the deep interior rotation periods of the ice giants. Moreover, if there is deep differential rotation within Uranus and Neptune no single solid-body rotation period could characterize the bulk rotation of the planets. We use wind and shape data to investigate the rotation of Uranus and Neptune. The shapes (flattening) of the ice giants are not measured, but only inferred from atmospheric wind speeds and radio occultation measurements at a single latitude. The inferred oblateness values of Uranus and Neptune do not correspond to bodies rotating with the Voyager rotation periods. Minimization of wind velocities or dynamic heights of the 1 bar isosurfaces, constrained by the single occultation radii and gravitational coefficients of the planets, leads to solid-body rotation periods of ∼16.58 h for Uranus and ∼17.46 h for Neptune. Uranus might be rotating faster and Neptune slower than Voyager rotation speeds. We derive shapes for the planets based on these rotation rates. Wind velocities with respect to these rotation periods are essentially identical on Uranus and Neptune and wind speeds are slower than previously thought. Alternatively, if we interpret wind measurements in terms of differential rotation on cylinders there are essentially no residual atmospheric winds.     

On the DC ratio in the atmosphere of Uranus

A method for deriving mixing ratios in the outer planets, mostly independent of scattering processes, is applied to Uranus. It is shown that scattering processes play a major role in the line formation in the atmospheres of Uranus and Neptune; consequently, abundance ratios derived from a reflecting-layer model can be questionable. Using our method, we derive for Uranus $\text{D}\text{C}\text{< 6 × 10}{}^{\mathrm{?3}}$, which is significantly smaller than our result on Jupiter. The simplest explanation implies a C/H enrichment by at least a factor of 6 relative to the solar value.