Engineering topics

Fok, M.-C.

Publications and source records attributed to Fok, M.-C..

At least 19 records

NASA NTRS2018

Imaging Plasma Density Structures in the Soft X-Rays Generated by Solar Wind Charge Exchange with Neutrals

Both heliophysics and planetary physics seek to understand the complex nature of the solar wind's interaction with solar system obstacles like Earth's magnetosphere, the ionospheres of Venus and Mars, and comets. Studies with this objective are frequently conducted with the help of single or multipoint in situ electromagnetic field and particle observations, guided by the predictions of both local and global numerical simulations, and placed in context by observations from far and extreme ultraviolet (FUV, EUV), hard X-ray, and energetic neutral atom imagers (ENA). Each proposed interaction mechanism (e.g., steady or transient magnetic reconnection, local or global magnetic reconnection, ion pick-up, or the Kelvin- Helmholtz instability) generates diagnostic plasma density structures. The significance of each mechanism to the overall interaction (as measured in terms of atmospheric/ionospheric loss at comets, Venus, and Mars or global magnetospheric/ionospheric convection at Earth) remains to be determined but can be evaluated on the basis of how often the density signatures that it generates are observed as a function of solar wind conditions. This paper reviews efforts to image the diagnostic plasma density structures in the soft (low energy, 0.1-2.0 keV) X-rays produced when high charge state solar wind ions exchange electrons with the exospheric neutrals surrounding solar system obstacles.

Sibeck, David G.↗

NASA NTRS2018

The Unknown Hydrogen Exosphere: Space Weather Implications

Recent studies suggest that the hydrogen (H) density in the exosphere and geocorona might differ from previously assumed values by factors as large as 2. We use the SAMI3 (Sami3 is Also a Model of the Ionosphere) and Comprehensive Inner Magnetosphere-Ionosphere models to evaluate scenarios where the hydrogen density is reduced or enhanced, by a factor of 2, relative to values given by commonly used empirical models. We show that the rate of plasmasphere refilling following a geomagnetic storm varies nearly linearly with the hydrogen density. We also show that the ring current associated with a geomagnetic storm decays more rapidly when H is increased. With respect to these two space weather effects, increased exosphere hydrogen density is associated with reduced threats to space assets during and following a geomagnetic storm.

Krall, J.↗

NASA NTRS2016

Inverse Energy Dispersion of Energetic Ions Observed in the Magnetosheath

We present a case study of energetic ions observed by the Energetic Particle Detector (EPD) on the Magnetospheric Multiscale spacecraft in the magnetosheath just outside the subsolar magnetopause that occurred at 1000 UT on 8 December 2015. As the magnetopause receded inward, the EPD observed a burst of energetic (approximately 50-1000 keV) proton, helium, and oxygen ions that exhibited an inverse dispersion, with the lowest energy ions appearing first. The prolonged interval of fast antisunward flow observed in the magnetosheath and transient increases in the H components of global ground magnetograms demonstrate that the burst appeared at a time when the magnetosphere was rapidly compressed. We attribute the inverse energy dispersion to the leakage along reconnected magnetic field lines of betatron-accelerated energetic ions in the magnetosheath, and a burst of reconnection has an extent of about 1.5 R(sub E) using combined Super Dual Auroral Radar Network radar and EPD observations.

Lee, S. H.↗

NASA NTRS2016

Simulation of a Rapid Dropout Event for Highly Relativistic Electrons with the RBE Model

A flux dropout is a sudden and sizable decrease in the energetic electron population of the outer radiation belt on the time scale of a few hours. We simulated a flux dropout of highly relativistic 2.5 MeV electrons using the Radiation Belt Environment model, incorporating the pitch angle diffusion coefficients caused by electromagnetic ion cyclotron (EMIC) waves for the geomagnetic storm events of 23-26 October 2002. This simulation showed a remarkable decrease in the 2.5 MeV electron flux during main phase of the storm, compared to those without EMIC waves. This decrease was independent of magnetopause shadowing or drift loss to the magnetopause. We suggest that the flux decrease was likely to be primarily due to pitch angle scattering to the loss cone by EMIC waves. Furthermore, the 2.5 MeV electron flux calculated with EMIC waves correspond very well with that observed from Solar Anomalous and Magnetospheric Particle EXplorer spacecraft. EMIC wave scattering is therefore likely one of the key mechanisms to understand flux dropouts. We modeled EMIC wave intensities by the Kp index. However, the calculated dropout is a several hours earlier than the observed one. We propose that Kp is not the best parameter to predict EMIC waves.

Kang, S-B.↗

NASA NTRS2014

The Comprehensive Inner Magnetosphere-Ionosphere Model

Simulation studies of the Earth's radiation belts and ring current are very useful in understanding the acceleration, transport, and loss of energetic particles. Recently, the Comprehensive Ring Current Model (CRCM) and the Radiation Belt Environment (RBE) model were merged to form a Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. CIMI solves for many essential quantities in the inner magnetosphere, including ion and electron distributions in the ring current and radiation belts, plasmaspheric density, Region 2 currents, convection potential, and precipitation in the ionosphere. It incorporates whistler mode chorus and hiss wave diffusion of energetic electrons in energy, pitch angle, and cross terms. CIMI thus represents a comprehensive model that considers the effects of the ring current and plasmasphere on the radiation belts. We have performed a CIMI simulation for the storm on 5-9 April 2010 and then compared our results with data from the Two Wide-angle Imaging Neutral-atom Spectrometers and Akebono satellites. We identify the dominant energization and loss processes for the ring current and radiation belts. We find that the interactions with the whistler mode chorus waves are the main cause of the flux increase of MeV electrons during the recovery phase of this particular storm. When a self-consistent electric field from the CRCM is used, the enhancement of MeV electrons is higher than when an empirical convection model is applied. We also demonstrate how CIMI can be a powerful tool for analyzing and interpreting data from the new Van Allen Probes mission.

Inner Magnetosphere↗

NASA NTRS2011

Rapid Rebuilding of the Outer Radiation Belt

Recent observations by the radiation monitor (RDM) on the spacecraft Akebono have shown several cases of greater than 2.5 MeV radiation belt electron enhancements occurring on timescales of less than a few hours. Similar enhancements are also seen in detectors on board the NOAA/POES and TWINS 1 satellites. These intervals are shorter than typical radial diffusion or wave-particle interactions can account for. We choose two so-called "rapid rebuilding" events that occur during high speed streams (4 September 2008 and 22 July 2009) and simulated them with the Space Weather Modeling Framework configured with global magnetosphere, radiation belt, ring current, and ionosphere electrodynamics model. Our simulations produce a weaker and delayed dipolarization as compared to observations, but the associated inductive electric field in the simulations is still strong enough to rapidly transport and accelerate MeV electrons resulting in an energetic electron flux enhancement that is somewhat weaker than is observed. Nevertheless, the calculated flux enhancement and dipolarization is found to be qualitatively consistent with the observations. Taken together, the modeling results and observations support the conclusion that storm-time dipolarization events in the magnetospheric magnetic field result in strong radial transport and energization of radiation belt electrons.

Glocer, A.↗

NASA NTRS2011

Simulation and Twins Observations of the 22 July 2009 Storm

TWINS is the first mission to perform stereo imaging of the Earth's ring current. The magnetic storm on 22 July 2009 was at the time the largest storm observed since TWINS began routine stereo imaging in June 2008. On 22 July 2009, the Dst dropped to nearly .80 nT at 0700 and 1000 UT. During the main phase, and at the peak of the storm, TWINS 1 and 2 were near apogee and moving between predawn and postdawn local time. The energetic neutral atom (ENA) imagers on the two spacecraft captured the storm intensification and the formation of the partial ring current. The peak of the high-altitude ENA emissions was seen in the midnight-to-dawn local time sector. The development of this storm has been simulated using the comprehensive ring current model (CRCM) to understand and interpret the observed signatures. We perform CRCM runs with constant and time-varying magnetic field. The model calculations are validated by comparing the simulated ENA and ion flux intensities with TWINS ENA images and in situ ion data from a THEMIS satellite. Simulation with a static magnetic field produces a strong shielding electric field that skews the ion drift trajectories toward dawn. The model's corresponding peak ENA emissions are always more eastward than those in the observed TWINS images. On the other hand, the simulation with a dynamic magnetic field gives better spatial agreement with both ENA and in situ particle data, suggesting that temporal variations of the geomagnetic field exert a significant influence upon global ring current ion dynamics.

Fok, M.-C.↗

NASA NTRS2010

Magnetic Flux Circulation During Dawn-Dusk Oriented Interplanetary Magnetic Field

Magnetic flux circulation is a primary mode of energy transfer from the solar wind into the ionosphere and inner magnetosphere. For southward interplanetary magnetic field (IMF), magnetic flux circulation is described by the Dungey cycle (dayside merging, night side reconnection, and magnetospheric convection), and both the ionosphere and inner magnetosphere receive energy. For dawn-dusk oriented IMF, magnetic flux circulation is not well understood, and the inner magnetosphere does not receive energy. Several models have been suggested for possible reconnection patterns; the general pattern is: dayside merging; reconnection on the dayside or along the dawn/dusk regions; and, return flow on dayside only. These models are consistent with the lack of energy in the inner magnetosphere. We will present evidence that the Dungey cycle does not explain the energy transfer during dawn-dusk oriented IMF. We will also present evidence of how magnetic flux does circulate during dawn-dusk oriented IMF, specifically how the magnetic flux reconnects and circulates back.

Mitchell, E. J.↗

NASA NTRS2010

Ring Current Dynamics in Moderate and Strong Storms: Comparative Analysis of TWINS and IMAGE/HENA Data with the Comprehensive Ring Current Model

We present a comparative study of ring current dynamics during strong and moderate storms. The ring current during the strong storm is studied with IMAGE/HENA data near the solar cycle maximum in 2000. The ring current during the moderate storm is studied using energetic neutral atom (ENA) data from the Two Wide-Angle Imaging Neutral- Atom Spectrometers (TWINS) mission during the solar minimum in 2008. For both storms, the local time distributions of ENA emissions show signatures of postmidnight enhancement (PME) during the main phases. To model the ring current and ENA emissions, we use the Comprehensive Ring Current Model (CRCM). CRCM results show that the main-phase ring current pressure peaks in the premidnight-dusk sector, while the most intense CRCM-simulated ENA emissions show PME signatures. We analyze two factors to explain this difference: the dependence of charge-exchange cross section on energy and pitch angle distributions of ring current. We find that the IMF By effect (twisting of the convection pattern due to By) is not needed to form the PME. Additionally, the PME is more pronounced for the strong storm, although relative shielding and hence electric field skewing is well developed for both events.

Buzulukova, N.↗

NASA NTRS2010

Global Response to Local Ionospheric Mass Ejection

We revisit a reported "Ionospheric Mass Ejection" using prior event observations to guide a global simulation of local ionospheric outflows, global magnetospheric circulation, and plasma sheet pressurization, and comparing our results with the observed global response. Our simulation framework is based on test particle motions in the Lyon-Fedder-Mobarry (LFM) global circulation model electromagnetic fields. The inner magnetosphere is simulated with the Comprehensive Ring Current Model (CRCM) of Fok and Wolf, driven by the transpolar potential developed by the LFM magnetosphere, and includes an embedded plasmaspheric simulation. Global circulation is stimulated using the observed solar wind conditions for the period 24-25 Sept 1998. This period begins with the arrival of a Coronal Mass Ejection, initially with northward, but later with southward interplanetary magnetic field. Test particles are launched from the ionosphere with fluxes specified by local empirical relationships of outflow to electrodynamic and particle precipitation imposed by the MIlD simulation. Particles are tracked until they are lost from the system downstream or into the atmosphere, using the full equations of motion. Results are compared with the observed ring current and a simulation of polar and auroral wind outflows driven globally by solar wind dynamic pressure. We find good quantitative agreement with the observed ring current, and reasonable qualitative agreement with earlier simulation results, suggesting that the solar wind driven global simulation generates realistic energy dissipation in the ionosphere and that the Strangeway relations provide a realistic local outflow description.

Moore, T. E.↗

NASA NTRS2010

Recent Developments in the Radiation Belt Environment Model

The fluxes of energetic particles in the radiation belts are found to be strongly controlled by the solar wind conditions. In order to understand and predict the radiation particle intensities, we have developed a physics-based Radiation Belt Environment (RBE) model that considers the influences from the solar wind, ring current and plasmasphere. Recently, an improved calculation of wave-particle interactions has been incorporated. In particular, the model now includes cross diffusion in energy and pitch-angle. We find that the exclusion of cross diffusion could cause significant overestimation of electron flux enhancement during storm recovery. The RBE model is also connected to MHD fields so that the response of the radiation belts to fast variations in the global magnetosphere can be studied.Weare able to reproduce the rapid flux increase during a substorm dipolarization on 4 September 2008. The timing is much shorter than the time scale of wave associated acceleration.

Fok, M.-C.↗

NASA NTRS2010

Dynamics of Ring Current and Electric Fields in the Inner Magnetosphere During Disturbed Periods: CRCM-BATS-R-US Coupled Model

We present simulation results from a one-way coupled global MHD model (Block-Adaptive-Tree Solar-Wind Roe-Type Upwind Scheme, BATS-R-US) and kinetic ring current models (Comprehensive Ring Current Model, CRCM, and Fok Ring Current, FokRC). The BATS-R-US provides the CRCM/FokRC with magnetic field information and plasma density/temperature at the polar CRCM/FokRC boundary. The CRCM uses an electric potential from the BATS-R-US ionospheric solver at the polar CRCM boundary in order to calculate the electric field pattern consistent with the CRCM pressure distribution. The FokRC electric field potential is taken from BATS-R-US ionospheric solver everywhere in the modeled region, and the effect of Region II currents is neglected. We show that for an idealized case with southward-northward-southward Bz IMF turning, CRCM-BATS-R-US reproduces well known features of inner magnetosphere electrodynamics: strong/weak convection under the southward/northward Bz; electric field shielding/overshielding/penetration effects; an injection during the substorm development; Subauroral Ion Drift or Polarization Jet (SAID/PJ) signature in the dusk sector. Furthermore, we find for the idealized case that SAID/PJ forms during the substorm growth phase, and that substorm injection has its own structure of field-aligned currents which resembles a substorm current wedge. For an actual event (12 August 2000 storm), we calculate ENA emissions and compare with Imager for Magnetopause-to-Aurora Global Exploration/High Energy Neutral Atom data. The CRCM-BATS-R-US reproduces both the global morphology of ring current and the fine structure of ring current injection. The FokRC-BATS-R-US shows the effect of a realistic description of Region II currents in ring current-MHD coupled models.

Buzulukova, N.↗

NASA NTRS2009

Self-Consistent Model of Magnetospheric Electric Field, Ring Current, Plasmasphere, and Electromagnetic Ion Cyclotron Waves: Initial Results

Further development of our self-consistent model of interacting ring current (RC) ions and electromagnetic ion cyclotron (EMIC) waves is presented. This model incorporates large scale magnetosphere-ionosphere coupling and treats self-consistently not only EMIC waves and RC ions, but also the magnetospheric electric field, RC, and plasmasphere. Initial simulations indicate that the region beyond geostationary orbit should be included in the simulation of the magnetosphere-ionosphere coupling. Additionally, a self-consistent description, based on first principles, of the ionospheric conductance is required. These initial simulations further show that in order to model the EMIC wave distribution and wave spectral properties accurately, the plasmasphere should also be simulated self-consistently, since its fine structure requires as much care as that of the RC. Finally, an effect of the finite time needed to reestablish a new potential pattern throughout the ionosphere and to communicate between the ionosphere and the equatorial magnetosphere cannot be ignored.

Gamayunov, K. V.↗

NASA NTRS2007

Self-consistent Model of Magnetospheric Electric Field, RC and EMIC Waves

Electromagnetic ion cyclotron (EMIC) waves are an important magnetospheric emission, which is excited near the magnetic equator with frequencies below the proton gyro-frequency. The source of bee energy for wave growth is provided by temperature anisotropy of ring current (RC) ions, which develops naturally during inward convection from the plasma sheet These waves strongly affect the dynamic s of resonant RC ions, thermal electrons and ions, and the outer radiation belt relativistic electrons, leading to non-adiabatic particle heating and/or pitch-angle scattering and loss to the atmosphere. The rate of ion and electron scattering/heating is strongly controlled by the Wave power spectral and spatial distributions, but unfortunately, the currently available observational information regarding EMIC wave power spectral density is poor. So combinations of reliable data and theoretical models should be utilized in order to obtain the power spectral density of EMIC waves over the entire magnetosphere throughout the different storm phases. In this study, we present the simulation results, which are based on two coupled RC models that our group has developed. The first model deals with the large-scale magnetosphere-ionosphere electrodynamic coupling, and provides a self-consistent description of RC ions/electrons and the magnetospheric electric field. The second model is based on a coupled system of two kinetic equations, one equation describes the RC ion dynamics and another equation describes the power spectral density evolution of EMIC waves, and self-consistently treats a micro-scale electrodynamic coupling of RC and EMIC waves. So far, these two models have been applied independently. However, the large-scale magnetosphere-ionosphere electrodynamics controls the convective patterns of both the RC ions and plasmasphere altering conditions for EMIC wave-particle interaction. In turn, the wave induced RC precipitation Changes the local field-aligned current distributions and the ionospheric conductances, which are crucial for a large-scale electrodynamics. The initial results from this new self-consistent model of the magnetospheric electric field, RC and EMIC waves will be shown in this presentation.

Gamayunov, K. V.↗

NASA NTRS2007

Energetic Neutral Atom Emissions From Venus: VEX Observations and Theoretical Modeling

Venus has almost no intrinsic magnetic field to shield itself from its surrounding environment. The solar wind thus directly interacts with the planetary ionosphere and atmosphere. One of the by-products of this close encounter is the production of energetic neutral atom (ENA) emissions. Theoretical studies have shown that significant amount of ENAs are emanated from the planet. The launch of the Venus Express (VEX) in 2005 provided the first light ever of the Venus ENA emissions. The observed ENA flux level and structure are in pretty good agreement with the theoretical studies. In this paper, we present VEX ENA data and the comparison with numerical simulations. We seek to understand the solar wind interaction with the planet and the impacts on its atmospheres.

Fok, M.-C.↗

NASA NTRS2006

Global Aspects of Heliosphere-Geosphere Coupling

The magnetosphere serves both to hold off the solar wind and to couple it selectively to the Earth through the auroral zones. It has long been understood that the plasmasphere consists of geogenic plasmas that expand out of the sunlit low latitude ionosphere, but it was initially assumed that the hot low density plasma beyond the plasmasphere is largely of solar wind origin with a minor admixture of ionospheric plasmas, discovered via mass spectrometric observations in the early 70's. Since then, elaborate simulation models have been developed that have taken us from an era of cartoon physics to a new era of quantitative global comparisons between observations and theory. In most current global circulation models of the magnetosphere, the ionospheric load on the system is taken to lie exclusively in the thin F layer of the ionosphere. This layer is coupled with solar wind and magnetospheric plasmas via Maxwell stresses communicated by field aligned current systems, and with the thermosphere via ion-neutral charge exchange and Coulomb collisions. However, recent observations have shown us that ionospheric plasmas flow sporadically in various forms into the plasmasphere and up into the high latitude circulation cells of the outer magnetosphere, as driven by the solar wind and its variable intensity and magnetic field. Under some conditions, ionospheric material is observed to be the principal component of plasmas at the dayside magnetopause and in the plasma sheet and ring current regions. Given a global model of magnetospheric circulation, it is relatively straightforward to investigate the behaviors of ionospheric plasmas in response to solar wind drivers, and we report the results of such efforts here. We find that ionospheric plasmas dominate the plasma pressure in the magnetosphere in some regions, particularly when the solar wind is especially intense and-or southward directed. This result violates the assumption that the ionospheric load is confined to the F layer, and shows that the ionosphere is often an important dynamic element of the solar wind-ionosphere system throughout the magnetosphere, especially during larger geospace storms. This means that future global circulation models must account for ionospheric plasma inertia, heat capacity, and pressure, to be quantitatively credible.

Moore, Thomas E.↗

NASA NTRS2002

Self-Consistent Magnetosphere-Ionosphere Coupling: Theoretical Studies

A self-consistent ring current (RC) model has been developed that couples electron and ion magnetospheric dynamics with the calculation of the electric field. Two new features were taken into account in order to close the self-consistent magnetosphere- ionosphere coupling loop. First, in addition to the RC ions, we solve an electron kinetic equation in our model. Second, using the relation of Galand and Richmond [2001], we calculate the height integrated ionospheric conductances as a function of the precipitated high energy magnetospheric electrons and ions that are produced by our code. To validate the results of our model we simulate the magnetic storm of May 2, 1986, a storm that has has been comprehensively studied, and compare our results with different theoretical approaches. The self-consistent inclusion of the hot electrons and, their effect on the conductance results in deeper penetration of the magnetospheric electric field. In addition, a slight westward rotation of the potential pattern (compared to previous self-consistent results) is evident in the inner magnetosphere. This changes the hot plasma distribution, especially by allowing increased access of plasma sheet ions and electrons to low L shells.

Khazanov, G. V.↗

NASA NTRS1998

IMAGE Mission Science

Simultaneous, global measurement of major magnetospheric plasma systems will be performed for the first time with the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) Mission. The ring current, plasmasphere, and auroral systems will be imaged using energetic neutral and ultraviolet cameras. Quantitative remote measurement of the magnetosheath, plasmaspheric, and magnetospheric densities will be obtained through radio sounding by the Radio Plasma Imager. The IMAGE Mission will open a new era in global magnetospheric physics, while bringing with it new challenges in data analysis. An overview of the IMAGE Theory and Modeling team efforts will be presented, including the state of development of Internet tools that will be available to the science community for access and analysis of IMAGE observations.

Gallagher, D. L.↗