Space erosion of meteorites and the secular variation of cosmic rays (over 109 years)

The space erosion of stony meteorites has been determined to be 650μm 10⁶y^?1, while that of iron meteorites has been determined to be 22 μm 10⁶y^?1. The erosion rates are based on flux and size distributions of small particles in the solar system, meteoroid orbitals and the relation, determined by laboratory experiments, between excavated volume due to a collision and the size and velocity of the impacting small particle. Neither multiple collision or space erosion can explain the difference in cosmic ray exposure ages based on ${}^{40}\text{K}$ and those based on ${}^{36}\text{Cl}\text{,}{}^{39}\text{Ar and}{}^{10}\text{Be}$. It is concluded that there is a long term cosmic ray variation.     

The angular momenta of solar system bodies: Implications for asteroid strengths

The angular momentum H is plotted versus mass M for the planets and for all asteroids with known rotation rates and shapes, primarily taken from D. C. McAdoo and J. A. Burns [Icarus18, 285–293 (1973)]. An asteroid's angular momentum is derived from its rotation rate as determined by the period of its lightcurve, its shape as indicated by the lightcurve amplitude, and where possible its size as given by polarimetry or radiometry. The asteroid is assumed to be rotating about its axis of maximum moment of inertia. As previously found by F. F. Fish [Icarus7, 251–256 (1967]) and W. K. Hartmann and S. M. Larson [Icarus7, 257–260 (1967)], H is approximately proportional to $\text{M}{}^{\mathrm{<mtext>5</mtext><mtext>3</mtext>}}$, which shows that the asteroids and most planets spin with nearly the same rate. The very smallest asteroids on the plot deviate from the above reaction, usually containing excess angular momentum. This suggests that collisions have transferred substantial angular momentum to the smallest asteroids, perhaps causing their internal stress states to be substantially modified by centrifugal effects.The forces produced by gravitation are then compared to centrifugal effects for a rotating, triaxial ellipsoid of density 3 g cm^?3. For all asteroids with known properties the gravitational attraction is shown to be larger than the centrifugal acceleration of a particle on the surface: thus the observed asteroid regoliths are gravitationally bound. Poisson's equation for the gravitational potential is investigated and it is shown by mathematical and physical arguments that any arbitrarily shaped ellipsoid with the attractive surface force boundary condition found above will have only attractive internal forces. Thus the internal stress states in asteroids are always compressive so that asteroids could be internally fractured without losing their integrity.     

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.     

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.     

On the exponential decay of western hemisphere solar particle events

Numerical solutions of the Fokker-Planck equation governing the transport of solar protons are obtained in three dimensions (time t, radial distance r, energy T) with the diffusion coefficient represented by κ = κ₀r^bT^a. The October 4, 1968, solar flare particle event is re-examined, and the rise and decay of the proton flux profiles for > 10, ;30 and > 60 MeV particles can be reasonably well reproduced with an instantaneous injection and a distant (10 AU) free escape boundary. The best fit is achieved with a diffusion coefficient $\text{κ = 1.4 × 10}{}^{20}\text{r}{}^{0.5}\text{T}{}^{0.75}\text{cm}{}^2\text{sec}$ where r is in AU and T in MeV.     

An interpretation of Akasofu's substorm parameter

It is proposed that Akasofu's substorm parameter $\text{P = l}{}^2\text{VB}{}^2\text{sin}{}^4\text{(}\text{θ}\text{2}\text{)}$ be interpreted as the electromagnetic power flux which enters the magnetosphere from the solar wind.     

The chemistry of metastable species in the Venusian ionosphere

Reactions of metastable species are important in determining the densities of minor ions in the Venusian ionosphere. Calculations are carried out in which the coupled continuity and momentum equations are solved for twelve ions and four neutral species in the dayside ionosphere, including $\text{O}{}^{\mathrm{+}}\text{(}{}^2\text{D),}\text{O}{}^{\mathrm{+}}\text{(}{}^2\text{P),}\text{N}\text{(}{}^2\text{D),}\text{and N}\text{(}{}^2\text{P)}$. Altitude profiles of these metastable species are presented. Their reactions are shown to be a significant source of several minor ions, especially N₂⁺, CO⁺, and N⁺. The discrepancies which existed between model and measured densities of these ions are resolved.     

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

Models of the formation of the solar nebula

Models of the collapse of a protostellar cloud and the formation of the solar nebula reveal that the size of the nebula produced will be the larger of $\text{R}{}_{\mathrm{<mtext>CF</mtext>}}\text{≡ J}{}^2\text{/k}{}^2\text{GM}{}^3\text{and}\text{R}{}_{\mathrm{<mtext>V</mtext>}}\text{≡ (GMv/2c}{}_{\mathrm{<mtext>c</mtext>}}{}^3\text{)}\text{1}\text{2}$ (where J, M, and c_s are the total angular momentum, total mass, and sound speed of the protosetellar material; G is the gravitational constant; k is a number of order unity; and v is the effective viscosity in the nebula). From this result it can be deduced that low-mass nebulas are produced if P ≡ (R_V/R_CF)² ? 1; “massive” nebulas result if P ? 1. Gravitational instabilities are expected to be important for the evolution of P ? 1 nebulas. The value of J distinguishes most current models of the solar nebula, since P ∝ J^?4. Analytic expressions for the surface density, nebular mass flux, and photospheric temperature distributions during the formation stage are presented for some simple models that illustrate the general properties of growing protostellar disks. It does not yet seem possible to rule out either P ? 1 or P < 1 for the solar nebula, but observed or possible heterogeneities in composition and angular-momentum orientation favor P < 1 models.     

The O(1S) quantum yield from O2+ dissociative recombination

Recent laboratory studies show that the $\text{O}\text{(}{}^1\text{S}\text{)}$ quantum yield, $\text{f}\text{(}{}^1\text{S}\text{)}$, from O₂⁺ dissociative recombination varies considerably with the degree r of vibrational excitation. However, the suggestion that the high values for $\text{f}\text{(}{}^1\text{S}\text{)}$ deduced from airglow and auroral observations can be explained by invoking vibrational excitation, creates a number of problems. Firstly, the rapid vibrational deactivation of O₂⁺ ions by collisions with O atoms will keep r too low to account for the magnitude of $\text{f}\text{(}{}^1\text{S}\text{)}$; secondly, r varies considerably from one atmospheric source to another but its relative values (which should be reliable) do not co-vary with those of $\text{f}\text{(}{}^1\text{S}\text{)}$; thirdly, because r increases markedly above the peak of the $\text{X5577}\text{A}\text{?}$ dissociative recombination layer, the fits which theorists have obtained to the observed volume emission rate profiles would have to be regarded as fortuitious. It is tentatively suggested that $\text{f}\text{(}{}^1\text{S}\text{)}$ is higher in the airglow and aurora than in the laboratory plasma studied by Zipf (1980) because of the electron temperature dependence of the $\text{O}\text{(}{}^1\text{S}\text{)}$ specific recombination coefficient for O₂⁺(v^' ? 3) ions.The repulsive ${}^1\text{Σ}{}_{\mathrm{<mtext>u</mtext>}}\text{[}{}^1\text{D}\text{+}{}^1\text{s}\text{]}$ state of O₂ does not provide a suitable channel for the dissociative recombination. A possible alternative is the bound ${}^3\text{Π}{}_{\mathrm{<mtext>u</mtext>}}\text{[}{}^5\text{S}\text{+}{}^3\text{s}\text{]}$ state with predissociation to the repulsive ${}^3\text{Π}{}_{\mathrm{<mtext>u</mtext>}}\text{[}{}^3\text{P}\text{+}{}^1\text{s}\text{]}$ state.     

Geopotential coefficients of order 29, from analysis of the orbit of satellite 1967-104B

The orbit of the satellite 1967-104B has been analysed as it passed through 29:2 resonance with the Earth's gravitational field between January 1977 and September 1978. From the changes in inclination and eccentricity the following lumped 29th-order geopotential harmonic coefficients were obtained: $\text{10}{}^9\text{C}\text{?}{}_{29}{}^{0.2}\text{= 4.1 ± 0.8}$, $\text{10}{}^9\text{S}\text{?}{}_{29}{}^{0.2}\text{= 10.3 ± 2.4}$, $\text{10}{}^9\text{C}\text{?}{}_{29}{}^{1.1}\text{= ? 160 ± 19}$, $\text{10}{}^9\text{S}\text{?}{}_{29}{}^{1.1}\text{= 79 ± 10}$, $\text{10}{}^9\text{C}\text{?}{}_{29}{}^{\mathrm{?1.3}}\text{= 38 ± 14}$, $\text{10}{}^9\text{S}\text{?}{}_{29}{}^{\mathrm{?1.3}}\text{= 19 ± 5}$. These values have been compared with existing comprehensive geopotential models: the best agreement is with the model of Rapp (1981).     

Theory for the critical ionization velocity phenomenon

A theoretical consideration is given to the critical ionization velocity. A macroscopic theory is constructed, assuming that (1) the newly born ions have an unstable velocity distribution and (2) the resultant fluctuating electric field heats the electrons. For the critical ionization velocity to operate in the case of $\text{m}{}_{\mathrm{a}}\text{V}{}^2\text{2}\text{> eφ}{}_{\mathrm{<mtext>ion</mtext>}}$, the initial values of high speed ion (or newly ionized ion) density and suprathermal electron energy density should be above a certain threshold. This threshold depends upon the velocity diffusion co-efficient and the collisional loss rate of the newly ionized ion.     

Branching ratios for the photoionization cross section of atomic oxygen

Branching ratios $\text{σ(}\text{O}{}^0{}^3\text{P}\text{→}\text{O}{}^{\mathrm{+}}{}^2\text{D}{}^0\text{)}\text{σ(}\text{O}{}^0{}^3\text{P}\text{→}\text{O}{}^{\mathrm{+}}{}^4\text{S}{}^0\text{)}$ and $\text{σ(}\text{O}{}^0{}^3\text{P}\text{→}\text{O}{}^{\mathrm{+}}{}^2\text{P}{}^0\text{)}\text{σ (}\text{O}{}^0{}^3\text{P}\text{→}{}^4\text{S}{}^0\text{)}$ are calculated at 584 Å and 304 A employing the close-coupling approximation to compute the photoionization cross section values. The coupled channels include the states dominated by the ground configuration 1s²2s²p³ of O⁺and the next excited configuration ls²2s2p⁴. It is found that the partial c section $\text{σ(}{}^2\text{D}{}^0\text{)}$ decreases more rapidly than $\text{σ(}{}^2\text{P}{}^0\text{)}$, and at the lower wavelength 304 Å, the ratio $\text{σ(}{}^2\text{D}{}^0\text{)}\text{σ(}{}^4\text{S}{}^0\text{)}\text{<}\text{σ(}{}^2\text{P}{}^0\text{)}\text{σ(}{}^4\text{S}{}^0\text{)}$. Present results at 304 Å differ considerably from previous work.     

Photoionization of barium clouds via the 3D metastable levels

The photoionization of optically thin barium clouds is analyzed and shown to occur primarily by a two-step process involving the ${}^3\text{D}$ metastable term as the ntermediate state. The equilibrium populations of the ${}^1\text{D}$ and ${}^3\text{D}$ metastable levels are calculated and found to differ significantly from the values now in the literature. These populations are combined with our newly available photoionization data to calculate the resulting photoionization rate in the upper atmosphere.     

A possible mechanism of formation of faculae

We propose a new heating mechanism of faculae. We think that the formation of faculae is a result of the Joule dissipation of the Hall current generated by the interaction of the convection field of granules in an active region and the inter-granular magnetic field. For a region to generate effectively Hall current, its characteristic length must be such that the magnetic Reynolds number is less than 1. The equation of energy balance in the facula region is

$\text{16σT}{}^3{}_{\mathrm{p}}\text{(T}{}_{\mathrm{l}}\text{? T}{}_{\mathrm{p}}\text{)}{}^{\mathrm{n}}\text{HP}{}_{\mathrm{s}}{}^{\mathrm{a}}\text{H}{}^{\mathrm{?}}\text{=}\text{Qn}{}_{\mathrm{s}}\text{m}{}_{\mathrm{i}}\text{u}{}_{\mathrm{x}}{}^2\text{(τ}{}^2{}_{\mathrm{in}}\text{ω}{}_{\mathrm{i}}\text{)}$

.For five observational models of faculae, we calculated the corresponding velocity fields, and the results are in basic agreement with the observed fields. The present mechanism explains the dependence of the facula brightness on the magnetic and velocity fields, the apparent distribution of the faculae on the solar disk and suggest a possible interpretation of the five structures of faculae.     

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.     

Direct measurement of conjugate photoelectrons and predawn 630 nm airglow

A sounding rocket was flown during the predawn on 17 January, 1976 from Uchinoura, Japan, to measure directly the behaviour of the conjugate photoelectrons at magnetically low latitudes. On board the rocket were an electron energy analyzer, 630 nm airglow photometer, and plasma probes to measure electron density and temperature. The incoming flux of the photoelectrons was measured in the altitude range between 210 and 340 km. The differential flux at the top of the atmosphere was determined to be $\text{F = (1.3 ± 0.4) × 10}{}^{11}\text{exp}\text{[?}\text{E(}\text{eV}\text{)}\text{12}\text{]}$ electron · m^?2 · sr^?1 · s^?1 in the energy range 10 ? E ? 50 eV. The emission rate of the 630 nm airglow was observed in the altitude range between 90 and 360 km. The apparent emission rate observed at 80 km was 32 ± 5 R. From a theoretical calculation of the optical excitation rate using the observed electron flux data along with a model distribution of atomic oxygen, it was estimated that more than 65% of the emission could be produced by direct impact of the photoelectrons with atomic oxygen in the thermosphere between 200 and 360 km. Using the observed electron density and the model distribution of oxygen molecules the residual of the emission was ascribed to the excitation of O(¹D) through dissociative recombination, $\text{O}{}_2{}^{\mathrm{+}}\text{+}\text{e}\text{→}\text{O}{}^1\text{+}\text{O}{}^7$. The direct collisional excitation by ambient electrons is estimated to be negligibly small at the level of observed electron temperature.     

Airborne spectroscopy and spacecraft radiometry of Venus in the far infrared

In March 1979, the spectrum of Venus was recorded in the far infrared from the G.P. Kuiper Airborne Observatory when the planet subtended a phase angle of 62°. The brightness temperature was observed to be 275°K near 110 cm^?1, dropping to 230°K near 270 cm^?1. Radiance calculations, using temperature and cloud structure formation from the Pioneer Venus mission and including gaseous absorption by the collision-induced dipole of CO₂, yield results consistently brighter than the observations. Supplementing the spectral data, Pioneer Venus OIR data at similar phase angles provide the constraint that any additional infrared opacity must be contained in the upper cloud, H₂SO₄ to the Pioneer-measured upper cloud structure serves to reconcile the model spectrum and the observations, but cloud microphysics strongly indicates that such a high particle density haze $\text{(N ? 1.6 × 10}{}^7\text{cm}{}^{\mathrm{?3}}\text{)}$ is implausible. The atmospheric environment is reviewed with regard to the far infrared opacity and possible particle distribution modifications are discussed. We conclude that the most likely possibility for supplementing the far-infrared opacity is a population of large particles $\text{(}\text{r}\text{? 1 μm)}$ in the upper cloud with number densities less than 1 particle cm^?3 which has remained undetected by in situ measurements.