Dynamics of magnetospheric ion Composition as observed by the GEOS mass spectrometer

After one year of operation the GEOS-1 Ion Composition Experiment has surveyed plasma composition at all local times in the L range 3 8 and the energy per charge range from thermal to 16 keV/e. From measurements made in the keV range during eleven magnetic storms we find that the percentage of heavy (M/Q > 1) ions present in the outer magnetosphere increases by a factor of 3 to 10 during disturbances. We conclude that two independent sources (solar wind, characterized by ⁴He²⁺, and ionosphere, characterized by O⁺) give on the average comparable contributions to injected populations, although in a single event one or the other source may dominate. However, in magnetically quiet periods protons are the dominant species with a few percent of heavy ions. With the help of special satellite manoeuvres magnetic field aligned fluxes of 0.05-3 keV/e H⁺, He⁺, O⁺ with traces of O²⁺ have been observed which may be related to ion beams found previously at lower altitudes in the auroral zone. At still lower energies ( 1 eV/e) the thermal plasma population is found to be made up of six ion species, three of which, D⁺, He²⁺ and O²⁺, were unknown in the magnetosphere prior to the GEOS-1 measurements. We present here a study of the evolution of doubly charged ions and their parent populations over four consecutive days. Various production mechanisms for doubly charged ions are discussed. We argue that ionization of singly charged ions by UV and energetic electrons and protons is the dominant process for plasmasphere production. Furthermore, the observed high concentrations of O²⁺ at high altitudes are a result of production in the upper ionosphere and plasmasphere combined with upward transport by thermal diffusion. Throughout the 1 year lifetime of GEOS-1 the ICE functioned perfectly and, because of its novel design, a short review of technical performance is included here.     

Pulsars as gamma ray sources: Nebular shocks and magnetospheric gaps

I summarize the results of recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is primarily electrons and positrons with an admixture of heavy ions. Shocks which contain heavy ions that are a minority constituent by number but which carry most of the energy density in the upstream medium put 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) E ^-2, where N(E)dE is the number of particles with energy between E and E+dE. Synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front, provides the mechanism of thermalization and non-thermal particle acceleration. The maximum energy achievable by the pairs is _± m _± c ² = m _i c ² ₁/Z _i, where ₁ is the Lorentz factor of the upstream flow and Z _i is the atomic number of the ions. The shock's spatial structure contains a series of overshoots in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. These overshoots provide a new interpretation of the structure of the inner regions of the Crab Nebula, in particular of the wisps, surface brightness enhancements near the pulsar. The wisps appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure ₁ 4 × 10⁶, the upstream magnetic field B ₁ to be 3 × 10^-5 Gauss, as well as to show that the total ion flow is 3 × 10³⁴ elementary charges/sec, in good agreement with the total current flow predicted by the early Goldreich and Julian (1969) model. The total pair outflow is shown to be about 5 × 10³⁷ pairs per second, in good agreement with the particle flux required to explain the nebular X—ray source.The energetics of particle acceleration within the magnetospheres of rotation powered pulsars and the consequences for pulsed gamma ray emission are also briefly discussed. The gamma ray luminosity above 100 MeV is shown to scale in proportion to _R ^1/2 , as is in accord with some of the simplest ideas about polar cap models. Models based on acceleration in the outer magnetosphere are also briefly discussed.     

Design by extrapolation: an evaluation of fault tolerant avionics

Over the past 30 years, safety-critical avionics systems such as Fly-By-Wire (FBW) flight controls, full-authority digital engine controls, and other systems have been introduced on many commercial and military airplanes and spacecraft. Early FBW systems, such as on the F-16 and Airbus A320, were considered revolutionary and introduced with extreme caution. These early systems and their successors all make use of redundant and fault-tolerant avionics to provide the required dependability and safety, but have used significantly different architectures. This paper examines the different levels of criticality and fault tolerance required by different types of avionics systems, establishes architectural categories of fault-tolerant architectures, and identifies the discriminating features of the varied approaches. Examples of discriminators include the level of redundancy, methods of engaging backup systems, protection from software errors, and the use of dissimilar hardware and software. The strengths and weaknesses of the approaches will be identified. The paper concludes with some speculation on trends for future systems based on this evaluation of previous systems     

Laboratory experiments of magnetospheric interest

Space-related laboratory experiments can play an important role as a complement to observations and active experiments in the magnetosphere. Excluding laboratory experiments for mere developing or testing of techniques for space experiments, we may distinguish between two major types: (1) partial scale model experiments and (2) experiments for clarifying basic plasma physical processes known or expected to be important in the magnetosphere (but without the ambition to simulate actual space configurations). The limitations and potentialities of both types are discussed and examples of experiments are given. It is concluded that there should be an increasing need for the experiments of the second type. In particular, they are needed for the clarification of the response of a thin plasma to electric fields and its ability to carry electric currents. This encompasses such key questions as the nature and role of anomalous resistivity (and electron runaway in its presence), the possible formation of double layers (and the acceleration processes associated with them) and rapid dissipation of magnetic-field energy.     

Modelling of the magnetic field of magnetospheric ring current as a function of interplanetary medium parameters

The models are examined which are proposed elsewhere for describing the magnetic field dynamics in ring-currentDR during magnetic storms on the basis of the magnetospheric energy balance equation. The equation parameters, the functions of injectionF and decay , are assumed to depend on interplanetary medium parameters (F and during the storm main phase) and on ring-current intensity ( during the recovery phase). The present-day models are shown to be able of describing theDR variations to within a good accuracy (the r.m.s. deviation 5 < < 15 nT, the correlation coefficient 0.85 <r < 1). The models describe a fraction of the geomagnetic field variation during a magnetic storm controlled by the geoeffective characteristic of interplanetary medium and, therefore responds directly to the variation of the latter. The fraction forms the basis of the geomagnetic field variations in low and middle latitudes. The shorter-term variations ofDR are affected by the injections into the inner magnetosphere during substorm intervals.During magnetic storms, the auroral electrojets shift to subauroral latitudes. When determining theAE indices, the data from the auroral-zone stations must be supplemented with the data from subauroral observatories. Otherwise, erratic conclusions may be obtained concerning the character of the relationships ofDR toAE or ofAE to interplanetary medium parameters. Considering this circumstance, the auroral electrojet intensity during the main phase is closely related to the energy flux supplied to the ring current. It is this fact that gives rise simultaneously to the intensification of auroral electrojets and to the large-scale decrease of magnetic field in low latitudes.The longitudinal asymmetry of magnetic field on the Earth's surface is closely associated with the geoeffective parameters of interplanetary medium, thereby making it possible to model-estimate the magnetic field variations during magnetic storms at given observatories. The inclusion of the field asymmetry due to the system of large-scale currents improves significantly the agreement between the predicted and model field variations at subauroral and midlatitude observatories. The first harmonic amplitude of field variation increases with decreasing latitude. This means that the long-period component of theD _st -variation asymmetry is due rather to the ring-current asymmetry, while the shorter-term fluctuations are produced by electrojets. The asymmetry correlates better with theAL indices (westward electrojet) than with theAU indices (eastward electrojet).The total ion energy in the inner magnetosphere during the storm main phase is sufficient for the magnetic field observed on the Earth's surface to be generated. The energy flux to the ring current is 15% of the -energy flux into the magnetosphere.     

Consequences of hydromagnetic waves on magnetospheric particle dynamics

Since long ULF waves are known to play a dominant role in the dynamics of energetic protons. Recent observations have shown that they also contribute to the heating and/or acceleration of both heavy ions and electrons. This review will deal with all the aspects of wave particle interactions that involve electromagnetic ULF waves in the equatorial region of the magnetosphere. We will consider successively Pc-3-4-5 pulsations, magnetosonic waves and ion cyclotron waves. In the latter case, for which both the experimental data and the theoretical interpretation are more advanced, we will discuss the linear, quasi-linear and nonlinear aspects of the interaction. Results of recent numerical simulations of this type of interaction will also be reported. A final section will be devoted to ULF waves observed in the vicinity of the Io torus.     

Study of magnetospheric dynamics using energetic solar particles

Information can be obtained from energetic particle measurements through the chemical composition, energy spectrum, directional anisotropy, temporal and spatial intensity variations. This is equivalent to saying that there is a distribution functionf _k(p,r,t) wherek corresponds to thekth particle species of momentump at positionr and timet.Particle transport is described by the Boltzmann equation, and because the densities are generally low in the case of cosmic rays or energetic solar flare particles, collective transport effects can be neglected. In the absence of magnetospheric motion it is relatively easy to treat the problems of particle transport as simple propagation of charged particles in a stationary magnetic field configuration using, for instance, trajectory calculations in model fields. The method here is to use correlated measurements of the particle distribution at two points along a dynamic trajectory, and in this way to learn something about the geomagnetic field. This approach provides a good basis from which to study magnetospheric dynamics. If the magnetosphere moves, large scale electric fields, turbulent electromagnetic fields and sources and sinks affect the propagation of energetic particles considerably. These effects change the distribution functionf _k(p,r,t) and can thus be detected.In this paper, we shall show the importance of the single particle approximation (trajectories in a reference field) in forming the basis of our understanding of the quiet-time penetration of cosmic rays into the magnetosphere, we shall consider the steady dynamics such as wave-particle inter-action and field line reconnection, which is believed to exist nearly all the time, and finally we shall review the work which has been done in the much more complex and less well-understood field of impulsive dynamics such as geomagnetic storms and substorms. This last topic is only just beginning to be investigated in detail, and it is hoped that the study of impulsive dynamics, using energetic particles, may be as successful as the study of the quiet magnetosphere and the steady dynamics.     

Laboratory experiments directed toward the investigation of magnetospheric phenomena

The penetration of fast electrons ( 5 keV) into an artificial magnetosphere and their precipitation on the terrella surface is investigated. These fast electrons act as radioactive tracers allowing the experimental determination of the global picture of plasma flow around the magnetosphere and its intrusion into the latter. Two different zones of precipitation are observed, distinctly separated on the day-side and merging into each other on the night-side. The high latitude penetration region on the day-side is not localized around the neutral points, but is stretched in longitude forming polar cusps toward dusk and dawn. The lower latitude precipitation zone, embracing the whole terrella is due to the particle precipitation from a radiation belt formed in this experiment. The source of these belt particles seems to be located in the plasmasheet on the night side. Besides the polar cusps, a plasma intrusion from the sides of the magnetosphere in the equatorial region is observed. This equatorial gap, originating on the day-side, is gradually transformed into the plasmasheet in the magnetospheric tail. On the basis of these experimental data a model of the magnetosphere is discussed.     

Observations of low energy magnetospheric plasma outside the plasmasphere

After some introductory discussions about morphological concepts and limitations of various measurement techniques, existing low energy plasma data, orginating primarily from the GEOS, Dynamics Explorer, and Prognoz spacecraft, is described and discussed. The plasmasphere measurements are not included (but for some observations of plasmasphere refilling). It is finally concluded that we are very far from a complete picture of the low-energy plasma component in the magnetosphere and that this problem has to be given high priority in planning payloads of future space plasma physics missions.     

Steady magnetospheric convection: A review of recent results

Theoretical pressure balance arguments have implied that steady convection is hardly possible in the terrestrial magnetotail and that steady energy input necessarily generates a cyclic loading-unloading sequence, i.e., repetitive substorms. However, observations have revealed that enhanced solar wind energy input to the magnetospheric system may either lead to substorm activity or enhanced but steady convection. This topic is reviewed with emphasis on several recent case studies of the Steady Magnetospheric Convection (SMC) events. In these cases extensive data sets from both satellite and ground-based instruments from various magnetospheric and ionospheric regions were available.Accurate distinction of the spatial and temporal scales of the magnetospheric processes is vital for correct interpretation of the observations during SMC periods. We show that on the large scale, the magnetospheric configuration and plasma convection are stable during SMC events, but that both reveal considerable differences from their quiet-time assemblies. On a shorter time scale, there are numerous transient activations which are similar to those found during substorms, but which presumably originate from a more distant tail reconnection process, and map to the poleward boundary of the auroral oval. The available observations and the unresolved questions are summarized here.The tail magnetic field during SMC events resembles both substorm growth and recovery phases in the neartail and midtail, respectively, but this configuration may remain stable for up to ten hours. Based on observations and model results we discuss how the magnetospheric system avoids pressure balance problems when the plasma convects earthward.Finally, the importance of further coordinated studies of SMC events is emphasized. Such studies may shed more light on the substorm dynamics and help to verify quantitatively the theoretical models of the convecting magnetosphere.     

Transport of energetic solar particles on closed magnetospheric field lines

Whereas the entry mechanism of energetic solar particles into the open field line region of the magnetosphere is now a rather well understood process, transport processes of solar particles in the closed field line region are still unclear and under dispute. The main difficulty lies not only in the fact that different field models predict different behavior of the particles in the quasi-trapping region (e.g. cut-off latitude), but that dynamic changes of the magnetosphere as geomagnetic storms and substorms greatly influence the particle distribution. The present review tries to summarize the status of knowledge regarding solar proton behavior on closed magnetospheric field lines. Together with a presentation of recent measurements in the closed field line region relevant theoretical problems are discussed. They fall either under the study of single particle motion in different static magnetospheric configurations (due to different field models or due to real, e.g. ring current induced changes), or under the study of resonant interaction processes as pitch angle scattering and radial diffusion.Invited Lecture, Second Meeting of the European Geophysical Society, September 1974, Trieste, Italy.     

Conjugate features of magnetospheric electron dynamics observed at balloon altitudes

Analysis of recent observations (from balloons, spacecraft, and surface observatories) demonstrate regional, shell, and nearpoint conjugacy at L ～ 7 during precipitative events which were characterized by local acceleration as well as release of gradient-drifted electrons injected during substorms. A number of new features of magnetospheric dynamics relating to substorm development and sudden-commencement effects, have been brought to light which, though poorly understood at present, may prove of considerable importance and are worthy of further investigation.

1. During the initial period of instability in substorm evolution, preceding the slower magnetotail convective injection, precipitation of waves of electrons in rapid polewards motion exhibit L-shell conjugacy near midnight.
2. Transient, large scale expansions of the magnetospheric electron population accompanied by temporally imbedded substorms display large scale regional conjugacy and are simultaneously observed as similarly transient intensity dropouts at balloon altitudes.
3. Precipitation from gradient-drifting electrons in the dayside magnetosphere exhibits near point-conjugacy, at least down to the order of 50 km and quite probably less.

Similarly tight conjugacy applies to the release of electrons showing a specific local response to sudden commencements.

1. Analysis of the approach to and attainment of spectral equilibrium in the precipitation observed from drifting electrons may provide information about either, or both, the source spectrum at injection and the process of local release.
2. The specific precipitation effect sometimes observed at the time of an SC remains a rather puzzling feature, although it seems clear now that the acceleration and/or release process responsible is of a highly local nature and works selectively at small pitch angles well within the magnetospheric boundary. Coupling of the interplanetary shock with the magnetosphere must be an important aspect, but the details are not clear as yet.
3. On at least one occasion, a large part (perhaps all) of the magnetospheric electron population varied in a nearly synchronous manner in response to solar wind induced distortions during the variable compressive phase of a sudden commencement geomagnetic storm.

In the ongoing effort to identify and understand acceleration and release mechanisms involved in magnetospheric dynamics, balloon-borne experiments will continue to be useful, providing essential information presently unattainable by other means.     

Extraction of ion distributions from magnetospheric ENA and EUV images

Energetic neutral atom (ENA) and extreme ultra-violet photon (EUV) imagers will soon be probing magnetospheric ion distributions from the NASA space missions IMAGE and TWINS. Although ENA and EUV images will differ greatly, the same basic mathematical approach can be applied to deducing the ion distributions: extracting the parameters of a model ion distribution in a model magnetic field (and, in the case of ENA, interacting with a model cold neutral population). The model ion distribution is highly non-linear in its many parameters (as many as 38 have been used) in order to describe the strong spatial gradients of ion intensities in the magnetosphere. We have developed several new computer algorithms to accomplish the extraction by minimizing the differences between a simulated (instrument-specific) image and an observed image (or set of images). Towards the goal of a truly automated `hands-off extraction algorithm, we have combined three algorithms into a Hierarchical Simplex Algorithm. At each step of the minimization, it first tries a sophisticated and efficient Adaptive Conjugate Gradient algorithm. Then, if the error function is not reduced, it defers to an intermediate Analytic Simplex algorithm, and (if this step also fails) it finally defaults to the robust but inefficient Downhill Simplex algorithm. Whenever a step is successful, the algorithm begins the next step at the top of the hierarchy. We also offer a completely different approach (without minimization) for the interpretation of sharp `edges in the images (e.g., the plasmapause in He ⁺ 30.4 nm EUV images of the plasmasphere). We demonstrate mathematically that the equatorial shape of the plasmapause can be constructed directly from the image using a simple graphical algorithm.     

Quantitative models of the magnetospheric magnetic field: Methods and results

The paper reviews various approaches to the problem of evaluation and numerical representation of the magnetic field distributions produced within the magnetosphere by the main electric current systems including internal Earth's sources, the magnetopause surface current, the tail plasma sheet, the large-scale systems of Birkeland current, the currents due to radiation belt particles, and the partial ring current circuit. Some basic physical principles as well as mathematical background for development of magnetospheric magnetic field models are discussed.A special emphasis is placed on empirical modeling based on datasets created from large bodies of spacecraft measurements. A review of model results on the average magnetospheric configurations and their dependence on the geomagnetic disturbance level and the state of interplanetary medium is given. Possibilities and perspectives for elaborating the instantaneous models capable of evaluating a current distribution of magnetic field and force line configuration based on a synoptic monitoring the intensity of the main magnetospheric electric current systems are also discussed. Some areas of practical use of magnetospheric models are reviewed in short. Magnetospheric plasma and energetic particle measurements are considered in the context of their use as an independent tool for testing and correcting the magnetic field models.     

The radial diffusion of trapped particles induced by fluctuating magnetospheric fields

Small fluctuations in magnetospheric electric and magnetic fields lead to random changes in the radial positions of trapped particles. The characteristics of this radial diffusion are described theoretically in terms of the statistical properties of the field fluctuations, in particular the power spectra of the various spatial components. A large body of trapped particle data demonstrates that diffusion with the predicted properties actually takes place. These data include the average radial and energy distributions and the time variations in particle fluxes. This radial diffusion has a major influence on the structure of the radiation belts; and since the net flow of particles is inward at most positions of the magnetosphere, the process acts as a strong source of trapped particles. Further experiments are needed to establish the importance of this mechanism relative to magnetic storm effects and to collective instabilities.     

The magnetogram inversion technique: Applications to the problem of magnetospheric substorms

The magnetogram inversion technique (MIT) is based upon recordings of geomagnetic variations at the worldwide network of ground-based magnetometers. MIT ensures a calculation of a global spatial distribution of the electric field, currents and Joule heating in the ionosphere. Variant MIT-2 provides, additionally, continuous monitoring of the following parameters: Poynting vector flux from the solar wind into the magnetosphere (); power, both dissipated and accumulated in the magnetosphere; magnetic flux in the open tail; and the magnetotail length (l _T) (distance between the dayside and nightside neutral points in the Dungey model). Using MIT-2 and data of direct measurements in the solar wind, an analysis is made of a number of substorms, and a new scenario of substorms is suggested. The scenario includes the convection model, the model with a neutral line and the model of magnetosphere-ionosphere coupling (outside the current sheet), i.e., the three known models. A brief review is given of these and some other substorms models. A new element in the scenario is the strong positive feedback in the primary generator circuit, which ensures growth of the ratio = / _Aby an order of magnitude or more during the substorms. Here _Ais the Pointing vector flux in the Akasofu-Perrault approximation, i.e., without the feedback taken into account. The growth of during the substorm is caused only by the feedback effect. It is assumed that the feedback arises due to an elongation of the magnetotail, i.e., a growth of l _Tby a factor of (23) during the substorm.In the active phase of substorm, a part (the first active phase) has been identified, where the principal role in the energetics is played by the feedback mechanism and the external energy source (although the internal source plus reconnection inside the plasma sheet make a marked contribution). In the second active phase (expansion) the external generator (solar wind) is switched off, and the main role is now played by the internal energy source (the tail magnetic field and ionospheric wind energy).Models of DP-2 DP-1 transitions are also considered, as well as the magnetospheric substorm-solar flare analogy.     

Small scale structure of magnetospheric electron density through on-line tracking of plasma resonances

The plasma resonance phenomena observed at f _pe, nf _ce, and f _qn by the GEOS-1 S-301 relaxation sounder are identified through a pattern recognition software process implemented in a mini-computer which receives on-line the compressed data. First, this processing system distributes in real time f _pe and f _ce measurements to the ground media. Second, it drives and controls automatically the S-301 on-board experiment by sending appropriate telecommands: the tracking of resonances is performed by shortening the frequency sweeps to a narrow range centered on the resonance location. Examples of such tracking sequences are presented, exhibiting sampling rates of the electron density measurements from once every 22 s (slowest rate) to once every 86 ms (highest rate available). The results give evidence of the existence of very small scale structures in the magnetospheric density, having characteristic sizes of the order of a few 10² m or/and a few 10^-1 s. The relative amplitude of these density fluctuations is typically 1%. Because of satellite spinning, fixed frequency sounding sequences allow to measure in a few seconds the directivity features of the plasma resonance signals. Examples of directional patterns in the plane perpendicular to the geomagnetic field are presented: the electrostatic nature of the waves received at f _pe, nf _ce, and f _qe being consistent with these patterns, the corresponding k vector orientations become available. The Bernstein modes properties are used to interpret the cf _ce and f _qe results.     

基于广义Richardson外插方法的阻力预测精度分析

利用计算流体力学(CFD)方法预测阻力是飞行器气动设计中的关键环节。采用广义Richardson外插方法分别对数值方法预测二维简单构型的压差阻力、摩擦阻力和三维复杂构型的总阻力的精度进行了分析。在NACA 0012翼型无黏绕流和平板湍流边界层两个算例中验证了NSAWET程序和广义Richardson外插方法,分别得到数值算法预测压差阻力和摩擦阻力能达到的名义精度。进而模拟三维通用研究模型(CRM)翼身组合体绕流,得到的阻力名义精确值在DPW 5的统计误差带范围之内;综合DPW 5的计算结果来看,不同CFD解算器的结果之间存在一定差别,阻力预测精度总体上不符合二阶。可见,标准Richardson方法采用的二阶精度假设难以普遍适用,有必要采用广义Richardson外插方法得到名义精度。针对不合理的名义精度,采用Roache建议的方法加以限制。广义Richardson外插方法有助于提高误差分析的合理性,可以进一步降低网格对阻力预测的影响。     

波长外推法在多波长温度测量中的应用

波长外推法作为多波长温度测量法一种数据处理方法,可以充分利用多个波长下测量得到的数据,经过最小二乘线性拟合得到温度测量结果,具有稳定性好,精确度高等优点。文章详细介绍了波长外推法的基本原理和使用方法,并搭建平台进行试验验证,结果表明这种方法对灰体以及在某一波段内近似灰体的目标能够得到很好的结果。