Interplanetary Physics in Mainland China

During the past two years (2016-2018), great achievements have been made in the Chinese research of interplanetary physics, with nearly 100 papers published in the academic journals. The achievements are including but not limited to the following topics:solar corona; solar wind and turbulence; filament/prominence and jets; solar flare; radio bursts; particle acceleration at coronal shocks; magnetic flux ropes; instability; instrument; Coronal Mass Ejections (CMEs) and their interplanetary counterparts; Magnetohydrodynamic (MHD) numerical modeling; solar energetic particles and cosmic rays. The progress further improves our understanding of the eruptions of solar activities, their evolutions and propagations in the heliosphere, and final geoeffects on our Earth. These results were achieved by the Chinese solar and space scientists independently or via international collaborations. This paper will give a brief review of these achievements.      

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.     

Evolution of Dst index,cosmic rays and global temperature during solar cycles 20–23

We have studied conditions in interplanetary space, which can have an influence on galactic cosmic ray (CR) and climate change. In this connection the solar wind and interplanetary magnetic field parameters and cosmic ray variations have been compared with geomagnetic activity represented by the equatorial Dst index from the beginning 1965 to the end of 2012. Dst index is commonly used as the solar wind–magnetosphere–ionosphere interaction characteristic. The important drivers in interplanetary medium which have effect on cosmic rays as CMEs (coronal mass ejections) and CIRs (corotating interaction regions) undergo very strong changes during their propagation to the Earth. Because of this CMEs, coronal holes and the solar spot numbers (SSN) do not adequately reflect peculiarities concerned with the solar wind arrival to 1 AU. Therefore, the geomagnetic indices have some inestimable advantage as continuous series other the irregular solar wind measurements. We have compared the yearly average variations of Dst index and the solar wind parameters with cosmic ray data from Moscow, Climax, and Haleakala neutron monitors during the solar cycles 20–23. The descending phases of these solar cycles (CSs) had the long-lasting solar wind high speed streams occurred frequently and were the primary contributors to the recurrent Dst variations. They also had effects on cosmic rays variations. We show that long-term Dst variations in these solar cycles were correlated with the cosmic ray count rate and can be used for study of CR variations. Global temperature variations in connection with evolution of Dst index and CR variations is discussed.     

New insights from cross-correlation studies between solar activity indices and cosmic-ray flux during Forbush decrease events

Observed galactic cosmic ray intensity can be subjected to a transient decrease. These so-called Forbush decreases are driven by coronal mass ejection induced shockwaves in the heliosphere. By combining in situ measurements by space borne instruments with ground-based cosmic ray observations, we investigate the relationship between solar energetic particle flux, various solar activity indices, and intensity measurements of cosmic rays during such an event. We present cross-correlation study done using proton flux data from the SOHO/ERNE instrument, as well as data collected during some of the strongest Forbush decreases over the last two completed solar cycles by the network of neutron monitor detectors and different solar observatories. We have demonstrated connection between the shape of solar energetic particles fluence spectra and selected coronal mass ejection and Forbush decrease parameters, indicating that power exponents used to model these fluence spectra could be valuable new parameters in similar analysis of mentioned phenomena. They appear to be better predictor variables of Forbush decrease magnitude in interplanetary magnetic field than coronal mass ejection velocities.     

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.     

Some characteristics of propagation of flare- or CME-associated interplanetary shock waves

Many interplanetary shock waves have a fast mode MHD wave Mach number between one and two and the ambient solar wind plasma and magnetic field are known to fluctuate. Therefore a weak, fast, MHD interplanetary shock wave propagating into a fluctuating solar wind region or into a solar wind stream will be expected to vary its strength.It is possible that an interplanetary shock wave, upon entering such a region will weaken its strength and degenerate into a fast-mode MHD wave. It is even possible that the shock may dissipate and disappear.A model for the propagation of a solar flare - or CME (Coronal Mass Ejections) - associated interplanetary shock wave is given. A physical mechanism is described to calculate the probability that a weak shock which enters a turbulent solar wind region will degenerate into a MHD wave. That is, the shock would disappear as an entropy-generate entity. This model also suggests that most interplanetary shock waves cannot propagate continuously with a smooth shock surface. It is suggested that the surface of an interplanetary shock will be highly distorted and that parts of the shock surface can degenerate into MHD waves or even disappear during its global propagation through interplanetary space. A few observations to support this model will be briefly described.Finally, this model of shock propagation also applies to corotating shocks. As corotating shocks propagate into fluctuating ambient solar wind regions, shocks may degenerate into waves or disappear.     

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.     

Prediction of expected global climate change by forecasting of galactic cosmic ray intensity time variation in near future based on solar magnetic field data

A method of prediction of expected part of global climate change caused by cosmic ray (CR) by forecasting of galactic cosmic ray intensity time variation in near future based on solar activity data prediction and determined parameters of convection-diffusion and drift mechanisms is presented. This gave possibility to make prediction of expected part of global climate change, caused by long-term cosmic ray intensity variation. In this paper, we use the model of cosmic ray modulation in the Heliosphere, which considers a relation between long-term cosmic ray variations with parameters of the solar magnetic field. The later now can be predicted with good accuracy. By using this prediction, the expected cosmic ray variations in the near Earth space also can be estimated with a good accuracy. It is shown that there are two possibilities: (1) to predict cosmic ray intensity for 1–6 months by using a delay of long-term cosmic ray variations relatively to effects of the solar activity and (2) to predict cosmic ray intensity for the next solar cycle. For the second case, the prediction of the global solar magnetic field characteristics is crucial. For both cases, reliable long-term cosmic ray and solar activity data as well as solar magnetic field are necessary. For solar magnetic field, we used results of two magnetographs (from Stanford and Kitt Peak Observatories). The obtained forecasting of long-term cosmic ray intensity variation we use for estimation of the part of global climate change caused by cosmic ray intensity changing (influenced on global cloudiness covering).     

Evidence for mass-loading of the Venus magnetosheath

The observed magnetic field configuration in the Venus magnetosheath contains information about the solar wind mass-loading processes occurring as a result of the extension of the neutral atmosphere into the magnetosheath. In this paper, magnetic field signatures of various mass-loading processes are discussed and experimental results from the Pioneer Venus Orbiter magnetometer experiment are examined for evidence of these signatures. The data suggest that the ?$\text{V}$X$\text{B}$ acceleration process, stochastic pickup of ionospheric ions, and $\text{J}$X$\text{B}$ force “scavenging” at the ionopause all occur at various times.     

Space storm measurements of the July 2005 solar extreme events from the low corona to the Earth

The Athens Neutron Monitor Data Processing (ANMODAP) Center recorded an unusual Forbush decrease with a sharp enhancement of cosmic ray intensity right after the main phase of the Forbush decrease on 16 July 2005, followed by a second decrease within less than 12 h. This exceptional event is neither a ground level enhancement nor a geomagnetic effect in cosmic rays. It rather appears as the effect of a special structure of interplanetary disturbances originating from a group of coronal mass ejections (CMEs) in the 13–14 July 2005 period. The initiation of the CMEs was accompanied by type IV radio bursts and intense solar flares (SFs) on the west solar limb (AR 786); this group of energetic phenomena appears under the label of Solar Extreme Events of July 2005. We study the characteristics of these events using combined data from Earth (the ARTEMIS IV radioheliograph, the Athens Neutron Monitor (ANMODAP)), space (WIND/WAVES) and data archives. We propose an interpretation of the unusual Forbush profile in terms of a magnetic structure and a succession of interplanetary shocks interacting with the magnetosphere.     

Interaction of shocks with a current sheet and reconnection in the solar corona

Interaction of shocks with a current sheet is investigated within a 2D MHD model based on an improved FCT numerical scheme. Basic parameters of the problem are chosen to correspond to situations in the solar corona with low plasma β and moderate shock strength. Slow and fast MHD shocks are introduced with shock normal parallel to magnetic field lines. The interaction with the current sheet causes distortion of the shock front and this distorts the magnetic field lines and generates electric current. Large current densities are generated especially when the fast MHD shock becomes the intermediate MHD shock at the current sheet. Then peak values of the current density are about 3–4 times larger than the initial undisturbed values in the current sheet.     

Cosmic rays in the dynamic heliosheath

Cosmic ray modulation in the outer heliosphere is discussed from a modeling perspective. Emphasis is on the transport and acceleration of these particles at and beyond the solar wind termination shock in the inner heliosheath region and how this changes over a solar cycle. We will show that by using numerical models, and by comparing results to spacecraft observations, much can be learned about the dependence of cosmic ray modulation on solar cycle changes in the solar wind and heliospheric magnetic field. While the first determines the heliospheric geometry and shock structure, the latter results in a time-dependence of the transport coefficients. Depending on energy, both these effects contribute to cosmic ray intensities in the inner heliosheath changing over a solar cycle.     

强磁场扰动对宇宙线调制的统计研究

对1978─1982太阳活动高年时发生的激波、强磁场扰动及激波与强磁场扰动共存这三类事件引起的宇宙线变化进行了统计研究,得到如下结果:(1)激波与强磁场扰动共存时引起的宇宙线强度下降最为显着;只有激波或强磁场扰动时,宇宙线的强度变化相对较小;(2)标志速度间断的激波是产生宇宙线Forbush下降的重要因素;(3)速度间断在强磁场扰动对宇宙线的调制中可能起一个触发的作用。      

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.     

The local time dependence of the anisotropic solar cosmic ray flux.

The distribution of the solar cosmic radiation flux over the earth is not uniform, but the result of complex phenomena involving the interplanetary magnetic field, the geomagnetic field and latitude and longitude of locations on the earth. The latitude effect relates to the geomagnetic shield; the longitude effect relates to local time. For anisotropic solar cosmic ray events the maximum particle flux is always along the interplanetary magnetic field direction, sometimes called the Archimedean spiral path from the sun to the earth. During anisotropic solar cosmic ray event, the locations on the earth viewing "sunward" into the interplanetary magnetic field direction will observe the largest flux (when adjustments are made for the magnetic latitude effect). To relate this phenomena to aircraft routes, for anisotropic solar cosmic ray events that occur during "normal quiescent" conditions, the maximum solar cosmic ray flux (and corresponding solar particle radiation dose) will be observed in the dawn quadrant, ideally at about 06 hours local time.     

Relativistic solar protons in the ground level event of 23 February 1956: New study

In this study we applied again to the outstanding solar particle event of 23 February 1956, the largest one in the entire history of observations of solar cosmic rays. Due to significant improvement of the analysis/modeling techniques and new understanding of physical processes in the solar atmosphere and interplanetary space, a possibility arises to interpret the old data in the light of modern concept of multiple particle acceleration at/near the Sun. In our new analysis the data of available then neutron monitors and muon telescopes are used. The technique of the analysis includes: (a) calculation of asymptotic cones of ground-based detectors; (b) modeling of cosmic ray detector responses at variable parameters of the flux of solar relativistic protons; (c) determination of primary solar proton parameters outside magnetosphere by comparison of computed responses with observations. Certain evidence was obtained that the flux of relativistic solar protons consisted of two distinct components: prompt and delayed ones. The prompt component with exponential energy spectrum caused a giant impulse-like increase at a number of European cosmic ray stations. The delayed component had a power-law spectrum and was a cause of gradual increase at cosmic ray stations in the North American region. A numerical simulation of the proton acceleration in the vicinity of the magnetic reconnection region brings to the proton spectrum with exponential dependence on energy. This agrees with observational data for the prompt component. It is also shown that the huge increase in ∼5000% on neutron monitors was due to the prompt component only with the exponential proton spectrum. The power-law spectrum of comparable intensity gave considerably smaller effect.     

Observed solar cycle modulation of galactic cosmic rays over a range of rigidities

We have studied the long-term, steady-state, solar cycle modulation of galactic cosmic ray intensity for seven cycles (17–23). Our analysis is based on the data obtained with a variety of detectors on earth (neutron monitors of the global network and muon detectors) as well as telescopes flown on high altitude balloons and on-board near-earth satellites. The median rigidity of response for these detectors to galactic cosmic ray spectrum lies in the range 1–70 GV. We correlate cosmic ray data to sunspot numbers, Ap, solar wind bulk speed (V), magnetic field (B), as well as to the cycle maximum (M), minimum (m), and the epochs of the solar polar field reversals. This enables us to derive the rigidity dependence of observations, and helps us to define the characteristics of the modulation function in the heliosphere.     

50–200 GV cosmic rays at various heliolatitudes

We consider regular motion of 50 – 200 GV particles in a large-scale interplanetary magnetic field model which contains a wavy neutral sheet responsible for the sector-structure. Numerical calculations based upon energy losses along various trajectories are carried out to obtain the predicted omni-directional density and anisotropy of cosmic rays at various solar latitudes. A marked difference is found between odd and even solar cycles. The post-1969 field configuration gives small radial and large latitudinal gradient: cosmic ray density increases toward the poles. The latitudinal gradient turns out smaller and of opposite sense for the pre-1969 epoch. Anisotropy changes dramatically as we move off the solar equator: corotation appears to be restricted to low latitudes.     

Cone model-based SEP event calculations for applications to multipoint observations

The problem of modeling solar energetic particle (SEP) events is important to both space weather research and forecasting, and yet it has seen relatively little progress. Most important SEP events are associated with coronal mass ejections (CMEs) that drive coronal and interplanetary shocks. These shocks can continuously produce accelerated particles from the ambient medium to well beyond 1 AU. This paper describes an effort to model real SEP events using a Center for Integrated Space weather Modeling (CISM) MHD solar wind simulation including a cone model of CMEs to initiate the related shocks. In addition to providing observation-inspired shock geometry and characteristics, this MHD simulation describes the time-dependent observer field line connections to the shock source. As a first approximation, we assume a shock jump-parameterized source strength and spectrum, and that scatter-free transport occurs outside of the shock source, thus emphasizing the role the shock evolution plays in determining the modeled SEP event profile. Three halo CME events on May 12, 1997, November 4, 1997 and December 13, 2006 are used to test the modeling approach. While challenges arise in the identification and characterization of the shocks in the MHD model results, this approach illustrates the importance to SEP event modeling of globally simulating the underlying heliospheric event. The results also suggest the potential utility of such a model for forcasting and for interpretation of separated multipoint measurements such as those expected from the STEREO mission.