skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Thursday, February 12 until 1:00 AM ET on Friday, February 13 due to maintenance. We apologize for the inconvenience.


Title: Nonlinear Three-Dimensional Simulations of the Gradient Drift and Secondary Kelvin–Helmholtz Instabilities in Ionospheric Plasma Clouds
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
Award ID(s):
1847905 2048422
PAR ID:
10493462
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Atmosphere
Volume:
14
Issue:
4
ISSN:
2073-4433
Page Range / eLocation ID:
676
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The polar and high latitude regions of the ionosphere are host to complex plasma processes involving Magnetosphere-Ionosphere (MI) coupling, plasma convection, and auroral dynamics. The magnetic field lines from the polar cusp down through the auroral region map out to the magnetosphere and project the footprint of the large-scale convective processes driven by the solar wind onto the ionosphere. This region is also a unique environment where the magnetic field is oriented nearly vertical, resulting in horizontal drifts along closed, localized, convection patterns, and where prolonged periods of darkness during the winter result in the absence of significant photoionization. This set of conditions results in unique ionospheric structures which can set the stage for the generation of the gradient drift instability (GDI). The GDI occurs when the density gradient and ExB plasma drift are in the same direction. The GDI is a source of structuring at density gradients and may give rise to ionospheric irregularities that impact over-the-horizon radars and GPS signals. While the plasma ExB drifts are supplied by magnetospheric convection and MI coupling, sharp density gradients in the polar regions will be present at polar holes. Since the GDI occurs where the density gradient and plasma drift are parallel, the ionospheric irregularities caused by the GDI should occur at the leading edge of the polar hole. If so, the resulting production of small-scale density irregularities may, if the density is high enough, give rise to scintillation of GNSS signals and backscatter on HF radars. In this study, we investigate whether these irregularities can occur at the edges of polar holes as detected by the HF radar scatter. We use the Ionospheric Data Assimilation 4-Dimentional (IDA4D) and Assimilative Mapping of Ionospheric Electrodynamics (AMIE) models to characterize the high latitude ionospheric density and ExB drift convective structures, respectively, for one of nine polar hole events identified using RISR-N incoherent scatter radar in Forsythe et al [2021]. The combined IDA4D and AMIE assimilative outputs indicate where the GDI could be triggered, e.g., locations where the density gradient and ExB drift velocity have parallel components and the growth rate is smaller than the characteristic time over which the convective pattern changes, in this case, ~1/15 min. The presence of decameter ionospheric plasma irregularities is detected using the Super Dual Auroral Radar Network (SuperDARN). SuperDARN radars are HF coherent scatter radars. The presence of ionospheric radar returns in regions unstable to GDI grown strongly suggest the GDI is producing decameter scale plasma irregularities. The statistical analyses conducted in the above investigation do not show a clear pattern of enhanced scatter with larger computed GDI growth rates. Further investigation must be conducted before concluding that the GDI does not cause irregularities detectable with HF radar at polar holes. 
    more » « less
  2. We address the sources and dynamics of vorticity and helicity and their relations in transitions to turbulence arising due to Kelvin–Helmholtz instability (KHI) “Tube” and “Knot” (T&K) events. Such events are common in the atmosphere and oceans, and initial numerical simulations reveal that T&K dynamics significantly accelerate turbulence transitions and enhance KHI peak and mean energy dissipation rates. KHI T&K events arise where emerging KH billows exhibit varying wavelengths, phases, amplitudes, and/or discontinuities along their axes. As the KH billows intensify, these regions evolve roughly orthogonal billow cores and induced vortex tubes in close proximity. Their mutual advection as they intensify induces large-amplitude Kelvin vortex waves, or “twist waves”, that arise where locally uniform vortices are distorted by axial or radial advection. The twist waves propagate along, and fragment, the vortex tubes and billow cores, thus accounting for the emergence of helicity and the down-scale energy, enstrophy, and helicity fluxes within the turbulence inertial range. We describe the results of four direct numerical simulations (DNS) addressing KHI T&K dynamics in large and idealized small domains. The large-domain vorticity fields reveal the character and diversity of KHI T&K dynamics, the emergence of twist waves at larger and smaller scales, and their driving of turbulence transitions. Two small-domain DNS exhibit idealized KHI T&K events arising from KH billows that are mis-aligned and that exhibit phase variability along their axes. A third examines the interactions of two vortex tubes in close proximity. These reveal that twist waves drive the character and evolutions of the vorticity and helicity fields. 
    more » « less
  3. Abstract The occurrence of plasma irregularities and ionospheric scintillation over the Caribbean region have been reported in previous studies, but a better understanding of the source and conditions leading to these events is still needed. In December 2021, three ground-based ionospheric scintillation and Total Electron Content monitors were installed at different locations over Puerto Rico to better understand the occurrence of ionospheric irregularities in the region and to quantify their impact on transionospheric signals. Here, the findings for an event that occurred on March 13–14, 2022 are reported. The measurements made by the ground-based instrumentation indicated that ionospheric irregularities and scintillation originated at low latitudes and propagated, subsequently, to mid-latitudes. Imaging of the ionospheric F-region over a wide range of latitudes provided by the GOLD mission confirmed, unequivocally, that the observed irregularities and the scintillation were indeed caused by extreme equatorial plasma bubbles, that is, bubbles that reach abnormally high apex heights. The joint ground- and space-based observations show that plasma bubbles reached apex heights exceeding 2600 km and magnetic dip latitudes beyond 28 ° . In addition to the identification of extreme plasma bubbles as the source of the ionospheric perturbations over low-to-mid latitudes, GOLD observations also provided experimental evidence of the background ionospheric conditions leading to the abnormally high rise of the plasma bubbles and to severe L-band scintillation. These conditions are in good agreement with the theoretical hypothesis previously proposed. Graphical Abstract 
    more » « less
  4. Small-scale ionospheric plasma structures can cause scintillation in radio signals passing through the ionosphere. The relationship between the scintillated signal and how plasma structuring develops is complex. We model the development of small-scale plasma structuring in and around an idealized polar cap patch observed by the Resolute Bay Incoherent Scatter Radars (RISR) with the Geospace Environment Model for Ion-Neutral Interactions (GEMINI). Then, we simulate a signal passing through the resulting small-scale structuring with the Satellite-beacon Ionospheric-scintillation Global Model of the upper Atmosphere (SIGMA) to predict the scintillation characteristics that will be observed by a ground receiver at different stages of instability development. Finally, we compare the predicted signal characteristics with actual observations of scintillation from ground receivers in the vicinity of Resolute Bay. We interpret the results in terms of the nature of the small-scale plasma structuring in the ionosphere and how it impacts signals of different frequencies and attempt to infer information about the ionospheric plasma irregularity spectrum. 
    more » « less
  5. Small-scale ionospheric plasma structures can cause scintillation in radio signals passing through the ionosphere. The relationship between the scintillated signal and how plasma structuring develops is complex. We model the development of small-scale plasma structuring in and around an idealized polar cap patch observed by the Resolute Bay Incoherent Scatter Radars (RISR) with the Geospace Environment Model for Ion-Neutral Interactions (GEMINI). Then, we simulate a signal passing through the resulting small-scale structuring with the Satellite-beacon Ionospheric-scintillation Global Model of the upper Atmosphere (SIGMA) to predict the scintillation characteristics that will be observed by a ground receiver at different stages of instability development. Finally, we compare the predicted signal characteristics with actual observations of scintillation from ground receivers in the vicinity of Resolute Bay. We interpret the results in terms of the nature of the small-scale plasma structuring in the ionosphere and how it impacts signals of different frequencies and attempt to infer information about the ionospheric plasma irregularity spectrum. 
    more » « less