磁赤道地區(qū)2007
—2013年COSMIC掩星反演和國(guó)際參考電離層模型輸出結(jié)果分析

黃智， 袁洪

1 江蘇師范大學(xué)物理與電子工程學(xué)院， 江蘇 徐州　221116　2 中國(guó)科學(xué)院光電研究院， 北京　100094

?

磁赤道地區(qū)2007
—2013年COSMIC掩星反演和國(guó)際參考電離層模型輸出結(jié)果分析

黃智1， 袁洪2

1 江蘇師范大學(xué)物理與電子工程學(xué)院， 江蘇 徐州2211162 中國(guó)科學(xué)院光電研究院， 北京100094

摘要利用2007—2013年電離層測(cè)高儀位于磁赤道觀測(cè)站Jicamarca (11.95°S，76.8°W，地磁緯度為1°N)的垂測(cè)數(shù)據(jù)和COSMIC掩星反演電離層資料，分析了不同太陽(yáng)活動(dòng)條件下兩種探測(cè)技術(shù)獲取電離層特征參數(shù)峰值密度NmF2和峰值高度HmF2的相關(guān)性，同時(shí)也探討了國(guó)際參考電離層模型IRI輸出參數(shù)與測(cè)高儀垂測(cè)數(shù)據(jù)的相關(guān)性.此外，進(jìn)一步分析了COSMIC掩星和IRI模型在不同地方時(shí)高估或低估垂測(cè)參數(shù)的分布特征.結(jié)果表明：(1)由COSMIC掩星反演和IRI模型輸出參數(shù)NmF2與測(cè)高儀垂測(cè)值NmF2得到的相關(guān)系數(shù)都在0.8以上.太陽(yáng)活動(dòng)低年COSMIC探測(cè)得到的NmF2相關(guān)性高于太陽(yáng)活動(dòng)高年得到的結(jié)果，但I(xiàn)RI模型在太陽(yáng)高年得到的NmF2相關(guān)性好于太陽(yáng)活動(dòng)低年的計(jì)算結(jié)果.(2) 由COSMIC掩星反演和IRI模型輸出的HmF2與測(cè)高儀垂測(cè)值HmF2在春秋季得到的相關(guān)性較高，夏季的相關(guān)性最弱.由COSMIC掩星探測(cè)HmF2得到的季節(jié)和時(shí)間相關(guān)系數(shù)大都集中在0.8和0.6以上，但由IRI模型得到的HmF2相關(guān)性降低，特別是太陽(yáng)活動(dòng)低年的夏季和傍晚其相關(guān)系數(shù)低于0.3.(3)太陽(yáng)低年COSMIC掩星和IRI模型白天時(shí)段大都高估、夜間至凌晨前后低估電離層參數(shù)NmF2和HmF2；但太陽(yáng)活動(dòng)高年NmF2在地方時(shí)午夜后則呈現(xiàn)高估的特點(diǎn).相關(guān)結(jié)果為未來(lái)IRI模型的進(jìn)一步完善以及低緯地區(qū)電離層同化模式研究提供參考.

關(guān)鍵詞垂測(cè)儀； COSMIC； IRI； 電離層特征參數(shù)； 相關(guān)性

1引言

位于地面60～1000 km的電離層是日地空間環(huán)境中承上啟下的重要環(huán)節(jié)和關(guān)鍵層次，電離層含有的帶電粒子成分顯著改變了大氣的運(yùn)動(dòng)行為，會(huì)對(duì)穿越其中的電波造成很大的影響.其中，電離層F2層峰值密度NmF2和峰值高度HmF2是表征電離層形態(tài)的重要參數(shù)，一直是電離層理論和工程應(yīng)用領(lǐng)域的重點(diǎn)研究對(duì)象.

利用電離層測(cè)高儀能夠測(cè)量給定觀測(cè)點(diǎn)電子密度的高度分布剖面，從而進(jìn)一步獲取電離層峰值密度NmF2和峰值高度HmF2.利用測(cè)高儀反演得到的電離層特征參數(shù)具有較高的精度，是電離層研究中廣泛采用的探測(cè)技術(shù)之一.但由于河流、雪域以及極區(qū)等自然條件的限制很可能導(dǎo)致某些地區(qū)無(wú)法布站，這會(huì)造成局部電離層形態(tài)描述的缺失.國(guó)內(nèi)外學(xué)者通過(guò)構(gòu)建電離層模式能夠有效補(bǔ)償觀測(cè)數(shù)據(jù)缺失的問(wèn)題，其中最具代表性的電離層經(jīng)驗(yàn)?zāi)Ｊ绞怯蓢?guó)際空間委員會(huì)和國(guó)際無(wú)線電聯(lián)盟共同發(fā)起的項(xiàng)目——國(guó)際參考電離層模型(IRI) (Bilitza, 2001; Bilitza et al., 2014).IRI利用測(cè)高儀、非相干散射雷達(dá)、衛(wèi)星觀測(cè)資料以及探空火箭的觀測(cè)數(shù)據(jù)，同時(shí)引入太陽(yáng)活動(dòng)和地磁指數(shù)，并融合多個(gè)大氣參數(shù)模型，給定年份、日期、時(shí)間、地理(磁)緯度和經(jīng)度等，能夠較為準(zhǔn)確地描述全球范圍內(nèi)特定時(shí)間和地點(diǎn)的電離層形態(tài)變化.目前IRI已廣泛應(yīng)用于電離層形態(tài)學(xué)和空間工程電波修正研究，但它畢竟是電離層的經(jīng)驗(yàn)?zāi)Ｊ?，如果缺乏某一特定時(shí)間和空間區(qū)域內(nèi)的觀測(cè)數(shù)據(jù)，模型的性能很可能會(huì)降低.

2006年4月，由美國(guó)和中國(guó)臺(tái)灣地區(qū)聯(lián)合開發(fā)的空間科學(xué)項(xiàng)目氣象、電離層和氣候衛(wèi)星聯(lián)合觀測(cè)系統(tǒng)COSMIC(F3/C)正式投入運(yùn)行.COSMIC發(fā)射6顆衛(wèi)星，分布于距離地面800 km高度、傾角為72°的軌道面上.每顆衛(wèi)星均裝有4根天線，接收GPS衛(wèi)星L1和L2頻段信號(hào).GPS衛(wèi)星信號(hào)穿越電離層和大氣層因電子密度、水汽和溫度等因素影響而發(fā)生折射，產(chǎn)生時(shí)間延遲，從而推算電離層和大氣層中與之相關(guān)的物理參數(shù).由于COSMIC不受海陸空間分布的限制，不僅為全球范圍的電離層觀測(cè)提供了有效的途徑，同時(shí)為電離層模式的完善補(bǔ)充數(shù)據(jù).但COSMIC掩星探測(cè)采用的反演方法是利用傳播路徑上積分總電子含量或彎曲角，基于局部球?qū)ΨQ、直線傳播、頂部梯度一階近似等假設(shè)，采用Abel積分變化得到電離層電子密度高度分布剖面(Lei et al., 2007).反演的電子密度剖面由于地理位置的差異，與真實(shí)分布情況會(huì)有不同程度的偏差.國(guó)內(nèi)外部分學(xué)者針對(duì)電離層特征參量對(duì)COSMIC觀測(cè)數(shù)據(jù)進(jìn)行了評(píng)估.Lei等(2007)利用Millstone Hill非相干散射雷達(dá)觀測(cè)數(shù)據(jù)和SPIDR發(fā)布的垂直探測(cè)數(shù)據(jù)比較分析了COSMIC反演電離層參數(shù)NmF2和HmF2的探測(cè)精度，結(jié)果發(fā)現(xiàn)地基觀測(cè)與天基掩星觀測(cè)具有較好的一致性.Liu 等(2010)利用電離層數(shù)字測(cè)高儀Jicamarca的垂測(cè)數(shù)據(jù)，比較分析了2007年COSMIC掩星反演電離層特征參數(shù)的性能，并解釋了白天高估垂測(cè)值NmF2的原因.Chuo 等(2011)針對(duì)Jicamarca觀測(cè)站太陽(yáng)低年2006年5月至2008年4月期間的垂測(cè)資料與COSMIC掩星反演得到的電離層參數(shù)B0、NmF2和HmF2之間的關(guān)系進(jìn)行了深入細(xì)致的分析與探討.利用2006—2008年期間巴西低緯及赤道地區(qū)三個(gè)觀測(cè)站數(shù)據(jù)，Ely 等(2012)對(duì)電離層垂測(cè)數(shù)據(jù)和COSMIC掩星兩種探測(cè)技術(shù)得到的電離層參數(shù)進(jìn)行了比較，結(jié)果發(fā)現(xiàn)二者之間盡管具有較好的相關(guān)性，但也存在明顯的差異，并分析了引起差異的原因.孫凌峰等(2014)分析了中國(guó)地區(qū)三個(gè)位于不同緯度觀測(cè)站的垂測(cè)數(shù)據(jù)與COSMIC探測(cè)參數(shù)在太陽(yáng)活動(dòng)高年2012的相關(guān)性，結(jié)果表明二者之間的相關(guān)性和緯度分布、季節(jié)變化密切相關(guān).Hu等(2014)利用中國(guó)地區(qū)120°E經(jīng)線上不同緯度的測(cè)高儀觀測(cè)數(shù)據(jù)詳細(xì)分析了COSMIC反演參數(shù)峰值密度和峰值高度的精度，結(jié)果發(fā)現(xiàn)峰值密度的誤差與緯度密切相關(guān)，峰值高度的相對(duì)誤差小于峰值密度.

上述研究取得了一些有意義的結(jié)果，為合理利用天基掩星開展電離層形態(tài)研究提供了有價(jià)值的參考.但上述研究大都集中在COSMIC系統(tǒng)運(yùn)行初期即太陽(yáng)活動(dòng)低年，缺乏長(zhǎng)期觀測(cè)數(shù)據(jù)的統(tǒng)計(jì)分析.COSMIC掩星自2006年4月成功運(yùn)行以來(lái)，迄今為止已積累了大量原始觀測(cè)數(shù)據(jù)，為電離層特征參量在不同太陽(yáng)活動(dòng)條件下的性能評(píng)估提供了豐富的數(shù)據(jù)資源.低緯地區(qū)的電離層結(jié)構(gòu)較為復(fù)雜，不僅具有較高的電子密度，且存在較大的時(shí)間和空間梯度，一直是電離層模式研究的關(guān)注對(duì)象.近年來(lái)，以電離層經(jīng)驗(yàn)?zāi)Ｊ交蚶碚撃Ｊ綖楸尘皥?chǎng)，吸收包含COSMIC在內(nèi)的多類觀測(cè)數(shù)據(jù)，構(gòu)建低緯同化模式成為電離層領(lǐng)域研究的重點(diǎn)和難點(diǎn)問(wèn)題.因此，有必要全面考察COSMIC掩星和IRI模型得到的電離層特征參數(shù)在低緯地區(qū)的精度，為未來(lái)電離層模式IRI進(jìn)一步完善以及利用IRI為背景場(chǎng)開展同化研究提供參考信息.由于測(cè)高儀垂測(cè)數(shù)據(jù)具有較高的精度，本文以磁赤道觀測(cè)站Jicamarca (11.95°S，76.8°W，地磁緯度為1°N)電離層測(cè)高儀垂直探測(cè)參量為真實(shí)值，評(píng)估2007—2013年COSMIC掩星和國(guó)際參考電離層模型IRI獲取的電離層特征參量NmF2和HmF2的精度.文中以下第二部分描述了數(shù)據(jù)資源和數(shù)據(jù)處理方法；第三部分探討了COSMIC掩星和IRI模型得到的電離層特征參量與垂測(cè)數(shù)據(jù)的相關(guān)性、不同太陽(yáng)活動(dòng)條件下的回歸分析結(jié)果，并進(jìn)一步統(tǒng)計(jì)分析了COSMIC掩星和IRI模型高估或低估垂測(cè)數(shù)據(jù)的變化特征.第四、五部分分別對(duì)結(jié)果進(jìn)行了討論和總結(jié).

2數(shù)據(jù)資源和數(shù)據(jù)處理方法

文中電離層垂測(cè)資料來(lái)自磁赤道觀測(cè)站Jicamarca,該臺(tái)站擁有長(zhǎng)期持續(xù)的數(shù)字測(cè)高儀觀測(cè)記錄，相關(guān)電離圖數(shù)據(jù)可以通過(guò)數(shù)字電離層圖數(shù)據(jù)中心下載獲取.利用SAO軟件包(http:∥ulcar.uml.edu/digisonde.html)能夠讀取測(cè)高儀原始探測(cè)數(shù)據(jù)，并對(duì)觀測(cè)結(jié)果進(jìn)行人工標(biāo)定，生成標(biāo)定后的數(shù)據(jù)文件，其時(shí)間分辨為15 min.由F2層臨界頻率foF2測(cè)量結(jié)果，根據(jù)以下數(shù)學(xué)表達(dá)式：

(1)

計(jì)算得到峰值密度NmF2的大小，單位為el·m-3.

COSMIC掩星探測(cè)資料來(lái)自于COSMIC數(shù)據(jù)存檔與分析中心CDAAC提供的后處理數(shù)據(jù)文檔ionProf.文檔包含單次掩星事件隨高度變化的電子密度、地理緯度和經(jīng)度、高度、世界時(shí)信息.文中分析高度在100 km以上的數(shù)據(jù)，以剔除電子密度剖面的負(fù)值觀測(cè)數(shù)據(jù).此外，由于地磁活動(dòng)或電離層中不規(guī)則體等因素可能造成電子密度剖面劇烈的擾動(dòng).為消除擾動(dòng)數(shù)據(jù)對(duì)分析結(jié)果的負(fù)面影響，COSMIC反演參數(shù)NmF2的選取采用以下標(biāo)準(zhǔn)(Guo et al., 2011)：

(2)

為全面分析COSMIC掩星和IRI模型在不同太陽(yáng)活動(dòng)條件下的性能，文中選取2007—2013年連續(xù)7年的觀測(cè)數(shù)據(jù).圖1描述了太陽(yáng)活動(dòng)指數(shù)F10.7的變化.由圖可見，2007—2009年為太陽(yáng)活動(dòng)低年，F(xiàn)10.7指數(shù)的最大值不超過(guò)100 SFU，平均值約為71 SFU，特別是2008—2009年F10.7指數(shù)明顯偏小，太陽(yáng)活動(dòng)進(jìn)入極小期；2011—2013年F10.7指數(shù)明顯升高，最大值約為190 SFU, 平均值約為117 SFU，太陽(yáng)活動(dòng)進(jìn)入活躍期.

圖1　太陽(yáng)輻射指數(shù)F10.7在2007—2013期間的變化Fig.1　Variation of F10.7 index during the period of 2007—2013

3結(jié)果分析

首先探討不同季節(jié)COSMIC掩星探測(cè)、IRI模型輸出特征參量與電離層垂直探測(cè)結(jié)果的相關(guān)性.其季節(jié)劃分如表1所示.圖2給出了2007—2013年電離層測(cè)高儀、COSMIC掩星和IRI模型不同季節(jié)電離層NmF2中值日變化的比較結(jié)果，同時(shí)分別給出了COSMIC掩星探測(cè)NmF2、IRI模型輸出NmF2與測(cè)高儀垂測(cè)參數(shù)NmF2的相關(guān)性.由圖可見，太陽(yáng)活動(dòng)低年，COSMIC掩星反演結(jié)果更接近于測(cè)高儀垂測(cè)值，相關(guān)性較高；IRI模型輸出NmF2與測(cè)高儀垂測(cè)參數(shù)NmF2的變化較為一致，但在1600 LT—2400 LT的輸出結(jié)果具有較為明顯的偏差(冬季除外).與太陽(yáng)活動(dòng)低年的觀測(cè)相比，太陽(yáng)活動(dòng)高年電離層峰值密度明顯增大，COSMIC探測(cè)值和測(cè)高儀垂測(cè)值的偏差有所增加，相關(guān)性降低；但I(xiàn)RI模型輸出NmF2與垂測(cè)參數(shù)NmF2的偏差減小，相關(guān)性略有增加.無(wú)論太陽(yáng)活動(dòng)極小年還是太陽(yáng)活動(dòng)高年，COSMIC反演參數(shù)NmF2、IRI輸出參數(shù)NmF2與垂測(cè)參數(shù)NmF2之間的相關(guān)系數(shù)大都在0.9以上.圖3與圖2類似，但給出的是電離層參數(shù)HmF2中值日變化的比較結(jié)果和相關(guān)性分析.COSMIC探測(cè)HmF2與垂測(cè)HmF2的變化基本一致，春秋兩季具有較高的精度，相關(guān)系數(shù)大都在0.9以上，但在夏冬兩季相關(guān)性降低，特別是太陽(yáng)活動(dòng)低年夏季午夜后(0000 LT—0500 LT)與垂測(cè)值的偏差增大，其相關(guān)系數(shù)約為0.66.與COSMIC掩星的分析結(jié)果類似，IRI模型輸出HmF2與垂測(cè)HmF2在春秋的相關(guān)性好于夏冬季，但白天時(shí)段與垂測(cè)結(jié)果有較大的偏差，以太陽(yáng)活動(dòng)低年的夏季最為顯著，相關(guān)系數(shù)降低至0.5以下.需要指出的是由于COSMIC某一時(shí)段觀測(cè)數(shù)據(jù)的缺失，導(dǎo)致相關(guān)系數(shù)不定值NaN的出現(xiàn)，但少量缺失數(shù)據(jù)應(yīng)當(dāng)不會(huì)影響文中的分析結(jié)果.

表1　季節(jié)劃分范圍

COSMIC掩星探測(cè)能夠給出全球范圍電離層特征參量的變化，但對(duì)于給定觀測(cè)點(diǎn)難以獲取連續(xù)的觀測(cè)很容易導(dǎo)致某些時(shí)刻數(shù)據(jù)的缺失.為全面分析COSMIC掩星和IRI模型在不同太陽(yáng)活動(dòng)條件下的性能，文中利用測(cè)高儀垂測(cè)參數(shù)分別對(duì)太陽(yáng)活動(dòng)低年2007—2009和太陽(yáng)活動(dòng)高年2011—2013不同季節(jié)所有COSMIC探測(cè)和IRI模型輸出參數(shù)與測(cè)高儀垂測(cè)值的中位數(shù)進(jìn)行回歸分析，其NmF2散點(diǎn)分布圖如圖4所示.圖中m為回歸直線(黑色)的斜率，斜率越接近于1，意味二者測(cè)量值越接近；c為相關(guān)系數(shù).太陽(yáng)活動(dòng)低年，COSMIC掩星和IRI模型得到的NmF2與測(cè)高儀垂測(cè)值相關(guān)系數(shù)的最大值約為0.98，最小為0.88.從圖上可以看出回歸直線的斜率都超過(guò)了1，意味著COSMIC掩星和IRI模型都不同程度高估電離層參數(shù)NmF2.與太陽(yáng)活動(dòng)低年的觀測(cè)相比，太陽(yáng)活動(dòng)高年COSMIC探測(cè)NmF2和測(cè)高儀垂測(cè)NmF2的相關(guān)性略有降低，但也都在0.8以上；IRI模型輸出參數(shù)NmF2與測(cè)高儀垂測(cè)NmF2的相關(guān)性則略有升高，所有季節(jié)都在0.9以上.二者的回歸直線斜率與太陽(yáng)低年的結(jié)果相比，其值明顯降低，位于0.76～1之間.與圖4類似，但圖5給出的是電離層參數(shù)HmF2的散點(diǎn)分布.由圖可見，COSMIC掩星和IRI模型在春秋兩季的精度高于夏冬兩季，在太陽(yáng)活動(dòng)高年得到的結(jié)果略好于太陽(yáng)活動(dòng)低年的結(jié)果.與COSMIC的探測(cè)結(jié)果相比，IRI模型輸出HmF2的精度偏低，特別是在太陽(yáng)活動(dòng)低年的夏季，其斜率和相關(guān)系數(shù)僅約為0.39和0.27.太陽(yáng)活動(dòng)低年，COSMIC掩星和IRI模型除夏季外同樣存在高估垂測(cè)參數(shù)HmF2的現(xiàn)象.

圖5　COSMIC掩星探測(cè)和IRI模型得到的HmF2與垂測(cè)HmF2在太陽(yáng)低年(上半圖)和太陽(yáng)高年(下半圖)的散點(diǎn)分布圖(m為回歸直線的斜率，c為相關(guān)系數(shù))Fig.5　Scatter plots of the HmF2 from COSMIC/ IRI and the HmF2 from ionosonde during the period of 2007—2013 (the upper pane denotes the solar minimum of 2007—2009, the lower pane denotes the solar maximum of 2011～2013; m is the slope of the regress line and c is the correlation coefficient)

進(jìn)一步將一天劃分為凌晨(0400 LT—0800 LT)、白天(0800 LT—1600 LT)、傍晚(1600 LT—2000 LT)和夜間(2100 LT—0400 LT)4個(gè)時(shí)段，分析不同時(shí)段內(nèi)COSMIC掩星探測(cè)和IRI模型輸出電離層參數(shù)與測(cè)高儀垂測(cè)參數(shù)的相關(guān)性.在太陽(yáng)活動(dòng)低年(2007—2009)和高年(2011—2013)條件下，峰值密度和峰值高度在各時(shí)段的相關(guān)性如表格2所示.從表中可見，COSMIC反演得到的NmF2與測(cè)高儀垂測(cè)值在太陽(yáng)活動(dòng)低年相關(guān)系數(shù)都高于0.9，太陽(yáng)活動(dòng)高年凌晨時(shí)段約為0.98，其他時(shí)段分布在0.82～0.89之間.IRI模型輸出NmF2與垂測(cè)NmF2在太陽(yáng)低年和高年的凌晨時(shí)間相關(guān)系數(shù)都在0.95以上，其他時(shí)段則分布在在0.83～0.89之間.除傍晚時(shí)段，太陽(yáng)活動(dòng)高年IRI模型得到的相關(guān)系數(shù)略高于太陽(yáng)活動(dòng)低年的計(jì)算結(jié)果.與NmF2的計(jì)算結(jié)果相比，COSMIC掩星和IRI模型得到的HmF2與垂測(cè)值HmF2之間的相關(guān)系數(shù)都明顯下降.COSMIC探測(cè)HmF2與垂測(cè)HmF2之間的最大相關(guān)系數(shù)發(fā)生在夜間，大小約為0.83，最小相關(guān)系數(shù)發(fā)生在白天或凌晨，大小約為0.50.太陽(yáng)高年IRI模型輸出參數(shù)HmF2與垂測(cè)HmF2得到的相關(guān)系數(shù)在夜間時(shí)段相關(guān)系數(shù)接近0.90，凌晨和白天在0.5以上；太陽(yáng)活動(dòng)低年除白天時(shí)段為0.51，其他時(shí)間都低于0.5.無(wú)論太陽(yáng)活動(dòng)低年還是高年，IRI模型在傍晚時(shí)段相關(guān)系數(shù)最低，大小約為0.22.

圖4和圖5中回歸直線的斜率在太陽(yáng)活動(dòng)低年大都明顯高于1，意味著COSMIC掩星和IRI模型都高估電離層特征參數(shù)NmF2和HmF2.分別利用太陽(yáng)活動(dòng)低年2007—2009和太陽(yáng)活動(dòng)高年2011—2013期間的數(shù)據(jù)，對(duì)不同地方時(shí)COSMIC掩星探測(cè)值、IRI模式輸出值與測(cè)高儀垂測(cè)值求差值，并統(tǒng)計(jì)正負(fù)差值的數(shù)量.圖6給出了COSMIC掩星和IRI模型得到的電離層參數(shù)NmF2高于垂測(cè)值或低于垂測(cè)值數(shù)量的條形分布圖，圖中黑色曲線為高于垂測(cè)值數(shù)量和低于垂測(cè)數(shù)量之間的差值.從圖6可以看出，太陽(yáng)活動(dòng)低年，COSMIC探測(cè)NmF2在夜間和凌晨前后的時(shí)間段內(nèi)大都低于垂測(cè)NmF2，白天(1000LT—1800LT)高于垂測(cè)NmF2，但冬季在夜間22時(shí)呈現(xiàn)稍微高估垂測(cè)值的現(xiàn)象.IRI模型輸出的NmF2約在午夜至凌晨前后低于垂測(cè)NmF2，其他時(shí)段(0800 LT—2300 LT)大都高于垂測(cè)值，但秋季0000LT—1400LT時(shí)間段內(nèi)呈現(xiàn)持續(xù)低估現(xiàn)象.從太陽(yáng)活動(dòng)高年期間的變化曲線可以看出，COSMIC掩星探測(cè)NmF2約在0600 LT—1000 LT低估垂測(cè)NmF2，但在午夜后(0000 LT—0500 LT) 與太陽(yáng)活動(dòng)低年正好相反，呈現(xiàn)高估NmF2的現(xiàn)象.此外，太陽(yáng)高年COSMIC掩星在春夏秋白天的不同時(shí)段內(nèi)探測(cè)值略微高于垂測(cè)值，夏季沒(méi)有明顯差別.IRI模型的估計(jì)情況與COSMIC探測(cè)有些相似，特別是在0000 LT—1000 LT時(shí)間內(nèi)都呈現(xiàn)了先高估后低估的現(xiàn)象.在夏冬兩季，兩者的變化也基本一致，且冬季1500 LT都略高估垂測(cè)值.但春秋季1000LT后的白天時(shí)段與COSMIC探測(cè)不同，IRI模型大都低估垂測(cè)值NmF2.與圖6類似，但圖7給出的是電離層參數(shù)HmF2的高估或低估變化.太陽(yáng)低年COSMIC掩星和IRI模型在午夜至凌晨前后均低估HmF2，其他時(shí)段大都呈現(xiàn)高估的現(xiàn)象.太陽(yáng)活動(dòng)高年COSMIC掩星和IRI模型在凌晨0500 LT左右低估HmF2，白天(1000 LT—1600 LT)則高估HmF2，在傍晚至夜間時(shí)段內(nèi)與垂測(cè)值的統(tǒng)計(jì)結(jié)果接近0.需要指出的是無(wú)論是電離層參數(shù)NmF2還是HmF2，太陽(yáng)活動(dòng)高年相對(duì)太陽(yáng)低年其數(shù)據(jù)量均偏少，這主要是由于COSMIC探測(cè)數(shù)據(jù)缺失造成的，未來(lái)有必要利用更豐富的探測(cè)數(shù)據(jù)驗(yàn)證結(jié)果的可靠性.

圖2　2007—2013年電離層峰值密度NmF2中值的季節(jié)變化和日變化(黑色實(shí)線為測(cè)高儀垂測(cè)值，叉號(hào)×為COSMIC探測(cè)值，圓圈○為IRI輸出值； C-S為COSMIC探測(cè)NmF2和測(cè)高儀垂測(cè)NmF2的相關(guān)系數(shù)，I-S為IRI輸出NmF2和測(cè)高儀垂測(cè)NmF2的相關(guān)系數(shù))Fig.2　Seasonal and diurnal variation of the ionospheric median NmF2 during the period of 2007—2013 (black solid line denotes digisonde measurements,× denotes COSMIC measurements and ○ denotes IRI outputs; C-S is the correlation coefficients between COSMIC-derived NmF2 and digisonde observations, I-S is the correlation coefficients between IRI-derived NmF2 and digisonde observations)

圖3　2007—2013年電離層峰值高度HmF2中值的季節(jié)變化和日變化(黑色實(shí)線為測(cè)高儀垂測(cè)值，叉號(hào)×為COSMIC探測(cè)值，圓圈○為IRI輸出值； C-S為COSMIC探測(cè)HmF2和測(cè)高儀垂測(cè)HmF2的相關(guān)系數(shù)，I-S為IRI輸出HmF2和測(cè)高儀垂測(cè)HmF2的相關(guān)系數(shù))Fig.3　Seasonal and diurnal variation of the ionospheric median HmF2 during the period of 2007—2013 (black solid line denotes digisonde measurements, × denotes COSMIC measurements and ○ denotes IRI outputs; C-S is the correlation coefficients between COSMIC-derived HmF2 and digisonde observations, I-S is the correlation coefficients between IRI-derived HmF2 and digisonde observations)

圖4　COSMIC掩星探測(cè)和IRI模型得到的NmF2與測(cè)高儀垂測(cè)NmF2在太陽(yáng)低年(上半圖)和太陽(yáng)高年 (下半圖)的散點(diǎn)分布圖 (m為回歸直線的斜率，c為相關(guān)系數(shù))Fig.4　Scatter plots of the NmF2 from COSMIC/IRI and the NmF2 from ionosonde during the period of 2007—2013 (the upper pane denotes the solar minimum of 2007—2009, the lower pane denotes the solar maximum of 2011—2013; m is the slope of the regress line and c is the correlation coefficient)

參數(shù)凌晨白天傍晚夜間COSMICNmF20.9798/0.98200.9099/0.81810.9494/0.86650.9416/0.8904IRINmF20.9548/0.97060.8261/0.86300.8781/0.83930.8705/0.8894COSMICHmF20.6764/0.50840.5045/0.70880.6633/0.64560.6109/0.8320IRIHmF20.3244/0.50470.5118/0.64160.2699/0.22650.4502/0.8997

圖6　COSMIC掩星和IRI模型得到的NmF2與垂測(cè)NmF2正負(fù)差值的數(shù)量分布(上半圖為太陽(yáng)活動(dòng)低年，下半圖為太陽(yáng)活動(dòng)高年, 黑色曲線為高于垂測(cè)值數(shù)量和低于垂測(cè)數(shù)量之間的差值；C-S為COSMIC探測(cè)值和測(cè)高儀垂測(cè)值之間的高估/低估結(jié)果，I-S為IRI輸出值和測(cè)高儀垂測(cè)值之間的高估/低估結(jié)果)Fig.6　Counts of the difference of the COSMIC/IRI-derived NmF2 and the ionosonde-derived NmF2 (the upper panel denotes the solar minimum, the lower panel denotes the solar maximum; black line denotes the differences between overestimation and underestimation; C-S denotes estimation between COSMIC anddigisonde measurements and I-S denotes estimation between IRI outputs and digisonde measurements)

圖7　COSMIC掩星和IRI模型輸出HmF2與垂測(cè)HmF2正負(fù)差值的數(shù)量分布(上半圖為太陽(yáng)活動(dòng)低年，下半圖為太陽(yáng)活動(dòng)高年，黑色曲線為高于垂測(cè)值數(shù)量和低于垂測(cè)數(shù)量之間的差值；C-S為COSMIC探測(cè)值和測(cè)高儀垂測(cè)值之間的高估/低估結(jié)果，I-S為IRI輸出值和測(cè)高儀垂測(cè)值之間的高估/低估結(jié)果)Fig.7　Counts of the difference of the COSMIC/IRI-derived HmF2 and the ionosonde-derived HmF2 (the upper panel denotes the solar minimum, the lower panel denotes the solar maximum; C-S denotesestimation between COSMIC and digisonde measurements and I-S denotes estimation between IRI outputs　and digisonde measurements)

4討論

與低緯地區(qū)觀測(cè)站Jicamarca垂測(cè)結(jié)果相比較，本文分析了不同太陽(yáng)活動(dòng)、不同季節(jié)、不同時(shí)間條件下COSMIC掩星和IRI模型得到的電離層特征參數(shù)NmF2和HmF2的精度，同時(shí)對(duì)COSMIC掩星和IRI模型高估或低估測(cè)高儀垂測(cè)值進(jìn)行統(tǒng)計(jì)分析.由COSMIC掩星和IRI模型得到的NmF2與垂測(cè)NmF2均具有較高的相關(guān)性，但COSMIC掩星和IRI模型得到的HmF2與垂測(cè)HmF2的相關(guān)性有所降低.與太陽(yáng)活動(dòng)低年的觀測(cè)相比，COSMIC掩星探測(cè)參數(shù)NmF2在太陽(yáng)高年與電離層垂測(cè)值得到的相關(guān)性略有降低，但由IRI模型輸出參數(shù)NmF2計(jì)算得到的相關(guān)性增大.COSMIC掩星和IRI模型凌晨前后低估垂測(cè)參數(shù)NmF2，白天時(shí)段則大都高估垂測(cè)參數(shù)NmF2，特別是太陽(yáng)低年尤為明顯.

Liu 等(2010)利用2007年Jicamarca垂測(cè)資料得到的結(jié)果表明COSMIC掩星探測(cè)參數(shù)NmF2和HmF2與測(cè)高儀垂測(cè)值的季節(jié)相關(guān)性分別在0.9和0.8以上；對(duì)比本文太陽(yáng)活動(dòng)低年分析結(jié)果，NmF2的相關(guān)性基本一致，HmF2的相關(guān)性略有降低，分布在0.7～0.9之間；文中給出的電離層參數(shù)在不同時(shí)段的相關(guān)性與Chuo等(2011)利用2006—2008觀測(cè)數(shù)據(jù)得到的結(jié)果大都吻合，但白天和傍晚時(shí)段HmF2稍低.Chuo等(2011)以Jicamarca的地理經(jīng)緯度為中心，選取空間距離不超過(guò)3°區(qū)域內(nèi)的COSMIC探測(cè)數(shù)據(jù)開展相關(guān)研究，文中同樣以觀測(cè)站位置為中心，但選取空間距離不超過(guò)5°的正方形區(qū)域內(nèi)的COSMIC掩星觀測(cè)數(shù)據(jù).此外，觀測(cè)數(shù)據(jù)的時(shí)間跨度也存在差別，因此相關(guān)結(jié)果的細(xì)微差異應(yīng)當(dāng)是合理的.COSMIC掩星基于球?qū)ΨQ假設(shè)條件，利用Abel積分反演得到的射線切點(diǎn)處的電子密度剖面實(shí)際上是整個(gè)掩星區(qū)域的平均值，當(dāng)電子密度的水平梯度較大時(shí)，容易造成較大的誤差(Yue et al., 2011).對(duì)于太陽(yáng)活動(dòng)高年，電離層電子濃度不僅顯著升高，時(shí)間和空間的梯度變化進(jìn)一步增大，致使COSMIC掩星反演NmF2的誤差增大，與垂測(cè)值之間的相關(guān)性降低.為分析不同太陽(yáng)活動(dòng)條件下COSMIC掩星探測(cè)電離層NmF2和HmF2與垂測(cè)結(jié)果的差異，進(jìn)一步統(tǒng)計(jì)分析了太陽(yáng)活動(dòng)低年(2007—2009)和太陽(yáng)活動(dòng)高年(2011—2013)電離層特征參數(shù)絕對(duì)偏差和相對(duì)誤差的分布，其結(jié)果如圖8所示.太陽(yáng)高年電離層參數(shù)NmF2位于±0.1×106el·cm-3范圍內(nèi)的絕對(duì)偏差和位于±30%范圍內(nèi)的相對(duì)誤差的百分比分別約為63%和51%，與太陽(yáng)活動(dòng)低年的結(jié)果相比分別下降29%和13%；均方根誤差RMSE和平均相對(duì)誤差增大，分別約為1.5×105el·cm-3和24.56%.太陽(yáng)高年電離層參數(shù)HmF2位于±25 km范圍內(nèi)的絕對(duì)誤差分布百分比和均方根誤差略遜于太陽(yáng)低年的結(jié)果，但位于±10%范圍內(nèi)的相對(duì)誤差分布和平均相對(duì)誤差則好于太陽(yáng)低年的結(jié)果.太陽(yáng)活動(dòng)高年COSMIC掩星反演參數(shù)NmF2的精度降低，但不同太陽(yáng)活動(dòng)水平對(duì)于HmF2的影響并不顯著.與Hu等(2014)在太陽(yáng)活動(dòng)高年(2011—2013)中國(guó)120°E地區(qū)的COSMIC掩星探測(cè)結(jié)果相比，NmF2和HmF2的精度高于低緯觀測(cè)站武漢和三亞的計(jì)算結(jié)果，其差異的主要原因可以歸結(jié)為中國(guó)低緯地區(qū)正好處于赤道北駝峰區(qū)域內(nèi)，電子濃度具有更大的梯度，可能導(dǎo)致較大的誤差.

圖8　COSMIC掩星探測(cè)參數(shù)和測(cè)高儀垂測(cè)參數(shù)的絕對(duì)誤差(左圖)和相對(duì)誤差(右圖)分布(灰色為太陽(yáng)活動(dòng)低年，黑色為太陽(yáng)活動(dòng)高年)Fig.8　Distribution of absolute error (left panel) and relative error (right panel) of ionospheric character parameters from COSMIC and ionosonde measurements (grey bar represents the solar minimum and black bar represents the solar maximum)

Yue等(2010)的研究指明白天時(shí)段Abel積分反演得到的射線切點(diǎn)正好位于赤道異常的駝峰區(qū)，由于球?qū)ΨQ假設(shè)可能會(huì)導(dǎo)致異常區(qū)內(nèi)(±10°—30°)的電子密度被低估，異常區(qū)附近區(qū)域(±10°)的電子密度被高估.從本文的結(jié)果可以看出，太陽(yáng)低年低估電離層參數(shù)NmF2主要發(fā)生在夜間和凌晨前后的時(shí)間，高估則發(fā)生在白天時(shí)段；這些現(xiàn)象與Yue等(2010)的分析結(jié)果基本一致，同時(shí)也與Liu等(2010)、Chuo等(2011)在太陽(yáng)低年的觀測(cè)結(jié)果相吻合.與太陽(yáng)活動(dòng)低年正好相反，太陽(yáng)活動(dòng)高年在0000 LT—0500 LT期間明顯高估電離層參數(shù)NmF2.太陽(yáng)活動(dòng)高年，太陽(yáng)輻射增強(qiáng)很可能導(dǎo)致凌晨時(shí)段赤道異常區(qū)內(nèi)的電子密度增強(qiáng)，從而導(dǎo)致NmF2的高估.

從IRI模型與電離層垂測(cè)參數(shù)的比較結(jié)果可以看出，NmF2的相關(guān)性高于HmF2，太陽(yáng)活動(dòng)高年得到的相關(guān)性高于太陽(yáng)低年的結(jié)果.與電離層垂測(cè)值相比，IRI模型大都高估電離層特征參數(shù)NmF2和HmF2，特別是太陽(yáng)活動(dòng)低年尤為顯著.2007—2009是太陽(yáng)活動(dòng)周期的低年，特別是2008/2009年不僅比以往低年的太陽(yáng)活動(dòng)程度要低，且持續(xù)時(shí)間長(zhǎng)，為1913年以來(lái)太陽(yáng)活動(dòng)的非尋常極小年(Bilitza et al., 2014).Lühr和Xiong(2010)比較了IRI2007和CHAMP/GRACE衛(wèi)星得到的電子密度剖面，結(jié)果發(fā)現(xiàn)IRI模型明顯高估兩種電子密度參數(shù)，其高估幅度在2008年和2009年分別約為50%和60%，在低緯白天時(shí)段更為顯著.Yue等(2013)比較了2008年IRI 2007和COSMIC掩星得到的斜向電子濃度總含量TEC, 結(jié)果發(fā)現(xiàn)在電離層F層高度范圍IRI模型均高估電子密度和電離層高度.Lühr和Xiong(2010)認(rèn)為在太陽(yáng)活動(dòng)極小年期間，F(xiàn)10.7指數(shù)不能完全準(zhǔn)確描述太陽(yáng)輻射情況導(dǎo)致IRI模型高估電離層參數(shù).IRI模型利用大量觀測(cè)數(shù)據(jù)，融合大氣物理模型參數(shù)，引入太陽(yáng)活動(dòng)指數(shù)和地磁指數(shù)，采用平滑、插值技術(shù)過(guò)濾了小尺度的擾動(dòng)，能夠較為準(zhǔn)確地反映電離層的“平均”變化.但I(xiàn)RI模型是基于數(shù)據(jù)驅(qū)動(dòng)的經(jīng)驗(yàn)?zāi)Ｊ?，如果缺乏某一區(qū)域或某一時(shí)段的數(shù)據(jù)，IRI模型輸出參數(shù)的精度就會(huì)降低.由于2008—2009期間太陽(yáng)活動(dòng)表現(xiàn)了不同以往低年的“異?！睒O小，電離層高度明顯下降，導(dǎo)致IRI模型高估電離層特征參數(shù).需要指出的是文中采用了IRI模型最新版本2012，由文中得到的分析結(jié)果，IRI 2012模型在太陽(yáng)活動(dòng)極小年仍然高估測(cè)高儀垂測(cè)參數(shù).與IRI 2007相比，IRI 2012主要在F1和F2中間層、極光E區(qū)、電子溫度和離子組成成分等方面進(jìn)行了有意義補(bǔ)充與改進(jìn)(Bilitza et al., 2014)，但低緯地區(qū)電離層NmF2和HmF2的估計(jì)精度仍需要進(jìn)一步提升.

5結(jié)論

利用2007—2013電離層測(cè)高儀磁赤道地區(qū)Jicamarca的垂測(cè)數(shù)據(jù)，探討了不同條件下COSMIC掩星探測(cè)和IRI 2012模型輸出參數(shù)NmF2和HmF2的性能，統(tǒng)計(jì)分析了COSMIC掩星和IRI模型在不同地方時(shí)高估或低估垂測(cè)參數(shù)的分布情況，得到了以下一些結(jié)論：

(1) COSMIC掩星探測(cè)NmF2、IRI模型輸出NmF2與測(cè)高儀垂測(cè)NmF2的季節(jié)、時(shí)段相關(guān)系數(shù)均不低于0.8，具有較好的相關(guān)性.

(2) COSMIC掩星探測(cè)HmF2、IRI模型輸出HmF2與測(cè)高儀垂測(cè)HmF2的相關(guān)性在春秋季節(jié)較高，其相關(guān)系數(shù)大都集中在0.8以上，但夏冬季節(jié)相關(guān)性則不同程度下降，其中由IRI模型在太陽(yáng)活動(dòng)低年夏季得到的相關(guān)系數(shù)最低，其值約為0.2.由COSMIC掩星計(jì)算得到的各時(shí)段HmF2相關(guān)系數(shù)大都在0.6以上，但由IRI模型計(jì)算得到的相關(guān)系數(shù)大多降至0.5以下，特別是傍晚時(shí)段相關(guān)系數(shù)約為0.2.

(3) 太陽(yáng)活動(dòng)低年COSMIC掩星探測(cè)NmF2與垂測(cè)NmF2之間的相關(guān)性高于太陽(yáng)活動(dòng)高年的相關(guān)結(jié)果，但太陽(yáng)活動(dòng)對(duì)于參數(shù)HmF2的相關(guān)結(jié)果并無(wú)明顯影響；太陽(yáng)活動(dòng)低年IRI模型輸出參數(shù)NmF2和HmF2與測(cè)高儀垂測(cè)參數(shù)的相關(guān)性低于太陽(yáng)活動(dòng)高年得到的結(jié)果.

(4) 太陽(yáng)活動(dòng)低年COSMIC掩星和IRI模型白天時(shí)段和夜間明顯高估電離層參數(shù)，凌晨前后則低估電離層參數(shù)；太陽(yáng)活動(dòng)高年兩者則在午夜后(0000 LT—0500 LT)明顯高估電離層NmF2，白天仍呈現(xiàn)高估HmF2的現(xiàn)象.

致謝感謝數(shù)字電離層圖數(shù)據(jù)中心和COSMIC數(shù)據(jù)存檔與分析中心CDAAC為本文提供數(shù)據(jù)資源.

References

Bilitza D. 2001. International reference ionosphere 2000. Radio Sci., 36(2): 261-275.

Bilitza D, Altadill D, Zhang Y L, et al. 2014. The international reference ionosphere 2012—a model of international collaboration. J. Space Weather Space Clim., 4: A07, doi: 10.1051/swsc/2014004.Chuo Y J, Lee C C, Chen W S, et al. 2011. Comparison between bottomside ionospheric profile parameters retrieved from FORMOSAT3 measurements and ground-based observations collected at Jicamarca. J. Atmos. Solar-Terr. Phys.,73(13):1665-1673. Ely C V, Batista I S, Abdu M A. 2012. Radio occultation electron density profiles from the FORMOSAT-3/COSMIC satellites over the Brazilian region: A comparison with Digisonde data. Adv. Space Res., 49(11): 1553-1562.

Guo P, Xu X, Zhang G X. 2011. Analysis of the ionospheric equivalent slab thickness based on ground-based GPS/TEC and GPS/COSMIC RO measurements. J. Atmos. Solar-Terr. Phys., 73(7-8): 839-846.Hu L H, Ning B Q, Liu L B, et al. 2014. Validation of COSMIC ionospheric peak parameters by the measurements of an ionosonde chain in China. Ann. Geophys., 32(10): 1311-1319. Lei J H, Syndergaard S, Burns A G, et al. 2007. Comparison of COSMIC ionospheric measurements with ground-based observations and model predictions: Preliminary results. J. Geophys. Res., 112(A7): A07308, doi: 10.1029/2006JA012240. Liu J Y, Lee C C, Yang J Y, et al. 2010. Electron density profiles in the equatorial ionosphere observed by the FORMOSAT-3/COSMIC and a digisonde at Jicamarca. GPS Solut., 14(1): 75-81, doi: 10.1007/s10291-009-0150-3.

Lühr H, Xiong C. 2010. IRI-2007 model overestimates electron density during the 23/24 solar minimum. Geophys. Res. Lett., 37(23): L23101, doi: 10.1029/2010GL045430.

Sun L F, Zhao B Q, Yue X A, et al. 2014. Comparison between ionospheric character parameters retrieved from FORMOSAT3 measurement and ionosonde observation over China. Chinese J. Geophys. (in Chinese), 57(11): 3625-3632, doi: 10.6038/cjg20141116.

Yue X A, Schreiner W S, Lei J H, et al. 2010. Error analysis of Abel retrieved electron density profiles from radio occultation measurements. Ann. Geophys., 28(1): 217-222, doi: 10.5194/angeo-28-217-2010.

Yue X A, Schreiner W S, Rocken C, et al. 2011. Evaluation of the orbit altitude electron density estimation and its effect on the Abel inversion from radio occultation measurements. Radio Sci., 46(1): RS1013, doi: 10.1029/2010RS004514.

Yue X A, Schreiner W S, Rocken C, et al. 2013. Validate the IRI2007 model by the COSMIC slant TEC data during the extremely solar minimum of 2008. Adv. Space Res., 51(4): 647-653.

附中文參考文獻(xiàn)

孫凌峰, 趙必強(qiáng), 樂(lè)新安等. 2014. 中國(guó)區(qū)域電離層垂測(cè)儀探測(cè)參量與COSMIC掩星反演結(jié)果比較研究. 地球物理學(xué)報(bào), 57(11): 3625-3632, doi: 10.6038/cjg20141116.

(本文編輯胡素芳)

基金項(xiàng)目國(guó)家自然科學(xué)基金(41104096)資助.

作者簡(jiǎn)介黃智,女,1973年生,副教授,博士,從事電波傳波及電離層物理相關(guān)研究.E-mail:wang86516@163.com

doi:10.6038/cjg20160701 中圖分類號(hào)P352

收稿日期2015-12-07，2016-01-11收修定稿

An analysis of ionospheric characteristic parameters retrieved from COSMIC and IRI at Jicamarca during the period of 2007—2013

HUANG Zhi1， YUAN Hong2

1SchoolofPhysicsandElectronicEngineering,JiangsuNormalUniversity,XuzhouJiangsu221116，China2AcademyofOpto-Electronics,ChineseAcademyofSciences,Beijing100094，China

AbstractTo evaluate comprehensively the performance of the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) and the International Reference Ionosphere (IRI), ionospheric characteristic parameters, the F2-peak density NmF2 and the F2-peak height HmF2from Digisonde over Jicamarca (11.95°S, 76.8°W, 1°N geomagnetic) are compared with data from the COSMIC radio occultation and IRI model during the different solar activity period of 2007—2013 in this study. The correlation of the digisonde observations and COSMIC/IRI measurements is analyzed emphatically, and then ionospheric parameters at Jicamarca being overestimated/underestimated by COSMIC and IRI are further discussed. The results are shown as follows: (1) Both the COSMIC-derived NmF2 and the IRI-derived NmF2 are in good agreement with the digisonde-derived measurements, and the correlation coefficients are more than 0.8. The correlation of the COSMIC-derived NmF2 in the solar minimum is greater than that in the solar maximum, and while the correlation of the IRI-derived NmF2 in the solar minimum is smaller than that in the solar maximum; (2) The HmF2 correlation between COSMIC/IRI and digisonde in spring and autumn is greater than that in summer and winter. Most seasonal and temporal correlation coefficients of COSMIC-derived HmF2 are more than 0.8 and 0.6, respectively. However, the correlation of the IRI-derived HmF2 decreased significantly, especially for the summer season in the solar minimum and dusk segment with a value of less than 0.3; (3) The ionospheric parameters are generally overestimated by COSMIC and IRI in the daytime and while they are underestimated from night to dawn in the solar minimum. Distinct from that in the solar minimum, the NmF2 are overestimated during the post-night (0000 LT—0500 LT) in the solar maximum. The current research could be useful for the ionospheric model improvement as well as data assimilation work based on the IRI model and COSMIC radio occultation data.

KeywordsDigisonde； COSMIC； IRI； Ionospheric character parameter； Correlation

黃智， 袁洪. 2016. 磁赤道地區(qū)2007—2013年COSMIC掩星反演和國(guó)際參考電離層模型輸出結(jié)果分析. 地球物理學(xué)報(bào)，59(7):2333-2343,doi:10.6038/cjg20160701.

Huang Z, Yuan H. 2016. An analysis of ionospheric characteristic parameters retrieved from COSMIC and IRI at Jicamarca during the period of 2007—2013. Chinese J. Geophys. (in Chinese),59(7):2333-2343,doi:10.6038/cjg20160701.