Mercury’s magnetosphere-solar wind interaction for northward and southward interplanetary magnetic field: Hybrid simulation results

Analysis of global hybrid simulations of Mercury’s magnetosphere-solar wind interaction is presented for northward and southward interplanetary magnetic field (IMF) orientations in the context of MESSENGER’s first two encounters with Mercury. The global kinetic simulations reveal the basic structure of this interaction, including a bow shock, ion foreshock, magnetosheath, cusp regions, magnetopause, and a closed ion ring belt formed around the planet within the magnetosphere. The two different IMF orientations induce different locations of ion foreshock and different magnetospheric properties: the dayside magnetosphere is smaller and cusps are at lower latitudes for southward IMF compared to northward IMF whereas for southward IMF the nightside magnetosphere is larger and exhibits a thin current sheet with signatures of magnetic reconnection and plasmoid formation. For the two IMF orientations the ion foreshock and quasi-parallel magnetosheath manifest ion-beam-driven large-amplitude oscillations, whereas the quasi-perpendicular magnetosheath shows ion-temperature-anisotropy-driven wave activity. The ions in Mercury’s belt remain quasi-trapped for a limited time before they are either absorbed by Mercury’s surface or escape from the magnetosphere. The simulation results are compared with MESSENGER’s observations.     

On the origin of energetic particles in the foreshock region of the Earth’s bow shock

We consider the processes related to the formation of the so-called foreshock region upstream of the Earth’s bow shock. We suggest a model based on the surfing of pick-up ions in the bow shock front in terms of which the ion acceleration mechanism in the front can be explained. We ascertain the physical conditions under which the accelerated ions lie upstream of the shock front and determine the direction of motion of the energetic ions. We conclude that it is this population of energetic ions (longitudinal beams) that plays a major role in forming the ion foreshock boundary.     

Magnetohydrodynamic study of shock waves in the current sheet of a solar coronal magnetic loop

A Lagrangian Remap (LareXd) Code is employed to investigate the shock wave formation in the current sheet of a solar coronal magnetic loop and its effect on the magnetic reconnection. We constructed the slow shock structure in the presence of viscosity and heat conduction parallel and perpendicular to the magnetic field and pairs of slow shocks propagate away from the central current sheet, the so-called diffusion region. Significant jumps in plasma density, pressure, velocity and magnetic field occur across the main shock while the temperature appears in the foreshock. In the presence of dissipative effects, the distinct jumps disappear and the shock profiles show smooth transition between the downstream and the upstream regions while the plasma density and the pressure show very narrow and a sharp decrease with time. These results can be applied to the heating of the solar corona, the structure of the magnetic reconnection and the solar wind.     

The Martian Electron Foreshock,Based on Mars Global Surveyor's Measurements

The Martian electron foreshock position and shape have been revealed using electron flux data obtained by the Electron Reflectometer (ER) experiment on board Mars Global Surveyor, in the time period between days 87 and 118 of 1998. From the electron flux measurements, 90 inbound and outbound electron foreshock crossings were detected. Apart from their spatial distribution, their temporal variability is also studied and investigated for possible correlation with the solar wind conditions, the solar radiation flux, and especially the interplanetary magnetic field configuration at the time. The outer boundary of the ion foreshock is identified as the inner limit of the electron foreshock. The energy distribution of the foreshock electrons is also discussed.     

Excitation of alfven turbulence in the solar wind ahead of the Earth bow shock by beams of high-velocity protons

One of the mechanisms of Alfven turbulence generation in the foreshock region is investigated by the example of the Earth bow shock. The effect of the temperature of high-velocity beams on characteristics of generated disturbances is examined. It is shown that the beam temperature has a significant impact on transverse scales of disturbances. The higher the temperature, the greater the limitations on transverse wavelengths. The development of instability in the propagation of reflected, intermediate, and diffusion proton beams in the foreshock region of the Earth bow shock is considered. Perturbation motion dynamics in foreshock region is analyzed.     

A hybrid simulation of Mercury’s magnetosphere for the MESSENGER encounters in year 2008

The latest measurements from the two encounters of the MESSENGER spacecraft in year 2008 have discovered several interesting features of the magnetosphere of Mercury. We have performed high-resolution 3D hybrid model calculations to simulate the solar wind interaction with the Hermean magnetosphere during the first two Mercury encounters of the MESSENGER spacecraft in 2008. It is found that the global structure of the Hermean magnetosphere is significantly controlled by the direction of the interplanetary magnetic field. The bow shock size and shape and the magnetotail configuration have very large differences in these two encounters with northward-pointing and southward-pointing interplanetary magnetic field, respectively. Comparisons are also given with the observed magnetic field profiles and the computational results. In general, good agreement can be found including the interesting feature of the relatively thick magnetopause current layer at outbound measurements. Our work shows that 3D hybrid simulation is a promising method to study in detail the Hermean magnetosphere in parallel with the post-MOI observations of the MESSENGER spacecraft and the Bepi-Colombo mission in future.     

Multi-spacecraft study of foreshock cavitons upstream of the quasi-parallel bow shock

In this work we perform the first multi-spacecraft analysis of two foreshock cavitons observed by the Cluster spacecraft. We also study the characteristics of their surrounding regions. Foreshock cavitons are a relatively new type of phenomena in the Earth's foreshock. They appear in regions deep inside the foreshock and are therefore always immersed in a sea of ULF waves and suprathermal particles. In the observational data the cavitons appear as simultaneous depressions of interplanetary magnetic field and plasma density. The two cavitons presented here have highly structured interiors and exhibit surface irregularities. They propagate sunwards in the reference frame of the solar wind plasma. Since their velocities are smaller than the solar wind velocity, the cavitons are convected towards the Earth by the solar wind flow. Their sizes are comparable to the size of the Earth. We show that the cavitons are different from other foreshock phenomena, such as cavities. The latter are thought to form by thermal expansion due to the excess of thermal pressure caused by intense flux of suprathermal ions in their interiors. Thermal pressure inside the cavitons is the same as in their surroundings, so they cannot form in this way. The proposed mechanism for the caviton formation includes nonlinear interactions between different types of ULF waves deep inside the foreshock.     

Lunar surface electric potential changes associated with traversals through the Earth's foreshock

We report an analysis of one year of Suprathermal Ion Detector Experiment (SIDE) Total Ion Detector (TID) “resonance” events observed between January 1972 and January 1973. The study includes only those events during which upstream solar wind conditions were readily available. The analysis shows that these events are associated with lunar traversals through the dawn flank of the terrestrial magnetospheric bow shock. We propose that the events result from an increase in lunar surface electric potential effected by secondary electron emission due to primary electrons in the Earth's foreshock region (although primary ions may play a role as well). This work establishes (1) the lunar surface potential changes as the Moon moves through the terrestrial bow shock, (2) the lunar surface achieves potentials in the upstream foreshock region that differ from those in the downstream magnetosheath region, (3) these differences can be explained by the presence of energetic electron beams in the upstream foreshock region and (4) if this explanation is correct, the location of the Moon with respect to the terrestrial bow shock influences lunar surface potential.     

Energetic Particle Fluxes during the Bastille Day Solar Eruption

We report on our observations of solar energetic particle fluxes of p, He, C, O, Ne, Mg, Si, and Fe ions measured by the Energetic and Relativistic Nucleon and Electron (ERNE) experiment associated with the Bastille Day solar flare and coronal mass ejection (CME) on 14 July 2000. We observed two clear maxima of the Fe/O ratio at the energies 8.5–15 MeV nucl⁻¹. The first Fe/O maximum occurred ∼ 3 hours after the beginning of the particle event, and the second maximum ∼ 22 hours after the first one at the arrival of the shock associated with the Bastille Day eruption. We also observed a change in the energy spectrum of oxygen concurrent with a change in the direction of the interplanetary magnetic field at the start of the second enhancement of the Fe/O ratio. We propose an interpretation of the particle event where observed interplanetary particle fluxes are associated with two different particle sources near the Sun and in interplanetary space. We suggest that heavy ions observed during the first period of the Fe/O enhancement were released when a coronal shock reached a magnetic foot point connected to 1 AU. The second maximum of Fe/O occurred when spacecraft encountered Fe-rich material stored in magnetic field flux tubes early in the event and was possibly reaccelerated by the interplanetary shock.     

Perturbative analysis of two-temperature radiative shocks with multiple cooling processes

The structure of the hot downstream region below a radiative accretion shock, such as that of an accreting compact object, may oscillate because of a global thermal instability. The oscillatory behaviour depends on the functional forms of the cooling processes, the energy exchanges of electrons and ions in the shock-heated matter, and the boundary conditions. We analyse the stability of a shock with unequal electron and ion temperatures, where the cooling consists of thermal bremsstrahlung radiation which promotes instability, plus a competing process which tends to stabilize the shock. The effect of transverse perturbations is considered also. As an illustration, we study the special case in which the stabilizing cooling process is of order 3/20 in density and 5/2 in temperature, which is an approximation for the effects of cyclotron cooling in magnetic cataclysmic variables. We vary the efficiency of the second cooling process, the strength of the electron–ion exchange and the ratio of electron and ion pressures at the shock, to examine particular effects on the stability properties and frequencies of oscillation modes.     

The Bastille day Shock and Merged Interaction Region at 63 au: Voyager 2 Observations

A transient flow system containing several streams and shocks associated with the Bastille Day 2000 solar event was observed by the WIND and ACE spacecraft at 1 AU. Voyager 2 (V2) at 63 AU observed this flow system after it moved through the interplanetary medium and into the distant heliosphere, where the interstellar pickup protons strongly influence the MHD structures and flow dynamics. We discuss the Voyager 2 magnetic and plasma observations of this event. Increases in the magnetic field strength B, density N, temperature T and speed V were observed at the front of a stream at V2, consistent with presence of a shock related to the Bastille Day shock at 1 AU. However, the jumps occurred in a 16.9-hour data gap, so that the shock was not observed directly, and the properties of the candidate shock cannot be determined precisely. The candidate shock was followed by a merged interaction region (MIR) that moved past V2 for at least 10 days. The first part of this MIR contains a structure that might be a magnetic cloud. Just ahead of the shock there was an abrupt increase in density associated with a decrease in temperature such that the solar wind thermal pressure was constant across it. Just behind the shock there was an abrupt decrease in density associated with a net increase in magnetic field strength. This appears to be a pressure balanced structure in which the interstellar pickup protons make a significant contribution.     

Field-aligned beams and reconnection in the jovian magnetotail

The release of plasma in the jovian magnetotail is observed in the form of plasmoids, travelling compression regions, field-aligned particle beams and flux-rope like events. We demonstrate that electrons propagate along the magnetic field lines in the plasma sheet boundary layer (PSBL), while close to the current sheet center the electron distribution is isotropic. The evidences of the counterstreaming electron beams in the PSBLs are also presented. Most of the field-aligned energetic ion beams are associated with the field-aligned electron beams and about half of them have the bipolar fluctuation of the meridional magnetic field component. Moreover they often show a normal velocity dispersion for the different species which fits well in the scenario of particle propagation from a single source. All features above are observed during jovian reconfiguration events which are typically bonded with plasma flow reversals. From all these characteristics, which are based on energetic particle and magnetic field measurements, we believe that the reconfiguration processes in the jovian magnetotail are associated with reconnection.     

The structure of radiative slow-mode shocks

We investigate the structure of slow-mode MHD shocks in a plasma where both radiation and thermal conduction are important. In such a plasma a slow shock dissociates into an extended foreshock, an isothermal subshock, and a downstream radiative cooling region. Our analysis, which is both numerical and analytical, focuses on the nearly switch-off shocks which are generated by magnetic reconnection in a strong magnetic field. These shocks convert magnetic energy into kinetic energy and heat, and we find that for typical flare conditions about f of the conversion occurs in the subshock while the remaining 1/3 occurs in the foreshock. We also find that no stable, steady-state solutions exist for radiative slow shocks unless the temperature in the radiative region downstream of the subshock falls below 10⁵ K. These results suggest that about 2/3 of the magnetic energy released in flare loops is released at the top of the loop, while the remaining 1/3 is released in the legs of the loop.     

The interaction of cosmic rays and magnetized plasma

A large flux of cosmic rays streaming through a magnetized plasma creates cavities of low plasma density and low magnetic field. The magnetic field focuses the cosmic ray trajectories into the cavities with the possible formation of filaments or beams of high-energy cosmic rays.     

Nonlinear Force-Free Reconstruction of the Global Solar Magnetic Field: Methodology

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday’s equation, give rise to a respective normal-field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modeled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal-field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition – the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time that a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but it depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.     

Acceleration of Pick-up Ions at the Solar Wind Termination Shock

It is generally accepted that pick-up ions act as a seed population for anomalous cosmic rays originating at the solar wind termination shock. We believe that the ion pre-acceleration process operating in the heliosphere up to the termination shock can be very important to inject the ions into the shock acceleration process. The pick-up ions pre-accelerated by solar wind turbulences have already a pronounced high energy tail when they reach the shock. Some fraction of these ions can experience further acceleration up to energies of anomalous cosmic rays by means of shock drift and diffusive acceleration. In the present paper the shock drift acceleration of pick-up ions suffering multiple reflection due to abrupt changes in both the strength and direction of the magnetic field through the shock is considered. The reflection process operates for high velocity particles different from the reflection by the electric cross-shock potential. During the first reflection the mean kinetic energy of pick-up ions increases by approximately a factor of 10. Reflected particles have highly anisotropic velocity distribution. Subsequent excursion of the particles in the turbulent upstream flow leads to diffusion in pitch-angle space and, as a result, the particles can return to the shock again suffering, thus, multiple encounters. In order to describe the motion of particles in the upstream and down streamparts of the flow we solve the Fokker-Plank transport equation for anisotropic velocity distribution function. This revised version was published online in July 2006 with corrections to the Cover Date.     

Flares in Sigmoidal Coronal Structures – a Case Study

We analyze radio observations, magnetograms and extrapolated field line maps, Hα filtergrams, and X-ray observations of two flare events (6 February 1992 in AR 7042 and 25 October 1994 in AR 7792) and study properties, evolution and energy release signatures of sigmoidal loop systems. During both events, the loop configuration seen in soft X-ray (SXR) images changes from a preflare sigmoidal shape to a relaxed post-flare loop system. The underlying magnetic field system consists of a quadrupolar configuration formed by a sheared arcade core and a remote field concentration. We demonstrate two possibilities: a sigmoidal SXR pattern can be due to a single continuous flux tube (the 1992 event). Alternatively, it can be due to a set of independent loops appearing like a sigmoid (the 1994 event). In both cases, the preflare and post-flare loops can be well reproduced by a linear force-free field and potential field, respectively, computed using preflare magnetograms. We find that thermal and non-thermal flare energy release indicators of both events become remarkably similar after applying spatial and temporal scale transformations. Using the spatial scaling between both events we estimated that the non-thermal energy release in the second event liberated about 1.7 times more energy per unit volume. A two-and-a-half times faster evolution indicates that the rate of the energy release per unit volume is more than four times higher in this event. A coronal type II burst reveals ignition and propagation of a coronal shock wave. In contrast, the first event, which was larger and released about a 10 times more energy during the non-thermal phase, was associated with a CME, but no type II burst was recorded. During both events, in addition to the two-ribbon flare process an interaction was observed between the flaring arcade and an emerging magnetic flux region of opposite polarity next to the dominant leading sunspot. The arcade flare seems to stimulate the reconnection process in an `emerging flux-type' configuration, which significantly contributes to the energy release. This regime is characterized by the quasiperiodic injection of electron beams into the surrounding extended field line systems. The repeated beam injections excite pulsating broadband radio emission in the decimetric-metric wavelength range. Each radio pulse is due to a new electron beam injection. The pulsation period (seconds) reflects the spatial scale of the emerging flux-type field configuration. Since broadband decimetric-metric radio pulsations are a frequent radio flare phenomenon, we speculate that opposite-polarity small-scale flux intrusions located in the vicinity of strong field regions may be an essential component of the energy release process in dynamic flares.     

Generation of low-frequency magnetic field disturbances in coronal loops by proton and electron beams

We consider the generation of low-frequency magnetic field disturbances in coronal loops when low density proton and electron beams propagate in them. Two mechanisms of low-frequency magnetic field perturbation generation are analyzed. The first mechanism is concerned with the longitudinal current generated by charged particles’ beams moving in the loop. It is shown that this mechanism of the Alfvén waves’ generation can lead to development of low-frequency perturbations even if the currents are very weak. It can facilitate the reconnection of magnetic fields and flare development. The second mechanism is not concerned with currents propagating in the coronal loop. It is shown that, in this case, the proton beams can cause an instability with significantly lower values of beam density. We found increments and criteria of the development of instabilities. Not only Alfvén-type perturbations can be generated as a result of development of those instabilities but also kinetic Alfvén-type perturbations can be generated.     

Transequatorial Filament Eruption and Its Link to a Coronal Mass Ejection

We revisit the Bastille Day flare/CME Event of 2000 July 14, and demonstrate that this flare/CME event is not related to only one single active region (AR). Activation and eruption of a huge transequatorial filament are seen to precede the simultaneous filament eruption and flare in the source active region, NOAA AR 9077, and the full halo-CME in the high corona. Evidence of reconfiguration of large-scale magnetic structures related to the event is illustrated by SOHO EIT and Yohkoh SXT observations, as well as, the reconstructed 3D magnetic lines of force based on the force-free assumption. We suggest that the AR filament in AR9077 was connected to the transequatorial filament. The large-scale magnetic composition related to the transequatorial filament and its sheared magnetic arcade appears to be an essential part of the CME parent magnetic structure. Estimations show that the filament-arcade system has enough magnetic helicity to account for the helicity carried by the related CMEs. In addition, rather global magnetic connectivity, covering almost all the visible range in longitude and a huge span in latitude on the Sun, is implied by the Nancay Radioheliograph (NRH) observations. The analysis of the Bastille Day event suggests that although the triggering of a global CME might take place in an AR, a much larger scale magnetic composition seems to be the source of the ejected magnetic flux, helicity and plasma. The Bastille Day event is the first described example in the literature, in which a transequatorial filament activity appears to play a key role in a global CME. Many tens of halo-CME are found to be associated with transequatorial filaments and their magnetic environment.     

Magnetospheric Cusps under Extreme Conditions: Cluster Observations and MHD Simulations Compared

MHD simulations are here applied to aid in the interpretation of three apparent cusp encounters by the Cluster 4 spacecraft in unusual places when the magnetosphere was under extreme solar wind and interplanetary magnetic field (IMF) conditions associated with the passage of magnetic clouds imbedded within fast ICMEs. At the time of each cusp encounter the IMF was very strong, generally northward in one case, generally equatorial in a second case, and generally southward in the third case. In the southward IMF case, the MHD models locate the origin of the cusp-like plasma by showing that the position of the spacecraft at the time of encounter was engulfed in a tongue of high-pressure plasma extending from the magnetopause into the magnetosphere. This tongue points to the northern-hemisphere cusp as the source of the feature. In the equatorial IMF case an elevated-pressure feature that apparently marked a cusp encounter in the computations coincided, however, with a passage in the solar wind of a dynamic pressure pulse, thus giving an alternative interpretation of the feature. However, Cluster data unambiguously identified the event as an encounter with magnetosheath-like plasma. Given that the Cluster observations classify the event as a true encounter with a cusp-like plasma feature (and not a compression event), the model simulations can be interpreted as identifying the origin of the feature to have been the northern-hemisphere cusp even though?—?and this is the interesting point?—?the observation point was in the southern hemisphere. In the northward IMF case, neither cusp (defined as a magnetic funnel linking the magnetopause to the Earth) was directly connected to the observation point. Instead, this encounter of magnetosheath-like plasma appears to be an instance of boundary-layer formation by means of the Song?–?Russell mechanism in which two-point magnetic reconnection entrains magnetosheath plasma on closed field lines when the IMF is northward.