Evolution of coronal mass ejections in the early stage

This work reports the investigation of two coronal mass ejections (CME) observed in white light, H, EUV and X-ray by various instruments both in space and on ground on February 18, 2003 and January 19, 2005, respectively. The white light coronal images show that the first CME began with the rarefaction of a region above the solar limb and was followed by the formation of its leading edge at the boundary of the rarefying region at altitude of 0.46 R from the solar surface. The rarefaction coincided the slow rising phase of the filament eruption, and the CME leading edge was observed to form as the filament eruption started to accelerate apparently. In the early stage of the second CME, a bright loop was first observed above the solar limb with height of 0.37 R in EUV images. We found that the more gradual CMEs initial process, the larger the timing difference between CMEs and their associated flares. The lower part of the filament brightened in H images as the filament rose to a certain height. These brightenings imply that the filament may be heated by magnetic reconnection below the filament in the early stage of the eruption. We suggest that the possible mechanism which led to the formation of the CME leading edge and cavity is magnetic reconnection which occurred under the filament when it reached a certain height.     

Two-velocity structure observed in the inner solar wind

Measurements of the motion of plasma density inhomogeneities in the inner solar wind are presented. The speeds were estimated using a cross-correlation analysis of radio frequency fluctuations of the Galileo spacecraft measured simultaneously at widely spaced ground stations. The radial projections of the correlation baselines on the pattern plane were of the order of several thousand kilometers. For cross-correlation functions calculated with comparatively short averaging times, we find that a pronounced two-velocity configuration is occasionally observed over the range of heliocentric distances 20 R < R < 40 R. The typical mean speed for such observations is about 300–400 km/s and the difference between the two predominant speeds is about 150–200 km/s. These results may indicate that the density fluctuations are associated with slow magnetosonic waves propagating in opposite directions at the local speed of sound in the reference frame moving with the mean solar wind speed. Quite reasonable estimates of the solar wind speed and speed of sound are obtained from this model. Another possible explanation of the two-velocity structures is that two independent solar wind streams are present simultaneously along different segments of the radio ray path.     

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.     

Relation between solar wind velocity and properties of its source region

Different kinds of coronal holes are sources of different kind of solar winds. A successful solar wind acceleration model should be able to explain all those solar winds. For the modeling it is important to find a universal relation between the solar wind physical parameters, such as velocity, and coronal physical parameters such as magnetic field energy. To clarify the physical parameters which control the solar wind velocity, we have studied the relation between solar wind velocity and properties of its source region such as photospheric/coronal magnetic field and the size of each coronal hole during the solar minimum. The solar wind velocity structures were derived by using interplanetary scintillation tomography obtained at Solar-Terrestrial Environment Laboratory, Japan. Potential magnetic fields were calculated to identify the source region of the solar wind. HeI 1083 nm absorption line maps obtained at Kitt Peak National Solar Observatory were used to identify coronal holes. As a result, we found a relation during solar minimum between the solar wind velocity and the coronal magnetic condition which is applicable to different kind of solar winds from different kind of coronal holes.     

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.     

一个基于观测的太阳风背景模型

建立由太阳光球磁场和日冕偏振亮度等观测约束的单流体太阳风模型,包括日冕和太阳风的等离子体密度、速度和磁场,温度还有待于以后处理.这里采用高山观测台(HAO)MKⅢ的日冕偏振亮度(pB)在1.36Rs上的观测概图,根据Guhathakurta在1996年发展的日冕电子密度反演模型确定日冕的电子密度分布.同时采用Wilcox太阳观测台(WSO)的光球磁场视向分量的观测概图作为底部边界,根据Zhao等在1994年发展的水平电流-电流片(HCCS)模型得到全球磁场.Phillips在1995年及McComas在2003年分别用Ulysses第一次和第二次跨极飞行的观测发现,归一化到1 AU的太阳风动量流密度除了在10°～30°的纬度范围内略低以外几乎不变.根据这一结论,结合已经得到的密度数据,就可以得到日冕和太阳风的速度.将上面的模型应用于1918卡林顿自转周稳态太阳风的研究,结果与太阳活动极小期的观测基本相符,但是与观测相比较低速高密度区偏大,因此密度模型还有待改进.      

Coronal radio-sounding detection of a CME during the 1997 Galileo solar conjunction

Frequency fluctuations of the Galileo S-band radio signal were recorded nearly continuously during the spacecraft’s solar conjunction from December 1996 to February 1997. A strong propagating disturbance, most probably associated with a coronal mass ejection (CME), was detected on 7 February when the radio ray path proximate point was on the west solar limb at about 54 solar radii from the Sun. The CME passage through the line of sight is characterized by a significant increase in the fluctuation intensity of the recorded frequency and by an increase in the plasma speed from about 234 km s⁻¹ up to about 755 km s⁻¹. These velocity estimates are obtained from a correlation analysis of frequency fluctuations recorded simultaneously at two widely-separated ground stations. The density turbulence power spectrum is found to steepen behind the CME front. The Galileo radio-sounding data are compared with SOHO/LASCO observations of the CME in the corona and with WIND spacecraft data near the Earth’s orbit.     

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.     

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

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

Periodic Variations of Drag Coefficient for the ANDE Spherical Satellites During its Lifetime

A drag coefficient (C_D) inversion method is introduced to study the variations of the drag coefficient for orbital satellites with spherical geometry. Drag coefficients of the four micro satellites in the Atmospheric Neutral Density Experiment (ANDE) are compiled out with this new method. The Lomb-Scargle Periodgram (LSP) analysis of the four ANDE satellites' C_D series has shown that there are obvious 5, 7, 9, and 27 days' period in those data. Interesting results are found through comparing the LSP analysis with series of the daily solar radio flux at 10.7 cm (F_10.7 index), the Ap index, and the daily averaged solar wind speed at 1AU. All series in the same time interval have an obvious period of about 27 days, which has already been explained as the association with the 27 days' solar rotation. The oscillating periods less than 27 days are found in series of C_D, Ap and solar wind speed at 1AU, e.g., the 5, 7, 9 days period. However, these short periods disappeared in the time series of F_10.7 index. The same periodicities of 5, 7, 9 days in Ap and solar wind are presented at the same time interval during the declining phase of solar cycle 23. While in the ascending phase of solar cycle 24, these short oscillations are not so obvious as that in the declining phase of solar cycle 23. These results provide definite evidence that the C_D variations with period of 5, 7 and 9 days are produced by a combination of space weather effects caused by the solar wind and geomagnetic activity.      

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.     

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.     

Radio sounding of the solar wind acceleration region with spacecraft signals

Data from coronal radio-sounding experiments carried out on various interplanetary spacecraft are used to derive the empirical radial dependence of solar wind velocity and density at heliocentric distances from 3 to 60 solar radii for heliolatitudes below 60° and for low solar activity. The radial dependencies of solar wind power and acceleration are derived from these results. Summaries of the radial behavior of characteristic parameters of the solar wind turbulence (e.g., the spectral index and the inner and outer turbulence scales), as well as the fractional density fluctuation, are also presented. These radio-sounding results provide a benchmark for models of the solar wind in its acceleration region.     

根据行星际闪烁(IPS)的双频互谱估计太阳风速度

本文讨论了观测频率为232MHz和327MHz时,利用互谱的方法估计太阳风速度的可能性、考虑路径积分效应,结果表明在17°≤ε≤50°范围内,对于几种可能的太阳风速度分布,由互谱所得结果与太阳风速度值差别小于15%.      

磁场强度对日冕定态结构的影响

以二维ＭＨＤ模型及时变方法为基础,内外边界完整的设影特征线边界条件,考察了太阳日冕大气的定态结构随偶极场强度的变化情况。模拟结果表明：随着偶极场强度的增加,磁场对太阳风的约束增强,低纬闭磁场打开程度减少,高纬与低纬区速度差增加,并且在阿尔文马赫数为１的点附近达到最大,速度过度区变陡;随着日心距离增加,低纬区宽度减小,速度过渡区变陡,可定性解释Ｕｌｙｓｓｅｓ飞船的新观测事实。     

Interplanetary shocks and sudden impulses during solar maximum (2000) and solar minimum (1995–1996)

In this work a study is performed on the correlation between fast forward interplanetary shock parameters at 1 Astronomical Unit and sudden impulse (SI) amplitudes in the H-component of the geomagnetic field, for periods of solar activity maximum (year 2000) and minimum (year 1995–1996). Solar wind temperature, density and speed, and total magnetic field, were taken to calculate the static pressures (thermal and magnetic) both in the upstream and downstream sides of the shocks. The variations of the solar wind parameters and pressures were then correlated with SI amplitudes. The solar wind speed variations presented good correlations with sudden impulses, with correlation coefficients larger than 0.70 both in solar maximum and solar minimum, whereas the solar wind density presented very low correlation. The parameter better correlated with SI was the square root dynamic pressure variation, showing a larger correlation during solar maximum (r = 0.82) than during solar minimum (r = 0.77). The correlations of SI with square root thermal and magnetic pressure were smaller than with the dynamic pressure, but they also present a good correlation, with r > 0.70 during both solar maximum and minimum. Multiple linear correlation analysis of SI in terms of the three pressure terms have shown that 78% and 85% of the variance in SI during solar maximum and minimum, respectively, are explained by the three pressure variations. Average sudden impulse amplitude was 25 nT during solar maximum and 21 nT during solar minimum, while average square root dynamic pressure variation is 1.20 and 0.86 nPa^1/2 during solar maximum and minimum, respectively. Thus on average, fast forward interplanetary shocks are 33% stronger during solar maximum than during solar minimum, and the magnetospheric SI response has amplitude 20% higher during solar maximum than during solar minimum. A comparison with theoretical predictions (Tsyganenko’s model corrected by Earth’s induced currents) of the coefficient of sudden impulse change with solar wind dynamic pressure variation showed excellent agreement, with values around 17 nT/nPa^1/2.     

Solar wind-driven electron radiation belt response functions at 100-min time scales

We present a simple yet numerically robust technique, using autoregressive linear filters, to remove unwanted “colored noise” from solar wind and radiation belt electron data at sub-daily resolution. The remaining signal is then studied using finite impulse response linear prediction filters to represent the driven portion of the linear dynamics that describe the coupling between solar wind speed and electron flux. Sub-daily resolution response profiles covering magnetic L-shells between 1.1 and 8.0 R_E are presented which are consistent with daily resolution response functions. Namely, while there is strong global coherence governing electron flux dynamics, there are at least two distinct responses. The first response is an immediate dropout of electrons between L = 4 and L = 7 that is at least a partly adiabatic effect associated with enhancements in the ring current. This is followed by a 1–2 day delayed enhancement across the same L-shells that is likely a result of increased radial diffusion. The second response is an immediate enhancement seen between L = 3 and L = 4 with a typical duration of less than one day. Plausible explanations for this second response are briefly discussed, but neither empirical nor theoretical evidence can establish conclusively a definite physical cause. Finally, the response profiles show significant solar cycle and seasonal dependencies, indicating that better model output might be achieved with: (1) additional simultaneous solar wind inputs; (2) more sophisticated dynamical model structures capable of incorporating non-linear feedback; and/or (3) time-adaptive linear filters that can track non-stationary dynamics in time.     

Solar wind turbulence during the solar cycle deduced from Galileo coronal radio-sounding experiments

Seven coronal radio-sounding campaigns were carried out during the active lifetime of the Galileo spacecraft in the years 1994–2002. The observational data analyzed in the present work are S-band frequency fluctuation measurements recorded during the solar conjunctions at different phases of solar activity cycle #23, specifically: periods near solar maximum (three conjunctions), near solar minimum (three conjunctions) and during the ascending phase (one conjunction). These data are all applicable to low heliographic latitudes, i.e. to the slow solar wind. The rms frequency fluctuation and power-law index of the frequency fluctuation temporal spectra are determined as a function of heliocentric distance. The turbulence power spectrum tends to be flatter inside ca. 20 solar radii during all phases of the solar cycle. This coincides with a transition in the flow from the inner acceleration region to the outer region of constant velocity. The radial falloff rate and absolute level of the rms frequency fluctuation are essentially invariant over the solar cycle.     

二维磁流体动力学平衡日冕(Ⅱ)

基于解析和数值相结合的方法,进一步讨论了非均匀引力场中日冕的二维磁流体动力学平衡。对临界点进行了比较仔细的处理。得到了包含闭场区、中性片和开场区的大尺度日冕磁场位形,闭场区和中性片构成冕流结构。在高纬和低纬地区几个太阳半径之外,等离子体径向流动速度超过了局地声速和局地Alfvén速度。在1AU处,太阳风速度可达到400kms_-1以上      

日冕扰动传播的数值模拟(Ⅲ)太阳风的影响

本文数值研究了二维磁流体动力学平衡基态下开场区日面冷物质径向喷射所引起的日冕动力学响应。结果表明:(1)在高密度环前方有一弱扰动区近似以Alfvén速度向外传播;(2)高密度环前缘移动速度随着径向距离而增加, 其增加值近似为局地太阳风速度;(3)高密度环中等离子体的最大径向速度约在4个太阳半径处趋于局地逃逸速度;(4)对于强开放场, 环形结构在θ方向上没有明显的扩张。这些结果可以更好地解释伴随有日珥的日冕物质抛射事件。