Fluctuations of cosmic rays and IMF in the vicinity of interplanetary shocks

Fluctuations of cosmic rays and interplanetary magnetic field upstream of interplanetary shocks are studied using data of ground-based polar neutron monitors as well as measurements of energetic particles and solar wind plasma parameters aboard the ACE spacecraft. It is shown that coherent cosmic ray fluctuations in the energy range from 10 keV to 1 GeV are often observed at the Earth’s orbit before the arrival of interplanetary shocks. This corresponds to an increase of solar wind turbulence level by more than the order of magnitude upstream of the shock. We suggest a scenario where the cosmic ray fluctuation spectrum is modulated by fast magnetosonic waves generated by flux of low-energy cosmic rays which are reflected and/or accelerated by an interplanetary shock.     

X-ray, γ-emission and energetic particles in near-Earth space as measured by CORONAS-F satellite: From maximum to minimum of the last solar cycle

The Russian solar observatory CORONAS-F was launched into a circular orbit on July 31, 2001 and operated until December 12, 2005. Two main aims of this experiment were: (1) simultaneous study of solar hard X-ray and γ-ray emission and charged solar energetic particles, (2) detailed investigation of how solar energetic particles influence the near-Earth space environment. The CORONAS-F satellite orbit allows one to measure both solar energetic particle dynamics and variations of the solar particle boundary penetration as well as relativistic electrons of the Earth’s outer radiation belt during and after magnetic storms. We have found that significant enhancements of relativistic electron flux in the outer radiation belt were observed not only during strong magnetic storms near solar maximum but also after weak storms caused by high speed solar wind streams. Relativistic electrons of the Earth’s outer radiation belt cause volumetric ionization in the microcircuits of spacecraft causing them to malfunction, and solar energetic particles form an important source of radiation damage in near-Earth space. Therefore, the present results and future research in relativistic electron flux dynamics are very important.     

Solar wind electron density and temperature over solar cycle 23: Thermal noise measurements on Wind

We present the solar wind plasma parameters obtained from the Wind spacecraft during more than nine years, encompassing almost the whole solar cycle 23. Since its launch in November 1994 Wind has frequently observed the in-ecliptic solar wind upstream of the Earth’s bow shock. The WIND/WAVES thermal noise receiver was specially designed to measure the in situ plasma thermal noise spectra, from which the electron density and temperature can be accurately determined. We present and discuss histograms of such measurements performed from 1994 to 2003. Using these large data sets, we study the density and core temperature variations with solar activity cycle and with different regimes of the solar wind. We confirm the anticorrelation of the electron density with the sunspot number, and obtain a positive correlation of the core temperature, with the sunspot number.     

Sharp boundaries of solar wind plasma structures and their relationship to solar wind turbulence

Sharp (<10 min) and large (>20%) solar wind ion flux changes are common phenomena in turbulent solar wind plasma. These changes are the boundaries of small- and middle-scale solar wind plasma structures which can have a significant influence on Earth’s magnetosphere. These solar wind ion flux changes are typically accompanied by only a small change in the bulk solar wind velocity, hence, the flux changes are driven mainly by plasma density variations. We show that these events occur more frequently in high-density solar wind. A characteristic of solar wind turbulence, intermittency, is determined for time periods with and without these flux changes. The probability distribution functions (PDF) of solar wind ion flux variations for different time scales are calculated for each of these periods and compared. For large time scales, the PDFs are Gaussian for both data sets. For small time scales, the PDFs from both data set are more flat than Gaussian, but the degree of flatness is much larger for the data near the sharp flux change boundaries.     

Solar wind interaction with comet Halley

In March 6 and 9, 1986 the spacecrafts ‘Vega-1’ and ‘Vega-2’ have flown through the coma of comet Halley and have carried measurements of plasma, energetic particles, magnetic field and plasma waves along its trajectory. A short review of these measurements and its comparison with theoretical models of solar wind interaction with comets are given.

The spacecrafts ‘Vega-1’ and ‘Vega-2’ have studied the solar wind loading by cometary ions, the structure of cometary bow shock and the processes in the inner coma of comet Halley. Exactly in this sequence we discuss the results of measurements and compare them with the theory.     

Energetic ion acceleration and transport in the upstream region of Jupiter: Voyager 1 and 2

Long-lived upstream energetic ion events at Jupiter appear to be very similar in nearly all respects to upstream ion events at earth. A notable difference between the two planetary systems is the enhanced heavy ion compositional signature reported for the Jovian events. This compositional feature has suggested that ions escaping from the Jovian magnetosphere play an important role in forming upstream ion populations at Jupiter. In contrast, models of energetic upstream ions at earth emphasize $\text{in situ}$ acceleration of reflected solar wind ions within the upstream region itself. Using Voyager 1 and 2 energetic (? 30 keV) ion measurements near the magnetopause, in the magnetosheath, and immediately upstream of the bow shock, we examine the compositional patterns together with typical energy spectra in each of these regions. We find characteristic spectral changes late in ion events observed upstream of the bow shock at the same time that heavy ion fluxes are enhanced and energetic electrons are present. A model involving upstream Fermi acceleration early in events and emphasizing energetic particle escape in the prenoon part of the Jovian magnetosphere late in events is presented to explain many of the features in the upstream region of Jupiter.     

行星际起伏向磁层顶的输运

时间尺度为分钟数量级的太阳风速度和行星际磁场大幅度扰动实际上始终存在于行星际空间的。这些扰动一直传输到紧贴磁层边界面外侧的区域。它们在磁鞘等离子体和磁层顶的相互作用过程中可能起很重要的作用。行星际起伏中的磁场分量在通过地球弓激波时首先经历一次跳跃,然后一部分扰动被带到磁层边界面处。在边界面附近磁场扰动幅度被大大地放大了。弓激波上游的太阳风条件控制了放大因子。本文所作的数值模拟研究结果表明,如果上游有大幅度的扰动,在边界面附近就有大幅度的Alfven起伏的磁场分量。当上游磁场接近垂直于日地联线时,放大因子变得相当大,而且放大因子随上游的等离子体β值和/或Alfven马赫数的增加而增加。上游各向异性对放大因子的影响不大。在磁层边界附近存在大幅度起伏表明这里不存在稳定的片流。      

Large geomagnetic storms of extreme solar event periods in solar cycle 23

During extreme solar events such as big flares or/and energetic coronal mass ejections (CMEs) high energy particles are accelerated by the shocks formed in front of fast interplanetary coronal mass ejections (ICMEs). The ICMEs (and their sheaths) also give rise to large geomagnetic storms which have significant effects on the Earth’s environment and human life. Around 14 solar cosmic ray ground level enhancement (GLE) events in solar cycle 23 we examined the cosmic ray variation, solar wind speed, ions density, interplanetary magnetic field, and geomagnetic disturbance storm time index (D_st). We found that all but one of GLEs are always followed by a geomagnetic storm with D_st  −50 nT within 1–5 days later. Most(10/14) geomagnetic storms have D_st index  −100  nT therefore generally belong to strong geomagnetic storms. This suggests that GLE event prediction of geomagnetic storms is 93% for moderate storms and 71% for large storms when geomagnetic storms preceded by GLEs. All D_st depressions are associated with cosmic ray decreases which occur nearly simultaneously with geomagnetic storms. We also investigated the interplanetary plasma features. Most geomagnetic storm correspond significant periods of southward B_z and in close to 80% of the cases that the B_z was first northward then turning southward after storm sudden commencement (SSC). Plasma flow speed, ion number density and interplanetary plasma temperature near 1 AU also have a peak at interplanetary shock arrival. Solar cause and energetic particle signatures of large geomagnetic storms and a possible prediction scheme are discussed.     

Cosmic ray anisotropies in the outer heliosphere

The observation of the directional distribution of energetic and cosmic ray particles has been done with the Voyager spacecraft over a long period. Since 2002, when the first flux enhancements of charged particles associated with the approach of Voyager 1 to the solar wind termination shock were observed, these anisotropy measurements have become of special interest. They play an important role to understand the magnetic field and shock structure and the basics of the modulation of cosmic ray and anomalous particles at and beyond the termination shock. They also serve as motivation to study the spatial behavior of galactic and anomalous cosmic ray anisotropies with numerical modulation models in order to illustrate how the radial anisotropy, at different energies, change from upstream to downstream of the termination shock. Observations made by Voyager 1 indicate that the termination shock is a complicated region than previously thought, hence the effects of the latitude dependence of the termination shock’s compression ratio and injection efficiency on the radial anisotropies of galactic and anomalous protons will be illustrated. We find that the magnitude and direction of the radial anisotropy strongly depends on the position in the heliosphere and the energy of particles. The effect of the TS on the radial anisotropy is to abruptly increase its value in the heliosheath especially in the A > 0 cycle for galactic protons and in both polarity cycles for anomalous protons. Furthermore, the global effect of the latitude dependence of the shock’s compression ratio is to increase the radial anisotropy for galactic protons throughout the heliosphere, while when combined with the latitude dependence of the injection efficiency this increase depends on modulation factors for anomalous protons and can even alter the direction of the radial anisotropy.     

Rapid and Asymmetric Response of the Earth's Bow Shock:Multipoint Observations

We analyze observations of three bow shock crossings which occurred during 2007, using upstream data from STEREO A/B, ACE and WIND, combined with multi-point THEMIS and Cluster data, and TC-1 data located near noon. During the crossing of 7 May 2007, we find that following a rapid reduction in solar wind ram pressure and subsequent pressure pulse seen by ACE and WIND upstream, the bow shock responds asymmetrically from dawn to dusk. Cluster data on the dawn-side suggest the bow shock is significantly flared and responds rapidly to the pulse arrival, while TC-1 at noon, and THEMIS on the dusk-side, are well matched to the model bow shock, but show a delayed response. The crossings observed on 21 May and 2 June show contrasting response matching the model boundary for northward Interplanetary Magnetic Field (IMF). The IMF and solar wind plasma data suggest that the bow shock crossing at dawn-dusk side and subsolar point were mainly caused by large and smaller scale features of the solar wind ram pressure rise rather than the influence of IMF.      

Comparison of heliospheric conditions near the earth during four recent solar maxima

Statistical properties of the daily averaged values of the solar activity (sunspot numbers, total solar irradiance and 10.7 cm radio emission indices), the solar wind plasma and the interplanetary magnetic field parameters near the Earth’s orbit are investigated for a period from 1964 to 2002 covering the maxima of four solar cycles from 20th to 23rd. Running half-year averages show significant solar cycle variations in the solar activity indices but only marginal and insignificant changes in comparison with background fluctuations for heliospheric bulk plasma and magnetic field parameters. The current 23rd cycle maximum is weaker than 21st and 22nd maxima, but slightly stronger than 20th cycle in most of solar and heliospheric manifestations.     

膨胀磁场和太阳风粒子相互作用三维数值模拟

采用三维模型,使用混合网格质点法HPIC(Hybrid Particle-in-Cell)对膨胀的磁场和太阳风相互作用过程进行数值模拟.研究了线圈产生的偶极子磁场在注入等离子体后和太阳风粒子的相互作用过程,并对以不同速度入射的等离子体引起的太阳风粒子的变化和磁场变化进行了比较.研究结果表明,偶极子磁场和太阳风作用时会产生弓形激波,此时磁压等于太阳风粒子的动压,当向线圈产生的偶极子磁场中注入高能等离子体时引起磁场膨胀,膨胀的磁场将会排斥太阳风粒子向外运动,从而引起弓形激波的变化,增大与太阳风相互作用的面积,并且粒子入射速度越大,磁场膨胀越明显,与太阳风相互作用愈强.      

An empirical model of the plasma environment around Mercury

Since the flyby observations by Mariner 10 in 1974 and 1975, Mercury has been one of the most interesting objects for space physics and planetary exploration. The MESSENGER and BepiColombo missions now plan to revisit this planet. In order to design plasma instruments for the BepiColombo mission, we have estimated electron and ion fluxes around Mercury with an empirical model, which has been developed for the Earth’s magnetotail. The solar wind data needed as input parameters are derived from Helios observations. The result shows that our predicted electron fluxes at aphelion agree well with the Mariner-10 data. It is also noted that ion instruments must cover a very wide dynamic range of proton fluxes. However, the applicability of the Earth’s magnetospheric model to Mercury is, in itself, an important issue for comparative magnetospheric studies.     

Global Hybrid Simulation of Magnetic Reconnection in the Magnetosheath

A three-dimensional (3-D) global hybrid simulation is carried out for the generation and structure of magnetic reconnection in the magnetosheath due to interaction of an interplanetary Tangential Discontinuity (TD) with the bow shock and magnetosphere. Runs are performed for solar wind TDs possessing different initial half-widths. As the TD propagates through the bow shock toward the magnetopause, it is greatly narrowed by a two-step compression processes, a "shock compression' followed by a subsequent ``convective compression'. In cases with a relatively thin solar wind TD, 3-D patchy reconnection is initiated in the transmitted TD, forming magnetosheath flux ropes. Multiple components of ion particles are present in the velocity distribution in the magnetosheath merging, accompanied by ion heating. For cases with a relatively wide initial TD, a dominant single X-line appears in the subsolar magnetosheath after the transmitted TD is narrowed. A shock analysis is performed for the detailed structure of magnetic reconnection in the magnetosheath. Rotational Discontinuity (RD)/Time-Dependent Intermediate Shock (TDIS) are found to dominate the reconnection layer, which and some weak slow shocks are responsible for the ion heating and acceleration.      

Shock acceleration of energetic particles in corotating interaction regions

A simple shock model for the acceleration of energetic particles in corotating interaction regions (CIR) in the solar wind is presented. Particles are accelerated at the forward and reverse shocks which bound the CIR by being compressed between the shock fronts and magnetic irregularities upstream from the shocks, or by being compressed between upstream irregularities and those downstream from the shocks. Particles also suffer adiabatic deceleration in the expanding solar wind, an effect not included in previous shock models for acceleration in CIR's. The model is able to account for the observed exponential spectra at earth, the observed behavior of the spectra with radial distance, the observed radial gradients in the intensity, and the observed differences in the intensity and spectra at the forward and reverse shocks.     

Energetic positive ions in the cometary “foreshock” region

Ions produced by ionization of the cometary neutrals interact with the solar wind protons to produce large amplitude oscillations of the ambient magnetic field. Such oscillations are convected towards the comet at the unperturbed solar wind speed far from the shock and at a lower speed closer to the shock (due to the solar wind mass loading); hence, they can energize the incoming ions by Fermi acceleration. The spatial extension of the acceleration region is of the order of 10⁶ km and the resulting energy spectrum is harder than in the Earth's bow shock case. The energization of cometary ions produces an additional deceleration of the solar wind. It is suggested that Comet Halley may be the most efficient “cosmic ray shock” in the solar system.     

Small-scale deformation of the bow shock

The prediction of the bow shock location is a proof of our understanding of the processes governing the solar wind – magnetosphere interaction. However, the models describing the bow shock location as a function of upstream parameters are based on a statistical processing of bow shock crossings observed by a single spacecraft. Such crossings locate the bow shock in motion, i.e., in a non-equilibrium state and this fact can be a source of significant errors. We have carefully analyzed a long interval of simultaneous observations of the bow shock and magnetopause and another interval of bow shock observations at two well-separated points. Our results suggest that often a small-scale deformation of the bow shock front due to magnetosheath fluctuations is the most appropriate interpretation of observations. Since the low-frequency magnetosheath variations exhibit largest amplitudes, a simultaneous bow shock displacement over a distance of 10–15 R_E can be observed. We suggest that bow shock models can be probably improved if the tilt angle would be implemented as a parameter influencing the bow shock location in high latitudes.     

Probing heliospheric diffusion coefficients with solar energetic particles

The dynamics of solar particle events provide a direct link to the understanding of properties of wave–particle interactions, and to the nature of the solar wind fluctuations. Depending on their energy, the often simultaneously observed electrons, protons and ions interact with different wavenumber ranges of the fluctuations, and are sensitive to various aspects of the dynamical nature of the solar wind turbulence. In general, the evolution of particle events is also sensitive to the spatial variation of the transport parameters between the Sun and a few AU. Together with in situ plasma and magnetic field observations this information can be used to extrapolate the properties of transport parameters into the more distant Heliosphere. Recent developments in the theory of parallel transport of energetic particles, and examples for the modelling of solar particle events and the derivation of transport parameters are considered. A dynamical quasi-linear theory is presented which gives special emphasis to the geometry and dynamic nature of the fluctuations, and which is able to provide particle mean free paths solely from observed plasma parameters, in good agreement with those derived by the modelling. Possibilities to apply the above results to the study of other energetic particle processes in the Heliosphere are discussed.     

Shock acceleration of energetic particles in wave heated corona

The solar wind wave heating models require substantial amount of wave power in order to efficiently heat and accelerate solar wind. The level of fluctuations is however limited by energetic particle observations. The simplest cyclotron sweep models result in convection-dominated transport, contradicting observations. However, models incorporating wave-wave -interactions, which cause wave energy to cascade in wavenumber, allow more reasonable energetic particle transport in the interplanetary space. The mean free path of the energetic particles remains still relatively short in the corona, providing favorable conditions for coronal mass ejection (CME) related shock acceleration. We study the consequences of this scenario on the energetic particle production related to CMEs. The role of self-generated waves is also discussed.     

Magnetic Field Fluctuations in the Solar Wind, Foreshock and Magnetosheath： Cluster Data Analysis

46 magnetosheath crossing events from the two years (2001.2-2003.1) of Cluster magnetic field measurements are identified and used to investigate the characters of the magnetic field fluctuations in the regions of undisturbed solar wind, foreshock, magnetosheath. The preliminary results indicate the properties of the plasma turbulence in the magnetosheath are strongly controlled by IMF orientation with respect to the bow shock normal. The amplitude of the magnetic field magnitude and direction variations behind quasi-parallel bow shock are larger than those behind quasi-perpendicular bow shock. Almost purely compressional waves are found in quasi-perpendicular magnetosheath.