留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

中国区域掩星观测与IRI-2016电离层峰值参数的比较

孙方方 罗佳 徐晓华 王涵

孙方方, 罗佳, 徐晓华, 王涵. 中国区域掩星观测与IRI-2016电离层峰值参数的比较[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
引用本文: 孙方方, 罗佳, 徐晓华, 王涵. 中国区域掩星观测与IRI-2016电离层峰值参数的比较[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
SUN Fangfang, LUO Jia, XU Xiaohua, WANG Han. Comparisons of Ionospheric Peak Parameters from Radio Occultation Observations and IRI-2016 Model Outputs over China[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
Citation: SUN Fangfang, LUO Jia, XU Xiaohua, WANG Han. Comparisons of Ionospheric Peak Parameters from Radio Occultation Observations and IRI-2016 Model Outputs over China[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391

中国区域掩星观测与IRI-2016电离层峰值参数的比较

doi: 10.13203/j.whugis20180391
基金项目: 

国家自然科学基金 41774032

国家自然科学基金 41774033

国家重点研发计划 2018YFC1503502

国家重点基础研究发展计划(973计划) 2013CB733302

详细信息

Comparisons of Ionospheric Peak Parameters from Radio Occultation Observations and IRI-2016 Model Outputs over China

Funds: 

The National Natural Science Foundation of China 41774032

The National Natural Science Foundation of China 41774033

the National Key Research and Development Program of China 2018YFC1503502

the National Key Basic Research and Development Program of China 2013CB733302

More Information
    Author Bio:

    SUN Fangfang, postgraduate, specializes in GNSS ionospheric radio occultation. E-mail:sunfangfang@whu.edu.cn

    Corresponding author: LUO Jia, PhD, associatet professor. E-mail: jialuo@whu.edu.cn
  • 摘要: 对2008—2014年由气象、电离层和气候卫星联合观测系统(Constellation Observing System for Meteorology Ionosphere and Climate,COSMIC)掩星观测与最新版国际电离层参考模型IRI-2016输出得到的电离层峰值参数(峰值密度NmF2与峰值高度hmF2)在中国区域进行了比较。IRI-2016模型输出值与COSMIC掩星反演值的相关性在太阳活动高年(2011—2014年)整体上高于太阳活动低年(2008—2010年)。在低年春秋季的当地时间(local time,LT)12:00—14:00,IRI-2016相对于COSMIC掩星在30°N~55°N区域内对NmF2和hmF2分别存在低估和高估现象,在15°N~30°N区域内则恰恰相反。对于NmF2,采用IRI_CCIR和IRI_URSI两种选项的模型输出值在中午时分均存在高估,在低年高估更为显著。对于hmF2,采用IRI_CCIR和IRI_AMTB两种选项的模型输出值在低年各季节均存在高估,且IRI_AMTB选项高估更显著,冬季最突出。结果表明,在中国区域由IRI-2016模型计算NmF2和hmF2时,分别推荐使用IRI_CCIR和IRI_Shubin选项。
  • 图  1  COSMIC反演的电子密度廓线和IRI-2016模型不同选项输出值的比较

    Figure  1.  Comparison Among Electron Density Profiles Derived from COSMIC RO and IRI-2016 Model with Different Model Options

    图  2  2008—2014年F10.7指数时间序列

    Figure  2.  Time Series of F10.7 Index During 2008—2014

    图  3  COSMIC掩星观测与IRI-2016模型采用推荐选项输出的NmF2与hmF2散点图

    Figure  3.  Scatter Plots of NmF2 and hmF2 from COSMIC RO and IRI-2016 Model with Recommended Model Options

    图  4  2008—2010年春秋季LT 12:00—14:00 COSMIC掩星与IRI-2016模型的NmF2与hmF2在中国区域的空间分布

    Figure  4.  Spatial Distributions of NmF2 and hmF2 over China Region at LT 12:00—14:00 During Equinox Months of 2008—2010 Derived from COSMIC RO and IRI-2016 Model

    图  5  不同季节COSMIC掩星与IRI-2016模型两种不同选项下NmF2的日变化图

    Figure  5.  Diurnal Variations of NmF2 from COSMIC RO and IRI-2016 Model Using Two Different Options for Different Seasons

    图  6  不同季节COSMIC掩星与IRI-2016模型不同选项下hmF2的日变化图

    Figure  6.  Diurnal Variations of hmF2 from COSMIC RO and IRI-2016 Model Using Different Options for Different Seasons

    表  1  2008—2014年NmF2和hmF2的MARE/%

    Table  1.   MARE of NmF2 and hmF2 During 2008―2014/%

    统计项 2008年 2009年 2010年 2011年 2012年 2013年 2014年
    NmF2 31.112 2 31.481 8 27.836 7 26.7882 21.888 1 24.443 1 25.204 4
    hmF2 11.689 0 11.944 2 9.747 9 8.460 8 7.071 9 8.021 6 7.344 6
    下载: 导出CSV
  • [1] Yue X, Schreiner W S, Rocken C, et al. Validate the IRI2007 Model by the COSMIC Slant TEC Data During the Extremely Solar Minimum of 2008[J].Advances in Space Research, 2013, 51(4):647-653 doi:  10.1016/j.asr.2011.08.011
    [2] Gordiyenko G I, Yakovets A F. Comparison of Midlatitude Ionospheric F Region Peak Parameters and Topside Ne Profiles from IRI2012 Model Prediction with Ground-Based Ionosonde and Alouette II Observations[J]. Advances in Space Research, 2017, 60(2):461-474 doi:  10.1016/j.asr.2017.01.006
    [3] Bilitza D, Altadill D, Zhang Y, et al. The International Reference Ionosphere 2012-a Model of International Collaboration[J]. Journal of Space Weather and Space Climate, 2014, 4(A07):1-12 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Doaj000003666511
    [4] Bilitza D, Altadill D, Truhlik V, et al. International Reference Ionosphere 2016:From Ionospheric Climate to Real-Time Weather Predictions[J]. Space Weather, 2017, 15:418-429 doi:  10.1002/2016SW001593
    [5] Altadill D, Magdaleno S, Torta J M, et al. Global Empirical Models of the Density Peak Height and of the Equivalent Scale Height for Quiet Conditions[J]. Advances in Space Research, 2013, 52(10):1756-1769 doi:  10.1016/j.asr.2012.11.018
    [6] Shubin V N. Global Median Model of the F2-Layer Peak Height Based on Ionospheric Radio-Occultation and Ground-Based Digisonde Observations[J]. Advances in Space Research, 2015, 56(5):916-928 doi:  10.1016/j.asr.2015.05.029
    [7] Lei J, Syndergaard S, Burns A G, et al. Comparison of COSMIC Ionospheric Measurements with Ground-Based Observations and Model Predictions:Preliminary Results[J]. Journal of Geophysical Research:Space Physics, 2007, 112(47):A07308 doi:  10.1029-2006JA012240/
    [8] Chu Y H, Su C L, Ko H T. A Global Survey of COSMIC Ionospheric Peak Electron Density and Its Height:A Comparison with Ground-Based Ionosonde Measurements[J]. Advances in Space Research, 2010, 46(4):431-439 doi:  10.1016/j.asr.2009.10.014
    [9] 孙凌峰, 赵必强, 乐新安, 等.中国区域电离层垂测仪探测参量与COSMIC掩星反演结果比较研究[J].地球物理学报, 2014, 57(11):3625-3632 http://d.old.wanfangdata.com.cn/Conference/8495752

    Sun L F, Zhao B Q, Yue X A, et al. Comparison Between Ionospheric Character Parameters Retrieved from FORMOSAT3 Measurements and Ionosonde Observation over China[J]. Chinese Journal of Geophysics, 2014, 57(11):3625-3632 http://d.old.wanfangdata.com.cn/Conference/8495752
    [10] 赵莹, 张小红. COSMIC掩星观测数据反演电离层电子密度廓线[J].武汉大学学报·信息科学版, 2010, 35(6):644-648 http://ch.whu.edu.cn/CN/abstract/abstract980.shtml

    Zhao Ying, Zhang Xiaohong. Inversion of Ionospheric Electron Density Profiles with COSMIC Occultation Data[J].Geomatics and Information Science of Wuhan University, 2010, 35(6):644-648 http://ch.whu.edu.cn/CN/abstract/abstract980.shtml
    [11] Zhao X K, Ning B Q, Zhang M L, et al. Comparison of the Ionospheric F-2 Peak Height Between Ionosonde Measurements and IRI2016 Predictions over China[J]. Advances in Space Research, 2017, 60(7):1524-1531 doi:  10.1016/j.asr.2017.06.056
    [12] Yang K F, Chu Y H, Su C L, et al. An Examination of FORMOSAT-3/COSMIC Ionospheric Electron Density Profile:Data Quality Criteria and Comparisons with the IRI Model[J]. Terrestrial, Atmospheric and Oceanic Sciences, 2009, 20(1):193-206 doi:  10.3319/TAO.2007.10.05.01(F3C)
    [13] 黄智, 袁洪.磁赤道地区2007-2013年COSMIC掩星反演和国际参考电离层模型输出结果分析[J].地球物理学报, 2016, 59(7):2333-2343 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201607001

    Huang Zhi, Yuan Hong. An Analysis of Ionospheric Characteristic Parameters Retrieved from COSMIC and IRI at Jicamarca During the Period of 2007-2013[J]. Chinese Journal of Geophysics, 2016, 59(7):2333-2343 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb201607001
    [14] Yue X, Schreiner W S, Lei J, et al. Error Analysis of Abel Retrieved Electron Density Profiles from Radio Occultation Measurements[J]. Annales Geophysicae, 2010, 28(1):217-222 doi:  10.5194/angeo-28-217-2010
    [15] Hu L, Ning B, Liu L, et al. Validation of COSMIC Ionospheric Peak Parameters by the Measurements of an Ionosonde Chain in China[J]. Annales Geophysicae, 2014, 32(10):1311-1319 doi:  10.5194/angeo-32-1311-2014
  • [1] 陈伟伟, 熊永良, 徐韶光, 高姝妹.  中国区域COSMIC数据掩星反射信号特征分析 . 武汉大学学报 ● 信息科学版, 2022, 47(2): 189-196. doi: 10.13203/j.whugis20190475
    [2] 王耀兴, 张秋昭, 沈震.  COSMIC掩星反演湿温廓线质量的时空分析 . 武汉大学学报 ● 信息科学版, 2021, 46(6): 887-894. doi: 10.13203/j.whugis20190051
    [3] 徐晓华, 朱洲宗, 罗佳.  利用IGRA2探空数据和COSMIC掩星资料对FY-3C掩星中性大气产品进行质量分析 . 武汉大学学报 ● 信息科学版, 2020, 45(3): 384-393. doi: 10.13203/j.whugis20180490
    [4] 徐晓华, 刘树纶, 罗佳.  利用COSMIC掩星折射指数分析全球大气边界层顶结构变化 . 武汉大学学报 ● 信息科学版, 2018, 43(1): 94-100. doi: 10.13203/j.whugis20160183
    [5] 徐晓华, 郭金城, 罗佳.  利用COSMIC掩星资料分析大气重力波参数的全球分布特征 . 武汉大学学报 ● 信息科学版, 2015, 40(11): 1493-1498. doi: 10.13203/j.whugis20130587
    [6] 屈小川, 安家春, 刘根.  利用COSMIC掩星资料分析南极地区对流层顶变化 . 武汉大学学报 ● 信息科学版, 2014, 39(5): 605-610. doi: 10.13203/j.whugis20130728
    [7] 瞿伟, 张勤, 李振洪, 王庆良.  山西清徐地裂缝构造活动参数反演 . 武汉大学学报 ● 信息科学版, 2013, 38(4): 421-425.
    [8] 张锦明, 游雄, 万刚.  径向基函数算法中插值参数对DEM精度的影响 . 武汉大学学报 ● 信息科学版, 2013, 38(5): 608-612.
    [9] 夏朋飞, 蔡昌盛, 戴吾蛟, 陈必焰.  地基GPS联合COSMIC掩星数据的水汽三维层析研究 . 武汉大学学报 ● 信息科学版, 2013, 38(8): 892-896.
    [10] 肖卫华, 符养, 杜晓勇, 高太长.  利用COSMIC掩星资料计算平流层大气风场方法分析 . 武汉大学学报 ● 信息科学版, 2012, 37(2): 199-204.
    [11] 许超钤, 史俊波, 郭际明, 徐晓华.  联合地基GPS和空基COSMIC掩星的可降水量研究 . 武汉大学学报 ● 信息科学版, 2011, 36(4): 467-470.
    [12] 赵莹, 张小红.  COSMIC掩星观测数据反演电离层电子密度廓线 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 644-648.
    [13] 徐晓华, ZHANG Kefei, 罗佳, FU Erjiang.  GPS掩星与探空观测统计比较中配对标准的比较研究 . 武汉大学学报 ● 信息科学版, 2009, 34(11): 1332-1335.
    [14] 徐晓华, 罗佳.  COSMIC掩星折射指数廓线的统计验证 . 武汉大学学报 ● 信息科学版, 2009, 34(2): 214-217.
    [15] 范磊, 符养, 杜晓勇, 李黄.  雾灵山山基掩星观测反演误差分析 . 武汉大学学报 ● 信息科学版, 2008, 33(1): 89-92.
    [16] 徐晓华, ZhangKefei, FuErjiang, 罗佳.  澳大利亚上空COSMIC掩星廓线的反演 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 800-804.
    [17] 徐晓华, 李征航, 罗佳.  GPS掩星数据反演中的地球扁率影响改正 . 武汉大学学报 ● 信息科学版, 2005, 30(6): 502-505.
    [18] 徐晓华, 李征航, 罗佳.  单颗LEO卫星轨道参数对GPS掩星事件分布和数量影响的模拟研究 . 武汉大学学报 ● 信息科学版, 2005, 30(7): 609-612.
    [19] 徐晓华, 李征航, 罗佳.  利用GPS/MET与CHAMP掩星资料反演地球中性大气参数的比较 . 武汉大学学报 ● 信息科学版, 2004, 29(3): 214-217.
    [20] 徐晓华, 李征航, 罗佳.  利用GPS掩星资料反演地球中性大气参数折射角方法研究 . 武汉大学学报 ● 信息科学版, 2003, 28(5): 589-592.
  • 加载中
图(6) / 表(1)
计量
  • 文章访问数:  1253
  • HTML全文浏览量:  163
  • PDF下载量:  82
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-30
  • 刊出日期:  2020-03-05

中国区域掩星观测与IRI-2016电离层峰值参数的比较

doi: 10.13203/j.whugis20180391
    基金项目:

    国家自然科学基金 41774032

    国家自然科学基金 41774033

    国家重点研发计划 2018YFC1503502

    国家重点基础研究发展计划(973计划) 2013CB733302

    作者简介:

    孙方方, 硕士生, 主要从事GNSS掩星电离层研究。sunfangfang@whu.edu.cn

    通讯作者: 罗佳, 博士, 副教授。jialuo@whu.edu.cn
  • 中图分类号: P228

摘要: 对2008—2014年由气象、电离层和气候卫星联合观测系统(Constellation Observing System for Meteorology Ionosphere and Climate,COSMIC)掩星观测与最新版国际电离层参考模型IRI-2016输出得到的电离层峰值参数(峰值密度NmF2与峰值高度hmF2)在中国区域进行了比较。IRI-2016模型输出值与COSMIC掩星反演值的相关性在太阳活动高年(2011—2014年)整体上高于太阳活动低年(2008—2010年)。在低年春秋季的当地时间(local time,LT)12:00—14:00,IRI-2016相对于COSMIC掩星在30°N~55°N区域内对NmF2和hmF2分别存在低估和高估现象,在15°N~30°N区域内则恰恰相反。对于NmF2,采用IRI_CCIR和IRI_URSI两种选项的模型输出值在中午时分均存在高估,在低年高估更为显著。对于hmF2,采用IRI_CCIR和IRI_AMTB两种选项的模型输出值在低年各季节均存在高估,且IRI_AMTB选项高估更显著,冬季最突出。结果表明,在中国区域由IRI-2016模型计算NmF2和hmF2时,分别推荐使用IRI_CCIR和IRI_Shubin选项。

English Abstract

孙方方, 罗佳, 徐晓华, 王涵. 中国区域掩星观测与IRI-2016电离层峰值参数的比较[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
引用本文: 孙方方, 罗佳, 徐晓华, 王涵. 中国区域掩星观测与IRI-2016电离层峰值参数的比较[J]. 武汉大学学报 ● 信息科学版, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
SUN Fangfang, LUO Jia, XU Xiaohua, WANG Han. Comparisons of Ionospheric Peak Parameters from Radio Occultation Observations and IRI-2016 Model Outputs over China[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
Citation: SUN Fangfang, LUO Jia, XU Xiaohua, WANG Han. Comparisons of Ionospheric Peak Parameters from Radio Occultation Observations and IRI-2016 Model Outputs over China[J]. Geomatics and Information Science of Wuhan University, 2020, 45(3): 403-410. doi: 10.13203/j.whugis20180391
  • 作为最为广泛使用的电离层经验模型之一,国际参考电离层模型(international reference ionosphere, IRI)对电离层研究和应用具有十分重要的意义[1]。IRI模型利用了垂测仪、非相干散射雷达、卫星观测及探空火箭数据,可以输出任意指定时间和经纬度的60~2 000 km的电子密度和总电子含量(total electron content, TEC)等多种电离层参数信息[2]。自1972年以来,随着电离层观测数据的积累、观测技术的革新和建模方法的改进,IRI模型已发布了多个版本,目前最新版本为IRI-2016 [3-4]

    峰值密度NmF2和峰值高度hmF2是表征电子密度廓线形态的重要参量。对于NmF2的计算,IRI-2016模型中包含IRI_CCIR和IRI_URSI两种模型选项。对于hmF2的计算,在旧版本模型使用的IRI_CCIR hmF2模型的基础上,IRI-2016模型中另外新增了IRI_AMTB和IRI_Shubin两种模型选项。hmF2的IRI_CCIR模型是基于1954—1958年全球分布的约150个垂测仪观测资料,利用hmF2与参数M(3 000)F2的强相关关系由M(3 000)F2计算得到。其中,M(3 000)F2为3 000 km高度处最大可用频率和F2层临界频率(foF2)的比值。这种间接计算hmF2的方法有诸多缺陷。IRI-2016模型中新增的两种计算hmF2的模型选项直接对hmF2进行建模,不依赖M(3 000)F2。其中,IRI_AMTB模型是根据全球分布的26个垂测仪(其中包括中国区域的武汉站、安阳站和三亚站)在1998—2006年的hmF2观测值,纯粹利用地基垂测仪观测值进行建模[5];IRI_Shubin模型的建立则使用了挑战性小卫星有效载荷(Challenging Minisatellite Payload,CHAMP)、重力恢复与气候实验(Gravity Recovery and Climate Experiment,GRACE)和气象、电离层和气候卫星联合观测系统(Constellation Observing System for Meteorology, Ionosphere and Climate,COSMIC)3个掩星任务的观测资料以及分布全球的62个垂测仪于1987—2012年的观测数据[6]

    COSMIC由美国国防部和中国台湾地区联合研制并于2006-04-15发射升空,在掩星观测任务中具有里程碑的意义[7]。COSMIC由6颗低轨卫星组成,轨道高度在800 km左右,轨道倾角为72°,并配备能够接收GPS信号的接收机。掩星事件发生时,GPS信号穿过电离层时受电子密度影响而发生折射,产生附加相位延迟。在信号直线传播和球对称假设等前提下,由传播路径上的总电子含量使用Abel积分变换方法可得到电子密度廓线。COSMIC掩星反演结果已经得到了多种独立观测手段的验证[8-10]

    中国区域(70°E~140°E,15°N~55°N)跨越中低纬地区,电离层变化幅度较大。文献[11]利用中国区域的4个垂测仪台站对IRI-2016模型在中国区域的性能进行了评估,但只能对垂测仪台站周围的电离层特征进行验证分析。COSMIC掩星观测具有全球覆盖、高精度、垂直分辨率高等特点,可用于对IRI-2016模型在整个中国区域的表现进行分析。本文的研究内容包括使用2008—2014年的COSMIC掩星数据对地磁平静条件下IRI-2016模型在中国区域的性能进行验证,并对IRI-2016模型在中国区域使用时不同NmF2和hmF2模型选项的输出结果进行统计分析,以确定IRI-2016模型在中国区域输出值的最优选项。

    • COSMIC掩星二级产品“ionPrf”数据可从COSMIC数据存档与分析中心(COSMIC Data Analysis and Archive Center, CDAAC)的网站http://cdaac-www.cosmic.ucar.edu下载得到。该数据包含了掩星事件发生的时间、位置、NmF2、hmF2和电子密度廓线等信息。在分析地磁平静条件下IRI模型的表现时,需要先剔除3 h ap指数大于15的时间段的COSMIC掩星数据[8]。另一方面,尽管CDAAC在反演电子密度廓线时采取了一定的质量检核措施,但通过其检核的产品中仍然存在问题廓线[8]。文献[12]提出根据顶部电离层电子密度廓线的斜率和电子密度波动的平均偏差(mean deviation, MD)来剔除这些问题廓线。本文选择COSMIC电子密度廓线的标准为在420~490 km的斜率小于-0.02;平均偏差不超过1.5。此外,峰值高度hmF2范围应在200~500 km之间。

    • IRI-2016模型输出值可从http://irimodel.org/网站获取。该模型中计算NmF2的选项包括IRI_CCIR模型和IRI_URSI模型。其中,IRI_URSI模型为推荐选项。计算峰值高度hmF2的选项包括原有的IRI_CCIR模型以及新增的IRI_AMTB模型和IRI_Shubin模型。其中,IRI_AMTB模型为推荐选项。模型的输入参数包括经纬度和时间。由于掩星观测具有水平漂移,本文选用COSMIC掩星观测的F2层峰值高度处的经纬度和发生时刻作为输入量,通过控制模型选项得到不同选项下输出的NmF2和hmF2。图 1(a)1(b)分别为当地时间(local time, LT)2018-02-03T13:35和2008-02-08T23:46的COSMIC掩星电子密度廓线和IRI-2016模型不同选项下输出结果的比较图。可以看出,COSMIC电子密度廓线和IRI-2016模型不同选项以及选项之间的输出结果存在着明显差异。

      图  1  COSMIC反演的电子密度廓线和IRI-2016模型不同选项输出值的比较

      Figure 1.  Comparison Among Electron Density Profiles Derived from COSMIC RO and IRI-2016 Model with Different Model Options

    • F10.7指数是指10.7 cm(2 800 MHz)波长的太阳辐射通量,单位为sfu,与太阳黑子数之间有强相关性,能很好地表征太阳活动强度。图 2是F10.7指数的月均值在第24太阳活动周2008—2014年期间的变化。可以看出,2008—2010年为太阳活动低年,尤其是2008、2009年,太阳活动达到了极低的水平;之后太阳活动水平逐渐提高,在2011—2014年太阳活动维持在较高水平,为太阳活动高年。

      图  2  2008—2014年F10.7指数时间序列

      Figure 2.  Time Series of F10.7 Index During 2008—2014

    • 为了了解IRI-2016模型在中国区域的表现,首先采用模型推荐的选项计算NmF2和hmF2,将结果与COSMIC掩星观测进行比较。图 3为2008年和2014年COSMIC掩星观测与IRI-2016模型输出的峰值密度NmF2与峰值高度hmF2的散点图。其中,红色实线代表COSMIC观测值等于IRI预测值,黑色虚线代表使用最小二乘方法得到的拟合线,RS分别表示两种方法得到的相关系数和拟合线斜率。由图 3(a)3(b)可以看出,尽管在太阳活动低年2008年和太阳活动高年2014年,COSMIC掩星观测和IRI模型输出的NmF2相关系数非常接近(分别为0.86和0.87),但在2014年数据对较为均匀地围绕红色实线分布,且拟合线和等值线非常接近,斜率高达0.90。而在2008年,数据对较为分散,且较多的数据对分布在等值线下侧,拟合直线与等值线之间存在较大偏差。拟合直线的斜率(0.70)表明在太阳活动低年2008年,相较于COSMIC掩星观测值,IRI-2016模型明显高估了NmF2。与峰值密度NmF2情形类似,图 3(c)图 3(d)表明, 在2014年IRI-2016模型峰值高度hmF2输出值与COSMIC掩星观测值更为接近,拟合直线几乎与等值线重合。而在2008年,拟合直线的斜率为0.76,远小于2014年的斜率0.97,说明2008年IRI-2016模型的hmF2输出值相对于COSMIC掩星观测存在明显高估。同时,由相关系数可以看出,无论在2008年还是2014年,NmF2的相关系数均高于hmF2,表明IRI-2016模型和COSMIC掩星观测输出的NmF2一致性优于hmF2。

      图  3  COSMIC掩星观测与IRI-2016模型采用推荐选项输出的NmF2与hmF2散点图

      Figure 3.  Scatter Plots of NmF2 and hmF2 from COSMIC RO and IRI-2016 Model with Recommended Model Options

      为了进一步定量分析太阳活动低年和高年COSMIC掩星观测值和IRI-2016模型输出值的差异,定义了相对中误差绝对值的均值(mean absolute relative error,MARE)M,计算公式为:

      $$ M=\frac{\mathop{\sum }^{}\left( \left| I-C \right|/C \right)}{N} $$

      式中,C表示由掩星数据反演的NmF2或者hmF2;I表示对应的由模型计算的NmF2或者hmF2;N为观测数据个数。

      表 1给出了2008—2014年的NmF2和hmF2的MARE数值。可以看出,无论是NmF2还是hmF2,太阳活动低年(2008—2010年)的MARE均高于太阳活动高年(2011—2014年),在太阳活动极小年(2008—2009年),IRI-2016模型输出值与COSMIC掩星观测值的差异最大。

      表 1  2008—2014年NmF2和hmF2的MARE/%

      Table 1.  MARE of NmF2 and hmF2 During 2008―2014/%

      统计项 2008年 2009年 2010年 2011年 2012年 2013年 2014年
      NmF2 31.112 2 31.481 8 27.836 7 26.7882 21.888 1 24.443 1 25.204 4
      hmF2 11.689 0 11.944 2 9.747 9 8.460 8 7.071 9 8.021 6 7.344 6

      因此,总体而言,对中国区域,无论是NmF2还是hmF2,以COSMIC掩星数据为参照,IRI-2016模型在太阳活动高年的表现都优于太阳活动低年。

    • 为了进一步分析中国区域COSMIC掩星反演的峰值参数与IRI-2016模型输出值间的差异,对比了差异更为明显的太阳活动低年(2008—2010年)峰值参数的空间分布。考虑到电离层活跃水平,分析的时间段选为最具代表性的春秋季的中午时分(LT 12:00—14:00)。在本文中,季节的划分标准为:3月、4月、9月、10月为春秋季,5月—8月为夏季,11月、12月、1月和2月为冬季。图 4为2008—2010年春秋季中午(LT 12:00—14:00)由COSMIC掩星与IRI-2016模型得到的峰值参数NmF2与hmF2的均值在中国区域的空间分布。图 4(a)~4(d)显示,COSMIC掩星观测与IRI模型输出的峰值参数空间分布形态均反映出低纬地区电离层更为活跃,同时观测值与模型输出值之间存在明显差异。图 4(e)4(g) 表明,对于NmF2,与COSMIC掩星反演结果对比,IRI模型输出值在15°N~30°N的低纬区域明显高估,相对偏差最高可达100%以上;而在30°N~55°N的中纬区域,IRI模型输出的NmF2普遍表现为低估。图 4(f)4(h)表明,对于hmF2,情况恰恰相反:在30°N~55°N中纬区域,IRI模型的hmF2输出值普遍高于COSMIC观测值,而在15°N~30°N低纬区域,模型输出值普遍低于COSMIC观测值。由图 4(f)可见,在太阳活动低年,IRI-2016模型与COSMIC在中国区域得到的hmF2的绝对偏差一般不超过50 km。

      图  4  2008—2010年春秋季LT 12:00—14:00 COSMIC掩星与IRI-2016模型的NmF2与hmF2在中国区域的空间分布

      Figure 4.  Spatial Distributions of NmF2 and hmF2 over China Region at LT 12:00—14:00 During Equinox Months of 2008—2010 Derived from COSMIC RO and IRI-2016 Model

    • 以上分析中所使用的IRI-2016模型输出值均基于模型的推荐选项得到。除了推荐选项,IRI-2016对NmF2和hmF2的计算均提供其他模型选项,而使用不同选项的模型输出结果存在差异。为了了解不同选项的模型在中国区域的表现,进一步以COSMIC掩星数据为参照,对IRI-2016模型中两种计算NmF2的选项(IRI_CCIR和IRI_URSI模型选项)和3种计算hmF2的选项(IRI_CCIR、IRI_AMTB和IRI_Shubin模型选项)的输出结果进行了比较。图 5给出了不同季节COSMIC掩星与IRI-2016模型两种不同选项得到的NmF2均值的日变化情况,其中图 5(a)~5(c)图 5(d)~5(f)分别为太阳活动低年和高年的结果。统计中如果掩星事件发生时间与某一整点时刻相差不超过30 min,则该时刻的掩星数据和IRI模型输出结果被用于计算该整点时刻的平均值。可以看出,无论采用哪个选项,NmF2在约LT 04:00达到最小值,之后快速升高,在约LT 14:00达到全天最大值,之后逐渐下降。春秋季NmF2高于其他季节,电离层最为活跃。基于两种不同选项的IRI-2016模型输出值与COSMIC观测值均有相似的日变化趋势,但分别基于两种选项的模型输出结果之间存在一定差异。图 5(a)~5(c)表明采用两种选项的IRI-2016模型输出结果非常接近,但在中午时分,IRI_CCIR模型和IRI_URSI模型均高估了峰值密度,其中IRI_URSI的高估程度更严重。在春秋季和冬季,IRI_CCIR模型和IRI_URSI模型在LT 17:00—05:00则低估了峰值密度,其中IRI_URSI低估程度更严重。图 5(d)~5(f)图 5(a)~5(c)的对比表明,随着太阳活动水平的增加,NmF2随之增大,但中午时分采用两种选项的IRI-2016模型输出值与COSMIC观测值之间的差异均有减小。特别是在太阳活动高年的夏季,IRI-2016模型和COSMIC观测值全天都符合得很好。而在春秋季LT 16:00—04:00,两种选项的IRI-2016模型输出值均小于COSMIC观测值,其中IRI_URSI的偏离度相对更大。

      图  5  不同季节COSMIC掩星与IRI-2016模型两种不同选项下NmF2的日变化图

      Figure 5.  Diurnal Variations of NmF2 from COSMIC RO and IRI-2016 Model Using Two Different Options for Different Seasons

      图 6给出了不同季节COSMIC掩星观测与IRI-2016模型3种不同选项输出的hmF2日变化情况。可以看出,3种选项对应的IRI-2016模型输出与COSMIC观测的hmF2呈现相似的日变化趋势,即在LT 06:00―07:00达到最小值,在LT 12:00达到小高峰,并在约LT 24:00LT达到全天最大值。在太阳活动低年,3种计算hmF2的模型选项中,IRI_Shubin模型输出值和COSMIC观测值最为接近,其他两种模型输出值相对于COSMIC观测值均存在高估。IRI_AMTB模型输出值相对于COSMIC观测值的高估现象更为突出,冬季最显著。与太阳活动低年相比较,在太阳活动高年,3种选项对应的hmF2模型输出值与COSMIC观测值的一致性均相对更好。但IRI_AMTB模型输出值在夜间的偏差仍较大。

      图  6  不同季节COSMIC掩星与IRI-2016模型不同选项下hmF2的日变化图

      Figure 6.  Diurnal Variations of hmF2 from COSMIC RO and IRI-2016 Model Using Different Options for Different Seasons

    • 本文以COSMIC掩星数据为参照,对IRI-2016模型在中国区域的电离层峰值参数NmF2和hmF2输出值随太阳活动强度、季节和昼夜的变化进行了对比分析。整体而言,模型输出值与COSMIC掩星观测值的相关性较高,但太阳活动低年的相关性低于太阳活动高年。两组数据拟合直线的斜率表明,IRI-2016模型输出值与COSMIC掩星观测值的一致性在太阳活动低年低于太阳活动高年。

      文献[13]使用位于磁赤道附近的Jicamarca站垂测仪数据与IRI模型输出值进行验证分析,结果发现模型输出值与垂测仪观测值在太阳活动高年的相关性高于太阳活动低年,这与本文的研究结果一致。这主要是由于太阳活动低年,特别是2008年、2009年为太阳活动非寻常极小年,太阳活动水平极低且持续时间较长,IRI模型输出值会产生较大误差,甚至出现不合理的输出结果[3]

      太阳活动低年COSMIC掩星与IRI-2016模型峰值参数在中国区域的空间分布比较情况显示,虽然两者的空间分布形态基本一致,但仍存在明显差异。在15°N~30°N范围内,IRI-2016模型往往高估NmF2而低估hmF2。在30°N~55°N范围内,IRI-2016模型则存在低估NmF2而高估hmF2的现象。COSMIC掩星与IRI-2016模型之间的差异一方面是由于IRI-2016模型本身的误差,另一方面则是由于COSMIC掩星反演结果存在的误差。在球对称假设的前提下,COSMIC掩星反演的电子密度实际上是掩星切点附近整个掩星区域的平均值。因此,在电子密度梯度较大的区域,COSMIC掩星反演的电子密度与真实大气状态相比可能存在较大差异[14]。中国区域低纬度地区位于赤道电离层异常的驼峰区域附近,电子密度梯度很大,导致COSMIC反演的电子密度廓线相较于实际值可能存在较大偏差。

      文献[15]使用中国区域位于同一经度带上的4个垂测仪台站(漠河站、北京站、武汉站和三亚站)数据对COSMIC掩星峰值参数进行了比较验证,发现在漠河站两者符合程度较高,但在武汉站和北京站中午时分COSMIC NmF2观测值高于垂测仪观测值,在三亚站中午时分COSMIC NmF2观测值低于垂测仪观测值。反演存在的误差可能使得在中纬度的武汉站(31.0°N,114.6°E)COSMIC反演的NmF2高于IRI-2016模型值,在低纬度的三亚站(18.3°N,109.6°E)COSMIC反演的NmF2低于IRI-2016模型值。文献[11]使用这4个垂测仪台站数据对IRI-2016模型计算hmF2的不同选项输出值结果进行比较分析,发现IRI-2016模型推荐的IRI_AMTB模型选项在4个台站均高估hmF2。本文结果也表明IRI-2016模型的hmF2输出值在30°N~55°N高于COSMIC掩星观测值。对于15°N~30°N范围,IRI-2016模型的hmF2输出值相对于掩星资料存在低估,原因可能是由于该范围位于赤道电离层异常驼峰附近,掩星反演的hmF2高于实际值[11]

    • 本文对2008—2014年IRI-2016模型在中国区域的电离层峰值密度NmF2和峰值高度hmF2输出值与COSMIC掩星数据进行了对比,并对IRI-2016模型中计算这两个峰值参数的不同选项进行了比较分析。结论如下:

      1)IRI-2016模型输出值与COSMIC掩星观测数据的一致性在太阳活动高年比太阳活动低年更好。模型NmF2输出值与COSMIC掩星观测值的相关性高于hmF2输出值与COSMIC掩星观测值的相关性。

      2)在太阳活动低年春秋季的中午时段,与COSMIC掩星数据相比,IRI-2016模型在30°N~55°N对NmF2和hmF2分别存在低估和高估现象;在15°N~30°N区域内则恰恰相反。

      3)以COSMIC掩星观测值为参照,IRI_CCIR模型和IRI_URSI模型均对NmF2在中午时分存在高估现象,该现象在太阳活动低年更为显著。而在春秋季的LT 17:00—04:00时间段,采用这两个模型选项得到的IRI-2016模型输出值均低于COSMIC观测值。

      4)在太阳活动低年,无论任何季节,IRI_CCIR模型和IRI_AMTB模型的hmF2输出值均高于COSMIC掩星观测值,其中IRI_AMTB模型输出值相对于掩星观测的偏差更大,尤其是在冬季。IRI_Shubin模型的hmF2输出值与COSMIC观测值最为接近,即使在太阳活动高年,两者的偏差也小于其他两种模型。这表明COSMIC掩星资料在IRI_Shubin模型改进hmF2输出值中发挥了重要作用。

      5)根据本文的研究结果,在中国区域采用IRI-2016计算NmF2和hmF2时,分别推荐使用IRI_CCIR和IRI_Shubin模型,而非官方推荐的IRI_URSI和IRI_AMTB模型。

参考文献 (15)

目录

    /

    返回文章
    返回