We describe, analyze, and validate the recently developed Alfven Wave Solar Model, a three-dimensional global model starting from the top of the chromosphere and extending into interplanetary space ...(out to 1-2 AU). This model solves the extended, two-temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low-frequency Alfven waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and to wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counterpropagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment in any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model; no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady state (in the corotating frame) properties of the solar environment are reproduced, using three adjustable parameters: the Poynting flux of chromospheric Alfven waves, the perpendicular correlation length of the turbulence, and a pseudoreflection coefficient. We compare model results for Carrington rotation 2063 (2007 November-December) with remote observations in the extreme-ultraviolet and X-ray ranges from the Solar Terrestrial Relations Observatory, Solar and Heliospheric Observatory, and Hinode spacecraft and with in situ measurements by Ulysses. The results are in good agreement with observations. This is the first global simulation that is simultaneously consistent with observations of both the thermal structure of the lower corona and the wind structure beyond Earth's orbit.
We present a new MHD model for simulating the large-scale structure of the solar corona and solar wind under "steady state" conditions stemming from the Wang-Sheeley-Arge empirical model. The ...processes of turbulent heating in the solar wind are parameterized using a phenomenological, thermodynamical model with a varied polytropic index. We employ the Bernoulli integral to bridge the asymptotic solar wind speed with the assumed distribution of the polytropic index on the solar surface. We successfully reproduce the mass flux from Sun to Earth, the temperature structure, and the large-scale structure of the magnetic field. We reproduce the solar wind speed bimodal structure in the inner heliosphere. However, the solar wind speed is in a quantitative agreement with observations at 1 AU for solar maximum conditions only. The magnetic field comparison demonstrates that the input magnetogram needs to be multiplied by a scaling factor in order to obtain the correct magnitude at 1 AU.
Context. Pre-equinox measurements of comet 67P/Churyumov-Gerasimenko with the mass spectrometer ROSINA/DFMS on board the Rosetta spacecraft revealed a strongly heterogeneous coma. The abundances of ...major and various minor volatile species were found to depend on the latitude and longitude of the nadir point of the spacecraft. The observed time variability of coma species remained consistent for about three months up to equinox. The chemical variability could be generally interpreted in terms of surface temperature and seasonal effects superposed on some kind of chemical heterogeneity of the nucleus. Aims. We compare here pre-equinox (inbound) ROSINA/DFMS measurements from 2014 to measurements taken after the outbound equinox in 2016, both at heliocentric distances larger than 3 AU. For a direct comparison we limit our observations to the southern hemisphere. Methods. We report the similarities and differences in the concentrations and time variability of neutral species under similar insolation conditions (heliocentric distance and season) pre- and post-equinox, and interpret them in light of the previously published observations. In addition, we extend both the pre- and post-equinox analysis by comparing species concentrations with a mixture of CO2 and H2O. Results. Our results show significant changes in the abundances of neutral species in the coma from pre- to post-equinox that are indicative of seasonally driven nucleus heterogeneity. Conclusions. The observed pre- and post-equinox patterns can generally be explained by the strong erosion in the southern hemisphere that moves volatile-rich layers near the surface.
The recent solar minimum with very low activity provides us a unique opportunity for validating solar wind models. During CR2077, the number of sunspots was near the absolute minimum of solar cycle ...23. For this solar rotation, we perform a multi-spacecraft validation study for the recently developed three-dimensional, two-temperature, Alfv'en-wave-driven global solar wind model. By using in situ observations from the Solar Terrestrial Relations Observatory (STEREO) A and B, Advanced Composition Explorer (ACE), and Venus Express, we compare the observed proton state (density, temperature, and velocity) and magnetic field of the heliosphere with that predicted by the model. The validation results suggest that most of the model outputs for CR2077 can fit the observations very well. Based on this encouraging result, we therefore expect great improvement for the future modeling of coronal mass ejections (CMEs) and CME-driven shocks.
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.
Coronal mass ejections (CMEs), a kind of energetic solar eruptions, are an integral subject of space weather research. Numerical magnetohydrodynamic (MHD) modeling, which requires powerful ...computational resources, is one of the primary means of studying the phenomenon. With increasing accessibility of such resources, grows the demand for user‐friendly tools that would facilitate the process of simulating CMEs for scientific and operational purposes. The Eruptive Event Generator based on Gibson‐Low flux rope (EEGGL), a new publicly available computational model presented in this paper, is an effort to meet this demand. EEGGL allows one to compute the parameters of a model flux rope driving a CME via an intuitive graphical user interface. We provide a brief overview of the physical principles behind EEGGL and its functionality. Ways toward future improvements of the tool are outlined.
Key Points
Increase visibility of a recently developed online tool, EEGGL
Provide brief overview of the basic principles EEGGL is based on
Outline the direction of future development of EEGGL
We have recently extended the global magnetohydrodynamic (MHD) model BATS‐R‐US to account for pressure anisotropy. Since the inner magnetosphere dynamics cannot be fully described even by anisotropic ...MHD, we coupled our anisotropic MHD model with two inner magnetospheric models: the Rice Convection Model (RCM) and the Comprehensive Ring Current Model (CRCM). The coupled models provide better representations of the near‐Earth plasma, especially during geomagnetic storms. In this paper, we present the two‐way coupling algorithms with both ring current models. The major difference between these two couplings is that the RCM assumes isotropic and constant pressures along closed field lines, while the CRCM resolves pitch angle anisotropy. For model validation, we report global magnetosphere simulations performed by the coupled models. The simulation results are compared to the results given by the coupled isotropic MHD and ring current models. We find that in the global MHD simulations coupled with ring current models, pressure anisotropy results in a thinner magnetosheath, a shorter tail, a much smaller Earthward plasma jet from the tail reconnection site, and is also important in controlling the magnetic field configuration. The comparisons with satellite data for the magnetospheric event simulations show improvements on reproducing the measured tail magnetic field and inner magnetospheric flow velocity when including pressure anisotropy in the ring current model coupled global MHD model.
Key Points
We implement the coupling between Anisotropic BATS‐R‐US and ring current models.
The coupled models are applied to global magnetospheric simulations.
We compare the anisotropic and isotropic MHD simulations.
•We compare multifluid MHD and Hybrid models of the solar wind interaction with a weak comet.•The multifluid MHD results show a distinct improvement over previous single fluid simulations.•Multifluid ...MHD reproduces major features of the particle gyration as shown by the Hybrid model.•We identify a repetitive pattern of Mach cones along the cometary ion tail.•Both models show similar magnetic field draping and the formation of the anti-sunward tail.
The interaction of a comet with the solar wind undergoes various stages as the comet’s activity varies along its orbit. For a comet like 67P/Churyumov–Gerasimenko, the target comet of ESA’s Rosetta mission, the various features include the formation of a Mach cone, the bow shock, and close to perihelion even a diamagnetic cavity. There are different approaches to simulate this complex interplay between the solar wind and the comet’s extended neutral gas coma which include magnetohydrodynamics (MHD) and hybrid-type models. The first treats the plasma as fluids (one fluid in basic single fluid MHD) and the latter treats the ions as individual particles under the influence of the local electric and magnetic fields. The electrons are treated as a charge-neutralizing fluid in both cases. Given the different approaches both models yield different results, in particular for a low production rate comet. In this paper we will show that these differences can be reduced when using a multifluid instead of a single-fluid MHD model and increase the resolution of the Hybrid model. We will show that some major features obtained with a hybrid type approach like the gyration of the cometary heavy ions and the formation of the Mach cone can be partially reproduced with the multifluid-type model.
Solar energetic particles (SEPs) are an important aspect of space weather. SEP events possess a high destructive potential, since they may cause disruptions of communication systems on Earth and be ...fatal to crew members on board spacecraft and, in extreme cases, harmful to people on board high-altitude flights. However, currently the research community lacks efficient tools to predict such a hazardous threat and its potential impacts. Such a tool is a first step for mankind to improve its preparedness for SEP events and ultimately to be able to mitigate their effects. The main goal of the presented research effort is to develop a computational tool that will have the forecasting capability and can serve as an operational system that will provide live information on the current potential threats posed by SEP based on the observations of the Sun. In the present paper we discuss the fundamentals of magnetohydrodynamical simulations to be employed as a critical part of the desired forecasting system.