Three-dimensional global simulation of multiple ICMEs’ interaction and propagation from the Sun to the heliosphere following the 25–28 October 2003 solar events

This study performs simulations of interplanetary coronal mass ejection (ICME) propagation in a realistic three-dimensional (3D) solar wind structure from the Sun to the Earth by using the newly developed hybrid code, HAFv.2+3DMHD. This model combines two simulation codes, Hakamada–Akasofu–Fry code version 2 (HAFv.2) and a fully 3D, time-dependent MHD simulation code. The solar wind structure is simulated out to 0.08 AU (18 Rs) from source surface maps using the HAFv.2 code. The outputs at 0.08 AU are then used to provide inputs for the lower boundary, at that location, of the 3D MHD code to calculate solar wind and its evolution to 1 AU and beyond. A dynamic disturbance, mimicking a particular flare’s energy output, is delivered to this non-uniform structure to model the evolution and interplanetary propagation of ICMEs (including their shocks). We then show the interaction between two ICMEs and the dynamic process during the overtaking of one shock by the other. The results show that both CMEs and heliosphere current sheet/plasma sheet were deformed by interacting with each other.     

Helios飞船观测到的行星际日冕物质抛射事件

根据异常低的质子温度判据，从Heliosl和2飞船等离子体观测数据中(0．3-1AU)识别出大约160个行星际日冕物质抛射事件(ICME)，并在统计意义上分析了ICME在内日球空间的传播和演化规律．     

Propagation and interaction of interplanetary transient disturbances. Numerical simulations

We study the heliocentric evolution of ICME-like disturbances and their associated transient forward shocks (TFSs) propagating in the interplanetary (IP) medium comparing the solutions of a hydrodynamic (HD) and magnetohydrodynamic (MHD) models using the ZEUS-3D code [Stone, J.M., Norman, M.L., 1992. Zeus-2d: a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. i – the hydrodynamic algorithms and tests. Astrophysical Journal Supplement Series 80, 753–790]. The simulations show that when a fast ICME and its associated IP shock propagate in the inner heliosphere they have an initial phase of about quasi-constant propagation speed (small deceleration) followed, after a critical distance (deflection point), by an exponential deceleration. By combining white light coronograph and interplanetary scintillation (IPS) measurements of ICMEs propagating within 1 AU [Manoharan, P.K., 2005. Evolution of coronal mass ejections in the inner heliosphere: a study using white-light and scintillation images. Solar Physics 235 (1–2), 345–368], such a critical distance and deceleration has already been inferred observationally. In addition, we also address the interaction between two ICME-like disturbances: a fast ICME 2 overtaking a previously launched slower ICME 1. After interaction, the leading ICME 1 accelerates and the tracking ICME 2 decelerates and both ICMEs tend to arrive at 1 AU having similar speeds. The 2-D HD and MHD models show similar qualitative results for the evolution and interaction of these disturbances in the IP medium.     

Evolution of solar wind structures from 0.72 to 1 AU

Understanding the evolution of solar wind structures in the inner heliosphere as they approach the Earth is important to space weather prediction. From the in situ solar wind plasma and magnetic field measurements of Pioneer Venus Orbiter (PVO) at 0.72 AU (1979–1988), and of Wind/Advanced Composition Explorer (ACE) missions at 1 AU (1995–2004), we identify and characterize two major solar wind structures, stream interaction regions (SIRs) and interplanetary coronal mass ejections (ICMEs). The average percentage of SIRs occurring with shocks increases significantly from 3% to 24% as they evolve from 0.72 to 1 AU. The average occurrence rate, radial extent, and bulk velocity variation of SIRs do not change from 0.72 to 1 AU, while peak pressure and magnetic field strength both decrease with the radial evolution of SIRs. Within the 0.28 AU distance from the orbit of Venus to that of Earth, the average fraction of ICMEs with shocks increases from 49% to 66%, and the typical radial extent of ICMEs expands by about a fraction of 1.4, with peak pressure and magnetic field strength decreasing significantly. The mean occurrence rate and expansion velocity of ICMEs do not change from 0.72 to 1 AU.     

Corotating events in the energy range 4–13 MeV as observed on board Helios 1 and 2 in 1975 and 1976

Intensity increases of 4 – 13 MeV protons correlated with recurrent high speed solar wind streams were observed on board Helios 1 and 2 at heliocentric distances between 0.3 and 1.0 AU. Study of events which occurred in the time period 1975 – 1976 showed that recurrent fast solar wind streams were sometimes not accompanied by corotating events and small events were more frequently observed than large corotating events, which occurred only occasionally. The explanation for these phenomena seems to be the variation of the number of suprathermal particles injected into the acceleration process. Radial gradients of corotating events from February till April 1976 were investigated. The gradient changed its sign in successive solar rotations and exhibited strong time dependence inside 1 AU. Clearly negative gradients were measured in corotating events which occurred in February and April 1976. These observations showed that probably there were no stationary corotating structures in interplanetary space. We suggest that this effect is related to local disturbances of solar origin, e. g. flare initiated shocks.     

太阳风速度及日球电流片对CME渡越时间的影响

基于1996-2005年88个引起重大地磁暴的CME(日冕物质抛射)事件、1996-2000年的47个CME事件以及1997-2002年的29个全晕状CME事件,结合ACE卫星在1AU处的太阳风和行星际磁场观测资料以及Wilcox Solar Observatory(WSO)天文台的太阳光球层磁图,分析了背景太阳风速度和日球电流片对CME到达1AU处渡越时间预报误差的影响.结果表明,背景太阳风速度与CME渡越时间误差并没有明显的相关性,在考虑了磁云通量管轴相对黄道面夹角的影响后相关性依然不明显.然而日球电流片对CME渡越时间却有明显的影响,对于初速度较小的异侧CME事件,其渡越时间大于同侧事件;而对于具有较大初速度的CME事件,异侧事件的渡越时间明显小于同侧事件.研究结果表明,CME与太阳风以及日球电流片的相互作用并不是简单的对流相互作用,造成高速CME异侧事件快于同侧事件到达地球的因素非常复杂,有待深入研究.      

太阳风粘性效应的重新考察

利用Helios飞船对太阳风高速流的等离子体和磁场观测,本文考察了0.3-1AU作用在高速流上各种体积力的径向分布.结果表明:普遍采用的在碰撞为主等离子体中导出的经典粘性系数表式不适用于太阳风高速流;为使惯性力与各体积力之和相等,必须有大小跟这一空间范围内的太阳重力相近的粘性减速力.本文提出一种新的方法,导出了太阳风高速流中实际粘性力的径向分布,并求出了实际粘性系数的径向分布.结果发现,相对于经典粘性系数,1AU处太阳风高速流的粘性系数的减小不小于十分之一.      

激波在行星际介质中的能量耗散

导出激波下游介质相对上游介质能流通量增量公式,并由HeliosA,B飞船太阳风观测资料得出不同流速太阳风流中各参量随日心距的幂律变化。以此作为背景值分别计算出磁能、内能、动能和总能在不同日心距离处的能量耗散率。结果指出激波后介质以动能增加为主,内能次之,磁能最少;总能耗率在近日处较大,但下降较快。从0.3-1.0AU,不同强度激波总能耗随初始Alfvén激波数A₁₀增大而增大,对A₁₀从2.0-10.0的计算结果与观测值一致。      

SEP events and multi-spacecraft observations: Constraints on theory

Fifteen solar energetic particle (SEP) events have been analyzed using proton flux data recorded by the Helios 1, Helios 2, and IMP 8 spacecraft in the energy range ∼4–40 MeV during 1974–1982. For each of the events at least two of the spacecraft have their nominal magnetic footpoint within 20° in heliocentric longitude from each other. The SEP events are sub-grouped as a function of their heliocentric longitudinal separation and heliocentric radial distance from the SEP associated solar flare and several case studies are presented in this paper. Main results concerning their usage in estimating the SEP radial dependence are given. Moreover, we investigate the behavior of the third not connected spacecraft in order to study the dependence of the proton flux as a function of flare location. It is found that the contribution of the longitudinal gradient in determining variations in the SEP proton flux is particularly relevant for spacecraft having their magnetic connection footpoint separated from the flare between 30° and 50°.     

A polytropic model for the solar wind

The solar wind fills the heliosphere and is the background medium in which coronal mass ejections propagate. A realistic modelling of the solar wind is therefore essential for space weather research and for reliable predictions. Although the solar wind is highly anisotropic, magnetohydrodynamic (MHD) models are able to reproduce the global, average solar wind characteristics rather well. The modern computer power makes it possible to perform full three dimensional (3D) simulations in domains extending beyond the Earth’s orbit, to include observationally driven boundary conditions, and to implement even more realistic physics in the equations. In general, MHD models for the solar wind often make use of additional source and sink terms in order to mimic the observed solar wind parameters and/or they hide the not-explicitly modelled physical processes in a reduced or variable adiabatic index. Even the models that try to take as much as possible physics into account, still need additional source terms and fine tuning of the parameters in order to produce realistic results. In this paper we present a new and simple polytropic model for the solar wind, incorporating data from the ACE spacecraft to set the model parameters. This approach allows to reproduce the different types of solar wind, where the simulated plasma variables are in good correspondence with the observed solar wind plasma near 1 AU.     

第24太阳活动周地球附近磁云与非磁云事件的统计分析

为研究第24太阳活动周中磁云（Magnetic Clouds,MC）与非磁云（Non-Magnetic Clouds,non-MC）的等离子体性质及其对空间天气的影响,使用1AU处的观测数据对2008-2015年168个ICME事件进行统计与分析,其中认证出磁云事件68个,占总数的40.48%.通过分析磁云与非磁云等离子体参数对空间天气环境的影响及与太阳活动的关系,整体性质的对比及在第23和24太阳活动周中性质的对比,可以发现:在第24太阳周中,磁云引起的磁暴强度普遍大于非磁云,南北向磁场分量是引起磁暴的重要参数;磁云数和太阳黑子数有很好的相关性,非磁云数与行星际日冕物质抛射总数及黑子数的相关性稍弱,磁云数在太阳周的不同阶段表现出不同的分布特性;磁云的磁场强度和南向磁场分量整体强于非磁云,两者质子温度、密度等参数差异不大.第24周磁云事件引起的地磁效应整体上弱于第23周磁云事件,这与第24周磁云事件最大南向磁场分量、传播速度以及质子温度整体小于第23周磁云事件有关.      

Magnetosphere response to the 2005 and 2006 extreme solar events as observed by the Cluster and Double Star spacecraft

The four identical Cluster spacecraft, launched in 2000, orbit the Earth in a tetrahedral configuration and on a highly eccentric polar orbit (4–19.6 R_E). This allows the crossing of critical layers that develop as a result of the interaction between the solar wind and the Earth’s magnetosphere. Since 2004 the Chinese Double Star TC-1 and TC-2 spacecraft, whose payload comprise also backup models of instruments developed by European scientists for Cluster, provided two additional points of measurement, on a larger scale: the Cluster and Double Star orbits are such that the spacecraft are almost in the same meridian, allowing conjugate studies. The Cluster and Double Star observations during the 2005 and 2006 extreme solar events are presented, showing uncommon plasma parameters values in the near-Earth solar wind and in the magnetosheath. These include solar wind velocities up to ∼900 km s⁻¹ during an ICME shock arrival, accompanied by a sudden increase in the density by a factor of ∼5 and followed by an enrichment in He⁺⁺ in the secondary front of the ICME. In the magnetosheath ion density values as high as 130 cm⁻³ were observed, and the plasma flow velocity there reached values even higher than the typical solar wind velocity. These resulted in unusual dayside magnetosphere compression, detection of penetrating high-energy particles in the magnetotail, and ring current development following several successive injections of energetic particles in the inner magnetosphere, which “washed out” the previously formed nose-like ion structures.     

A new parameter to define interplanetary coronal mass ejections

The interplanetary manifestations of coronal mass ejections, ICMEs, have many signatures in the solar wind but none of these signatures in the velocity, density, temperature, magnetic field, plasma composition or energetic particles uniquely and unambiguously identifies the occurrence of an ICME. Different investigators identify different events when confronted with the same data. Herein, we present a single physical parameter that combines information from multiple plasma components and that holds the promise of defining a beginning and an end of the region of influence ICME and an indication of the location of the encounter with the ICME relative to its central meridian. This parameter is the total plasma pressure perpendicular to the magnetic field, consisting of the sum of the magnetic pressure and plasma kinetic or thermal pressure. It provides a vehicle for classifying the nature of the ICME encounter and, in many cases, provides an unambiguous start and stop time of the event. However, it does not provide a start and stop time for any embedded flux rope. This identification depends on examination of the magnetic field.     

How to address the accuracy of empirical magnetic field models?

Empirical magnetic field models are compared with high-altitude magnetic field measurements and results from an MHD simulation. Comparison of the T96 model and observations from GOES-8 and GOES-9 shows that if the observed solar wind and IMF parameters are used to compute the model field, the model field is more stretched than the observed field. On the other hand, if measurements made by one spacecraft are used to find the model parameters that give a best-fit field at that location, the RMS error can be reduced also at the other spacecraft four hours away in local time. Comparison of T96 model and MHD simulation results shows that the empirical models have a thinner current sheet than the MHD simulation, but that the lobe field values are quite similar to each other. Furthermore, a comparison of an event-oriented, modified T89 model and MHD simulation by Pulkkinen et al. [2000] reveals that if the empirical model is constructed by fitting to in-situ measurements, the resulting model is very similar to the MHD simulation magnetic field. These results indicate that an efficient method of utilizing the present-day empirical models is to select model parameters based on measurements from a few individual points.     

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.     

太阳风质子束流密度的径向变化

利用Helios2飞船的数据,对太阳风速度分布中质子束流部分与整个质子的密度之比随日心距离的变化做了分析.为了排除碰撞因素的影响,有针对性地分析了太阳风高速流(600相似文献   

行星际日冕物质抛射地磁效应研究的支持向量机方法初步研究

行星际日冕物质抛射（Interplanetary Coronal Mass Ejection,ICME）与地球磁层相互作用并带来地磁暴等地磁扰动.从Richardson和Cane提供的近地球ICME列表中筛选出ICME事件集,基于ICME扰动期间的行星际等离子体与磁场数据提取出特征.通过计算各特征的费舍尔分值（Fisher Score）,对这些特征进行选择,发现行星际磁场南北向分量持续时间小于-10nT且激波等扰动所带来的ICME扰动开始时,太阳风速度的增量等特征与ICME事件的地磁效应密切相关.这与现有的传统统计研究结果一致.以这些特征为基础,训练得到的径向基函数支持向量机能够以0.78±0.08的准确率判断ICME事件是否会产生中等及以上强度的地磁暴（Dst ≤-50nT）.      

Time-dependent cosmic ray modulation

Time-dependent cosmic ray modulation is calculated over multiple solar cycles using our well established two-dimensional time-dependent modulation model. Results are compared to Voyager 1, Ulysses and IMP cosmic ray observations to establish compatibility. A time-dependence in the diffusion and drift coefficients, implicitly contained in recent expressions derived by , ,  and , is incorporated into the cosmic ray modulation model. This results in calculations which are compatible with spacecraft observations on a global scale over consecutive solar cycles. This approach compares well to the successful compound approach of Ferreira and Potgieter (2004). For both these approaches the magnetic field magnitude, variance of the field and current sheet tilt angle values observed at Earth are transported time-dependently into the outer heliosphere. However, when results are compared to observations for extreme solar maximum, the computed step-like modulation is not as pronounced as observed. This indicates that some additional merging of these structures into more pronounced modulation barriers along the way is needed.     

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.