Particle energy fluxes in an unstable plasma with vortex structures in an inhomogeneous geomagnetic field in the topside ionosphere

Plasma inhomogeneities extending along geomagnetic field lines in the ionosphere and magnetosphere can have a vortex structure. Electromagnetic waves can propagate in plasma inhomogeneities in the waveguide channel mode. It has been indicated that energy and particle fluxes related to the development of small-scale electrostatic turbulence in a magnetized plasma with an unstable electron component promotes an increase in plasma density gradients in the walls of waveguide channels and an enhancement in plasma vortices. At low L shells in the region of the geomagnetic equator, the development of plasma electrostatic instability and the damping of drifting plasma vortices in the inhomogeneous geomagnetic field in the topside ionosphere can be the main mechanism by which large-scale (∼1000 km) regions with a decreased plasma density are formed.     

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.     

Effect of magnetospheric convection on the energy distribution of protons from the Earth radiation belts

The structure of energetic protons from the Earth radiation belts, averaged for a magnetically quiet period, can be explained by the equilibrium between the radial diffusion transfer, loss due to Coulomb collisions, charge exchange with ambient neutral hydrogen of the geocorona, and drift of protons under the influence of magnetospheric convection. By transfer we mean diffusion owing to fluctuations related to substorms in the large-scale electric and magnetic fields. Equatorially mirroring protons with energies of 1–750 keV have been considered, and the theoretical predictions of the proton energy spectra for L = 1.0−6.6 have been compared with the observations on several satellites.     

High-latitude HF Doppler observations of ULF waves. 1. Waves with large spatial scale sizes

A quantitative study of observations of the ionospheric signatures of magnetospheric ultra low frequency (ULF) waves by a high-latitude (geographic: 69.6°N 19.2°E) high-frequency Doppler sounder has been undertaken. The signatures, which are clearly correlated with pulsations in ground magnetometer data, exhibit periods in the range 100–400 s and have azimuthal wave numbers in the range 3–8. They are interpreted here as local field line resonances. Phase information provided by O- and X-mode Doppler data support the view that these are associated with field line resonances having large azimuthal scale sizes. The relative phases and amplitudes of the signatures in the Doppler and ground magnetometer data are compared with a model for the generation of Doppler signatures from incident ULF waves. The outcome suggests that the dominant mechanism involved in producing the Doppler signature is the vertical component of an E × B bulk motion of the local plasma caused by the electric field perturbation of the ULF wave.     

The dynamic cusp at low altitudes: a case study utilizing Viking,DMSP-F7, and Sondrestrom incoherent scatter radar observations

Coincident multi-instrument magnetospheric and ionospheric observations have made it possible to determine the position of the ionospheric footprint of the magnetospheric cusp and to monitor its evolution over time. The data used include charged particle and magnetic field measurements from the Earth-orbiting Viking and DMSP-F7 satellites, electric field measurements from Viking, interplanetary magnetic field and plasma data from IMP-8, and Sondrestrom incoherent scatter radar observations of the ionospheric plasma density, temperature, and convection. Viking detected cusp precipitation poleward of 75.5○ invariant latitude. The ionospheric response to the observed electron precipitation was simulated using an auroral model. It predicts enhanced plasma density and elevated electron temperature in the upper E- and F-regions. Sondrestrom radar observations are in agreement with the predictions. The radar detected a cusp signature on each of five consecutive antenna elevation scans covering 1.2 h local time. The cusp appeared to be about 2○ invariant latitude wide, and its ionospheric footprint shifted equatorward by nearly 2○ during this time, possibly influenced by an overall decrease in the IMF B_z component. The radar plasma drift data and the Viking magnetic and electric field data suggest that the cusp was associated with a continuous, rather than a patchy, merging between the IMF and the geomagnetic field.     

Turbulent mixing in the Northern North Sea: a numerical model study

The aim of this study is to intervalidate observations and numerical simulation results for the turbulent dissipation rate under strong wind conditions in the Northern North Sea during one week in October 1998. The observations were obtained by spatially and temporally averaging measurements of small-scale shear with a free-falling shear probe. The 1D numerical model used for this study is based on a state-of-the-art two-equation k− turbulence model with an algebraic second-moment closure scheme. It is discussed by means of annual and seasonal model simulations how the influence of heat and salt advection and internal waves can be accounted for. After these precautions, the agreement between observations and simulations of the turbulent dissipation rate are fairly good. Remaining differences cannot only be explained by problems such as undersampling and noise level, but also by idealising model assumptions.     

ULF wave measurements onboard the Interball auroral probe

The IESP experiment implemented onboard the Interball auroral probe measures the six components (3B, 3E) of the waves in the ULF range: 0.1–10 Hz and from time to time 0–30 Hz. Two different kinds of waves have been observed in the auroral region at altitudes between 10 000 and 20 000 km: (1) electrostatic emissions which consist of quasi-monochromatic structures with frequencies above the oxygen gyrofrequency, superimposed on a wide band signal interpreted as a Doppler broadening, (2) electromagnetic wide band spectrum fluctuations. These emissions are interpreted as current-driven electromagnetic or electrostatic ion cyclotron waves. The electromagnetic/electrostatic character is controlled by the plasma parameter _i and by the O⁺ concentration.     

On the current-voltage relationship in auroral breakups and westwards-travelling surges

Auroral precipitating electrons pass through an acceleration region before entering the atmosphere. Regardless of what produces it, a parallel electric field is assumed to cause the acceleration. It is well known that from kinetic theory an expression for the corresponding upward field-aligned current can be calculated, which under certain assumptions can be linearized to j = KV. The K constant, referred to as the Lyons-Evans-Lundin constant, depends on the source density and thermal energy of the magnetospheric electrons; it is an important parameter in magnetosphere-ionosphere coupling models. However, the K parameter is still rather unknown, and values are found in a wide range of 10^–8–10^–10Sim^–2. In this study, we investigated how the type of auroral structure affects the K values. We look at onset and westwards-travelling surge (WTS) events and make comparisons with earlier results from observations of more stable auroral arcs. A new analysis technique for studying those magnetospheric parameters using ground-based measurements is introduced. Electron density measurements are taken with the EISCAT radar, and through an inversion technique the flux-energy spectra are calculated. Source densities, thermal energies and potential drops are estimated from fittings of accelerated Maxwellian distributions. With this radar technique we have the possibility to study the changes of the mentioned parameters during the development of onsets and the passage of surges over EISCAT. The study indicates that the linearization of the full Knight formulation holds even for the very high potential drops and thermal temperatures found in the dynamic onset and WTS events. The values of K are found to be very low, around 10^–11Sm^–2 in onset cases as well as WTS events. The results may establish a new technique where ionospheric measurements are used for studying the ionosphere-magnetosphere coupling processes.     

On the role of collective interactions in asymmetric ring current formation

The contribution of resonant wave-particle interactions to the formation and decay of the magnetospheric ring current is analysed in the framework of a self-consistent set of equations which take into account azimuthal plasmasphere asymmetry. It is shown that the cyclotron interaction of westward drifting energetic protons with Alfven waves in the evening-side plasmaspheric bulge region leads to the formation of a ring current asymmetry located near 18:00 MLT. The time-scale of this asymmetry is determined by the proton drift time through the plasmaspheric bulge and is about 1 - 3 h. A symmetrical ring current decays mainly due to charge exchange processes. The theory is compared with known experimental data on ions and waves in the ring current and on low-latitude magnetic disturbances. New low-latitude magnetometer data on the magnetic storm of 24 - 26 July 1986 are also discussed. The model presented explains the observed localization of an asymmetrical ring current loop in the evening sector and the difference in relaxation time-scales of the asymmetry and the D_st index. It also explains measured wave turbulence levels in the evening-side plasmasphere and wave observation statistics.     

Excitation of ion-acoustic waves in the high-latitude ionosphere

The broadband electrostatic turbulence generally observed in the high-latitude ionosphere is a superposition of nonlocal waves of ion-acoustic and ion-cyclotron types. In the presence of a shear of ion parallel velocity, ion-acoustic modes can be induced by an instability emerging due to an inhomogeneous distribution of energy density. This paper is devoted to the studies of excitation of oblique ion-acoustic wave in background configurations with inhomogeneous profiles of both electric field and ion parallel velocity. A numerical algorithm has been developed, and instability was simulated at various parameters of background plasma. The general possibility of oblique ion-acoustic wave generation by a gradient of ion parallel velocity is shown. In this case, the wave spectrum is found to be broadband, which agrees with satellite observations.     

Magnetosheath Interaction with the High Latitude Magnetopause

We present both statistical and case studies of magnetosheath interaction with the high-latitude magnetopause on the basis of Interball-1 and other ISTP spacecraft data. We discuss those data along with recently published results on the topology of cusp-magnetosheath transition and the roles of nonlinear disturbances in mass and energy transfer across the high-latitude magnetopause. For sunward dipole tilts, a cusp throat is magnetically open for direct interaction with the incident flow that results in the creation of a turbulent boundary layer (TBL) over an indented magnetopause and downstream of the cusp. For antisunward tilts, the cusp throat is closed by a smooth magnetopause; demagnetized ‘plasma balls’ (with scale ∼ few R_E, an occurrence rate of ∼25% and trapped energetic particles) present a major magnetosheath plasma channel just inside the cusp. The flow interacts with the ‘plasma balls’ via reflected waves, which trigger a chaotization of up to 40% of the upstream kinetic energy. These waves propagate upstream of the TBL and initiate amplification of the existing magnetosheath waves and their cascade-like decays during downstream passage throughout the TBL. The most striking feature of the nonlinear interaction is the appearance of magnetosonic jets, accelerated up to an Alfvenic Mach number of 3. The characteristic impulsive local momentum loss is followed by decelerated Alfvenic flows and modulated by the TBL waves; momentum balance is conserved only on time scales of the Alfvenic flows (1/f_A ∼12 min). Wave trains at f_A∼1.3 mHz are capable of synchronizing interactions throughout the outer and inner boundary layers. The sonic/Alfvenic flows, bounded by current sheets, control the TBL spectral shape and result in non-Gaussian statistical characteristics of the disturbances, indicating the fluctuation intermittency. We suggest that the multi-scale TBL processes play at least a comparable role to that of macro-reconnection (remote from or in the cusp) in solar wind energy transformation and population of the magnetosphere by the magnetosheath plasma. Secondary micro-reconnection constitutes a necessary chain at the small-scale (∼ion gyroradius) edge of the TBL cascades. The thick TBL transforms the flow energy, including deceleration and heating of the flow in the open throat, ‘plasma ball’ and the region downstream of the cusp.     

Simulation studies of processes associated with stimulated electromagnetic emission (SEE) in the ionosphere

Results of numerical simulation studies of processes associated with Stimulated Electromagnetic Emission (SEE) produced during ionospheric heating experiments are presented. A one-dimensional magnetized electrostatic Particle-In-Cell (PIC) simulation model with uniform plasma density is used to investigate electrostatic wave generation in the region where the pump frequency ω₀ approximately equals the upper hybrid frequency ω_uh. In particular, the simulation plasma is driven with a uniform oscillating electric field to represent the long wavelength pump wave and power spectra of the electrostatic waves produced are taken. The pump wave frequency and amplitude are varied to consider the effects on the simulation power spectrum. The upper hybrid frequency in the model is varied through harmonics of the electron cyclotron frequency Ω_ce to consider the effects of stepping the pump frequency through cyclotron harmonics. The power spectrum from the simulation plasma is richly structured. The resulting power spectra show sidebands upshifted and downshifted from the pump frequency by multiples of the lower hybrid frequency ω_lh. The structure of the spectrum is highly sensitive to the proximity of the upper hybrid frequency to the cyclotron harmonic frequencies.     

High-latitude HF Doppler observations of ULF waves: 2. Waves with small spatial scale sizes

The DOPE (Doppler Pulsation Experiment) HF Doppler sounder located near Tromsø, Norway (geographic: 69.6°N 19.2°E; L = 6.3) is deployed to observe signatures, in the high-latitude ionosphere, of magnetospheric ULF waves. A type of wave has been identified which exhibits no simultaneous ground magnetic signature. They can be subdivided into two classes which occur in the dawn and dusk local time sectors respectively. They generally have frequencies greater than the resonance fundamentals of local field lines. It is suggested that these may be the signatures of high-m ULF waves where the ground magnetic signature has been strongly attenuated as a result of the scale size of the waves. The dawn population demonstrate similarities to a type of magnetospheric wave known as giant (Pg) pulsations which tend to be resonant at higher harmonics on magnetic field lines. In contrast, the waves occurring in the dusk sector are believed to be related to the storm-time Pc5s previously reported in VHF radar data. Dst measurements support these observations by indicating that the dawn and dusk classes of waves occur respectively during geomagnetically quiet and more active intervals.     

Weimer电场模式在地球磁层内的特征

地球磁层中的电场是研究磁层物理的重要参数,目前常用的对流电场有均匀晨昏电场和投影电场.电离层电场可以看做磁层电场沿磁力线在电离层的投影,本文选取的电离层电场模型为Weimer(2001模式)电场.利用T96磁场模式,沿磁力线将电离层电场投影到磁层空间,得到一个新的磁层电场模式,并讨论了磁暴、行星际磁场(IMF)、太阳风参数和亚暴等对磁层电场的影响.利用该模型计算的电场结果与卫星探测结果相符.     

Nonlinear wave structures in collisional plasma of auroral E-region ionosphere

Studies of the auroral plasma with small-scale inhomogenieties producing the VHF-radar reflections (radar aurora) when observed in conditions of the saturated Farley-Buneman instability within the auroral E region, show strong nonlinear interactions and density fluctuations of 5–15%. Such nonlinearity and high fluctation amplitudes are inconsistent with the limitations of the weak turbulence theory, and thus a theory for arbitrary amplitudes is needed. To this end, a nonlinear theory is described for electrostatic MHD moving plasma structures of arbitrary amplitude for conditions throughout the altitude range of the collisional auroral E region. The equations are derived, from electron and ion motion selfconsistent with the electric field, for the general case of the one-dimensional problem. They take into account nonlinearity, electron and ion inertia, diffusion, deviation from quasi-neutrality, and dynamical ion viscosity. The importance of the ion viscosity for dispersion is stressed, while deviation from the quasi-neutrality can be important only at rather low plasma densities, not typical for the auroral E region. In a small amplitude limit these equations have classical nonlinear solutions of the type of electrostatic shock wave or of knoidal waves. In a particular case these knoidal waves degrade to a dissipative soliton. A two-dimensional case of a quasi-neutral plasma is considered in the plane perpendicular to the magnetic field by way of the Poisson brackets, but neglecting the nonlinearity and ion inertia. It is shown that in these conditions an effective saturation can be achieved at the stationary turbulence level of order of 10%.     

Properties and origin of energetic particles at the duskside of the Earth’s magnetosheath throughout a great storm

We study an interval of 56 h on January 16 to 18, 1995, during which the GEOTAIL spacecraft traversed the duskside magnetosheath from X ≅ −15 to −40 R_E and the EPIC/ICS and EPIC/STICS sensors sporadically detected tens of energetic particle bursts. This interval coincides with the expansion and growth of a great geomagnetic storm. The flux bursts are strongly dependent on the magnetic field orientation. They switch on whenever the B_z component approaches zero (B_z ≅ 0 nT). We strongly suggest a magnetospheric origin for the energetic ions and electrons streaming along these “exodus channels”. The time profiles for energetic protons and “tracer” O⁺ ions are nearly identical, which suggests a common source. We suggest that the particles leak out of the magnetosphere all the time and that when the magnetosheath magnetic field connects the spacecraft to the magnetotail, they stream away to be observed by the GEOTAIL sensors. The energetic electron fluxes are not observed as commonly as the ions, indicating that their source is more limited in extent. In one case study the magnetosheath magnetic field lines are draped around the magnetopause within the YZ plane and a dispersed structure for peak fluxes of different species is detected and interpreted as evidence for energetic electrons leaking out from the dawn LLBL and then being channelled along the draped magnetic field lines over the magnetopause. Protons leak from the equatorial dusk LLBL and this spatial differentiation between electron and proton sources results in the observed dispersion. A gradient of energetic proton intensities toward the Z_GSM= 0 plane is inferred. There is a permanent layer of energetic particles adjacent to the magnetosheath during this interval in which the dominant component of the magnetic field was B_z.     

Intermittency in midlatitude Pi2 pulsations observed during magnetospheric substorms and in the absence of these phenomena

The results of studying the distribution character of the amplitudes and time intervals between wave packets of Pi2 geomagnetic pulsations, observed during the nighttime development of magnetospheric substorms and in the absence of these phenomena, have been presented. The analog records from the midlatitude Borok Geophysical Observatory (geographic coordinates φ = 58.03°; λ = 38.97°) for 1995 and 1997 have been used to analyze Pi2 pulsations. Three groups of pulsations have been analyzed: Pi2 observed during sub-storms related to the external impact on the magnetosphere, Pi2 spontaneously originating during substorms, and Pi2 observed in the absence of substorms on the nightside of the magnetosphere. Interplanetary magnetic field B _y and B _z components and the solar wind dynamic pressure (ρV ²) have been considered as possible triggers of magnetospheric substorms. It has been indicated that the distributions of the amplitude and the duration of time intervals between Pi2 bursts are approximated by the power and exponential functions, respectively, which is typical of intermittent processes. The hypothesis that the processes of magnetospheric plasma turbulization can be among the Pi2 pulsation burst sources has been put forward. It is assumed that the obtained characteristic values can be used to qualitatively estimate the degree of plasma turbulence on the nightside of the magnetosphere when a sequence of Pi2 wave packets is excited.     

大功率高频电波加热电离层中参量衰减不稳定性研究(1):数值模拟

参量衰减不稳定性(Parametric Decay Instability,PDI)在大功率高频(High frequency,HF)电波与电离层等离子体相互作用的过程中扮演着十分重要的角色,本文采用广义Zakharov方法对常规的等离子体流体力学方程组进行相应处理后,并在近似实际的电离层背景和电波传播模型下,构建了高频电波加热电离层激发PDI的数值计算模型.模拟结果发现:在毫秒量级的时间尺度上,大功率高频电波在寻常波(Ordinary wave,O波)反射点高度附近激发出了朗缪尔波(Langmuir wave)和离子声波(Ion-Acoustic wave)两种等离子体静电波模,模拟中产生的朗缪尔波和离子声波相应波数为5~11rad·m~(-1),结果与利用色散关系求出的理论值4~7 rad·m~(-1)近似一致,密度扰动幅值从10~6m~(-3)量级指数级增长到了10~(10)m~(-3)量级,直至能显著影响与"低频"密度背景相关的等离子体频率后,出现了等离子体"空穴"结构以及朗缪尔波被"俘获"现象,在扰动空间内的小尺度静电场幅值最高能达到100 V·m~(-1)量级,最终造成一种强烈的局地化"空穴"湍流现象.本文的研究有助于深入理解PDI的物理机制,对研究大功率高频电波与电离层等离子体之间复杂的非线性相互作用也有着非常重要的意义.     

Auroral electrojet dynamics during magnetic storms, connection with plasma precipitation and large-scale structure of the magnetospheric magnetic field

Effect of the equatorward shift of the eastward and westward electrojets during magnetic storms main phase is analyzed based on the meridional chains of magnetic observatories EISCAT and IMAGE and several Russian observatories (geomagnetic longitude ≈110°, corrected geomagnetic latitudes 74°>φ>51°.) Magnetic storms of various Dst index intensity where the main phase falls on 1000 UT - 2400 UT interval were selected so that one of the observatory chains was located in the afternoon - near midnight sector of MLT. The eastward electrojet center shifts equatorward with Dst intensity increase: when Dst ≈ −50 nT the electrojet center is located at φ ≈ 62°, when Dst ≈ −300 nT it is placed at φ ≈ 54°. The westward electrojet center during magnetic storms main phase for intervals between substorms shifts equatorward with Dst increase: at φ ≈ 62° when Dst ≈ −100 nT and at φ ≈ 55° when Dst ≈ −300 nT. During substorms within the magnetic storms intervals the westward electrojet widens poleward covering latitudes φ ≈ 64°–65°. DMSP (F08, F10 and F11) satellite observations of auroral energy plasma precipitations at upper atmosphere altitudes were used to determine precipitation region structure and location of boundaries of various plasma domains during magnetic storms on May 10–11, 1992, February 5–7 and February 21–22, 1994. Interrelationships between center location, poleward and equatorward boundaries of electrojets and characteristic plasma regions are discussed. The electrojet center, poleward and equatorward boundaries along the magnetic observatories meridional chain were mapped to the magnetosphere using the geomagnetic field paraboloid model. The location of auroral energy oxygen ion regions in the night and evening magnetosphere is determined. Considerations are presented on the mechanism causing the appearance in the inner magnetosphere during active intervals of magnetic storms of ions with energy of tens KeV. In the framework of the magnetospheric magnetic field paraboloid model the influence of the ring current and magnetospheric tail plasma sheet currents on large-scale magnetosphere structure is considered.     

Geospheric effects of the solar flare of December 13, 2006

Disturbances in the magnetic field and magnetospheric and ionospheric plasma registered on December 14–16, 2006, during a strong magnetic storm caused by a solar flare of 4B/X3.4 class are studied. It is shown that in the north of Yakutia, interactions between the Earth’s magnetosphere and the region of high dynamic pressure of the solar wind led to the formation of sporadic layers in the ionospheric E and F regions, depletion of the critical frequency of the F2 layer, and total absorption. At the end of the magnetic storm’s main phase, anomalously high values of foF2 exceeding the quiet level by a factor of 1.5–1.7 were detected. It was found that the disturbances detected by ground-based observatories had developed on the background of changes in the temperature, density, and the pitch-angle distribution of particles at the geostationary orbit manifesting radial shifts of magnetospheric structures (magnetopause, cusp/cleft, and plasma sheet) relative to the observation points. A conclusion is drawn that in this case, changes in the near-Earth plasma and magnetic field manifest the dynamics of the physical conditions at the magnetospheric boundary and diurnal rotation of the Earth.