2. 中国电波传播研究所, 新乡 453007
2. China Research Institute of Radio Wave Propagation, Xinxiang 453007, China
太阳耀斑的研究和预报工作一直是各国天文工作者研究的重点.看似平静的太阳,实际上每时每刻都在发生剧烈的运动,太阳耀斑正是太阳最剧烈的活动.太阳耀斑突然爆发后,喷射出大量的可见光、X射线、伽马射线、紫外线、带电粒子等,将会破坏电离层,影响地球空间环境,干扰地球磁场和高空电离层,使它失去反射无线电电波的功能,造成短波无线通信、广播质量严重下降或完全中断,手机和互联网崩溃,全球定位导航系统失灵,电力系统瘫痪等,影响人们正常的生产和生活.伴随太阳耀斑喷射的高能粒子1~2天后到达地球附近,还有可能会引起一系列的地球物理现象,如电离层暴、地磁暴等,会给在轨运行的卫星、宇宙飞船上的设备、特别是出舱活动的航天员带来非常大的危害.为了能有效预报太阳耀斑以及利用太阳耀斑爆发所释放的能量,几十年来国内外天文学家一直都致力于研究与太阳耀斑相关的工作(栗志等,1999;王水和魏奉思,2007;朱岗崑等,2009).例如,近年来提出的基于序列数据的太阳耀斑预测方法以及利用GPS对太阳耀斑进行监测和预报(马宇倩等,2004;孔庆颜等,2010;李强等,2012).但这些方法都还存在着经济成本较高、数据处理繁琐等问题.因此,本文在了解VLF相位在极区的传播规律和南极中山站的电离层特征的基础上(贺龙松等,2000;沈长寿等,2005;王健等,2006;徐中华等,2006),通过观测1991-02-25南极地区特大太阳耀斑爆发时,以及随后粒子沉降发生过程中VLF信号的相位变化情况,分析其相位变化与太阳耀斑爆发及粒子沉降的相关性,为直接预测太阳耀斑及粒子沉降事件发生的时间和级别提供依据.
1 观测原理和观测系统
VLF信号的频率在3~30 kHz之间,波长长,衰减低,稳定性好,因此在矿产资源探测、海底通讯、地震预测以及天文预报等方面都有广泛的应用.近几十年来,大量的研究资料表明,VLF信号相位对太阳活动辐射出的X射线、伽马射线、高能粒子十分敏感,因此,研究VLF信号相位变化与太阳耀斑之间相关性,对预报太阳耀斑,减小其对人类的破坏,具有非常重要的意义.VLF信号在地球表面和低电离层之间的波导内传播,通常用波导模(Budden,1961)的概念来解释,一般比较复杂,日出和日落时有其特殊的变化规律(Crombie,1966;Schoute-Vannneck,1974;牛有田等,2009;王虎等,2011),但对于远距离的白天情况下,VLF信号只有一阶模,相位变化比较平稳.当太阳发生耀斑时,太阳突然喷射大量的X射线及高能带电粒子都会引起低电离层D层电离浓度突然增加,D层的导电率增大,D层对VLF信号的有效反射高度突然降低,VLF信号的相位突然超前,即发生了所谓的SPA现象.
图 1甚低频信号观测系统框图.单鞭天线接收由Omega甚低频无线电导航系统G台或E台发射的频率为10.2 kHz的VLF信号相位,送入到VLF信号锁相接收机.精度优于1×10-11的5 MHz的正弦波由铷原子频标提供,作为接收机的本地标准频率.记录仪记录VLF信号锁相接收机输出的相位信号.系统测相精度为±1 μs.
图 2为1991-02-25及1991-02-26南极中山10.2 kHz VLF相位日变化曲线.实线表示1991-02-25有太阳耀斑爆发的Omega系统G台——南极中山站之间VLF信号相位变化实测曲线;虚线表示1991-02-26未有太阳耀斑发生的G台——中山站VLF信号相位正常变化实测曲线.
表 1给出了1991-02-25特大太阳耀斑爆发时至1991-02-27粒子沉降发生后,这一时间段在南极中山站观测记录的VLF相位变化的部分数据.
太阳耀斑爆发后喷射出X射线以光速传播,经过8分18秒到达地球附近,即可引起SPA现象.从图 2中可以看出,1991-02-25发生太阳耀斑时VLF信号相位变化曲线与1991-02-26未发生太阳耀斑VLF信号的相位变化相比,在1991-02-25 13:09LT至16:00LT之间明显发生了相位超前现象.我们推断这次太阳耀斑从1991-02-25 13:09LT开始,到13:21LT达到峰值,说明这次太阳耀斑经过约12分钟的时间,对D层的电离达到高峰.之后,X射线流量逐渐减小,D层的电离率减小,离子的复合加快.到16:00LT,X射线流量恢复常态,D层中空气分子的电离和离子的复合又达到了新的动态平衡状态,即这次太阳耀斑对VLF信号传播的扰动基本上结束了,总共持续了大约3个小时.
根据1—8X射线的峰值流量一般可分为X、M、C、B等级别.设X射线峰值流量为F,则当1×10-4≤F<1×10-3 erg/(cm2·s)时,太阳耀斑为B级;当1×10-3≤F<1×10-2 erg/(cm2·s)时,太阳耀斑为C级;当1×10-2≤F<1×10-1 erg/(cm2·s)时,太阳耀斑为M级;当F≥1×10-1 erg/(cm2·s)时,太阳耀斑为X级.在图 2中可以看到,1991-02-25 VLF相位变化量为44 μs(1 μs=4.9°),文献(王素琴等,1988)中给出D层有效反射高度的降低Δh与VLF相位Δφ之间的关系式,修正后可计算出Δh=15.5 km.修正公式为
其中,信号频率f=10.2 kHz,λ为信号波长,地球半径a=6369 km,G台——中山站的路径距离D=5314.6 km,白天电离层高度h=70 km.文献(刘万通,1987)给出了太阳X射线通量密度F与电离层降低高度Δh之间的拟合方程.根据文献中提出的方法,将数据再次整理拟合,可得出拟合方程为
由公式(2)计算出X射线的通量为1.1×10-1erg/(cm2·s),根据太阳耀斑的划分级别,推测此次太阳耀斑可达到X1.1级.从NGDC发布的报告中可查到,1991-02-25 8:00 UT~11:00 UT期间发生特大太阳耀斑,级别达到X1.2.由于地方时比世界时早5 h,LT≈UT+5,经过换算后发现,根据VLF信号相位异常推测出的太阳耀斑爆发时间和级别与卫星观测数据基本吻合.
粒子沉降事件大概太阳耀斑爆发后的1~2天到达地球附近.在持续的观测中发现,VLF信号相位在1991-02-26无异常变化,从1991-02-27 22:57LT开始,VLF相位又开始发生超前现象,直至23:45LT才恢复正常,历时48分钟,相位超前达幅度28 μs.我们推测在这段时间发生了粒子沉降事件.另外,从观测数据中看到,太阳耀斑爆发造成VLF相位发生异常的开始时间是1991-02-25 13:09LT,随之而来的粒子沉降事件使VLF相位异常超前的开始时间是1991-02-27 22:57LT,此期间共经历了57.8小时.日地距离大约为149597870千米,因此,可计算出太阳风速度为718.5 km/s.经对比发现,观测VLF相位来推测和计算的高能粒子沉降事件发生时间及太阳风速度,与NGDC发布的卫星观测结果吻合性较好.
4 结 论
以上分析表明,VLF相位对太阳耀斑爆发及高能粒子沉降非常敏感.太阳耀斑一旦发生,VLF相位便会出现异常超前现象,并且超前的幅度与太阳耀斑的级别之间存在直接关系.因此,通过观测VLF信号的相位变化情况,可以推测是否发生太阳耀斑.根据VLF相位变化的幅度,结合相应的拟合公式,能够推算太阳耀斑的级别.同样可以由VLF相位的变化来确定粒子沉降事件发生的时间及太阳风速度.本文中VLF相位异常的时间与卫星观测的数据还存在一定误差,所以仍然需要对观测系统、VLF信号频率等可能会给观测结果带来影响的因素进行调整,通过多次试验观测分析,使试验结果达到最佳.另外,电离层的降低高度和通量密度的拟合公式也还需几千次甚至上万次的观测数据来进行分析修正,只有这样才能使预测结果足够精确.如果这两个问题都能解决,则可以仅仅通过观测VLF信号的相位变化,来确定太阳耀斑级别,准确预测粒子沉降时间,使地球上居住的人类及时采取防范措施,把太阳耀斑爆发的带来的危害降到最低.
[1] | Budden K G. 1961. The Wave-guide Mode Theory of Wave Propagation [M]. London: Cambridge University Press. |
[2] | Crombie D D. 1966. Further observations of sunrise and sunset fading of very-low-frequency signals[J]. Radio Science, 1(1): 47-51. |
[3] | He L S, Liu R Y, Liu S L, et al. 2000. Overall properties of F region around solar minimum at Zhongshan station, Antarctica[J]. Chinese J. Geophys. (in Chinese), 43(3): 289-295. |
[4] | Kong Q Y, Feng J, Liu W, et al. 2010. Development on short-term prediction of iononspheric parameters[J]. Progress in Geophysics (in Chinese), 25(6): 1968-1976, doi: 10. 3969/j. issn. 1004-2903. 2010. 06. 011. |
[5] | Li Q, Ning B Q, Zhao B Q, et al. 2012. Applications of the CMONOC based GNSS data in monitoring and investigation of ionospheric space weather [J]. Chinese J. Geophys. (in Chinese), 55(7): 2193-2202, doi: 10. 6038/j. issn. 0001-5733. 2012. 07. 006. |
[6] | Li Z, Yang X N, Zhou W H. 1999. Some new methods of the automatic observation, forecast, data processing and research for solar activity [J]. Progress in Geophysics (in Chinese), 14(S1): 135-141. |
[7] | Liu W T. 1987. Correlative analysis between sudden phase anomalies of VLF signals and solar X-ray events [J]. Chinese Journal of Space Science (in Chinese), 7(3): 20-24. |
[8] | Ma Y Q, Wang H Y, Zhang C M, et al. 2004. A conjunctive study of solar flare 20010402 and related solar proton events by the observation of SZ2/XD and ZY1/CBMC [J]. Chinese J. Geophys. (in Chinese), 47(5): 737-742. |
[9] | Niu Y T, Chen J F, Hao H Z, et al. 2009. Observation and analysis of total solar eclipse of 1997 March 9 using VLF signal [J]. Journal of Henan Normal University (Natural Science) (in Chinese), 37(2): 64-66, 69. |
[10] | Schoute-Vanneck C A. 1974. VLF radio transmissions at sunrise [J]. Journal of Atmospheric and Terrestrial Physics, 36(3): 489-499. |
[11] | Shen C S, Zi M J, Wang J S, et al. 2005. Features of the polar ionosphere over the Zhongshan station in Antarctic [J]. Chinese J. Geophys. (in Chinese), 48(1): 1-6. |
[12] | Tschu K K, Wu Y Y, Du A M, et al. 2009. Review on solar-terrestrial activity and associated indices[J]. Progress in Geophysics (in Chinese), 24(4): 1167-1175, doi: 10. 3969/j. issn. 1004-2903. 2009. 04. 002. |
[13] | Wang H, Wang J X, Bai G X, et al. 2011. Real-time monitoring of ionosphere changes in the Yangtze River Delta region based on GPS technology during the total solar eclipse of 22 July 2009[J]. Chinese J. Geophys. (in Chinese), 2011, 54(7): 1718-1726, doi: 10. 3969/j. issn. 0001-5733. 2011. 07. 004. |
[14] | Wang J, Chen E P, Zhao X Z. 2006. Observation and study on phase of very low frequency propagation over polar area [J]. Chinese Journal of Radio Science (in Chinese), 21(3): 386-390. |
[15] | Wang S, Wei F S. 2007. Progress of space weather research in China [J]. Progress in Geophysics (in Chinese), 22(4): 1025-1029. |
[16] | Wang S Q, Pan S D, Ma L X. 1988. Study on the correlation between the solar X-ray burst and the VLF phase anomaly[J]. Publications of the Shaanxi Astronomical Observatory (in Chinese), 11(Z1): 52-59. |
[17] | Xu Z H, Liu R Y, Liu S L, et al. 2006. Variations of the ionospheric F2 layer critical frequency at Zhongshan station, Antarctica [J]. Chinese J. Geophys. (in Chinese), 49(1): 1-8. |
[18] | 贺龙松, 刘瑞源, 刘顺林,等. 2000. 太阳活动低年南极中山站电离层F层的平均特性[J]. 地球物理学报, 43(3): 289-295. |
[19] | 孔庆颜, 冯静, 柳文,等. 2010. 电离层参数短期预报研究进展[J]. 地球物理学进展, 25(6): 1968-1976, doi: 10. 3969/j. issn. 1004-2903. 2010. 06. 011. |
[20] | 李强, 宁百齐, 赵必强,等. 2012. 基于陆态网络GPS数据的电离层空间天气监测与研究[J]. 地球物理学报, 55(7): 2193-2202, doi: 10. 6038/j. issn. 0001-5733. 2012. 07. 006. |
[21] | 栗志, 杨晓宁, 周卫红. 1999. 对太阳活动区的自动监测、资料处理及研究的一些新方法[J]. 地球物理学进展, 14(S1): 135-141. |
[22] | 刘万通. 1987. VLF相位突然异常与太阳X射线事件相关性的分析[J]. 空间科学学报, 7(3): 20-24. |
[23] | 马宇倩, 王焕玉, 张承模,等. 2004. 2001年4月2日太阳耀斑及其太阳质子事件的观测结果研究[J]. 地球物理学报, 47(5): 737-742. |
[24] | 牛有田, 陈建锋, 郝好贞,等. 2009. 1997-03-09日全食的甚低频观测与分析[J]. 河南师范大学学报(自然科学版), 37(2): 64-66, 69. |
[25] | 沈长寿, 资民筠, 王劲松,等. 2005. 南极中山站电离层的极区特征[J]. 地球物理学报, 48(1): 1-6. |
[26] | 朱岗崑, 吴迎燕, 杜爱民,等. 2009. 关于日地关系主要元素与各种指数——综述及展望[J]. 地球物理学进展, 24(4): 1167-1175, doi: 10. 3969/j. issn. 1004-2903. 2009. 04. 002. |
[27] | 王虎, 王解先, 白贵霞,等. 2011. 基于GPS技术实时监测2009年7月22日日全食期间长三角区域电离层TEC变化[J]. 地球物理学报, 54(7): 1718-1726, doi: 10. 3969/j. issn. 0001-5733. 2011. 07. 004. |
[28] | 王健, 陈恩平, 赵协中. 2006. 甚低频横跨极区传播相位的观测与研究[J]. 电波科学学报, 21(3): 386-390. |
[29] | 王水, 魏奉思. 2007. 中国空间天气研究进展[J]. 地球物理学进展, 22(4): 1025-1029. |
[30] | 王素琴, 潘铄德, 麻兰香. 1988. 太阳X射线爆发与甚低频传播相位异常关系的研究[J]. 陕西天文台台刊, 11(Z1): 52-59. |
[31] | 徐中华, 刘瑞源, 刘顺林,等. 2006. 南极中山站电离层F2层临界频率变化特征[J]. 地球物理学报, 49(1): 1-8. |