Plasma clouds associated with Comet P/Borrelly dust impacts

The NASA DS1 spacecraft encountered Comet P/Borrelly on September 22, 2001 at a distance of ∼2171 km on the sunward side of the comet. The flyby speed was ∼16.5 km s⁻¹. Using high temporal resolution (50 μs) absolute electric field amplitude measurements from a ∼1 m dipole antenna, new features of plasma clouds created by cometary dust impacts have been detected. The pulses have 1/e exponential decays of ∼650 μs duration, exponentially shaped overshoots with rise times of ∼2 ms, and exponential-shaped overshoot decay times of ∼10 ms. Assuming a plasma temperature of 10⁴ K, these pulse features have been explained as plasma cloud space charge effects from the electron, proton and heavy ion portions of the clouds passing the antenna. Complex pulse shapes were also detected. These are believed to be due to either plasma cloud scattering off of the spacecraft, or to secondary impacts. Small electric pulses of duration 10-15 ms of cometary origin were detected but are presently unexplained. The electric component of the plasma wave spectra at closest approach had an f^−2.4 power law shape from 10 Hz to 1 kHz. The electron cyclotron frequency was approximately 1 kHz. One possible explanation of the wave spectrum is that whistler mode waves associated with phase steepened cometary plasma waves are dispersed, leading to the broad spectrum. Finally, based on the present results, a new type of low-cost, large-area dust detector is proposed.     

Trapping and dynamical evolution of interplanetary dust particles in Earth’s quasi-satellite resonance

We used numerical simulations to model the orbital evolution of interplanetary dust particles (IDPs) evolving inward past Earth’s orbit under the influence of radiation pressure, Poynting–Robertson light drag (PR drag), solar wind drag, and gravitational perturbations from the planets. A series of β values (where β is the ratio of the force from radiation pressure to that of central gravity) were used ranging from 0.0025 up to 0.02. Assuming a composition consistent with astronomical silicate and a particle density of 2.5 g cm⁻³ these β values correspond to dust particle diameters ranging from 200 μm down to 25 μm. As the dust particle orbits decay past 1 AU between 4% (for β = 0.02, or 25 μm) and 40% (for β = 0.0025, or 200 μm) of the population became trapped in 1:1 co-orbital resonance with Earth. In addition to traditional horseshoe type co-orbitals, we found about a quarter of the co-orbital IDPs became trapped as so-called quasi-satellites. Quasi-satellite IDPs always remain relatively near to Earth (within 0.1–0.3 AU, or 10–30 Hill radii, R_H) and undergo two close-encounters with Earth each year. While resonant perturbations from Earth halt the decay in semi-major axis of quasi-satellite IDPs their orbital eccentricities continue to decrease under the influence of PR drag and solar wind drag, forcing the IDPs onto more Earth-like orbits. This has dramatic consequences for the relative velocity and distance of closest approach between Earth and the quasi-satellite IDPs. After 10⁴–10⁵ years in the quasi-satellite resonance dust particles are typically less than 10R_H from Earth and consistently coming within about 3R_H. In the late stages of evolution, as the dust particles are escaping the 1:1 resonance, quasi-satellite IDPs can have deep close-encounters with Earth significantly below R_H. Removing the effects of Earth’s gravitational acceleration reveals that encounter velocities (i.e., velocities “at infinity”) between quasi-satellite IDPs and Earth during these close-encounters are just a few hundred meters per second or slower, well below the average values of 2–4 km s⁻¹ for non-resonant Earth-crossing IDPs with similar initial orbits. These low encounter velocities lead to a factor of 10–100 increase in Earth’s gravitationally enhanced impact cross-section (σ_grav) for quasi-satellite IDPs compared to similar non-resonant IDPs. The enhancement in σ_grav between quasi-satellite IDPs and cometary Earth-crossing IDPs is even more pronounced, favoring accretion of quasi-satellite dust particles by a factor of 100–3000 over the cometary IDPs. This suggests that quasi-satellite dust particles may dominate the flux of large (25–200 μm) IDPs entering Earth’s atmosphere. Furthermore, because quasi-satellite trapping is known to be directly correlated with the host planet’s orbital eccentricity the accretion of quasi-satellite dust likely ebbs and flows on 10⁵ year time scales synchronized with Earth’s orbital evolution.     

Decreased values of cosmic dust number density estimates in the Solar System

Experiments to investigate the effect of impacts on side-walls of dust detectors such as the present NASA/ESA Galileo/Ulysses instrument are reported. Side walls constitute 27% of the internal area of these instruments, and increase field of view from 140° to 180°. Impact of cosmic dust particles onto Galileo/Ulysses Al side walls was simulated by firing Fe particles, 0.5-5 μm diameter, 2-50 km s⁻¹, onto an Al plate, simulating the targets of Galileo and Ulysses dust instruments. Since side wall impacts affect the rise time of the target ionization signal, the degree to which particle fluxes are overestimated varies with velocity. Side-wall impacts at particle velocities of 2-20 km s⁻¹ yield rise times 10-30% longer than for direct impacts, so that derived impact velocity is reduced by a factor of ∼2. Impacts on side wall at 20-50 km s⁻¹ reduced rise times by a factor of ∼10 relative to direct impact data. This would result in serious overestimates of flux of particles intersecting the dust instrument at velocities of 20-50 km s⁻¹. Taking into account differences in laboratory calibration geometry we obtain the following percentages for previous overestimates of incident particle number density values from the Galileo instrument [Grün et al., 1992. The Galileo dust detector. Space Sci. Rev. 60, 317-340]: 55% for 2 km s⁻¹ impacts, 27% at 10 km s⁻¹ and 400% at 70 km s⁻¹. We predict that individual particle masses are overestimated by ∼10-90% when side-wall impacts occur at 2-20 km s⁻¹, and underestimated by ∼10-¹⁰² at 20-50 km s⁻¹. We predict that wall impacts at 20-50 km s⁻¹ can be identified in Galileo instrument data on account of their unusually short target rise times. The side-wall calibration is used to obtain new revised values [Krüger et al., 2000. A dust cloud of Ganymede maintained by hypervelocity impacts of interplanetary micrometeoroids. Planet. Space Sci. 48, 1457-1471; 2003. Impact-generated dust clouds surrounding the Galilean moons. Icarus 164, 170-187] of the Galilean satellite dust number densities of 9.4×¹⁰⁻⁵, 9.9×¹⁰⁻⁵, 4.1×¹⁰⁻⁵, and 6.8×¹⁰⁻⁵ m⁻³ at 1 satellite radius from Io, Europa, Ganymede, and Callisto, respectively. Additionally, interplanetary particle number densities detected by the Galileo mission are found to be 1.6×¹⁰⁻⁴, 7.9×¹⁰⁻⁴, 3.2×¹⁰⁻⁵, 3.2×¹⁰⁻⁵, and 7.9×¹⁰⁻⁴ m⁻³ at heliocentric distances of 0.7, 1, 2, 3, and 5 AU, respectively. Work by Burchell et al. [1999b. Acceleration of conducting polymer-coated latex particles as projectiles in hypervelocity impact experiments. J. Phys. D: Appl. Phys. 32, 1719-1728] suggests that low-density “fluffy” particles encountered by Ulysses will not significantly affect our results—further calibration would be useful to confirm this.     

Three-body capture of irregular satellites: Application to Jupiter

We investigate a new theory of the origin of the irregular satellites of the giant planets: capture of one member of a ∼100-km binary asteroid after tidal disruption. The energy loss from disruption is sufficient for capture, but it cannot deliver the bodies directly to the observed orbits of the irregular satellites. Instead, the long-lived capture orbits subsequently evolve inward due to interactions with a tenuous circumplanetary gas disk.We focus on the capture by Jupiter, which, due to its large mass, provides a stringent test of our model. We investigate the possible fates of disrupted bodies, the differences between prograde and retrograde captures, and the effects of Callisto on captured objects. We make an impulse approximation and discuss how it allows us to generalize capture results from equal-mass binaries to binaries with arbitrary mass ratios.We find that at Jupiter, binaries offer an increase of a factor of ∼10 in the capture rate of 100-km objects as compared to single bodies, for objects separated by tens of radii that approach the planet on relatively low-energy trajectories. These bodies are at risk of collision with Callisto, but may be preserved by gas drag if their pericenters are raised quickly enough. We conclude that our mechanism is as capable of producing large irregular satellites as previous suggestions, and it avoids several problems faced by alternative models.     

Formation of the regular satellites of giant planets in an extended gaseous nebula I: subnebula model and accretion of satellites

We model the subnebulae of Jupiter and Saturn wherein satellite accretion took place. We expect each giant planet subnebula to be composed of an optically thick (given gaseous opacity) inner region inside of the planet’s centrifugal radius (where the specific angular momentum of the collapsing giant planet gaseous envelope achieves centrifugal balance, located at rCJ ∼ 15R_J for Jupiter and rCS ∼ 22R_S for Saturn) and an optically thin, extended outer disk out to a fraction of the planet’s Roche-lobe (R_H), which we choose to be ∼R_H/5 (located at ∼150 R_J near the inner irregular satellites for Jupiter, and ∼200R_S near Phoebe for Saturn). This places Titan and Ganymede in the inner disk, Callisto and Iapetus in the outer disk, and Hyperion in the transition region. The inner disk is the leftover of the gas accreted by the protoplanet. The outer disk may result from the nebula gas flowing into the protoplanet during the time of giant planet gap-opening (or cessation of gas accretion). For the sake of specificity, we use a solar composition “minimum mass” model to constrain the gas densities of the inner and outer disks of Jupiter and Saturn (and also Uranus). Our model has Ganymede at a subnebula temperature of ∼250 K and Titan at ∼100 K. The outer disks of Jupiter and Saturn have constant temperatures of 130 and 90 K, respectively.Our model has Callisto forming in a time scale ∼10⁶ years, Iapetus in 10⁶-10⁷ years, Ganymede in 10³-10⁴ years, and Titan in 10⁴-10⁵ years. Callisto takes much longer to form than Ganymede because it draws materials from the extended, low density portion of the disk; its accretion time scale is set by the inward drift times of satellitesimals with sizes 300-500 km from distances ∼100R_J. This accretion history may be consistent with a partially differentiated Callisto with a ∼300-km clean ice outer shell overlying a mixed ice and rock-metal interior as suggested by Anderson et al. (2001), which may explain the Ganymede-Callisto dichotomy without resorting to fine-tuning poorly known model parameters. It is also possible that particulate matter coupled to the high specific angular momentum gas flowing through the gap after giant planet gap-opening, capture of heliocentric planetesimals by the extended gas disk, or ablation of planetesimals passing through the disk contributes to the solid content of the disk and lengthens the time scale for Callisto’s formation. Furthermore, this model has Hyperion forming just outside Saturn’s centrifugal radius, captured into resonance by proto-Titan in the presence of a strong gas density gradient as proposed by Lee and Peale (2000). While Titan may have taken significantly longer to form than Ganymede, it still formed fast enough that we would expect it to be fully differentiated. In this sense, it is more like Ganymede than like Callisto (Saturn’s analog of Callisto, we expect, is Iapetus). An alternative starved disk model whose satellite accretion time scale for all the regular satellites is set by the feeding of planetesimals or gas from the planet’s Roche-lobe after gap-opening is likely to imply a long accretion time scale for Titan with small quantities of NH₃ present, leading to a partially differentiated (Callisto-like) Titan. The Cassini mission may resolve this issue conclusively. We briefly discuss the retention of elements more volatile than H₂O as well as other issues that may help to test our model.     

Dust ring formation due to sublimation of dust grains drifting radially inward by the Poynting-Robertson drag: An analytical model

Dust particles exposed to the stellar radiation and wind drift radially inward by the Poynting-Robertson (P-R) drag and pile up at the zone where they begin to sublime substantially. The reason they pile up or form a ring is that their inward drifts due to the P-R drag are suppressed by stellar radiation pressure when the ratio of radiation pressure to stellar gravity on them increases during their sublimation phases. We present analytic solutions to the orbital and mass evolution of such subliming dust particles, and find their drift velocities at the pileup zone are almost independent of their initial semimajor axes and masses. We derive analytically an enhancement factor of the number density of the particles at the outer edge of the sublimation zone from the solutions. We show that the formula of the enhancement factor reproduces well numerical simulations in the previous studies. The enhancement factor for spherical dust particles of silicate and carbon extends from 3 to more than 20 at stellar luminosities L_?=0.8-500L_⊙, where L_⊙ is solar luminosity. Although the enhancement factor for fluffy dust particles is smaller than that for spherical particles, sublimating particles inevitably form a dust ring as long as their masses decrease faster than their surface areas during sublimation. The formulation is applicable to dust ring formation for arbitrary shape and material of dust in dust-debris disks as well as in the Solar System.     

The nature of the volcanic activity at Loki: Insights from Galileo NIMS and PPR data

Loki is the largest patera and the most energetic hotspot on Jupiter's moon Io, in turn the most volcanically active body in the Solar System, but the nature of the activity remains enigmatic. We present detailed analysis of Galileo Near-Infrared Mapping Spectrometer (NIMS) and PhotoPolarimeter/Radiometer (PPR) observations covering the 1.5-100 μm wavelength range during the I24, I27, and I32 flybys. The general pattern of activity during these flybys is consistent with previously proposed models of a resurfacing wave periodically crossing a silicate lava lake. In particular our analysis of the I32 NIMS observations shows, over much of the observed patera, surface temperatures and implied ages closely matching those expected for a wave advancing counterclockwise at 0.94-1.38 km/day. The age pattern is different than other published analyses which do not show as clearly this azimuthal pattern. Our analysis also shows two additional distinctly different patera surfaces. The first is located along the inner and outer margins where components with a 3.00-4.70-μm color temperature of 425 K exist. The second is located at the southwestern margin where components with a 550-K color temperature exist. Although the high temperatures could be caused by disruption of a lava lake crust, some additional mechanism is required to explain why the southwest margin is different from the inner or outer ones. Finally, analysis of the temperature profiles across the patera reveal a smoothness that is difficult to explain by simple lava cooling models. Paradoxically, at a subpixel level, wide temperature distributions exist which may be difficult to explain by just the presence of hot cracks in the lava crust. The resurfacing wave and lava cooling models explain well the overall characteristics of the observations. However, additional physical processes, perhaps involving heat transport by volatiles, are needed to explain the more subtle features.     

Interacting ellipsoids: a minimal model for the dynamics of rubble-pile bodies

A simple mechanical model is formulated to study the dynamics of rubble-pile asteroids, formed by the gravitational re-accumulation of fragments after the collisional breakup of a parent body. In this model, a rubble-pile consists of N interacting fragments represented by rigid ellipsoids, and the equations of motion explicitly incorporate the minimal degrees of freedom necessary to describe the attitude and rotational state of each fragment. In spite of its simplicity, our numerical examples indicate that the overall behavior of our model is in line with several known properties of collisional events, like the energy and angular momentum partition during high velocity impacts. Therefore, it may be considered as a well defined minimal model.     

Possible long-term decline in impact rates: 2. Lunar impact-melt data regarding impact history

Crater counts at lunar landing sites with measured ages establish a steep decline in cratering rate during the period ∼3.8 to ∼3.1 Gyr ago. Most models of the time dependence suggest a roughly constant impact rate (within factor ∼2) after about 3 Gyr ago, but are based on sparse data. Recent dating of impact melts from lunar meteorites, and Apollo glass spherules, clarifies impact rates from ∼3.2 to ∼2 Gyr ago or less. Taken together, these data suggest a decline with roughly 700 Myr half-life around 3 Gyr ago, and a slower decline after that, dropping by a factor ∼3 from about ∼2.3 Gyr ago until the present. Planetary cratering involved several phases with different time behaviors: (1) rapid sweep-up of most primordial planetesimals into planets in the first hundred Myr, (2) possible later effects of giant planet migration with enhanced cratering, (3) longer term sweep-up of leftover planetesimals, and finally (4) the present long-term “leakage” of asteroids from reservoirs such as the main asteroid belt and Kuiper belt. In addition, at any given point on the Moon, a pattern of “spikes” (sharp maxima of relatively narrow time width) will appear in the production rate of smaller craters (?500 m?), not only from secondary debris from large primary lunar impacts at various distances from the point in question, but also from asteroid breakups dotted through Solar System history. The pattern of spikes varies according to type of sample being measured (i.e., glass spherules vs impact melts). For example, several data sets show an impact rate spike ∼470 Myr ago associated with the asteroid belt collision that produced the L chondrites (see Section 3.6 below). Such spikes should be less prominent in the production record of craters of D? few km. These phenomena affect estimates of planetary surfaces ages from crater counts, as discussed in a companion paper [Quantin, C., Mangold, N., Hartmann, W.K., Allemand, P., 2007. Icarus 186, 1-10]. Fewer impact melts and glass spherules are found at ∼3.8 Gyr than at ∼3.5 Gyr ago, even though the impact rate itself is known to have been higher at 3.8 Gyr ago than 3.5 Gyr. This disproves the assertion by Ryder [Ryder, G., 1990. EOS 71, 313, 322-323] and Cohen et al. [Cohen, B.A., Swindle, T.D., Kring, D.A., 2000. Science 290, 1754-1756] that ancient impact melts are a direct proxy for ancient impact (cf. Section 3.3). This result raises questions about how to interpret cratering history before 3.8 Gyr ago.