ICME期间磁暴的太阳风参数和极光沉降能量

基于1995-2004年ICME驱动的强烈磁暴(SA型)、强磁暴(SB型)和延迟型主相暴(SC型)三种磁暴类型,对1AU处太阳风动压、太阳风速度、行星际磁场、EK-L电场以及极光沉降能量进行时序叠加分析,并分别与-vBz耦合函数和Newell耦合函数进行对比.结果表明,三种磁暴在ICME到达前期的太阳风动压较稳定,背景太阳风、极光沉降能量、行星际磁场和磁层存在相对平静期. ICME到达前期SA型磁暴的背景太阳风速度、行星际磁场南向分量以及极光沉降能量的均值高于另外两种磁暴类型,这说明大型日冕物质抛射在ICME到达前就对行星际磁场、背景太阳风和HP产生了影响.磁暴急始后,SC型磁暴的EK-L电场斜率小,峰值延后且行星际磁场北向分量增强,这些都是磁暴主相延迟的表现,极光沉降能量随着行星际磁场转为南向而增加.     

异常SSC事件的行星际原因分析

磁暴急始(SSC)是强烈太阳风动压或行星际激波与磁层相互作用的结果.通常SSC事件的上升时间在4～10 min,我们把上升时间超过15 min的SSC事件称为异常SSC事件.本文利用地磁SYM-H指数鉴别出了5个有地磁观测历史以来发生的上升时间大于15 min的异常SSC事件,并利用Wind,ACE,IMP 8,Goes,Geotail多点卫星太阳风观测数据和地磁观测数据,分析了异常SSC事件的行星际原因.结果表明,异常SSC事件通常都是强烈行星际扰动引起的,5个异常SSC事件有4个对应于行星际激波,有3个对应于多步太阳风动压跃变,有1个对应于行星际电场大幅度变化;由行星际激波产生的异常SSC事件,其上升时间依赖于行星际激波的方向,方向相对于日地连线越偏,上升时间越长;异常SSC事件上升时间与行星际磁场方向关系不明显.      

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

Modeling and data analysis of a Forbush decrease

We study the Forbush decrease of the galactic cosmic ray intensity observed in 9–25 September 2005 using the experimental data and a newly developed time-dependent three dimensional modeling. We analyze neutron monitors and muon telescopes, and the interplanetary magnetic field data. We demonstrate a clear relationship between the rigidity (R) spectrum exponent (γ) of the Forbush decrease and the exponent (ν) of the power spectral density of the components of the interplanetary magnetic field in the frequency range of ∼ 10⁻⁶–10 ⁻⁵ Hz. We confirm that an inclusion of the time-dependent changes of the exponent ν makes the newly developed nonstationary three dimensional model of the Forbush decrease compatible with the experimental data. Also, we show that the changes of the rigidity spectrum exponent γ does not depend on the level of convection of the galactic cosmic rays stream by solar wind; depending on the changes of the exponent ν, i.e. on the state of the turbulence of the interplanetary magnetic field.     

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.     

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.     

行星际闪烁功率谱的Fresnel现象与太阳风速度

讨论了弱散射条件下行星际闪烁(IPS)功率谱的Fresnel现象。数值计算结果表明,无行星际扰动存在时,IPS功率谱的膝点频率仍是Fresnel频率,在日张角。ε≤35°的范围内,Fresnel振荡现象仍然存在,极小值所在频率与太阳风电子密度为一维结构薄屏近似的解析结果相差不大。行星际扰动存在时,功率谱的膝点频率并未改变,扰动区太阳风速度决定着Fresnel振荡出现的频率。      

行星际激波传播至1AU的渡越时间预测

考虑了激波爆发源角宽度、能量、驱动时间、激波速度及其与背景太阳风之间的相互作用,利用流体力学扰动方程建立起一个激波扰动传播模型,用于研究激波从太阳传播到地球轨道附近(1 AU处)所需要的时间(渡越时间)问题.为印证扰动传播模型的适用性,利用1979-1989年间的27个激波事件,以及1997年2月到2000年1月间的68个激波事件,对激波到达地球轨道附近的渡越时间进行了预测,并将结果与STOA和ISPM预报模型结果进行了比较.实验表明,该模型在所有95个事件中,渡越时间相对误差小于10%的事件数占总事件数的25.26%;相对误差小于20%的占总事件数的50.53%;相对误差小于30%的占总事件的65.26%.      

Peculiarities of galactic cosmic ray Forbush effects during October–November 2003

Features of two successive Forbush effects of the galactic cosmic ray intensity in October–November 2003 have been studied based on the neutron monitors data. The rigidity spectrum of the galactic cosmic ray intensity in the course of the first Forbush effect (22–27 October) is gradually hardening, while the rigidity spectrum of the second Forbush effect (28 October–10 November) from the starting moment is very hard. As far, the energy range of the turbulence of the interplanetary magnetic field is in general responsible for the diffusion of galactic cosmic ray particles of the energy 5–50 GeV (to which neutron monitors are sensitive), we postulate that the gradually hardening (from day to day) of the rigidity spectrum of the first Forbush effect is associated with the enhancement of the power spectral density in the energy range of the interplanetary magnetic field turbulence caused by the large scale irregularities generated due to the interaction of the extending high speed disturbances with the background solar wind. The very hard rigidity spectrum (from the starting moment) of the second Forbush effect is generally associated with the well established new structure of the energy range of the interplanetary magnetic field turbulence enriched by the already created large scale irregularities. The gradually softening of the rigidity spectrum during the recovery phase of the second Forbush effect confirms that the disturbed interplanetary magnetic field turbulence step by step returns to the initial state.     

Gyrosynchrotron Radiation of Interplanetary CMEs

CMEs (Coronal Mass Ejections) are an important means of energy release in the solar corona. Solar Polar Orbit Radio Telescope (SPORT) is a mission being proposed for observing the propagation of interplanetary CMEs from solar polar orbit. The main payload onboard SPORT is a synthetic aperture interferometric radiometer, which receives radio emission of interplanetary CMEs. It is identified that there are mainly three radio emission mechanisms of CMEs, i.e., bremsstrahlung, gyrosynchrotron emission and plasma emission. Among these emission types, bremsstrahlung emission is the main emission mechanism of the high-density plasma clouds of interplanetary CMEs. Gyrosynchrotron emission is the continuous emission generated by high-energy electrons from CMEs, while plasma emission is the main mechanism of transient radio bursts from CMEs. In this paper, the gyrosynchrotron emission of interplanetary CMEs is focused on. Firstly, the mechanism of gyrosynchrotron emission is reviewed. Secondly, a review of the physical parameter models of background solar wind and interplanetary CMEs is presented. After these, the brightness temperature and polarization of gyrosynchrotron emission of interplanetary CMEs are calculated and analyzed. Finally, the detectability of gyrosynchrotron emission of interplanetary CMEs by radio meters is discussed briefly.      

Comment on estimating the solar proton environment that may affect Mars missions.

Estimates of the energetic proton environment for a Mars mission are generally extrapolated from the solar proton observations at 1 AU. We find that solar particle events may be divided into two general classes. Events dominated by a near-sun injection of particles onto interplanetary magnetic field lines leading to the spacecraft position represent the "classical" solar particle event associated with solar activity. This class of event will scale in radial distance by the classical power law extrapolation. The extended-interplanetary-shock source generates a maximum flux as the shock passes the detection location. This class of event typically generates maximum fluence, but in this case, the flux and fluence will not scale in the classical manner with radial distance.     

Comparing proton fluxes of central meridian SEP events with those predicted by SOLPENCO

We have developed an operational code, SOLPENCO, that can be used for space weather prediction schemes of solar energetic particle (SEP) events. SOLPENCO provides proton differential flux and cumulated fluence profiles from the onset of the event up to the arrival of the associated traveling interplanetary shock at the observer’s position (either 1.0 or 0.4 AU). SOLPENCO considers a variety of interplanetary scenarios where the SEP events develop. These scenarios include solar longitudes of the parent solar event ranging from E75 to W90, transit speeds of the associated shock ranging from 400 to 1700 km s⁻¹, proton energies ranging from 0.125 to 64 MeV, and interplanetary conditions for the energetic particle transport characterized by specific mean free paths. We compare the results of SOLPENCO with flux measurements of a set of SEP events observed at 1 AU that fulfill the following four conditions: (1) the association between the interplanetary shock observed at 1 AU and the parent solar event is well established; (2) the heliolongitude of the active region site is within 30° of the Sun–Earth line; (3) the event shows a significant proton flux increase at energies below 96 MeV; (4) the pre-event intensity background is low. The results are discussed in terms of the transit velocity of the shock and the proton energy. We draw conclusions about both the use of SOLPENCO as a prediction tool and the required improvements to make it useful for space weather purposes.     

Diffusive shock acceleration at interplanetary perpendicular shock waves: Influence of the large scale structure of turbulence on the maximum particle energy

We calculate the maximum energy that a particle can obtain at perpendicular interplanetary shock waves by the mechanism of diffusive shock acceleration. The influence of the energy range spectral index of the two-dimensional modes of the interplanetary turbulence is explored. We show that changes in this parameter lead to energies that differ in at least one order of magnitude. Therefore, the large scale structure of the turbulence is a key input if the maximum particle energy is calculated.     

从属函数在地磁扰动预报研究中的初步应用

根据1966-1982年期间有关太阳耀斑、行星际激波和地磁扰动的观测资料而建立的从属函数,对1984-1985年间的行星际闪烁观测中能证认出的耀斑-激波所引起的地磁扰动作了预报试验。结果表明:(1)磁扰开始时间预报的相对误差,δT/T≤10%的事件数为20个,占总事件数的50%,δT/T≤20%的事件占总事件数的70%以上;(2)磁扰幅度(ΣKp)大小的预报,其相对误差δΣKp/ΣKp≤30%的事件数为32个,占总事件数的80%,而δΣKp/ΣKp≥60%仅占15%.本文方法显示了一定潜力,有待从聚类分析方面进一步深入。      

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.     

多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

Contribution of geometry of interaction between interplanetary and terrestrial magnetic fields into global magnetospheric state and geomagnetic activity

The paper presents results of our study of dependence of geomagnetic activity from geoeffective parameters taking into account mutual orientation of the interplanetary magnetic field, electric field of the solar wind and geomagnetic moment. We attract a reconnection model elaborated by us made allowance for changes of geometry of the solar wind–magnetosphere interaction during annual and diurnal motions of the Earth. We take as our data base the interplanetary magnetic field and solar wind velocity measured at 1 a.u. at ecliptic plane for the period of 1963–2005 and K_p, D_st, a_m indices. Taken as a whole a geoeffective parameter suggested by us explains 95% of observed variations of the indices. Changes of the geometric factor determined by mutual orientation of the solar wind electric field and geomagnetic moment explain larger than 75% of observed statistical variations of D_st and a_m indices. Based on our results we suggest a new explanation of semi-annual and UT variation of geomagnetic activity.     

Solar flare particle injection spectra and spectra of flux maxima at the point of observation

Numerical models of impulsive solar flare particle events usually assume the radial diffusion coefficient to be independent of energy per nucleon, T, although the observations indicate a T^0.5 dependence (constant mean free path). The assumption of a constant diffusion coefficient results in a preservation of a power law injection spectrum at all radial distances throughout the event. We investigate the effect of an energy dependent diffusion coefficient on the spectrum of flux maxima at a fixed point in interplanetary space. This spectrum is harder than that of initial differential number densities close to the sun. Furthermore, the spectrum hardens with increasing radial distance which seems to be at variance with observations.     

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.     

Orientations of LASCO Halo CMEs and their connection to the flux rope structure of interplanetary CMEs

Coronal mass ejections (CMEs) observed near the Sun via LASCO coronographic imaging are the most important solar drivers of geomagnetic storms. ICMEs, their interplanetary, near-Earth counterparts, can be detected in situ, for example, by the Wind and ACE spacecraft. An ICME usually exhibits a complex structure that very often includes a magnetic cloud (MC). They can be commonly modelled as magnetic flux ropes and there is observational evidence to expect that the orientation of a halo CME elongation corresponds to the orientation of the flux rope. In this study, we compare orientations of elongated CME halos and the corresponding MCs, measured by Wind and ACE spacecraft. We characterize the MC structures by using the Grad–Shafranov reconstruction technique and three MC fitting methods to obtain their axis directions. The CME tilt angles and MC fitted axis angles were compared without taking into account handedness of the underlying flux rope field and the polarity of its axial field. We report that for about 64% of CME–MC events, we found a good correspondence between the orientation angles implying that for the majority of interplanetary ejecta their orientations do not change significantly (less than 45 deg rotation) while travelling from the Sun to the near-Earth environment.