Zakharov simulations of beam-induced turbulence in the auroral ionosphere

Recent detections of strong incoherent scatter radar echoes from the auroral F region, which have been explained as the signature of naturally produced Langmuir turbulence, have motivated us to revisit the topic of beam-generated Langmuir turbulence via simulation. Results from one-dimensional Zakharov simulations are used to study the interaction of ionospheric electron beams with the background plasma at the F region peak. A broad range of beam parameters extending by more than 2 orders of magnitude in average energy and electron number density is considered. A range of wave interaction processes, from a single parametric decay, to a cascade of parametric decays, to formation of stationary density cavities in the condensate region, and to direct collapse at the initial stages of turbulence, is observed as we increase the input energy to the system. The effect of suprathermal electrons, produced by collisional interactions of auroral electrons with the neutral atmosphere, on the dynamics of Langmuir turbulence is also investigated. It is seen that the enhanced Landau damping introduced by the suprathermal electrons significantly weakens the turbulence and truncates the cascade of parametric decays

Do magnetospheric shear Alfvén waves generate sufficient electron energy flux to power the aurora?

Using a self-consistent drift-kinetic simulation code, we investigate whether electron acceleration owing to shear Alfvén waves in the plasma sheet boundary layer is sufficient to cause auroral brightening in the ionosphere. The free parameters used in the simulation code are guided by in situ observations of wave and plasma parameters in the magnetosphere at distances >4 RE from the Earth. For the perpendicular wavelength used in the study, which maps to ∼4 km at 110 km altitude, there is a clear amplitude threshold which determines whether magnetospheric shear Alfvén waves above the classical auroral acceleration region can excite sufficient electrons to create the aurora. Previous studies reported wave amplitudes that easily exceed this threshold; hence, the results reported in this paper demonstrate that auroral acceleration owing to shear Alfvén waves can occur in the magnetosphere at distances >4 RE from the Earth

Ionospheric Plasma Transport and Loss in Auroral Downward Current Regions

A detailed study of the effects of auroral current systems on thermal ionospheric plasma transport and loss is conducted using a new ionospheric model. The mathematical formulation of the model is a variation on the 5‐moment approximation which describes the temporal evolution of density, drift, and temperature for five different ion species in two spatial dimensions. The fluid system is closed through a 2‐D electrostatic treatment of the auroral currents. This model is used to examine the interplay between ion heating, perpendicular transport, molecular ion generation, and type‐1 ion upflows in a self‐consistent way for the first time. Simulations confirm that the depletion of E‐region plasma due to current closure occurs on extremely fast time scales (5–30 s), and that it is dependent on current system scale size. Near the F‐region peak, the loss is mostly due to enhanced recombination from the conversion of the plasma to molecular ions. The F‐region loss process is fairly slow (120–300 s) by comparison to lower altitude processes and is highly electric field dependent. On similar time scales, transient ion upflows from frictional heating move plasma from the near topside ionosphere (∼500 km) to higher regions, leaving depletions and enhancing plasma densities at very high altitudes. Results indicate the existence of large molecular ion upflows near the F‐region peak and may shed some light on ionospheric source regions for outflowing molecular ions. Neutral atmospheric winds and densities are also shown to play an important role in modulating molecular ion densities, frictional heating, and currents

On the Source of Energetic Electron Precipitation during Auroral Substorms

Precipitating auroral electrons are believed to originate mainly from parallel electric fields set up at the auroral acceleration region (AAR) extending up to 20,000 km altitude. However, electrons of energy greater than 100 keV are probably generated by acceleration processes beyond the AAR. Observational evidence for the source location of these energetic electrons are hard to come by. In our current work, we present simultaneous magnetically conjugate measurements of energetic electron spectra estimated at the ionosphere using the Poker Flat Incoherent Scatter Radar (PFISR) and measured at the inner plasma sheet by the THEMIS spacecraft. The flux of precipitating electrons of energy greater than 100 keV demonstrate a striking spatio-temporal correlation with that of the inner plasma sheet electrons. This suggests that the source of the energetic electrons lie at or beyond the inner plasma sheet, and that the acceleration processes within the auroral acceleration zone don't contribute substantially to their energization. Using simultaneous THEMIS measurements of wave power, we speculate that the electromagnetic ion cyclotron (EMIC) and Chorus waves are likely candidates for electron acceleration within the inner plasma sheet apart from the usual candidates of betatron and fermi acceleration. However, between the ionosphere and the plasma sheet, electrons of energy less than 100 keV show significant differences in their energy spectra after the substorm onset suggesting an active AAR

A diagnosis of the plasma waves responsible for the explosive energy release of substorm onset

During geomagnetic substorms, stored magnetic and plasma thermal energies are explosively converted into plasma kinetic energy. This rapid reconfiguration of Earth’s nightside magnetosphere is manifest in the ionosphere as an auroral display that fills the sky. Progress in understanding of how substorms are initiated is hindered by a lack of quantitative analysis of the single consistent feature of onset; the rapid brightening and structuring of the most equatorward arc in the ionosphere. Here, we exploit state-of-the-art auroral measurements to construct an observational dispersion relation of waves during substorm onset. Further, we use kinetic theory of high-beta plasma to demonstrate that the shear Alfven wave dispersion relation bears remarkable similarity to the auroral dispersion relation. In contrast to prevailing theories of substorm initiation, we demonstrate that auroral beads seen during the majority of substorm onsets are likely the signature of kinetic Alfven waves driven unstable in the high-beta magnetotail

Dynamics of Density Cavities Generated by Frictional Heating: Formation, Distortion, and Instability

A simulation study of the generation and evolution of mesoscale density cavities in the polar ionosphere is conducted using a time-dependent, nonlinear, quasi-electrostatic model. The model demonstrates that density cavities, generated by frictional heating, can form in as little as 90 s due to strong electric fields of ∼120 mV/m, which are sometimes observed near auroral zone and polar cap arcs. Asymmetric density cavity features and strong plasma density gradients perpendicular to the geomagnetic field are naturally generated as a consequence of the strong convection and finite extent of the auroral feature. The walls of the auroral density cavities are shown to be susceptible to large-scale distortion and gradient-drift instability, hence indicating that arc-related regions of frictional heating may be a source of polar ionospheric density irregularities

Incoherent Scatter Radar Estimation of F Region Ionospheric Composition During Frictional Heating Events

A method is developed for estimating F region ion composition from incoherent scatter radar (ISR) measurements during times of frictional ion heating. The technique addresses ion temperature‐mass ambiguities in the IS spectra by self‐consistently modeling ion temperature profiles, including the effects of ion temperature anisotropies and altitude‐independent neutral winds. The modeled temperature profiles are used in a minimization procedure to estimate ion composition consistent with the recorded IS spectra. The proposed method is applicable to short‐integration (min) data sets from either single‐beam or multiple‐beam experiments. Application of the technique to Sondrestrom ISR measurements shows increases in F region molecular ions in response to frictional heating, a result consistent with previous theoretical and observational work. Estimates of ion composition are shown to be relatively insensitive to moderate variations in the neutral atmospheric model, which serves as input to the method. The technique developed in this work is uniquely qualified for studying highly variable ion composition near auroral arcs and associated processes such as molecular ion upflows. It also addresses a systematic source of error in standard ISR analysis methods when they are applied in such situations

The Optical Manifestation of Dispersive Field‐Aligned Bursts in Auroral Breakup Arcs

High‐resolution optical observations of a substorm expansion show dynamic auroral rays with surges of luminosity traveling up the magnetic field lines. Observed in ground‐based imagers, this phenomenon has been termed auroral flames, whereas the rocket signatures of the corresponding energy dispersions are more commonly known as field‐aligned bursts. In this paper, observations of auroral flames obtained at 50 frames/s with a scientific‐grade Complementary Metal Oxide Semiconductor (CMOS) sensor (30° × 30° field of view, 30 m resolution at 120 km) are used to provide insight into the nature of the precipitating electrons similar to high‐resolution particle detectors. Thanks to the large field of view and high spatial resolution of this system, it is possible to obtain a first‐order estimate of the temporal evolution in altitude of the volume emission rate from a single sensor. The measured volume emission rates are compared with the sum of modeled eigenprofiles obtained for a finite set of electron beams with varying energy provided by the TRANSCAR auroral flux tube model. The energy dispersion signatures within each auroral ray can be analyzed in detail during a fraction of a second. The evolution of energy and flux of the precipitation shows precipitation spanning over a large range of energies, with the characteristic energy dropping from 2.1 keV to 0.87 keV over 0.2 s. Oscillations at 2.4 Hz in the magnetic zenith correspond to the period of the auroral flames, and the acceleration is believed to be due to Alfvenic wave interaction with electrons above the ionosphere

Reconstruction of Fine Scale Auroral Dynamics

We present a feasibility study for a high frame rate, short baseline auroral tomographic imaging system useful for estimating parametric variations in the precipitating electron number flux spectrum of dynamic auroral events. Of particular interest are auroral substorms, characterized by spatial variations of order 100 m and temporal variations of order 10 ms. These scales are thought to be produced by dispersive Alfvén waves in the near-Earth magnetosphere. The auroral tomography system characterized in this paper reconstructs the auroral volume emission rate to estimate the characteristic energy and location in the direction perpendicular to the geomagnetic field of peak electron precipitation flux using a distributed network of precisely synchronized ground-based cameras. As the observing baseline decreases, the tomographic inverse problem becomes highly ill-conditioned; as the sampling rate increases, the signal-to-noise ratio degrades and synchronization requirements become increasingly critical. Our approach to these challenges uses a physics-based auroral model to regularize the poorly-observed vertical dimension. Specifically, the vertical dimension is expanded in a low-dimensional basis consisting of eigenprofiles computed over the range of expected energies in the precipitating electron flux, while the horizontal dimension retains a standard orthogonal pixel basis. Simulation results show typical characteristic energy estimation error less than 30% for a 3 km baseline achievable within the confines of the Poker Flat Research Range, using GPS-synchronized Electron Multiplying CCD cameras with broad-band BG3 optical filters that pass prompt auroral emissions.National Science Foundation Atmosphere and Geospace Directorate, Grants 1216530, 123737