The reversal in the direction of the ionospheric electric field in the region known as the Harang discontinuity is known to play an important role in the large-scale dynamics of the ionosphere and magnetosphere. This project uses the Rice Convection Model (RCM) to simulate the Harang reversal in order to evaluate the physical processes responsible for its formation and characteristics, as well as its evolution during substorms. The RCM can appropriately take into account the physics of the plasma drift transport within the plasma sheet, and it solves for the ionospheric electric potential self-consistently with the plasma sheet dynamics through their coupling by field-aligned currents.
This project will further develop simulations to be more realistic. A magnetic field solver will provide magnetic fields that are in force balance with the RCM plasma pressures. Particle boundary conditions will be established using observations from the Geotail satellite. This combined model will be used to provide quantitative evaluation of the physical processes believed to govern the Harang reversal. The quantitative results from the simulations will be compared with ground-based radar observations. To evaluate the underlying physical processes, individual physical parameters in the model will be varied to see how the variations effect the Harang reversal. By varying a single paramter while keeping others fixed it will be possible to determine the effect of each parameter on the formation and characteristics of the Harang reversal.
To understand the relationship between the Harang reversal and magnetic substorms and to test a substorm theory, this project will evaluate the evolution of the Harang reversal and its correlations to the plasma sheet during the substorm growth phase and during the initial stage after a strong decrease in plasma sheet convection. The project will examine whether or not a decrease in the convective electric field leads towards increasingly rapid changes and instability as predicted by the theory.
The simulations will be the first to provide quantitative assessment of the Harang reversal, its correlations with the plasma sheet, its evolution and response to variations of interplanetary conditions during substorms, and its potential use to predict substorms, which will be crucial to understanding magnetosphere-ionosphere coupling within the plasma sheet and substorm triggering.
The Harang reversal is a region in the nightside ionosphere at auroral latitudes where the north-south component of electric fields reverses direction. Observations have shown that the Harang reversal is highly related to space weather disturbances. Its formation is a result of coupling of physical processes between the ionosphere and magnetosphere (M-I coupling). It enhances during the substorm growth phase and decreases dramatically at onset and it location during the substorm growth phase controls the location of a future onset. In this project, we used Rice Convection Model (RCM) with a force-balanced magnetic field solver to understand the underlying M-I coupling physics of the Harang reversal and to evaluate how the characteristics of the Harang reversal changes in response to the solar wind-magnetosphere interaction that leads to the substorm growth phase and geomagnetic storms. From the simulations, we have determined that the existence of an overlap in local time of Region 2 upward and downward field-aligned current (FAC) is necessary for the formation of the Harang reversal. The distribution of FAC is mainly determined by particles of different energy with the downward FAC in this overlap region associated with electrons and low energy ions and the upward FAC associated with high-energy ions. We have found as the plasma sheet density and temperature change with the interplanetary magnetic field (IMF), a colder and denser plasma sheet during northward IMF results in a shaper Harang reversal and at higher latitudes than does a hotter and less dense plasma sheet during southward IMF. There are often mesoscale plasma sheet disturbances during substorms and storms. Our results revealed that the Harang reversal responses quite differently to different types of mesoscale disturbances. A bubble (lower density) disturbance enhances the Harang reversal while a blob (higher density) disturbance weakens the reversal. During the storm main phase, the Harang reversal expands wider in magnetic local times (longitudes) and extends to lower latitudes than during the substorm growth phase. The Harang reversal disappears following a sharp decrease in the coupling between solar wind-and magnetosphere, such as that happens when IMF turns from southward to northward. These important findings have greatly advanced our understanding of "coupling between the solar wind and the magnetosphere and between the magnetosphere and ionosphere", provided explanations and predictions to the observed changes in the Harang reversal, and provided improvement to the predictions of current space weather models. The above scientific achievements were obtained through the tasks of developing numerical model and simulations, as well as analyzing satellite observations and simulation results. The majority of these tasks were conducted by a female Ph.D. student, through which she has learned and acquired the essential skills of science research. She has received her Ph.D. degree in 2012 and is now an independent researcher with many collaborations with other researchers and continues to make scientific contributions to the space physics community.