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.     

Deep Space 1 at comet 19P/Borrelly: Magnetic field and plasma observations

On September 22, 2001 the Deep Space 1 spacecraft performed a flyby at comet 19P/Borrelly at a solar distance of 1.36 AU leading the Earth by 74^° in longitude. The spacecraft-comet distance at closest approach was 2171 km. The bow shock had a magnetic compression ratio of 2.5 at a distance of 147 100 km from the nucleus. Deep Space 1 first entered the sheath region essentially from the north polar region. Fluctuations from the cometary ion pickup were present throughout the sheath region and even well upstream of the shock, as expected. The magnetic field pileup region had a peak field strength of 83 nT and was shown to be consistent with a pressure equal to the solar wind ram pressure. The peak field location was offset from the time of closest approach. It is uncertain whether this is a spatial or temporal variation. Draping of magnetic fields around the nucleus was sought, but evidence for this was not apparent in the data. A possible explanation is that the interplanetary solar wind was composed of turbulent short-scale fields, and thus the fields were not symmetric about the point of closest approach. During the flyby phase there were in general few intervals of ACE data where there were large scale Parker spiral fields. With the addition of plasma data, the shock properties are investigated. The characteristics of magnetic draping, pileup and fluctuations are explored. These comet 19P/Borrelly results are contrasted with other cometary flyby results.     

Multi-Spacecraft Study of the 21 January 2005 ICME

We examine the near-Earth Interplanetary Coronal Mass Ejection (ICME) apparently related to the intense Solar Energetic Particle (SEP) event of 20 January 2005. Our purpose is to contribute to the understanding of the macroscopic structure, evolution and dynamics of the solar corona and heliosphere. Using Cluster, ACE and Wind data in the solar wind, and Geotail data in the magnetosheath, we perform a multi-spacecraft analysis of the ICME-driven shock, post-shock magnetic discontinuities and ejecta. Traversals by the well-separated near-Earth spacecraft provide a coherent picture of the ICME geometry. Following the shock, the ICME sequence starts with a hot pileup, i.e.,? a sheath, followed by a fast ejecta characterised by a non-compressive density enhancement (NCDE), which is caused essentially by an enrichment in helium. The plasma and magnetic observations of the ejecta are consistent with the outskirts of a structure in strong expansion, consisting of nested magnetic loops still connected to the Sun. Within the leading edge of the ejecta, we establish the presence of a tilted current sheet substructure. An analysis of the observations suggests that the tilted current sheet is draped within the overlying cloud canopy, ahead of a magnetic cloud-like structure. The flux rope interpretation of this structure near L1, confirmed by observations of the corresponding magnetic cloud, provided by Ulysses at 5.3 AU and away from the Sun?–?Earth line, indicates that the bulk of the cloud is in the northwest sector as seen from the Earth, with its axis nearly perpendicular to the ecliptic. This is consistent with the primary direction of travel of the fast halo CME observed at the Sun. Moreover, the NCDE and helium enrichment are consistent with the position near the streamer belt of the flaring active region NOAA 10720 associated with the CME. However, differences between interplanetary and solar observations indicate a large rotation of the erupting filament and overlying arcade, which can be attributed to the flux rope being subject to the helical kink instability.     

Solar Energetic Particle Events in the 23rd Solar Cycle: Interplanetary Magnetic Field Configuration and Statistical Relationship with Flares and CMEs

We study the influence of the large-scale interplanetary magnetic field configuration on the solar energetic particles (SEPs) as detected at different satellites near Earth and on the correlation of their peak intensities with the parent solar activity. We selected SEP events associated with X- and M-class flares at western longitudes, in order to ensure good magnetic connection to Earth. These events were classified into two categories according to the global interplanetary magnetic field (IMF) configuration present during the SEP propagation to 1 AU: standard solar wind or interplanetary coronal mass ejections (ICMEs). Our analysis shows that around 20 % of all particle events are detected when the spacecraft is immersed in an ICME. The correlation of the peak particle intensity with the projected speed of the SEP-associated coronal mass ejection is similar in the two IMF categories of proton and electron events, ≈?0.6. The SEP events within ICMEs show stronger correlation between the peak proton intensity and the soft X-ray flux of the associated solar flare, with correlation coefficient r=0.67±0.13, compared to the SEP events propagating in the standard solar wind, r=0.36±0.13. The difference is more pronounced for near-relativistic electrons. The main reason for the different correlation behavior seems to be the larger spread of the flare longitude in the SEP sample detected in the solar wind as compared to SEP events within ICMEs. We discuss to what extent observational bias, different physical processes (particle injection, transport, etc.), and the IMF configuration can influence the relationship between SEPs and coronal activity.     

The mechanism of magnetospheric substorm and the MHD waves of the solar wind

A mechanism of the Earth's magnetospheric substorm is proposed. It is suggested that the MHD waves may propagate across the magnetopause from the magnetosheath into the magnetotail and will be dissipated in the plasma sheet, heating the plasma and accelerating the particles. When the solar wind parameters change, the Poynting flux of the waves transferred from the magnetosheath into the tail, may be greater than 10¹⁸ erg s^?1. The heated plasma and accelerated particles in the plasma sheet will be injected into the inner magnetosphere, and this may explain the process of the ring current formation and auroral substorm.The Alfvén wave can only propagate along the magnetic force line into the magnetosphere in the open magnetosphere, but the magnetosonic wave can propagate in both the open and closed magnetosphere. When the IMF turns southward, the configuration of the magnetosphere will change from a nearly closed model into some kind of open one. The energy flux of Alfvén waves is generally larger than that of the magnetosonic wave. This implies that it is easy to produce substorms when the interplanetary magnetic field (IMF) has a large southward component, but the substorm can also be produced even if the IMF is directed northward.     

On one mechanism for the generation of a reverse solar wind shock in the magnetosheath in front of the Earth’s magnetosphere

A possible mechanism for the generation of a reverse fast shock in the magnetosheath in the solar wind flow around the Earth’s magnetosphere is considered. It is shown that such a shock can emerge through the breaking of a nonlinear fast magnetosonic compression wave reflected from the magnetopause toward the bow shock rear. In this case, the magnetopause is represented as a tangential discontinuity with a zero normal magnetic field component at it and the mechanism under consideration is assumed to be secondary with respect to the sudden disturbance of the bow shock-Earth’s magnetosphere system by a nonstationary solar wind shock. A possible confirmation of the process under study by in-situ SC3 experimental observations of the bow shock front motion on the Cluster spacecraft is pointed out.     

Tracking a Ulysses High-latitude ICME Event Back to Its Solar Origins

High-latitude interplanetary mass ejections (ICMEs) observed beyond 1 AU are not studied very often. They are useful for improving our understanding of the 3D heliosphere. As there are only few such events registered by the Ulysses spacecraft, the task of detecting their solar counterparts is a challenge, especially during high solar activity periods, because there are dozens coronal mass ejections (CMEs) registered by SOHO that might be chosen as candidates. We analyzed a high-latitude ICME registered by the Ulysses spacecraft on 18 January 2002. Our investigation focused on the correlation between various plasma parameters that allow the identification to be made of the ICME and its components such as the forward shock, the magnetic cloud and the reverse shock.     

Sector-structured interplanetary magnetic field associated with the fast plasma streams in 1985–1996

An analysis of 373 well-defined high-speed solar-wind streams observed at 1 AU during the years 1985–1996 is outlined. The distribution of the occurrence of these streams as a function of Bartels rotation days using the dominant polarity of the interplanetary magnetic field (IMF) associated with the referred fast streams shows that a four-sector pattern for the positive IMF polarity and a two-sector pattern for the negative IMF polarity are the dominant features in the investigated period. The high-speed streams seem to occur at preferred Bartels days: positive polarity streams are most frequent near Bartels days 5 and 18, while negative polarity streams are most frequent in days 14 and 23. Moreover, the corotating streams with positive IMF polarity prefer to occur in days 5 and 18 of the Bartels rotation period, whereas flare-generated streams with negative IMF polarity occur in days 14 and 23. The observed distribution of Bartels days is probably related to the distribution of the solar sources of high-speed solar wind streams as the solar wind carries with it the photospheric magnetic polarity of the solar source region. In addition, the distribution of the streams reveals a similar behaviour during the ascending and the declining phase of the last solar cycle (22nd) in contrast to the previous one where it has an opposite appearance. Determined differences in the characteristics of the sector structured IMF associated with the fast streams of the last cycle with the previous one (21st) and some similarities with the alternate solar cycle (20th) seem to be attributed to the 22-year magnetic cycle and to the polarity reversals of the polar magnetic field of the Sun. As the magnetic sectors are due to multiple crossings of the solar equatorial plane by a large-scale, warped heliospheric current sheet, it is suggested that the two-sector pattern arises from a tilted solar magnetic dipole component and the more commonly observed four-sector pattern from a quadrupole component of the solar interplanetary magnetic field.     

Effects of the Weak Polar Fields of Solar Cycle 23: Investigation Using OMNI for the STEREO Mission Period

The current solar cycle minimum seems to have unusual properties that appear to be related to weak solar polar magnetic fields. We investigate signatures of this unusual polar field in the ecliptic near-Earth interplanetary magnetic field (IMF) for the STEREO period of observations. Using 1 AU OMNI data, we find that for the current solar cycle declining phase to minimum period the peak of the distribution for the values of the ecliptic IMF magnitude is lower compared to a similar phase of the previous solar cycle. We investigate the sources of these weak fields. Our results suggest that they are related to the solar wind stream structure, which is enhanced by the weak polar fields. The direct role of the solar field is therefore complicated by this effect, which redistributes the solar magnetic flux at 1 AU nonuniformly at low to mid heliolatitudes.     

On the relationship of cosmic ray intensity with solar,interplanetary, and geophysical parameters

The flux rate of cosmic rays incident on the Earth’s upper atmosphere is modulated by the solar wind and the Earth’s magnetic field. The amount of solar wind is not constant due to changes in solar activity in each solar cycle, and hence the level of cosmic ray modulation varies with solar activity. In this context, we have investigated the variability and the relationship of cosmic ray intensity with solar, interplanetary, and geophysical parameters from January 1982 through December 2008. Simultaneous observations have been made to quantify the exact relationship between the cosmic ray intensity and those parameters during the solar maxima and minima, respectively. It is found that the stronger the interplanetary magnetic field, solar wind plasma velocity, and solar wind plasma temperature, the weaker the cosmic ray intensity. Hence, the lowest cosmic ray intensity has good correlations with simultaneous solar parameters, while the highest cosmic ray intensity does not. Our results show that higher solar activity is responsible for a higher geomagnetic effect and vice versa.     

Interplanetary Coronal Mass Ejections Observed by <Emphasis Type="Italic">Ulysses</Emphasis> Through Its Three Solar Orbits

An extended Ulysses interplanetary coronal mass ejections (ICMEs) list is used to statistically study the occurrence rate of ICMEs, the interaction of ICMEs with solar wind, and the magnetic field properties in ICMEs. About 43% of the ICMEs are identified as magnetic clouds (MCs). It is found that the occurrence rate of ICMEs approximately follows the solar activity level, except for the second solar orbit; the rate is higher in the southern heliolatitude than in the northern heliolatitude; and it roughly decreases with the increase of ICME speeds. Our results show that the speed difference between the ICME and the solar wind in front of it shows a slight decrease with increasing heliocentric distance for ICMEs preceded by a shock, whereas no such dependence is found for the ICMEs without shock association. It is also found that approximately 23% of the ICMEs are associated with radial field events, in which the interplanetary magnetic field with near-radial direction lasts for many hours, in the Ulysses detected range, and these associated events are not necessarily confined to fast ICMEs or the trailing portions of ICMEs. Nearly all these associated events occur during the period of declining solar wind speed and most of them occur at low heliolatitudes.     

An indication for an intrinsically different coronal temperature

As is well known, the interplanetary magnetic field is the extension of the solar magnetic field, and the solar wind is the extension of the solar corona (Biermann, 1951; Parker, 1958). Consequently, knowledge of these two interplanetary structures reflects a knowledge of the solar corona itself.Calculations of radial proton temperature gradients within magnetic sector boundaries, observed by Helios-1 and Helios-2 heliocentric satellites between 0.3 and 1.0 AU, revealed lower temperatures in the solar corona in the region where these boundaries merge.     

Study of the Cosmic-Ray Modulation During the Passage of ICMEs and CIRs

We compare the cosmic-ray response to interplanetary coronal mass ejections (ICMEs) and corotating interaction regions (CIRs) during their passage in near-Earth space. We study the relative importance of various structures/features identified during the passage of the ICMEs and CIRs observed during Cycle 23 (1995?–?2009). The identified ICME structures are the shock front, the sheath, and the CME ejecta. We isolate the shock arrival time, the passage of the sheath region, the arrival of ejecta, and the end time of their passage. Similarly, we isolate the CIR arrival, the associated forward shock, the stream interface, and the reverse shock during the passage of a CIR. For the cosmic-ray intensity, we utilize the data from high counting rate neutron monitors. In addition to neutron monitor data, we utilize near-simultaneous and same time-resolution data of interplanetary plasma and field, namely the solar-wind velocity, the interplanetary magnetic field (IMF) vector, and its variance. Further, we also utilize some derived interplanetary parameters. We apply the method of the superposed-epoch analysis. As the plasma and field properties are different during the passage of different structures, both in ICMEs and CIRs, we systematically vary the epoch time in our superposed-epoch analysis one by one. In this way, we study the role and effects of each of the identified individual structures/features during the passage of the ICMEs and CIRs. Relating the properties of various structures and the corresponding variations in plasma and field parameters with changes of the cosmic-ray intensity, we identify the relative importance of the plasma/field parameters in influencing the amplitude and time profiles of the cosmic-ray intensity variations during the passage of the ICMEs and CIRs.     

The Bastille day Magnetic Clouds and Upstream Shocks: Near-Earth Interplanetary Observations

The energetic charged particle, interplanetary magnetic field, and plasma characteristics of the `Bastille Day' shock and ejecta/magnetic cloud events at 1 AU occurring over the days 14–16 July 2000 are described. Profiles of MeV (WIND/LEMT) energetic ions help to organize the overall sequence of events from the solar source to 1 AU. Stressed are analyses of an outstanding magnetic cloud (MC2) starting late on 15 July and its upstream shock about 4 hours earlier in WIND magnetic field and plasma data. Also analyzed is a less certain, but likely, magnetic cloud (MC1) occurring early on 15 July; this was separated from MC2 by its upstream shock and many heliospheric current sheet (HCS) crossings. Other HCS crossings occurred throughout the 3-day period. Overall this dramatic series of interplanetary events caused a large multi-phase magnetic storm with min Dst lower than −300 nT. The very fast solar wind speed (≥ 1100 km s⁻¹) in and around the front of MC2 (for near average densities) was responsible for a very high solar wind ram pressure driving in the front of the magnetosphere to geocentric distances estimated to be as low as ≈ 5 R _E, much lower than the geosynchronous orbit radius. This was consistent with magnetic field observations from two GOES satellites which indicated they were in the magnetosheath for extended times. A static force-free field model is used to fit the two magnetic cloud profiles providing estimates of the clouds' physical and geometrical properties. MC2 was much larger than MC1, but their axes were nearly antiparallel, and their magnetic fields had the same left-handed helicity. MC2's axis and its upstream shock normal were very close to being perpendicular to each other, as might be expected if the cloud were driving the shock at the time of observation. The estimated axial magnetic flux carried by MC2 was 52×10²⁰ Mx, which is about 5 times the typical magnetic flux estimated for other magnetic clouds in the WIND data over its first 4 years and is 17 times the flux of MC1. This large flux is due to both the strong axially-directed field of MC2 (46.8 nT on the axis) and the large radius (R ₀=0.189 AU) of the flux tube. MC2's average speed is consistent with the expected transit time from a halo-CME to which it is apparently related.     

Radial deformation of the solar current sheet as a cause of geomagnetic storms

It is suggested that the solar current sheet, extending from a coronal streamer, develops a large-scale radial deformation, at times with a very steep gradient at the Earth's distance. The associated magnetic field lines (namely, the interplanetary magnetic field (IMF) lines) are expected to have also a large gradient in the vicinity of the current sheet. It is also suggested that some of the major geomagnetic storms occur when the Earth is located in the region where IMF field lines have a large dip angle with respect to the ecliptic plane for an extended period (6–48 h), as a result of a steep radial deformation of the current sheet.     

Asymmetry of plasma fluxes at Mars. ASPERA-3 observations and hybrid simulations

The asymmetry of fluxes of solar wind and planetary ions is studied by using the ASPERA-3 observations onboard the Mars Express spacecraft in February 2004 to March 2006. Due to the small scale of the Martian magnetosphere and its induced origin, the flow pattern near Mars is sensitive to the directions of the interplanetary magnetic and electric (-V×B) fields. Asymmetry of the magnetic field draping produces an asymmetry in plasma flows in the plane containing the IMF. The crustal magnetic fields on Mars also influence the flow pattern. Scavenging of planetary ions is less efficient in the regions of strong crustal magnetization and therefore the escape fluxes of planetary ions in the southern hemisphere are smaller. The results of the observations are compared to simulations based on a 3D hybrid model with several ion species.     

The magnetic field near Mars: A comparison between a hybrid model, Mars Global Surveyor and Mars Express observations

Mars Express (MEX) does not carry its own magnetometer which complicates interpretation of ASPERA-3/MEX ion measurements. The direction of the interplanetary magnetic field (IMF) is especially important because it, among other things, determines the direction of the convective electric field and orientation of the cross tail current sheet and tail lobes. In this paper we present a case study to show the properties of the magnetic field near Mars in a quasi-neutral hybrid (QNH) model at the orbits where the Mars Global Surveyor (MGS) has made measurements, present a method to derive the IMF clock angle by comparing fields in a hybrid model and the direction of the magnetic field measured by MGS by deriving the IMF clock angle. We also use H⁺ ring velocity distribution observations upstream of the bow shock measured by the IMA/ASPERA-3 instrument on board MEX spacecraft. These observations are used to indirectly provide the orientation of the IMF. We use a QNH model (HYB-Mars) where ions are modeled as particles while electrons form a mass-less charge neutralizing fluid. We found that the direct MGS and non-direct IMA observations of the orientation magnetic field vectors in non-crustal magnetic field regions are consistent with the global magnetic field draping pattern predicted by the global model.     

On the association of predominant polarity of North-South component of IMF in the GSE system near 1 AU with solar polar magnetic fields during 1967–2020

The skewness of the monthly distribution of GSE latitudinal angles of Interplanetary Magnetic Field (IMF) observed near the Earth (Sk) is found to show anti-correlation with sunspot activity during the solar cycles 20–24. Sk can be considered as a measure of the predominant polarity of north-south component of IMF (B_z component) in the GSE system near 1 AU. Sk variations follow the magnitude of solar polar magnetic fields in general and polarity of south polar fields in particular during the years 1967–2020. Predominant polarity of Sk is found to be independent of the heliographic latitude of Earth. Sk basically reflects the variations of the solar dipolar magnetic field during a sunspot cycle. It is also found that IMF sector polarity variation is not a good indicator of the magnitude changes in solar polar magnetic fields during a sunspot cycle. This is possibly due to the influence of non-dipolar components of the solar magnetic field and the associated north-south asymmetries in the heliospheric current sheet.     

A global 3-D simulation of interplanetary dynamics in June 1991

The global dynamics of the solar wind and interplanetary magnetic field in June 1991 is simulated based on a fully three-dimensional, time-dependent numerical MHD model. The numerical simulation includes eight transient disturbances associated with the major solar flares of June 1991. The unique features of the present simulation are: (i) the disturbances are originated at the coronal base (1R _s) and their propagation through inhomogeneous ambient solar wind is simulated out to 1.5 AU; (ii) as a background for the transients, the global steady-state solar wind structure inferred from the 3-D steady-state model (Usmanov, 1993c) is used. The parameters of the initial pulses are prescribed in terms of the near-Sun shock velocities (as inferred from the metric Type II radio burst observations) relative to the preshock steady-state flow parameters at the flare sites. The computed parameters at the Earth's location for the period 1–18 June, 1991 are compared with the available observations of the interplanetary magnetic field, solar wind velocity, density, and with variation of the geomagnetic activityK _pindex.     

Deformation of ICMEs/MCs along their path

Interplanetary coronal mass ejections (ICMEs) and their subset, magnetic clouds (MCs), are important manifestations of solar activity which have substantial impact on the geomagnetic field. We re-analyze events already identified in Wind and Voyager 2 data and estimate changes of their geometry along the path from the Sun. The analysis is based on the thickness of the sheath between a shock and a particular ICME or MC which is proportional to the apparent curvature radius of ICMEs/MCs. We have found that this apparent radius of curvature increases with the Mach number and this effect is attributed to the larger deformation of the fast ICME/MC. Further, the relative sheath thickness that is proportional to the flux rope oblateness decreases with the magnetic field intensity inside the ICME/MC and increases with the heliospheric distance.