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

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

太阳风速度的源表面二维结构

以IPS观测数据为基础，在考虑了加速剖面对太阳风源位置和速度大小的影响下，探讨了1984年1745-1755Carrington周太阳风速度在源表面上的二维分布.结果显示：观测和理论结果均显示出此期间北半球太阳风速度明显高于南半球，特别是经度0一30°，300-360°N的高速流均有不同程度的从极区向南侵的现象.低速带与磁中性线有一定程度的类似即均为单峰结构.但理论和观测给出的低速带的峰值位置与磁中性线的峰值位置略有差异，而且磁中性线的报幅大于观测和理论的结果.     

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.     

The Compton–Getting effect of energetic particles with an anisotropic pitch-angle distribution

This paper presents a simulation of anisotropy measurements by the low-energy charged particle (LECP) experiment on Voyager 1 for cases when the particle pitch-angle distribution function in the solar wind plasma reference frame is not isotropic. The model includes both the Compton–Getting anisotropy and perpendicular diffusion anisotropy that possibly exists in the upstream region of the termination shock. The results show that the Voyager 1 data cannot rule out either the model with zero solar wind speed or the one with a finite speed on qualitative basis. The determination of solar wind speed using the Compton–Getting effect is affected by the assumption of the magnetic field direction and perpendicular diffusion anisotropy. Because the pitch-angle distribution anisotropy is so large, a small uncertainty in the magnetic field direction can produce very different solar wind speeds ranging from ∼0 to >400 km/s. In fact, if the magnetic field is chosen to be in the Parker spiral direction, which is consistent with the magnetometer measurement on Voyager 1, the derived solar wind speed is still close to the supersonic value. Only the two lowest-energy channels of the LECP instrument may give a definitive answer to the solar wind speed. However, because these channels contain a very high level of cosmic ray background, an uncertainty of just a few percent in the background can entirely hamper the estimate of solar wind speed.     

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

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

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.     

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.     

WSA太阳风经验模型及其应用

Wang-Sheely-Arge （WSA）模型是对准稳态太阳风的经验和物理相结合的描述,其利用观测的日面磁图作为输入,可以提前3到4天预测L1点处的太阳风速度和行星际磁场极性.WSA模型是在WS模型的基础上经过若干改进形成的实时预报模式,之后又借鉴Distance from the Coronal Hole Boundary （DCHB）模型的参数,进一步改进了太阳风速度关系式,形成了目前常见的形式.WSA经验模型由日冕磁场模型、太阳风速度关系式和一维运动学模型三部分组成.在实际应用过程中,基本步骤包括观测磁图预处理、日冕三维磁场反演、计算日冕磁场参数、计算太阳风的速度分布和将太阳风映射到1AU等环节.在发展过程中,WSA模型经历了一些细节上的调整变化,例如观测磁图数据的来源、日冕磁场模型的类型、经验速度关系中自由系数的取值等.许多研究对如何改善模型的预报效果进行了探索.      

Ionospheric storm due to solar Coronal mass ejection in September 2017 over the Brazilian and African longitudes

Coronal mass ejection (CME) occurs when there is an abrupt release of a large amount of solar plasma, and this cloud of plasma released by the Sun has an intrinsic magnetic field. In addition, CMEs often follow solar flares (SF). The CME cloud travels outward from the Sun to the interplanetary medium and eventually hits the Earth’s system. One of the most significant aspects of space weather is the ionospheric response due to SF or CME. The direction of the interplanetary magnetic field, solar wind speed, and the number of particles are relevant parameters of the CME when it hits the Earth’s system. A geomagnetic storm is most geo-efficient when the plasma cloud has an interplanetary magnetic field southward and it is accompanied by an increase in the solar wind speed and particle number density. We investigated the ionospheric response (F-region) in the Brazilian and African sectors during a geomagnetic storm event on September 07–10, 2017, using magnetometer and GPS-TEC networks data. Positive ionospheric disturbances are observed in the VTEC during the disturbed period (September 07–08, 2017) over the Brazilian and African sectors. Also, two latitudinal chains of GPS-TEC stations from the equatorial region to low latitudes in the East and West Brazilian sectors and another chain in the East African sector are used to investigate the storm time behavior of the equatorial ionization anomaly (EIA). We noted that the EIA was disturbed in the American and African sectors during the main phase of the geomagnetic storm. Also, the Brazilian sector was more disturbed than the African sector.     

太阳风和地球磁层相互作用的两种可能类型

本文对太阳活动20周不同活动期间的太阳风参数与地磁活动性指数分别进行了相关分析,并进一步对太阳活动极大和极小年分别对B_z和太阳风参数V、T、N的时均值日方差作了分析比较。结果指出,除目前普遍认为的IMF与地磁场重联导致的磁扰外,还有一类与B_z无关,而是由高温、高速、热不均匀太阳风等离子体导致的地磁扰动类型。      

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

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

Effect of the Interplanetary Electric Field on the Magnetopause From Global MHD Simulations

The north-south component B_z of the Interplanetary Magnetic Field (IMF) and solar wind dynamic pressure P_d are generally treated as the two main factors in the solar wind that determine the geometry of the magnetosphere. By using the 3D global MHD simulations, we investigate the effect of the Interplanetary Electric Field (IEF) on the size and shape of magnetopause quantitatively. Our numerical experiments confirm that the geometry of the magnetopause are mainly determined by P_dB_z, as expected. However, the dawn-dusk IEFs have great impact on the magnetopause erosion because of the magnetic reconnection, thus affecting the size and shape of the magnetopause. Higher solar wind speed with the same B_z will lead to bigger dawn-dusk IEFs, which means the higher reconnection rate, and then results in more magnetic flux removal from the dayside. Consequently, the dayside magnetopause moves inward and flank magnetopause moves outward.      

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.     

Energetic positive ions in the cometary “foreshock” region

Ions produced by ionization of the cometary neutrals interact with the solar wind protons to produce large amplitude oscillations of the ambient magnetic field. Such oscillations are convected towards the comet at the unperturbed solar wind speed far from the shock and at a lower speed closer to the shock (due to the solar wind mass loading); hence, they can energize the incoming ions by Fermi acceleration. The spatial extension of the acceleration region is of the order of 10⁶ km and the resulting energy spectrum is harder than in the Earth's bow shock case. The energization of cometary ions produces an additional deceleration of the solar wind. It is suggested that Comet Halley may be the most efficient “cosmic ray shock” in the solar system.     

On the long-term variability of the heliosphere–magnetosphere environment

Annual means of measured and reconstructed solar, heliospheric, and magnetospheric parameters are used to infer solar activity signatures at the Hale and Gleissberg cycles timescales. Available open solar flux, modulation strength, cosmic ray flux, total solar irradiance data, reconstructed back to 1700, solar wind parameters (speed and density) and the magnitude of the heliospheric magnetic field at 1 AU, reconstructed back to 1870, as well as the time series of geomagnetic activity indices (aa, IDV, IHV), going back to 1870, have been considered. Simple filtering procedures (successive 11-, 22-, and 88-year running averages and differences between them) and scaling by the standard deviation from the average value for the common interval covered by the data show that the long-discussed variation in the 20th century (a pronounced increase since ∼1900, followed by a depression in the ‘60s and a new, slower, increase) seen in the 11-year averages of parameters such as geomagnetic activity indices and reconstructed heliospheric magnetic field strength, solar wind speed, open solar flux, is a result of the superposition in data of solar activity signatures at Hale and Gleissberg cycles timescales. The Hale and Gleissberg signals were characterized and similarities and differences in the temporal behavior of the analyzed parameters at these timescales are discussed. The similarities in the studied parameters point to a common pacing source, the solar dynamo.     

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⁻².     

嵌套闭磁场结构内CME产生和传播的数值模拟

给出了特殊类型的日冕物质抛射(CME)数值模拟定性结果,这种CME核心闭磁场结构前半部分磁力线的方向与太阳整体偶极场磁力线的方向相反．计算结果表明,这种CME核心闭磁场结构磁力线与太阳整体偶极场反向磁力线之间存在过渡磁场结构,在向外传播时过渡磁场结构所占的面积逐渐增大．这一结果可以用来解释飞船为什么能够观测到一类双极磁云,这类磁云前半部分磁场方向与太阳整体偶极场方向相反．为了模拟这一数值结果,强调需要采用包含嵌套闭磁场的冕流背景结构,并在合适的位置触发CME．     

Modeling and experimental study of the 27-day variation of galactic cosmic-ray intensity for a solar wind velocity depending on heliolongitude

We develop a three-dimensional (3-D) model of the 27-day variation of galactic cosmic-ray (GCR) intensity with a spatial variation of the solar wind velocity. A consistent, divergence-free interplanetary magnetic field is derived by solving the corresponding Maxwell equations with a variable solar wind speed, which reproduces in situ observed experimental data for the time interval to be analyzed (24 August 2007–28 February 2008). We perform model calculations for the GCR intensity using the variable solar wind and the corresponding magnetic field. Results are compatible with experimental data; the correlation coefficient between our model predictions and observed 27-day GCR variation is 0.80 ± 0.05.     

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.     

膨胀磁场和太阳风粒子相互作用三维数值模拟

采用三维模型,使用混合网格质点法HPIC(Hybrid Particle-in-Cell)对膨胀的磁场和太阳风相互作用过程进行数值模拟.研究了线圈产生的偶极子磁场在注入等离子体后和太阳风粒子的相互作用过程,并对以不同速度入射的等离子体引起的太阳风粒子的变化和磁场变化进行了比较.研究结果表明,偶极子磁场和太阳风作用时会产生弓形激波,此时磁压等于太阳风粒子的动压,当向线圈产生的偶极子磁场中注入高能等离子体时引起磁场膨胀,膨胀的磁场将会排斥太阳风粒子向外运动,从而引起弓形激波的变化,增大与太阳风相互作用的面积,并且粒子入射速度越大,磁场膨胀越明显,与太阳风相互作用愈强.