Deuteron whistler and trans-equatorial propagation of the ion cyclotron whistler

New ion cyclotron whistlers which have the asymptotic frequency of one half the local proton gyrofrequency, $\text{G}{}_{\mathrm{p}}\text{2}$, and the minimum (or equatorial) proton gyrofrequency, G_pm, along the geomagnetic field line passing through the satellite have been found in the low-latitude topside ionosphere from the spectrum analysis of ISIS VLF electric field data received at Kashima, Japan. Ion cyclotron whistlers with asymptotic frequency of G_pm or $\text{G}{}_{\mathrm{pm}}\text{2}$ are observed only in the region of $\text{B}{}_{\mathrm{m}}\text{>}\text{B}\text{2}$ or rarely $\text{B}{}_{\mathrm{m}}\text{>}\text{B}\text{4}$, where B is the local magnetic field and B_m is the mini magnetic field along the geomagnetic field line passing through the satellite.The particles with one half the proton gyrofrequency may be the deuteron or alpha particle. Theoretical spectrograms of the electron whistlers (R-mode) and the ion cyclotron whistlers (L-mode) propagating along the geomagnetic field lines are computed for the appropriate distributions of the electron density and the ionic composition, and compared with the observed spectrograms.The result shows that the ion cyclotron whistler with the asymptotic frequency of $\text{G}{}_{\mathrm{p}}\text{2}$ is the deuteron whistler, and that the ion cyclotron whistlers with the asymptotic frequency of G_pm or $\text{G}{}_{\mathrm{pm}}\text{2}$ are caused by the trans-equatorial propagation of the proton or deuteron whistler from the other hemisphere.     

Comment on the group velocity of surface waves on the plasma sheet boundary

In the recent estimation by Maltsev and Lyatsky (1984) of the group velocity of surface waves on the inner boundary of the plasma sheet, the effect of the curvature of the field lines of the ambient magnetic field of the Earth on the spectrum has been assessed. The authors have not accounted for the fact, however, that the group velocity of the compressional surface magnetohydrodynamic waves itself is nonzero transverse to the magnetic field—a characteristic which has been omitted in the spectrum of Chen and Hasegawa (1974), being used by Maltsev and Lyatsky.This characteristic of compressional surface MHD waves is inherent for the spectrum $\text{ω = (}\text{k}{}_6\text{k}\text{)V}{}_{\mathrm{A}}\text{(k}{}^2{}_6\text{+ 2k}{}^2{}_{\mathrm{\perp }}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, obtained by Nenovski (1978) in the cold plasma limit V_A ? V_S(V_A is Alfvén velocity, and V_S, sound velocity). A comment has been made on the restrictions, proceeding from the approximation, used by Maltsev and Lyatsky. The estimation of the velocities for movements of auroral riometer absorption bays have been reviewed.     

The influence of ethane and acetylene upon the thermal structure of the Jovian atmosphere

Ethane and acetylene, both of which possess more efficient emission bands than methane, have been incorporated into a thermal structure model for the atmosphere of Jupiter. Choosing for illustrative purposes the mixing ratios $\text{[}\text{C}{}_2\text{H}{}_6\text{]}\text{[}\text{H}{}_2\text{]}\text{= 10}{}^{\mathrm{?5}}$ and $\text{[}\text{C}{}_2\text{H}{}_2\text{]}\text{[}\text{H}{}_2\text{]}\text{= 5 × 10}{}^{\mathrm{?7}}$, it is found that these hydrocarbon gases lower the atmospheric temperature within the thermal inversion region by as much as 20 K, subsequently reducing the emission intensity of the 7.7 μm CH₄ band below the observed result. It is qualitatively shown, however, that this cooling by C₂H₆ and C₂H₂ could be compensated by aerosol heating resulting from a uniformily mixed aerosol which absorbs 15% of the incident solar radiation. Such aerosol heating has been suggested by uv albedo observations.     

The role of protons of the radiation belts in the generation of Pc 3–5

A new theory of the Alfvén wave generation in inhomogeneous finite β two component plasma is developed ($\text{β =}\text{8πρ}\text{β}{}_0{}^2$, ρ and B₀ are plasma pressure and unperturbed magnetic field, respectively). The analysis was carried out for these waves both for long wave approximation kρ_i ? 1 as well as for $\text{kρ}{}_{\mathrm{i}}\text{? 1}$ (k and ρ_i are wave vector and larmor radius of protons). The influence of the loss-cone on the development of the instability is considered. The theory is applied to explain the generation mechanism of Pc 3–5.     

Love numbers of the giant planets

We calculate the Love numbers k_n for n = 2 to 10, and determine the “gravitational noise” from tides. The new values k₂ for Jupiter, Saturn, and Uranus yield new estimates for the planetary dissipation functions: $\text{Q}{}_{\mathrm{<mtext>J</mtext>}}\text{? 2.5 × 10}{}^4$, $\text{Q}{}_{\mathrm{<mtext>S</mtext>}}\text{? 1.4 × 10}{}^4\text{, Q}{}_{\mathrm{<mtext>U</mtext>}}\text{? 5 × 10}{}^3$.     

Analysis of HI content in galaxies

We have collected data on 241 galaxies from 13 sources and made a statistical analysis after reduction to a uniform system. We found that the Hubble sequence is one of increasing $\text{〈}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{M}{}_{\mathrm{<mtext>T</mtext>}}\text{〉}$ and $\text{〈}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{L}{}_{\mathrm{<mtext>B</mtext>}}\text{〉}$, these mean values increasing monotonically from .0016 and .024 at E to .084 and .83 at Im, but the dispersion is large.The HI content in barred spiral is greater than that in ordinary spirals, and this is consistent with their statistics of angular momentum and colour.The HI content is related to colour; it is greater in bluer systems. The large dispersion suggests that it also depends on some other factors, but these are smoothed out when averaged over each type, resulting in a linear relation between $\text{〈log(}\text{M}{}_{\mathrm{<mtext>M</mtext>}}\text{M}{}_{\mathrm{<mtext>T</mtext>}}\text{〉}$ and 〈(B ? V^O_T)〉. Unlike the colour-colour diagram, the large dispersion on the $\text{log}\text{(}\text{M}{}_{\mathrm{<mtext>H</mtext>}}\text{L}{}_{\mathrm{<mtext>B</mtext>}}\text{) ? (B ? V}{}^0{}_{\mathrm{<mtext>T</mtext>}}\text{)}$ is not related to peculiar galaxies.     

The photochemistry of hydrocarbons in the atmosphere of Titan

Detailed photochemical models are constructed for two model atmospheres: (1) 100% CH₄ and (2) 50% H₂, 50% CH₄. Both models predict large column densities of C₂H₂ and C₂H₆ (～1 cm atm) for eddy mixing rates ～10⁵ cm² sec^?, which are comparable to rates appropriate for Jupiter. These column densities vary inversely with the eddy diffusion coefficient. The models confirm the interpretation by Danielson et al. (1973) of the 12μ feature in the spectra of Gillet et al. (1973) as emission by C₂H₆ in a thermal inversion region. The $\text{C}{}_2\text{H}{}_6\text{C}{}_2\text{H}{}_2$ mixing ratio is sensitive to the net escape rate of H atoms from the exobase.     

Characteristics of localized resonance coupling oscillations of the slow magnetosonic wave in a non-uniform plasma

We analyze linear resonance oscillations in a non-uniform one-fluid finite-β plasma, which is oversimplified to understand easily fundamental characteristics of the resonance oscillations. A linear resonance oscillation of localized slow magnetosonic mode $\text{[ω}{}^2{}_{\mathrm{s}}\text{=}\text{ω}{}^2{}_{\mathrm{A}}\text{(1 +}\text{V}{}^2{}_{\mathrm{A}}\text{V}{}^2{}_{\mathrm{s}}\text{)}\text{]}$, which has the diamagnetic property in a uniform plasma, is newly found to be excited in the radially non-uniform plasma. The localized slow resonance indicates a radially polarized compressional oscillation (δB_∥ ? δB_H ? δB_D). The sense of the Alfvénic polarizations in the H-D plane near the resonant point is a function of both the propagation in the azimuthal direction and the slope of wave amplitude in the radial direction, whereas the sense of the resonant slow magnetosonic polarizations changes in accordance only with the switch in the azimuthal propagation direction. Further multi-satellite studies are necessary to establish the resonant structures of the slow magnetosonic waves in the magnetosphere.     

Low-frequency upstream wave as a probable source of low-latitude Pc 3–4 magnetic pulsations

Daytime Pc 3–4 pulsation activities observed at globally coordinated low-latitude stations [SGC (L = 1.8,λ = 118.0°W), EWA(1.15,158.1°W), ONW(1.3,141.5°E)] are evidently controlled by the cone angle θ_XB of the IMF observed at ISEE 3. Moreover, the Pc 3–4 frequencies ($\text{?}$) at the low latitudes and high latitude (COL; L = 5.6 and λ = 147.9°W) on the ground and that of compressional waves at geosynchronous orbit (GOES 2; L = 6.67 and λ = 106.7°W) are also correlated with the IMFmagnitude(B_IMF).The correlation of $\text{?}$ of the compressional Pc 3–4 waves at GOES 2 against B_IMF is higher than those of the Pc 3–4 pulsations at the globally coordinated ground stations, i.e., γ = 0.70 at GOES 2, and (0.36,0.60,0.66,0.54) at (COL, SGC, EWA, ONW), respectively. The standard deviation ($\text{σ}{}_{\mathrm{n}}\text{= ± Δ?}\text{mHz}$) of the observed frequencies from the form $\text{?}$ (mHz) = 6.0 × B_IMF (nT) is larger at the ground stations than at GOES 2, i.e., $\text{Δ? = ± 6.6}\text{mHz at}$GOES 2, and ±(13.9, 9.1, 10.7, 12.1) mHz at (COL, SGC, EWA, ONW), respectively. The correlations between the IMF magnitude B_IMF and Pc 3–4 frequencies at the low latitudes are higher than that at the high latitude on the ground, which can be interpreted by a “filtering action” of the magnetosphere for daytime Pc 3–4 magnetic pulsations. The scatter plots of pulsation frequency $\text{?}$ against the IMF magnitude B_IMF for the compressional Pc 3–4 waves at GOES 2 are restricted within the forms $\text{? = 4.5 × B}{}_{\mathrm{<mtext>IMF</mtext>}}\text{and}\text{? = 7.5 × B}{}_{\mathrm{<mtext>IMF</mtext>}}$. The frequency distribution is in excellent agreement with the speculation ($\text{<ω}{}_{\mathrm{sc}}\text{Ω}{}_{\mathrm{i}}\text{= 0.3 ～ 0.5}$) of the spacecraft frame frequency of the magnetosonic right-hand waves excited by the anomalous ion cyclotron resonance with reflected ion beams with V₆ = 650 ～ 1150 km s^?1 in the solar wind frame observed by the ISEE satellite in the Earth's foreshock. These observational results suggest that the magnetosonic right-handed waves excited by the reflected ion beams in the Earth's foreshock are convected through the magnetosheath to the magnetopause, transmitted into the magnetosphere without significant changes in spectra, and then couple with various HM waves in the Pc 3–4 frequency range at various locations in the magnetosphere.     

Venus: Cloud optical depth and surface albedo from Venera 8

We have used the measurements of the solar flux obtained by the Venera 8 spacecraft inside the atmosphere of Venus and the values of the Venus spherical albedo to deduce the characteristics of the clouds and of the ground. The method used is the exponential kernel approximation and the results have been tested by exact computations with the spherical harmonics method.A cloud layer with an optical thickness $\text{τ}{}_1\text{? 144}$, an albedo for single scattering $\text{ω}{}_0\text{= 0.9998}$ in the rear infrared, above a Rayleigh layer between 0 and 32 km and a ground of reflectivity ? = 0.4, gives a good agreement with the experimental results. A model with two cloud layers is also discussed.     

Convex profiles from asteroid lightcurves

Integral geometry is used to solve a two-dimensional simplification of the three-dimensional lightcurve inversion problem, and a method is introduced for obtaining a convex profile, $\text{P}$, from asteroid lightcurve data. Whenever four ideal conditions are satisfied, $\text{P}$ is an estimator for the asteroid's “mean cross section,” $\text{C}$, a convex set defined as the average of all cross sections $\text{C}$ cut by planes a distance z above the asteroid's equatorial plane. $\text{C}$ is therefore a two-dimensional average of the asteroid's three-dimensional shape. The ideal conditions are that (A) each curve $\text{C}\text{(z)}$ is convex. (B) the asteroid's scattering law is uniform and geometric, (C) the astrocentric declinations of the Sun and Earth are zero, and (D) the solar phase angle θ ≠ 0. If Condition C is known to hold, the extend to which the lightcurve can be accounted for by a geometrically scattering convex object can be quantified in terms of an appropriate “goodness-of-fit” static. If the solar phase angle is zero, as for radar “lightcurve,” then (i) method yields a profile $\text{P}{}_{\mathrm{<mtext>s</mtext>}}$ the symmetrization $\text{C}{}_{\mathrm{<mtext>s</mtext>}}$ of $\text{C}$; (ii) Condition A need not hold and if it does not, then the inversion yields the symmetrization of the asteroid's mean convex hull; and (iii) Fourier analysis of the lightcurve can reveal violation of Condition B. Doppler-frequency resolution of radar echoes at several rotational phases adds information by constraining the convex hull $\text{H}{}_{\mathrm{<mtext>p</mtext>}}$ of the asteroid's (not necessarily convex) polar silhouette. Estimation of a convex profile from a photoelectric or radar lightcurve is a problem in weighted-least-squares optimization subject to inequality constraints. The solution uses a recursive quadratic programming algorithm to derive a Fourier parameterization for $\text{P}$ from the coefficients in the lightcurve's Fourier expansion. The method has been tested by inverting analytically generated lightcurves for geometrically scattering ellipsoids with semiaxes a ? b ? c, and the inversion yields $\text{P}\text{=}\text{P}{}_{\mathrm{<mtext>s</mtext>}}\text{?}\text{C}\text{=}\text{C}{}_{\mathrm{<mtext>s</mtext>}}\text{=}\text{H}{}_{\mathrm{<mtext>p</mtext>}}$ when the viewing geometry (Condition C) is close to ideal. For situations when the asteroid's pole direction is unknown, a test is offered of the hyphothesis that a given lightcurve can be due to a geometrically scattering ellipsoid with $\text{a}\text{c}\text{? ?}$, where ? is an priori upper bound on the maximum axis ratio. Convex profiles are presented for 15 Eunomia, 118 Peitho, 246 Asporina, 281 Lucretia. 790 Pretoria, 1685 Toro, and 1978 CA.     

The Moon's zonal tidal effect in the variation of the earth's spin rate

Using the time observations obtained by 8 instruments in the Chinese Joint System during the years 1966–1980, we analyse the Moon's zonal tidal effect. The results show that the effects of the M_f and M_m waves are obvious. From this, the parameters $\text{K}\text{C}$ of the zonal tide are estimated and the weighted averages of the 8 instruments are $\text{(}\text{K}\text{C}\text{)}\text{M}{}_{\mathrm{f}}\text{= 0.909 ± 0.114}\text{and}\text{(}\text{K}\text{C}\text{)}\text{M}{}_{\mathrm{m}}\text{= 0.905 ± 0.083}$ respectively.     

Generation of magnetic-field aligned currents,parallel electric fields,and inverted-V structures by plasma pressure inhomogeneities in the magnetosphere

Magnetic-field aligned currents driven by plasma pressure inhomogeneities (plasma clouds) in the distant magnetosphere are analyzed quantitatively. A parallel potential drop is found to be established in the upward current region whenever a spatial scale D₀ for the pressure gradient in the equatorial magnetosphere is smaller than $\text{≈ 3g}{}_0\text{B}{}_{\mathrm{i}}\text{B}{}_0$, where g₀ is a hot electron gyroradius in the equatorial magnetic field B₀ (B_i denotes the magnetic induction in the ionosphere). A theoretical derivation is given for the experimentally observed linear relation T = AE_p + T₀ between the characteristic energy T of precipitating magnetospheric electrons and the peak energy E_p in inverted-V electron spectra. Three-dimensional potential structures accelerating electrons earthward are shown to be established beneath some model clouds which could correspond to a large scale inverted-V structure and to a thin (~ 1 km) auroral arc.     

An estimator of the underlying size distribution of overlapping impact-craters

The stochastic model of lunar type impact-crater formation which assumes (a) random impacts, (b) circular craters, each obliterating any portions of earlier craters lying within, and (c) a probability P_i(t) that a newly formed crater (primary or secondary) has an area a_i is analyzed to develop a method of estimating P_i from the final overlapping pattern. It is found that if each crater is weighted by the fraction of the rim which is visible and which lies in an observation area A, then the expected value of the weighted sum $\text{Ω}{}_{\mathrm{i}}$ of craters of area a_i is simply proportional to P_i for any degree of coverage under several conditions, including (a) constant P_i for all i, and (b) P_i stepping from a constant early value to zero (for some i's) with otherwise arbitrary bombardment. Furthermore, in the general case, the expected value of the contribution $\text{ΔΩ}{}_{\mathrm{i}}\text{(t}{}_0\text{)}\text{to}\text{Ω}{}_{\mathrm{i}}$ produced during t₀ ± Δt/2 is found to be proportional to P_i(t₀). Thus measurement of $\text{Ω}{}_{\mathrm{i}}$ in the first two cases, or of $\text{ΔΩ}{}_{\mathrm{i}}$ if crater age data is available in the last case, provides an estimate of the desired P_i. Therefore the $\text{Ω}{}_{\mathrm{i}}$ introduce the correct weighting factors that just compensate for the effect of overlap.Expressions for the variances of $\text{Ω}{}_{\mathrm{i}}\text{and}\text{Ω = Σ}{}_{\mathrm{i}}\text{Ω}{}_{\mathrm{i}}$ are derived from which it is shown that under the above conditions, $\text{Ω}{}_{\mathrm{i}}\text{/Ω}\text{or}\text{ΔΩ}{}_{\mathrm{i}}\text{/ΔΩ}$ are consistent estimators of P_i. Formal evaluation of the variances is carried out in the special case of constant P_i and no secondary cratering. A criterion for the degree of coverage is given; in particular it is shown that the expectation of $\text{σ = Σ}{}_{\mathrm{i}}\text{a}{}_{\mathrm{i}}\text{Ω}{}_{\mathrm{i}}$ at saturation is just A.     

Role of Coulomb collisions in the equatorial F-region plasma instabilities

The effect of collisions on electrostatic instabilities driven by gravity and density gradients perpendicular to the ambient magnetic field is studied. Electron collisions tend to stabilize the short wavelength (k_y?_i ? 1, where k_y is the perpendicular wavenumber of the instability and ?_i is the ion Larmor radius) kinetic interchange mode. In the presence of weak ion-ion collisions, this mode gets converted into an unmagnetized ion interchange mode which has maximum growth rate one order smaller than that of the collisionless mode. On the other hand, electron collisions can excite a long wavelength resistive interchange mode in a wide wavenumber regime ($\text{10}{}^{\mathrm{?}}\text{3 ? k}{}_{\mathrm{y}}\text{?}{}_{\mathrm{i}}\text{? 0.3}$) with growth rates comparable to that of the collisional Rayleigh-Taylor mode. The results may be relevant to some of the spread F irregularities.     

Theoretical interpretation of the ground-based photometry of Saturn's B ring

New photographic photometry at small tilt angles during the 1979 and 1981 apparitions is combined with earlier data to yield several physical parameters for Saturn's B ring in red and blue colors. Phase curves are obtained for a mean tilt angle B ? 6°. The value of the volume density D is 0.020±0.004 with no indication of dependence on either the color or the tilt angle for 6°<B<26°. This conclusion is not altered significantly if the individual ring particles have a phase function similar to the phase curves of bright solar system objects. For the geometric albedo of a single particle we derive 0.61±0.04 (red) and 0.41±0.03 (blue), which are superior to earlier estimates because of the additional data now available. These values and the derived amount of multiple scattering as a function of tilt angle constrain the particle phase function in the red to be moderately backscattering. Inferred values of the particle single-scattering albedo are $\text{0.7≤}\text{ω}{}_0\text{(}\text{red}\text{) ≤0.92}$ and $\text{0.5≤}\text{ω}{}_0\text{(}\text{blue}\text{) ≤0.7}$, depending on the choice of phase function. No indication was found that the particle photometric properties might depend on the vertical distance from the central plane. Our results show that the ground-based photometry is entirely consistent with the classical, many-particle-thick ring model.     

The relationship of the boundary layer of the magnetosphere to IPRP events

A model is proposed in which a mixture of hot solar wind and cold atmospheric plasma flowing in the dayside equatorial boundary layer towards the dawn-dusk plane generates hydromagnetic waves near the frequency $\text{ω = ω}{}_{\mathrm{Bi}}\text{¦1 ? T}{}_{\mathrm{\parallel }}\text{¦T}{}_{\mathrm{\perp }}\text{¦}$ where ω_Bi is the ion gyrofrequency and T_∥, T_⊥ are the temperatures of the solar wind plasma, parallel and perpendicular respectively to the magnetic field B. The model accounts for the properties of IPRP events, i.e. intervals of geomagnetic pulsations of periods rising on average from about 2 s to about 7 s over an interval of about 5 min. The diagnostic potential of this phenomenon for study of the boundary layer is indicated.     

Photometric properties of zodiacal light particles

From published ground-base, spacecraft, and rocket photometry and polarimetry of the zodiacal light, a number of optical and physical parameters have been derived. It was assumed that the number density, mean particle size, and albedo vary with heliocentric distance, and shown that average individual interplanetary particles have a small but definite opposition effect, a mean single-scattering albedo in the V band at 1-AU heliocentric distance of 0.09 ± 0.01, and a zero-phase geometric albedo of 0.04. Modeled by a power law, both albedos decrease with increasing heliocentric distance as r^?0.54. The corresponding exponents for changes in mean particle size and number density are related in a simple way. The median orbital inclination of zodiacal light particles with respect to the ecliptic is 12°, close to the observed median value for faint asteroids and short-period comets. Furthermore, the color of dust particles and its variation with solar phase angle closely resemble those of C asteroids. These findings are, at least, consistent with the zodiacal cloud originating primarily from collisions among asteroids. Finally, a value of ?10¹⁸?_ErmE g was derived for the mass of the zodiacal cloud, where ?_E is the mean particle radius (in micrometers) at 1-AU-heliocentric distance. For extinction in the ecliptic, $\text{Δm = 10}{}^{\mathrm{?5}}\text{?}{}^{\mathrm{<mtext>?1</mtext><mtext>2</mtext>}}\text{mag}$ was obtained, where ? is the solar elongation in degrees.     

Editorial

The Galilean satellites Io, Europa, and Ganymede interact through several stable orbital resonances where $\text{λ}{}_1\text{? 2λ}{}_2\text{+}\text{ω}{}_1\text{= 0, λ}{}_1\text{? 2λ}{}_2\text{+}\text{ω}{}_2\text{= 180°, λ}{}_2\text{? 2λ}{}_3\text{+}\text{ω}{}_2\text{= 0}$ and λ₁ ? 3λ₂ + 2λ₃ = 180°, with λ_i being the mean longitude of the ith satellite and $\text{ω}{}_{\mathrm{i}}$ the longitude of the pericenter. The last relation involving all three bodies is known as the Laplace relation. A theory of origin and subsequent evolution of these resonances outlined earlier (C. F. Yoder, 1979b, Nature279, 747–770) is described in detail. From an initially quasi-random distribution of the orbits the resonances are assembled through differential tidal expansion of the orbits. Io is driven out most rapidly and the first two resonance variables above are captured into libration about 0 and 180° respectively with unit probability. The orbits of Io and Europa expand together maintaining the 2:1 orbital commensurability and Europa's mean angular velocity approaches a value which is twice that of Ganymede. The third resonance variable and simultaneously the Laplace angle are captured into libration with probability ～0.9. The tidal dissipation in Io is vital for the rapid damping of the libration amplitudes and for the establishment of a quasi-stationary orbital configuration. Here the eccentricity of Io's orbit is determined by a balance between the effects of tidal dissipation in Io and that in Jupiter, and its measured value leads to the relation $\text{k}{}_1\text{?}{}_1\text{/Q}{}_1\text{≈ 900k}{}_{\mathrm{<mtext>J</mtext>}}\text{/Q}{}_{\mathrm{<mtext>J</mtext>}}$ with the k's being Love numbers, the Q's dissipation factors, and f a factor to account for a molten core in Io. This relation and an upper bound on Q₁ deduced from Io's observed thermal activity establishes the bounds 6 × 10⁴ < Q_J < 2 × 10⁶, where the lower bound follows from the limited expansion of the satellite orbits. The damping time for the Laplace libration and therefore a minimum lifetime of the resonance is 1600 Q_J years. Passage of the system through nearby three-body resonances excites free eccentricities. The remnant free eccentricity of Europa leads to the relation $\text{Q}{}_2\text{/?}{}_2\text{? 2 × 10}{}^{\mathrm{?4}}\text{Q}{}_{\mathrm{<mtext>J</mtext>}}$ for rigidity μ₂ = 5 × 10¹¹ dynes/cm². Probable capture into any of several stable 3:1 two-body resonances implies that the ratio of the orbital mean motions of any adjacent pair of satellites was never this large.A generalized Hamiltonian theory of the resonances in which third-order terms in eccentricity are retained is developed to evaluate the hypothesis that the resonances were of primordial origin. The Laplace relation is unstable for values of Io's eccentricity e₁ > 0.012 showing that the theory which retains only the linear terms in e₁ is not valid for values of e₁ larger than about twice the current value. Processes by which the resonances can be established at the time of satellite formation are undefined, but even if primordial formation is conjectured, the bounds established above for Q_J cannot be relaxed. Electromagnetic torques on Io are also not sufficient to relax the bounds on Q_J. Some ideas on processes for the dissipation of ideal energy in Jupiter yield values of Q_J within the dynamical bounds, but no theory has produced a Q_J small enough to be compatible with the measurements of heat flow from Io given the above relation between Q₁ and Q_J. Tentative observational bounds on the secular acceleration of Io's mean motion are also shown not to be consistent with such low values of Q_J. Io's heat flow may therefore be episodic. Q_J may actually be determined from improved analysis of 300 years of eclipse data.