Titan: Far-infrared and microwave remote sensing of methane clouds and organic haze

It is shown that Titan's surface and plausible atmospheric thermal opacity sources—gaseous N₂, CH₄, and H₂, CH₄ cloud, and organic haze—are sufficient to match available Earth-based and Voyager observations of Titan's thermal emission spectrum. Dominant sources of thermal emission are the surface for wavelenghts λ ? 1 cm, atmospheric N₂ for 1 cm ? λ ? 200 μm,, condensed and gaseous CH₄ for 200 μm ? λ ? 20 μm, and molecular bands and organic haze for λ ? 20 μm. Matching computed spectra to the observed Voyager IRIS spectra at 7.3 and 52.7° emission angles yields the following abundances and locations of opacity sources: CH₄ clouds: 0.1 g cm^? at a planetocentric radius of 2610–2625 km, 0.3 g cm^?2 at 2590–2610 km, total 0.4 ± 0.1 g cm^–2 above 2590 km; organic haze: $\text{4 ± 2 × 10}{}^{\mathrm{?6}}\text{,}\text{g cm}\text{,}{}^{\mathrm{?2}}$ above 2750 km; tropospheric H₂: 0.3 ± 0.1 mol%. This is the first quantitative estimate of the column density of condensed methane (or CH₄/C₂H₆) on Titan. Maximum transparency in the middle to far IR occurs at 19 μm where the atmospheric vertical absorption optical depth is $\text{?0.6}$ A particle radius $\text{r}\text{? 2 μ}\text{m}$ in the upper portion of the CH₄ cloud is indicated by the apparent absence of scattering effects.     

The 63 μm radiation field in the earth's thermosphere and its influence on the atomic hydrogen temperature

The thermal escape of hydrogen from the Earth's atmosphere is strongly affected by its temperature at the exobase. It has been suggested recently that the hydrogen temperature might be significantly lower than the thermospheric temperature as a result of a collisional exchange of energy with atomic oxygen. The tendency is to cool the hydrogen since the energy of the excited ${}^3\text{P}{}_1$ level of oxygen can be lost from the atmosphere via magnetic dipole emission of the 63 μm line (${}^3\text{P}{}_2\text{?}{}^3\text{P}{}_1$). We present here a detailed calculation of the net cooling effect as a function of altitude throughout the thermosphere. The calculations have been performed for both day and night conditions and for periods of maximum and minimum solar activity conditions. It is found that its effect on ΔT/T varies from a very small value to a maximum of ～3%. We also provide the theoretical framework for describing deviations of the 63 μm emission from local thermodynamic equilibrium and show that these effects can cause the emission to be reduced by as much as 40% near 500 km.     

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.     

Absolute spectrophotometry of Neptune: 3390 to 7800 Å

Absolute spectrophotometry of Neptune from 3390 to 7800 Å, with spectral resolution of 10 Å in the interval 3390–6055 and 20 Å in the interval 6055–7800 Å, is reported. The results are compared with filter photometry (Appleby, 1973; Wamsteker, 1973; Savage et al., 1980) and with synthetic spectra computed on the basis of a parameterization proposed by Podolak and Danielson (1977) for aerosol scattering and absorption. A CH₄/H₂ ratio of is derived for the convectively mixed part of Neptune's atmosphere, and constrains optical properties of hypothetical aerosol layers.     

Radiation forces on small particles in the solar system

We present a new and more accurate expression for the radiation pressure and Poynting-Robertson drag forces; it is more complete than previous ones, which considered only perfectly absorbing particles or artificial scattering laws. Using a simple heuristic derivation, the equation of motion for a particle of mass m and geometrical cross section A, moving with velocity v through a radiation field of energy flux density S, is found to be (to terms of order $\text{v}\text{c}$)

$\text{m}\text{v}\text{?}\text{= (}\text{SA}\text{c}\text{)Q}{}_{\mathrm{<mtext>pr</mtext>}}\text{[(1 ?}\text{r}\text{?}\text{c}\text{)}\text{S}\text{?}\text{?}\text{v}\text{c}\text{]}$

, where ? is a unit vector in the direction of the incident radiation, $\text{r}\text{?}$ is the particle's radial velocity, and c is the speed of light; the radiation pressure efficiency factor Q_pr ≡ Q_abs + Q_sca(1 ? 〈cos α〉), where Q_abs and Q_sca are the efficiency factors for absorption and scattering, and 〈cos α〉 accounts for the asymmetry of the scattered radiation. This result is confirmed by a new formal derivation applying special relativistic transformations for the incoming and outgoing energy and momentum as seen in the particle and solar frames of reference. Q_pr is evaluated from Mie theory for small spherical particles with measured optical properties, irradiated by the actual solar spectrum. Of the eight materials studied, only for iron, magnetite , and graphite grains does the radiation pressure force exceed gravity and then just for sizes around 0.1 μm; very small particles are not easily blown out of the solar system nor are they rapidly dragged into the Sun by the Poynting-Robertson effect. The solar wind counterpart of the Poynting-Robertson drag may be effective, however, for these particles. The orbital consequences of these radiation forces-including ejection from the solar system by relatively small radiation pressures-and of the Poynting-Robertson drag are considered both for heliocentric and planetocentric orbiting particles. We discuss the coupling between the dynamics of particles and their sizes (which diminish due to sputtering and sublimation). A qualitative derivation is given for the differential Doppler effect, which occurs because the light received by an orbiting particle is slightly red-shifted by the solar rotation velocity when coming from the eastern hemisphere of the Sun but blue-shifted when from the western hemisphere; the ratio of this force to the Poynting-Robertson force is $\text{(}\text{R}{}_{\mathrm{\odot }}\text{r}\text{)}{}^2\text{[(}\text{w}{}_{\mathrm{\odot }}\text{n}\text{) ? 1]}$, where R_⊙ and w_⊙ are the solar radius and spin rate, and n is the particle's mean motion. The Yarkovsky effect, caused by the asymmetry in the reradiated thermal emission of a rotating body, is also developed relying on new physical arguments. Throughout the paper, representative calculations use the physical and orbital properties of interplanetary dust, as known from various recent measurements.     

Total current to cylindrical collectors in collisionless plasma flow

A theory is presented for charged-particle collection by a cylindrical conducting object, such as a spacecraft or an electrostatic probe, which is moving transversely through a collisionless plasma, such as those in the upper atmosphere and space. The calculation is approximate, using symmetric potential profiles which are exact for the infinite-cylinder stationary case. Theoretical current predictions are presented for ratios of collector potential to electron thermal energy eφ_c/kT_e from 0 to ?25, for ion-to-electron temperature ratios T_i/T_c = 1 and 0.5, ratio of collector radius to electron Debye length r_c/λ_D from 0 to 100, and ratio of flow speed to ion thermal speed $\text{S}{}_{\mathrm{i}}\text{= U/(2kT}{}_{\mathrm{i}}\text{/m}{}_{\mathrm{i}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{)}$ from 0 to 10. Comparisons with existing exact calculations by other authors show that none of these fulfil all of the requirements for nontrivial comparison. Appropriate parameter ranges for future exact calculations are thereby suggested. These are as follows: (a) r_c/λ_D should be large enough that the collector not be in or near orbit-limited conditions; (b) the ratio $\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,\; i}}\text{¦}$ of ion directed energy to potential energy change in the sheath, should be close to unity or if

$\text{S}{}_{\mathrm{i}}{}^2\text{/¦χ}{}_{\mathrm{c,i}}\text{¦? 1,}\text{then}\text{S}{}_{\mathrm{i}}\text{? 1}$

.     

Small-scale turbulence in the atmosphere of Venus

Previous studies based on radio scintillation measurements of the atmosphere of Venus have identified two regions of small-scale temperature fluctuations located in the vicinity of 45 and 60 km. A global study of the fluctuations near 60 km, which are consistent with wind-shear-generated turbulence, was conducted using the Pioneer Venus measurements. The structure constants of refractive index fluctuations c_n² and temperature fluctuations c_T² increase poleward, peak near 70° latitude, and decrease over the pole; c_n² varies from 2 × 10^?15 to 1.5 × 10^?14m^{$\text{2}\text{3}$} and c_T² from 4 × 10^?3 to 7 × 10^?2°K²m^{$\text{?2}\text{3}$}. These results indicate greater turbulent activity at the higher latitudes. In the region near 45 km the refractive index fluctuations and the corresponding temperature fluctuations are substantially lower. Based on the analysis of one representative occultation measurement, $\text{c}{}_{\mathrm{<mtext>n</mtext>}}{}^2\text{= 2 × 10}{}^{\mathrm{?16}}\text{m}{}^{\mathrm{<mtext>?2</mtext><mtext>3</mtext>}}\text{and}\text{c}{}_{\mathrm{<mtext>T</mtext>}}{}^2\text{= 7.3 × 10}{}^{\mathrm{?4}}\text{°}\text{K}{}^2\text{m}{}^{\mathrm{<mtext>?2</mtext><mtext>3</mtext>}}$ in the 45-km region. The fluctuations in this region also appear to be consistent with wind-shear-generated turbulence. The turbulence level is considerably weaker than that at 60 km; the energy dissipation rate ε is 4.9 × 10^?5m²sec^?3 and the small-scale eddy diffusion coefficient K is 2 × 10³ cm² sec^?1.     

Habitable zones about main sequence stars

Calculations show that a main sequence star which is less massive than the Sun has a continuously habitable zone about it which is not only closer in than the corresponding zone about the Sun, but is also relatively narrower. Let L(t) represent the luminosity after t billion years of a main sequence star of mass M, and let r_inner and r_outer represent the boundaries of the continuously habitable zone about such a star—that is, the zone in which an Earthlike planet will undergo neither a runaway greenhouse effect in the early stages of its history nor runaway glaciation after it develops an oxidizing atmosphere. Then our computer results indicate that $\text{r}{}_{\mathrm{<mtext>outer</mtext>}}\text{r}{}_{\mathrm{<mtext>inner</mtext>}}$ is roughly proportional to $\text{[}\text{L(3.5)}\text{L(1.0)}\text{]}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$. This ratio is smaller for stars less massive than the Sun (because they evolve more slowly), and the width of the continuously habitable zone about a main sequence star is therefore a strong function of the initial stellar mass. Our calculations show that r_inner = r_outer for $\text{M～0.83M?}$ (i.e., K1 stars), and it therefore appears that there is no continuously habitable zone about most K stars, nor any about M stars.     

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.     

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.     

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.     

Residual mass from atmospheric ablation of small meteoroids

Numerical solutions of the equations of meteor ablation in the Earth's atmosphere have been obtained using a variable step size Runge-Kutta technique in order to determine the size of the residual mass resulting from atmospheric flight. The equations used include effects of meteoroid heat capacity and thermal radiation, and a realistic atmospheric density profile. Results were obtained for initial masses in the range 10^?7–10^?2 g, and for initial velocities less than 24 km s^?1 (results indicated no appreciable residual mass for meteors with velocities above 24 km s^?1 in this mass range). The following function has been obtained to provide the logarithm of the ratio of the residual mass following atmospheric ablation to the original preatmospheric mass

$\text{log r = 4.7 ?0.33v}{}_{\mathrm{\infty }}\text{?0.013v}{}_{\mathrm{\infty }}{}^2\text{+ 1.2 log m}{}_{\mathrm{\infty }}\text{+ 0.08 log}{}^2\text{m}{}_{\mathrm{\infty }}\text{?0.083v}{}_{\mathrm{\infty }}\text{log m}{}_{\mathrm{\infty }}\text{M}$

The pre-atmospheric mass and velocity are represented by m_∞ and v_∞.When the results are expressed in terms of the size of the residual mass following atmospheric ablation as a function of the initial mass and velocity, it is found that the final residual mass is almost independent of the original mass of the meteoroid, but very strongly dependent on the original velocity. For example, the residual mass is very nearly 10^?7 g for a meteoroid with velocity 18 kms^?1 for initial masses from 10^?7 to 10^?3 g. On the other hand, a slight change in the initial velocity to 20 km s^?1 will shift the residual mass to approx. 10^?8 g. This strong velocity dependence coupled with the weak dependence on the original mass has important consequences for the sampling of ablation product micrometeorites.     

Electric fields and potential patterns in the high-latitude ionosphere for different situations in interplanetary space

The potential ? of the electric field at high latitudes has been obtained by solving numerically the second order differential equation in spherical coordinates:

$\text{?}\text{1}\text{2}\text{(rσ}{}^{\mathrm{H}}\text{?}{}_{\mathrm{\theta }}\text{)}{}_{\mathrm{\theta }}\text{+}\text{1}\text{r}\text{(σ}{}^{\mathrm{H}}\text{?λ)}{}_{\mathrm{\lambda }}\text{+}\text{1}\text{r}\text{(σ}{}^{\mathrm{P}}\text{?}{}_{\mathrm{\lambda }}\text{)}{}_{\mathrm{\theta }}\text{?(σ}{}^{\mathrm{P}}\text{?}{}_{\mathrm{\theta }}\text{)}{}_{\mathrm{\lambda }}\text{=}\text{1}\text{r}\text{(rψ}{}_{\mathrm{\theta }}\text{)}{}_{\mathrm{\theta }}\text{+}\text{1}\text{r}{}^2\text{ψ}{}_{\mathrm{\lambda \lambda }}$

, where θ is colatitude, λ is longitude, σ^H and σ^P are the height-integrated Hall and Perdersen ionospheric conductivities, r = sinθ, and ψ is the current function. The boundary condition is ? = 0 on the geomagnetic parallel θ = 34°. Values of ψ are determined from geomagnetic field variations at the Earth's surface from geomagnetic field variations at the Earth's surface for various conditions in interplanetary space. σ^P and σ^H are taken to vary with season, local time, tilt of the geomagnetic dipole axis (UT), and intensity of corpuscular precipitation (the model proposed by Wallis and Budzinski, 1981). The model distributions of ?^M and E^M = -▽?^m so obtained are compared with observational results. The feasibility has been demonstrated of interpreting the statistical results and individual measurement data in terms of a unified dynamic model of ionospheric electric fields. The model makes allowance for the changes of electromagnetic “weather” in interplanetary space.     

Possibility of the investigation of interplanetary shock interaction with the earth's bow shock,magnetopause and plasmapause by means of a non-coherent response method

The problem of interaction between the interplanetary shock of 8 March, 1970 and the Earth's bow shock, magnetopause and plasmapause is considered. Estimates are made using existing models of the moments of initial impulsive interaction of interplanetary shocks with the bow shock and of the secondary interaction of the resulting split discontinuities with the magnetopause, plasmapause and a modified bow shock. Using computed data on the plasma's concentration jumps at discontinuities and on the latters' velocities, estimates have been carried out of remote sounding and the response signals' phase difference change rates Δf (which were found to be of the order of ～ 10^?3?10^?2Hz) appearing on the radio path with a non-coherent response near the subsolar region. It has been ascertained that the non-coherent response method permits, by using generators with a stability of $\text{ε =}\text{δ}{}_{\mathrm{<mtext>r</mtext>}}\text{f}\text{f}{}_{\mathrm{<mtext>O</mtext>}}\text{= 10}{}^{\mathrm{?11}}\text{?10}{}^{\mathrm{?10}}$, effective investigation (with a good time resolution) of the impulsive interaction of interplanetary shocks with the plasma discontinuities of the bow shockmagnetopause-plasmapause system.     

Determination of magnetic field direction in a comet's head

The Stokes parameters of resonance radiation scattered by a Na atom with the angular momentum F aligned by directed unpolarized radiation in a magnetic field H ～ 10^?5?10^?1 Oe are presented. An influence of the orientation of the magnetic field on these parameters are studied; the intensity ratio $\text{I(D}{}_2\text{)}\text{I(D}{}_1\text{)}$ changes within ±5%, and the polarization degree P(D₂) within ±25%. Measurements of $\text{I(D}{}_2\text{)}\text{I(D}{}_1\text{)}$ and P(D₂), if the geometry of scattering is known, may give information on the direction of the magnetic field in the sodium atmospheres of comets, as well as Io's sodium cloud or man-made cosmic clouds.     

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.     

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$.     

On fourier-analysis of geomagnetic pulsations pc-1, “two-fold” and “three-fold” pulsations pc-1 and fundamental frequencies of the magnetospheric cavity—I

The paper gives the results of detailed studies of the frequency spectra $\text{S}{}_{\mathrm{s}}\text{(?)}$ of the chain of the wave packets F_s(t) of geomagnetic pulsations PC-1 recorded at the Novolazarevskaya station. The bulk of the energy of F_s(t) is concentrated in the vicinity of the central frequencies $\text{?}{}_{\mathrm{s0}}$ of spectra—the carrier frequencies of the signals. The velocity $\text{V}{}_0\text{≌ 6.10}{}^3\text{km s}{}^{\mathrm{?1}}$ of the flux of protons generating these signals correspond to them. The spectra of the signals have oscillations—“satellites” irregularly distributed in frequency. These satellites, as the authors believe, testify to the presence of the individual groups of protons of low concentration whose velocities vary within 10³–10⁴ km s^?1.Their energy is only of the order of 10^?2–10^?3 of the energy of the main proton flux. Clearly pronounced maxima on double and triple frequencies $\text{? = 2?}{}_{\mathrm{s0}}\text{and}\text{3?}{}_{\mathrm{s0}}$ are detected. They show that the generation of pulsations PC-1 is accompanied by the generation on the overtones of wave packets called in this paper “two-fold” and “three-fold” pulsations PC-1. Intensive symmetrical satellites of a modulation character have been discovered on frequencies $\text{?}{}^{\mathrm{\pm }}{}_{\mathrm{sK}}$. Frequency differences $\text{Δ?}{}_{\mathrm{sK}}{}^{\mathrm{\pm }}\text{=}\text{¦?}{}_{\mathrm{s0}}\text{? ?}{}_{\mathrm{sK}}{}^{\mathrm{\pm }}\text{¦}\text{= (0.011,0.022}\text{and}\text{0.035)}\text{Hz}$ correspond to them. The authors believe that the values of Δ?^±_sK are resonance frequencies of the magnetospheric cavity in which geomagnetic pulsations PC-1 are generated. It is established that the values of Δ?^±_sK coincide closely with the carrier frequencies of geomagnetic pulsations PC-3 and PC-4 generated in the magnetosphere. This leads to the conclusion that the resonance oscillations of the magnetospheric cavity are their source. Thus, the generation of geomagnetic pulsations of different types and resonance oscillations in the magnetosphere are integrated into a unified process. The importance of the results obtained and the necessity to check further their trustworthiness and universality, using experimental data gathered in different conditions, is stressed.     

Limb darkening of meteorites and asteroids

This paper presents the results of a laboratory study of the limb darkening, near opposition, of the carbonaceous chondrites Orgueil (C1), Murchison (C2), and Allende (C3), the ordinary chondrite Bruderheim (L6), and a stainless-steel powder. These materials represent possible analogs for the surface materials of C, S, and M asteroids respectively. At low phase angles, the limb-darkening behavior of all materials studied is well represented by Minnaert's law. For carbonaceous chondrites, the Minnaert limb-darkening parameter k is nearly independent of wavelength for wavelengths between 0.4 and 0.9 μm, with a typical value of k = 0.55. The reflectance parameter, B₀, varies from 0.045 to 0.065 over the same range of wavelengths. Both k and B₀ are larger for the stainless-steel powder and the ordinary chondrite, due to the increased importance of multiple scattering in the surface layer. If no limb darkening were present, k would equal $\text{1}\text{2}$ and the geometric albedo (p) of an asteroid would equal the normal reflectance $\text{(r}{}_{\mathrm{<mtext>n</mtext>}}\text{? B}{}_0\text{)}$ of its surface material. For bodies whose surface material is appreciably limb darkened, the geometric albedo measured at the telescope will be lower than the true normal reflectance of surface material; we estimate that for S and M objects r_n ? 1.05 p. In the case of nonspherical asteroids, because the distribution of incidence and emission angles varies as the asteroid rotates, the geometric albedo must change with aspect. If limb darkening is not considered when interpreting asteroid light curves, the values of b/a derived will be too extreme. This effect is probably too small to be observed for C asteroids, because of their intrinsically low reflectances, but could be appreciable for S and M objects.     

Distribution functions for energetic oxygen atoms in the earth's lower atmosphere

Direct photolysis of O₃ and quenching of $\text{O}\text{(}{}^1\text{D}\text{)}$ by N₂ provide abundant sources of fast oxygen atoms for the Earth's lower atmosphere. The concentration of atoms with energy above 0.7 eV may exceed the concentration of $\text{O}\text{(}{}^1\text{D}\text{)}$ for all altitudes below 18 km and these atoms may play an important role in lower atmospheric chemistry. Distribution functions for $\text{O}\text{(}{}^3\text{P}\text{)}$ are given for the energy interval 0.1-1.3 eV, for a range of altitudes from 0 to 62 km.