留言板

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

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

NTCM广播电离层模型性能评估

赵齐乐 杨超 耿涛 马福建 张强

赵齐乐, 杨超, 耿涛, 马福建, 张强. NTCM广播电离层模型性能评估[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
引用本文: 赵齐乐, 杨超, 耿涛, 马福建, 张强. NTCM广播电离层模型性能评估[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
ZHAO Qile, YANG Chao, GENG Tao, MA Fujian, ZHANG Qiang. Performance Analysis of NTCM Ionospheric Broadcast Models[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
Citation: ZHAO Qile, YANG Chao, GENG Tao, MA Fujian, ZHANG Qiang. Performance Analysis of NTCM Ionospheric Broadcast Models[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348

NTCM广播电离层模型性能评估

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

国家自然科学基金 41674004

详细信息
    作者简介:

    赵齐乐, 博士, 教授, 主要从事GNSS精密定轨定位数据处理方法研究。zhaoql@whu.edu.cn

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

Performance Analysis of NTCM Ionospheric Broadcast Models

Funds: 

The National Natural Science Foundation of China 41674004

More Information
    Author Bio:

    ZHAO Qile, PhD, professor, specializes in GNSS precision orbit and positioning. E-mail:zhaoql@whu.edu.cn

    Corresponding author: GENG Tao, PhD, associate professor. E⁃mail:gt_gengtao@whu.edu.cn
  • 摘要: 对于全球卫星导航系统单频用户而言,广播电离层模型是改正电离层误差的主要方法之一。对NTCM-BC模型、MNTCM-BC模型、NTCM-Klobuchar模型和传统Klobuchar模型的性能进行了对比分析。结果表明,在一个太阳活动周期内,NTCM-BC和MNTCM-BC模型计算的垂直总电子含量精度相当,且精度最高,NTCM-Klobuchar模型次之,Klobuchar模型最低。在低纬度测站,NTCM系列模型较Klobuchar模型最高能减少20 TECU左右的误差。与Klobuchar模型相比,采用NTCM-BC和MNTCM-BC模型进行伪距单点定位,高程精度提高约38%,特别是太阳活动高峰期间,低纬度测站最高能提高3 m左右。
  • 图  1  2005—2016年电离层模型的RMS值时序图

    Figure  1.  RMS Values of Ionospheric Models forthe Year of 2005—2016

    图  2  2005—2016年每天电离层模型RMS值随F10.7指数的变化

    Figure  2.  Daily RMS Values of the Ionospheric Models for the Year of 2005—2016 with F10.7 Index

    图  3  北纬4个IGS测站的电离层模型每小时RMS值时序图

    Figure  3.  Hourly RMS Values of VTEC of the Ionospheric Models at 4 IGS Stations in Northern Latitude

    图  4  南纬3个IGS测站的电离层模型每小时RMS值时序图

    Figure  4.  Hourly RMS Values of VTEC of the Ionospheric Models at 3 IGS Stations in Southern Latitude

    表  1  IGS测站信息

    Table  1.   IGS Station Coordinates

    测站名 地理纬度 地理经度 城市 国家
    YAKT 62.030 8°N 129.680 2°E 雅库茨克 俄罗斯
    BJFS 39.608 6°N 115.892 2°E 北京 中国
    TCMS 24.797 7°N 120.987 2°E 新竹 中国
    PIMO 14.635 5°N 121.077 5°E 奎松市 菲利宾
    DARW 12.843 6°S 131.132 5°E 达尔文 澳大利亚
    PERT 31.801 7°S 115.885 0°E 珀斯 澳大利亚
    MAC1 54.499 4°S 158.935 6°E 麦夸里岛 澳大利亚
    下载: 导出CSV

    表  2  2009年12月9日的定位精度/m

    Table  2.   RMS of Position Error on December 9th, 2009/m

    测站名 Klobuchar NTCM-Klobuchar NTCM-BC MNTCM-BC
    水平 高程 水平 高程 水平 高程 水平 高程
    YAKT 1.42 3.17 1.19 1.83 1.27 1.78 1.28 1.77
    BJFS 1.09 3.71 1.23 2.82 1.04 2.04 1.04 2.04
    TCMS 1.27 3.35 1.09 3.45 1.17 2.01 1.23 2.03
    PIMO 1.04 3.17 1.01 3.57 1.05 2.36 1.06 2.39
    DARW 0.99 2.93 0.92 2.18 0.87 1.47 0.88 1.50
    PERT 1.31 3.15 1.39 2.21 1.43 2.01 1.48 2.04
    MAC1 1.12 2.07 1.04 1.85 1.04 1.66 1.06 1.67
    平均值 1.18 3.08 1.12 2.56 1.12 1.90 1.15 1.92
    下载: 导出CSV

    表  3  2014年9月12日的定位精度/m

    Table  3.   RMS of Position Error on September 12th, 2014/m

    测站名 Klobuchar NTCM-Klobuchar NTCM-BC MNTCM-BC
    水平 高程 水平 高程 水平 高程 水平 高程
    YAKT 1.28 3.08 1.22 1.48 1.26 1.56 1.15 1.65
    BJFS 1.92 2.20 1.44 1.85 1.47 2.13 1.38 2.50
    TCMS 6.18 4.77 5.40 3.21 5.52 3.29 5.59 3.24
    PIMO 2.45 5.76 2.66 3.07 2.43 3.11 2.46 3.14
    DARW 5.00 6.90 3.41 3.62 3.56 3.64 3.47 3.65
    PERT 1.12 3.04 1.04 2.09 1.07 2.04 1.00 2.10
    MAC1 1.27 2.10 1.29 1.80 1.30 1.99 1.30 1.94
    平均值 2.75 3.98 2.35 2.45 2.37 2.54 2.34 2.60
    下载: 导出CSV
  • [1] 章红平, 韩文慧, 黄玲, 等.地基GNSS全球电离层延迟建模[J].武汉大学学报·信息科学版, 2012, 37(10):1186-1189 http://ch.whu.edu.cn/CN/Y2012/V37/I10/1186

    Zhang Hongping, Han Wenhui, Huang Ling, et al. Modeling Global Ionospheric Delay with IGS Ground-Based GNSS Observations[J]. Geomatics and Information Science of Wuhan University, 2012, 37(10):1186-1189 http://ch.whu.edu.cn/CN/Y2012/V37/I10/1186
    [2] Klobuchar J A. Ionospheric Time Delay Algorithm for Single Frequency GPS Users[J]. IEEE Transactions on Erospace and Electronic Systems, 1987, 23(3): 325-331 doi:  10.1109-TAES.1987.310829/
    [3] Hochegger G, Nava B, Radicella S M, et al. A Family of Ionospheric Models for Different Uses[J]. Phys Chem Earth Part C, 2000, 25 (4): 307-310 http://www.sciencedirect.com/science/article/pii/S1464191700000222
    [4] Wu X, Hu X, Wang G, et al. Evaluation of COMPASS Ionospheric Model in GNSS Positioning[J]. Adv Space Res, 2013, 51 (6): 959-968 doi:  10.1016/j.asr.2012.09.039
    [5] Jakowski N, Hoque M M, Mayer C. A New Global TEC Model for Estimating Transionospheric Radio Wave Propagation Errors[J]. Journal of Geodesy, 2011, 85 (12):965-974 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ea2b79b6d9f5d5b30ebc6699a0f4d07e
    [6] Hoque M M, Jakowski N, Jens B. An Ionosphere Broadcast Model for Next Generation GNSS[C]. ION GNSS, Florida, USA, 2015
    [7] Hoque M M, Jakowski N. An Alternative Ionospheric Correction Model for Global Navigation Satellite Systems[J]. Journal of Geodesy, 2015, 89(4): 391-406 doi:  10.1007/s00190-014-0783-z
    [8] Zhang X H, Ma F J, Ren X D, et al. Evaluation of NTCM-BC and a Proposed Modification for Single-Frequency Positioning[J]. GPS Solutions, 2017, 21(4):1535-1548 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=7ef1c09d16a29a6f7d2b534b5bf275dd
    [9] Hoque M M, Jakowski N, Berdermann J. Ionospheric Correction Using NTCM Driven by GPS Klobuchar Coefficients for GNSS Applications[J]. GPS Solutions, 2017, 21(4):1563-1572 doi:  10.1007/s10291-017-0632-7
  • [1] 赵齐乐, 许小龙, 马宏阳, 刘经南.  GNSS实时精密轨道快速计算方法及服务 . 武汉大学学报 ● 信息科学版, 2018, 43(12): 2157-2166. doi: 10.13203/j.whugis20180374
    [2] 张强, 赵齐乐.  OpenMP并行计算在全球电离层反演中的应用 . 武汉大学学报 ● 信息科学版, 2018, 43(2): 227-233, 240. doi: 10.13203/j.whugis20160065
    [3] 刘帅, 贾小林, 孙大伟.  GNSS星载原子钟性能评估 . 武汉大学学报 ● 信息科学版, 2017, 42(2): 277-284. doi: 10.13203/j.whugis20150344
    [4] 陈俊平, 周建华, 严宇, 陈倩, 王彬.  GNSS数据处理时空参数的相关性 . 武汉大学学报 ● 信息科学版, 2017, 42(11): 1649-1657. doi: 10.13203/j.whugis20170278
    [5] 徐克科, 伍吉仓, 雷伟伟.  基于GNSS观测网络的断层滑移时空反演 . 武汉大学学报 ● 信息科学版, 2016, 41(10): 1391-1397. doi: 10.13203/j.whugis20140400
    [6] 杨元喜, 曾安敏, 景一帆.  函数模型和随机模型双约束的GNSS数据融合及其性质 . 武汉大学学报 ● 信息科学版, 2014, 39(2): 127-131. doi: 10.13203/j.whugis20130378
    [7] 魏二虎, 万丽华, 金双根, 刘经南.  联合GNSS和SLR观测对地球自转参数的解算与分析 . 武汉大学学报 ● 信息科学版, 2014, 39(5): 581-585. doi: 10.13203/j.whugis20120213
    [8] 陈俊平, 张益泽, 谢益炳, 周旭华.  超大观测网络及多GNSS系统的快速数据处理 . 武汉大学学报 ● 信息科学版, 2014, 39(3): 253-257. doi: 10.13203/j.whugis20120712
    [9] 吴寒, 姚宜斌, 陈鹏, 肖勇.  GNSS电离层层析成像算法研究 . 武汉大学学报 ● 信息科学版, 2013, 38(12): 1405-1408.
    [10] 李英冰, 闫景仙, 熊程波, 陈中新.  GNSS观测值的压缩方法研究 . 武汉大学学报 ● 信息科学版, 2012, 37(7): 818-822.
    [11] 刘志平, 赵自强, 郭广礼.  电离层总电子含量时空特征分析及分区建模 . 武汉大学学报 ● 信息科学版, 2012, 37(11): 1360-1363.
    [12] 章红平, 韩文慧, 黄玲, 耿长江.  地基GNSS全球电离层延迟建模 . 武汉大学学报 ● 信息科学版, 2012, 37(10): 1186-1189.
    [13] 李伟强, 杨东凯, 李明里, 张其善.  面向遥感的GNSS反射信号接收处理系统及实验 . 武汉大学学报 ● 信息科学版, 2011, 36(10): 1204-1208.
    [14] 卢鋆, 杨强文.  960MHz至1610MHz频段全球业务及GNSS频谱兼容分析 . 武汉大学学报 ● 信息科学版, 2011, 36(10): 1172-1176.
    [15] 刘志平, 何秀凤, 郭广礼, 查剑锋.  GNSS模糊度降相关算法及其评价指标研究 . 武汉大学学报 ● 信息科学版, 2011, 36(3): 257-261.
    [16] 蔡华, 赵齐乐, 孙汉荣, 胡志刚.  GNSS实时数据质量控制 . 武汉大学学报 ● 信息科学版, 2011, 36(7): 820-824.
    [17] 刘经南, 赵莹, 张小红.  GNSS无线电掩星电离层反演技术现状与展望 . 武汉大学学报 ● 信息科学版, 2010, 35(6): 631-635.
    [18] 张飞舟, 杨东凯, 陈嘉, 杨伯钢.  GNSS软件接收机的实现与测试分析 . 武汉大学学报 ● 信息科学版, 2009, 34(5): 577-580.
    [19] 李征航, 陈锴, 刘万科, 王文丽.  GNSS电离层延迟模型的数学统一与方法扩展 . 武汉大学学报 ● 信息科学版, 2007, 32(8): 699-703.
    [20] 耿涛, 丁志辉, 谢新, 吕逸飞.  基于载波相位差分的多频多GNSS测速精度评估 . 武汉大学学报 ● 信息科学版, 0, 0(0): 0-0. doi: 10.13203/j.whugis20200226
  • 加载中
图(4) / 表(3)
计量
  • 文章访问数:  1187
  • HTML全文浏览量:  117
  • PDF下载量:  131
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-01-31
  • 刊出日期:  2019-05-05

NTCM广播电离层模型性能评估

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

    国家自然科学基金 41674004

    作者简介:

    赵齐乐, 博士, 教授, 主要从事GNSS精密定轨定位数据处理方法研究。zhaoql@whu.edu.cn

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

摘要: 对于全球卫星导航系统单频用户而言,广播电离层模型是改正电离层误差的主要方法之一。对NTCM-BC模型、MNTCM-BC模型、NTCM-Klobuchar模型和传统Klobuchar模型的性能进行了对比分析。结果表明,在一个太阳活动周期内,NTCM-BC和MNTCM-BC模型计算的垂直总电子含量精度相当,且精度最高,NTCM-Klobuchar模型次之,Klobuchar模型最低。在低纬度测站,NTCM系列模型较Klobuchar模型最高能减少20 TECU左右的误差。与Klobuchar模型相比,采用NTCM-BC和MNTCM-BC模型进行伪距单点定位,高程精度提高约38%,特别是太阳活动高峰期间,低纬度测站最高能提高3 m左右。

English Abstract

赵齐乐, 杨超, 耿涛, 马福建, 张强. NTCM广播电离层模型性能评估[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
引用本文: 赵齐乐, 杨超, 耿涛, 马福建, 张强. NTCM广播电离层模型性能评估[J]. 武汉大学学报 ● 信息科学版, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
ZHAO Qile, YANG Chao, GENG Tao, MA Fujian, ZHANG Qiang. Performance Analysis of NTCM Ionospheric Broadcast Models[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
Citation: ZHAO Qile, YANG Chao, GENG Tao, MA Fujian, ZHANG Qiang. Performance Analysis of NTCM Ionospheric Broadcast Models[J]. Geomatics and Information Science of Wuhan University, 2019, 44(5): 655-660. doi: 10.13203/j.whugis20170348
  • 在全球卫星导航系统(Global Navigation Satellite System,GNSS)定位应用中,双频用户主要利用不同频率的载波相位、伪距进行线性组合来消除电离层误差[1]。但对于单频用户而言,通常采用广播电离层模型或者经验模型来改正电离层误差。

    目前,GNSS系统播发不同的广播电离层模型。美国GPS使用Klobuchar模型[2],其电离层误差改正率大概在50%左右;欧洲Galileo使用NeQuick模型[3],其电离层误差改正效果优于Klobuchar模型;北斗系统使用北斗Klobuchar模型,模型改正率能达到65%[4]。文献[5]设计了一种新的电离层改正方案,考虑了时间、地磁纬度、电离层峰值和太阳活动等因素,称之为NTCM经验模型。文献[6]把NTCM经验模型的12参数简化到9参数,称之为下一代广播电离层模型,即NTCM-BC模型。研究表明[7],NTCM-BC模型与NeQuick模型精度相当,但其模型结构更简单。文献[8]在此基础上,更加细致地刻画了电离层双峰区域,提出了MNTCM-BC模型。同期,文献[9]在NTCM经验模型的基础上,提出由Klobuchar模型驱动的NTCM-Klobuchar模型。

    本文从时域、空域和定位结果3个方面对Klobuchar模型和NTCM系列模型(NTCM-BC模型、MNTCM-BC模型和NTCM-Klobuchar模型)这4种广播电离层模型的性能进行分析评估。

    • Klobuchar模型是由导航电文播发的8个参数(α0~α3, β0~β3)构成的3阶表达式。其将晚间的电离层时延设置为常数D(5×10-9 s),白天倾斜方向上的时延Δτ设置成当地时间的余弦函数[9]

      $$ {\rm{\Delta }}\tau = \left[ {D + A{\rm{cos}}\left( {\frac{{2{\rm{ \mathsf{ π} }}\left( {t - {T_P}} \right)}}{P}} \right)} \right] \cdot {M_F} $$ (1)

      式中,$ A = \mathop {\mathop \sum\limits_{i = 0} }\limits^3 {\alpha _i}\phi _m^i;P = \mathop {\mathop \sum\limits_{i = 0} }\limits^3 {\beta _i}\phi _m^i$为穿刺点的地磁纬度;t为穿刺点的当地时间;TP=14;MF为映射函数,表达式为:

      $$ {M_F} = 1 + 16{(0.53 - {e_l})^3} $$ (2)

      式中,el为高度角。

    • NTCM-BC模型是由当地时间、地理纬度和地磁纬度、低纬度双峰方程构成的非线性表达式[6]

      $$ T_{{\rm{model}}}^{{\rm{vert}}} = {F_1} \cdot {F_2} \cdot {F_3} $$ (3)
      $$ \left\{ {\begin{array}{*{20}{l}} {{F_1} = {\rm{cos}}{\chi ^{{\rm{***}}}} + {\rm{cos}}{\chi ^{{\rm{**}}}}({c_1}{\rm{cos}}{V_D} + {c_2}{\rm{cos}}{V_{SD}} + }\\ {{c_3}{\rm{sin}}{V_{SD}} + {c_4}{\rm{cos}}{V_{TD}} + {c_5}{\rm{sin}}{V_{TD}})}\\ {{F_2} = 1 + {c_6}{\rm{cos}}{\phi _m}}\\ {{F_3} = {c_7} + {c_8}{\rm{exp}}\left( {{E_{{c_1}}}} \right) + {c_9}{\rm{exp}}\left( {{E_{{c_2}}}} \right)} \end{array}} \right. $$ (4)

      其中,

      $$ {\rm{cos}}{\chi ^{***}} = {\rm{cos}}\left( {\phi - \delta } \right) + 0.4 $$ (5)
      $$ {\rm{cos}}{\chi ^{**}} = {\rm{cos}}\left( {\phi - \delta } \right) - \frac{2}{{\rm{ \mathsf{ π} }}} \cdot \phi \cdot {\rm{sin}}\delta $$ (6)
      $$ {E_{{c_1}}} = - \frac{{{{({\phi _m} - {\phi _{{c_1}}})}^2}}}{{2\sigma _{{c_1}}^2}}, {\phi _{{c_1}}} = {16^ \circ }{\rm{N}}, {\sigma _{{c_1}}} = {12^ \circ } $$ (7)
      $$ {E_{{c_2}}} = - \frac{{{{({\phi _m} - {\phi _{{c_2}}})}^2}}}{{2\sigma _{{c_2}}^2}}, {\phi _{{c_2}}} = {10^ \circ }{\rm{S}}, {\sigma _{{c_2}}} = {13^ \circ } $$ (8)
      $$ \delta = 23.45{\rm{sin}}\left( {2{\rm{ \mathsf{ π} }} \cdot \frac{{284 + d}}{{365}}} \right) \cdot \frac{{\rm{ \mathsf{ π} }}}{{180}} $$ (9)

      式(3)~(9)中,Tmodelvert为电离层垂直总电子含量(vertical total electron content,VTEC);c1~c9为模型参数;VD=2π·(LT-14)/24;VSD=2π·LT/12;VTD=2π·LT/8;ϕϕm分别为地理纬度和地磁纬度;d为年积日;LT为当地时间。

      MNTCM-BC模型计算VTEC的方程与NTCM-BC模型基本相同,但式(7)和式(8)表达为σc1=σc2=18°-8°cosVD[8]

      NTCM-Klobuchar模型加上了季节变化方程和太阳活动方程[9]

      $$ T_{{\rm{model}}}^{{\rm{vert}}} = {F_1} \cdot {F_2} \cdot {F_3} \cdot {F_4} \cdot {F_5} $$ (10)

      其中,

      $$ \left\{ \begin{array}{l} {F_1} = {\rm{cos}}{\chi ^{{\rm{***}}}} + {\rm{cos}}{\chi ^{{\rm{**}}}}(0.879{\rm{}}09{\rm{cos}}{V_D} + \\ 0.174{\rm{}}66{\rm{cos}}{V_{SD}} - 0.025{\rm{sin}}{V_{SD}} + \\ 0.060{\rm{}}58{\rm{cos}}{V_{TD}} + 0.007{\rm{}}14{\rm{sin}}{V_{TD}})\\ {F_2} = 1 + 0.199{\rm{}}2{\rm{cos}}{V_A} - 0.026{\rm{}}34{\rm{cos}}{V_{SA}}\\ {F_3} = 1 - 0.318{\rm{}}36{\rm{cos}}{\phi _m}\\ {F_4} = 1 + 0.992{\rm{}}21{\rm{exp}}\left( {{E_{{c_1}}}} \right) + 1.016{\rm{}}12{\rm{exp}}\left( {{E_{{c_2}}}} \right)\\ {F_5} = 2.594{\rm{}}18 + 0.351{\rm{}}27K \end{array} \right. $$ (11)

      式中,

      $$ \left\{ {\begin{array}{*{20}{l}} {{V_A} = 2{\rm{ \mathsf{ π} }}\left( {d - 18} \right)/365.25}\\ {{V_{SA}} = 4{\rm{ \mathsf{ π} }}\left( {d - 6} \right)/365.25}\\ {\delta = 23.44{\rm{sin}}\left( {0.985{\rm{}}6 \cdot \left( {d - 80.7} \right) \cdot \frac{{\rm{ \mathsf{ π} }}}{{180}}} \right) \cdot \frac{{\rm{ \mathsf{ π} }}}{{180}}} \end{array}} \right. $$

      式(11)中的K是由Klobuchar模型计算的A点(10°N, 90°W)和B点(10°S, 90°W)的VTEC之和:

      $$ K = {T_A} + {T_B} $$ (12)
    • 本文分别从时域、空域和定位结果这3个角度来分析电离层模型的性能。在时域方面,采用欧洲定轨中心(The Center for Orbit Determination in Europe,CODE)提供的一个太阳活动周期(2005—2016年)的电离层格网产品作为参考值。采用前一天的CODE产品估算NTCM-BC或MNTCM-BC模型的9个参数,估算策略采用迭代非线性最小二乘拟合[7]。NTCM-Klobuchar模型参数采用文献[9]研究所得的经验系数。

      在空域方面,本文选取了全球不同纬度分布的7个国际GNSS服务(International GNSS Service,IGS)测站(见表 1),测站均坐落在东经120°左右,当地时间比世界时晚8 h左右。测站的电离层参考值采用双线性内插CODE电离层格网产品得到。

      表 1  IGS测站信息

      Table 1.  IGS Station Coordinates

      测站名 地理纬度 地理经度 城市 国家
      YAKT 62.030 8°N 129.680 2°E 雅库茨克 俄罗斯
      BJFS 39.608 6°N 115.892 2°E 北京 中国
      TCMS 24.797 7°N 120.987 2°E 新竹 中国
      PIMO 14.635 5°N 121.077 5°E 奎松市 菲利宾
      DARW 12.843 6°S 131.132 5°E 达尔文 澳大利亚
      PERT 31.801 7°S 115.885 0°E 珀斯 澳大利亚
      MAC1 54.499 4°S 158.935 6°E 麦夸里岛 澳大利亚

      评价指标采用均方根(root mean square,RMS)R,计算公式为:

      $$ R = \sqrt { < {{({V_{{\rm{mod}}}} - {V_{{\rm{ref}}}})}^2} > } $$ (13)

      式中,Vmod为4个模型的电离层计算值;Vref为电离层参考值;符号 < > 表示对累积和取平均。

      为了忽略卫星轨道和钟差误差的影响,以评估电离层改正模型的定位性能,本文采用精密星历对上述7个测站进行伪距单点定位。观测文件、精密星历、钟差等文件均从网站ftp://cddis.gsfc.nasa.gov下载,测站的参考坐标均从IGS产品sinex文件中获取。

    • 图 1给出了2005—2016年电离层模型RMS(10 d取一均值)的统计结果,其中N-Klo代表NTCM-Klobuchar模型,N-BC代表NTCM-BC模型,MN-BC代表MNTCM-BC模型。由图 1可知,在太阳活动平静期(2005—2010年),4个模型的RMS值在10 TECU以下,在2009年低至7 TECU;在太阳活动活跃期(2011—2016年),模型精度降低,特别是太阳活动高峰年(2014年),Klobuchar模型的RMS值能高达20 TECU。

      图  1  2005—2016年电离层模型的RMS值时序图

      Figure 1.  RMS Values of Ionospheric Models forthe Year of 2005—2016

      总体来看,NTCM-BC模型和MNTCM-BC模型精度相当,且精度最高,Klobuchar模型最低,NTCM-Klobuchar模型介于二者之间。NTCM-BC和MNTCM-BC模型精度在2005—2010年比Klobuchar模型提高了3~4 TECU,在2011—2016年提高了3~10 TECU。

      从以上分析可以看出,电离层模型精度受太阳活动影响。为了更好地分析太阳活动与模型精度的关系,图 2给出了12 a期间每天电离层模型精度随太阳活动指数F10.7(单位为sfu,即太阳辐射通量)的变化关系。整体来看,4个模型的误差与太阳活动水平呈现一致性的变化特征,即随着F10.7指数的变大(即太阳活动增强),RMS值也随之变大,模型精度变差。

      图  2  2005—2016年每天电离层模型RMS值随F10.7指数的变化

      Figure 2.  Daily RMS Values of the Ionospheric Models for the Year of 2005—2016 with F10.7 Index

      图 2进一步发现,Klobuchar模型的斜率最大,NTCM-Klobuchar模型次之,NTCM-BC和MNTCM-BC模型最小。总体来看,F10.7指数在60~100之间时,NTCM-BC和MNTCM-BC模型的RMS值比Klobuchar模型小3~5 TECU,NTCM-Klobuchar模型的RMS值比Klobuchar模型小1~2 TECU;F10.7指数大于120时,Klobuchar模型的RMS值为5~22 TECU,NTCM-Klobuchar模型的RMS值为3~20 TECU,NTCM-BC和MNTCM-BC模型的RMS值均在15 TECU之下;在F10.7指数大于160时MNTCM-BC模型的RMS值比NTCM-BC模型低0.5 TECU左右。

    • 纬度是影响电离层误差的重要因素之一。为了进一步评估电离层模型在不同太阳活动期间和不同纬度分布的改正精度,本文分别在太阳活动低峰年和高峰年选取一天,图 3图 4中左、右列分别代表低峰年2009年12月9日和高峰年2014年9月12日7个测站每小时VTEC的RMS值。

      图  3  北纬4个IGS测站的电离层模型每小时RMS值时序图

      Figure 3.  Hourly RMS Values of VTEC of the Ionospheric Models at 4 IGS Stations in Northern Latitude

      图  4  南纬3个IGS测站的电离层模型每小时RMS值时序图

      Figure 4.  Hourly RMS Values of VTEC of the Ionospheric Models at 3 IGS Stations in Southern Latitude

      图 3图 4可知,电离层模型的RMS值与当地时间密切相关。在太阳活动平静期(2009年12月9日),NTCM-BC和MNTCM-BC模型的RMS值在5 TECU以内,且随当地时间变化平稳,而NTCM-Klobuchar和Klobuchar模型在PIMO和TCMS测站的RMS值可达10~15 TECU。在太阳活动活跃期(2014年9月12日),4个模型RMS值在当地时间10—20时普遍大于其他时刻,且相比于太阳活动平静期明显增大。

      图 3图 4还可知,电离层模型在低纬度测站的改正误差大于中、高纬度测站。在太阳活动活跃期间,NTCM系列模型能显著降低电离层误差。Klobuchar模型在TCMS和PIMO测站的RMS值高达45 TECU,而NTCM系列模型在同样时间段仅有20 TECU。4个模型在低纬度测站的误差值呈现双峰现象,且NTCM系列模型刻画得更为明显。NTCM-Klobuchar模型在太阳活动平静期的TCMS和PIMO测站表现较差,可能是由于该模型采用经验系数所致[9]

    • 本文基于4个广播电离层模型,分别选取2009年12月9日(太阳活动平静期)和2014年9月12日(太阳活动活跃期)7个IGS测站(见表 1)的观测数据进行单历元伪距单点定位,并对全天解算的所有历元的水平和高程精度进行统计,结果如表 2表 3所示。

      表 2  2009年12月9日的定位精度/m

      Table 2.  RMS of Position Error on December 9th, 2009/m

      测站名 Klobuchar NTCM-Klobuchar NTCM-BC MNTCM-BC
      水平 高程 水平 高程 水平 高程 水平 高程
      YAKT 1.42 3.17 1.19 1.83 1.27 1.78 1.28 1.77
      BJFS 1.09 3.71 1.23 2.82 1.04 2.04 1.04 2.04
      TCMS 1.27 3.35 1.09 3.45 1.17 2.01 1.23 2.03
      PIMO 1.04 3.17 1.01 3.57 1.05 2.36 1.06 2.39
      DARW 0.99 2.93 0.92 2.18 0.87 1.47 0.88 1.50
      PERT 1.31 3.15 1.39 2.21 1.43 2.01 1.48 2.04
      MAC1 1.12 2.07 1.04 1.85 1.04 1.66 1.06 1.67
      平均值 1.18 3.08 1.12 2.56 1.12 1.90 1.15 1.92

      表 3  2014年9月12日的定位精度/m

      Table 3.  RMS of Position Error on September 12th, 2014/m

      测站名 Klobuchar NTCM-Klobuchar NTCM-BC MNTCM-BC
      水平 高程 水平 高程 水平 高程 水平 高程
      YAKT 1.28 3.08 1.22 1.48 1.26 1.56 1.15 1.65
      BJFS 1.92 2.20 1.44 1.85 1.47 2.13 1.38 2.50
      TCMS 6.18 4.77 5.40 3.21 5.52 3.29 5.59 3.24
      PIMO 2.45 5.76 2.66 3.07 2.43 3.11 2.46 3.14
      DARW 5.00 6.90 3.41 3.62 3.56 3.64 3.47 3.65
      PERT 1.12 3.04 1.04 2.09 1.07 2.04 1.00 2.10
      MAC1 1.27 2.10 1.29 1.80 1.30 1.99 1.30 1.94
      平均值 2.75 3.98 2.35 2.45 2.37 2.54 2.34 2.60

      表 2表 3中可以发现:在水平方向上,采用NTCM系列模型的定位精度差异在0.3 m以内,略优于Klobuchar模型。在高程方向上,在太阳活动平静期,采用NTCM-BC和MNTCM-BC模型的定位精度较Klobuchar模型能提高约1 m,较NTCM-Klobuchar模型能提高约0.6 m;在太阳活动活跃期,采用NTCM系列模型的定位精度相当,在低纬度测站比Klobuchar模型最高能提高约3 m。

    • 本文介绍了最新的广播电离层改正模型,并从时域、空域、定位等不同角度评估了Klobuchar模型、NTCM-BC模型、MNTCM-BC模型和NTCM-Klobuchar模型的性能。通过实验得出:

      1)4个电离层模型计算的VTEC误差值均表现出随太阳活动变化的特征。

      2)从一个太阳活动周期的电离层模型精度来看,NTCM-BC和MNTCM-BC模型精度相当,且精度最高,NTCM-Klobuchar模型次之,Klobuchar模型最低。

      3)4个模型均呈现低纬度地区误差峰值特征。与Klobuchar模型相比,NTCM系列模型的RMS值在低纬度测站最高能减少20 TECU。

      4)从采用精密星历的伪距单点定位结果来看,在水平方向上,采用NTCM系列模型的定位精度相差不大,略优于Klobuchar模型;但在高程方向,采用NTCM-BC和MNTCM-BC模型的定位精度较Klobuchar模型提高了约38%,特别是太阳活动高峰期的低纬度测站,最高能提高3 m左右。

参考文献 (9)

目录

    /

    返回文章
    返回