Numerical MHD modeling of the solar driver for the destabilization of a coronal helmet streamer

We investigate the forms of the solar driver which cause the destabilization of helmet streamers. Two forms of solar drivers are considered; (i) emergence of a flux-rope from sub-photospheric levels and (ii) application of a photospheric shear motion to a streamer-flux rope system. Numerical results showed that both cases exhibit the characteristics of commonly observed coronal mass ejections (CMEs), but the propagation speed of the CME is higher than the background solar wind speed when the solar driver is the emerging magnetic flux and is the same as the solar wind speed when the photospheric shear is used as the solar driver. The energy constraint allowing the magnetic field transition from a closed to an open configuration is also addressed.     

Turbulence,complexity, and solar flares

The issue of predicting solar flares is one of the most fundamental in physics, addressing issues of plasma physics, high-energy physics, and modelling of complex systems. It also poses societal consequences, with our ever-increasing need for accurate space weather forecasts. Solar flares arise naturally as a competition between an input (flux emergence and rearrangement) in the photosphere and an output (electrical current build up and resistive dissipation) in the corona. Although initially localised, this redistribution affects neighbouring regions and an avalanche occurs resulting in large scale eruptions of plasma, particles, and magnetic field. As flares are powered from the stressed field rooted in the photosphere, a study of the photospheric magnetic complexity can be used to both predict activity and understand the physics of the magnetic field. The magnetic energy spectrum and multifractal spectrum are highlighted as two possible approaches to this.     

Observational study of magnetic chirality of solar active regions

We present the evolution of magnetic field and relationship with the magnetic (current) helicity in solar active regions from a series of photospheric vector magnetograms obtained at Huairou Solar Observing Station near Beijing, and also longitudinal magnetograms by MDI of SOHO, white light and 171 Å images by TRACE and soft X-ray images by Yohkoh.The conclusions in the analysis of the formation process of complex and delta magnetic configuration in some super active regions are the following: (1) The magnetic shear and gradient provide the non-potentiality of the magnetic field of active regions reflecting the existence of electric current. (2) Some of large-scale delta active regions could be due to the emergence of highly sheared non-potential magnetic flux bundles from the subatmosphere with amount of magnetic helicity, in addition to the emergence of twisted magnetic ropes. (3) We also present some results on the study of the magnetic (current) helicity in solar active regions.     

Flux emergence,flux imbalance,magnetic free energy and solar flares

Emergence of complex magnetic flux in the solar active regions lead to several observational effects such as a change in sunspot area and flux embalance in photospheric magnetograms. The flux emergence also results in twisted magnetic field lines that add to free energy content. The magnetic field configuration of these active regions relax to near potential-field configuration after energy release through solar flares and coronal mass ejections. In this paper, we study the relation of flare productivity of active regions with their evolution of magnetic flux emergence, flux imbalance and free energy content. We use the sunspot area and number for flux emergence study as they contain most of the concentrated magnetic flux in the active region. The magnetic flux imbalance and the free energy are estimated using the HMI/SDO magnetograms and Virial theorem method. We find that the active regions that undergo large changes in sunspot area are most flare productive. The active regions become flary when the free energy content exceeds 50% of the total energy. Although, the flary active regions show magnetic flux imbalance, it is hard to predict flare activity based on this parameter alone.     

What are the physical mechanisms of eruptions and CMEs?

CMEs are due to physical phenomena that drive both, eruptions and flares in active regions. Eruptions/CMEs must be driven from initially force-free current-carrying magnetic field. Twisted flux ropes, sigmoids, current lanes and pattern in photospheric current maps show a clear evidence of currents parallel to the magnetic field. Eruptions occur starting from equilibria which have reached some instability threshold. Revisiting several data sets of CME observations we identified different mechanisms leading to this unstable state from a force free field. Boundary motions related to magnetic flux emergence and shearing favor the increase of coronal currents leading to the large flares of November 2003. On the other hand, we demonstrated by numerical simulations that magnetic flux emergence is not a sufficient condition for eruptions. Filament eruptions are interpreted either by a torus instability for an event occurring during the minimum of solar activity either by the diffusion of the magnetic flux reducing the tension of the restraining arcade. We concluded that CME models (tether cutting, break out, loss of equilibrium models) are based on these basic mechanisms for the onset of CMEs.     

Photospheric flows in the active regions (asymmetric and localized Doppler velocities)

We observed 10 active regions through their disk passage during June 25–August 25, 1988, with the Tower Vector Magnetograph (TVM) of Marshall Space Flight Center. The TVM was used in scanning mode to measure the photospheric Doppler velocities with the Line-Center-Magnetogram (LCM) technique in the spectral line of FeI 5250.2 Å. In this paper we present the result of a subset of observations obtained while the active regions were situated away from the solar limb. A wide range of magnetic complexity and associated chromospheric activity characterized these active regions. It was found that the value of zero-crossing wavelength of the integrated Stokes-V profile of two opposite magnetic polarities were different, corresponding to Doppler velocities ranging from ∼100 m s⁻¹ to ∼1475 m s⁻¹. The measurements of relative velocities between different locations, connected by magnetic flux tubes as inferred from YOHKOH soft X-ray and TRACE 171 Å Fe IX images, showed widely different values of dominant localized flows. The region of parasite polarity, which showed recurrent chromospheric activity, was blue shifted with respect to the main “magnetic element” of the same polarity. Some of them were also the sites of sheared magnetic field configuration. The magnitude of the relative velocity between the leading and following polarity is more for the active regions of higher “field asymmetry”.     

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

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

Coronal index as a solar activity index applied to space weather

All possible changes of the solar activity can be expressed by the coronal index of solar activity that represents the averaged daily power of the green corona emitted from the Sun’s visible hemisphere. The representative character of this index allows us to study long-term, intermediate and short-term variations of the Sun as a star. This index can be expressed well as a function of other solar indices. As green line reflects the distribution of the photospheric magnetic fields in the solar corona, the dependence of this index on the solar magnetic field is confirmed by means of statistical analysis of these two parameters. Daily values of the coronal index, as well as of the magnetic field data obtained from the Wilcox Solar Observatory, has been analysed by Fast Fourier analysis and Wavelet Transform analysis for the time period 1966–1998 covering more than three solar cycles. Periodicities of 11.4, 3.2, 2.3, 1.7, 1, 0.29, 0.07 and 0.04 years have been found in both parameters that means once again that the coronal index is probably related to the underlying photospheric magnetic fields and can be used as a global index of solar activity useful for Space Weather studies.     

2011年8月4日螺旋状日冕物质抛射爆发事件

利用多卫星多波段的综合观测数据,通过追踪光球表面等离子体速度分析计算了耀斑爆发前后磁螺度的变化,发现耀斑爆发前活动区中光球表面存在强的水平剪切运动,活动区磁螺度的注入主要由这种剪切运动所产生;使用CESE-MHD-NLFFF重建了耀斑爆发前后活动区的磁场位形,推测出耀斑过程中存在磁绳结构的抛射.基于这些分析,给出了这一螺旋状抛射结构的形成机制:爆发前暗条西侧足点的持续剪切运动驱动磁通量绳增加扭转,高度扭缠的通量绳与东侧足点附近的开放磁力线重联并与东侧足点断开,进而向外抛出并伴随解螺旋运动.另外,利用1AU处WIND卫星的观测数据在对应的行星际日冕物质抛射中找到典型磁云的观测特征.这表明除了传统上双足点均在太阳表面的磁云模型,这种单足点固定于太阳表面的磁通量绳爆发图景同样可能在行星系际空间形成磁云结构.研究结果对进一步认识磁云结构具有重要意义.      

Reconstruction of solar irradiance using the Group sunspot number

We present a reconstruction of total solar irradiance since 1610 to the present based on variations of the surface distribution of the solar magnetic field. The latter is calculated from the historical record of the Group sunspot number using a simple but consistent physical model. Our model successfully reproduces three independent data sets: total solar irradiance measurements available since 1978, total photospheric magnetic flux from 1974 and the open magnetic flux since 1868 (as empirically reconstructed from the geomagnetic aa-index). The model predicts an increase in the total solar irradiance since the Maunder Minimum of about 1.3 Wm⁻².     

太阳光球磁场特征物理量在质子事件短期预报中的应用

利用SOHO/MDI全日面纵向磁图, 计算了三个描述太阳活动区磁场复杂性和非势性的特征物理量, 即纵向磁场最大水平梯度Bz, 强梯度中性线长度L, 孤立奇点数目η. 为检验太阳光球磁场特征在质子事件短期预报中是否有效, 采用BP神经网络方法, 建立了基于这三个磁场特征物理量简单的太阳质子事件短期(24h)预报模型. 模 型在对2002年和2003年连续两年的样本检测中, 有很高的准确率(2002年和2003年 分别为90 %, 87.54 %)和较高的 质子事件报准率(2002年和2003年分别为60 %, 75 %),从而为光球磁场特征物理量作为质子事件预报的有效因子提供了依据.      

Correlation between solar flare productivity and photospheric vector magnetic fields

Studying the statistical correlation between the solar flare productivity and photospheric magnetic fields is very important and necessary. It is helpful to set up a practical flare forecast model based on magnetic properties and improve the physical understanding of solar flare eruptions. In the previous study ([Cui, Y.M., Li, R., Zhang, L.Y., He, Y.L., Wang, H.N. Correlation between solar flare productivity and photospheric magnetic field properties 1. Maximum horizontal gradient, length of neutral line, number of singular points. Sol. Phys. 237, 45–59, 2006]; from now on we refer to this paper as ‘Paper I’), three measures of the maximum horizontal gradient, the length of the neutral line, and the number of singular points are computed from 23990 SOHO/MDI longitudinal magnetograms. The statistical relationship between the solar flare productivity and these three measures is well fitted with sigmoid functions. In the current work, the three measures of the length of strong-shear neutral line, total unsigned current, and total unsigned current helicity are computed from 1353 vector magnetograms observed at Huairou Solar Observing Station. The relationship between the solar flare productivity and the current three measures can also be well fitted with sigmoid functions. These results are expected to be beneficial to future operational flare forecasting models.     

A theoretical model for the magnetic helicity of solar active regions

Active regions on the solar surface are known to possess magnetic helicity, which is predominantly negative in the northern hemisphere and positive in the southern hemisphere. Choudhuri et al. [Choudhuri, A.R. On the connection between mean field dynamo theory and flux tubes. Solar Phys. 215, 31–55, 2003] proposed that the magnetic helicity arises due to the wrapping up of the poloidal field of the convection zone around rising flux tubes which form active regions. Choudhuri [Choudhuri, A.R., Chatterjee, P., Nandy, D. Helicity of solar active regions from a dynamo model. ApJ 615, L57–L60, 2004] used this idea to calculate magnetic helicity from their solar dynamo model. Apart from getting broad agreements with observational data, they also predict that the hemispheric helicity rule may be violated at the beginning of a solar cycle. Chatterjee et al. [Chatterjee, P., Choudhuri, A.R., Petrovay, K. Development of twist in an emerging magnetic flux tube by poloidal field accretion. A&A 449, 781–789, 2006] study the penetration of the wrapped poloidal field into the rising flux tube due to turbulent diffusion using a simple 1-d model. They find that the extent of penetration of the wrapped field will depend on how weak the magnetic field inside the rising flux tube becomes before its emergence. They conclude that more detailed observational data will throw light on the physical conditions of flux tubes just before their emergence to the photosphere.     

The creation of fine structure by magnetic fields

The solar plasma is strongly structured by the presence of magnetic field. This structuring is manifest in the photosphere in the form of flux tubes, from the readily visible sunspots to the sub-telescopic intense tubes, so that the atmosphere is divided into strong-field media or field-free media. In the corona, by contrast, the magnetic field permeates the whole of the atmosphere and structuring consists principally of density and temperature inhomogeneities. We discuss some of the causes of magnetic structuring, including kinematic concentration, convective collapse and magnetoconvection for photospheric tubes, spicules in the chromosphere, and thermal instability for coronal loops.     

The characteristics of sharp (small-scale) boundaries of solar wind plasma and magnetic field structures

We investigate properties of large (>20%) and sharp (<10 min) solar wind ion flux changes using INTERBALL-1 and WIND plasma and magnetic field measurements from 1996 to 1999. These ion flux changes are the boundaries of small-scale and middle-scale solar wind structures. We describe the behavior of the solar wind velocity, temperature and interplanetary magnetic field (IMF) during these sudden flux changes. Many of the largest ion flux changes occur during periods when the solar wind velocity is nearly constant, so these are mainly plasma density changes. The IMF magnitude and direction changes at these events can be either large or small. For about 55% of the ion flux changes, the sum of the thermal and magnetic pressure are in balance across the boundary. In many of the other cases, the thermal pressure change is significantly more than the magnetic pressure change. We also attempted to classify the types of discontinuities observed.     

Source investigation of impulsive He-rich particle events

We have investigated the source characteristic and coronal magnetic field structure of six impulsive solar energetic particle (SEP) events selected from Wang et al. [Wang, Y.-M., Pick, M., Mason, G.M. Coronal holes, jets, and the origin of ³He-rich particle events. ApJ 639, 495, 2006] and Pick et al. [Pick, M., Mason, G.M., Wang, Y.-M., Tan, C., Wang, L. Solar source regions for ³He-rich solar energetic particle events identified using imaging radio, optical, and energetic particle observations. ApJ 648, 1247, 2006]. Some results are obtained: first, 2 events are associated with wide (≈100°) CMEs (hereafter wide CME events), another 4 events are associated with narrow (?40°) CMEs (hereafter narrow CME events); second, the coronal magnetic field configuration of narrow CME events appear more simple than that of the wide CME events; third, the photospheric magnetic field evolutions of all these events show new emergence of fluxes, while one case also shows magnetic flux cancellation; fourth, the EUV jets usually occurred very close to the footpoint of the magnetic field loop, while meter type III bursts occurred near or at the top of the loop and higher than EUV jets. Furthermore, the heights of type III bursts are estimated from the result of the coronal magnetic field extrapolations.     

Solar atmosphere wave dynamics generated by solar global oscillating eigenmodes

The solar atmosphere exhibits a diverse range of wave phenomena, where one of the earliest discovered was the five-minute global acoustic oscillation, also referred to as the p-mode. The analysis of wave propagation in the solar atmosphere may be used as a diagnostic tool to estimate accurately the physical characteristics of the Sun’s atmospheric layers.In this paper, we investigate the dynamics and upward propagation of waves which are generated by the solar global eigenmodes. We report on a series of hydrodynamic simulations of a realistically stratified model of the solar atmosphere representing its lower region from the photosphere to low corona. With the objective of modelling atmospheric perturbations, propagating from the photosphere into the chromosphere, transition region and low corona, generated by the photospheric global oscillations the simulations use photospheric drivers mimicking the solar p-modes. The drivers are spatially structured harmonics across the computational box parallel to the solar surface. The drivers perturb the atmosphere at 0.5?Mm above the bottom boundary of the model and are placed coincident with the location of the temperature minimum. A combination of the VALIIIC and McWhirter solar atmospheres are used as the background equilibrium model.We report how synthetic photospheric oscillations may manifest in a magnetic field free model of the quiet Sun. To carry out the simulations, we employed the magnetohydrodynamics code, SMAUG (Sheffield MHD Accelerated Using GPUs).Our results show that the amount of energy propagating into the solar atmosphere is consistent with a model of solar global oscillations described by Taroyan and Erdélyi (2008) using the Klein-Gordon equation. The computed results indicate a power law which is compared to observations reported by Ireland et al. (2015) using data from the Solar Dynamics Observatory/Atmospheric Imaging Assembly.     

Model-independent large-scale magnetohydrodynamic quantities in magnetic clouds

Magnetic clouds are the interplanetary manifestation of coronal mass ejections, which are transient expulsions of major quantities of magnetized plasma, from the Sun toward the heliosphere. The magnetic flux and helicity are two key physical magnitudes to track solar structures from the photosphere-corona to the interplanetary medium. To determine the content of flux and helicity in magnetic clouds, we have to know their 3D structure. However, since spacecrafts register data along a unique direction, several aspects of their global configuration cannot be observed. We present a method to estimate the magnetic flux and the magnetic helicity per unit length in magnetic clouds, directly from in situ magnetic observations, assuming only a cylindrical symmetry for the magnetic field configuration in the observed cross-section of the cloud. We select a set of 20 magnetic clouds observed by the spacecraft Wind and estimate their magnetic flux and their helicity per unit length. We compare the results obtained from our direct method with those obtained under the assumption of a helical linear force-free field. This direct method improves previous estimations of helicity in clouds.     

Solar X-ray and EUV corona and their relation to the magnetic activity

We have studied the solar magnetic cycle in corona using X-ray data from YOHKOH and Extreme Ultraviolet data from SOHO/EIT. Soft X-ray data last the period from after the maximum cycle 22 to the maximum cycle 23 (1991–2001). The SOHO/EIT Extreme Ultraviolet data are used for the period from 1996 to 2003. These data provide us a unique opportunity to look at the solar corona on the solar disc and to compare with the magnetic activity, directly.Our studies reveal a close relationship between the coronal emissions and the photospheric magnetic field in the axisymmetrical case. The evolution of coronal structures in X-ray and EUV can be considered as a proxy of the coronal magnetic field and demonstrates a development of the solar magnetic cycle in corona. It is shown that the most important feature of the coronal cycle is the forming of giant loops structure visible in X-ray and, partially, in EUV (284A) on the solar disk.     

The role of magnetic field shear in solar flares

We present observational results and their physical implications garnered from the deliberations of the FBS Magnetic Shear Study Group on magnetic field shear in relation to flares. The observed character of magnetic shear and its involvement in the buildup and release of flare energy are reviewed and illustrated with emphasis on recent results from the Marshall Space Flight Center vector magnetograph. It is pointed out that the magnetic field in active regions can become sheared by several processes, including shear flow in the photosphere, flux emergence, magnetic reconnection, and flux submergence. Modeling studies of the buildup of stored magnetic energy by shearing are reported which show ample energy storage for flares. Observational evidence is presented that flares are triggered when the field shear reaches a critical degree, in qualitative agreement with some theoretical analyses of sheared force-free fields. Finally, a scenario is outlined for the class of flares resulting from large-scale magnetic shear; the overall instability driving the energy release results from positive feedback between reconnection and eruption of the sheared field.