The role of the bow shock in solar wind-magnetosphere coupling

In this paper we examine the role of the bow shock in coupling solar wind energy to the magnetosphere using global magnetohydrodynamic simulations of the solar wind-magnetosphere interaction with southward IMF. During typical solar wind conditions, there are two significant dynamo currents in the magnetospheric system, one in the high-latitude mantle region tailward of the cusp and the other in the bow shock. As the magnitude of the (southward) IMF increases and the solar wind becomes a low Mach number flow, there is a significant change in solar wind-magnetosphere coupling. The high-latitude magnetopause dynamo becomes insignificant compared to the bow shock and a large load appears right outside the magnetopause. This leaves the bow shock current as the only substantial dynamo current in the system, and the only place where a significant amount of mechanical energy is extracted from the solar wind. That energy appears primarily as electromagnetic energy, and the Poynting flux generated at the bow shock feeds energy back into the plasma, reaccelerating it to solar wind speeds. Some small fraction of that Poynting flux is directed into the magnetosphere, supplying the energy needed for magnetospheric dynamics. Thus during periods when the solar wind flow has a low Mach number, the main dynamo in the solar wind-magnetosphere system is the bow shock

Coupling of Coronal and Heliospheric Magnetohydrodynamic Models: Solution Comparisons and Verification

Two well-established magnetohydrodynamic (MHD) codes are coupled to model the solar corona and the inner heliosphere. The corona is simulated using the MHD algorithm outside a sphere (MAS) model. The Lyon–Fedder–Mobarry (LFM) model is used in the heliosphere. The interface between the models is placed in a spherical shell above the critical point and allows both models to work in either a rotating or an inertial frame. Numerical tests are presented examining the coupled model solutions from 20 to 50 solar radii. The heliospheric simulations are run with both LFM and the MAS extension into the heliosphere, and use the same polytropic coronal MAS solutions as the inner boundary condition. The coronal simulations are performed for idealized magnetic configurations, with an out-of-equilibrium flux rope inserted into an axisymmetric background, with and without including the solar rotation. The temporal evolution at the inner boundary of the LFM and MAS solutions is shown to be nearly identical, as are the steady-state background solutions, prior to the insertion of the flux rope. However, after the coronal mass ejection has propagated through the significant portion of the simulation domain, the heliospheric solutions diverge. Additional simulations with different resolution are then performed and show that the MAS heliospheric solutions approach those of LFM when run with progressively higher resolution. Following these detailed tests, a more realistic simulation driven by the thermodynamic coronal MAS is presented, which includes solar rotation and an azimuthally asymmetric background and extends to the Earth’s orbit

Predicting magnetospheric dynamics with a coupled sun-to-Earth model: challenges and first results

Results from the first Sun-to-Earth coupled numerical model developed at the Center for Integrated Space Weather Modeling are presented. The model simulates physical processes occurring in space spanning from the corona of the Sun to the Earth's ionosphere, and it represents the first step toward creating a physics-based numerical tool for predicting space weather conditions in the near-Earth environment. Two 6- to 7-d intervals, representing different heliospheric conditions in terms of the three-dimensional configuration of the heliospheric current sheet, are chosen for simulations. These conditions lead to drastically different responses of the simulated magnetosphere-ionosphere system, emphasizing, on the one hand, challenges one encounters in building such forecasting tools, and on the other hand, emphasizing successes that can already be achieved even at this initial stage of Sun-to-Earth modeling

Geotail and LFM comparisons of plasma sheet climatology: 2. Flow variability

[1] We characterize the variability of central plasma sheet bulk flows with a 6-year Geotail data set and a 2-month Lyon-Fedder-Mobarry (LFM) global MHD simulation at two spatial resolutions. Comparing long databases of observed and simulated parameters enable rigorous statistical tests of the model\u27s ability to predict plasma sheet properties during routine driving conditions and represent a new method of global MHD validation. In this study, we use probability density functions (PDFs) to compare the statistics of plasma sheet velocities in the Geotail observations with those in the LFM simulations. We find that the low-resolution model grossly underestimates the occurrence of fast earthward and tailward flows. Increasing the simulation resolution inherently changes plasma sheet mass transport in the model, allowing the development of fast, bursty flows. These flows fill out the wings of the velocity distribution and bring the PDF into closer agreement with observations

Do we know the actual magnetopause position for typical solar wind conditions?

We compare predicted magnetopause positions at the subsolar point and four reference points in the terminator plane obtained from several empirical and numerical MHD models. Empirical models using various sets of magnetopause crossings and making different assumptions about the magnetopause shape predict significantly different magnetopause positions (with a scatter >1 RE) even at the subsolar point. Axisymmetric magnetopause models cannot reproduce the cusp indentations or the changes related to the dipole tilt effect, and most of them predict the magnetopause closer to the Earth than nonaxisymmetric models for typical solar wind conditions and zero tilt angle. Predictions of two global nonaxisymmetric models do not match each other, and the models need additional verification. MHD models often predict the magnetopause closer to the Earth than the nonaxisymmetric empirical models, but the predictions of MHD simulations may need corrections for the ring current effect and decreases of the solar wind pressure that occur in the foreshock. Comparing MHD models in which the ring current magnetic field is taken into account with the empirical Lin et al. model, we find that the differences in the reference point positions predicted by these models are relatively small for Bz=0. Therefore, we assume that these predictions indicate the actual magnetopause position, but future investigations are still needed.Key PointsEmpirical models predict significantly different magnetopause positions even at the subsolar pointAxisymmetric empirical models predict the magnetopause closer to the Earth than nonaxisymmetric empirical models for zero tilt angleResults of MHD models with the ring current magnetic field lie close to results of the nonaxisymmetric Lin et al. modelPeer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/134087/1/jgra52758_am.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/134087/2/jgra52758.pd

Geotail and LFM comparisons of plasma sheet climatology: 1. Average values

[1] We compare the statistics of central plasma sheet properties from 6 years of Geotail observations with 2 months of Lyon-Fedder-Mobarry (LFM) global MHD simulations. This statistical validation effort represents an inherently new method of systematically characterizing and quantifying global MHD model performance. For our comparison, we identify the central plasma sheet in the observations and simulation by identical criteria and select the simulation interval to ensure statistically similar distributions of solar wind conditions in both studies. After verifying our plasma sheet selection by inspecting the magnetic signatures of both studies, we compare the resultant number densities, thermal pressures, thermal energies, and bulk flows as functions of position across the equatorial plane. We find that the LFM model successfully reproduces the gross features of the global plasma sheet in a statistical sense. However, our comparison also reveals certain systematic discrepancies between the model and the observations. The LFM predicts a plasma sheet which is too dense, too cool, and exhibits faster globally averaged bulk flows than the observed plasma sheet. By quantifying the LFM overestimate of ionospheric transpolar potential and showing that ΦPC correlates with plasma sheet flow speed, we demonstrate that 15% of the plasma sheet velocity discrepancy is reflected in a ΦPC overestimate. This statistical validation effort represents an essential first step toward the rigorous, quantitative evaluation of a global MHD model in the plasma sheet

Properties of Neutral Charmed Mesons in Proton--Nucleus Interactions at 70 GeV

The results of treatment of data obtained in the SERP-E-184experiment "Investigation of mechanisms of the production of charmed particles in proton-nucleus interactions at 70 GeV and their decays" by irradiating the active target of the SVD-2 facility consisting of carbon, silicon, and lead plates, are presented. After separating a signal from the two-particle decay of neutral charmed mesons and estimating the cross section for charm production at a threshold energy {\sigma}(c\v{c})=7.1 \pm 2.4(stat.) \pm 1.4(syst.) \mub/nucleon, some properties of D mesons are investigated. These include the dependence of the cross section on the target mass number (its A dependence); the behavior of the differential cross sections d{\sigma}/dpt2 and d{\sigma}/dxF; and the dependence of the parameter {\alpha} on the kinematical variables xF, pt2, and plab. The experimental results in question are compared with predictions obtained on the basis of the FRITIOF7.02 code.Comment: 9 pages, 9 figures,3 table