盛傳貞,張京奎,張寶成
( 1.衛(wèi)星導(dǎo)航系統(tǒng)與裝備技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,石家莊 050081;2.中國科學(xué)院精密測量科學(xué)與技術(shù)創(chuàng)新研究院,武漢 430000 )
基于全球衛(wèi)星導(dǎo)航系統(tǒng)(GNSS)建立的全球電離層格網(wǎng)產(chǎn)品(GIM)是全球電離層總電子含量(TEC)研究和應(yīng)用的重要數(shù)據(jù)資源,不僅是反演電離層變化規(guī)律、揭示電離層物理演變規(guī)律的基礎(chǔ)數(shù)據(jù)產(chǎn)品[1-2],而且也可作為單頻導(dǎo)航定位研究中電離層延遲修正的重要技術(shù)手段[3-4].國際GNSS服務(wù)(IGS)的電離層綜合中心(IAAC)和全球連續(xù)監(jiān)測評估系統(tǒng)(iGMAS)綜合中心發(fā)布了預(yù)報(bào)GIM、快速GIM 和最終GIM 三種類型的電離層產(chǎn)品[5].由于不同機(jī)構(gòu)解算不同類型GIM 產(chǎn)品時(shí)采用了不同的數(shù)學(xué)模型及GNSS跟蹤站數(shù)據(jù),解算策略的不同導(dǎo)致不同機(jī)構(gòu)產(chǎn)生的不同類型GIM 產(chǎn)品精度存在差異[6].
系統(tǒng)深入評估和分析了不同GIM產(chǎn)品在不同區(qū)域的應(yīng)用效果,可為不同區(qū)域的電離層TEC探測、環(huán)境和災(zāi)害監(jiān)測及導(dǎo)航用戶電離層延遲修正服務(wù)提供參考[7].國內(nèi)外學(xué)者針對不同GIM產(chǎn)品的應(yīng)用性能開展了大量的精度評估工作,然而目前已開展的研究絕大部分都是針對不同分析中心的單一類型GIM 產(chǎn)品在全球范圍內(nèi)的平均應(yīng)用性能[8-12].由于電離層活動具有顯著的區(qū)域特性,特別是我國低緯地區(qū)受到電離層赤道異常影響[13],不同GIM產(chǎn)品在全球范圍內(nèi)的評估結(jié)果不能直接適用于中國區(qū)域.文獻(xiàn)[12]和文獻(xiàn)[14]討論了不同類型的GIM 產(chǎn)品在全球范圍的內(nèi)符合精度,文獻(xiàn)[8]討論了四種最終GIM 產(chǎn)品在中國區(qū)域的精度.然而,目前就不同機(jī)構(gòu)不同類型的GIM產(chǎn)品在中國區(qū)域的精度分析工作還沒有開展.本文基于不同太陽活動水平下GNSS跟蹤站和Jason 測高衛(wèi)星實(shí)測電離層TEC對不同機(jī)構(gòu)不同類型的GIM產(chǎn)品在中國及周邊海域區(qū)域的精度進(jìn)行評估,并分析了不同GIM產(chǎn)品對中國區(qū)域單頻導(dǎo)航定位精度的影響.
目前,發(fā)布GIM 產(chǎn)品的IGS的電離層聯(lián)合分析中心(IAACs)主要包括美國噴氣推進(jìn)實(shí)驗(yàn)室(JPL)、歐洲定軌中心(CODE)、歐洲太空局(ESA/ESOC)、西班牙加泰羅尼亞理工大學(xué)(UPC)、中國科學(xué)院(CAS)和武漢大學(xué)衛(wèi)星導(dǎo)航定位技術(shù)研究中心(WHU)等機(jī)構(gòu).我國建立的iGMAS綜合中心也發(fā)布了綜合國內(nèi)各分析中心產(chǎn)品后的快速GIM 及最終GIM 產(chǎn)品.根據(jù)GIM 產(chǎn)品時(shí)效性,GIM 產(chǎn)品主要分為快速GIM(延遲1天)、最終GIM(延遲11天)以及預(yù)報(bào)GIM(提前1天或2天)產(chǎn)品.
如表1所示,本文評估的GIM 產(chǎn)品包含iGMAS綜合中心提供的快速GIM 和最終GIM 以及IAACs的分析中心CODE、UPC、ESA、CAS、JPL提供的預(yù)報(bào)GIM、快速GIM 和最終GIM 產(chǎn)品.
表1 不同機(jī)構(gòu)發(fā)布的GIM 產(chǎn)品
本文將GNSS參考站(后文統(tǒng)稱為“跟蹤站”)實(shí)測數(shù)據(jù)獲取的電離層TEC變化量和測高衛(wèi)星實(shí)測電離層TEC為參考,分析不同機(jī)構(gòu)不同類型GIM 產(chǎn)品在中國及周邊區(qū)域的電離層延遲改正精度及其對中國區(qū)域?qū)Ш接脩舳ㄎ痪鹊挠绊?
1)與GNSS實(shí)測數(shù)據(jù)獲取的電離層傾斜總電子含量(STEC)歷元間變化進(jìn)行對比
在未發(fā)生周跳的連續(xù)觀測弧段內(nèi),基于GNSS跟蹤站的無幾何組合的載波相位觀測值歷元間變化可以計(jì)算衛(wèi)星和接收機(jī)視線方向上電離層TEC的歷元間變化,其精度優(yōu)于0.1 TECU.由于衛(wèi)星高度角最大時(shí)由天頂方向總電子含量(VTEC)轉(zhuǎn)換為STEC時(shí)受到的投影函數(shù)誤差最小,實(shí)際計(jì)算電離層歷元間變化量時(shí),通常選擇衛(wèi)星高度角最大時(shí)刻的載波相位觀測值為參考,具體計(jì)算公式[15]如下:
由式(1)可知,采用電離層TEC歷元間變化方法進(jìn)行GIM 精度驗(yàn)證時(shí),無需額外計(jì)算衛(wèi)星和接收機(jī)的差分碼偏差(DCB)參數(shù).因此,與直接利用扣除DCB影響的GNSS電離層TEC測量值進(jìn)行精度驗(yàn)證相比,電離層TEC歷元間變化方法的優(yōu)勢在于消除了衛(wèi)星和接收機(jī)DCB解算誤差對于驗(yàn)證結(jié)果的不利影響.電離層TEC歷元間變化方法是IGS電離層工作組評估各電離層分析中心GIM 產(chǎn)品精度采用最廣泛的方法之一.
2)與測高衛(wèi)星電離層TEC對比
測高衛(wèi)星通過向海洋表面發(fā)射的雙頻信號可實(shí)現(xiàn)信號傳播路徑上的電離層TEC反演.測高衛(wèi)星電離層TEC是目前獨(dú)立于GNSS之外覆蓋性較廣、連續(xù)觀測時(shí)間最長的電離層觀測數(shù)據(jù),也是IGS電離層工作組用來評估不同分析中心GIM 產(chǎn)品性能的有效數(shù)據(jù)之一.需要說明的是:由于測高衛(wèi)星僅能測量從海洋表面測高衛(wèi)星軌道高度以下的電離層電子密度,其測量的電離層TEC與GIM 電離層TEC通常存在量級為2~5 TECU 的系統(tǒng)性偏差.
3)GIM 產(chǎn)品對定位精度影響評估
對于單頻衛(wèi)星導(dǎo)航用戶,可以采用GIM 產(chǎn)品來削弱電離層延遲對導(dǎo)航定位結(jié)果的影響.為評估和掌握不同GIM產(chǎn)品在定位中的應(yīng)用效果,利用GNSS跟蹤站觀測數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)單點(diǎn)定位,評估不同GIM 產(chǎn)品對定位的影響.以IGS發(fā)布的精密周解坐標(biāo)產(chǎn)品作為坐標(biāo)參考“真值”,基于動態(tài)標(biāo)準(zhǔn)單點(diǎn)定位估計(jì)測站坐標(biāo)及接收機(jī)鐘差參數(shù).標(biāo)準(zhǔn)單點(diǎn)定位數(shù)據(jù)處理中,采用的數(shù)據(jù)為GPS和北斗衛(wèi)星導(dǎo)航系統(tǒng)(BDS)觀測數(shù)據(jù)采樣率為30 s,衛(wèi)星截止高度角為15°,衛(wèi)星軌道及鐘差采用廣播星歷的軌道及鐘差參數(shù),對流層誤差采用UNB3模型改正.
本文采用不同電離層模型的TEC計(jì)算值與電離層TEC參考值之間的平均偏差(BIAS)、標(biāo)準(zhǔn)差(STD)、均方根(RMS)和相對誤差(PER)作為精度分析的指標(biāo),其計(jì)算公式為
為了合理全面的評估不同機(jī)構(gòu)不同類型GIM產(chǎn)品在中國區(qū)域的實(shí)際應(yīng)用效果,選取不同太陽活動水平條件下的數(shù)據(jù)開展性能評估工作.由于iGMAS從2015年開始才穩(wěn)定提供綜合快速GIM 和最終GIM 產(chǎn)品,本文選取的數(shù)據(jù)測試時(shí)段為最近一個(gè)太陽活動周期中2015年(太陽活動中年)和2020年(太陽活動低年)各一個(gè)月的數(shù)據(jù)(年積日第100—130天).為評估不同GIM 產(chǎn)品的外符合精度,本文選擇中國地殼運(yùn)動觀測網(wǎng)絡(luò)(CMONC)跟蹤站以及Jason 測高衛(wèi)星獲取的電離層TEC數(shù)據(jù)進(jìn)行GIM產(chǎn)品精度驗(yàn)證.為充分利用IGS發(fā)布的精密坐標(biāo)產(chǎn)品,在研究GIM 產(chǎn)品對定位精度影響時(shí)選擇中國區(qū)域的IGS跟蹤站進(jìn)行實(shí)驗(yàn).如圖1所示,本文選取的GNSS數(shù)據(jù)由59 個(gè)CMONC跟蹤站和8個(gè)IGS跟蹤站組成.
圖1 選取的GNSS檢測站及Jason 穿刺點(diǎn)軌跡分布圖
需要說明的是:由于獲取不到2020年的Jason測高衛(wèi)星數(shù)據(jù),選擇2019年度的Jason-3測高衛(wèi)星數(shù)據(jù)進(jìn)行實(shí)驗(yàn).此外,部分機(jī)構(gòu)的GIM 產(chǎn)品存在某些時(shí)段缺失的現(xiàn)象,如iGMAS分別缺失2015年年積日第101天、102天和124天的快速GIM 產(chǎn)品;ESA 分別缺失2020年的預(yù)報(bào)產(chǎn)品、2015年年積日第101—103和120—122天的快速GIM產(chǎn)品;UPC分別缺失2019年的快速GIM 產(chǎn)品和2020年度的預(yù)報(bào)GIM 產(chǎn)品.
圖2、圖3分別給出了不同機(jī)構(gòu)GIM 產(chǎn)品電離層TEC的歷元間變化RMS精度時(shí)間序列及其在各CMONC跟蹤站的分布情況,不同太陽活動水平下的精度分布繪制在虛線的左右兩側(cè).表2~表4分別給出了電離層TEC的歷元間變化平均精度統(tǒng)計(jì)情況.從上述圖表可知,不同機(jī)構(gòu)的GIM 產(chǎn)品在中國區(qū)域的電離層TEC的歷元間變化精度具有如下特點(diǎn):
表2 不同機(jī)構(gòu)預(yù)測GIM 產(chǎn)品電離層TEC的歷元間變化精度統(tǒng)計(jì)TECU
表3 不同機(jī)構(gòu)快速GIM 產(chǎn)品電離層TEC的歷元間變化精度統(tǒng)計(jì)TECU
表4 不同機(jī)構(gòu)最終GIM 產(chǎn)品電離層TEC的歷元間變化精度統(tǒng)計(jì)TECU
圖2 2015年(虛線左側(cè))和2020年(虛線右側(cè))不同機(jī)構(gòu)GIM 產(chǎn)品的電離層TEC的歷元間變化精度時(shí)間序列
圖3 2015年(虛線左側(cè))和2020年(虛線右側(cè))不同機(jī)構(gòu)GIM產(chǎn)品的電離層TEC的歷元間變化精度在各跟蹤站的分布情況
1)各機(jī)構(gòu)不同類型GIM產(chǎn)品的精度由高到低分別是最終GIM、快速GIM 和預(yù)報(bào)GIM 產(chǎn)品;各GIM產(chǎn)品在太陽活動中年的精度通常低于太陽活動低年,且在低緯度跟蹤站的精度通常低于中緯度跟蹤站,這與電離層在不同太陽活動條件及不同地理緯度的變化復(fù)雜程度密切相關(guān).
2)在太陽活動中年,iGMAS綜合的快速和最終GIM 在同類產(chǎn)品中的精度最高,其RMS均值為2.65 TECU;除iGMAS綜合的GIM 產(chǎn)品外,UPC的最終GIM 和快速GIM 產(chǎn)品精度優(yōu)于其他機(jī)構(gòu)的同類產(chǎn)品,其RMS值分別為2.69 TECU 和2.70 TECU;CAS的快速GIM 和ESA 的最終GIM 產(chǎn)品精度略低于其他機(jī)構(gòu)同類產(chǎn)品.
3)在太陽活動低年,JPL 的快速GIM 和最終GIM 在同類產(chǎn)品中精度最低;盡管iGMAS的快速GIM 產(chǎn)品精度在某些天存在異?,F(xiàn)象,與其他機(jī)構(gòu)的快速GIM 產(chǎn)品精度差異仍小于0.11 TECU;CODE和iGMAS的最終產(chǎn)品精度略優(yōu)于其他機(jī)構(gòu)的最終GIM 產(chǎn)品.
4)CODE 和UPC的預(yù)報(bào)GIM 產(chǎn)品精度略優(yōu)于ESA 的預(yù)報(bào)產(chǎn)品.CODE 預(yù)報(bào)1天與預(yù)報(bào)2天的GIM產(chǎn)品精度大致相當(dāng).
圖4給出了Jason-2和Jason-3測高衛(wèi)星電離層TEC和GIM 電離層TEC之差的均值和標(biāo)準(zhǔn)差,表5給出了快速和最終GIM電離層TEC的RMS精度統(tǒng)計(jì)情況.其中,Jason-2和Jason-3測高衛(wèi)星有效觀測數(shù)據(jù)時(shí)間分別為2015年年積日第100—130天和2019年年積日第102—130天.從圖4可以看出,以Jason 測高衛(wèi)星獲取的TEC為參考,不同機(jī)構(gòu)的GIM產(chǎn)品在中國周邊海洋區(qū)域的精度有以下特點(diǎn):
圖4 不同GIM產(chǎn)品相對于Jason VTEC的BIAS和STD
表5 不同機(jī)構(gòu)快速和最終GIM 產(chǎn)品相對于Jason VTEC的RMS精度統(tǒng)計(jì)TECU
1) Jason 測高衛(wèi)星電離層TEC與不同機(jī)構(gòu)各GIM的電離層TEC的系統(tǒng)性偏差均小于0,這與測高衛(wèi)星觀測數(shù)據(jù)不包含測高衛(wèi)星軌道高度至GNSS衛(wèi)星軌道高度處的電離層TEC是一致的.
2)測高衛(wèi)星和不同機(jī)構(gòu)GIM 電離層TEC的差異在電離層活動中年(2015年)時(shí)的差異較低年(2019年)要明顯偏大,這與電離層在活動水平存在一定的相關(guān)性.不同機(jī)構(gòu)的預(yù)報(bào)GIM、快速GIM 和最終GIM與測高衛(wèi)星TEC差異的標(biāo)準(zhǔn)差在太陽活動中年與低年BIAS分別為3.56 TECU、4.36 TECU 和4.63 TECU.
3)對比太陽活動中年不同參考基準(zhǔn)獲取的GIM標(biāo)準(zhǔn)差結(jié)果發(fā)現(xiàn),以實(shí)測GNSS電離層TEC為基準(zhǔn)獲取的GIM精度約為測高衛(wèi)星電離層TEC為基準(zhǔn)的精度的2倍,表明不同機(jī)構(gòu)的GIM 產(chǎn)品在中國周邊海洋區(qū)域的精度顯著低于GNSS數(shù)據(jù)覆蓋的中國陸地區(qū)域.
4)從系統(tǒng)偏差指標(biāo)來看,JPL 的快速和最終GIM產(chǎn)品在同類型產(chǎn)品中與Jason 測高衛(wèi)星電離層TEC的系統(tǒng)偏差最大,平均差值均為5.5 TECU,其他機(jī)構(gòu)的快速和最終GIM產(chǎn)品與Jason 測高衛(wèi)星TEC的BIAS分別在1.88~2.84 TECU 和2.68~4.07 TECU 之間.
5)從標(biāo)準(zhǔn)差指標(biāo)來看,除JPL 和UPC的快速和最終GIM產(chǎn)品精度大致相當(dāng)外,其他機(jī)構(gòu)均表現(xiàn)為同一機(jī)構(gòu)最終GIM 精度略高于快速和預(yù)報(bào)GIM產(chǎn)品;CODE 和UPC的預(yù)報(bào)GIM 產(chǎn)品精度大致相當(dāng),ESA 的預(yù)報(bào)產(chǎn)品精度在三種預(yù)報(bào)產(chǎn)品中精度最差.
6)從RMS指標(biāo)來看,在太陽活動中年,UPC的快速GIM 和最終GIM 產(chǎn)品在同類產(chǎn)品中精度最高;在太陽活動低年,iGMAS的最終GIM和CAS的快速GIM 在同類產(chǎn)品中精度最高;JPL 的快速GIM和最終GIM 產(chǎn)品在不同太陽活動水平下均顯著低于其他結(jié)構(gòu)的同類型GIM 產(chǎn)品.
以2015年和2020年年積日第100 —130天選取的中國區(qū)域IGS測站,對用戶采用不同GIM模型的單頻標(biāo)準(zhǔn)單點(diǎn)定位定位結(jié)果作進(jìn)一步分析.圖5和圖6分別給出了所有跟蹤站在天頂和水平方向上的平均定位誤差的1天內(nèi)時(shí)間序列圖.表6至表8分別給出了試驗(yàn)期間不同GIM 產(chǎn)品在所有測試站的天頂和水平方向定位誤差的RMS統(tǒng)計(jì)情況.對比上述圖表可知,不同機(jī)構(gòu)GIM產(chǎn)品獲取的定位精度存在如下特點(diǎn):
表6 不同快速GIM 產(chǎn)品在各測試站天頂和水平方向上定位誤差RMS統(tǒng)計(jì)m
表7 不同最終GIM 產(chǎn)品在各跟蹤站天頂和水平方向上定位誤差RMS統(tǒng)計(jì)m
表8 不同預(yù)報(bào)GIM 產(chǎn)品和不施加電離層修正在各跟蹤站定位誤差RMS統(tǒng)計(jì)m
圖5 2015年(虛線左側(cè))和2020年(虛線右側(cè))不同GIM產(chǎn)品獲取的水平方向定位精度時(shí)間序列
圖6 2015年(虛線左側(cè))和2020年(虛線右側(cè))不同GIM產(chǎn)品獲取的高程方向定位精度時(shí)間序列
1)采用不同GIM產(chǎn)品進(jìn)行標(biāo)準(zhǔn)單點(diǎn)定位定位與不采用電離層延遲修正相比在天頂和水平方向的定位精度均有所提升,且不同GIM產(chǎn)品獲取的定位精度差異主要體現(xiàn)在天頂方向上.與未對電離層延遲進(jìn)行修正的單點(diǎn)定位結(jié)果相比,不同機(jī)構(gòu)的預(yù)報(bào)GIM、快速GIM和最終GIM 產(chǎn)品在太陽活動中年在天頂方向上的定位精度分別提升65.8%~73.9%、72.7%~74%和62.5%~71.8%,在水平方向上的定位精度分別提升28.5%~44.8%、36.1%~46.2%和34.7%~41%.
2)不同機(jī)構(gòu)的GIM產(chǎn)品在太陽活動低年獲取的定位精度差異較小,在太陽活動中年差異較為顯著,且最終GIM 產(chǎn)品獲取的定位精度差異較快速GIM產(chǎn)品小.在太陽活動中年,UPC的最終GIM 和快速GIM 在同類GIM產(chǎn)品中獲取的水平方向定位精度最高;除ESA 產(chǎn)品外,不同機(jī)構(gòu)的快速GIM、最終GIM 和預(yù)報(bào)GIM 產(chǎn)品獲取的天頂方向定位誤差最大差異分別為0.17 m、0.06 m 和0.06 m,獲取的水平方向定位誤差最大差異分別為0.4 m、0.11 m 和0.18 m.在太陽活動低年,不同機(jī)構(gòu)的預(yù)報(bào)GIM、快速GIM產(chǎn)品獲取的天頂方向定位誤差最大差異分別為0.05 m、0.06 m 和0.01 m,獲取的水平方向定位誤差最大差異分別為0.04 m、0.05 m 和0 m.
3)CODE的兩種預(yù)報(bào)GIM 和UPC的預(yù)報(bào)GIM產(chǎn)品精度大致相當(dāng),二者獲取的天頂方向上的定位精度要顯著優(yōu)于ESA 的預(yù)報(bào)GIM 產(chǎn)品.
4)根據(jù)太陽活動中年各跟蹤站上獲取的定位時(shí)間天內(nèi)時(shí)間序列可知,白天的定位精度比晚上定位精度低,這與電離層日間變化趨勢是一致的.
本文基于中國區(qū)域均勻分布的59個(gè)陸態(tài)網(wǎng)跟蹤站數(shù)據(jù)以及8個(gè)IGS跟蹤站對不同機(jī)構(gòu)發(fā)布的預(yù)報(bào)GIM、快速GIM 和最終GIM 產(chǎn)品在中國及周邊海域區(qū)域的應(yīng)用精度進(jìn)行了評估.通過對比GIM、跟蹤站實(shí)測GNSS數(shù)據(jù)和Jason 測高衛(wèi)星的電離層TEC,并分析不同GIM 產(chǎn)品對中國區(qū)域?qū)Ш接脩魡晤l定位精度的影響,針對不同GIM 產(chǎn)品在中國及周邊海域區(qū)域的應(yīng)用精度,得到以下結(jié)論:
1)不同類型GIM 產(chǎn)品在中國區(qū)域的精度由高到低分別為最終GIM、快速GIM 和預(yù)報(bào)GIM 產(chǎn)品;
2)不同機(jī)構(gòu)不同類型GIM 產(chǎn)品的精度差異在太陽活動中年大于太陽活動低年,在低緯度區(qū)域高于中緯度區(qū)域;
3)當(dāng)采用不同機(jī)構(gòu)不同類型的GIM 產(chǎn)品進(jìn)行中國區(qū)域定位試驗(yàn)時(shí),除ESA 產(chǎn)品獲取的定位精度較差外,其他機(jī)構(gòu)的GIM產(chǎn)品獲取的定位精度大致相當(dāng);
4)總體而言,iGMAS和UPC的快速和最終GIM產(chǎn)品精度優(yōu)于其他機(jī)構(gòu)的GIM 產(chǎn)品,基于GIM 數(shù)據(jù)進(jìn)行中國區(qū)域電離層探測時(shí),建議用戶使用時(shí)優(yōu)先選擇iGMAS和UPC的快速和最終GIM產(chǎn)品;
5) CODE 預(yù)報(bào)1天的GIM產(chǎn)品與預(yù)報(bào)2天的GIM產(chǎn)品精度相當(dāng),優(yōu)于ESA 的預(yù)報(bào)GIM 產(chǎn)品精度.
隨著不同分析中心GIM 解算策略的不斷改進(jìn),不同類型GIM產(chǎn)品的精度和穩(wěn)定性將會進(jìn)一步提升,未來需要采用更長時(shí)段的數(shù)據(jù)對不同機(jī)構(gòu)不同類型的GIM 產(chǎn)品在中國區(qū)域的應(yīng)用精度進(jìn)行進(jìn)一步分析.此外,由于實(shí)時(shí)GIM 產(chǎn)品數(shù)據(jù)較少,本文未對實(shí)時(shí)GIM產(chǎn)品的精度進(jìn)行評估,未來也需進(jìn)一步分析實(shí)時(shí)GIM 產(chǎn)品在中國區(qū)域的應(yīng)用精度.