联合太阳和行星际物理参数预测行星际激波能否到达地球

选取第23太阳活动周（1997—2006年）期间542例由太阳爆发活动驱动的行星际激波事件,分析确定了太阳源头和行星际空间中影响行星际激波能否到达地球轨道的关键物理参数;在此基础上,建立了预测行星际激波能否到达地球的新预报模型（EdEaSPM）. 回溯预报结果表明,EdEaSPM模型的预报成功率约为66%,略高于国际一流预报模型的预报成功率;EdEaSPM模型的虚报率未超过50%,改善了当前国际主流模型虚报率较大的情况;对于偏度指标,虽然当前所有模型的偏度值均大于1,但EdEaSPM模型的偏度值最接近于1且明显小于其他模型的偏度值;EdEaSPM模型的其他评价指标也都高于国际主流模型的相应指标. 此外,选取2012年期间的激波事件对EdEaSPM模型进行了预报检验,预测结果与实际情况吻合. EdEaSPM模型不仅能够提前约1～3天进行预报,而且预报效果与国际一流模型具有可比性,尤其是在提高预报成功率及降低虚报率方面具有一定优势.      

不确定理论下竞争失效系统的可靠性分析

针对有新研发产品，故障数据较少的复杂系统，提出了不确定环境下自然退化和外部冲击相互独立的竞争失效模型。考虑了系统同时遭受自然退化和外部冲击，连续的自然退化用一个不确定过程刻画，冲击到达的时间间隔和每次冲击对系统造成的损坏量分别用2个不同的不确定变量来刻画。运用不确定理论，分别在极端冲击模型、累积冲击模型、δ冲击模型下，研究了系统的确信可靠度，结果表明:在有新研发的产品、故障数据较少的复杂系统，用不确定理论的方法来描述模型更合适，并通过数值分析显示了模型的有效性。      

An empirical model to predict the 1-AU arrival of interplanetary shocks

We extend the empirical coronal mass ejection (CME) arrival model of Gopalswamy et al. [Gopalswamy, N. et al. Predicting the 1-AU arrival times of coronal mass ejections, J. Geophys. Res. 106, 29207, 2001] to predict the 1-AU arrival of interplanetary (IP) shocks. A set of 29 IP shocks and the associated magnetic clouds observed by the Wind spacecraft are used for this study. The primary input to this empirical shock arrival model is the initial speed of white-light CMEs obtained using coronagraphs. We use the gas dynamic piston–shock relationship to derive the ESA model which provides a simple means of obtaining the 1-AU speed and arrival times of interplanetary shocks using CME speeds.     

运用一维太阳风MHD CE/SE模型预报行星际激波到达时间

在一维太阳风磁流体(MHD)数值模拟中, 应用时空守恒元和解元(Conservation Element/Solution Element, CE/SE)方法, 建立了一个行星际激波扰动传播模型(1D-MHD (CE/SE)模型), 用来预报行星际激波到达时间. 收集了1997年2月至2002年8月间的137个激波事件, 对激波到达地球轨道附近的传播时间进行了预测, 并将预报结果与STOA, ISPM, HAFv.2以及SPM模型所得结果进行比较. 对于相同的样本事件, 1D-MHD (CE/SE)模型给出的渡越时间平均绝对值误差并不大于其他4个模型, 且该模型预报的相对误差小于10 %的事件占25.6 %, 小于30 %的事件占69.3 %, 小于50 %的事件占87.6 %, 其预报精度与其他模型相比基本相当. 这表明该模型在空间天气的激波到达时间预报方面有潜在的应用价值.      

Propagation of interplanetary shocks into the Earth’s magnetosphere

This paper is devoted to the study of propagation of disturbances caused by interplanetary shocks (IPS) through the Earth’s magnetosphere. Using simultaneous observations of various fast forward shocks by different satellites in the solar wind, magnetosheath and magnetosphere from 1995 till 2002, we traced the interplanetary shocks into the Earth’s magnetosphere, we calculated the velocity of their propagation into the Earth’s magnetosphere and analyzed fronts of the disturbances. From the onset of disturbances at different satellites in the magnetosphere we obtained speed values ranging from 500 to 1300 km/s in the direction along the IP shock normal, that is in a general agreement with results of previous numerical MHD simulations. The paper discusses in detail a sequence of two events on November 9th, 2002. For the two cases we estimated the propagation speed of the IP shock caused disturbance between the dayside and nightside magnetosphere to be 590 km/s and 714–741 km/s, respectively. We partially attributed this increase to higher Alfven speed in the outer magnetosphere due to the compression of the magnetosphere as a consequence of the first event, and partially to the faster and stronger driving interplanetary shock. High-time resolution GOES magnetic field data revealed a complex structure of the compressional wave fronts at the dayside geosynchronous orbit during these events, with initial very steep parts (10 s). We discuss a few possible mechanisms of such steep front formation in the paper.     

电离层等效电流体系对行星际激波的响应研究

行星际激波是导致地球磁层-电离层系统发生扰动的重要原因之一,其可以通过对磁层-电离层系统电流体系的改变来影响地磁变化.本文采用全球三维磁流体力学数值模拟方法,分析了行星际激波作用下电离层等效电流体系的即时响应.模拟结果表明,在激波作用下伴随着异常场向电流对的产生,电离层在午前午后出现一对反向的等效电流涡.这对涡旋一边向极侧和夜侧运动,一边经历强度增强和减弱直至消失的过程.激波过后等效电流体系图像逐渐演化为激波下游行星际条件控制的典型图像.这个响应过程与行星际激波强度有关,激波强度越强,则反向的等效电流涡旋强度越大,寿命也就越短.      

Propagation of normal and faster CMEs in the interplanetary medium

We have analyzed 101 Coronal Mass Ejection (CME) events and their associated interplanetary CMEs (ICMEs) and interplanetary (IP) shocks observed during the period 1997–2005 from the list given by Mujiber Rahman et al. (2012). The aim of the present work is to correlate the interplanetary parameters such as, the speeds of IP shocks and ICMEs, CME transit time and their relation with CME parameters near the Sun. Mainly, a group of 10 faster CME events (V_INT > 2200 km/s) are compared with a list of 91 normal events of Manoharan et al. (2004). From the distribution diagrams of CME, ICME and IP shock speeds, we note that a large number of events tends to narrow towards the ambient (i.e., background) solar wind speed (∼500 km/s) in agreement with the literature. Also, we found that the IP shock speed and the average ICME speed measured at 1 AU are well correlated. In addition, the IP shock speed is found to be slightly higher than the ICME speed. While the normal events show CME travel time in the range of ∼40–80 h with a mean value of 65 h, the faster events have lower transit time with a mean value of 40 h. The effect of solar wind drag is studied using the correlation of CME acceleration with interplanetary (IP) acceleration and with other parameters of ICMEs. While the mean acceleration values of normal and faster CMEs in the LASCO FOV are 1 m/s², 18 m/s², they are −1.5 m/s² and −14 m/s² in the interplanetary medium, respectively. The relation between CME speed and IP acceleration for normal and faster events are found to agree with that of  and  except slight deviations for the faster events. It is also seen that the faster events with less travel time face higher negative acceleration (>−10 m/s²) in the interplanetary medium up to 1 AU.     

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

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

Solar and interplanetary parameters of CMEs with and without type II radio bursts

We have analyzed 101 CMEs, and their associated ICMEs and interplanetary (IP) shocks observed during the period 1997–2005. The main aim of the present work is to study the interplanetary characteristics of metric and DH type II associated CMEs such as, shock strength, IP shock speed, ICME speed, stand off distance and transit time. Among these 101 CMEs, 38 events show both metric and DH type II bursts characteristics. There are no metric and DH type II association for 52 events. While DH type II alone is found in 7 cases, metric type II alone is found in 4 events. It is found that the mean speeds of CMEs increase progressively from CMEs without type II events to CMEs associated with metric and DH type IIs as suggested by Gopalswamy et al. (2005). In addition, we found that the speeds of ICMEs and IP shocks progressively increase in the following order: events without metric and DH type IIs, events with metric alone, events with DH alone and events with both metric and DH type IIs. Similarly the Mach number is found to increase in the same order. While there is not much change in the stand-off distance among these cases, it is minimum (∼18 R_⊙) for CMEs with speed greater than 2200 km/s. The above results confirm that more energetic CMEs can produce both metric and DH type IIs for which the interplanetary parameters such as mean values of ICME speed and IP shock speed and Mach number are found to be higher.     

具有变点的不确定竞争失效退化系统的可靠性建模

针对具有变点的竞争失效退化系统的可靠性进行了建模，系统同时遭受连续内部磨损退化和外部冲击，内部磨损退化由一个不确定过程表示，外部冲击过程是一个更新回报过程，即外部冲击到达的时间间隔和每次冲击的损坏量都是随机变量。外部环境条件的变化导致随机次数的冲击后损坏量服从的分布发生变化，通过运用不确定理论和机会理论，分别推导出系统在3种不同冲击模式下的确信可靠度公式，并通过数值算例验证了模型的正确性和有效性。      

How do interplanetary shock impact angles control the size of the geoeffective magnetosphere?

In this paper, we investigate temporal and spatial magnetosphere response to the impact of interplanetary (IP) shocks with different inclinations and speeds on the Earth’s magnetosphere. A data set with more than 500 IP shocks is used to identify positive sudden impulse (SI⁺) events as expressed by the SuperMAG partial ring current index. The SI⁺ rise time (RT), defined as the time interval between compression onset and maximum SI⁺ signature, is obtained for each event. We use RT and a model suggested by Takeuchi et al. (2002) to calculate the geoeffective magnetospheric distance (GMD) in the shock propagation direction as a function of shock impact angle and speed for each event. GMD is a generalization of the geoeffective magnetosphere length (GML) suggested by Takeuchi et al. (2002), defined from the subsolar point along the X line toward the tail. We estimate statistical GMD and GML values which are then reported for the first time. We also show that, similarly to well-known results for RT, the highest correlation coefficient for the GMD and impact angle is found for shocks with high speeds and small impact angles, and the faster and more frontal the shock, the smaller the GMD. This result indicates that the magnetospheric response depends heavily on shock impact angle. With these results, we argue that the prediction and forecasting of space weather events, such as those caused by coronal mass ejections, will not be accurately accomplished if the disturbances’ angles of impact are not considered as an important parameter within model and observation scheme capabilities.     

Relationships between energetic storm particle events and interplanetary shocks driven by full and partial halo coronal mass ejections

We have analysed energetic storm particle (ESP) events in 116 interplanetary (IP) shocks driven by front-side full and partial halo coronal mass ejections (CMEs) with speeds $>$400 km s^?1during the years 1996–2015. We investigated the occurrence and relationships of ESP events with several parameters describing the IP shocks, and the associated CMEs, type II radio bursts, and solar energetic particle (SEP) events. Most of the shocks (57 %) were associated with an ESP event at proton energies $>$1 MeV.The shock transit speeds from the Sun to 1 AU of the shocks associated with an ESP event were significantly greater than those of the shocks without an ESP event, and best distinguished these two groups of shocks from each other. The occurrence and maximum intensity of the ESP events also had the strongest dependence on the shock transit speed compared to the other parameters investigated. The correlation coefficient between ESP peak intensities and shock transit speeds was highest (0.73 $\pm$ 0.04) at 6.2 MeV. Weaker dependences were found on the shock speed at 1 AU, Alfvénic and magnetosonic Mach numbers, shock compression ratio, and CME speed. On average all these parameters were significantly different for shocks capable to accelerate ESPs compared to shocks not associated with ESPs, while the differences in the shock normal angle and in the width and longitude of the CMEs were insignificant.The CME-driven shocks producing energetic decametric–hectometric (DH) type II radio bursts and high-intensity SEP events proved to produce also more frequently ESP events with larger particle flux enhancements than other shocks. Together with the shock transit speed, the characteristics of solar DH type II radio bursts and SEP events play an important role in the occurrence and maximum intensity of ESP events at 1 AU.     

A first step towards proton flux forecasting

We present a preliminary version of a potential tool for real time proton flux prediction which provides proton flux profiles and cumulative fluence profiles at 0.5 and 2 MeV of solar energetic particle events, from their onset up to the arrival of the interplanetary shock at the spacecraft position (located at 1 or 0.4 AU). Based on the proton transportation model by Lario et al. [Lario, D., Sanahuja, B., Heras, A.M. Energetic particle events: efficiency of interplanetary shocks as 50 keV E < 100 MeV proton accelerators. Astrophys. J. 509, 415–434, 1998] and the magnetohydrodynamic shock propagation model of Wu et al. [Wu, S.T., Dryer, M., Han, S.M. Non-planar MHD model for solar flare-generated disturbances in the Heliospheric equatorial plane. Sol. Phys. 84, 395–418, 1983], we have generated a database containing “synthetic” profiles of the proton fluxes and cumulative fluences of 384 solar energetic particle events. We are currently validating the applicability of this code for space weather forecasting by comparing the resulting “synthetic” flux profiles with those of several real events.     

Relative fluxes of shock and prompt peaks in SEP event time profiles

Peak fluxes are an important property of gradual solar energetic particle (SEP) event time profiles from both astro/heliophysical and applications perspectives. However, the peak flux in an event may occur at the event onset, or at the time of the interplanetary shock arrival (the ESP or energetic storm particles). This makes an important difference in the interpretation of the peak flux, and in any attempts to characterize or model it. This paper describes a study of SEP data sets from ACE, IMP-8 and GOES toward determining the relative properties of these peak fluxes for protons with energies near 1, 10, and 50 MeV. The results suggest that for gradual events with both peaks, the ESP peak often dominates at 1 MeV energies and is dominant about half the time at 10 MeV. Moreover, the prompt peak fluxes can be used to estimate the shock peak (ESP event) up to days ahead, especially in the lower energy range.     

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.     

Solar wind low-energy energetic ion enhancements: A tool to forecast large geomagnetic storms

Predicting the occurrence of large geomagnetic storms more than an hour in advance is an important, yet difficult task. Energetic ion data show enhancements in flux that herald the approach of interplanetary shocks, usually for many hours before the shock arrival. We present a technique for predicting large geomagnetic storms (Kp  7) following the arrival of interplanetary shocks at 1 AU, using low-energy energetic ions (47–65 keV) and solar wind data measured at the L1 libration point. It is based on a study of the relationship between energetic ion enhancements (EIEs) and large geomagnetic storms by Smith et al. [Smith, Z., Murtagh, W., Smithtro, C. Relationship between solar wind low-energy energetic ion enhancements and large geomagnetic storms. J. Geophys. Res. 109, A01110, 2004. doi:10.1029/ 2003JA010044] using data in the rise and maximum of solar cycle 23 (February 1998–December 2000). An excellent correlation was found between storms with Kp  7 and the peak flux of large energetic ion enhancements that almost always (93% of time in our time period) accompany the arrival of interplanetary shocks at L1. However, as there are many more large EIEs than large geomagnetic storms, other characteristics were investigated to help determine which EIEs are likely to be followed by large storms. An additional parameter, the magnitude of the post-shock total magnetic field at the L1 Lagrangian point, is introduced here. This improves the identification of the EIEs that are likely to be followed by large storms. A forecasting technique is developed and tested on the time period of the original study (the training data set). The lead times, defined as the times from the arrival of the shock to the start of the 3-h interval of maximum Kp, are also presented. They range from minutes to more than a day; the average for large storms is 7 h. These times do not include the extra warning time given when the EI flux cross the high thresholds ahead of the shock. Because the data-stream used in the original study is no longer available, we extended the original study (1998–2000) to 2001, in order to: (a) investigate EIEs in 2001; (b) present a validation of the technique on an independent data set; (c) compare the results based on the original (P1) energy channel to those of the replacement (P1′) and (d), determine new EIE thresholds for forecasting geomagnetic storms using P1′ data. The verification of this P1′ training data set is also presented, together with lead times.     

异常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事件上升时间与行星际磁场方向关系不明显.      

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.     

基于LSTM的TTE网络速率约束流量预测

时间触发以太网（TTE）中的速率约束（RC）流量为事件触发流量，在RC流量动态调度的应用场景下，若能预测未来短时间内数条RC流量到达交换节点的序列，使交换节点提前进行调度决策，以减小RC流量时延，提高网络吞吐量。对RC流量到达序列预测问题进行了研究，建立了RC流量的到达序列模型，提出了基于长短期记忆网络（LSTM）算法的RC流量预测算法。利用OMNET++工具进行TTE网络仿真，得到多组混合关键性配置下RC流量的传输数据；以此作为输入样本对预测算法进行训练和测试。实验结果显示，LSTM算法在RC流量预测问题的准确率达到了70%以上。通过对比实验说明所提算法适用于RC流量预测场景。      

Response of the low latitude geomagnetic field to the major proton event of November 2001

A major solar flare eruption occurred at 16:20 UT on 4 November 2001, followed by strong solar radiation storm and proton event recorded by the SOHO and other interplanetary satellites. Coronal mass ejection associated with the flare event triggered an interplanetary shock, which impacted the geomagnetic field after about 33 h. The shock impact was quite intense to produce a SSC magnitude of 80 nT in the low latitude ground magnetic records followed by sharp and deep main phase (Dst −300 nT) in the first stage, following the density (N_p) enhancement. High time resolution digital magnetic field data from the equatorial and low latitude stations in India are analyzed to study the influence of various IP parameters on the intensity and duration of the magnetic storm. A double step storm was found to be in progress caused by the multiple injections. During the period of recovery, after a period of 8 h, a third stage of depression in the ground magnetic field was set in, which corresponded to the southward directed Bz. The energy transfer processes associated with the event is presented.