skip to main content

Title: Simulations of Secondary Farley‐Buneman Instability Driven by a Kilometer‐Scale Primary Wave: Anomalous Transport and Formation of Flat‐Topped Electric Fields

Since the 1950s, high frequency and very high frequency radars near the magnetic equator have frequently detected strong echoes caused ultimately by the Farley‐Buneman instability (FBI) and the gradient drift instability (GDI). In the 1980s, coordinated rocket and radar campaigns made the astonishing observation of flat‐topped electric fields coincident with both meter‐scale irregularities and the passage of kilometer‐scale waves. The GDI in the daytimeEregion produces kilometer‐scale primary waves with polarization electric fields large enough to drive meter‐scale secondary FBI waves. The meter‐scale waves propagate nearly vertically along the large‐scale troughs and crests and act as VHF tracers for the large‐scale dynamics. This work presents a set of hybrid numerical simulations of secondary FBIs, driven by a primary kilometer‐scale GDI‐like wave. Meter‐scale density irregularities develop in the crest and trough of the kilometer‐scale wave, where the total electric field exceeds the FBI threshold, and propagate at an angle near the direction of total Hall drift determined by the combined electric fields. The meter‐scale irregularities transport plasma across the magnetic field, producing flat‐topped electric fields similar to those observed in rocket data and reducing the large‐scale wave electric field to just above the FBI threshold value. The self‐consistent reduction in driving electric field helps explain why echoes from the FBI propagate near the plasma acoustic speed.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Space Physics
Page Range / eLocation ID:
p. 734-748
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Thermal effects in the ionospheric plasma instabilities including the Farley‐Buneman instability (FBI), the gradient‐drift instability (GDI), and the ion‐thermal instability (ITI) are analyzed, focusing on analytic analysis of the previously obtained general expression of the combined instability growth rate. It is shown that thermal effects lead to generally nonmonotonic behavior of the growth rate with the background electric fieldor, equivalently with the plasma drift speed, in contrast with the conventional quadratic (FBI) or linear (GDI) dependence. The threshold electric field for FBI is demonstrated to exist at all altitudes, although at higher altitudes and/or at longer wavelength the threshold fields become too high to be observed in the ionosphere. The GDI growth rate is shown to be significantly modified by the thermal effects at long wavelengths. In the absence of thermal effects, the growth rate is proportional to a cosine of the flow angle, the angle between the differential plasma drift and the wavevector. Thermal effects result in an additional phase shift that maximizes at an altitude around 120 km, using representative high‐latitude ionospheric parameters, and in an increase in the proportionality coefficient by up to a factor of 2. The study highlights the usefulness of an analytic treatment for revealing additional insights into the instability behavior in the broad range of ionospheric parameters that may remain hidden using only numerical treatment.

    more » « less
  2. A newly developed three-dimensional electrostatic fluid model solving continuity and current closure equations aims to study phenomena that generate ionospheric turbulence. The model is spatially discretized using a pseudo-spectral method with full Fourier basis functions and evolved in time using a four-stage, fourth-order Runge Kutta method. The 3D numerical model is used here to investigate the behavior and evolution of ionospheric plasma clouds. This problem has historically been used to study the processes governing the evolution of the irregularities in the F region of the ionosphere. It has been shown that these artificial clouds can become unstable and structure rapidly (i.e., cascade to smaller scales transverse to the ambient magnetic field). The primary mechanism which causes this structuring of ionospheric clouds is the E×B, or the gradient drift instability (GDI). The persistence and scale sizes of the resulting structures cannot be fully explained by a two-dimensional model. Therefore, we suggest here that the inclusion of three-dimensional effects is key to a successful interpretation of mid-latitude irregularities, as well as a prerequisite for a credible simulation of these processes. We investigate the results of 2D and 3D nonlinear simulations of the GDI and secondary Kelvin–Helmholtz instability (KHI) in plasma clouds for three different regimes: highly collisional (≈200 km), collisional (≈300 km), and inertial (≈450 km). The inclusion of inertial effects permits the growth of the secondary KHI. For the three different regimes, the overall evolution of structuring of plasma cloud occurs on longer timescales in 3D simulations. The inclusion of three-dimensional effects, in particular, the ambipolar potential in the current closure equation, introduces an azimuthal “twist“ about the axis of the cloud (i.e., the magnetic field B). This azimuthal “twist” is observed in the purely collisional regime, and it causes the perturbations to have a non-flute-like character (k‖≠0). However, for the 3D inertial simulations, the cloud rapidly diffuses to a state in which the sheared azimuthal flow is substantially reduced; subsequently, the cloud becomes unstable and structures, by retaining the flute-like character of the perturbations (k‖=0).

    more » « less
  3. Abstract

    Observations of 30‐MHz coherent backscatter from sporadic‐Eionization layers were obtained with a VHF imaging radar located in Ithaca, New York. The volume probed by the radar lies at relatively high magnetic latitudes, on the northern edge of the mid‐latitude region and underneath the ionospheric trough. Banded, quasi‐periodic (QP) echoes observed from Ithaca are similar to those found in lower midlatitude regions. The Doppler shifts observed are smaller and, so far, do not appear to reach the threshold for Farley‐Buneman instability. However, many of the echoes exhibit fine‐scale structure, with secondary bands or braids oriented obliquely to the primary bands. Secondary bands have been seen only rarely at lower middle latitudes. In previous observations, the QP scattering has been linked to unstable neutral wind shears. Neutral wind shear commonly found in the lower thermosphere could play a key role in the formation of these irregularities and explain some morphological features of the resulting plasma density irregularities and the radar echoes. We consider whether neutral instability and turbulence in the lower thermosphere is the likely cause for some of the structuring in the sporadic‐Elayers. Results of 3D numerical simulations of atmospheric dynamics in the mesosphere to lower thermosphere support the proposition. In particular, we focus on Ekman‐type instabilities that, like the more common Kelvin‐Helmholtz instabilities, are inflection point instabilities, although specifically associated with turning shears, and result in convective rolls aligned close to the mean wind direction, with smaller‐scale secondary waves aligned normal to the primary structures.

    more » « less
  4. Abstract

    Behavior of unstable plasma waves generated by the Farley‐Buneman instability (FBI) and the gradient drift instability (GDI) is analyzed in the transitional valley region near 120 km in altitude. The analysis is based on the expression for the FBI/GDI growth rateγthat has been recently generalized to include ion inertia effects for arbitrary altitude and wavelength, within the limits imposed by the fluid and local approaches. It is found that the ion inertia leads to a different instability behavior when the convection component is between the two critical values determined by the ion acoustic speedCsand the ratioribetween the ion collision and gyrofrequency. The most interesting case occurs near 120 km, just below whereri=1. From analysis of electron density gradientsG=n/nthat result in marginal instability conditionγ=0 (i.e., critical gradientsG0), there exists a critical scale whereG0=0 and below which all waves are unstable to FBI. Above this scale,G0>0 and gradients need to be sufficiently strongG>G0for the plasma to become unstable through GDI. There also exists a maximum in dependence, which refers to the least unstable scale and gradient. For convection outside of the specified range, no critical or least unstable scale exists, which is a typical situation outside of the transitional valley region. Overall, this analysis shows that the FBI convection thresholds and the GDI critical gradients are modified by the ion inertia and that the effects are most pronounced in the transitional valley region near 120 km.

    more » « less
  5. Abstract

    Radar and sounding rocket observations of plasma irregularities in theF‐region ionosphere acquired on 19 June 2019 during NASA experiment Too WINDY on Kwajalein Atoll are presented. The experiment was conducted near local midnight during a period of low solar flux and quiet geomagnetic conditions. Plasma density irregularities were seen by the rocket and also in the incoherent scatter radar data to emerge and persist mainly in the topside. Density irregularities in the bottomside remained very small by comparison throughout the observations. Zonal plasma drifts measured by the rocket were highly structured in the topside. Patches of coherent scatter entrained in the large‐scale topside density irregularities appeared to propagate slowly westward in a narrow flow channel detected by the rocket. Broadband ELF emissions were also detected in the topside. Some of the characteristics of the topside irregularities are typical of postsunset equatorialF‐region irregularities observed frequently by coherent scatter radars, and some of the common features in the coherent scatter database are reviewed. Two scenarios that have been proposed to account for postmidnight spreadFare tested computationally. One involves unseasonably large background zonal electric fields, and the other involves forcing from below by neutral waves and turbulence. Neither scenario appears to be able to account for the Too WINDY observations and the preponderance of topside irregularities without bottomside precursors in particular.

    more » « less