Engineering Papers⌕ Search

Engineering topics

Raeder, Joachim

Publications and source records attributed to Raeder, Joachim.

Explosive Magnetotail Activity

Modes and manifestations of the explosive activity in the Earth’s magnetotail, as well as its onset mechanisms and key pre-onset conditions are reviewed. Two mechanisms for the generation of the pre-onset current sheet are discussed, namely magnetic flux addition to the tail lobes, or other high-latitude perturbations, and magnetic flux evacuation from the near-Earth tail associated with dayside reconnection. Reconnection onset may require stretching and thinning of the sheet down to electron scales. It may also start in thicker sheets in regions with a tailward gradient of the equatorial magnetic field B$_z$; in this case it begins as an ideal-MHD instability followed by the generation of bursty bulk flows and dipolarization fronts. Indeed, remote sensing and global MHD modeling show the formation of tail regions with increased B$_z$, prone to magnetic reconnection, ballooning/interchange and flapping instabilities. While interchange instability may also develop in such thicker sheets, it may grow more slowly compared to tearing and cause secondary reconnection locally in the dawn-dusk direction. Post-onset transients include bursty flows and dipolarization fronts, micro-instabilities of lower-hybrid-drift and whistler waves, as well as damped global flux tube oscillations in the near-Earth region. They convert the stretched tail magnetic field energy into bulk plasma acceleration and collisionless heating, excitation of a broad spectrum of plasma waves, and collisional dissipation in the ionosphere. Collisionless heating involves ion reflection from fronts, Fermi, betatron as well as other, non-adiabatic, mechanisms. Ionospheric manifestations of some of these magnetotail phenomena are discussed. Explosive plasma phenomena observed in the laboratory, the solar corona and solar wind are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron Physics in 3‐D Two‐Fluid 10‐Moment Modeling of Ganymede's Magnetosphere

Abstract We studied the role of electron physics in 3‐D two‐fluid 10‐moment simulation of Ganymede's magnetosphere. The model captures nonideal physics like the Hall effect, electron inertia, and anisotropic, nongyrotropic pressure effects. A series of analyses were carried out: (1) The resulting magnetic field topology and electron and ion convection patterns were investigated. The magnetic fields were shown to agree reasonably well with in situ measurements by the Galileo satellite. (2) The physics of collisionless magnetic reconnection were carefully examined in terms of the current sheet formation and decomposition of generalized Ohm's law. The importance of pressure anisotropy and nongyrotropy in supporting the reconnection electric field is confirmed. (3) We compared surface “brightness” morphology, represented by surface electron and ion pressure contours, with oxygen emission observed by the Hubble Space Telescope. The correlation between the observed emission morphology and spatial variability in electron/ion pressure was demonstrated. Potential extension to multi‐ion species in the context of Ganymede and other magnetospheric systems is also discussed.

Wang, Liang↗

Simulation of magnetic holes formation in the magnetosheath

Magnetic holes have been frequently observed in the Earth's magnetosheath and are believed to be the consequence of the nonlinear evolution of the mirror instability. Mirror mode perturbations mainly form as magnetic holes in regions where the plasma is marginally mirror stable with respect to the linear instability criterion. We present an expanding box particle-in-cell simulation to mimic the changing conditions in the magnetosheath as the plasma is convected through it that produces mirror mode magnetic holes. We show that in the initial nonlinear evolution, where the plasma conditions are mirror unstable, the magnetic peaks are dominant, while later, as the plasma relaxes toward marginal stability, the fluctuations evolve into deep magnetic holes. Finally, while the averaged plasma parameters in the simulation remain close to the mirror instability threshold, the local plasma in the magnetic holes is highly unstable to mirror instability and locally mirror stable in the magnetic peaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reply to comment by Remya et al. on “Effects of electron temperature anisotropy on proton mirror instability evolution”

In the comment to our paper, Remya et al. (2017) state that we conclude that their theory is incorrect; however, no such conclusion is in our paper. In fact, as stated in their paper, we agree with their theory that shows the impact of heavy ions and electron temperature anisotropy on the competition of the ion anisotropy instabilities. While their linear theory is correct, our paper focused on the nonlinear evolution, where one needs to be careful in assuming a given electron anisotropy, because electrons themselves can be unstable to the electron whistler instability, which quickly lowers the anisotropy to levels where, in the absence of heavy ions, it is not sufficient to significantly change the balance between proton cyclotron and mirror mode. We agree that the electron whistler instability will not lead to complete isotropization of the electrons but only lower it to the instability threshold. In the parameter regime addressed, this limited isotropization will still eliminate the dominance of mirror mode and restore the usual dominance of the proton cyclotron mode, so our point still stands. Our simulations showed an isotropization of the electrons beyond the electron whistler threshold. In this paper, we will show that there are two contributing reasons: The nonlinear evolution of the mirror instability affects the electron anisotropy, as does unphysical numerical heating due to the limited resolution of a particle-in-cell simulation. Finally, we further discuss the coexistence of electron whistler instability and mirror instability, and we agree that both instabilities can be present in the magnetosheath.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the Ionosphere-Thermosphere Response to a Geomagnetic Storm Using Physics-based Magnetospheric Energy Input: OpenGGCM-CTIM Results

The magnetosphere is a major source of energy for the Earth's ionosphere and thermosphere (IT) system. Current IT models drive the upper atmosphere using empirically calculated magnetospheric energy input. Thus, they do not sufficiently capture the storm-time dynamics, particularly at high latitudes. To improve the prediction capability of IT models, a physics-based magnetospheric input is necessary. Here, we use the Open Global General Circulation Model (OpenGGCM) coupled with the Coupled Thermosphere Ionosphere Model (CTIM). OpenGGCM calculates a three-dimensional global magnetosphere and a two-dimensional high-latitude ionosphere by solving resistive magnetohydrodynamic (MHD) equations with solar wind input. CTIM calculates a global thermosphere and a high-latitude ionosphere in three dimensions using realistic magnetospheric inputs from the OpenGGCM. We investigate whether the coupled model improves the storm-time IT responses by simulating a geomagnetic storm that is preceded by a strong solar wind pressure front on August 24, 2005. We compare the OpenGGCM-CTIM results with low-earth-orbit satellite observations and with the model results of Coupled Thermosphere-Ionosphere-Plasmasphere electrodynamics (CTIPe). CTIPe is an up-to-date version of CTIM that incorporates more IT dynamics such as a low-latitude ionosphere and a plasmasphere, but uses empirical magnetospheric input. OpenGGCMCTIM reproduces localized neutral density peaks at approx. 400 km altitude in the high-latitude dayside regions in agreement with in situ observations during the pressure shock and the early phase of the storm. Although CTIPe is in some sense a much superior model than CTIM, it misses these localized enhancements. Unlike the CTIPe empirical input models, OpenGGCM-CTIM more faithfully produces localized increases of both auroral precipitation and ionospheric electric fields near the high-latitude dayside region after the pressure shock and after the storm onset, which in turn effectively heats the thermosphere and causes the neutral density increase at 400 km altitude.

Connor, Hyunju K.↗

Effects of electron temperature anisotropy on proton mirror instability evolution

Proton mirror modes are large amplitude nonpropagating structures frequently observed in the magnetosheath. It has been suggested that electron temperature anisotropy can enhance the proton mirror instability growth rate while leaving the proton cyclotron instability largely unaffected, therefore causing the proton mirror instability to dominate the proton cyclotron instability in Earth's magnetosheath. In this paper we use particle-in-cell simulations to investigate the electron temperature anisotropy effects on proton mirror instability evolution. Contrary to the hypothesis, electron temperature anisotropy leads to excitement of the electron whistler instability. Finally, our results show that the electron whistler instability grows much faster than the proton mirror instability and quickly consumes the electron-free energy so that there is no electron temperature anisotropy left to significantly impact the evolution of the proton mirror instability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

GEM-CEDAR Challenge: Poynting Flux at DMSP and Modeled Joule Heat

Poynting flux into the ionosphere measures the electromagnetic energy coming from the magnetosphere. This energy flux can vary greatly between quiet times and geomagnetic active times. As part of the Geospace Environment Modeling-coupling energetics and dynamics of atmospheric regions modeling challenge, physics-based models of the 3-D ionosphere and ionospheric electrodynamics solvers of magnetosphere models that specify Joule heat and empirical models specifying Poynting flux were run for six geomagnetic storm events of varying intensity. We compared model results with Poynting flux values along the DMSP-15 satellite track computed from ion drift meter and magnetic field observations. Although being a different quantity, Joule heat can in practice be correlated to incoming Poynting flux because the energy is dissipated primarily in high latitudes where Poynting flux is being deposited. Within the physics-based model group, we find mixed results with some models overestimating Joule heat and some models agreeing better with observed Poynting flux rates as integrated over auroral passes. In contrast, empirical models tend to underestimate integrated Poynting flux values. Modeled Joule heat or Poynting flux patterns often resemble the observed Poynting flux patterns on a large scale, but amplitudes can differ by a factor of 2 or larger due to the highly localized nature of observed Poynting flux deposition that is not captured by the models. In addition, the positioning of modeled patterns appear to be randomly shifted against the observed Poynting flux energy input. This study is the first to compare Poynting flux and Joule heat in a large variety of models of the ionosphere.

Rastaetter, Lutz↗

GEM-CEDAR Challenge: Comparing Ionospheric Models with Poynting Flux from DMSP Observations

As part to the GEM-CEDAR challenge we are extending the model-data comparisons to electrodynamic in-situ measurements in low-Earth orbit. We use DMSP observations of electric and magnetic fields to compute Poynting Flux values along the satellite track in high latitudes including the auroral zones and the polar cap. Models of the ionosphere that include electrodynamic parameters have been run for five events selected for the GEM-CEDAR modeling challenge for which DMSP data are available for comparison. Combined with a magnetic field model we use the modeled electric fields to compute Poynting Flux and Joule Dissipation values from outputs of CTIPe, TIE-GCM, the ionospheric electrodynamics solvers of the SWMF, LFM and OpenGGCM magnetosphere-ionosphere coupled models, and the Weimer electric field model. The online metrics analysis tool at the Community Coordinated Modeling Center (CCMC) has been updated to handle the analysis of separate short segments of available data (high-latitude sections of the satellite orbit) with model outputs to analyze how well auroral patterns are being reproduced by the models. We present initial results from the new analysis tool in terms of model yields (ratio of the difference between maximum and minimum values of model results to the observation), timing/location errors of local maxima in the inbound and outbound auroral crossings as well as cross-correlations for individual passes. We collect the information for many DMSP passes and present an analysis for model performance during quiet and geomagnetically disturbed time periods using half-orbit integrated values as well.

Rastaetter, Lutz↗

Dst Index in the 2008 GEM Modeling Challenge - Model Performance for Moderate and Strong Magnetic Storms

The GEM 2008 modeling challenge efforts are expanding beyond comparing in-situ measurements in the magnetosphere and ionosphere to include the computation of indices to be compared. The Dst index measures the largest deviations of the horizontal magnetic field at 4 equatorial magnetometers from the quiet-time background field and is commonly used to track the strength of the magnetic disturbance of the magnetosphere during storms. Models can calculate a proxy Dst index in various ways, including using the Dessler-Parker Sckopke relation and the energy of the ring current and Biot-Savart integration of electric currents in the magnetosphere. The GEM modeling challenge investigates 4 space weather events and we compare models available at CCMC against each other and the observed values of Ost. Models used include SWMF/BATSRUS, OpenGGCM, LFM, GUMICS (3D magnetosphere MHD models), Fok-RC, CRCM, RAM-SCB (kinetic drift models of the ring current), WINDMI (magnetosphere-ionosphere electric circuit model), and predictions based on an impulse response function (IRF) model and analytic coupling functions with inputs of solar wind data. In addition to the analysis of model-observation comparisons we look at the way Dst is computed in global magnetosphere models. The default value of Dst computed by the SWMF model is for Bz the Earth's center. In addition to this, we present results obtained at different locations on the Earth's surface. We choose equatorial locations at local noon, dusk (18:00 hours), midnight and dawn (6:00 hours). The different virtual observatory locations reveal the variation around the earth-centered Dst value resulting from the distribution of electric currents in the magnetosphere during different phases of a storm.

Rastaetter, Lutz↗

Origins and Transport of Ions during Magnetospheric Substorms

We investigate the origins and the transport of ions observed in the near-Earth plasma sheet during the growth and expansion phases of a magnetospheric substorm that occurred on November 24, 1996. Ions observed at Geotail were traced backward in time in time-dependent magnetic and electric fields to determine their origins and the acceleration mechanisms responsible for their energization. Results from this investigation indicate that, during the growth phase of the substorm, most of the ions reaching Geotail had origins in the low latitude boundary layer (LLBL) and had alread@, entered the magnetosphere when the growth phase began. Late in the growth phase and in the expansion phase a higher proportion of the ions reaching Geotail had their origin in the plasma mantle. Indeed, during the expansion phase more than 90% of the ions seen by Geotail were from the mantle. The ions were accelerated enroute to the spacecraft; however, most of the ions' energy gain was achieved by non-adiabatic acceleration while crossing the equatorial current sheet just prior to their detection by Geotail. In general, the plasma mantle from both southern and northern hemispheres supplied non-adiabatic ions to Geotail, whereas the LLBL supplied mostly adiabatic ions to the distributions measured by the spacecraft.

Ashour-Abdalla, Maha↗

Using Global Simulations of the Magnetosphere for Multi-Satellite Mission Planning and Data Analysis

We use global simulations of Earth's magnetosphere to assess the scientific return from a multi-satellite mission in the magnetosphere. We examine 4 different scenarios with 20, 40, 80, and 160 satellites, respectively. The satellite orbits are randomized with perigee distances ranging from 2 to 5R(sub E), apogee distances between 10 and 50 R(sub E), and within +/-5R(sub E) of the geocentric solar ecliptic (GSE) equator. For each of these satellite configurations we examine the expected observations during a typical substorm by using time traces obtained from a global simulation at the satellite positions. The 160 satellite configuration yields sufficient information to distinguish between different substorm models without any temporal/spatial ambiguities. An 80 satellite configuration still provides sufficient information for this task, however for fewer events with good satellite conjunctions and with less statistical certainty. For constellations with fewer than 40 satellites time-space ambiguities are likely to remain in the observation. However, any multi-satellite constellation would be a quantum leap in magnetospheric research because of the unprecedented coverage of other regions, because it would enable new measurement techniques that are unique to multi-satellite missions, and because it would enable the use of data assimilation techniques in global models for the first time.

Raeder, Joachim↗

Development of a Model for the Night Side Magnetopause using Global Simulations

During the final year of this investigation we have finished several event studies that we considered necessary for the development of a tail magnetopause model and for the calibration of our simulation code. We have not reached the ultimate goal of the project, i.e., the development of an analytical tail magnetopause model. In the course of the investigation we have learned that such a model would be much more complex than we had anticipated. However, the investigations that we conducted towards this goal have led to significant results and discoveries that are of considerable value for understanding the tail magnetopause. These are summarized in the following sections.

Raeder, Joachim↗

Visualization of numerical simulation results: An interactive tool for mission oriented theory

With the advent of multispacecraft missions, it became obvious that theoretical models are needed, not only to interpret local spacecraft observations, but also to link single point measurements. Models that address space plasma phenomena on local and global scales, that involve nonlinear processes, and that use self-consistent approaches, were developed. An interactive system for visualizing and analyzing numerical simulation results was also developed. This system is based on visualization modules which use the application visualization system and the NCAR graphic packages. Examples from recent simulations illustrate how these modules can be used for displaying and manipulating simulation results.

Berchem, Jean↗

Magnetic flux ropes at the high-latitude magnetopause

We examine the consequences of magnetic reconnection at the high-latitude magnetopause using a three-dimensional global magnetohydrodynamic simulation of the solar wind interaction with the Earth's magnetosphere. Magnetic field lines from the simulation reveal the formation of magnetic flux ropes during periods with northward interplanetary magnetic field. These flux ropes result from multiple reconnection processes between the lobes field lines and draped magnetosheath field lines that are convected around the flank of the magnetosphere. The flux ropes identified in the simulation are consistent with features observed in the magnetic field measured by Hawkeye-1 during some high-latitude magnetopause crossings.

Berchem, Jean↗

A simulation study of multiple ion wave generation downstream of low Mach number quasiperpendicular shocks

The behavior of minor ions just downstream of a low Mach number quasi-perpendicular shock is investigated both theoretically and by computer simulations. Because all ions see the same cross shock electric field their deceleration depends on their charge to mass ratio, yielding different downstream velocities. It is shown that these differences in velocity can lead to coherent wave structures in the downstream region of quasi-perpendicular shocks with a narrow transition layer. These waves are shown to be multi ion hybrid waves in contrast to mirror waves and ion cyclotron waves. Under favorable conditions these waves should be observable both at interplanetary shocks and at planetary bowshocks.

Motschmann, Uwe↗

A global magnetohydrodynamic simulation of magnetospheric dynamics when the IMF is southward: Mapping to the auroral zone

A high resolution global magnetohydrodynamic simulation model is used to investigate magnetospheric dynamics during intervals with southward interplanetary magnetic field (IMF). When the southward IMF reaches the dayside magnetopause reconnection begins and magnetic flux is convected into the tail lobes. After about 35 m, reconnection begins within the plasma sheet near midnight at x = -14 RE. Later the x-line moves towards the magnetopause. The reconnection occurs just tailward of the region where the tail attaches onto the dipole dominated inner magnetosphere. Later when all the plasma sheet field lines have reconnected a plasmoid moves down the tail. The region of the ionosphere where the energy flux from the magnetosphere is greatest is calculated. The energy flux is confined to a region which approximates the auroral oval.

Walker, Raymond J.↗