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

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

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.     

MHD simulation of magnetic field configuration above the active region NOAA 10365

The current sheet (CS) creation before a flare in the vicinity of a singular line above the active region NOAA 10365 is shown in numerical experiments. Such a way the possibility of energy accumulation for a solar flare is demonstrated. These data and results of observation confirm the electrodynamical solar flare model that explains solar flares and CME appearance during CS disruption. The model explains also all phenomena observed in flares. For correct reproduction of the real boundary conditions the magnetic flux between spots should be taken into account. The full system of 3D MHD equations are solved using the PERESVET code. For setting the boundary conditions the method of photospheric magnetic maps is used. Such a method permits to take into account all evolution of photospherical magnetic field during several days before the flare.     

Inferring plasma flow velocities from photospheric vector magnetic field observations for the investigation of flare onsets

The amount of emergence and submergence of magnetized plasma and the horizontal motion of the footpoints of flux tubes might be crucial for the dynamics of the solar atmosphere. Although the rate of flux emergence and submergence can be observationally determined near the polarity inversion line (Chae et al., 2004), the same is not true for regions away from the PIL. Also, the horizontal motions cannot be directly measured in the solar photosphere. In this sense, the evolution of the photospheric magnetic field provides valuable information which can be used to estimate photospheric plasma flows since magnetic field and plasma are closely associated (frozen-in-condition). We used three methods to estimate the photospheric plasma motion from magnetic field observations. The methods were applied to photospheric vector magnetic field data of active region NOAA 9077, observed by the Huairou Solar Observing Station (HSOS) of the National Astronomical Observatories of China before and after the ‘Bastille Day’ flare on July 13th and 14th, 2000.     

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.     

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

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

结合光球磁场特征物理量的质子事件短期预报

利用描述太阳活动区光球磁场复杂性和非势性特征的三个物理量(纵向磁场最大水平梯度|▽_hB_z |_m, 强梯度中性线长度L, 孤立奇点数目η)建立了质子事件短期预报模型, 验证了磁场特征物理量对质子事件短期预报的有效性. 目前已建立或使用的太阳质子事件短期预报模型中仍然没有正式将磁场特征物理量作为预报因子. 由于太阳质子事件是小概率事件, 其物理产生机制尚不完全清楚, 这些预报模型往往存在虚报率偏高或报准率偏低的问题. 本文试图将原有质子事件模型所用的传统因子与磁场特征物理量结合起来, 利用神经网络方法建立一个更为有效的质子事件短期预报模型. 利用1997--2001年的训练数据集1871个样本建立了输入层为传统预报因子的模型A以及输入层为传统预报因子和磁场特征物理量的模型B. 通过对2002--2003年973个样本的测试数据集进行模拟预报发现, 模型A与B在具有相同质子事件报准率的情况下, 模型B的虚报率明显降低. 这进一步验证了磁场特征物理量在质子事件短期预报中的作用, 进而可以加强对太阳质子事件的实际预报能力.      

太阳宁静区色球磁场观测

利用中国科学院北京天文台怀柔太阳磁场望远镜,对日面中心宁静区光球和色球磁场进行了长时间的积分观测。通过对光球、色球以及色球不同层次的长时间积分的观测发现,网络磁元从光球到色球扩展不大,并且部分内网络磁元升到了色球。这些结果对描述太阳磁场的两大模型提出了挑战。      

Motion of the HXR sources in solar flares: Yohkoh images and statistics

Using the Yohkoh Hard X-Ray Telescope (HXT) data, we have examined motions of the hard X-ray (HXR) sources during 72 solar flares occurred from 1991 September to 2001 December. In these flares, we have found 198 intense sources that are presumably the chromospheric footpoints (FPs) of flare loops. The average velocity V and the velocity dispersion σ were determined by a linear regression for these sources. For 80% of them, the ratio of V to 3σ is larger than 1, strongly suggesting that the regular motions of the HXR sources dominate their chaotic motions.For 43 of 72 flares, coalignment of the HXT images with the photospheric magnetograms allows us to consider the HXR sources located on the both sides of the photospheric neutral line (NL) as the FP sources, and to distinguish between three main types of the FP motions. The type I is the motions of the HXR sources preferentially away from and nearly perpendicular to the NL. Less than 5% of the flares show this pattern of motion. In the type II, the sources move mainly along the NL in anti-parallel directions. Such motions have been found in 26% of flares. The type III involves a similar pattern of motions as the type II but all the HXR sources move in the same direction along the NL. Flares of this type constitute 30% of the flares. About 19% of flares can be described as a combination of these basic types. The remaining 20% of flares seem to be more complicated or less regular in the motion scale under consideration. An interpretation of results is suggested.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

Modeling large solar flares

The basic ideas to model the large solar flares are reviewed and illustrated. Some fundamental properties of potential and non-potential fields in the solar atmosphere are recalled. In particular, we consider a classification of the non-potential fields or, more exactly, related electric currents, including reconnecting current layers. The so-called ‘rainbow reconnection’ model is presented with its properties and predictions. This model allows us to understand main features of large flares in terms of reconnection. We assume that in the two-ribbon flares, like the Bastille-day flare, the magnetic separatrices are involved in a large-scale shear photospheric flow in the presence of reconnecting current layers generated by a converging flow.     

太阳质子耀斑的一个统计性质

本文研究太阳质子耀斑相对于太阳光球大尺度平均磁场中性线的分布, 给出一个新得到的质子耀斑的统计性质。研究太阳质子事件及其源耀斑的统计性质, 是太阳物理和空间物理学的重要前沿课题。从太阳活动预报及地球空间环境预报研究的角度看, 一个重要的问题是质子耀斑在      

A study of a long duration B9 flare-CME event and associated shock

We present and discuss here the observations of a small long duration GOES B-class flare associated with a quiescent filament eruption, a global EUV wave and a CME on 2011 May 11. The event was well observed by the Solar Dynamics Observatory (SDO), GONG H$\alpha$, STEREO and Culgoora spectrograph. As the filament erupted, ahead of the filament we observed the propagation of EIT wave fronts, as well as two flare ribbons on both sides of the polarity inversion line (PIL) on the solar surface. The observations show the co-existence of two types of EUV waves, i.e., a fast and a slow one. A type II radio burst with up to the third harmonic component was also associated with this event. The evolution of photospheric magnetic field showed flux emergence and cancellation at the filament site before its eruption.     

Physics of ion acceleration in the solar flare on 2005 September 7 determines γ-ray and neutron production

Relativistic neutrons were observed by the neutron monitors at Mt. Chacaltaya and Mexico City and by the solar neutron telescopes at Chacaltaya and Mt. Sierra Negra in association with an X17.0 flare on 2005 September 7. The neutron signal continued for more than 20 min with high statistical significance. Intense emissions of γ-rays were also registered by INTEGRAL, and during the decay phase by RHESSI. We analyzed these data using the solar-flare magnetic-loop transport and interaction model of Hua et al. [Hua, X.-M., Kozlovsky, B., Lingenfelter, R.E. et al. Angular and energy-dependent neutron emission from solar flare magnetic loops, Astrophys. J. Suppl. Ser. 140, 563–579, 2002], and found that the model could successfully fit the data with intermediate values of loop magnetic convergence and pitch-angle scattering parameters. These results indicate that solar neutrons were produced at the same time as the γ-ray line emission and that ions were continuously accelerated at the emission site.     

Detailed analysis of dynamic evolution of three Active Regions at the photospheric level before flare and CME occurrence

We present a combined analysis of the applications of the weighted horizontal magnetic gradient (denoted as ${\mathit{WG}}_M$ in Korsós et al. (2015)) method and the magnetic helicity tool (Berger and Field, 1984) employed for three active regions (ARs), namely NOAA AR 11261, AR 11283 and AR 11429. We analysed the time series of photospheric data from the Solar Dynamics Observatory taken between August 2011 and March 2012. During this period the three ARs produced a series of flares (eight M- and six X-class) and coronal mass ejections (CMEs). AR 11261 had four M-class flares and one of them was accompanied by a fast CME. AR 11283 had similar activities with two M- and two X-class flares, but only with a slow CME. Finally, AR 11429 was the most powerful of the three ARs as it hosted five compact and large solar flare and CME eruptions. For applying the ${\mathit{WG}}_M$ method we employed the Debrecen sunspot data catalogue, and, for estimating the magnetic helicity at photospheric level we used the Space-weather HMI Active Region Patches (SHARP’s) vector magnetograms from SDO/HMI (Solar Dynamics Observatory/Helioseismic and Magnetic Imager). We followed the evolution of the components of the ${\mathit{WG}}_M$ and the magnetic helicity before the flare and CME occurrences. We found a unique and mutually shared behaviour, called the U-shaped pattern, of the weighted distance component of ${\mathit{WG}}_M$ and of the shearing component of the helicity flux before the flare and CME eruptions. This common pattern is associated with the decreasing-receding phases yet reported only known to be a necessary feature prior to solar flare eruption(s) but found now at the same time in the evolution of the shearing helicity flux. This result leads to the conclusions that (i) the shearing motion of photospheric magnetic field may be a key driver for solar eruption in addition to the flux emerging process, and that (ii) the found decreasing-approaching pattern in the evolution of shearing helicity flux may be another precursor indicator for improving the forecasting of solar eruptions.     

Emission heights of coronal bright points on Fe XII radiance map

The study of coronal bright points (BPs) is important for understanding coronal heating and the origin of the solar wind. Previous studies indicated that coronal BPs have a highly significant tendency to coincide with magnetic neutral lines in the photosphere. Here we further studied the emission heights of the BPs above the photosphere in the bipolar magnetic loops that are apparently associated with them. As BPs are seen in projection against the disk their true emission heights are unknown. The correlation of the BP locations on the Fe XII radiance map from EIT with the magnetic field features (in particular neutral lines) was investigated in detail. The coronal magnetic field was determined by an extrapolation of the photospheric field (derived from 2-D magnetograms obtained from the Kitt Peak observatory) to different altitudes above the disk. It was found that most BPs sit on or near a photospheric neutral line, but that the emission occurs at a height of about 5 Mm. Some BPs, while being seen in projection, still seem to coincide with neutral lines, although their emission takes place at heights of more than 10 Mm. Such coincidences almost disappear for emissions above 20 Mm. We also projected the upper segments of the 3-D magnetic field lines above different heights, respectively, on to the tangent x–y plane, where x is in the east–west and y in the south–north direction. The shape of each BP was compared with the respective field-line segment nearby. This comparison suggests that most coronal BPs are actually located on the top of their associated magnetic loops. Finally, we calculated for each selected BP region the correlation coefficient between the Fe XII intensity enhancement and the horizontal component of the extrapolated magnetic field vector at the same x–y position in planes of different heights, respectively. We found that for almost all the BP regions we studied the correlation coefficient, with increasing height, increases to a maximal value and then decreases again. The height corresponding to this maximum was defined as the correlation height, which for most bright points was found to range below 20 Mm.     

Evolution of flare ribbons and energyrelease

We examined the relation between the evolutions of the H flare ribbons and the released magnetic energiesat a solar flare which occurred on 2001 April 10. This is the first study to evaluate the released energy quantitatively, based on the magnetic reconnection model, and by using the data obtained with the multi wavelength observation. We measured the, photospheric magnetic field strengths and the separation speeds of the fronts of the H flare ribbon, and compared them the nonthermal behaviors observed in HXRs and microwaves. Those nonthermal radiation sources tell us when and where large energy releases occur. Then, by using the photospheric and chromospheric features, we estimated the released magnetic energy at the flare. The estimated energy release rates at the H kernels associated with the HXR sources are locally large enough to explain the difference between the spatial distribution the H kernels and the HXR sources. Their temporal evolution of the energy release rates also shows peaks corresponding to HXR bursts.     

太阳活动区磁场扭绞与耀斑产率

本文以1972年10月的太阳活动区McMath 12094为范例, 研究了活动区磁场扭绞与耀斑产率的关系.先在常α无力场模型假定下, 以观测到的活动区光球磁场为边值, 对活动区在日面中心附近4天(10月28—31日), 推算出代表活动区磁场平均扭绞程度的无力因子α, 从而外推出活动区在这4天的三维磁力线形态.然后以这些资料为基础, 进一步讨论了活动区磁场演化特征, 磁场扭绞与耀斑产率的关系, 并且近似用单极场模型估算了通过活动区前导大黑子A的电流、电流密度以及因大黑子逆时针旋转造成磁场扭绞所贮存的能量.本文主要结论为:(1)活动区McMath 12094从10月27日起保持较强扭绞, 10月30日达到极大, 10月31日后扭绞减弱.活动区磁场扭绞的主要原因是光球中的磁流体力学作用所导致的前导大黑子A的逆时针旋转。(2)代表活动区磁场平均扭绞程度的无力因子α与活动区耀斑产率同步变化, 表明活动区磁场扭绞与耀斑产率成正相关.(3)通过活动区前导大黑子A的本影电流为4.3—6.6×1012A, 因扭绞产生的自由能贮存为0.44—1.11×1032erg.活动区中的电流密度达到0.96—1.47×10A·m-2.这样高的电流密度可能是该活动区高耀斑产率的重要原因.      

Fluorescent excitation of photospheric Fe K-α emission during solar flares

The Bent Crystal Spectrometer on the NASA Solar Maximum Mission satellite provides high spectral and temporal resolution observations of the Fe Kα lines. We have analysed spectra from almost 50 solar flares that occurred during 1980. These data strongly support fluorescent excitation of photospheric iron by photons of E > 7.11 keV emitted by the hot coronal plasma produced during the flare. After comparison of the data with a model, we discuss the observed Kα line widths, estimates of the size of the emitting region, the height of the coronal source and the photospheric iron abundance.     

Reorganization of the solar corona following a C4.7 flare

Yohkoh X-ray images, multifrequency two-dimentional observations of the Nancay Radioheliograph, Kitt Peak and Mees magnetograms provide a unique set of data with which to study a C4.7 long-duration flare that was observed close to the equator (S07, W11) on 25 Oct. 1994 at 09:49 UT. Linear force-free field extrapolations indicate a very high degree of non-potentiality in the active region. The X-ray flare started with the expansion of spectacular twisted loops. Fifteen minutes after the flare onset sporadic radio (type III) bursts were observed spreading over an area of almost 1/3 of the solar disc and two remote X-ray brightenings appeared over quiet regions of opposite magnetic polarity located in on opposite hemispheres of the Sun. In the close vicinity of these remote brightenings two coronal holes formed. The timing and location of these events combined with the overall magnetic configuration provide evidence for a large-scale magnetic reconnection occurring between the expanding twisted loops and the overlying huge loops which inter-connect quiet solar regions.