Anomalous resistivity in the geomagnetic tail region

The beam cyclotron instability and electron acoustic instability, driven by cross-tail current and inhomogeneity in density and magnetic field, are found to be unstable in the earth's magnetic tail region. The anomalous resistivities due to these instabilities are found to be of the order of $\text{(10}{}^{\mathrm{?1}}\text{?10}{}^{\mathrm{?3}}\text{)Ω}{}_{\mathrm{e}}{}^{\mathrm{?1}}$ (Ω_e being the electron gyro frequency). It is also suggested that the non-linear saturation of the beam cyclotron instability may lead to conditions favourable for exciting ion acoustic instability.     

Numerical studies of the transport of solar protons in interplanetary space

Numerical solutions of the Fokker-Planck equation governing the transport of solar protons are obtained using the Crank-Nicholson technique with the diffusion coefficient represented by K_r=K₀r^b where r is radial distance from the Sun and b can take on positive or negative values. As b ranges from +1 to ?3, the time to the observation of peak flux decreases by a factor of 5 for 1 MeV protons when $\text{V}\text{K}{}_0\text{=}\text{3 AU}{}^{\mathrm{b?1}}$ where V is the solar wind speed. The time to peak flux is found to be very insensitive to assumptions concerning the solar and outer scattering boundary conditions and the presence of exponential time decay in the flux does not depend on the existence of an outer boundary. At $\text{V}\text{K}{}_0\text{? 15 AU}{}^{\mathrm{b?1}}$, 1 MeV particles come from the Sun by an almost entirely convective process and suffer large adiabatic deceleration at b?0 but for b=+1, large Fermi acceleration is possible at all reasonable $\text{V}\text{K}{}_0$ values. Implications of this result for the calculation and measurement of particle diffusion coefficients is discussed. At $\text{b?0}$, the pure diffusion approximation to transport overestimates by a factor 2 or more the time to peak flux but as b becomes more negative, the additional effects of convection and energy loss become less important.     

The influence of a flow field on the stability of the force-free magnetic field

In order to analyse the convective instability of the force-free magnetic field, an exact solution of the MHD equation for the magnetic field (1) together with the flow field (2) of constant speed V₀ making an angle θ with the magnetic field, was chosen as the unperturbed state. The stability of the fields between two parallel conducting walls of seperation d was studied by a linear perturbation method, which led to the eigenvalue problem (12), X being given by (13). It was shown by an approximate variational method that instability will set in by the flow field if V₀ is greater than 1/ $\text{3}$ times Alfven velocity V_A. For , the stability of the force-free field (1) is not influenced by the flow field, which may still be significant in other respects. Perturbations transverse to the magnetic field were found to be the most unstable modes.     

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.     

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

Interactions of plasmas with magnetic field boundaries

The magnetopause, the boundary layer, or current sheath, which separates the magnetosphere from the solar wind, is the particular interaction considered in this paper.The collision free electron skin depth, $\text{ξ}{}_{\mathrm{e}}\text{=}\text{c}\text{ω}{}_{\mathrm{pe}}$, where c is the velocity of light and ω_pe, is the plasma frequency, gives a classical measure of the penetration depth of a collisionless plasma by an electromagnetic field. This penetration depth is small compared with the dimensions of the magnetosphere and hence the boundary layer may be conveniently considered in one dimension.In General all one dimensional solutions lie within an order of magnitude of the value of ξ_e, the only exception being the important one, in which the electric field perpendicular to the current sheath plane is not present, either due to a particular trapped particle distribution or due to a short circuiting end effect. For this exception the thickness is increased by the factor ($\text{m}{}_{\mathrm{i}}\text{i/m}{}_{\mathrm{e}}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$.The current sheath solutions discussed are equilibrium solutions but not necessarily stable equilibrium solutions.The extension of the models to three dimensions has a larger effect than might at first be expected. The effect may be intuitively understood as a consequence of flux conservation in the sheath. The one dimensional solutions then correspond to the current sheath profiles at the thinnest point of the three dimensional sheath.     

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.     

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.     

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.     

A note on the asymptotic density-potential relation of a spiral galaxy with a finite thickness

It is shown that the asymptotic σ₁(r) and ψ₁(r) relations can be derived very simply by using the method of double series expansion, where σ₁, ψ₁(r,0) and ψ₁ are the surface density perturbation, the gravitational potential perturbation at the symmetric plane Z=0 and the average potential perturbation respectively. The results are accurate to the order of both m²(kr)^?2 and k〈∣z∣〉, where m is the number of spiral arms, k is the radial wave number, r is the distance from the centre of the galaxy, and 〈∣z∣〉 is the average vertical distance of a star from the Symmetrie plane Z=0. Such an accuracy is usually sufficient for the discussion of spiral modes in a spiral galaxy of small but finite disk thickness. It is pointed out that ψ₁(r,0)～(σ₁(r) relation can be expressed in a unified form for different vertical density profiles if 〈∣z∣〉 is adopted as the thickness scale, and that ψ₁(r,0)～(σ₁(r) can be expressed in a unified form for different vertical density profiles if 〈∣z?z∣〉 the average vertical separation between two stars, is adopted as the thickness scale. Only the value of the ratio $\text{〈|z?z′|〉z}\text{〈|z|〉}$ is a functional of the vertical density profile. However, for the usual physically meaningful profiles, these values are very close to each other: It is $\text{2}$ for the Gaussian profile, $\text{1}\text{Ln}\text{2}\text{= 1.443}$ for the $\text{rmsech}{}^2\text{(}\text{z}\text{z}{}_1\text{(r)}\text{)}$ profile, and 1.5 for the $\text{exp}\text{[}\text{?|z|}\text{z}{}_1\text{(r)}\text{]}$ profile.     

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.     

Analysis of the orbital inclination of Tansei 3 rocket 1977-12B

The orbit of Tansei 3rocket(1977-12B) has been determined at 47 epochs between 1 October 1977 and 19 March 1979 using over 1700 observations and the RAE orbit refinement program PROP6. The rate of change of the inclination was examined to evaluate values of the atmospheric rotation rate, Λ rev day^?1. Analysis yielded the value Λ = 1.1 ± 0.05 at height 315 ± 30 km, average conditions; or alternatively Λ = 1.1 ± 0.1 at height 347 ± 12 km, slight winter bias and Λ = 1.07 ± 0.1 at height 270 ± 18 km, average conditions, supplying further evidence of a decrease in rotation rates from the 1960s to the 1970s.Analysis of the inclination at 15th-order resonance yielded the lumped harmonic values

$\text{10}{}^9\text{C}{}^{\mathrm{0,1}}{}_{15}\text{= 13.4 ± 6.2, 10}{}^9\text{S}{}^{\mathrm{0,1}}{}_{15}\text{= 0.7 ± 13.3}$

for inclination 65.485°.     

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.     

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.     

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.     

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.     

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.     

On the nature of “overshoot” in a collisionless shock

In this paper, by comparing experimental data on bow shock with MHD-relationships on a flat shock discontinuity, allowing for the presence behind the front of turbulent electrostatic oscillations and of an ion beam, an analysis is made of the nature of the “overshoot” of magnetic field (density) behind the front of a collisionless shock wave. It is shown that the large value of plasma compression in the overshoot region ($\text{n}{}_2{}^{\mathrm{ov}}\text{n}{}_1\text{) ～ 6}$, in excess of the maximum allowable value of density jump ($\text{n}{}_2\text{n}{}_1\text{)|}{}_{\mathrm{<mtext>max</mtext>}}\text{= (γ +}\text{1}\text{γ ? 1}\text{)|}{}_{\mathrm{\gamma }}\text{=}{}_{\mathrm{<mtext>5</mtext><mtext>3</mtext>}}\text{= 4}\text{at a Mach number}\text{M → ∞}$, is attributable to the presence in the “overshoot” of a high level of lowerhybrid electrostatic oscillations with an energy density W ? nT.     

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

     

Application of methane band-model parameters to the visible and near-infrared spectrum of Uranus

Laboratory band-model absorption coefficients of CH₄ are used to calculate the Uranus spectrum from 5400 to 10,400 Å. A good fit of both strong and weak bands for the Uranus spectrum over the entire wavelength interval is achieved for the first time. Three different atmospheric models are employed: a reflecting layer model, a homogeneous scattering layer model, and a clear atmosphere sandwiched between two scattering layers. The spectrum for the reflecting layer model exhibits serious discrepancies but shows that large amounts of CH₄ (5–10 km-am) are necessary to reproduce the Uranus spectrum. Both scattering models give reasonably good fits. The homogeneous model requires a particle scattering albedo $\text{(}\text{g}\text{?}\text{w}{}_{\mathrm{<mtext>p</mtext>}}\text{) ? 0.998}$ and an abundance per scattering mean free path $\text{(}\text{a}\text{?}\text{)}\text{of}\text{a}\text{?}\text{1}\text{km-am}$. The parameters derived from the sandwich layer model are: forsb the upper scattering layer a continuum single scattering albedo ($\text{g}\text{?}\text{w}{}_0$) of 0.995 and a scattering optical depth variable with wavelength consistent with Rayleigh scattering; for the clear layer they are a CH₄ abundance (a) of 2.2 km-am and an effective pressure (p) ? 0.1 atm; for the lower cloud deck a Lambert reflectivity (L) of 0.9 resulted. A severe depletion of CH₄ in the upper scattering layer is required. An enrichment of CH₄/H₂ over the solar ratio by a factor of 4–14 in the lower atmosphere is, however, indicated.