skip to main content

This content will become publicly available on December 1, 2023

Title: An electrodynamics model for Data Interpretation and Numerical Analysis of ionospheric Missions and Observations (DINAMO)
Abstract We introduce a new numerical model developed to assist with Data Interpretation and Numerical Analysis of ionospheric Missions and Observations (DINAMO). DINAMO derives the ionospheric electrostatic potential at low- and mid-latitudes from a two-dimensional dynamo equation and user-specified inputs for the state of the ionosphere and thermosphere (I–T) system. The potential is used to specify the electric fields and associated F -region E × B plasma drifts. Most of the model was written in Python to facilitate the setup of numerical experiments and to engage students in numerical modeling applied to space sciences. Here, we illustrate applications and results of DINAMO in two different analyses. First, DINAMO is used to assess the ability of widely used I–T climatological models (IRI-2016, NRLMSISE-00, and HWM14), when used as drivers, to produce a realistic representation of the low-latitude electrodynamics. In order to evaluate the results, model E × B drifts are compared with observed climatology of the drifts derived from long-term observations made by the Jicamarca incoherent scatter radar. We found that the climatological I–T models are able to drive many of the features of the plasma drifts including the diurnal, seasonal, altitudinal and solar cycle variability. We also identified discrepancies between more » modeled and observed drifts under certain conditions. This is, in particular, the case of vertical equatorial plasma drifts during low solar flux conditions, which were attributed to a poor specification of the E -region neutral wind dynamo. DINAMO is then used to quantify the impact of meridional currents on the morphology of F -region zonal plasma drifts. Analytic representations of the equatorial drifts are commonly used to interpret observations. These representations, however, commonly ignore contributions from meridional currents. Using DINAMO we show that that these currents can modify zonal plasma drifts by up to ~ 16 m/s in the bottom-side post-sunset F -region, and up to ~ 10 m/s between 0700 and 1000 LT for altitudes above 500 km. Finally, DINAMO results show the relationship between the pre-reversal enhancement (PRE) of the vertical drifts and the vertical shear in the zonal plasma drifts with implications for equatorial spread F. « less
Award ID(s):
Publication Date:
Journal Name:
Progress in Earth and Planetary Science
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. In March 2014 an all-sky imager (ASI) was installed at the Jicamarca Radio Observatory (11.95°S, 76.87°W; 0.3°S MLAT). We present results of equatorial spread F (ESF) characteristics observed at Jicamarca and at low latitudes. Optical 6300 and 7774Å airglow observations from the Jicamarca ASI are compared with other collocated instruments and with ASIs at El Leoncito, Argentina (31.8°S, 69.3°W; 19.8°S MLAT), and Villa de Leyva, Colombia (5.6°N, 73.52°W; 16.4°N MLAT). We use Jicamarca radar data, in incoherent and coherent modes, to obtain plasma parameters and detect echoes from irregularities. We find that ESF depletions tend to appear in groups with a group-to-group separation around 400–500km and within-group separation around 50–100km. We combine data from the three ASIs to investigate the conditions at Jicamarca that could lead to the development of high-altitude, or topside, plumes. We compare zonal winds, obtained from a Fabry–Pérot interferometer, with plasma drifts inferred from the zonal motion of plasma depletions. In addition to the ESF studies we also investigate the midnight temperature maximum and its effects at higher latitudes, visible as a brightness wave at El Leoncito. The ASI at Jicamarca along with collocated andmore »low-latitude instruments provide a clear two-dimensional view of spatial and temporal evolution of ionospheric phenomena at equatorial and low latitudes that helps to explain the dynamics and evolution of equatorial ionospheric/thermospheric processes.

    « less
  2. Variations of vertical atmospheric electric field E z have been attributed mainly to meteorological processes. On the other hand, the theory of electromagnetic waves in the atmosphere, between the bottom ionosphere and earth’s surface, predicts two modes, magnetic H (TE) and electric E (TH) modes, where the E-mode has a vertical electric field component, E z . Past attempts to find signatures of ULF (periods from fractions to tens of minutes) disturbances in E z gave contradictory results. Recently, study of ULF disturbances of atmospheric electric field became feasible thanks to project GLOCAEM, which united stations with 1 sec measurements of potential gradient. These data enable us to address the long-standing problem of the coupling between atmospheric electricity and space weather disturbances at ULF time scales. Also, we have reexamined results of earlier balloon-born electric field and ground magnetic field measurements in Antarctica. Transmission of storm sudden commencement (SSC) impulses to lower latitudes was often interpreted as excitation of the electric TH 0 mode, instantly propagating along the ionosphere–ground waveguide. According to this theoretical estimate, even a weak magnetic signature of the E-mode ∼1 nT must be accompanied by a burst of E z well exceeding the atmospheric potential gradient. We havemore »examined simultaneous records of magnetometers and electric field-mills during >50 SSC events in 2007–2019 in search for signatures of E-mode. However, the observed E z disturbance never exceeded background fluctuations ∼10 V/m, much less than expected for the TH 0 mode. We constructed a model of the electromagnetic ULF response to an oscillating magnetospheric field-aligned current incident onto the realistic ionosphere and atmosphere. The model is based on numerical solution of the full-wave equations in the atmospheric-ionospheric collisional plasma, using parameters that were reconstructed using the IRI model. We have calculated the vertical and horizontal distributions of magnetic and electric fields of both H- and E-modes excited by magnetospheric field-aligned currents. The model predicts that the excitation rate of the E-mode by magnetospheric disturbances is low, so only a weak E z response with a magnitude of ∼several V/m will be produced by ∼100 nT geomagnetic disturbance. However, at balloon heights (∼30 km), electric field of the E-mode becomes dominating. Predicted amplitudes of horizontal electric field in the atmosphere induced by Pc5 pulsations and travelling convection vortices, about tens of mV/m, are in good agreement with balloon electric field and ground magnetometer observations.« less
  3. 630.0 nm all-sky imaging data are used to detect airglow depletions associated with equatorial spread F. Pairs of imagers located at geomagnetically conjugate locations in the American sector at low and mid-latitudes provide information on the occurrence rate and zonal motion of airglow depletions. Airglow depletions are seen extending to magnetic latitudes as high as 25°. An asymmetric extension is observed with structures in the northern hemisphere reaching higher latitudes. By tracking the zonal motion of airglow depletions, zonal plasma drifts in the thermosphere can be inferred and their simultaneous behavior in both hemispheres investigated. Case studies using El Leoncito and Mercedes imagers in the southern hemisphere, and the respective magnetically conjugate imagers at Villa de Leyva and Arecibo, provide consistent evidence of the influence of the South Atlantic Magnetic Anomaly on the dynamics and characteristics of the thermosphere–ionosphere system at low and mid-latitudes.
  4. Context. Coronal mass ejections (CMEs) are large eruptions of magnetised plasma from the Sun that are often accompanied by solar radio bursts produced by accelerated electrons. Aims. A powerful source for accelerating electron beams are CME-driven shocks, however, there are other mechanisms capable of accelerating electrons during a CME eruption. So far, studies have relied on the traditional classification of solar radio bursts into five groups (Type I–V) based mainly on their shapes and characteristics in dynamic spectra. Here, we aim to determine the origin of moving radio bursts associated with a CME that do not fit into the present classification of the solar radio emission. Methods. By using radio imaging from the Nançay Radioheliograph, combined with observations from the Solar Dynamics Observatory, Solar and Heliospheric Observatory, and Solar Terrestrial Relations Observatory spacecraft, we investigate the moving radio bursts accompanying two subsequent CMEs on 22 May 2013. We use three-dimensional reconstructions of the two associated CME eruptions to show the possible origin of the observed radio emission. Results. We identified three moving radio bursts at unusually high altitudes in the corona that are located at the northern CME flank and move outwards synchronously with the CME. The radio bursts correspondmore »to fine-structured emission in dynamic spectra with durations of ∼1 s, and they may show forward or reverse frequency drifts. Since the CME expands closely following an earlier CME, a low coronal CME–CME interaction is likely responsible for the observed radio emission. Conclusions. For the first time, we report the existence of new types of short duration bursts, which are signatures of electron beams accelerated at the CME flank. Two subsequent CMEs originating from the same region and propagating in similar directions provide a complex configuration of the ambient magnetic field and favourable conditions for the creation of collapsing magnetic traps. These traps are formed if a CME-driven wave, such as a shock wave, is likely to intersect surrounding magnetic field lines twice. Electrons will thus be further accelerated at the mirror points created at these intersections and eventually escape to produce bursts of plasma emission with forward and reverse drifts.« less
  5. Context. Spectroscopic observations of the emission lines formed in the solar transition region commonly show persistent downflows on the order of 10−15 km s −1 . The cause of such downflows, however, is still not fully clear and has remained a matter of debate. Aims. We aim to understand the cause of such downflows by studying the coronal and transition region responses to the recently reported chromospheric downflowing rapid redshifted excursions (RREs) and their impact on the heating of the solar atmosphere. Methods. We have used two sets of coordinated data from the Swedish 1 m Solar Telescope, the Interface Region Imaging Spectrograph, and the Solar Dynamics Observatory for analyzing the response of the downflowing RREs in the transition region and corona. To provide theoretical support, we use an already existing 2.5D magnetohydrodynamic simulation of spicules performed with the Bifrost code. Results. We find ample occurrences of downflowing RREs and show several examples of their spatio-temporal evolution, sampling multiple wavelength channels ranging from the cooler chromospheric to the hotter coronal channels. These downflowing features are thought to be likely associated with the returning components of the previously heated spicular plasma. Furthermore, the transition region Doppler shifts associated with them aremore »close to the average redshifts observed in this region, which further implies that these flows could (partly) be responsible for the persistent downflows observed in the transition region. We also propose two mechanisms – (i) a typical upflow followed by a downflow and (ii) downflows along a loop –from the perspective of a numerical simulation that could explain the ubiquitous occurrence of such downflows. A detailed comparison between the synthetic and observed spectral characteristics reveals a distinctive match and further suggests an impact on the heating of the solar atmosphere. Conclusions. We present evidence that suggests that at least some of the downflowing RREs are the chromospheric counterparts of the transition region and lower coronal downflows.« less