留言板

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

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

单站区域电离层TEC建模及精度分析

李涌涛 赵昂 李建文 车通宇 潘林 陈晨

李涌涛, 赵昂, 李建文, 车通宇, 潘林, 陈晨. 单站区域电离层TEC建模及精度分析[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
引用本文: 李涌涛, 赵昂, 李建文, 车通宇, 潘林, 陈晨. 单站区域电离层TEC建模及精度分析[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
LI Yongtao, ZHAO Ang, LI Jianwen, CHE Tongyu, PAN Lin, CHEN Chen. Regional Ionospheric TEC Modeling and Accuracy Analysis Based on Observations from a Station[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
Citation: LI Yongtao, ZHAO Ang, LI Jianwen, CHE Tongyu, PAN Lin, CHEN Chen. Regional Ionospheric TEC Modeling and Accuracy Analysis Based on Observations from a Station[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286

单站区域电离层TEC建模及精度分析

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

国家自然科学基金 41904042

详细信息
    作者简介:

    李涌涛,硕士,工程师,主要研究方向为电离层TEC建模与GNSS数据处理。641594920@qq.com

    通讯作者: 李建文,博士,教授。ZZLJW@126.com
  • 中图分类号: P237

Regional Ionospheric TEC Modeling and Accuracy Analysis Based on Observations from a Station

Funds: 

The National Natural Science Foundation of China 41904042

More Information
    Author Bio:

    LI Yongtao, master, engineer, specializes in ionospheric TEC modeling and GNSS data analysis. E-mail: 641594920@qq.com

    Corresponding author: LI Jianwen, PhD, professor. E-mail: ZZLJW@126.com
  • 摘要: 为了分析单站区域电离层总电子含量(total electron content, TEC)模型的适用范围和精度,基于2~‍15阶次球谐函数,分别建立了欧洲区域16个单站区域电离层TEC模型,生成了区域格网TEC,并与欧洲定轨中心(Center for Orbit Determination in Europe,CODE)、国际全球导航卫星系统服务组织(International Global Navigation Satellite System Service, IGS)和全球GNSS监测评估系统(International GNSS Monitoring and Assessment System, iGMAS)等全球电离层产品比较分析。结果表明,基于低阶(2×2阶或3×3阶)球谐函数建立的单站区域电离层TEC模型,以单站为中心,在经纬度10°×10°范围以内(半径小于600 km),电离层TEC精度与CODE、iGMAS和IGS等全球电离层产品的TEC精度相当,约为1.0 TECU,实现了在一定区域内(半径小于600 km)利用单站建立电离层TEC模型替代全球电离层TEC模型,高效地为区域内单频用户提供高精度的电离层延迟改正。
  • 图  1  格网经纬度范围

    Figure  1.  Latitude and Longitude Range of Grid

    图  2  欧洲区域16个GNSS跟踪站分布图

    Figure  2.  Map of 16 GNSS Tracking Stations in the European Region

    图  3  2017年DOY第357~363天太阳F10.7指数

    Figure  3.  Solar F10.7 Index from DOY 357 to 363 in 2017

    图  4  2017年DOY第357~363天地磁Kp指数

    Figure  4.  Geomagnetic Kp Index from DOY 357 to 363 in 2017

    图  5  16个测站0°×0°范围7天的TEC平均值

    Figure  5.  TEC Mean Values of 16 Stations in the Range of 0°×0° for 7 Days

    图  6  bzrg站0°×0°范围7天的TEC平均值和RMS平均值

    Figure  6.  Mean Values of TEC and RMS for 7 Days in the Range of 0°×0° at bzrg Station

    图  7  wtzr站0°×0°范围7天的TEC平均值和RMS平均值

    Figure  7.  Mean Values of TEC and RMS for 7 Days in the Range of 0°×0° at wtzr Station

    图  8  16个测站不同范围7天的TEC平均值

    Figure  8.  TEC Mean Values of 16 Stations in Different Ranges for 7 Days

    图  9  pots站不同范围7天的TEC和RMS平均值

    Figure  9.  Mean Values of TEC and RMS of Different Ranges for 7 Days at pots Station

    图  10  ganp测站不同范围7天的TEC和RMS平均值

    Figure  10.  Mean Values of TEC and RMS of Different Ranges for 7 Days at ganp Station

    图  11  pots站观测的GPS卫星星下点轨迹图

    Figure  11.  GPS Sub Satellite Point Trajectory Map Observed at pots Station

    图  12  IGS和ganp单站2~3阶次区域电离层TEC云图

    Figure  12.  TEC Cloud Images of Second-and-Third-Order Ionosphere of IGS and ganp Stations

    图  13  IGS和pots单站2~3阶次区域电离层TEC云图

    Figure  13.  TEC Cloud Images of Second-and-Third-Order Ionosphere of IGS and pots Stations

  • [1] 毛文飞. 电离层TEC预测模型的研究及应用[D]. 抚州: 东华理工大学, 2018

    Mao Wenfei. Research and Application of Prediction Model of Total Electron Content in Ionosphere[D]. Fuzhou: East China University of Technology, 2018
    [2] 王成, 毛大智, 施闯, 等. 全球电离层模型的分布式并行解算[J]. 武汉大学学报·信息科学版, 2018, 43(8): 1207-1213 doi:  10.13203/j.whugis20160274

    Wang Cheng, Mao Dazhi, Shi Chuang, et al. Distributed Parallel Estimation for Global Ionospheric Modeling[J]. Geomatics and Information Science of Wuhan University, 2018, 43(8): 1207-1213 doi:  10.13203/j.whugis20160274
    [3] 薛军琛, 宋淑丽, 朱文耀. 基于BDS/GPS双系统的全球电离层建模[J]. 中国科学: 物理学力学天文学, 2015, 45(7): 28-37 https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK201507005.htm

    Xue Junchen, Song Shuli, Zhu Wenyao. Global Ionosphere Model Based on BDS/GPS Dual-System Observations[J]. Scientia Sinica(Physica, Mechanica and Astronomica), 2015, 45(7): 28-37 https://www.cnki.com.cn/Article/CJFDTOTAL-JGXK201507005.htm
    [4] 李征航, 黄劲松. GPS测量与数据处理[M]. 武汉: 武汉大学出版社, 2005

    Li Zhenghang, Huang Jinsong. GPS Surveying and Data Processing[M]. Wuhan: Wuhan University Press, 2005
    [5] 袁运斌. 基于GPS的电离层监测及延迟改正理论与方法的研究[D]. 武汉: 中国科学院研究生院, 2002

    Yuan Yunbin. Study on Theories and Methods of Correcting Ionospheric Delay and Monitoring Ionosphere Based on GPS[D]. Wuhan: Graduate School of Chinese Academy of Sciences, 2002
    [6] 刘磊, 姚宜斌, 孔建, 等. 附加GIM约束的全球电离层建模[J]. 大地测量与地球动力学, 2017, 37 (1): 67-71 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201701015.htm

    Liu Lei, Yao Yibin, Kong Jian, et al. Global Ionosphere Modeling with GIM Constraint[J]. Journal of Geodesy and Geodynamics, 2017, 37(1): 67-71 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201701015.htm
    [7] 李子申, 王宁波, 李敏, 等. 国际GNSS服务组织全球电离层TEC格网精度评估与分析[J]. 地球物理学报, 2017, 60(10): 3718-3729 doi:  10.6038/cjg20171003

    Li Zishen, Wang Ningbo, Li Min, et al. Evaluation and Analysis of the Global Ionospheric TEC Map in the Frame of International GNSS Services[J]. Chinese Journal of Geophysics, 2017, 60(10): 3718-3729 doi:  10.6038/cjg20171003
    [8] 袁运斌, 霍星亮, 张宝成. 近年来我国GNSS电离层延迟精确建模及修正研究进展[J]. 测绘学报, 2017, 46(10): 1364-1378 doi:  10.11947/j.AGCS.2017.20170349

    Yuan Yunbin, Huo Xingliang, Zhang Baocheng. Research Progress of Precise Models and Correction for GNSS Ionospheric Delay in China over Recent Years[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1364-1378 doi:  10.11947/j.AGCS.2017.20170349
    [9] 杨业, 谢益炳, 陆亚峰, 等. 基于陆态网数据低阶球谐函数电离层区域建模[J]. 测绘地理信息, 2015, 40 (3): 45-47

    Yang Ye, Xie Yibing, Lu Yafeng, et al. Modeling the Low-Order Spherical Harmonic Function Ionospheric Region Based on the Land State Network Data[J]. Journal of Geomatics, 2015, 40(3): 45-47
    [10] 柳景斌, 王泽民, 王海军, 等. 利用球冠谐分析方法和GPS数据建立中国区域电离层TEC模型[J]. 武汉大学学报·信息科学版, 2008, 33(8): 792-795 http://ch.whu.edu.cn/article/id/1682

    Liu Jingbin, Wang Zemin, Wang Haijun, et al. Modeling Regional Ionosphere Using GPS Measurements over China by Spherical Cap Harmonic Analysis Methodology[J]. Geomatics and Information Science of Wuhan University, 2008, 33(8): 792-795 http://ch.whu.edu.cn/article/id/1682
    [11] 阿尔察, 刘四清, 黄文耿, 等. 中国电离层TEC同化现报系统[J]. 地球物理学报, 2018, 61(6): 2186-2197 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201806002.htm

    A Ercha, Liu Siqing, Huang Wengeng, et al. Ionospheric TEC Data Assimilation and Now-Casting System over China[J]. Chinese Journal of Geophysics, 2018, 61(6): 2186-2197 https://www.cnki.com.cn/Article/CJFDTOTAL-DQWX201806002.htm
    [12] 付伟, 姚顽强, 陈鹏, 等. 基于多项式的区域电离层建模精度分析与验证[J]. 北京测绘, 2018, 32(10): 1219-1225 https://www.cnki.com.cn/Article/CJFDTOTAL-BJCH201810024.htm

    Fu Wei, Yao Wanqiang, Chen Peng, et al. Polynomial Based Accuracy Analysis and Validation of Regional Ionospheric Modeling[J]. Beijing Surveying and Mapping, 2018, 32(10): 1219-1225 https://www.cnki.com.cn/Article/CJFDTOTAL-BJCH201810024.htm
    [13] 杨帆, 李宁. 区域电离层延迟模型研究[J]. 测绘与空间地理信息, 2018, 41(7): 15-18 doi:  10.3969/j.issn.1672-5867.2018.07.005

    Yang Fan, Li Ning. Study on Regional Ionospheric Delay Model[J]. Geomatics & Spatial Information Technology, 2018, 41(7): 15-18 doi:  10.3969/j.issn.1672-5867.2018.07.005
    [14] Farzaneh S, Forootan E. Reconstructing Regional Ionospheric Electron Density: A Combined Spherical Slepian Function and Empirical Orthogonal Function Approach[J]. Surveys in Geophysics, 2018, 39 (2): 289-309 doi:  10.1007/s10712-017-9446-y
    [15] 周要宗, 匡翠林, 窦邵华. 高阶电离层延迟对中国区域双差定位影响研究[J]. 大地测量与地球动力学, 2019, 39(5): 506-510 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201905013.htm

    Zhou Yaozong, Kuang Cuilin, Dou Shaohua. Analysing Higher-Order Ionospheric Effects on GNSS Double-Differenced Positioning in China Region[J]. Journal of Geodesy and Geodynamics, 2019, 39 (5): 506-510 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201905013.htm
    [16] 陆仲连. 球谐函数[M]. 广州: 解放军出版社, 1988

    Lu Zhonglian. Spherical Harmonic Function[M]. Guangzhou: Liberation Army Publishing House, 1988
    [17] 王健, 党亚民, 王虎. 全球电离层延迟建模及精度分析[J]. 测绘工程, 2019, 28(2): 17-21 https://www.cnki.com.cn/Article/CJFDTOTAL-CHGC201902004.htm

    Wang Jian, Dang Yamin, Wang Hu. Modeling Global Ionospheric Delay and Its Accuracy Analysis[J]. Engineering of Surveying and Mapping, 2019, 28(2): 17-21 https://www.cnki.com.cn/Article/CJFDTOTAL-CHGC201902004.htm
    [18] 畅鑫, 张伟. 实时区域电离层TEC建模与单频PPP实验[J]. 大地测量与地球动力学, 2013, 33(5): 39-43 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201305009.htm

    Chang Xin, Zhang Wei. Real Time Regional Ionospheric TEC Modeling and Single-Frequency PPP Experiment[J]. Journal of Geodesy and Geodynamics, 2013, 33(5): 39-43 https://www.cnki.com.cn/Article/CJFDTOTAL-DKXB201305009.htm
    [19] International GNSS Service. IGS Products[EB/OL]. [2007-04-15]. http://www.igs.org/products
    [20] International GNSS Monitoring and Assessment System. iGMAS Products[EB/OL]. [2014-12-13]. http://www.igmas.org/Product/Cpdetail/detail/nav_id/4/cate_id/36.html
  • [1] 刘立龙, 陈雨田, 黎峻宇, 田祥雨, 贺朝双.  活跃期区域电离层总电子短期预报及适用性分析 . 武汉大学学报 ● 信息科学版, 2019, 44(12): 1757-1764. doi: 10.13203/j.whugis20180145
    [2] 王成, 施闯, 章红平.  全球电离层模型的OpenMP多线程并行解算 . 武汉大学学报 ● 信息科学版, 2018, 43(6): 813-819. doi: 10.13203/j.whugis20150223
    [3] 王成, 毛大智, 施闯, 章红平.  全球电离层模型的分布式并行解算 . 武汉大学学报 ● 信息科学版, 2018, 43(8): 1207-1213. doi: 10.13203/j.whugis20160274
    [4] 郑凯, 郭博峰, 张小红.  接收机钟跳对单站GPS测速的影响及改正方法 . 武汉大学学报 ● 信息科学版, 2017, 42(3): 304-308, 327. doi: 10.13203/j.whugis20150119
    [5] 邹璇, 李宗楠, 陈亮, 宋伟伟, 王成, 唐卫明.  一种历元间差分单站单频周跳探测与修复方法 . 武汉大学学报 ● 信息科学版, 2017, 42(10): 1406-1410. doi: 10.13203/j.whugis20150805
    [6] 沈迎春, 闫昊明, 彭鹏, 白希选, 田代恒.  质量负荷引起地表形变的格林函数和球谐函数方法对比研究 . 武汉大学学报 ● 信息科学版, 2017, 42(7): 1008-1014. doi: 10.13203/j.whugis20150201
    [7] 舒宝, 刘晖, 张明, 吴丹.  北斗系统硬件延迟解算及精度分析 . 武汉大学学报 ● 信息科学版, 2016, 41(2): 279-284. doi: 10.13203/j.whugis20140058
    [8] 谢益炳, 伍吉仓, 陈俊平, 刘伟洲.  GPS和GLONASS组合的全球实测电离层TEC建模 . 武汉大学学报 ● 信息科学版, 2014, 39(8): 930-934. doi: 10.13203/j.whugis20130128
    [9] 曾艳艳, 于锦海, 万晓云.  任意球冠下Stokes-Neumman混合边值问题的球谐函数有限逼近方法 . 武汉大学学报 ● 信息科学版, 2014, 39(1): 65-69.
    [10] 聂文锋, 胡伍生, 潘树国, 宋玉兵.  利用GPS双频数据进行区域电离层TEC提取 . 武汉大学学报 ● 信息科学版, 2014, 39(9): 1022-1027. doi: 10.13203/j.whugis20130046
    [11] 王成, 王解先, 段兵兵.  附有国际参考电离层约束的全球电离层模型 . 武汉大学学报 ● 信息科学版, 2014, 39(11): 1340-1346.
    [12] 谢益炳, 陈俊平, 伍吉仓, 胡丛玮.  不同约束条件对电离层电子含量和硬件延迟的影响 . 武汉大学学报 ● 信息科学版, 2014, 39(7): 799-803.
    [13] 王建强, 李建成, 王正涛, 赵国强.  球谐函数变换快速计算扰动引力 . 武汉大学学报 ● 信息科学版, 2013, 38(9): 1039-1043.
    [14] 刘志平, 赵自强, 郭广礼.  电离层总电子含量时空特征分析及分区建模 . 武汉大学学报 ● 信息科学版, 2012, 37(11): 1360-1363.
    [15] 苏勇, 范东明, 游为.  缔合Legendre函数二阶导数的快速稳定递推算法 . 武汉大学学报 ● 信息科学版, 2012, 37(12): 1409-1412.
    [16] 章红平, 韩文慧, 黄玲, 耿长江.  地基GNSS全球电离层延迟建模 . 武汉大学学报 ● 信息科学版, 2012, 37(10): 1186-1189.
    [17] 张小红, 程世来, 李星星, 郭斐.  单站GPS载波平滑伪距精密授时研究 . 武汉大学学报 ● 信息科学版, 2009, 34(4): 463-465.
    [18] 李征航, 陈锴, 刘万科, 黄欢.  顾及f~3项的电离层延迟模型 . 武汉大学学报 ● 信息科学版, 2007, 32(2): 139-143.
    [19] 蔡昌盛, 高井祥, 李征航.  利用GPS监测电离层总电子含量的季节性变化 . 武汉大学学报 ● 信息科学版, 2006, 31(5): 451-453.
    [20] 张传定, 吴晓平.  非心摄动引力的快速计算方法研究 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 87-90.
  • 加载中
图(13)
计量
  • 文章访问数:  659
  • HTML全文浏览量:  151
  • PDF下载量:  77
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-24
  • 刊出日期:  2022-01-05

单站区域电离层TEC建模及精度分析

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

    国家自然科学基金 41904042

    作者简介:

    李涌涛,硕士,工程师,主要研究方向为电离层TEC建模与GNSS数据处理。641594920@qq.com

    通讯作者: 李建文,博士,教授。ZZLJW@126.com
  • 中图分类号: P237

摘要: 为了分析单站区域电离层总电子含量(total electron content, TEC)模型的适用范围和精度,基于2~‍15阶次球谐函数,分别建立了欧洲区域16个单站区域电离层TEC模型,生成了区域格网TEC,并与欧洲定轨中心(Center for Orbit Determination in Europe,CODE)、国际全球导航卫星系统服务组织(International Global Navigation Satellite System Service, IGS)和全球GNSS监测评估系统(International GNSS Monitoring and Assessment System, iGMAS)等全球电离层产品比较分析。结果表明,基于低阶(2×2阶或3×3阶)球谐函数建立的单站区域电离层TEC模型,以单站为中心,在经纬度10°×10°范围以内(半径小于600 km),电离层TEC精度与CODE、iGMAS和IGS等全球电离层产品的TEC精度相当,约为1.0 TECU,实现了在一定区域内(半径小于600 km)利用单站建立电离层TEC模型替代全球电离层TEC模型,高效地为区域内单频用户提供高精度的电离层延迟改正。

English Abstract

李涌涛, 赵昂, 李建文, 车通宇, 潘林, 陈晨. 单站区域电离层TEC建模及精度分析[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
引用本文: 李涌涛, 赵昂, 李建文, 车通宇, 潘林, 陈晨. 单站区域电离层TEC建模及精度分析[J]. 武汉大学学报 ● 信息科学版, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
LI Yongtao, ZHAO Ang, LI Jianwen, CHE Tongyu, PAN Lin, CHEN Chen. Regional Ionospheric TEC Modeling and Accuracy Analysis Based on Observations from a Station[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
Citation: LI Yongtao, ZHAO Ang, LI Jianwen, CHE Tongyu, PAN Lin, CHEN Chen. Regional Ionospheric TEC Modeling and Accuracy Analysis Based on Observations from a Station[J]. Geomatics and Information Science of Wuhan University, 2022, 47(1): 69-78. doi: 10.13203/j.whugis20190286
  • 电离层大气分子是受太阳高能辐射及宇宙射线激发电离产生的高层大气,距地高度约60~1 000 km,电离层是日地空间的重要组成部分和研究对象,电离层总电子含量(total electron content,TEC)是描述电离层性质的重要物理参数[1-2]。电磁波信号穿过电离层时,受电离层影响,传播路径和速度发生改变,成为卫星导航中重要的误差源之一[3-4]。为了提高导航定位精度,需要对观测量进行电离层延迟改正,以国际全球导航卫星系统服务组织(International Global Navigation Satellite System Service,IGS)为代表的全球电离层云图(global ionosphere maps,GIM)是用于电离层延迟改正的常用方法之一[5-6],GIM以格网的形式给出了全球范围内电离层总电子含量TEC及相应的均方根(root mean square,RMS)值,适用于全球范围电离层延迟改正。对于区域用户而言,如以1 000 km半径作为活动范围的用户,GIM中区域以外的TEC信息没有用到而且也没有必要,相比于GIM,利用尽可能少的测站建立区域内电离层TEC模型具有更高的效率和更适用的范围,因此,建立单站区域电离层TEC模型,分析其适用范围和精度对区域内用户提高导航定位精度尤为重要。研究人员分析了IGS及欧洲定轨中心(Center for Orbit Determination in Europe,CODE)、中国科学院、美国喷气推进实验室(Jet Propulsion Laboratory,JPL)和西班牙加泰罗尼亚理工大学等分析中心的GIM精度,以及近年来中国GNSS电离层延迟建模的进展,为电离层TEC精度分析提供了一定的参考[7-8];基于球谐函数和球冠谐函数建立了中国区域电离层TEC模型,提供了一定的电离层区域建模参考[9-10];实现了中国电离层TEC同化现报系统[11];基于多项式曲面拟合建立并分析了区域电离层的模型精度[12]。以上研究的区域范围较广,所用数据模型也是高阶次模型,而基于单站低阶次球谐函数建立区域电离层TEC模型的相关研究较少。

    本文基于2~15不同阶次球谐函数分别建立了欧洲区域16个GNSS跟踪站单站区域电离层TEC模型,并与CODE、iGMAS和IGS的全球电离层TEC格网产品进行了对比,分析了单站区域电离层TEC模型的适用范围和精度。

    • 在GNSS信号传播过程中,假设沿信号传播路径上的所有电子都集中在无厚度的理想薄层上,薄层高度取450 km,并基于薄层假设利用单站全球定位系统(global positioning system,GPS)双频观测数据,采用载波相位平滑伪距的方法计算传播斜路径上的倾斜电离层总电子含量(slant total electronic content,STEC)值,根据投影映射函数将STEC转化为测站上方的垂直电离层总电子含量(vertical total electron content,VTEC),并与球谐函数模型联立,求解出相应的球谐函数系数,最终将球谐函数系数转化为区域电离层TEC格网文件[13]

      双频接收机的测距码和载波相位观测方程[914]为:

      Pk,ij=ρij+c×(dti-dtj)+Ik,ij+Tij+ c×rk,i+c×skj+εP,k,ijLk,ij=ρij+c×(dti-dtj)-Ik,ij+Tij+ εp,k,ij-λkNk,ij

      式中,k为频率编号;ji分别为卫星编号和接收机编号;PL分别为测距码和载波相位观测值;ρ为接收机和卫星间的几何距离;c为真空中的光速;dtidtj分别为接收机钟差和卫星钟差;TI分别为信号传播路径上的对流层延迟和电离层延迟;rs分别为接收机延迟偏差和卫星硬件延迟偏差;ε为测距码的观测噪声;λ为波长;N为整周模糊度。分别对P1P2L1L2进行差分,可得P4L4组合:

      P4=P1,ij-P2,ij= (I1,ij-I2,ij)+cDj+cDiL4=L1,ij-L2,ij=-(I1,ij-I2,ij)-λ(N1,ij-N2,ij)

      式中,Dj=s1j-s2jDi=r1,i-r2,i分别为卫星码偏差和接收机码偏差。采用相位观测量L4对伪码观测量P4进行平滑,伪距测量中电离层延迟量为:

      I=40.28f2δSTEC

      式中,I为电离层延迟量;f为载波相位频率;δSTEC为观测路径上倾斜电离层TEC。将式(3)代入式(2)可得:

      P4=40.28(1f12-1f22)δSTEC+cDj+cDiL4=-40.28(1f12-1f22)δSTEC-λ(N1,ij-N2,ij)

      为了提高伪距观测值的精度,采用Hatch滤波进行载波相位平滑伪距,对于连续观测弧段,平滑方式为:

      P4,st=k-1k[P4,st-1+L4t- L4t-1]+1kP4t,t>1P4,s1=P41,t=1

      式中,P4,s为平滑后的伪距观测值;t为历元;k为从平滑开始的总历元数。

      P4,s代入式(4)中替代P4可得:

      P4,s=40.28(1f12-1f22)δSTEC+cDj+cDi

      式中,f1f2分别为载波相位L1L2的频率。卫星和接收机硬件延迟偏差在一天内可视为无变化,基于减少待估参数和提高建模精度考虑,本文采用CODE解算的卫星硬件延迟偏差参与建模。

      将STEC通过单层模型投影函数最终转化为测站上空垂直电离层VTEC[15],转化过程为:

      F(z)=δSTECδVTEC=11-(ReRe+Hi·sinz)2δVTEC=f12f2240.28(ReRe+Hi·sinz)1-(ReRe+Hi·sinz)2× (ρ2'-ρ1')

      式中,F(z)为穿刺点处的投影函数;δVTEC为测站上空垂直电离层总电子含量;z为接收机天顶距;Re为地球半径;Hi为薄层高度。

      球谐函数是研究全球范围场力的重要数学模型工具[16],基于电离层薄层假设利用球谐函数进行单站区域电离层TEC建模,即用较少的参数很好地拟合区域电离层TEC的整体情况[417]。球谐函数为:

      δVTEC(φ,λ)=n=0nmaxm=0nP˜nmsinφ×A˜nmcos(mλ)+ B˜nmsin(mλ)

      式中,φλ分别为穿刺点处的纬度和经度;P˜nmsinφ为正规化后的nm次勒让德函数;A˜mnB˜mn分别为球谐函数模型系数。球谐函数阶次越低,所求系数越少,解算效率越高,但建模范围也有限,即低阶球谐函数在进行大范围建模时拟合结果较差,系数个数Np为:

      Np=nmax+12-nmax-mmaxnmax-mmax+1

      将式(7)中第2式和式(8)联立,利用最小二乘拟合求解球谐函数模型系数A˜mnB˜mn,通过函数转换,将球谐系数转换为通用电离层TEC格网格式[18]

    • CODE发布的电离层TEC最终产品精度为2~8 TECU,时间滞后约7 d[7],是利用全球300多个IGS跟踪站,基于15×15阶球谐函数建立全球电离层TEC模型解算得到的;IGS发布的电离层TEC最终产品精度为2~8 TECU,时间滞后为12 d[19],是综合了CODE和JPL两个分析中心的电离层产品得到的;国内iGMAS发布的电离层TEC最终产品精度为2~8 TECU,时间滞后约12 d[20],是由8个iGMAS分析中心的电离层产品加权得到的。因为这3个机构发布的电离层产品比较权威,所以用其作为电离层TEC精度比较分析的基准。

      利用单个测站GNSS数据建立电离层TEC模型,并生成电离层TEC格网转换(ionosphere map exchange format,IONEX)格式文件,与CODE、IGS和iGMAS的GIM格式相同,即将测区按纬度和经度分辨率为2.5°×5°格网化,给出了每个格网点的电离层TEC值和相应的RMS值,每天分为采样间隔为2 h的13幅电离层TEC格网图。取以观测站为中心,格网经纬度范围分别为0°×0°、10°×5°、10°×10°、10°×15°、20°×15°和20°×20°内所有格网点上电离层TEC的平均值和RMS的平均值,并分别与CODE、iGMAS和IGS的电离层GIM产品中相应的TEC的平均值和RMS的平均值进行对比分析,格网经纬度范围如图 1所示。

      图  1  格网经纬度范围

      Figure 1.  Latitude and Longitude Range of Grid

    • 本文研究区域范围为欧洲5°W~45°E、30°N~60°N。选取区域中的16个IGS跟踪站(gope、ganp、polv、mat1、wtzr、riga、onsa、morp、ebre、ankr、opmt、pots、leij、bzrg、glsv、mdvj)如图 2所示。观测时间为2017-12-23—2017-12-29,年积日(day of year,DOY)为第357~363天,共7 d连续24 h观测数据,采样间隔为30 s,卫星截止高度角为15°。以每个测站GNSS观测数据为背景,分别建立单站区域的电离层TEC模型,并生成相应的单站电离层TEC格网IONEX格式文件。

      图  2  欧洲区域16个GNSS跟踪站分布图

      Figure 2.  Map of 16 GNSS Tracking Stations in the European Region

    • 电离层TEC的变化受到太阳活动剧烈程度和磁暴的影响较大,相比于空间活动平静时期,太阳剧烈活动和发生磁暴的概率较小且影响时间较短,因此,充分考虑这二者可能带来的不确定影响,本文仅选取太阳活动平静和未发生磁暴的时段进行建模和数据分析。

      观测时段内2017年DOY第357~363天的太阳活动F10.7指数和地磁Kp指数分别如图 3图 4所示。

      图  3  2017年DOY第357~363天太阳F10.7指数

      Figure 3.  Solar F10.7 Index from DOY 357 to 363 in 2017

      图  4  2017年DOY第357~363天地磁Kp指数

      Figure 4.  Geomagnetic Kp Index from DOY 357 to 363 in 2017

      图 3图 4可知,太阳活动F10.7指数均小于100×10-22 W/(m2·Hz),说明太阳活动处于平静时期;地磁Kp指数均小于5,说明未发生磁暴,因此,在此期间解算的电离层TEC具有一般性。

    • 1)分别基于2~15阶次球谐函数建立16个测站的单站电离层TEC模型,并生成电离层格网IONEX格式产品。分别提取每天每个测站对应的电离层格网产品中距离每个测站最近的格网点(即测站0°×0°范围)的TEC值和RMS值,并求7 d的平均值,与CODE、IGS和iGMAS相对应的值进行比较。16个测站2~15阶次球谐函数模型电离层、CODE、IGS和iGMAS在每个测站0°×0°范围7 d的TEC平均值和RMS平均值如图 5所示。图 5中,codg、iscg和igsg分别表示CODE、iGMAS和IGS电离层产品的结果。

      图  5  16个测站0°×0°范围7天的TEC平均值

      Figure 5.  TEC Mean Values of 16 Stations in the Range of 0°×0° for 7 Days

      图 5(a)可知,基于2~15阶次球谐函数建立的16个单站区域电离层TEC模型在每个测站0°×0°范围7 d的TEC平均值约为4~5 TECU,与CODE和iGMAS相当,与IGS的差值小于1 TECU,说明2~15阶次球谐函数建立的单站区域电离层TEC模型在适用性上具有一致性,即在单站单点0°×0°范围内,区域电离层TEC模型的建立上差别不大,但随着球谐函数阶次数的不断增加,其需要估计的系数也不断增加,2×2阶次和15×15阶次球谐函数需要估计的球谐系数分别为9个和256个,后者的估计系数是前者的28倍,在解算效率上,前者明显优于后者,相比之下,可以选择低阶(2×2或3×3)球谐函数进行单站区域电离层TEC建模。

      图 5(b)可知,在RMS方面,基于2~15阶次球谐函数建立的16个单站区域电离层TEC模型在每个测站0°×0°范围7 d的RMS平均值相差较小,且与IGS相当,约为0.9 TECU,与CODE和iGMAS最大相差0.8 TECU。图 5(b)中显示的所有电离层产品的RMS均小于1.0 TECU,均优于CODE、iGMAS或者IGS电离层产品的官方精度(2~8 TECU),说明在这一区域所有电离层产品均具有较高的精度,相比于高阶次球谐函数,可以选择低阶(2×2或3×3)球谐函数进行单站区域电离层TEC建模。

      本文展示了bzrg和wtzr两个单站解算的电离层产品在测站0°×0°范围7 d的TEC平均值和RMS平均值,分别如图 6图 7所示。

      图  6  bzrg站0°×0°范围7天的TEC平均值和RMS平均值

      Figure 6.  Mean Values of TEC and RMS for 7 Days in the Range of 0°×0° at bzrg Station

      图  7  wtzr站0°×0°范围7天的TEC平均值和RMS平均值

      Figure 7.  Mean Values of TEC and RMS for 7 Days in the Range of 0°×0° at wtzr Station

      图 6图 7可知,基于2~15不同阶次球谐函数模型建立的bzrg和wtzr两个单站区域电离层TEC模型,在0°×0°范围内TEC值与CODE、iGMAS和IGS相当,RMS均小于1.0 TECU,与16个测站平均值的结果一致。

    • 以基于2~15不同阶次球谐函数分别建立的16个单站的区域电离层TEC格网文件为基础,取以每个观测站为中心,格网经纬度分别为0°×‍0°、10°×5°、10°×10°、10°×15°、20°×15°和20°×20°范围内所有格网点上电离层TEC的平均值和RMS的平均值,并分别与CODE、iGMAS和IGS的电离层GIM产品中相应的TEC的平均值和RMS的平均值进行对比分析,16个测站不同范围7天的TEC平均值和RMS平均值如图 8所示。

      图  8  16个测站不同范围7天的TEC平均值

      Figure 8.  TEC Mean Values of 16 Stations in Different Ranges for 7 Days

      图 8(a)整体反映了16个单站建立的区域电离层TEC模型在不同区域的TEC均值情况,其中,2~15不同阶次的球谐函数电离层TEC模型的TEC值整体处于CODE、iGMAS和IGS之间,表明若以CODE、iGMAS和IGS的电离层产品为比较基准,则在以单个测站为中心,经纬度范围20°×20°的区域内,本文的方法建立的单站区域电离层模型TEC值与CODE、iGMAS和IGS的电离层产品的TEC值精度相当。

      图 8(b)整体反映了16个单站建立的区域电离层TEC模型在不同区域的RMS均值情况。由图 8(b)可知,以单测站为中心,同一取值范围区域(以20°×20°为例)内,随球谐函数模型阶次的增加,RMS随之增大,且与IGS的RMS相差在3 TECU以内,2×2或3×3阶次球谐函数的RMS则与IGS的RMS相当。以单测站为中心,在10°×10°范围内,RMS约为0.8~1.5 TECU,且与IGS的差值约为0.5 TECU。超过这一个范围,如10°×15°、20°×15°和20°×20°以内,其RMS值与IGS的RMS相差较大,约为2~3 TECU。

      综合图 8(a)图 8(b)可知,就其TEC而言,单站建立的区域电离层TEC模型经纬度范围20°×20°的区域内TEC值与IGS的TEC值相当;但就其RMS而言,经纬度范围10°×10°以外的RMS与IGS的RMS相差较大,经纬度范围10°×‍10°以内的球谐函数阶次越低,其模型RMS越小,因此,以单测站为中心,在10°×10°范围区域内,基于低阶(2×2或3×3)球谐函数建立的单站区域电离层TEC模型的精度与IGS的TEC模型的精度相当,即单站区域电离层TEC模型的适用半径约为600 km,精度约为1 TECU。

      pots站、ganp站解算的区域电离层TEC产品以测站为中心,格网经纬度分别为0°×0°、10°×5°、10°×10°、10°×15°、20°×15°和20°×20°不同范围内所有格网点7天的TEC平均值和RMS平均值,分别如图 9图 10所示。

      图  9  pots站不同范围7天的TEC和RMS平均值

      Figure 9.  Mean Values of TEC and RMS of Different Ranges for 7 Days at pots Station

      图  10  ganp测站不同范围7天的TEC和RMS平均值

      Figure 10.  Mean Values of TEC and RMS of Different Ranges for 7 Days at ganp Station

      图 9可知,基于2~15阶次球谐函数,以pots测站建立的单站区域电离层TEC模型,在距离测站中心不同区域范围内,其TEC值比CODE和iGMAS的TEC值高约0.6 TECU,比IGS的TEC值低约0.2 TECU,整体位于CODE、iGMAS和IGS电离层产品之间。RMS随球谐函数阶次和取值区域范围的增大而增大,在10°×10°范围以内,单站模型的RMS与IGS的RMS相当。

      图 10可知,基于2~15阶次球谐函数,以ganp测站建立的单站区域电离层TEC模型,在距离测站中心不同区域范围内,其TEC值与CODE和iGMAS的TEC值相当,比IGS小约0.4 TECU,整体位于CODE、iGMAS和IGS电离层产品之间。RMS随球谐函数阶次和取值区域范围的增大而增大,在10°×10°范围以内,单站模型的RMS与IGS的RMS相当。

      pots站和ganp站所观测的GPS卫星星下点轨迹如图 11所示。由图 11可知,pots和ganp两个测站的GPS卫星星下点轨迹范围均为经度约40°,纬度约25°,即经纬度范围约40°×25°,其范围远大于以上的数据分析的最大范围20°×20°,在进一步的定量数据分析和精度比较评估中,将单站区域电离层TEC的适用范围缩小至经纬度10°×‍10°(半径小于600‍ km),在精度和范围上均足以满足对某一城市单频用户电离层延迟改正的需求,既提高了电离层产品的实效性,精度也与CODE、iGMAS和IGS电离层产品相当,约为1 TECU,因此,对区域小范围内(半径小于600 km)为单频用户提供高精度的电离层延迟改正具有重要意义,同时,在单频导航定位综合服务中具有广阔的应用前景。

      图  11  pots站观测的GPS卫星星下点轨迹图

      Figure 11.  GPS Sub Satellite Point Trajectory Map Observed at pots Station

      本文展示了部分单站电离层TEC云图,2017-12-25(DOY第359天)协调世界时(coordinate universal time,UTC)10:00、12:00、14:00,IGS和基于2×2和3×3阶次球谐函数以ganp单站观测数据建立电离层TEC模型的电离层TEC云图如图 12所示。

      图  12  IGS和ganp单站2~3阶次区域电离层TEC云图

      Figure 12.  TEC Cloud Images of Second-and-Third-Order Ionosphere of IGS and ganp Stations

      2017-12-25(DOY第359天)UTC 10:00、12:‍00、14:00,IGS和基于2×2和3×3阶次球谐函数以pots单站观测数据建立电离层TEC模型的电离层TEC云图如图 13所示。

      图  13  IGS和pots单站2~3阶次区域电离层TEC云图

      Figure 13.  TEC Cloud Images of Second-and-Third-Order Ionosphere of IGS and pots Stations

      图 12图 13可知,以ganp或pots测站为中心,格网经纬度在10°×10°范围以内,以ganp或pots建立的单站区域电离层TEC模型与IGS全球电离层TEC模型具有高度一致性,表明基于低阶(2×2或3×3)球谐函数建立的单站区域电离层TEC模型的精度与IGS全球电离层TEC模型精度相当。

    • 为了分析单站区域电离层TEC建模的适用范围和精度,本文基于2~15阶次球谐函数,分别建立了欧洲区域16个单站区域电离层TEC模型,并生成了区域电离层TEC格网。通过与相同格式的CODE、iGMAS和IGS等权威电离层TEC产品的比较分析,结果表明,单站建立的区域电离层TEC模型经纬度范围20°×20°的区域内TEC值与IGS的TEC值相当,相应的RMS随球谐函数阶次和取值区域范围的增大而增大。顾及TEC值和RMS,整体而言,在以单站为中心、经纬度范围10°×10°以内(半径小于600 km),基于低阶(2×2或3×3)球谐函数建立的单站区域电离层TEC模型的精度与CODE、iGMAS和IGS等的全球电离层TEC模型相当,其精度约为1.0 TECU。因此,可以基于低阶(2×2或3×3)球谐函数建立单站区域电离层TEC模型,在半径600 km范围内代替全球电离层TEC模型,极大地提高了电离层TEC建模的解算效率和产品实效性,对区域单频用户提供及时、高精度的电离层延迟改正信息具有重要意义。

参考文献 (20)

目录

    /

    返回文章
    返回