Alex Glocer

Adaptive numerical algorithms in space weather modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different relevant physics in different domains. A multi-physics system can be modeled by a software framework comprising several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamic (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics-based space weather modeling and even short-term forecasting.

Simulation of the 23 July 2012 Extreme Space Weather Event: What if This Extremely Rare CME Was Earth Directed?

Extreme space weather events are known to cause adverse impacts on critical modern day technological infrastructure such as high-voltage electric power transmission grids. On 23 July 2012, NASAs Solar Terrestrial Relations Observatory-Ahead (STEREO-A) spacecraft observed in situ an extremely fast coronal mass ejection (CME) that traveled 0.96 astronomical units (∼1 AU) in about 19 h. Here we use the Space Weather Modeling Framework (SWMF) to perform a simulation of this rare CME. We consider STEREO-A in situ observations to represent the upstream L1 solar wind boundary conditions. The goal of this study is to examine what would have happened if this Rare-type CME was Earth-bound. Global SWMF-generated ground geomagnetic field perturbations are used to compute the simulated induced geoelectric field at specific ground-based active INTERMAGNET magnetometer sites. Simulation results show that while modeled global SYM-H index, a high-resolution equivalent of the Dst index, was comparable to previously observed severe geomagnetic storms such as the Halloween 2003 storm, the 23 July CME would have produced some of the largest geomagnetically induced electric fields, making it very geoeffective. These results have important practical applications for risk management of electrical power grids.

The Comprehensive Inner Magnetosphere-Ionosphere Model

Simulation studies of the Earths 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.

Modeling Extreme 'Carrington-Type' Space Weather Events Using Three-dimensional Global MHD Simulations

There is a growing concern over possible severe societal consequences related to adverse space weather impacts on man-made technological infrastructure. In the last two decades, significant progress has been made toward the first-principles modeling of space weather events, and three-dimensional (3-D) global magnetohydrodynamics (MHD) models have been at the forefront of this transition, thereby playing a critical role in advancing our understanding of space weather. However, the modeling of extreme space weather events is still a major challenge even for the modern global MHD models. In this study, we introduce a specially adapted University of Michigan 3-D global MHD model for simulating extreme space weather events with a Dst footprint comparable to the Carrington superstorm of September 1859 based on the estimate by Tsurutani et. al. (2003). Results are presented for a simulation run with “very extreme” constructed/idealized solar wind boundary conditions driving the magnetosphere. In particular, we describe the reaction of the magnetosphere-ionosphere system and the associated induced geoelectric field on the ground to such extreme driving conditions. The model setup is further tested using input data for an observed space weather event of Halloween storm October 2003 to verify the MHD model consistence and to draw additional guidance for future work. This extreme space weather MHD model setup is designed specifically for practical application to the modeling of extreme geomagnetically induced electric fields, which can drive large currents in ground-based conductor systems such as power transmission grids. Therefore, our ultimate goal is to explore the level of geoelectric fields that can be induced from an assumed storm of the reported magnitude, i.e., Dst∼=−1600 nT.

Extended magnetohydrodynamics with embedded particle‐in‐cell simulation of Ganymede's magnetosphere

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHD-EPIC) algorithm is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the interaction between Jupiters magnetospheric plasma and Ganymedes magnetosphere. We compare the MHD-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the importance of kinetic effects in controlling the configuration and dynamics of Ganymedes magnetosphere. We find that the Hall MHD and MHD-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular, the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHD-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-D structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHD-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

Magnetosphere‐ionosphere energy interchange in the electron diffuse aurora

[1] The diffuse aurora has recently been shown to be a major contributor of energy flux into the Earth’s ionosphere. Therefore, a comprehensive theoretical analysis is required to understand its role in energy redistribution in the coupled ionosphere-magnetosphere system. In previous theoretical descriptions of precipitated magnetospheric electrons (E � 1 keV), the major focus has been the ionization and excitation rates of the neutral atmosphere and the energy deposition rate to thermal ionospheric electrons. However, these precipitating electrons will also produce secondary electrons via impact ionization of the neutral atmosphere. This paper presents the solution of the Boltzman-Landau kinetic equation that uniformly describes the entire electron distribution function in the diffuse aurora, including the affiliated production of secondary electrons (E < 600 eV) and their ionosphere-magnetosphere coupling processes. In this article, we discuss for the first time how diffuse electron precipitation into the atmosphere and the associated secondary electron production participate in ionosphere-magnetosphere energy redistribution.

How Hospitable Are Space Weather Affected Habitable Zones? The Role of Ion Escape

Atmospheres of exoplanets in the habitable zones around active young G-K-M stars are subject to extreme X-ray and EUV (XUV) fluxes from their host stars that can initiate atmospheric erosion. Atmospheric loss affects exoplanetary habitability in terms of surface water inventory, atmospheric pressure, the efficiency of greenhouse warming, and the dosage of the UV surface irradiation. Thermal escape models suggest that exoplanetary atmospheres around active K-M stars should undergo massive hydrogen escape, while heavier species including oxygen will accumulate forming an oxidizing atmosphere. Here, we show that non-thermal oxygen ion escape could be as important as thermal, hydrodynamic H escape in removing the constituents of water from exoplanetary atmospheres under supersolar XUV irradiation. Our models suggest that the atmospheres of a significant fraction of Earth-like exoplanets around M dwarfs and active K stars exposed to high XUV fluxes will incur a significant atmospheric loss rate of oxygen and nitrogen, which will make them uninhabitable within a few tens to hundreds of Myr, given a low replenishment rate from volcanism or cometary bombardment. Our non-thermal escape models have important implications for the habitability of the Proxima Centauris terrestrial planet.

The electric wind of Venus: A global and persistent “polar wind”‐like ambipolar electric field sufficient for the direct escape of heavy ionospheric ions

Understanding what processes govern atmospheric escape and the loss of planetary water is of paramount importance for understanding how life in the universe can exist. One mechanism thought to be important at all planets is an “ambipolar” electric field that helps ions overcome gravity. We report the discovery and first quantitative extraterrestrial measurements of such a field at the planet Venus. Unexpectedly, despite comparable gravity, we show the field to be five times stronger than in Earths similar ionosphere. Contrary to our understanding, Venus would still lose heavy ions (including oxygen and all water-group species) to space, even if there were no stripping by the solar wind. We therefore find that it is possible for planets to lose heavy ions to space entirely through electric forces in their ionospheres and such an “electric wind” must be considered when studying the evolution and potential habitability of any planet in any star system.

THE INTERACTION OF VENUS-LIKE, M-DWARF PLANETS WITH THE STELLAR WIND OF THEIR HOST STAR

We study the interaction between the atmospheres of Venus-like, non-magnetized exoplanets orbiting an M-dwarf star, and the stellar wind using a multi-species Magnetohydrodynaic (MHD) model. We focus our investigation on the effect of enhanced stellar wind and enhanced EUV flux as the planetary distance from the star decreases. Our simulations reveal different topologies of the planetary space environment for sub- and super-Alfvenic stellar wind conditions, which could lead to dynamic energy deposition in to the atmosphere during the transition along the planetary orbit. We find that the stellar wind penetration for non-magnetized planets is very deep, up to a few hundreds of kilometers. We estimate a lower limit for the atmospheric mass-loss rate and find that it is insignificant over the lifetime of the planet. However, we predict that when accounting for atmospheric ion acceleration, a significant amount of the planetary atmosphere could be eroded over the course of a billion years.

The role of the Hall effect in the global structure and dynamics of planetary magnetospheres: Ganymede as a case study

We present high-resolution Hall MHD simulations of Ganymedes magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moons surface, where they may contribute to Ganymedes aurora. Within the magnetopause and magnetotail current sheets, Hall J × B forces accelerate ions to the local Alfven speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymedes magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian “double magnetopause” structure, an Alfvenic ion jet extending across the upstream magnetopause, and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves.