Nongyrotropic particle distributions in space plasmas

In nonstationary, strong inhomogeneous or open plasmas particle orbits are rather complicated. If the nonstationary time scale is smaller than the gyration period, if the inhomogeneity scale is smaller than the gyration radius, i.e. at magnetic plasma boundaries, or if the plasma has sources and sinks in phase space, then nongyrotropic distribution functions occur. The stability of such plasma configurations is studied in the framework of linear dispersion theory. In an open plasma nongyrotropy drives unstable waves parallel and perpendicular to the background magnetic field, whereas in the gyrotropic limit the plasma is stable. In nonstationary plasmas nongyrotropy drives perpendicular unstable waves only. Temporal modulation couples a seed mode with its side lobes and thus it renders unstable wave growth more difficult. As an example of an inhomogeneous plasma a magnetic halfspace is discussed. In a layer with thickness of the thermal proton gyroradius a nongyrotropic distribution is formed which may excite unstable parallel and perpendicular propagating waves.     

On the dispersion of two coexisting nongyrotropic ion species

Space observations in the solar wind and simulations of high Mach number bow-shocks have detected particle populations with two coexisting nongyrotropic ion species. We investigate the influence of these two sources of free energy on the stability of parallel (with respect to the ambient magnetic field) and perpendicular propagation. For parallel modes, we derive their dispersion equation in a magnetoplasma with protons and alpha particles that may exhibit stationary nongyrotropy (SNG) and discuss the characteristics of its solutions. Kinetic simulations study the behaviour of perpendicular electrostatic (Bernstein-like) waves in a plasma whose ion populations (positrons and fictitious singly-charged particles with twice the electron mass, for the sake of simulation feasability) can be time-varying nongyrotropic (TNG). The results show that the coexistence of two gyrophase bunched species does not significantly enhance the parallel SNG instability already found for media with only one nongyrotropic species, whereas it strongly intensifies the growth of Bernstein-like modes in TNG plasmas.     

Sausage mode instability of thin current sheets as a cause of magnetospheric substorms

Observations have shown that, prior to substorm explosions, thin current sheets are formed in the plasma sheet of the Earth’s magnetotail. This provokes the question, to what extent current-sheet thinning and substorm onsets are physically, maybe even causally, related. To answer this question, one has to understand the plasma stability of thin current sheets. Kinetic effects must be taken into account since particle scales are reached in the course of tail current-sheet thinning. We present the results of theoretical investigations of the stability of thin current sheets and about the most unstable mode of their decay. Our conclusions are based upon a non-local linear dispersion analysis of a cross-magnetic field instability of Harris-type current sheets. We found that a sausage-mode bulk current instability starts after a sheet has thinned down to the ion inertial length. We also present the results of three-dimensional electromagnetic PIC-code simulations carried out for mass ratios up to M_i/m_e = 64. They verify the linearly predicted properties of the sausage mode decay of thin current sheets in the parameter range of interest.     

A mode filter for plasma waves in the Hall-MHD approximation

A filter method is presented which allows a qualitative and quantitative identification of wave modes observed with plasma experiments on satellites. Hitherto existing mode filters are based on the MHD theory and thus they are restricted to low frequencies well below the ion cyclotron frequency. The present method is generalized to cover wave modes up to the characteristic ion frequencies. The spectral density matrix determined by the observations is decomposed using the eigenvectors of the linearized Hall-MHD equations. As the wave modes are dispersive in this formalism, a precise determination of the κ-vectors requires the use of multi-point measurements. Therefore the method is particularly relevant to multi-satellite missions. The method is tested using simulated plasma data. The Hall-MHD filter is able to identify the modes excited in the model plasma and to assign the correct energetic contributions. By comparison with the former method it is shown that the simple MHD filter leads to large errors if the frequency is not well below the ion cyclotron frequency. Further the range of validity of the linear theory is examined rising the simulated wave amplitudes.     

WIND observations of coherent electrostatic waves in the solar wind

The time domain sampler (TDS) experiment on WIND measures electric and magnetic wave forms with a sampling rate which reaches 120 000 points per second. We analyse here observations made in the solar wind near the Lagrange point L1. In the range of frequencies above the proton plasma frequency f_pi and smaller than or of the order of the electron plasma frequency f_pe, TDS observed three kinds of electrostatic (e.s.) waves: coherent wave packets of Langmuir waves with frequencies f ≃ f_pe, coherent wave packets with frequencies in the ion acoustic range f_pi ≥ f ≥ f_pe, and more or less isolated non-sinusoidal spikes lasting less than 1 ms. We confirm that the observed frequency of the low frequency (LF) ion acoustic wave packets is dominated by the Doppler effect: the wavelengths are short, 10 to 50 electron Debye lengths λ_D. The electric field in the isolated electrostatic structures (IES) and in the LF wave packets is more or less aligned with the solar wind magnetic field. Across the IES, which have a spatial width of the order of ≃25_D, there is a small but finite electric potential drop, implying an average electric field generally directed away from the Sun. The IES wave forms, which have not been previously reported in the solar wind, are similar, although with a smaller amplitude, to the weak double layers observed in the auroral regions, and to the electrostatic solitary waves observed in other regions in the magnetosphere. We have also studied the solar wind conditions which favour the occurrence of the three kinds of waves: all these e.s. waves are observed more or less continuously in the whole solar wind (except in the densest regions where a parasite prevents the TDS observations). The type (wave packet or IES) of the observed LF waves is mainly determined by the proton temperature and by the direction of the magnetic field, which themselves depend on the latitude of WIND with respect to the heliospheric current sheet.     

Heating of beam ions by ion acoustic waves

Satellite measurements show that ion beams above the auroral acceleration region are heated to hundreds of eV in a direction perpendicular to the magnetic field. We show that ion acoustic waves may be responsible for much of this heating. Even in the absence of a positive slope in the velocity distribution of the beam ions, ion acoustic waves can be generated by a fan instability. We present analytical estimates of the wave growth rate and ion beam heating rate. These estimates, which are confirmed by particle simulations, indicate that the perpendicular temperature of the beam ions will increase by 30 eV/s, or by 1 eV in 20–25 km. From the simulations we also conclude that the heating saturates at a perpendicular temperature around 200 eV, which is consistent with observations.     

Waveform and packet structure of lion roars

The Equator-S magnetometer is very sensitive and has a sampling rate of normally 128 Hz. The high sampling rate allows for the first time fluxgate magnetometer measurements of ELF waves between the ion cyclotron and the lower hybrid frequencies in the equatorial dayside magnetosheath. The so-called lion roars, typically seen by the Equator-S magnetometer at the bottom of the magnetic troughs of magnetosheath mirror waves, are near-monochromatic packets of electron whistler waves lasting for a few wave cycles only, typically 0.25 s. They are right-hand circularly polarized waves with typical amplitudes of 0.5–1 nT at around one tenth of the electron gyrofrequency. The cone angle between wave vector and ambient field is usually smaller than 1.5°.     

A generation mechanism for chorus emission

A chorus generation mechanism is discussed, which is based on interrelation of ELF/VLF noise-like and discrete emissions under the cyclotron wave-particle interactions. A natural ELF/VLF noise radiation is excited by the cyclotron instability mechanism in ducts with enhanced cold plasma density or at the plasmapause. This process is accompanied by a step-like deformation of the energetic electron distribution function in the velocity space, which is situated at the boundary between resonant and nonresonant particles. The step leads to the strong phase correlation of interacting particles and waves and to a new backward wave oscillator (BWO) regime of wave generation, when an absolute cyclotron instability arises at the central cross section of the geomagnetic trap, in the form of a succession of discrete signals with growing frequency inside each element. The dynamical spectrum of a separate element is formed similar to triggered ELF/VLF emission, when the strong wavelet starts from the equatorial plane. The comparison is given of the model developed using some satellite and ground-based data. In particular, the appearance of separate groups of chorus signals with a duration 2–10 s can be connected with the preliminary stage of the step formation. BWO regime gives a succession period smaller than the bounce period of energetic electrons between the magnetic mirrors and can explain the observed intervals between chorus elements.     

Parametric decay of beam-driven Langmuir wave and enhanced ion-acoustic fluctuations in the ionosphere: a weak turbulence approach

We present a model that describes the decay of beam generated Langmuir waves into ion-acoustic waves in the topside ionosphere. This calculation is done within the frame of the weak turbulence approximation. We study the spectral signature of such a process as seen by a VHF incoherent scatter radar. An incoherent scatter (IS) spectrum is characterized by two maxima at k_radar and −k_radar, the right and left ion lines respectively. It is shown that, for reasonable beam parameters, the parametric decay of beam-generated Langmuir waves can enhance either the right, the left or both ion lines simultaneously. The shape of the spectrum can change drastically on time scale of about 0.1 to 1 s. The role of the beam parameter as well as the ionospheric parameters is also investigated. For a given beam number density, the beam energy or the background density are important to trigger either the left or the right ion line. A large energy spread of the beam or low electron collision frequencies can explain the simultaneous observations of the left and the right ion line. The importance of the electron collision frequency can explain the altitude distribution of the coherent echoes observed by incoherent scatter radars.     

Concerning the generation of geomagnetic giant pulsations by drift-bounce resonance ring current instabilities

Giant pulsations are nearly monochromatic ULF-pulsations of the Earth’s magnetic field with periods of about 100 s and amplitudes of up to 40 nT. For one such event ground-magnetic observations as well as simultaneous GEOS-2 magnetic and electric field data and proton flux measurements made in the geostationary orbit have been analysed. The observations of the electromagnetic field indicate the excitation of an odd-mode type fundamental field line oscillation. A clear correlation between variations of the proton flux in the energy range 30–90 keV with the giant pulsation event observed at the ground is found. Furthermore, the proton phase space density exhibits a bump-on-the-tail signature at about 60 keV. Assuming a drift-bounce resonance instability as a possible generation mechanism, the azimuthal wave number of the pulsation wave field may be determined using a generalized resonance condition. The value determined in this way, m = −21±4, is in accord with the value m = −27±6 determined from ground-magnetic measurements. A more detailed examination of the observed ring current plasma distribution function f shows that odd-mode type eigenoscillations are expected for the case ∂f/∂W ≥ 0, much as observed. This result is different from previous theoretical studies as we not only consider local gradients of the distribution function in real space, but also in velocity space. It is therefore concluded that the observed giant pulsation is the result of a drift-bounce resonance instability of the ring current plasma coupling to an odd-mode fundamental standing wave. The generation of the bump-on-the-tail distribution causing ≥f/≥W ≥ 0 can be explained due to velocity dispersion of protons injected into the ring current. Both this velocity dispersion and the necessary substorm activity causing the injection of protons into the nightside magnetosphere are observed.     

Transport of thermal plasma above the auroral ionosphere in the presence of electrostatic ion-cyclotron turbulence

The electron component of intensive electric currents flowing along the geomagnetic field lines excites turbulence in the thermal magnetospheric plasma. The protons are then scattered by the excited electromagnetic waves, and as a result the plasma is stable. As the electron and ion temperatures of the background plasma are approximately equal each other, here electrostatic ion-cyclotron (EIC) turbulence is considered. In the nonisothermal plasma the ion-acoustic turbulence may occur additionally. The anomalous resistivity of the plasma causes large-scale differences of the electrostatic potential along the magnetic field lines. The presence of these differences provides heating and acceleration of the thermal and energetic auroral plasma. The investigation of the energy and momentum balance of the plasma and waves in the turbulent region is performed numerically, taking the magnetospheric convection and thermal conductivity of the plasma into account. As shown for the quasi-steady state, EIC turbulence may provide differences of the electric potential of δ V ≈ 1–10 kV at altitudes of 500 < h < 10 000 km above the Earth’s surface. In the turbulent region, the temperatures of the electrons and protons increase only a few times in comparison with the background values.     

Magnetospheric lion roars

The Equator-S magnetometer is very sensitive and has a sampling rate normally of 128 Hz. The high sampling rate for the first time allows detection of ELF waves between the ion cyclotron and the lower hybrid frequencies in the equatorial dawnside magnetosphere. The characteristics of these waves are virtually identical to the lion roars typically seen at the bottom of the magnetic troughs of magnetosheath mirror waves. The magnetospheric lion roars are near-monochromatic packets of electron whistler waves lasting for a few wave cycles only, typically 0.2 s. They are right-hand circularly polarized waves with typical amplitudes of 0.5 nT at around one tenth of the electron gyrofrequency. The cone angle between wave vector and ambient field is nearly always smaller than 1°.     

A case study of a radially polarized Pc4 event observed by the Equator-S satellite

A 16 mHz Pc4 pulsation was recorded on March 17, 1998, in the prenoon sector of the Earths magnetosphere by the Equator-S satellite. The event is strongly localized in radial direction at approximately L = 5 and exhibits properties of a field line resonance such as an ellipticity change as seen by applying the method of the analytical signal to the magnetic field data. The azimuthal wave number was estimated as m 150. We discuss whether this event can be explained by the FLR mechanism and find out that the change in ellipticity is more a general feature of a localized Alfvén wave than indicative of a resonant process.     

Study of the magnetic turbulence in a corotating interaction region in the interplanetary medium

We study the geometry of magnetic fluctuations in a CIR observed by Pioneer 10 at 5 AU between days 292 and 295 in 1973. We apply the methodology proposed by Bieber et al. to make a comparison of the relative importance of two geometric arrays of vector propagation of the magnetic field fluctuations: slab and two-dimensional (2D). We found that inside the studied CIR this model is not applicable due to the restrictions imposed on it. Our results are consistent with Alfvenic fluctuations propagating close to the radial direction, confirming Mavromichalaki et al.’s findings. A mixture of isotropic and magnetoacoustic waves in the region before the front shock would be consistent with our results, and a mixture of slab/2D and magnetoacoustic waves in a region after the reverse shock. We base the latter conclusions on the theoretical analysis made by Kunstmann. We discuss the reasons why the composite model can not be applied in the CIR studied although the fluctuations inside it are two dimensional.     

A simple oblique dip model for geomagnetic micropulsations

It is pointed out that simple models adopted so far have tended to neglect the obliquity of the magnetic field lines entering the Earths surface. A simple alternative model is presented, in which the ambient field lines are straight, but enter wedge shaped boundaries at half a right-angle. The model is illustrated by assuming an axially symmetric, compressional, impulse type disturbance at the outer boundary, all other boundaries being assumed to be perfectly conducting. The numerical method used is checked from the instant the excitation ceases, by an analytical method. The first harmonic along field lines is found to be of noticeable size, but appears to be mainly due to coupling with the fundamental, and with the first harmonic across field lines.     

Observations of the upper frequency cutoffs of the auroral kilometric radiation

Intense auroral kilometric radiation (AKR) is being frequently observed with POLRAD from the Auroral Probe (Interball-2). Observations of the abrupt upper frequency cutoffs (UFCs) in the spectra of AKR are reported. The UFCs can be observed at a frequency range from 300 to 700 kHz, corresponding to AKR generation altitudes from approximately 4800 to 2100 km, and are distributed in magnetic local time (MLT) hours similarly to the AKR events, with a maximum at 1 h MLT. The observed frequency extent of the UFCs is 12 kHz, and is often determined by the instrumental resolution (4 kHz). It is suggested that the UFC may be associated with an abrupt switching on of the generation mechanism, when the electron density becomes sufficiently low inside a plasma depletion at an altitude where the ratio of f_pe/f_ce crosses some threshold value. The steepness of the UFCs can imply a non-linear process of generation. The estimated distance of the e-folding field aligned wave amplification is between 3 and 8 km. The UFCs are sometimes, though very seldom (10%), accompanied by narrow band (less than 4 kHz) ridges of radiation observed at the cutoff frequency. They are smoothly drifting in frequncy for several minutes. The power density of radiation in the ridge can be up to 2 orders of magnitude stronger than in the accompanying wide band emission of AKR. The ridge at UFC can imply either energy concentration at the source bottom, or focusing, if specific conditions for the escape of the radiation are assumed.     

High-beta plasma blobs in the morningside plasma sheet

Equator-S frequently encountered, i.e. on 30%0of the orbits between 1 March and 17 April 1998, strong variations of the magnetic field strength of typically 5–15-min duration outside about 9R_E during the late-night/early-morning hours. Very high-plasma beta values were found, varying between 1 and 10 or more. Close conjunctions between Equator-S and Geotail revealed the spatial structure of these “plasma blobs” and their lifetime. They are typically 5–10° wide in longitude and have an antisymmetric plasma or magnetic pressure distribution with respect to the equator, while being altogether low-latitude phenomena (<15°). They drift slowly sunward, exchange plasma across the equator and have a lifetime of at least 15–30 min. While their spatial structure may be due to some sort of mirror instability, little is known about the origin of the high-beta plasma. It is speculated that the morningside boundary layer somewhat further tailward may be the source of this plasma. This would be consistent with the preference of the plasma blobs to occur during quiet conditions, although they are also found during substorm periods. The relation to auroral phenomena in the morningside oval is uncertain. The energy deposition may be mostly too weak to generate a visible signature. However, patchy aurora remains a candidate for more disturbed periods.     

Nonstationary pearl pulsations as a signature of magnetospheric disturbances

We analyse long-lasting (several hours) Pc1 pearl pulsations with decreasing, increasing or constant central frequencies. We show that nonstationary pearl events (those with either decreasing or increasing central frequency) are observed simultaneously with increasing auroral magnetic activity at the nightside magnetosphere while the stationary events (constant central frequency) correspond to quiet magnetic conditions. Events with decreasing central frequency are observed mostly in the late morning and daytime whereas events with increasing central frequency appear either early in the morning or in the afternoon. We explain the diurnal distribution of the nonstationary pearl pulsations in terms of proton drifts depending on magnetic activity, and evaluate the magnetospheric electric field based on the variation of the central frequency of pearl pulsations.     

On the nonlinear triggering of VLF emissions by power line harmonic radiation

VLF ground data from Porojarvi in N. Finland has been presented. Spectrograms reveal frequent occurrence of power line harmonic radiation (PLHR), originating from the Finnish power system and from heavy industrial plant. The radiation is seen to penetrate the magnetosphere since numerous occurrences of PLHR triggered emissions are seen. Risers predominate but fallers and hooks are also observed. A well-established 1D Vlasov simulation code has been used to simulate these emissions, using plausible magnetospheric data for a range of L values from L = 4 to L = 5.5. The code is able to reproduce risers fallers and hooks in close agreement with observations. The results shed considerable insight into the generation structure of both risers and fallers.