國際電離層模型磁暴模式的適用性分析

張 嘯，聞德保，2，湯佳明

國際電離層模型磁暴模式的適用性分析

張 嘯1，聞德保1，2，湯佳明1

(1.長沙理工大學交通運輸工程學院，長沙 410004；2.中國科學院大地測量與地球動力學國家重點實驗室，武漢 430077)

針對IRI2012磁暴模式在中國大陸強太陽活動程度和磁擾條件下的適用性問題，利用該模型模擬了電離層TEC及自由電子密度在2003年8月磁暴發生前后的時空變化，并與GPS-TEC 、測高儀數據及GPS電離層層析重構結果作了比較，結果表明：IRI2012磁暴模式能夠反映強太陽活動程度下TEC的分布特征，模擬精度可達±2.3 TECU～±7.4 TECU，并隨緯度降低而變大，TEC周日峰值較實測數據提前2 h出現，TEC模擬值在白天偏高，在夜間則偏低；F2層自由電子密度模擬值在峰值密度和峰值高度上與武漢測高儀數據有較好的一致性，但在磁擾發生的條件下，IRI2012磁暴模式由于不能有效的反映電子密度剖面赤道異常區域的北向移動以及赤道異常核在緯度方向的收縮和在高度方向的擴展，使得TEC模擬精度也出現較大衰減，最大達±14.4 TECU。

IRI2012；磁暴模式；電離層；差異性

0 引言

國際電離層參考模型(the international reference ionosphere，IRI)是由空間研究委員會(the committee on space research，COSPAR)和國際無線電科學聯合會(the international union of radio science，URSI)于二十世紀六十年代末共同建立的電離層經驗模型，可計算距地面50～2000 km空間范圍內電子密度、電子溫度、離子構成、離子溫度等電離層物理參量的月平均狀態，目前已經更新到2012版本，并且包含磁暴模式。隨著以無線電測量為主要手段的衛星導航、遙感、遙測的廣泛應用，電離層時空物理特性對電離層延遲的建模和預測顯得越來越重要，IRI作為主要電離層經驗模型之一，匯集了全球數百個電離層觀測站及衛星觀測資料以及多個大氣參數模型，反映了電離層在平均狀態下的重要物理特征[1～3]。國內外學者針對IRI模型作了大量的工作[4～6]，IRI模型也被廣泛應用于電離層研究當中，文獻[7]利用IRI2007預測電子密度剖面，得出太陽活動高年和活動低年F2層峰值密度在白天和夜晚的高度范圍；文獻[8]利用非相干散射雷達觀測數據同IRI2001進行比較，表明太陽活動低年IRI2001能較好的反映電離層電子密度隨時間和高度的變化，而在太陽活動高年500 km以下總電子含量( total electron content，TEC)達到最高值以后IRI2001計算的電離層頂電子密度過大，使得TEC明顯高于雷達數據；文獻[9]比較了IRI、NeQuick和Klubchar三種電離層模式，認為三種模式均能較好的模擬出電離層TEC背景特征，而且在電離層延遲改正中IRI精度較其他兩者高；文獻[10]利用IRI2012，IRI plas計算了歐洲中緯度區域的垂直總電子含量(vertical total electron content，VTEC)并與由全球定位系統(global positioning system，GPS)獲取的總電子含量(total electron content by GPS，GPS-TEC)作了比較分析；文獻[11]利用IRI2007作為背景場，將單基站TEC地圖精度在±0.25 個總電子含量單位(total electron content unit，TECU；相當于1016el/m2)范圍內的比例提高到70%以上；文獻[12]利用IRI2012模擬的電子密度作為初值，加快電離層層析迭代收斂速率，以拉布拉斯算子抑制結果對初值精度的依賴，實現了一種新的電離層層析算法。

鑒于IRI模型的研究和應用需求，本文針對IRI2012磁暴模式在中國區域強太陽活動程度和磁擾條件下的適用性問題，模擬了2003年8月磁暴中國區域及武漢上空的TEC分布、TEC周日變化及自由電子密度剖面，并與GPS電離層層析結果、GPS-TEC及測高儀觀測的電子密度(Ionosondes -electron density，Ionosondes-Ne)作了對比，分析了IRI2012磁暴模式與電離層實況的差異特征，以期利用該差異探求一種以IRI模型為背景的電離層磁暴異常探測模式。

1 數據準備

本文所采用的數據為2003年8月磁暴期間電離層自由電子總含量以及電離層自由電子密度，主要包括：IRI模型計算的總電子含量(total electron content calculated by IRI，IRI-TEC)是IRI自由電子密度(electron density calculated by IRI，IRI-Ne)在垂直路徑上的積分值，IRI2012自由電子密度在不同層別則通過不同的模型計算得到，本文中頂部電子濃度選用NeQuick模型，底部厚度參數選擇B0ABT，F2層峰值模型選用CCIR，并且開啟磁暴模式，此外，太陽和電離層指數(ig_rz.dat)、地磁場系數(dgrfyear.dat計算，igrfyear.dat預測)、F10.7(日平均、81日滑動平均、365日滑動平均)及地磁指數(apf107.dat)、F2層峰值模型系數(CCIR%%.ASC)均采取文件形式讀入，并在igrf.for、irifun.for、irisub.for中修改相應文件存放路徑。GPS-TEC由國際全球衛星導航系統服務組織(International GNSS Service，IGS)提供的TEC數據，非格網點采用二維雙線性內插得到。

圖1 2003年8月磁暴前后地磁擾動情況

電離層層析-電子密度(computerized ionospheric tomography-electron density，CIT-Ne)是基于中國陸態網GPS數據，利用組合電離層層析算法[13]計算得到的電離層自由電子密度。Ionosondes-Ne為電離層觀測站數字測高儀數據。

2 IRI2012模擬數據與電離層實測數據的比較

2.1 磁暴前后電離層TEC的空間分布特征

圖2為中國區域電離層IRI-TEC與GPS-TEC在磁暴前一天(17日)的分布圖，其間F10.7日平均值達到122太陽通量單位(solar flux unit，sfu)，屬于強太陽活動程度，IRI-TEC能較好的反映出TEC北半球赤道異常區域(20°N左右)的結構、大小和隨太陽高度角變化的TEC峰值位置，以及中國區域TEC梯度南向逐漸增大的變化特性，但IRI-TEC整體上較GPS-TEC在白天時段偏高，在夜間時段則偏低。

由于赤道電急流的加劇，產生與地磁場相反的強電磁場，造成全球性地磁擾動，改變了赤道異常區域F2層自由電子的沉降位置，加之太陽高能離子流加熱兩極大氣，形成赤道向中性風[14]，造成赤道異常區域南移或北移，在局部區域便體現為正相或負相暴[15]；18日的磁暴在中國區域是正相暴，赤道異常區域北移，中國中部南部GPS-TEC升高，梯度變大，北部中高緯區域GPS-TEC則呈現出負向擾動。而IRI-TEC除了在赤道異常區域較前日增大以外，未出現上述GPS-TEC的分布變化特征。

圖2 地磁平靜日IRI-TEC與GPS-TEC分布的比較

圖3 磁暴日IRI-TEC與GPS-TEC分布的比較

2.2 IRI-TEC與GPS-TEC的周日變化差異

為比較IRI-TEC與GPS-TEC周日變化差異，本文選取中國區域相鄰經度鏈上，緯度從高而低的北京、武漢、廣州3個測站為研究對象，如圖4(a)所示，在太陽活動程度較高但地磁較為平靜的條件下，IRI-TEC能夠反映電離層TEC周日變化趨勢，但TEC峰值較GPS-TEC提前約2 h左右，且峰值大小在北京區域較GPS-TEC低，在武漢區域則與之相當，而在廣州則偏高；IRI-TEC的模擬精度隨著緯度降低而逐漸增大，依次為±2.3 TECU、±4.1 TECU、±7.4 TECU；如圖4(b)，當磁暴發生，由于赤道異常區域北移，GPS-TEC急劇增大，峰值出現時間與IRI-TEC基本一致，但IRI-TEC的模擬精度也隨磁暴發生而變低，分別為±9.3 TECU、±14.4 TECU、±8.1 TECU，基本失去了電離層TEC的預報功能。

圖4 IRI-TEC與GPS-TEC的周日差異

2.3 IRI-Ne剖面與測高儀數據的比較

圖5展示了武漢站電離層測高儀觀測和IRI2012模擬的自由電子密度在白天(01 UT)和夜晚(13 UT)隨高度變化的變化情況其中，UT為世界時(universal time)。在太陽活動程度較高但地磁平靜的條件下，IRI2012模型白天F2層峰值密度約6.4×1011el/m3，峰值高度在250 km左右，與測高儀觀測的結果有較好的一致性(如圖5(a))；在夜晚IRI2012與測高儀在峰值高度上保持一致，在峰值密度上略有差異(如圖5(b))，總體上能較為真實的反映電離層的電子密度變化。隨著磁暴的爆發，IRI2012的F2層層頂計算模型數值偏低，且峰值高度在250 km左右，低于測高儀350 km的測量值(如圖5(c))；在夜晚，雖然峰值高度與測高儀保持一致300 km左右(如圖5(d))，但峰值密度明顯小于測高儀觀測值，這主要與IRI2012中F2層峰值高度模型和峰值密度模型相關。

圖5 地磁平靜和擾動條件下IRI-Ne剖面與測高儀數據的比較

2.4 IRI-Ne經度剖面與CIT-Ne的比較

圖6展示了在地磁平靜條件下，利用電離層層析組合算法得到的電子密度重構結果與IRI2012的模擬結果在114°經度鏈上的時序剖面圖，從圖6各子圖可看出，IRI2012的模擬結果與GPS電離層層析得到的結果具有較好的一致性，赤道異常核01 UT開始高度為350 km，07 UT增加到400 km，13UT又降到350 km；從北緯10°到北緯50°的區域，赤道異常核的層頂具有明顯的傾斜結構，且IRI-Ne較CIT-Ne要緩，這與所用的地磁模型有關。

圖7展示的是磁擾條件下在114°經度鏈上電子密度的時序剖面圖，由對應子圖可知CIT-Ne相對于IRI-Ne有明顯的異常特征，如圖4(a)，赤道異常區域的傾斜結構被破壞，在北緯30°到40°之前開始向高低兩處擴展，至07 UT(如圖4(b))，赤道異常核收縮并向北移動到而北緯20°左右處，之后不斷向高低兩個方向延伸(如圖4(c))，擴大到100～700 km的高度范圍，這也說明了3.2中IRI-TEC的模擬精度在武漢區域最差的原因；然而IRI-Ne則與前一日并無明顯差異。

圖6 地磁平靜日IRI電子密度剖面與GPS電離層層析的比較(1012 el/m3)

圖7 地磁平靜日IRI電子密度剖面與GPS電離層層析的比較(1011 el/m3)

3 結束語

利用IRI2012磁暴模式模擬了強太陽活動程度下，磁暴期間TEC及自由電子密度的分布和變化，通過與GPS-TEC、Ionosondes-Ne、CIT-Ne的對比分析，得出如下結論：

1)在地磁平靜條件下，IRI2012能夠反映強太陽活動程度下TEC分布隨太陽高度角變化的變化特性及經向梯度變化特征，包括北半球赤道異常區域的位置和大小；在TEC峰值模擬計算上比GPS-TEC實測數據提前2 h，計算精度最高為±2.3 TECU，并隨緯度的減小而降低至±7.4 TECU；模擬值在白天較GPS-TEC高，夜間則偏低；

2)IRI2012在 F2層峰值高度及峰值密度的模擬計算上與武漢測高儀數據具有較好的一致性，但在磁暴條件下，由于赤道異常區移動即正負向擾動的不確定性，IRI2012模式F2層頂模型不能很好的描述Ne白天在400～1 000 km的異常變化，夜晚F2層峰值高度與測高儀一致，但峰值密度偏低4.7×1011el/m3；

3)在自由電子密度剖面模擬方面，IRI2012能夠反映強太陽活動程度下的自由電子密度的時序變化，磁暴時CIT-Ne赤道異常核的北向移動，緯度向收縮及高度方向的擴張則未能有效反映，造成IRI-TEC較大的精度衰減。

[1] ROBERT F，PFAFF R F.The Near-Earth Plasma a Environment[J].Space Science Reviews，2012，168(1-4)：23-112.

[2] FORBES J M，PALO S E，ZHANG Xiao-li.Variability of the Ionosphere[J].Journal of Atmospheric and Solar-terrestrial Physics，2000，62(8)：685-693.

[3] 蔡超軍，黃江，鄧柏昌，等.IRI-2007預測TEC在廣州區域的適用性分析[J].空間科學學報，2013，33(3)：277-285.

[4] BILITZA D.International Reference Ionosphere 2000[J].Radio Science，2001，36 (2)，261-275.

[5] BILITZA D.International Reference Ionosphere 2002：Examples of Improvements and New Features[J].Advance in Space Research，2003，31(3)，757-767.

[6] BILITZA D，REINISCH B W.International Reference Ionosphere 2007.Improvements and New Features[J].Advance in Space Research，2008，42(4)：599-609.

[7] OLADIPO O A，ADENIYI J O，RADICELLA S M，et al.Variability of the Ionospheric Electron Density at Fixed Heights and Validation of IRI-2007 Profile’s Prediction at Ilorin[J].Advance in Space Research，2011，47(3)：496-505.

[8] 蔡紅濤，馬淑英，SCHLEGEL K.高緯電離層氣候學特征研究—EISCAT雷達觀測及與IRI模式的比較[J].地球物理學報，2005，48(3)：471-479.

[9] 方涵先，翁利斌，楊升高.IRI、NeQuick和Klobuchar模式比較研究[J].地球物理學報，2005，48(3)：471-479.

[10]KENPANKHO P，SUPNITHI P，NAGATSUMA T.Vertical TEC Representation by IRI 2012 and IRI Plas Models for European Midlatitudes[J].Advance in Space Research，2013，52(10)：1820-1826.

[11]翁利斌，方涵先，楊升高，等.基于IRI背景場的單站電離層TEC地圖重構技術[J].空間科學學報，2011，31(4)：453-458.

[12]聞德保，呂慧珠，張嘯.電離層層析重構的一種新算法[J].地球物理學報，2014，57(11)：3611-3616.

[13]聞德保.基于GNSS的電離層層析重構算法及其應用[M].北京：測繪出版社，2013.

[14]FULLER-ROWELL T J，RICHMOND A D，MARUYAMA N.Global Modeling of Storm-time Thermospheric Dynamics and Electrodynamics[C]//KINTNER Jr PM，COSTER A J，FULLER-ROWELL T J，et al.Midlatitude Ionospheric Dynamics and Disturbances.New York：John Wiley & Sons，Inc，2008：187-200.doi：10.1029/181GM06.

[15]XIA Chun-liang，WAN Wei-xing，YUAN Hong，et al.Analysis of the Intense Magnetic Storm of July，2000 and of October，2003 Using the Technique for Nowcasting of GPS TEC Data[J].Chinese Journal of Space Science，2005，25(4)：259-266.

The Applicability Analysis of International Reference Ionosphere in Storm Mode

ZHANGXiao1，WENDe-Bao1，2，TANGJia-Ming1

(1.School of Traffic and Transportation Engineering，Changsha University of Science & Technology， Changsha 410004，China；2.State Key Laboratory of Geodesy and Earth’s Dynamics，Chinese Academy of Science，Wuhan 430077，China)

For the applicability over China of IRI2012(the international reference ionosphere )in storm mode with strong solar activity or magnetic disturbance，the model is used to simulate the temporal and spatial variation of the ionospheric TEC(total electron content) and free electron density during magnetic storms in August 2003，and compared with the GPS-TEC，ionosondes data and reconstruction of ionospheric electron density by CIT(computerized ionospheric tomograph)，the results showed that：IRI2012 in storm mode can reflect the distribution of TEC over China with strong solar activity，the simulation accuracy is up to ± 2.3 TECU ～ ± 7.4 TECU increasing with latitude decline.In addition，TEC analog values is higher in daytime and lower at night compared with the measured data，its diurnal peak appears two hours earlier than GPS-TEC.The simulation of peak density and peak height of electron density in F2 layer is in well agreement with the ionosondes in Wuhan，but in the condition that magnetic disturbance occurs，the analog values of electron density profiles by IRI2012 in storm mode are unable to effectively reflect such changes of the area of equatorial anomaly that integral movement toward the north，contraction of nuclear in the latitude direction and expansion in the height direction，as a result the simulation accuracy of TEC make larger attenuation over this region，which is up to ± 14.4 TECU.

IRI2012；Storm mode；Ionosphere；Difference

張嘯，聞德保，湯佳明.國際電離層模型磁暴模式的適用性分析[J].導航定位學報,2015,3(3):89-94.ZHANG Xiao，WEN De-Bao，TANG Jia-Ming.The Applicability Analysis of International Reference Ionosphere in Storm Mode[J].Journal of Navigation and Positioning,2015,3(3):89-94.

10.16547/j.cnki.10-1096.20150318.

2015-05-18

國家自然科學基金(41174001)；湖南省杰出青年基金(14JJ1021)；教育部科學技術類項目(213028A)；湖南省教育廳重點項目(12A002)；湖南省研究生科研創新基金資助(CX2015B350)。

張嘯(1989—)，男(苗族)，貴州安順人，碩士生，主要從事電離層層析成像與電離層天氣的研究。

聞德保(1974—)，河南光山人，教授，博士，主要從事電離層層析成像、GNSS氣象學和無線電掩星技術等方面的研究。

P228

A

2095-4999(2015)-03-0089-06