Ion and neutral sources and sinks within Saturn's inner magnetosphere: Cassini results

Using ion-electron fluid parameters derived from Cassini Plasma Spectrometer (CAPS) observations within Saturn's inner magnetosphere as presented in Sittler et al. [2006a. Cassini observations of Saturn's inner plasmasphere: Saturn orbit insertion results. Planet. Space Sci., 54, 1197-1210], one can estimate the ion total flux tube content, N_IONL², for protons, H⁺, and water group ions, W⁺, as a function of radial distance or dipole L shell. In Sittler et al. [2005. Preliminary results on Saturn's inner plasmasphere as observed by Cassini: comparison with Voyager. Geophys. Res. Lett. 32(14), L14S04), it was shown that protons and water group ions dominated the plasmasphere composition. Using the ion-electron fluid parameters as boundary condition for each L shell traversed by the Cassini spacecraft, we self-consistently solve for the ambipolar electric field and the ion distribution along each of those field lines. Temperature anisotropies from Voyager plasma observations are used with (T_⊥/T_∥)_W⁺∼5 and (T_⊥/T_∥)_H⁺∼2. The radio and plasma wave science (RPWS) electron density observations from previous publications are used to indirectly confirm usage of the above temperature anisotropies for water group ions and protons. In the case of electrons we assume they are isotropic due to their short scattering time scales. When the above is done, our calculation show N_IONL² for H⁺ and W⁺ peaking near Dione's L shell with values similar to that found from Voyager plasma observations. We are able to show that water molecules are the dominant source of ions within Saturn's inner magnetosphere. We estimate the ion production rate S_ION∼10²⁷ ions/s as function of dipole L using N_H⁺, N_W⁺ and the time scale for ion loss due to radial transport τ_D and ion-electron recombination τ_REC. The ion production shows localized peaks near the L shells of Tethys, Dione and Rhea, but not Enceladus. We then estimate the neutral production rate, S_W, from our ion production rate, S_ION, and the time scale for loss of neutrals by ionization, τ_ION, and charge exchange, τ_CH. The estimated source rate for water molecules shows a pronounced peak near Enceladus’ L shell L∼4, with a value S_W∼2×10²⁸ mol/s.     

Derivation of density and temperature from the Cassini–Huygens CAPS electron spectrometer

In this paper we present two methods to derive electron fluid parameters from the CAPS–ELS spectrometer on board the Cassini spacecraft currently in orbit around Saturn. In the first part of the paper we give a basic overview of the instrument and describe the challenges inherent in the derivation of density and temperature values using these techniques. We then describe a method to calculate electron moments by integrating the particle distribution function. We also describe a second technique in which we fit the electron energy spectrum with a Gaussian curve and use the peak energy of this curve to derive density and temperature values. We then compare the two methods with particular emphasis on their application to Cassini SOI observations in the saturnian environment and point out the limitations of the two techniques. We will show that results from the two very different methods are in agreement when the physical properties of the environment and of the observed electron populations have been inferred from inspection of the raw data. Finally we will suggest future developments that will remove these limitations.     

Adiabatic chaos in the Prometheus–Pandora system

The chaotic orbital motion of Prometheus and Pandora, the 16th and 17th satellites of Saturn, is studied. Chaos in their orbital motion, as found by Goldreich & Rappaport and Renner & Sicardy, is due to interaction of resonances in the resonance multiplet corresponding to the 121:118 commensurability of the mean motions of the satellites. It is shown rigorously that the system moves in adiabatic regime. The Lyapunov time (the 'time horizon of predictability' of the motion) is calculated analytically and compared to the available numerical–experimental estimates. For this purpose, a method of analytical estimation of the maximum Lyapunov exponent in the perturbed pendulum model of non-linear resonance is applied. The method is based on the separatrix map theory. An analytical estimate of the width of the chaotic layer is made as well, based on the same theory. The ranges of chaotic diffusion in the mean motion are shown to be almost twice as big compared to previous estimates for both satellites.     

Core formation in giant gaseous protoplanets

Sedimentation rates of silicate grains in gas giant protoplanets formed by disk instability are calculated for protoplanetary masses between 1 MSaturn to 10 MJupiter. Giant protoplanets with masses of 5 MJupiter or larger are found to be too hot for grain sedimentation to form a silicate core. Smaller protoplanets are cold enough to allow grain settling and core formation. Grain sedimentation and core formation occur in the low mass protoplanets because of their slow contraction rate and low internal temperature. It is predicted that massive giant planets will not have cores, while smaller planets will have small rocky cores whose masses depend on the planetary mass, the amount of solids within the body, and the disk environment. The protoplanets are found to be too hot to allow the existence of icy grains, and therefore the cores are predicted not to contain any ices. It is suggested that the atmospheres of low mass giant planets are depleted in refractory elements compared with the atmospheres of more massive planets. These predictions provide a test of the disk instability model of gas giant planet formation. The core masses of Jupiter and Saturn were found to be ∼0.25 M_⊕ and ∼0.5 M_⊕, respectively. The core masses of Jupiter and Saturn can be substantially larger if planetesimal accretion is included. The final core mass will depend on planetesimal size, the time at which planetesimals are formed, and the size distribution of the material added to the protoplanet. Jupiter's core mass can vary from 2 to 12 M_⊕. Saturn's core mass is found to be ∼8 M_⊕.     

An Experimental Evidence Of Non-Ducted Propagation During Daytime From The Esd Whistler Observations At Low Latitude Ground Station Jammu

The propagationmechanism of low latitude daytime whistlers is investigated on the basis of ground measurements made continuously during daytime in North India at Jammu (geomag. lat. 22°26°N;L = 1.17). On February 14, 1998 extremely small dispersion (ESD) whistlers with dispersion varying from 5–10 sec^1/2 in surprisingly large numbers were recorded at Jammu during daytime in the late afternoon. The results of a study of the characteristics of ESD whistlers are presented and the discussion indicates that ESD whistlers recorded are the VLF waves radiated from the return stroke of the lightning discharge launched at the ionosphere with different initial wave normal angles, propagated upwards under eitherquasi-longitudinal conditions or pro-longitudinal whistler mode, turned around at different heights due to quasi-transverse propagation and received at Jammu with the dispersion of the order of 5–10 sec^1/2. The validity of this suggestion has been tested by performing actualray-tracing computations in thepresence of equatorial anomaly model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Pc3 pulsations during variable IMF conditions

Pc3 geomagnetic field fluctuations detected at low latitude (L’Aquila, Italy) during the passage of a high velocity solar wind stream, characterized by variable interplanetary magnetic field conditions, are analyzed. Higher frequency resonant fluctuations and lower frequency phenomena are simultaneously observed; the intermittent appearance and the variable frequency of the longer period modes can be well interpreted in terms of the variable IMF elements; moreover their polarization characteristics are consistent with an origin related to external waves propagating in antisunward direction. A comparison with simultaneous observations performed at Terra Nova Bay (Antarctica) provides additional evidence for a clear relationship between the IMF and Pc3 pulsations also at very high latitudes.     

1989年3月特强磁暴期间的电离层暴

用欧亚大陆地面电离层垂测站资料考察1989年3月12～16日磁暴期间的电离层暴形态及其发展变化. 特强磁暴引发的电离层暴是全球性的,但自磁层沉降的高能粒子对热层低部的加热程度及区域分布不同,因而各经度链区域内电离层暴的特征也有所差异. 本文研究表明,与理论推断对照,欧洲地区内F2层最大电子密度NmaxF2(或f0F2)并不出现正暴现象,而负暴自高纬向低纬的发展则与典型的热层环流结果相符. 此外,此磁暴过程期间在中低纬区存在明显的波动过程. 在亚洲高纬地区,磁暴初期13日有约10 h的正暴,而负暴过程则与欧洲地区类似,但不太清晰;且无波动现象. 磁暴期间,同一经度链的中低纬地区,夜间常发生多站同时的h′F突增. 本文再次证实这是一般磁暴期间常出现的普遍现象.     

Ground-based radar observations of Titan: 2000-2008

We have observed Titan with the Arecibo Observatory’s 12.6 cm wavelength radar system during the last eight oppositions of the Saturn system with sufficient sensitivity to characterize its scattering properties as a function of sub-Earth longitude. In a few sessions the Green Bank Telescope was used as the receiving instrument in a bistatic configuration to boost sub-radar track length and integration time. Radar echo spectra have been obtained for a total of 92 viewing geometries with sub-Earth locations scattered through all longitudes and at latitudes between 7.6°S and 26.3°S, close to the maximum southern excursion of the sub-Earth track. We find Titan to have globally average radar albedos at this wavelength of 0.161 in the opposite circular polarization sense as that transmitted (OC) and 0.074 in the same sense (SC), giving a polarization ratio SC/OC of 0.46. These values are intermediate between lower reflectivity rocky surfaces and higher reflectivity clean icy surfaces. The variations with longitude in general mirror the surface brightness variations seen through the infrared atmospheric windows. Xanadu Regio’s radar reflectivity and polarization ratio are higher than the global averages, and suggest that its composition is relatively cleaner water ice or, possibly, some other material with low propagation loss at radio wavelengths. For all echo spectra most of the power is in a broad diffuse component but with a specular component whose strength and narrowness is highly variable as a function of surface location. For all data we fit a sum of the standard Hagfors scattering law describing the specular component and an empirical diffuse radar scattering model to extract bulk parameters of the surface. Many areas exhibit very narrow specular reflections implying terrain that are quite flat on centimeter to meter scales over spans of tens to perhaps hundreds of kilometers. The proportion of spectra showing these narrow specular echoes has fallen significantly over the observational time span, indicating either a latitudinal effect related to terrain differences or changing surface conditions over the past several years. A few radar tracks, especially those from the 2008 session, overlap some high resolution Cassini RADAR imagery swaths to allow a direct comparison with terrain.     

Regolith depth growth on an icy body orbiting Saturn and evolution of bidirectional reflectance due to surface composition changes

Using a Markov chain model, we consider the regolith growth on a small body in orbit around Saturn, subject to meteoritic bombardment, and assuming all impact ejecta are re-collected. We calculate the growth of regolith and the fractional pollution, assuming an initial pure ice body and amorphous carbon as a pollutant. We extend the meteorite flux of Cuzzi and Estrada (Cuzzi, J., Estrada, P. [1998]. Icarus 132, 1-35) to larger sizes to consider the effect of disruption of the moonlet on other moonlets in the ensemble. This is a relatively small effect, completely negligible for moonlets of 1 m radius. For the given impact model, fractional pollution reaches 22% for 1 m bodies, but only 3% for 10 m bodies, 1.7% for 20 m bodies, and 1% for 30 m bodies after 4 byr. By considering an ensemble of moonlets, which have identical cross-sections for releasing and capturing ejecta, this analysis can be extended to a model of particles in Saturn’s rings, where the calculated spectra can be compared to observed ring spectra. The measured spectral reflectance of Saturn’s rings from Cassini observations therefore constrains the size and age of the ring particles. The comparison between 1 m, 10 m, 20 m, and 30 m particles confirms that for larger ring mass, the current rings would be less polluted; for the largest particles, we expect negligible changes in the UV spectrum after 4 byr of meteoritic bombardment. We consider two end members for mixing of the meteoritic material: areal and intimate. Given the uncertainties in the actual mixing of the meteoritic infall and in its composition (as a worst case, we assume the meteoritic material is 100% amorphous carbon, intimately mixed) initially pure ice 30 m ring particles would darken after 4 byr of exposure by 15%.     

Dipolarization fronts in the magnetotail plasma sheet

We present a THEMIS study of a dipolarization front associated with a bursty bulk flow (BBF) that was observed in the central plasma sheet sequentially at X=−20.1, −16.7, and −11.0R_E. Simultaneously, the THEMIS ground network observed the formation of a north-south auroral form and intensification of westward auroral zone currents. Timing of the signatures in space suggests earthward propagation of the front at a velocity of 300 km/s. Spatial profiles of current and electron density on the front reveal a spatial scale of 500 km, comparable to an ion inertial length and an ion thermal gyroradius. This kinetic-scale structure traveled a macroscale distance of 10R_E in about 4 min without loss of coherence. The dipolarization front, therefore, is an example of space plasma cross-scale coupling. THEMIS observations at different geocentric distances are similar to recent particle-in-cell simulations demonstrating the appearance of dipolarization fronts on the leading edge of plasma fast flows in the vicinity of a reconnection site. Dipolarization fronts, therefore, may be interpreted as remote signatures of transient reconnection.