Since the 1960s, the Farley-Buneman instability has played an important role in probing the E-region ionosphere. The intervening years have seen significant progress in the linear theory of this instability, its relation to other instabilities, and some of its observational signatures. However, the saturation mechanism and nonlinear behavior remain open topics because of their role in controlling energy flow in the E-region plasma. This paper explores the saturated state of the Farley-Buneman instability in 2-D and 3-D kinetic simulations of the high-latitude ionosphere, at three different simulated altitudes: 107, 110, and 113 km. These simulations show irregularity amplitude growth and saturation in all runs, but irregularity growth takes much longer in 2-D than in 3-D. Once the simulations reach saturation, wave power in the meter-scale regime falls off as a power law below the wavelength of peak growth, but the power law index is larger in 3-D than in 2-D. At longer wavelengths, the 3-D spectrum is much flatter than the 2-D spectrum. This implies that purely 2-D simulations of the Farley-Buneman instability may overestimate irregularity amplitudes at decameter scales and may also underestimate the efficiency of ion Landau damping at the ion mean-free-path scale. From a physical perspective, the relatively flat spectra above the wavelength of peak growth in 3-D simulations imply a wavelength-independent saturation mechanism across a range of altitudes. Finally, both 2-D and 3-D simulations demonstrate the importance of accounting for zeroth-order ion drift when estimating the flow angle of density irregularities.
Plain Language SummaryAt approximately 90-120 km in altitude, the Earth's atmosphere comprises a gas containing mostly neutral molecules, and a small number of positively charged ions and electrons. This sort of gas is called a weakly ionized plasma. Aeronomy researchers have known for decades that the weakly ionized plasma around 90-120 km in Earth's atmosphere can grow unstable and reflect radio waves, and they have used radio-wave reflection as one technique to probe this region of geospace. However, the aeronomy community still does not understand the nature of the fully developed plasma instabilities in this region. This work addresses that problem with computer simulations. It finds that 2-D simulations do not capture all the characteristics of 3-D simulations, which should mimic reality more closely. It also finds that whatever mechanism creates the fully developed 3-D instability should apply to wavelengths from a few meters to tens of meters. Finally, this work suggests that the direction along which the unstable waves travel should change in a predictable way as altitude increases. These conclusions will help future researchers estimate atmospheric parameters based on radar and rocket measurements.