准平行无碰撞激波上游能量粒子的观测研究

通过Cluster卫星在2005年3月16日观测到的一个准平行激波观测事例,研究了准平行激波上游低频等离子体波动与能量离子之间的关系.卫星观测结果表明,在准平行激波上游,离子微分能通量受到了非线性波动的调制.在磁场强度较小区域,离子微分能通量较高.产生这种现象的可能原因是准平行激波上游的非线性波动可以捕获离子,被捕获的离子在波动中来回弹跳并被电场加速,从而导致磁场强度较小区域离子微分能通量较高.这一观测结果与已有的混合模拟结果相吻合.      

Positive ion distributions in the morning auroral zone: Local acceleration and drift effects

We report on the typical structure of the large scale ion precipitation in the morning sector of the auroral zone and associated low frequency electromagnetic waves. Data obtained during near radial passes of the AUREOL-3 satellite point to a distinction between two main precipitation regions: 1) In the poleward part of the auroral zone the latitudinal variation of the average energy (or temperature) of the precipitated ions (mainly H⁺) indicate that they are adiabatically accelerated in the outer magnetosphere. This “high energy” (? 3 to > 20 keV) precipitation is usually associated with a low energy (E < 110 eV) upward flowing 0⁺ and H⁺ component, and 2) near the boundary between discrete and diffuse electron aurorae a drastic change in the ion characteristics is observed. The flux of energetic precipitated H⁺ ions is sharply reduced, which suggests the formation of an Alfvén layer. However, intense fluxes of precipitated H⁺, O⁺, and He⁺ ions with energies < 3 keV are observed equatorward of the Alfvén layer, in coincidence with the diffuse aurora and in association with quasi-monochromatic electromagnetic waves with frequencies around the proton gyrofrequency. As the characteristic convection and bounce times of the low energy upward flowing ion component are comparable (τ > 3 hours) we suggest that the precipitation of ionospheric ions inside the diffuse aurora results from convection and corotation of the ions accelerated to suprathermal energies at higher latitudes.     

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.     

Cometary kilometric radio waves and plasma waves correlated with ion pick-up effect at comet Halley

From the discrete spectra of the emissions from the comet in the frequency range from 30 to 195 kHz named CKR (Cometary Kilometric Radiation), movements of the bow shock at comet Halley are concluded, i.e., the observed CKR emissions can be interpreted as being generated and propagating from the moving shock. The motion of the shocks are possibly associated with time variation of the solar wind and of the cometary outgassings. By in-situ plasma waves observations using PWP (Plasma Wave Probe) onboard the Sakigake spacecraft, the characteristic spectra of the electrostatic electron plasma waves, the electron cyclotron harmonic waves, and the ion sound waves have been detected during the interval of the Halley's comet fly-by. Compared with the results of a Faraday cup observation and a magnetometer, it is concluded that these plasma wave phenomena are the manifestation of the ion pick-up processes. The ion pick-up processes are taking place even in the remote region within a distance range from 7×10⁶ to 10⁷ km from the cometary nucleus.     

太阳风中磁场重联离子扩散区的大尺度特征

利用ACE和WIND卫星2007年1月6日的联合探测, 在1AU附近发现了一个等离子体密度极低的Petschek-like重联喷流区. 该喷流区内部出现了非常明显的Hall双极磁场、等离子体密度下降区以及与Hall电流相符的低能段电子投掷角分布. 这些特征与重联离子扩散区的Hall效应非常吻合, 说明很可能在太阳风中观测到了一个离子扩散区. 分析表明, 与之相关的磁场重联为准稳态快速完全反向重联, 其扩散区以一对慢模波为边界, 空间尺度达到80个离子惯性长度, 表现出了大尺度重联的特征.      

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.     

Geomagnetic activity associated with magnetic clouds,magnetic cloud-like structures and interplanetary shocks for the period 1995–2003

Using nine years (1995–2003) of solar wind plasma and magnetic field data, solar sunspot number, and geomagnetic activity data, we investigated the geomagnetic activity associated with magnetic clouds (MCs), magnetic cloud-like structures (MCLs), and interplanetary shock waves. Eighty-two MCs and one hundred and twenty-two MCLs were identified by using solar wind and magnetic field data from the WIND mission, and two hundred and sixty-one interplanetary shocks were identified over the period of 1995–2003 in the vicinity of Earth. It is found that MCs are typically more geoeffective than MCLs or interplanetary shocks. The occurrence frequency of MCs is not well correlated with sunspot number. By contrast, both occurrence frequency of MCLs and sudden storm commencements (SSCs) are well correlated with sunspot number.     

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.     

Observations of solar radio bursts from meter to kilometer wavelengths

Using the Clark Lake Radioheliograph data we present direct evidence that type III electron streams propagate in dense coronal streamers. We also present imaging observations of meter-decameter microbursts, which appear to be similar to those observed in hard X-rays. At meter-decameter wavelengths, these microbursts appear to be due to plasma radiation. From observations made with ISSE-3, we discuss the characteristics of hectometer and kilometer wavelength radio bursts. In particular, we show that from studies of type III storms that the exciter electrons propagate along spiral structures, where the density is enhanced and that there is an acceleration of the solar wind. We discuss type II bursts at kilometer wavelengths, compare them with meter type II bursts and discuss their association with interplanetary shocks. We show that the interaction between type III electron streams and shocks at kilometer wavelengths can provide information on the interplanetary shock geometry. Finally, we discuss the possibility that some shock associated (SA) events may be emissions caused by electrons accelerated lower in the atmosphere rather than high in the corona in type II shocks.Recent advances in solar research have resulted from new work on plasma radiation theory, new observations of active regions and flares across the electromagnetic spectrum and the availability of spacecraft in situ measurements of solar ejecta. In this paper, we review some results obtained with the Clark Lake multifrequency radioheliograph at meter-decameter wavelengths and from satellite multifrequency directive observations at hectometer and kilometer wavelengths. We present evidence that type III electrons propagate in dense coronal streamers, and that frequently observed microbursts (presumably of type III) at meter-decameter wavelengths are due to plasma radiation. We discuss observations of hectometer and kilometer type III radio storms which reveal information about active region structures, interplanetary magnetic field configuration, and solar wind acceleration. We also discuss kilometer type II bursts, interactions between type III electrons and interplanetary shocks, and present some new results on shock associated (SA) events.     

Ion relaxation processes in the heliospheric interface: how perturbed are ion distribution functions?

Charge-exchange processes between interstellar H-/O-atoms and protons of the bulk of the interstellar plasma flow downstream of the outer bowshock in the heliospheric interface induce secondary ions leading to non-relaxated velocity distribution functions. The relaxation of these freshly induced ions towards an equilibrium distribution occurs due to Coulomb interactions and wave–particle interactions with the background turbulence. Since Coulomb interactions are of low relevance, we study here in detail the effect of wave–particle interactions. To find the turbulence levels in the interface we consider the MHD-wave transformation at the outer shock surface between the interface and the local interstellar plasma. The turbulence in the outer interface region is shown to be dominated by incompressible Alfvén waves both for cases of quasiparallel and quasiperpendicular shocks. Also we show that waves propagating towards the shock are more intensive than those propagating away from it. The level of Alfvén turbulence in the interface is estimated using the recent data on local interstellar turbulence deduced from observations of interstellar scintillations of distant radiosources. Two proton relaxation processes are considered: quasilinear resonant interactions with Alfvén waves and non-linear self-induced wave–particle scattering. The corresponding diffusion coefficients are estimated, and typical time periods for protons and oxygen ions relaxation are shown to be of the same order of magnitude as H-/O-atoms passage time over the extent of the interface. This indicates that perturbed ion distribution functions must be expected there.     

Excitation of flare-induced waves in coronal loops and the effects of radiative cooling

EUV imaging observations from several space missions (SOHO/EIT, TRACE, and SDO/AIA) have revealed a presence of propagating intensity disturbances in solar coronal loops. These disturbances are typically interpreted as slow magnetoacoustic waves. However, recent spectroscopic observations with Hinode/EIS of active region loops revealed that the propagating intensity disturbances are associated with intermittent plasma upflows (or jets) at the footpoints which are presumably generated by magnetic reconnection. For this reason, whether these disturbances are waves or periodic flows is still being studied. This study is aimed at understanding the physical properties of observed disturbances by investigating the excitation of waves by hot plasma injections from below and the evolution of flows and wave propagation along the loop. We expand our previous studies based on isothermal 3D MHD models of an active region to a more realistic model that includes full energy equation accounting for the effects of radiative losses. Computations are initialized with an equilibrium state of a model active region using potential (dipole) magnetic field, gravitationally stratified density and temperature obtained from the polytropic equation of state. We model an impulsive injection of hot plasma into the steady plasma outflow along the loops of different temperatures, warm (～1 MK) and hot (～6 MK). The simulations show that hot jets launched at the coronal base excite slow magnetoacoustic waves that propagate to high altitudes along the loops, while the injected hot flows decelerate rapidly with heights. Our results support that propagating disturbances observed in EUV are mainly the wave features. We also find that the effect of radiative cooling on the damping of slow-mode waves in 1–6 MK coronal loops is small, in agreement with the previous conclusion based on 1D MHD models.     

Dispersive shock waves in partially ionised plasmas

Compressional waves propagating in the partially ionised solar lower atmospheric plasmas can easily steepen into nonlinear waves, including shocks. Here we investigate the effect of weak dispersion generated by Hall currents perpendicular to the ambient magnetic field on the characteristics of shock waves. Our study will also focus on the interplay between weak dispersion and partial ionisation of the plasma. Using a multiple scale technique we derive the governing equation in the form of a Korteweg-de Vries-Burgers equation. The effect of weak dispersion on shock waves is obtained using a perturbation technique. The secular behaviour of second order terms is addressed with the help of a renormalization technique. Our results show that dispersion modifies the characteristics of shock waves and this change is dependent also on the ionisation degree of the plasma. Dispersion can create short lived oscillations in the shocked plasma. The shock fronts become wider with the increase in the number of neutrals in the plasma.     

多重激波的数值模拟研究

本文用一维混合粒子模拟Code研究了包括中间激波在内的多重激波.模拟了四种情形,可以分为两类:(1)由快激波和中间激波构成的两重激波,(2)快激波、中间激波和慢激波构成的三重激波.结果表明:多重激波是不稳定的,它趋向于发展成磁流体旋转间断和MHD波,左旋圆偏振波逐渐在上游区内发展起来.文章对导致多重激波不稳定性的可能原因进行了简单的讨论.     

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.     

Nonplanar electron acoustic shock waves

The properties of cylindrical and spherical electron acoustic shock waves (EASWs) in an unmagnetized plasma consisting of cold electrons, immobile ions and Boltzmann distributed hot electrons are investigated by employing the reductive perturbation method. A Korteweg–de Vries Burgers (KdVB) equation is derived and its numerical solution is obtained. The effects of several parameters and ion kinematic viscosity on the basic features of EA shock waves are discussed in nonplanar geometry. It is found that nonplanar EA shock waves behave quite differently from their one-dimensional planar counterpart.     

Forward and reverse CIR shocks at 4–5 AU: Ulysses

In this article, we study fast shocks at CIR boundaries during an extended interval of 15 consecutive major high speed solar wind streams in 1992–1993. Ulysses was 4–5 AU from the sun. The Abraham-Schrauner shock normal method and the Rankine-Hugoniot relations were used to determine fast shock directions and speeds. Out of 33 potential CIR shocks, 14 were determined to be fast forward shocks (FSs) and 14 were fast reverse shocks (RSs). Of the remaining 5 events, 2 were forward waves and 3 were reverse waves. CIR edges at latitudes below ∼30^o were, for the most part, bounded by fast magnetosonic shocks. The forward shocks were generally quasi-perpendicular (average θ_nBo = 67^o). The reverse shocks were more oblique (average θ_nBo = 52^o), but they extended to all angles. Both FSs and RSs had magnetosonic Mach numbers ranging from 1 to 5 or 6. The average Mach numbers were 2.4 and 2.6 for FSs and RSs, respectively. The shock Mach numbers were noted to generally decrease with increasing latitude. The non-shock events or waves were noted to occur preferentially at high (∼−30° to −35°) heliolatitudes where stream-stream interactions were presumably weaker. These results are consistent with expectations, indicating the general accuracy of the Abraham-Schrauner technique.     

Acceleration and transport of energeticparticles at CME-driven shocks

Historically, solar energetic particle (SEP) events are classified in two classes as “impulsive” and “gradual”. Whether there is a clear distinction between the two classes is still a matter of debate, but it is now commonly accepted that in large “gradual” SEP events, Fermi acceleration, also known as diffusive shock acceleration, is the underlying acceleration mechanism. At shock waves driven by coronal mass ejections (CMEs), particles are accelerated diffusively at the shock and often reach > MeV energies (and perhaps up to GeV energies). As a CME-driven shock propagates, expands and weakens, the accelerated particles can escape ahead of the shock into the interplanetary medium. These escaping energized particles then propagate along the interplanetary magnetic field, experiencing only weak scattering from fluctuations in the interplanetary magnetic field (IMF). In this paper, we use a Monte-Carlo approach to study the transport of energetic particles escaping from a CME-driven shock. We present particle spectra observed at 1 AU. We also discuss the particle “crossing number” at 1AU and its implication to particle anisotropy. Based on previous models of particle acceleration at CME-driven shocks, our simulation allows us to investigate various characteristics of energetic particles arriving at various distances from the sun. This provides us an excellent basis for understanding the observations of high-energy particles made at 1 AU by ACE and WIND.     

Acceleration theory for 5–40 keV ions at interplanetary shocks

Power-law spectra f(E)∝E^?2.7 of < 40 keV suprathermal ions within ～10⁷ km of propagating interplanetary shocks are explained by diffusive scattering near a plane shock. The theory fits the 25 November 1977 event with a mean free path perpendicular to the shock with is 0.01 AU in front of the shock and less than .0003 AU behind it, for 1 keV ions. The theory predicts a steepening spectrum at higher energies, of the form $\text{f(v)∝v}{}^{\mathrm{?4}}\text{exp(??λdv/ur)}\text{where}\text{u = (ΔV)}{}^2\text{/2V}{}_{\mathrm{W}}$ depends on the plasma velocity jump ΔV and the plasma speed V_W and mean free path λ in front of the shock     

Magnetic reconnection coupled with heat conduction

Magnetic reconnection coupled with anisotropic nonlinear heat conduction is studied by using an MHD simulation. Due to the heat conduction effect, the adiabatic slow-mode MHD shocks emanating from the neutral point are dissociated into conduction fronts and isothermal shocks.     

Sources and losses of ring current ions: An update

The cleft ion fountain has been identified as a prodigious source of upflowing suprathermal ionospheric plasma. Modeling efforts have traced the path of these ions from the polar ionosphere along trajectories where the ions are energized to keV energies and deposited in the near earth plasma sheet. Mass and energy dispersion of these ions accounts in a natural way for the observed variation in heavy ion content of the plasma sheet. Observations of ion composition in the plasma sheet by the AMPTE and ISEE spacecraft establish that ionospheric ions dominate in the near earth plasma sheet but solar wind ions become significant tailward. The heavy ion content of the plasma sheet increases with both solar cycle and magnetic activity. Direct injection of ionospheric ions into the ring current has been observed in the outer plasmasphere. Several mechanisms for the direct injection of ions from the plasmasphere and ionosphere into the ring current have appeared. Estimation of ionospheric source strengths and residence times have led to an estimate of the magnetospheric densities that would result solely from an ionospheric outflow populating the magnetosphere. Estimated densities were quite reasonable even without inclusion of a solar wind source of ions. Ring current ions decay primarily via charge exchange with the hydrogen geocorona, however, the roles of pitch angle diffusion and Coulomb collisions in this decay process are being clarified.

Modeling and observations of ENA by the 1SEE1 spacecraft has led to a re-affirmation of the dominant role of charge exchange in ring current decay. Ion cyclotron waves contribute to ring current decay in the dusk bulge region. The role of low frequency. (< 1 Hz) ion cyclotron waves in the plasmasphere is still unclear. Other wave modes may be responsible for the pitch angle diffusion and subsequent loss of ring current ions. Coulomb collisional energy losses from ring current O⁺ to thermal electrons are sufficient to power SAR arcs and represent an energy sink for ring current O⁺ within the plasmasphere. Coulomb collisions may be important for decay of low energy (< 10 KeV) ring current ions in the plasmasphere.