基于修正Kp指数的中纬区单站电离层暴预报方法研究(2)

来源:地球物理学进展 作者:高 琴1,2, 刘立 发表于:2010-01-13 20:58  点击:
【关健词】犓p,中纬,磁暴,电离层暴,预测,STORM
结果作对比,我们主要检验了模型在夏季和春秋季 条件下对电离层扰动的预报效果.值得一提的是,在 犓f 的计算中,我们采用犓p 2 的27 天滑动均值犓qm 作为地磁指数的背景值,而在应用到实际预测时,由 于当

结果作对比,我们主要检验了模型在夏季和春秋季
条件下对电离层扰动的预报效果.值得一提的是,在
犓f 的计算中,我们采用犓p
2 的27 天滑动均值犓qm
作为地磁指数的背景值,而在应用到实际预测时,由
于当前日期后13 天的犓p 值未知,犓qm 的计算应该
相应地作些修改.在这里,我们只取当前日期前13
天的犓2
p 均值作为犓qm.
磁暴发生时间分别为1991年的4月4~7日、5
月13~16日、6月17~20日、7月8~11日、8月12
~15日和9月9~12日.各个磁暴的犓p 变化情况
见图3.本文模型及STORM 对这6 次磁暴下电离
层扰动的预报结果见图4.
与STORM 相比,本文模型有一个共同特点:
所预测的δ犳狅犉2 比STORM 先达到负最小值,且从
负最小值恢复到平静状态的速度更快.5 月13~16
日和7月8~11日的两次磁暴比较特殊,犓p 都有连
续的两个峰值(见图3).这反映的是两次磁层能量
注入过程,通常电离层响应会出现叠加效应[28].在
图4 中可看到,δ犳狅犉2 观测值有分两个阶段下降的
现象,而本文模型也预测到两次连续的负相扰动,与
观测值较接近,但STORM 模型只反映出一次
δ犳狅犉2 下降.值得一提的是,7月8~11日磁暴中犓p
峰值的时间间隔比5 月13~16 日的短,因此对应
δ犳狅犉2 分段下降的时间间隔也相应不同,而本文模
型很好地表现了这一差别.
图3 1991年6次磁暴的犓p 变化情况
Fig.3 Geomagneticactivities(犓p)for6stormsin1991
1946
 6期高 琴,等:基于修正犓p 指数的中纬区单站电离层暴预报方法研究
图4 1991年6次磁暴对应Yakutsk站的δ犳狅犉2 实际观测值及预测值
灰线表示观测值,黑色实线表示本文模型预测值,黑色点划线表示STORM 预测值
Fig.4 DataandpredictionofourmodelandSTORMatYakutskforthe6stormsin1991
Thegraylineistheobservation,theblacksolidlineshowstheoutputofourmodel
andtheblackdashedlineistheoutputofSTORM.
  对于5月、6月、7月及8月的4次暴,Yakutsk
站都表现为先发生幅度较大且持续时间较长的电离
层正暴,然后发生大的负暴.从图中可看到,
STORM 完全没有预报正暴,而我们的模型都对正
暴作出了反映.虽然所得到的正暴幅度还不够理想,
但效果比STORM 好许多,尤其是对正扰动持续时
间的预报与实际值很吻合.
下面,我们对两个模型的预报效果作一下定量
的比较与分析.采用绝对误差作为检验标准,其计算
式如下
RMSE = δ犳狅犉2mod -δ犳狅犉2obs , (4)
其中δ犳狅犉2mod为模型预测值,δ犳狅犉2obs为实际值.根
据这6 次磁暴样本,我们分别计算了本文模型及
STOMR 模型的预报值相对于实际值的平均绝对偏
差,如图5所示.从图中可看到,在δ犳狅犉2<40%范
图5 Yakutsk站模型预报绝对误差分析
Fig.5 Errorsofmodelprediction
andSTORMatYakutsk
1947
地 球 物 理 学 进 展24卷 
围内,本文模型的绝对误差明显比STORM 模型
小;在δ犳狅犉2>40%情况下,两个模型的误差存在较
大的起伏,可能是由于样本中满足条件的数据太少
的缘故. AraujoPradere and FullerRowell
(2002)[13]提出了一个衡量模型改善程度的因子:
 %improvement=
 
(RMSE(STORM)-RMSE(model))
RMSE(STORM)×100
. (5)
我们把两个模型在δ犳狅犉2<40%情况下的误差单独
作比较分析,得到本文模型的预测效果相对
STORM 提高了约35%.
此外,为验正模型方法的可靠性,我们改变训练
样本的时间范围.分别取1957~1980 年、1957~
1970年和1957~1960 年这三种情况,并将这三种
情况的计算结果与之前的1957~1990年作对比.所
得结果如图6所示.可以看到,除了1957~1960 年
的情况,样本范围取1957~1970 或1957~1980 年
时所得到的犪和犫与1957~1990年数据所得的计
算结果没有很大的偏离.这说明,只要所取的训练样
本时间范围足够长(如一个太阳周期),所得的模型
系数都能用于下一步电离层暴预报.
图6 不同训练样本下的模型系数a和b
Fig.6 Comparisonofaandbfordifferenttrainingsamples
  综上所述,本文提出的模型方法可靠性较高,建
立电离层暴预报模型在中纬Yakutsk 站的应用效
果比STORM 模型好.
进一步地,我们对更多的中纬区电离层台站作
了相同的分析与检验,发现本文模型方法在许多台
站都得到了比STORM 更好的预报效果,但同时也
有许多台站的预报效果并不理想.由于篇幅限制,这
些台站的分析结果在这里就不一一列出.对于模型
并不能适用于所有中纬区台站,我们认为其主要原
因在于,本文模型方法假设电离层暴主要受地磁犓p
变化的影响,假设电离层暴的季节性效应占主导地
位,而由于电离层暴的复杂性,这两个假设在某些台
站并不能成立,因而导致模型的预报效果不够好.
3 小 结
本文主要借鉴文献[9]中构建修正犓p 指数犓f
的方法,统计分析中纬Yakutsk 站δ犳狅犉2 与犓f 的
关系发现两者有非常好的线性关系,且季节性特征
明显,最终建立了一个简易的中纬区单站电离层暴
预报模型,其模型系数只随月份不同.
1948
 6期高 琴,等:基于修正犓p 指数的中纬区单站电离层暴预报方法研究
模型检验表明,本文模型在Yakutsk 站达到了
比STORM 模型更好的预报效果.模型的平均绝对
误差明显比STORM 的小,对δ犳狅犉2<40% 的扰动
的预测效果也比STORM 平均提高了35%左右.特
别是模型能很好地反映磁暴连续发生时对应δ犳狅犉2
分阶段下降的形态.与STORM 不能预测电离层正
相扰动相比,本文模型在这方面改善了很多,可以对
正相扰动的幅度和持续时间作出一定的反映.
文中还用该模型在其它中纬台站作了检验,预
报效果有好有差.虽然并不能在所有台站都达到比
STORM 更好的预报效果,但由于模型系数是用单(责任编辑:南粤论文中心)转贴于南粤论文中心: http://www.nylw.net(南粤论文中心__代写代发论文_毕业论文带写_广州职称论文代发_广州论文网)
顶一下
(0)
0%
踩一下
(0)
0%


版权声明:因本文均来自于网络,如果有版权方面侵犯,请及时联系本站删除.