A model of Jupiter's magnetospheric magnetic field with variable magnetopause flaring

We present a new model of the jovian magnetosphere in which the flaring of the magnetopause boundary can be varied. Magnetopause flaring is expected to vary due to changing conditions in the upstream interplanetary medium, related both to the dynamic pressure of the solar wind, and to changes in the direction of the interplanetary magnetic field. The model includes a tilted dipole field, which is screened by the magnetopause, a tail field current system, and the field of a screened equatorial current disc.     

Keck near-infrared observations of Saturn's E and G rings during Earth's ring plane crossing in August 1995

We present near-infrared (1.24-2.26 μm) images of Saturn's E and G rings which were taken with the W.M. Keck telescope in 1995 August 9-11, during the period that Earth crossed Saturn's ring plane. Our data confirm that the E ring is very blue. Its radial and vertical structure are found to be remarkably similar to that apparent in the HST ringplane crossing data at visible wavelengths, reinforcing models of the ring's peculiar narrow or very steep particle size distribution. Our data show unambiguously that the satellite Tethys is a secondary source of material for the E ring. The G ring is found to be distinctly red, similar in color to Jupiter's main ring, indicative of a (more typical) broad particle size distribution.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

Long term evolution and chaotic diffusion of the insolation quantities of Mars

As the obliquity of Mars is strongly chaotic, it is not possible to give a solution for its evolution over more than a few million years. Using the most recent data for the rotational state of Mars, and a new numerical integration of the Solar System, we provide here a precise solution for the evolution of Mars' spin over 10 to 20 Myr. Over 250 Myr, we present a statistical study of its possible evolution, when considering the uncertainties in the present rotational state. Over much longer time span, reaching 5 Gyr, chaotic diffusion prevails, and we have performed an extensive statistical analysis of the orbital and rotational evolution of Mars, relying on Laskar's secular solution of the Solar System, based on more than 600 orbital and 200,000 obliquity solutions over 5 Gyr. The density functions of the eccentricity and obliquity are specified with simple analytical formulas. We found an averaged eccentricity of Mars over 5 Gyr of 0.0690 with standard deviation 0.0299, while the averaged value of the obliquity is 37.62° with a standard deviation of 13.82°, and a maximal value of 82.035°. We find that the probability for Mars' obliquity to have reached more than 60° in the past 1 Gyr is 63.0%, and 89.3% in 3 Gyr. Over 4 Gyr, the position of Mars' axis is given by a uniform distribution on a spherical cap limited by the obliquity 58.62°, with the addition of a random noise allowing a slow diffusion beyond this limit. We can also define a standard model of Mars' insolation parameters over 4 Gyr with the most probable values 0.068 for the eccentricity and 41.80° for the obliquity.     

The gravitational instability in the dust layer of a protoplanetary disk:: axisymmetric solutions for nonuniform dust density distributions in the direction vertical to the midplane

The gravitational instability in the dust layer of a protoplanetary disk with nonuniform dust density distributions in the direction vertical to the midplane is investigated. The linear analysis of the gravitational instability is performed. The following assumptions are used: (1) One fluid model is adopted, that is, difference of velocities between dust and gas are neglected. (2) The gas is incompressible. (3) Models are axisymmetric with respect to the rotation axis of the disk. Numerical results show that the critical density at the midplane is higher than the one for the uniform dust density distribution by Sekiya (1983, Prog. Theor. Phys. 69, 1116-1130). For the Gaussian dust density distribution, the critical density is 1.3 times higher, although we do not consider this dust density distribution to be realistic because of the shear instability in the dust layer. For the dust density distribution with a constant Richardson number, which is considered to be realized due to the shear instability, the critical density is 2.85 times higher and is independent of the value of the Richardson number. Further, if a constant Richardson number could decrease to the order of 0.001, the gravitational instability would be realized even for the dust to gas surface density ratio with the solar abundance. Our results give a new restriction on planetesimal formation by the gravitational instability.     

Relative velocities among accreting planetesimals in binary systems: The circumprimary case

We investigate classical planetesimal accretion in a binary star system of separation ab?50 AU by numerical simulations, with particular focus on the region at a distance of 1 AU from the primary. The planetesimals orbit the primary, are perturbed by the companion and are in addition subjected to a gas drag force. We concentrate on the problem of relative velocities Δv among planetesimals of different sizes. For various stellar mass ratios and binary orbital parameters we determine regions where Δv exceed planetesimal escape velocities vesc (thus preventing runaway accretion) or even the threshold velocity vero for which erosion dominates accretion. Gaseous friction has two crucial effects on the velocity distribution: it damps secular perturbations by forcing periastron alignment of orbits, but at the same time the size-dependence of this orbital alignment induces a significant Δv increase between bodies of different sizes. This differential phasing effect proves very efficient and almost always increases Δv to values preventing runaway accretion, except in a narrow eb?0 domain. The erosion threshold Δv>vero is reached in a wide (ab,eb) space for small <10-km planetesimals, but in a much more limited region for bigger ?50-km objects. In the intermediate vesc<Δv<vero domain, a possible growth mode would be the type II runaway growth identified by Kortenkamp et al. [Kortenkamp, S., Wetherill, G., Inaba, S., 2001. Science 293, 1127-1129].     

Near-infrared adaptive optics imaging of the satellites and individual rings of Uranus

We present the first Earth-based images of several of the individual faint rings of Uranus, as observed with the adaptive optics system on the W.M. Keck II telescope on four consecutive days in October 2003. We derive reflectivities based on multiple measurements of 8 minor moons of Uranus as well as Ariel and Miranda in filters centered at wavelengths of 1.25(J), 1.63(H), and 2.1(Kp) μm. These observations have a phase angle of 1.84°-1.96°. We find that the small satellites are somewhat less bright than in observations made by the HST at smaller phase angles, confirming an opposition surge effect. We calculate albedoes for the ring groups and for each ring separately. We find that the ε ring particles, as well as the particles in the three other ring groups, have albedoes near 0.043 at these phase angles. The equivalent depths of some of the individual rings are different than predicted based upon ring widths from occultation measurements (assuming a constant particle ring brightness); in particular the γ ring is fainter and the η ring brighter than expected. Our results indicate that q, the ratio of ε ring intensity at apoapse vs. periapse, is close to 3.2±0.16. This agrees well with a model that has a filling factor for the ε ring of 0.06 (Karkoschka, 2001, Icarus 151, 78-83). We also determine values of the north to south brightness ratio for the individual rings and find that in most cases they are close to unity.     

Thermal response of Iapetus to an eclipse by Saturn's rings

Measurements of Iapetus as seen at 20 and 2.2 μm in the shadow of Saturn's ring are given, providing the thermal response to a rapidly varying heat input. The 20 μm thermal emission follows the 2.2 μm flux input closely. The observations, plus a simple diffusion calculation, imply that the surface of Iapetus is made of material having a very small thermal inertia, probably .