Airglow from the inner comas of comets

The intensities of radiation from the inner comas of comets which are composed primarily of water and carbon monoxide have been calculated. Only “airglow” emissions initiated by the absorption of extreme ultraviolet radiation have been considered. The photoionizations of H₂O, CO, CO₂, and N₂ are the most important emission sources, although photoelectron excitation is also considered. Among the emission features for which intensities were calculated are H₂O⁺ ($\text{A}\text{?}{}^2\text{A}{}_{\mathrm{1?}}\text{X}\text{?}{}^2\text{B}{}_1$), CO⁺ (first negative), CO (fourth positive), CO (Cameron), CO₂⁺ ($\text{B}\text{?}{}^2\text{?}{}_{\mathrm{u}}\text{?}\text{X}\text{?}{}^2\text{II}{}_{\mathrm{g}}$), N₂ (Vegard-Kaplan), N₂⁺ (first negative), and OI (1304 Å). In the inner coma (collision region) these airglow mechanisms are shown to be possible competitors with the usually assumed resonance scattering and flourescence excitation mechanisms which are appropriate for the outer coma and tail.     

New high-speed photometry of EH Lib and a binary interpretation of its period change

Six times of maxima of the ultrashort-period cepheid variable EH Librae were measured in 1980 May to June and in 1981 January, with a three-channel photocounting high-speed photoelectric photometer. These, together with all the photoelectric times of maxima over the past 30 years, are used to re-examine the nature of the change of the period. We found that we can fix the times of maxima by the following formula

$\text{T}{}_{\max }\text{= T}{}_0\text{+P}{}_0\text{E+}\text{1}\text{2}\text{βE}{}^2\text{+A}\text{sin}\text{2π}\text{EP}{}_0\text{E}{}_0$

where T₀ = HJD 2433438.6088 and P₀ = 0.0884132445 d represent the initial maximum epoch and the pulsation period, β = ?2.8 × 10^?8/yr; A = 0.0015 d, P₀ = 6251 d = 17.1 yr are the semi-amplitude and the period of the sine curve, and E is the number of periods elapsed since T₀, and (E₀ = 70700).If we interpret this 17.1 year periodicity as a modulation of the phase of maximum by binary motion, then the semi-amplitude of the orbital radial velocity variation is K = 2πasini/E₀ = 0.45 km/s and the mass function is

$\text{f(m)=}\text{m}{}^3{}_2\text{sin}{}^3\text{i}\text{(m}{}_1\text{m}{}_2\text{)}{}^2\text{=}\text{(a}\text{sin}\text{i)}{}^3\text{E}{}^2{}_0\text{=6 x 10}{}^{\mathrm{?5}}\text{M}{}_{\mathrm{\odot }}$

     

Ring current simulation in connection with interplanetary space conditions

A ring current model has been obtained which permits calculations ofD_st variations on the Earth's surface during magnetic storms. The changes in D_st are described by the equation

$\text{d}\text{d}\text{t}\text{D}{}_{\mathrm{sto}}\text{= F(EM)?}\text{D}{}_{\mathrm{sto}}\text{tau;}$

where $\text{D}{}_{\mathrm{sto}}\text{= D}{}_{\mathrm{st}}\text{-bp}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{+~tc}$; p = mnv² is solar wind pressure; F(EM) is the function, controlled by the electromagnetic parameters of interplanetary medium, of injection into ring current ; τ is the constant of ring current decay. $\text{C}\text{= C}\text{u}\text{τ?=18}\text{nT}$, where C is the level of the D_st-variation field measurements; ? is the injection function characterizing the quasisteady-state injection of energy into the ring-current region. The constant Ç is determined from the condition that the change of the ring current energy from magnetic storm commencement to end should equal the difference between the injected and dissipated energy throughout the storm. The values of the factors b and τ were found by the method of minimizing the sum of the quadratic deviations of the calculated D_st from the values observed throughout the storm : b = 0.23 nT/(eV cm^?3)^{$\text{1}\text{2}$}, τ = 8.2 h at D_st? ? 55 nT and τ = 5.8 h at -120 ? D_st ? — 55 nT. The injection function F(EM) is of the form F(EM) = d(E_y? A) at the values of the azimuthal component of the solar wind electric field E_y ? A, and F(EM) =0 at A?E_y.d = ? 1.2 × 10^?3 Ts^?1 (mV/m)^?1 and A = ? 0.9 mV m^?1 . Thus, the injection to ring current is possible at the northward B_z component of the IMF.     

Outlying plasmasphere structure detected by whistlers

Whistlers recorded at Eights ($\text{L ? 4}$) and Byrd ($\text{f ? 7}$), Antarctica have been used to study large-scale structure in equatorial plasma density at geocentric distances $\text{?3–6}\text{R}{}_{\mathrm{E}}$. The observations were made during conditions of magnetic quieting following moderate disturbance. The structures were detected by a “scanning” process involving relative motion, at about one tenth of the Earth's angular velocity or greater, between the observed density features and the observing whistler station or stations. Three case studies are described, from 26 March 1965, 11 May 1965 and 29 August 1966. The cases support satellite results by showing outlying high density regions at $\text{?4–6}\text{R}{}_{\mathrm{E}}$ that are separated from the main plasmasphere by trough-like depressions ranging in width from $\text{?0.2}$ to 1 R_E. The structures evidently endured for periods of 12 hr or more. In the cases of deepest quieting their slow east-west motions with respect to the Earth are probably of dynamo origin. The cases observed during deep quieting (11 May 1965 and 29 August 1966) suggest the approximate rotation with the Earth of structure formed during previous moderate disturbance activity in the dusk sector. The third case, from 26 March 1965, may represent a structure formed near local midnight. The reported structures appear to be closely related to the bulge phenomenon. The present work supports other experimental and theoretical evidence that the dusk sector is one of major importance in the generation of outlying density structure. It is inferred that irregularities of the type reported here regularly develop near 4–5 R_E during moderate substorm activity. This research suggests that at least a major class of the density structures that develop near 4 R_E are tail-like in nature, joined to the main body of the plasmasphere. The apparent disagreement with Chappell's results from OGO 5, which are interpreted as showing regions of “detached” plasma beyond 5 R_E, may be related to the pronounced spatial structure of electric fields observed in high-latitude ionospheric regions that are conjugate to the magnetospheric regions in which the OGO-5 observations were made.     

On the photodissociation of water vapour in the mesosphere

The photodissociation of water vapour in the mesosphere depends on the absorption of solar radiation in the region (175–200 nm) of the O₂ Schumann-Runge band system and also at H-Lyman alpha. The photodissociation products are OH + H, OH + H, O + 2H and H₂ + O at Lyman alpha; the percentages for these four channels are 70, 8, 12 and 10%, respectively, but OH + H is the only channel between 175 and 200 nm. Such proportions lead to a production of H atoms corresponding to practically the total photodissociation of H₂O, while the production of H₂ molecules is only 10% of the H₂O photodissociation by Lyman alpha.The photodissociation frequency (s^?1) at Lyman alpha can be expressed by a simple formula

$\text{J}{}_{\mathrm{<mtext>Ly</mtext><mtext>\alpha </mtext>}}\text{H}{}_2\text{O=4.5 ×10}{}^{\mathrm{?6}}\text{1+0.2}\text{F}{}_{10.7}\text{?65}\text{100}\text{exp}\text{[?4.4 ×10}{}^{\mathrm{?19}}\text{N}{}^{0.917}\text{]}$

where F_{10.7 cm} is the solar radioflux at 10.7 cm and N the total number of O₂ molecules (cm^?2), and when the following conventional value is accepted for the Lyman alpha solar irradiance at the top of the Earth's atmosphere ($\text{Δλ = 3.5}\text{A}\text{?}$) q_Lyα,∞ = 3 × 10¹¹ photons cm^?2 s^1?.The photodissociation frequency for the Schumann-Runge band region is also given for mesospheric conditions by a simple formula

$\text{J}{}_{\mathrm{<mtext>SRB</mtext>}}\text{(}\text{H}{}_2\text{O}\text{) = J}{}_{\mathrm{<mtext>SRB</mtext><mtext>,\infty </mtext>}}\text{(}\text{H}{}_2\text{O}\text{) exp [?10}{}^{\mathrm{?7}}\text{N}{}^{0.35}\text{]}$

where J_SRB,∞(H₂O) = 1.2 × 10^?6 and 1.4 × 10^?6 s^?1 for quiet and active sun conditions, respectively.The precision of both formulae is good, with an uncertainty less than 10%, but their accuracy depends on the accuracy of observational and experimental parameters such as the absolute solar irradiances, the variable transmittance of O₂ and the H₂O effective absorption cross sections. The various uncertainties are discussed. As an example, the absolute values deduced from the above formulae could be decreased by about 25-20% if the possible minimum values of the solar irradiances were used.     

Impact fragmentation experiments of basalts and pyrophyllites

Results of impact fragmentation experiments for basalts and pyrophyllites are reported. Aluminum cylindrical projectiles were impacted on cubic basalt and pyrophyllite targets at velocities of 70 to 990 m/sec. The targets and projectiles were 20 g to 3.3 kg and 2 to 20 g in weight respectively. Weights of the fragments produced by impacts were measured and the size distributions of fragments were examined. Data of the largest fragment mass (m_L) normalized to the original target mass (M_t), m_L/M_t, correlate better with the nondimensional impact stress, P_I, a new scaling parameter introduced by H. Mizutani, Y. Takagi, and S. Kawakami (1984, in preparation) than the conventional projectile's kinetic energy per unit target mass, E/M_t, used in the previous studies. All the m_L/M_t data for basalts obtained in the present study are summarized by m_L/M_t = 2.95 × 10^?2P_I^?1 where P_I = P₀L³/YR³, P₀ = peak shock pressure, L = projectile size, R = target size and Y = material strength of target. For aluminum targets, however, the m_L/M_t is 2.5 orders of magnitude larger than that for brittle targets at impacts with the same P_I. Size distributions of fragments expressed in a log N - log (m/M_t) diagram divided into three regimes bounded by two inflection points. In each regime the curve is expressed by $\text{N (>}\text{m}\text{M}{}_{\mathrm{t}}\text{) = A (}\text{m}\text{M}{}_{\mathrm{t}}\text{)}{}^{\mathrm{?a}}$. The slopes, a, of the log $\text{N -}\text{log}\text{(}\text{m}\text{M}{}_{\mathrm{t}}\text{)}$ curves in the regimes of a large and a medium size range are positively correlated with the nondimensional impact stress, P_I, and expressed as a = C₃ + a₃log P_I. The slopes, a, in the smallest size range are, on the other hand, nearly constant and have values of $\text{0.5}\text{to}\text{0.7 (}\text{1}\text{2}\text{?}\text{2}\text{3}\text{)}$. Present results indicate that the impact fragmentation is scaled well by the new scaling parameter, P_I, of Mizutani, Takagi, and Kawakami and that the present experimental data may shed new light on planetary impact processes.     

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.     

Non-adiabatic processes in a two-fluid plasma and energization of the plasma sheet plasma

The change of energy of a collisionless, two-fluid plasma consists of the adiabatic gain or loss of energy, which is due to the work done by the electromagnetic forces, and of the non-adiabatic change associated with the presence of the “rest” field $\text{E}{}^1\text{=}\text{E}\text{+ (}\text{1}\text{c}\text{)}\text{V}\text{×}\text{B}$. The non-adiabatic gain or loss of energy per unit ti may be expressed by the relation

$\text{Q=}\text{E}{}_{\mathrm{\parallel }}\text{·}\text{i}{}_{\mathrm{\parallel }}\text{+}\text{c}\text{eNB}{}^2\text{B·}\text{f}\text{?}\text{×f}$

where i is the density of conductive current, N the ion number-density, and f (f?) the sum of inertia and pressure divergence of ions (electrons). Symbols of parallelism refer to the direction of B.A special case of non-adiabatic energization of a slowly convecting plasma sheet plasma is discussed in some detail. Regardless of the value of V, the non-adiabatic energization may significantly exceed any conceivable energization associated with the electric field $\text{?(}\text{1}\text{c}\text{)}\text{V}\text{×}\text{B}$.     

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.     

Latitudinal structures of discrete arcs resulting from viscous interaction between sheared plasma flows

Latitudinal structures of discrete arcs are modelled as a consequence of the quasi-steady magnetosphere-ionosphere coupling involving viscous interaction between sunward and anti-sunward plasma flows in the magnetosphere. The quasi-steady state in the magnetosphere and ionosphere coupling is described by the magnetospheric and ionospheric current conservation and the field-aligned currentpotential relation assuming adiabatic electron motion along field lines. The upward and downward fieldaligned currents are assumed to be stably maintained by vorticity-induced space charges in the region of plasma flow reversal, where divergence of the magnetospheric electric field E_⊥ is negative and positive, respectively. By introducing the effective conductance Σ_dc arising from the anomalous viscosity, a specific relation between the dc field-aligned current density J_∥ and the magnetospheric electric field E_⊥ is derived as $\text{J}{}_{\mathrm{\parallel }}\text{=−Σ}{}_{\mathrm{dc}}\text{div}\text{E}{}_{\mathrm{\perp }}$. Sufficiently large potential drops to accelerate auroral electrons are shown to exist along the auroral field lines originating from the flow reversal region with div $\text{E}{}_{\mathrm{\perp }}\text{< 0}$. It is shown that the latitudinal structure of a discrete arc is primarily determined by the magnetospheric potential structure and the characteristic width is on the order of 10 km at the ionospheric altitude.     

Damping of geomagnetic pulsations by the ionosphere

A significant sink of geomagnetic pulsation energy is due to Joule dissipation in the ionosphere. To investigate this we have computed the damping experienced by standing Alfvén waves in a dipole magnetic field. Both the uncoupled poloidal and toroidal modes are considered with Joule dissipation being introduced through a boundary condition which relates the electric and magnetic field strengths at the ionosphere, viz: $\text{4πΣ}{}_{\mathrm{p}}\text{E}\text{c}\text{= b,}\text{where}\text{Σ}{}_{\mathrm{p}}$ is the height integrated Pederson conductivity. The damping rates are strongly dependent on the ionospheric conductivity and we find that typically the normalized damping rate, $\text{γ}\text{ω}\text{,}\text{is}\text{～0.1}$ for nightside values of conductivity and ～0.01 for the dayside. This would account for the observed scale of bandwidths in pulsation signals. Away from regions of extreme damping we find γ ∝ L^?1Σ_p^?1.     

The semi-detached binary system RZ draconis

The 1978 photoeletric observations of the late type close binary RZ Dra were reanalyzed with the Wilson and Devinney approach. Photometric parameters were determined (in Tab. 1.). The system is found to be semi-detached where the less massive component fills its Roche surface, whereas the other component almost does so. The configuration of the system is shown in Fig. 1. The absolute dimensions of the system are found to be M₁ = 0.61M⊙, $\text{M}{}_2\text{= 0.41}\text{M}\text{⊙}\text{R}{}_1\text{= 1.15}\text{R}\text{,}\text{?}\text{and}\text{R}{}_2\text{= 0.96}\text{R}\text{⊙}$. Both components appear to be overluminous and oversized for their masses and spectral types. Its evolutionary stage is also discussed. The variability in the brightness of the primary mlnimum(Fig. 4) indicates mass loss from the vicinity of L₂, which would be mainly responsible for the long-term decrease in its period.     

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.     

Estimates of quasipolar absorption by methane in the atmosphere of Titan

The basis for “quasipolar” absorption (QPA) by CH₄ is the existence of a small electric dipole moment in its ground state. The integrated intensity α_QPA at a temperature of 90K is calculated to be between 4.8 × 10^?5 and 1.9 × 10^?2 cm^?2 atm^?1. With an assumed mean pressure of 0.1 atm and a relative abundance of $\text{[}\text{CH}{}_4\text{]}\text{[}\text{H}{}_2\text{]}\text{= 1}$, it is estimated that the ratio of quasipolar to pressure-induced absorption (PIA) is 0.05 ? α_QPA/α_PIA ? 18 for the spectral range from 0 to 300 cm^?1. This result suggests that quasipolar absorption may contribute to a weak, CH₄-induced greenhouse in the atmosphere of Titan.     

Cross-section for the dissociative photoionization of hydrogen by 584 Å radiation: The formation of protons in the Jovian ionosphere

The cross-section for dissociative photoionization of hydrogen by 584 Å radiation has been measured, yielding a value of 5 × 10^?20 cm². The process can be explained as a transition from the $\text{X}{}^1\text{Σ}{}_{\mathrm{g}}{}^{\mathrm{+}}$ ground state to a continuum level of the $\text{X}{}^2\text{Σ}{}_{\mathrm{g}}{}^{\mathrm{+}}$ ionized state of H₂ The branching ratio for proton (H⁺) vs molecular ion (H₂⁺) production at this energy is 8 × 10^?3. This process is quite likely an important source of protons in the Jovian ionosphere near altitudes where peak ionization rates are found.     

The locating of hydromagnetic whistler sources and determination of their generating proton spectra

Results are given of the calculations of the group delay time propagating τ(ω, φ₀) of hydromagnetic whistlers, using outer ionospheric models closely resembling actual conditions. The τ(ω, φ₀) dependencies were compared with the experimental data of τ_exp(ω, φ₀) obtained from sonagrams. The sonagrams were recorded in the frequency range $\text{? ? (0.5?2.5) Hz}$ at observation points located at geomagnetic latitudes φ₀ = (53?66)° and in the vicinity of the geomagnetic poles. This investigation has led us to new and important conclusions.The wave packets (W.P.) forming hydromagnetic whistlers (H.W.) are mainly generated in the plasma regions at L = 3.5?4.0. This is not consistent with ideas already expressed in the literature that their generation region is $\text{L ? 3?10}$. The overwhelming majority of the τ_exp values differ considerably from the times at which wave packets would, in theory, propagate along the magnetic field lines corresponding to those of the geomagnetic latitudes φ₀ of the observation points. The second important fact is that the W.P. frequency ω is less than $\text{Ω}{}_{\mathrm{H}}$ everywhere along its propagation trajectory, including the apogee of the magnetic force line ($\text{Ω}{}_{\mathrm{H}}$ is the proton gyrofrequency). Proton flux spectra $\text{E ? (30?120)}$ keV, responsible for H.W. generation, were determined. Comparison of the Explorer-45 and OGO-3 measurements published in the literature, with our data, showed that the proton flux density energy responsible for the H.W. excitation $\text{N}{}_{\mathrm{p}}\text{(}\text{MV}{}_6{}^2\text{2}\text{) ? (5 × 10}{}^{\mathrm{?3}}\text{?10}{}^{\mathrm{?1}}\text{)}\text{H}{}_{\mathrm{a}}{}^2\text{8π}$ where H_a is the magnetic field force in the generation region of these W.P. The electron concentration is $\text{N}{}_{\mathrm{a}}\text{? (10}{}^2\text{?10}{}^3\text{) cm}{}^{\mathrm{?3}}$. The values given in the literature are $\text{N}{}_{\mathrm{a}}\text{? (10}{}^{\mathrm{?10}}\text{?10}{}^3\text{) cm}{}^{\mathrm{?3}}$. The e data considered also leads to the conclusion that the generating mechanism of the W.P. studied probably always co-exists with the mechanism of their amplification.     

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.     

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.     

The Rosseland mean opacity of stellar atmospheres

Rosseland mean opacities $\text{a}\text{?}{}_{\mathrm{R}}\text{(T,}\text{lgP}{}_{\mathrm{e}}\text{)}$ [see Table 2] per unit hydrogen particle are calculated in this paper for the atmosphere of the solar type [see Table 1].The sources of the continuous opacities under consideration in this paper are as follows: (1) H^?, HI, H^?₂, H⁺₂ He^?, HeI, HeII, CI, CII, CIII, NI, NII, NIII, OI, OII, NaI, MgI, MgII, AlI, AlII, SiI, SiII, CII, KI, CaII in the form of absorption. (2) HI, HeI, CI, NI, OI, H₂ in the form of Rayleigh scattering. (3) Free electrons in the form of Thomson scattering.     

Thinning of the near-earth (10∼15 RE) plasma sheet preceding the substorm expansion phase

The timing of the plasma-sheet thinning relative to the onset of the expansion phase of substorms is examined by the analysis of the OGO 5 electron (79 ± 23 keV) and proton (100～150 keV) data with the aid of simultaneous magnetic field observations. It is found that the timing of the thinning is significantly dependent on the distance. At $\text{x}{}^2\text{+ y}{}^2\text{? 15 R}{}_{\mathrm{E}}$ the thinning often starts before the onset, while at $\text{x}{}^2\text{+ y}{}^2\text{? 15 R}{}_{\mathrm{E}}$ it tends to occur after the onset, where x and y refer to solar magnetospheric coordinates. The thinning that precedes the expansion-phase onset has been found to reduce the thickness to ～1 R_E, and further thinning may occur in a spatially limited region. Hence it is conceivable that the formation of the neutral line characterizing the substorm expansion phase is the consequence of the thinning of the plasma sheet in the near-Earth region.