Tamas I. Gombosi

Space Weather Modeling Framework: A new tool for the space science community

[1] The Space Weather Modeling Framework (SWMF) provides a high-performance flexible framework for physics-based space weather simulations, as well as for various space physics applications. The SWMF integrates numerical models of the Solar Corona, Eruptive Event Generator, Inner Heliosphere, Solar Energetic Particles, Global Magnetosphere, Inner Magnetosphere, Radiation Belt, Ionosphere Electrodynamics, and Upper Atmosphere into a high-performance coupled model. The components can be represented with alternative physics models, and any physically meaningful subset of the components can be used. The components are coupled to the control module via standardized interfaces, and an efficient parallel coupling toolkit is used for the pairwise coupling of the components. The execution and parallel layout of the components is controlled by the SWMF. Both sequential and concurrent execution models are supported. The SWMF enables simulations that were not possible with the individual physics models. Using reasonably high spatial and temporal resolutions in all of the coupled components, the SWMF runs significantly faster than real time on massively parallel supercomputers. This paper presents the design and implementation of the SWMF and some demonstrative tests. Future papers will describe validation (comparison of model results with measurements) and applications to challenging space weather events. The SWMF is publicly available to the scientific community for doing geophysical research. We also intend to expand the SWMF in collaboration with other model developers.

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.

Global three-dimensional MHD simulation of a space weather event: CME formation, interplanetary propagation, and interaction with the magnetosphere

A parallel adaptive mesh refinement (AMR) finite-volume scheme for predicting ideal MHD flows is used to simulate the initiation, structure, and evolution of a coronal mass ejection (CME) and its interaction with the magnetosphere-ionosphere system. The simulated CME is driven by a local plasma density enhancement on the solar surface with the background initial state of the corona and solar wind represented by a newly devised “steady state” solution. The initial solution has been constructed to provide a reasonable description of the time-averaged solar wind for conditions near solar minimum: (1) the computed magnetic field near the Sun possesses high-latitude polar coronal holes, closed magnetic field flux tubes at low latitudes, and a helmet streamer structure with a neutral line and current sheet; (2) the Archimedean spiral topology of the interplanetary magnetic field is reproduced; (3) the observed two-state nature of the solar wind is also reproduced with the simulation yielding fast and slow solar wind streams at high and low latitudes, respectively; and (4) the predicted solar wind plasma properties at 1 AU are consistent with observations. Starting with the generation of a CME at the Sun, the simulation follows the evolution of the solar wind disturbance as it evolves into a magnetic cloud and travels through interplanetary space and subsequently interacts with the terrestrial magnetosphere-ionosphere system. The density-driven CME exhibits a two-step release process, with the front of the CME rapidly accelerating following the disruption of the near-Sun closed magnetic field line structure and then moving at a nearly constant speed of ∼560 km/s through interplanetary space. The CME also produces a large magnetic cloud (> 100 RS across) characterized by a magnetic field that smoothly rotates northward and then back again over a period of ∼2 days at 1 AU. The cloud does not contain a sustained period with a strong southward component of the magnetic field, and, as a consequence, the simulated CME is somewhat ineffective in generating strong geo-magnetic activity at Earth. Nevertheless, the simulation results illustrate the potential, as well as current limitations, of the MHD-based space weather model for enhancing the understanding of coronal physics, solar wind plasma processes, magnetospheric physics, and space weather phenomena. Such models will provide the foundation for future, more comprehensive space weather prediction tools.

ERUPTION OF A BUOYANTLY EMERGING MAGNETIC FLUX ROPE

We present a three-dimensional numerical magnetohydrodynamic simulation designed to model the emergence of a magnetic flux rope passing from below the photosphere into the corona. For the initial state, we prescribe a plane-parallel atmosphere that comprises a polytropic convection zone, photosphere, transition region, and corona. Embedded in this system is an isolated horizontal magnetic flux rope located 10 photospheric pressure scale heights below the photosphere. The flux rope is uniformly twisted, with the plasma temperature inside the rope reduced to compensate for the magnetic pressure. Density is reduced in the middle of the rope, so that this section buoyantly rises. The early evolution proceeds with the middle of the rope rising to the photosphere and expanding into the corona. Just as it seems the system might approach equilibrium, the upper part of the flux rope begins to separate from the lower, mass-laden part. The separation occurs through stretching of the field, which forms a current sheet, where reconnection severs the field lines to form a new system of closed flux. This flux then erupts into the corona. Essential to the eruption process are shearing motions driven by the Lorentz force, which naturally occur as the rope expands in the pressure-stratified atmosphere. The shearing motions transport axial flux and energy to the expanding portion of the magnetic field, driving the eruption.

67P/Churyumov-Gerasimenko, a Jupiter family comet with a high D/H ratio

The provenance of water and organic compounds on Earth and other terrestrial planets has been discussed for a long time without reaching a consensus. One of the best means to distinguish between different scenarios is by determining the deuterium-to-hydrogen (D/H) ratios in the reservoirs for comets and Earth’s oceans. Here, we report the direct in situ measurement of the D/H ratio in the Jupiter family comet 67P/Churyumov-Gerasimenko by the ROSINA mass spectrometer aboard the European Space Agency’s Rosetta spacecraft, which is found to be (5.3 ± 0.7) × 10−4—that is, approximately three times the terrestrial value. Previous cometary measurements and our new finding suggest a wide range of D/H ratios in the water within Jupiter family objects and preclude the idea that this reservoir is solely composed of Earth ocean–like water.

Time variability and heterogeneity in the coma of 67P/Churyumov-Gerasimenko

Comets contain the best-preserved material from the beginning of our planetary system. Their nuclei and comae composition reveal clues about physical and chemical conditions during the early solar system when comets formed. ROSINA (Rosetta Orbiter Spectrometer for Ion and Neutral Analysis) onboard the Rosetta spacecraft has measured the coma composition of comet 67P/Churyumov-Gerasimenko with well-sampled time resolution per rotation. Measurements were made over many comet rotation periods and a wide range of latitudes. These measurements show large fluctuations in composition in a heterogeneous coma that has diurnal and possibly seasonal variations in the major outgassing species: water, carbon monoxide, and carbon dioxide. These results indicate a complex coma-nucleus relationship where seasonal variations may be driven by temperature differences just below the comet surface.

Three‐dimensional multispecies MHD studies of the solar wind interaction with Mars in the presence of crustal fields

interaction of the solar wind with Mars. The three ions considered are H + ,O 2 , and O + , representing the solar wind and the two major ionospheric ion species, respectively. The calculations indicate that the presence of a hot oxygen corona does not, within the resolution and accuracy of the model, lead to any significant effect on the dayside bow shock and ionopause positions. Next the trans-terminator fluxes and escape fluxes down the tail were calculated neglecting the effects of the crustal magnetic field. The calculated flux values are consistent with the measured escape fluxes and the calculated limiting fluxes from the dayside ionosphere. Finally, a 60-order harmonic expansion model of the measured magnetic field was incorporated into the model. The crustal magnetic field did not cause major distortions in the bow shock but certainly had an important effect within the magnetosheath and on the apparent altitude of the ionopause. The model results also indicated the presence of ‘‘minimagnetocylinders,’’ consistent with the MGS observations. We also recalculated the trans-terminator and escape fluxes, for the nominal solar wind case, in the presence of the crustal magnetic field and found, as expected, that there is a decrease in the calculated escape flux; however, it is still reasonably close to the value estimated from the Phobos-2 observations. INDEX TERMS: 2780 Magnetospheric Physics: Solar wind interactions with unmagnetized bodies; 2459 Ionosphere: Planetary ionospheres (5435, 5729, 6026, 6027, 6028); 5440 Planetology: Solid Surface Planets: Magnetic fields and magnetism; 2728 Magnetospheric Physics: Magnetosheath; KEYWORDS: Mars, MHD, bow shock, escape flux, solar wind interaction, crustal magnetic field Citation: Ma, Y., A. F. Nagy, K. C. Hansen, D. L. DeZeeuw, T. I. Gombosi, and K. G. Powell, Three-dimensional multispecies MHD studies of the solar wind interaction with Mars in the presence of crustal fields, J. Geophys. Res., 107(A10), 1282, doi:10.1029/2002JA009293, 2002.

Heliosphere in the magnetized local interstellar medium : Results of a three-dimensional MHD simulation

The results of a three-dimensional adaptive magnetohydrodynamic (MHD) model of the interaction of a magnetized solar wind with a magnetized very local interstellar medium in the presence of neutral interstellar hydrogen are presented. The interplanetary magnetic field is approximated by the Parker spiral, and the direction of the interstellar magnetic field is taken to be arbitrary. It is demonstrated that magnetic field interaction has a very pronounced effect on the structure of the global heliosphere. In particular, it is shown that the interaction of the interstellar wind with the shocked solar wind significantly depends on the direction of the interstellar magnetic field. This effect is mainly manifested in the distances to the heliospheric boundaries and the shape of the heliosphere. Depending on the orientation of the interstellar magnetic field the upstream location of the termination shock is expected to be at 80 ± 10 AU. The termination shock is predicted to be weak in agreement with the available body of observations. It is found that under quiet solar conditions the spiral structure of the interplanetery magnetic field remains imprinted in the solar wind far beyond the termination shock. Numerical simulations indicate that magnetic fields have a stabilizing effect on the heliopause.

Solar wind stagnation near comets

The nature of the solar wind flow near comets is examined analytically in this paper. In particular, typical values for the stagnation pressure and magnetic barrier strength are estimated, taking into account magnetic field line tension and change-exchange cooling of the mass-loaded solar wind. A knowledge of the strength of the magnetic barrier is required in order to determine the location of the ionopause surface which separates the contaminated solar wind plasma from the outflowing plasma of the cometary ionosphere.

Three-dimensional multiscale MHD model of cometary plasma environments

First results of a three-dimensional multiscale MHD model of the interaction of an expanding cometary atmosphere with the magnetized solar wind are presented. The model starts with a supersonic and super-Alfvenic solar wind far upstream of the comet (25 Gm upstream of the nucleus) with arbitrary interplanetary magnetic field orientation. The solar wind is continuously mass loaded with cometary ions originating from a 10-km size nucleus. The effects of photoionization, electron impact ionization, recombination, and ion-neutral frictional drag are taken into account in the model. The governing equations are solved on an adaptively refined unstructured Cartesian grid using our new multiscale upwind scalar conservation laws-type numerical technique (MUSCL). We have named this the multiscale adaptive upwind scheme for MHD (MAUS-MHD). The combination of the adaptive refinement with the MUSCL-scheme allows the entire cometary atmosphere to be modeled, while still resolving both the shock and the diamagnetic cavity of the comet. The main findings are the following: (1) Mass loading decelerates the solar wind flow upstream of the weak cometary shock wave (M ≈ 2, MA ≈ 2), which forms at a subsolar standoff distance of about 0.35 Gm. (2) A cometary plasma cavity is formed at around 3 × 103 km from the nucleus. Inside this cavity the plasma expands outward due to the frictional interaction between ions and neutrals. On the nightside this plasma cavity considerably narrows and a relatively fast and dense cometary plasma beam is ejected into the tail. (3) Inside the plasma cavity a teardrop-shaped inner shock is formed, which is terminated by a Mach disk on the nightside. Only the region inside the inner shock is the “true” diamagnetic cavity. (4) The model predicts four distinct current systems in the inner coma: the density peak current, the cavity boundary current, the inner shock current, and finally the cross-tail current. (5) The calculated plasma parameters (magnetic field, plasma density, speed, and temperature) are in very good agreement with published Giotto observations.