电离层模型IRI-2016的NmF2和HmF2参数在平静期的误差分析

电离层模型IRI-2016的NmF2和HmF2参数在平静期的误差分析 Error Analysis on NmF2 and HmF2 Parameters of IRI-2016 Ionospheric Model during Magnetically Quiet Periods 1 2 3 first-author 佘丽丽，女，硕士生，现从事GNSS电离层和大气层掩星探测数据处理研究. E-mail:last_romeo@163.com 佘丽丽，女，硕士生，现从事GNSS电离层和大气层掩星探测数据处理研究. E-mail: last_romeo@163.com last_romeo@163.com 1 2 3 corresp 白伟华, E-mail: chenzz@whu.edu.cn 白伟华, E-mail: chenzz@whu.edu.cn chenzz@whu.edu.cn 1 2 3 1 2 3 1 3 1 3 1 3 1 3 1 3 1 3 1 3 中国科学院 国家空间科学中心, 北京 100190 中国科学院 国家空间科学中心, 北京 100190 National Space Science Center, The Chinese Academy of Sciences, Beijing 100190, China National Space Science Center, The Chinese Academy of Sciences, Beijing 100190, China 中国科学院大学，北京 100190 中国科学院大学，北京 100190 University of Chinese Academy of Sciences, Beijing 100190, China University of Chinese Academy of Sciences, Beijing 100190, China 北京市天基空间环境探测重点实验室，北京 100190 北京市天基空间环境探测重点实验室，北京 100190 Beijing Key Laboratory of Space Environment Explorations, Beijing 100190, China Beijing Key Laboratory of Space Environment Explorations, Beijing 100190, China 国际参考电离层模型IRI(the international reference ionosphere model)是利用全球电离层测站以及卫星观测数据建立的电离层经验模型, 并广泛用于电离层研究领域.本文选用2014年地磁暴平静期的IRI-2016模式和电离层垂测仪的数据，对电离层F2层峰值电子密度NmF2和峰值高度HmF2两个物理量，通过IRI模型最新版本(IRI-2016) 在地磁环境平静期的误差精度进行评估.根据不同磁倾角范围分别选择5个台站(Juliusruh (54°N，13°E)磁倾角：68.54°, Okinawa(26°N，128°E)磁倾角：37.26°, Jicamarca(12°S，77°W)磁倾角：-0.26°, Port_Stanley(52°S，58°W)磁倾角：-51.20°, Hermanus(34°S，19°E)磁倾角：-64.24°)进行NmF2和HmF2参数对比分析.5个台站中Port_Stanley站的NmF2相关性最差(相关系数0.68)，Okinawa站HmF2相关性最差(相关系数0.52)，Juliusruh站两个参数的相关性最高(相关系数0.87).以Juliusruh站为例进行白昼、夜晚以及不同季节的分析，可以看出白昼误差明显低于夜晚，相较于分季(春或秋)和冬季，IRI模型在夏季预测效果最好.5个台站的研究结果表明，IRI模型与垂测仪数据具备一致性，IRI模型的值高于实测数据值. The international reference ionosphere (IRI) model has been widely utilized in ionospheric research. It is an empirical ionosphere model based on worldwide available data of ground-based and satellite observations. In this paper, the NmF2 and HmF2 data of IRI-2016 model and ionosonde data during geomagnetic quiet periods in 2014 have been selected for a comparison study to estimate the accuracy and precision of the IRI-2016 model. This paper chooses five stations(Juliusruh(54°N, 13°E):68.54°, Okinawa(26°N, 128°E):37.26°, Jicamarca(12°S, 77°W):-0.26°, Port_Stanley(52°S, 58°W):-51.20°, Hermanus(34°S, 19°E):-64.24°) from different geomagnetic dip zones to compare and analyze NmF2 and HmF2 parameters. NmF2 data showed the worst correlation ( r =0.68) in Port_Stanley of the five stations and HmF2 data showed the worst correlation ( r =0.52) in Okinawa of the five stations. Juliusruh had the highest correlation ( r =0.87), both for NmF2 and HmF2. In the case of Juliusruh, analysis on day and night and seasonal characteristics indicated that the error in daytime was obviously lower than that in nighttime. Compared with winter and equinoxes, summer has the best agreement in IRI model. The results of five stations showed that IRI predictions had a good agreement with ionosonde data both for the NmF2 and HmF2, and IRI predictions systematically were higher than the observed results. 国际参考电离层模型 垂测仪 NmF2 HmF2 the international reference ionosphere model ionosonde NmF2 HmF2 1 国际参考电离层模型(the international reference ionosphere model, IRI模型)是一个由国际空间研究委员会(committee on space research, COSPAR)和国际无线电科学联盟(international union of radio science, URSI)联合发起的国际项目，其目的就是建立一个电离层经验标准模型.IRI模型是一个能够计算给定位置和给定时间电离层的电子密度、电子温度、离子温度、离子成分等参数的经验模型.IRI模型很好的描述了电离层各物理量的分布情况，故其对理论研究以及实际应用都有非常大的价值.IRI模型作为背景电离层主要用于无线电磁波传播，以及通过大地测量方法获得的电离层参数，进行可靠性、精度验证和数据分析等方面的科学研究 [ 1 1 ~ 8 8 1~8 ] ；层析成像技术和无线电掩星反演技术也都以IRI作为背景场进行数据分析和算法简化 [ 8 8 ~ 11 11 8~11 ] ；IGS(the international GNSS service)还将IRI模型作为二阶改正用于导航和定位中修正电离层延迟误差 [ 8 8 , 12 12 ~ 14 14 12~14 ] , 很多学者也都根据IRI模型推导出了GPS双频的高阶改正量 [ 8 8 , 15 15 ] ；此外，IRI模型还可用于对空间天气的预报 [ 8 8 , 16 16 ] . 目前已经有一些学者对IRI-2012及以下版本模型与垂测仪数据的对比进行了研究.相关文献指出IRI-2012及以下版本模型预测精度与季节相关，其中冬季的误差较大，夏季的误差较小 [ 17 17 , 18 18 ] ，同时预测还受纬度的影响 [ 19 19 ] ；文献[ 20 20 ]通过对比垂测仪的数据，发现IRI-2012及以下版本模型对白昼的预测误差要小于晚上. IRI模型主要是针对平静期电离层特性进行描述 [ 21 21 ] , 所以本文选择IRI-2016版本(IRI模式最新版本)对2014年地磁暴平静期的模式数据和电离层垂测仪的数据进行分析研究.针对电离层F2层峰值电子密度NmF2和峰值高度HmF2两个物理量，进行时间和空间对应的IRI模式与垂测仪数据对比，进而完成IRI模型最新版本(IRI-2016) 在地磁环境平静期NmF2和HmF2参数的误差精度评估. IRI模型介绍 1 IRI模型是利用全球遍布的电离层垂测站以及卫星观测数据建立的电离层经验模型.IRI模型主要的数据来源于全球电离层探测仪网络、非相干散射雷达、火箭探空数据、地面接收的卫星数据和电离层探测卫星数据等 [ 1 1 , 8 8 , 13 13 ] .IRI软件包内包含FORTRAN子程序、模型系数、IRI-2016模型的索引文件、README文件和许可文件.用户可通过调用IRI源代码进行相关参数的计算，如输入地点、时间、高度范围、模型选择和部分参数等，IRI即可调用相关子程序进行计算得到参数计算结果 [ 22 22 ] . 图 1 图 1 是根据IRI模型画出的2014年7月1日的全球电子密度峰值NmF2分布图.IRI模型是根据大量实测数据不断修正的一种经验模型，并且根据最新的观测数据进行不断的修正.目前已经发布了几个版本，并在不断改进.IRI官方网站可以提供IRI模型(2016、2012、2007、2001版本)的Fortran源代码、太阳和地磁的指数文件、实时的临界频率( f 0 F2) 和HmF2数据、IRI-2012的Matlab版本以及IRI-2001的PC windows版本，网站还提供在线计算和画图功能(基于IRI 2012、2007和2001版本).虽然目前的IRI版本已经更新到2016版本，但是大量的数据对比都是基于IRI-2012及之前版本，而最新版本IRI-2016目前还没有相应的比较，本文基于最新的IRI-2016版本进行了数据比对分析. (a) 00:00; (b) 06:00; (c) 12:00; (d) 18:00 (a) 00:00; (b) 06:00; (c) 12:00; (d) 18:00 数据选择和处理方法 1 数据时间选择 2 地磁暴是地球磁场全球性的剧烈扰动现象，以地磁指数来表征. K p 指数(the planetary three-hour-range K p index)是由全球地磁台网中13个地磁台站的K指数计算得到的，用于表示全球地磁活动性，每3 h一个值，其值逐渐增大表示地磁扰动逐渐增强.一般把 K p < 4称为平静期， K p =5，6称为中小磁暴， K p =7，8，9称为大磁暴 [ 23 23 , 24 24 ] .由于2015年SPIDR(space physics interactive data resource)网站提供的垂测仪数据不完整，只能获取前半年数据，且地磁暴平静期没有2014年多，故本文选择2014年的电离层垂测仪的数据进行分析研究. 图 2 图 2 是2014年的 K p 指数图，红线表示 K p =4，由 图 2 图 2 可知，2014年平静期占到全年的98%，本文过滤掉了不属于平静期的2%的数据，只选取2014年平静期数据进行分析. 电离层垂测仪数据选择 2 从SPIDR网站下载各台站垂测仪数据，由于设备、天气等客观原因，可能会导致数据的错误和缺失.在全球61个台站的 f 0 F2和HmF2两个物理量的数据可用率高，磁倾角为-80°~-40°，-40°~0°，0°附近，0°~40°，40°~80°的范围内分别选择5个台站的垂测仪实测数据进行分析.5个台站分别是Juliusruh站(位置(54°N，13°E)，磁倾角68.54°)、Okinawa站(位置(26°N，128°E)，磁倾角37.26°)、Jicamarca站(位置(12°S，77°W)，磁倾角-0.26°)、Port_Stanley站(位置(52°S，58°W)，磁倾角-51.20°)和Hermanus站(位置(34°S，19°E)，磁倾角-64.24°)，由于在-40°~0°范围内没有数据可用率达到50%以上的台站，所以选择Port_Stanley(磁倾角-51.20°)代替.全球61个台站在地图上的位置如 图 3 图 3 所示.蓝色圆圈是其他56个台站，红色五角星是本文选定的5个台站. 表 1 表 1 为5个台站垂测仪的经纬度、数据可用率、磁倾角等信息，可以看出5个台站的数据可用率都达到了60%以上，保证了数据分析的可靠性. 2 台站 2 经纬度 2 磁倾角/(°) 2 数据可用率/% f 0 F2 HmF2 Juliusruh 54°N，13°E 68.54 86 86 Okinawa 26°N，128°E 37.26 73 73 Jicamarca 12°S，77°W -0.26 82 82 Port_Stanley 52°S，58°W -51.20 86 86 Hermanus 34°S，19°E -64.24 84 84 数据处理 2 电离层垂测仪下载的原始数据包括F2层临界频率 f 0 F2和峰值高度HmF2, 其中 f 0 F2需要转换为峰值电子密度NmF2，计算公式如下： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459111&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459111&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459111&type=middle (1) 式中, f 0 F2表示F2层的临界频率(MHz)，NmF2表示峰值电子密度(cm -3 ).本文对IRI数据与电离层垂测仪数据进行对比分析，主要计算了误差的期望和标准差.其中，NmF2和HmF2的参数误差分别采用相对误差和绝对误差进行评估. NmF2的相对误差计算公式为： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459116&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459116&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459116&type=middle HmF2的绝对误差公式为： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459121&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459121&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459121&type=middle (2) 和(3) 式中 x 分别为NmF2和HmF2的IRI模式数据， x 0 分别为NmF2和HmF2的观测数据. 误差期望的计算公式为： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459126&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459126&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459126&type=middle (4) 式中 n 表示数据个数， e 表示误差，可以是相对误差 γ 也可以是绝对误差 δ ，下角标 i 代表对应的数据标号， e i 代表相应误差，误差标准差的计算公式如下: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459130&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459130&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=4459130&type=middle 误差统计结果 1 分别对5个台站画出IRI模型与实测HmF2和NmF2数据的日均变化散点图，通过计算相关系数，对NmF2和HmF2分别进行相对误差统计分析和绝对误差统计分析.然后以Juliusruh站为例进行季节性的误差分析和昼夜的误差分析，对IRI模型完成一致性的有效评估. F2层峰值电子密度NmF2分析 2 5个台站2014年NmF2分析 3 图 4 图 4 是2014年5个台站全年平均每天的NmF2统计分析图. 图 4(a) 图 4(a) ~ (e) (e) 是以IRI模型数据为 x 轴，垂测仪数据为 y 轴的相关性散点图，红线为对数据点做的线性拟合，绿线为 y = x 线，表示如果IRI预测值和观测值相同，那么该数据点应分布在 y = x 线上， r 为散点图的相关系数.从 图 4 图 4 中可以看出Port_Stanley站相关性最差, 相关系数为0.68，Okinawa站次之，其他几个站相关性都比较好，相关系数均大于0.8. (a)(f)(k):Juliusruh站; (b)(g)(l):Okinawa站; (c)(h)(m):Jicamara站; (d)(i)(n):Port_Stanley站; (e)(j)(o):Hermanus站; (a)~(e):数据相关性散点图; (f)~(j):IRI模型和垂测仪数据的日均变化图; (k)~(o):日均变化的相对误差的期望和标准差图 (a)(f)(k):Juliusruh station; (b)(g)(l):Okinawa station; (c)(h)(m):Jicamara station; (d)(i)(n): Port_Stanley station; (e)(j)(o): Hermanus station; (a)~(e):data correlation scatter diagram; (f)~(j):daily change between IRI model and ionosonde data; (k)~(o):expectation and standard deviation of daily change relative error 图 4(f) 图 4(f) ~ (j) (j) 是2014年5个台站IRI模型和垂测仪数据的日均变化图，虚线是垂测仪实测数据，实线是IRI模型数据.Jicamarca站在当地时间08:00和16:00出现了双峰结构，在05:00时有极小值，日出和夜晚偏差较大，其他时间一致性良好.而其他几个站大部分时间段吻合度都比较高，在中午出现一个峰值，中午偏差比较大，其中Okinawa站中午偏差达到最大. 图 4(k) 图 4(k) ~ (o) (o) 是2014年5个台站日均变化的相对误差的期望和标准差图，其中虚线是标准差，实线是期望.Jicamara站和Okinawa站在当地时间00:00到08:00标准差增大，最大达到109%，且出现震荡情况，之后迅速下降，保持平缓；Juliusruh站标准差变化比较平缓，在20%~50%之间；而Hermanus站和Port_Stanley站的标准差则在00:00~17:00间曲线较为平缓，在18:00时出现峰值，达到148%.期望与标准差趋势基本保持一致，只有Okinawa站与其标准差趋势不同，全天都比较平缓，期望维持在10%-40%之间.综上所述，5个台站中最大标准差出现在Hermanus站，18:00时达到148%；最大期望出现在Port_Stanley站，达到112%. 表 2 表 2 是2014年5个台站全天的NmF2的相对误差期望和标准差结果，可以看出误差的期望最接近0的是Juliusruh站和Hermanus站，标准差最小的则是Juliusruh站和Jicamarca站，所以IRI模型对Juliusruh站的NmF2的预测效果最好. 台站 经纬度 磁倾角/(°) 相对误差的期望/% 相对误差的标准差/% Juliusruh 54°N，13°E 68.54 2.24 34.69 Okinawa 26°N，128°E 37.26 23.98 63.33 Jicamarca 12°S，77°W -0.26 29.94 53.01 Port_Stanley 52°S，58°W -51.20 49.41 83.54 Hermanus 34°S，19°E -64.24 1.65 59.97 均值 21.44 58.91 以Juliusruh站为例的NmF2分析 3 由于Juliusruh站数据可用率高，相关性好，所以以此台站NmF2数据为例进行昼夜以及四季的对比分析，白昼的当地时间范围是08:00~18:00，夜晚的当地时间范围是18:00~06:00. 表 3 表 3 是Juliusruh站NmF2全天、白昼、夜晚相对误差的期望和标准差的统计表格.由 表 3 表 3 可知，白昼的相对误差的期望和标准差都低于夜晚. 台站 经纬度 磁倾角/(°) 时间/h 相对误差的期望/% 相对误差的标准差/% 3 Juliusruh 3 54°N，13°E 3 68.54 全天(00:00-24:00) 2.24 34.69 白昼(08:00-18:00) -1.57 29.51 夜晚(18:00-06:00) 6.92 37.85 图 5 图 5 是2014年Juliusruh站夏季、分季、冬季的NmF2统计分析图. 图 5(a) 图 5(a) ~ (c) (c) 是相关性散点图，可以看出冬季相关性最好，分季次之，夏季相关性最差，但数据点比较集中，3个季节相关性都达到0.7以上. 图 5(d) 图 5(d) ~ (f) (f) 是2014年Juliusruh站IRI模型和垂测仪数据在3个季节的日均变化图，可以看出夏季两条曲线趋势非常相近，误差最小，而分季夜晚曲线也很相近，但白昼偏差较大.冬季则在当地时间07:00到12:00偏差略大，其他时间段曲线非常相似，冬季出现了很陡的峰，其值很大，这也就是“冬季异常”现象 [ 25 25 ] . 图 5(g) 图 5(g) ~ (i) (i) 是2014年Juliusruh站3个季节的日均变化的相对误差的期望和标准差图.3个季节的标准差都比较平缓，变化不大，但是非常明显的是夏季标准差最小，在20%~30%之间，标准差最大出现在冬季，最大值达到66%；相对误差的期望在冬季变剧烈，在08:00时出现最小值-38%，在18:00时出现峰值69%，而分季和夏季期望变化不明显.综上可以看出IRI模型与垂测仪实测数据Juliusruh站夏季的相关性虽然较低，但是夏季的误差最小，所以IRI模型在夏季预测是十分准确的. (a)(d)(g):夏季; (b)(e)(h):分季; (c)(f)(i):冬季; (a)~(c):数据相关性散点图; (d)~(f):IRI模型和垂测仪数据的日均变化图; (g)~(i):日均变化的相对误差的期望和标准差图 (a)(d)(g):summer; (b)(e)(h):equinoxes; (c)(f)(i):winter; (a)~(c):data correlation scatter diagram; (d)~(f):daily change between IRI model and ionosonde data; (g)~(i):expectation and standard deviation of daily change relative error 表 4 表 4 是Juliusruh站NmF2夏季、分季、冬季相对误差的标准差、期望的统计表格，由表可知，夏季的标准差和期望都是最小的，这也验证了上一段的结论. 台站 经纬度 磁倾角/(°) 季节 相对误差的期望/% 相对误差的标准差/% 3 Juliusruh 3 54°N，13°E 3 68.54 夏季 2.09 23.01 分季 -5.12 37.09 冬季 17.70 46.80 F2层峰值高度HmF2分析 2 5个台站2014年HmF2分析 3 图 6 图 6 是2014年5个台站全年平均每天(即每一天的0点到24点)的HmF2统计分析图. 图 6(a) 图 6(a) ~ (e) (e) 是2014年5个台站的HmF2相关性散点图，从 图 6 图 6 中可以看出Okinawa站相关性最差，此站的NmF2相关系数也很低，其他几个站相关性都比较好，都达到0.7以上，其中Juliusruh站相关性最高. (a)(f)(k):Juliusruh站; (b)(g)(l):Okinawa站; (c)(h)(m):Jicamara站; (d)(i)(n):Port_Stanley站; (e)(j)(o):Hermanus站; (a)~(e):数据相关性散点图; (f)~(j):IRI模型和垂测仪数据的日均变化图; (k)~(o):日均变化的绝对误差的期望和标准差图 (a)(f)(k):Juliusruh station; (b)(g)(l):Okinawa station; (c)(h)(m):Jicamara station; (d)(i)(n): Port_Stanley station; (e)(j)(o): Hermanus station; (a)~(e):data correlation scatter diagram; (f)~(j):daily change between IRI model and ionosonde data; (k)~(o):expectation and standard deviation of daily change relative absolute error 图 6(f) 图 6(f) ~ (j) (j) 是2014年5个台站IRI模型和垂测仪数据的日均变化图.Jicamarca站IRI模型的预测出现峰值提前的现象，实测数据峰值出现在当地时间19:00，而IRI模型峰值出现在12:00，白昼的偏差很大，一致性不是很好；Hermanus站和Okinawa站都出现了两个波谷，不同的是，Hermanus站IRI预测值略高于实测数据的值，Okinawa站IRI预测值低于实测数据的值；Juliusruh站和Port_Stanley站HmF2从00:00开始不断下降，到达当地时间08:00趋于平缓，17:00逐渐上升，Juliusruh站IRI模型与实测数据十分吻合，Port_Stanley站IRI预测值略高于实测值. 图 6(k) 图 6(k) ~ (o) (o) 是2014年5个台站日均变化的绝对误差的期望和标准差图，除Jicamarca站外其他4个站期望曲线平缓，取值范围小，较集中，表明IRI模型在这几个站的模型数据是十分有效的.然而Jicamarca站的变化起伏很大，第一个极小值出现在当地时间07:00，再迅速上升，在11:00有最大值71 km，之后又快速下降，在20:00有最小值-49 km.Juliusruh站的标准差最接近0，这个台站的预测在HmF2也是很有效的.Okinawa站的标准差曲线略有震荡，其他几个站标准差变化也不大，标准差最大出现在02:00的Okinawa站，达到73 km. 表 5 表 5 是2014年5个台站全天的HmF2统计表格，统计了5个台站的绝对误差期望和标准差的具体数值，可以看出期望最接近0的是Juliusruh站，标准差最小的则是Juliusruh站和Hermanus站，所以IRI模型对NmF2的预测效果最好的也是Juliusruh站. 台站 经纬度 磁倾角/(°) 绝对误差的期望/km 绝对误差的标准差/km Juliusruh 54°N，13°E 68.54 2.08 20.28 Okinawa 26°N，128°E 37.26 -15.58 45.72 Jicamarca 12°S，77°W -0.26 14.34 45.26 Port_Stanley 52°S，58°W -51.20 15.72 32.04 Hermanus 34°S，19°E -64.24 15.30 24.21 均值 6.37 33.50 以Juliusruh站为例的HmF2分析 3 依然选择Juliusruh站进行HmF2数据的昼夜以及四季的对比分析. 表 6 表 6 是Juliusruh站NmF2全天、白昼、夜晚绝对误差的标准差、期望的统计表格.由 表 6 表 6 可知，白昼的绝对误差的标准差低于夜晚，但白昼的期望高于夜晚. 台站 经纬度 磁倾角/(°) 时间 绝对误差的期望/km 绝对误差的标准差/km 3 Juliusruh 3 54°N，13°E 3 68.54 全天 2.08 20.28 白昼 3.29 15.46 夜晚 -0.02 24.34 图 7 图 7 是2014年Juliusruh站夏季、分季、冬季的HmF2统计分析图. 图 7(a) 图 7(a) ~ (c) (c) 是相关性散点图，可以看出3个季节的相关系数一样，相关性都比较高. 图 7(d) 图 7(d) ~ (f) (f) 是2014年Juliusruh站IRI模型和垂测仪数据在3个季节的日均变化图，可以看出3个季节IRI模型和实测数据曲线趋势非常相近，一致性都很好.其中，夏季在中午偏差较大，而分季和冬季在日出时分偏差较大. 图 7(g) 图 7(g) ~ (i) (i) 是2014年Juliusruh站3个季节的日均变化的绝对误差的期望和标准差图，3个季节的绝对误差的期望震荡趋势非常相似，最大值出现在冬季达到28 km；3个季节的标准差趋势也很相近，本地时00:00到20:00都很平缓，之后略有上升，在24:00有最大值33 km. (a)(d)(g):夏季; (b)(e)(h):分季; (c)(f)(i):冬季; (a)~(c):数据相关性散点图; (d)~(f):IRI模型和垂测仪数据的日均变化图; (g)~(i):日均变化的绝对误差的期望和标准差图 (a)(d)(g):summer; (b)(e)(h):equinoxes; (c)(f)(i):winter; (a)~(c):data correlation scatter diagram; (d)~(f):daily change between IRI model and ionosonde data; (g)~(i):expectation and standard deviation of daily change absolute error 表 7 表 7 是Juliusruh站在夏季、分季、冬季HmF2参数绝对误差的标准差、期望的统计表格.由 表 7 表 7 可知，夏季的标准差最小，分季的期望最小.从 图 7 图 7 可以看出，相对而言，HmF2在夏季曲线最为平缓，所以HmF2在夏季的预测也是最为准确的. 台站 经纬度 磁倾角/(°) 季节 绝对误差的期望/km 绝对误差的标准差/km 3 Juliusruh 3 54°N，13°E 3 68.54 夏季 4.04 19.31 分季 -1.56 21.64 冬季 4.27 22.37 结论 1 本文选择2014年全年地磁暴平静期的IRI-2016模式数据和电离层垂测仪的数据，完成了NmF2和HmF2两个物理参数的误差精度评估.通过对2014年全年5个台站的NmF2和HmF2统计分析发现，IRI模型与实测数据的相关性较高，Juliusruh站一致性最好，Okinawa站和Port_Stanley一致性最差；两个参数的日均变化趋势一致，IRI模型值高于实测数据的值且白昼偏差较大，NmF2最大偏差出现在Okinawa站的中午，HmF2的最大偏差出现在Jicamarca站的中午；统计NmF2的相对误差的最大标准差是Hermanus站，最大期望出现在Port_Stanley站；HmF2的绝对误差的最大期望出现在Jicamarca站，最大标准差出现在Okinawa站.选择与IRI模型吻合度最高的Juliusruh站为例进行单点分析，对于NmF2而言，白昼的相对误差的期望和标准差都低于夜晚，在季节方面夏季误差最小，最大标准差出现在冬季(可达到66%)；对于HmF2而言，白昼的绝对误差的标准差低于夜晚，但白昼的期望高于夜晚，在季节方面没有N mF2明显，3个季节相关性差异不大，IRI模型与实测数据的一致性也都很高，但夏季的绝对误差小于其他两季节，所以IRI模型对夏季的预测较其他两个季节更为准确一些. 总之，对于NmF2和HmF2两个参数而言，在平静期IRI模型的数据是比较准确的，5个台站平均NmF2相对误差期望和平均HmF2的绝对误差都很小.以上分析初步获取了IRI-2016模式误差的统计特性，为模式的应用提供了先验误差信息，而且为模式的进一步改进奠定了基础.由于垂测仪数据数量和质量的限制，本研究在统计意义上尚不够完善，下一步计划获取更多的垂测仪数据(如子午工程的台站数据)，并且对其数据进行质量评估和控制，进而获取更具统计意义的IRI模式误差特性. 参考文献 BILITZA D, REINISCH B W. International reference ionosphere 2007: Improvements and new parameters[J]. Advances in Space Research, 2008, 42 (4):599-609. DOI:10.1016/j.asr.2007.07.048. MCKINNELL L A, FRIEDRICH M, STEINER R J. A new approach to modelling the daytime lower ionosphere at auroral latitudes[J]. Advances in Space Research, 2004, 34 (9):1943-1948. DOI: 10.1016/j.asr.2004.05.005. MCKINNELL L A, FRIEDRICH M. Results from a new auroral lower ionosphere model[J]. Advances in Space Research, 2006, 37 (5):1045-1050. DOI: 10.1016/j.asr.2004.11.008. HERNANDEZ-PAJARES M, JUAN J, SANZ J, et al. Combining GPS measurements and IRI model values for space weather specification[J]. Advances in Space Research, 2002, 29 (6):949-958.DOI: 10.1016/S0273-1177(02)00051-0. NIRANJAN K, SRIVANI B, GOPIKRISHNA S, et al. Spatial distribution of ionization in the equatorial and low-latitude ionosphere of the Indian sector and its effect on the Pierce point altitude for GPS applications during low solar activity periods[J].Journal of Geophysical Research, 2007, 112 (A5):A05304, DOI:10.1029/2006JA011989. HOCKE K, IGARASHI K. Structure of the earth's lower ionosphere observed by GPS/MET radio occultation[J].Journal of Geophysical Research, 2002, 107 (A5), SIA4-1. DOI:10.1029/2001JA900158. DEAR R, MITCHELL C. GPS interfrequency biases and total electron content errors in ionospheric imaging over Europe[J]. Radio Science, 2006, 41 : RS6007, DOI:10.1029/2005RS00 3269. BILITZA D, MCKINNELLL A, REINISCH B W, et al. The international reference ionosphere today and in the future[J]. Journal of Geodesy, 2001, 85 (12): 909-920, DOI:10.1007/s00 190-010-0427-x. BUST G S, GARNER T W, GAUSSIRAN II T L. Ionospheric data assimilation three-dimensional (IDA3D): A global, multisensory, electron density specification algorithm[J]. Journal of Geophysical Research, 2004, 109 : A11312. DOI:10.1029/2003JA010234. LEE J K, KAMALABADI F, MAKELA J J. Three-dimensional tomography of ionospheric variability using a dense GPS receiver array[J]. Radio Science, 2008, 43 :RS3001, DOI:10.1029/ 2007RS003716. GARCIA R, CRESPON F. Radio tomography of the Ionosphere:Analysis of an underdetermined, ill-posed inverse problem, and regional application[J]. Radio Science, 2008, 43 :RS2014, DOI:10.1029/2007RS00314. BILITZA D, KOBLINSKY C, BECKLEY B, et al. Using IRI for the computation of ionospheric corrections for altimeter data analysis[J]. Advances in Space Research, 1995, 15 (2):113-120.DOI:10.1016/S0273-1177(99)80032-5. BILITZA D. Internationalreference ionosphere 2000[J]. Radio Science, 2001, 36 (2):261-275. DOI: 10.1029/2000RS002432. ORùS R, HERNÁNDEZ-PAJARES M, JUAN J M, et al. Performance of different TEC models to provide GPS ionospheric corrections[J]. Journal of Atmospheric and Solar-terrestrial Physics, 2002, 64 (18): 2055-2062.DOI:10.1016/S1364-6826(02)00224-9. DATTA-BARUA S, WALTER T, BLANCH J, et al. Bounding higher-order ionospheric errors for the dual-frequency GPS user[J]. Radio Science, 2008, 43 :RS5010.DOI:10.1029/2007RS003772. GARNER T W, GAUSSIRAN T L II, TOLMAN B W, et al. Total electron content measurements in ionospheric physics[J]. Advances in Space Research, 2008, 42(4):720-726, DOI:10.1016/j.asr.2008.02.025. BILITZA D, STEVEN A B, MATHEW Y W, et al. Measurements and IRI model predictions during the recent solar minimum[J]. Journal of Atmospheric and Solar-terrestrial Physics, 2012, 86 :99-106.DOI: 10.1016/j.jastp.2012.06.010. RIOS V H, MEDINA C F, ALVAREZ P. Comparison between IRI predictions and digisonde measurements at Tucuman[J].Journal of Atmospheric and Solar-terrestrial Physics, 2007, 69 (4-5): 569-577.DOI: 10.1016/j.jastp.2006.10.004. 黄为权, 肖赛冠, 肖佐. IRI-2007对扩展F的预测与观测比较研究[J].空间科学学报, 2009, 29 (3):2 75-280.DOI: 10.11728/cjss2009.03.275. HUANG W Q, XIAO S G, XIAO Z. Comparison study of IRI-2007 spread-F occurrence predictions and observations[J]. Chinese Journal of Space Science, 2009, 29 (3):275-280.1111111155555555 MURTAZA G, IQBAL S, AMEEN M A. Comparing IRI and regional ionospheric model with ionosonde measurements in pakistan[J]. Advances in Space Research, 2008, 42 (4):682-690.DOI: 10.1109/IconSpace.2015.7283755. 万德焕, 黄江, 邓柏昌, 等.基于IRI-2012模型广州地区f 0 F 2 实测与预测的对比分析[J].空间科学学报, 2015, 35 (2): 166-177.DOI: 10.11728/cjss2015.02.166. WAN D H, HUANG J, DENG B C, et al. Comparison of f 0 F 2 between observation and the prediction based on IRI-2012 over Guangzhou[J]. Chinese Journal of Space Science, 2015, 35 (2): 166-177. DOI: 10.11728/cjss2015.02.166. 王成, 王解先.国际参考电离层在电离层电子密度特征分析中的应用[J].大地测量与地球动力学, 2012, 32 (2): 55-59. WANG C, WANG J X. International reference ionosphere and its application features analysis of ionosphere electron density [J]. Journal of Geodesy and Geodynamics, 2012, 32 (2): 55-59(Ch). BARTELS J.The Standardized Index Ks and the Planetary Index Kp[M]. Paris: IUGG Publication Office, 1949. BARTELS J, HECK N H, JOHNSTON H F. The three-hour-range index measuring geomagnetic activity [J]. Journal of Geophysical Research, 1939, 44 (4):411-454.DOI: 10.1029/TE044i004p00411. RISHBETH H, MULLER-WODARG I C F, ZOU L. Annual and semiannual variations in the ionospheric F2-layer: Ⅱ. Physical discussion[C]// Ann Geophys, Berlin: Springer, 2000, 18 :945-956. DOI: 10.1007/s00585-000-0945-6 http://dx.doi.org/10.1007/s00585-000-0945-6 . 10.14188/j.1671-8836.2017.03.001 P352 whdxxblxb-63-3-189 国家自然科学基金青年科学基金项目(41505030, 41405039, 41405040和41606206) 资助 2016-10-20 2017-06-24 2020-12-09 版权所有© 2017 《武汉大学学报(理学版)》编辑部 Copyright © 2017 Journal of Wuhan University(Natural Science Edition). All rights reserved. 2017 武汉大学学报(理学版) 3 63