Search NASA⌕ Search

Engineering topics

Raeder, Joachim

Publications and source records attributed to Raeder, Joachim.

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.↗

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.↗