基于地基GPS研究臺風引起電離層TEC的變化特征

許九靖,柯福陽，趙興旺

(1.安徽理工大學 測繪學院,安徽 淮南 232001;2.南京信息工程大學 遙感與測繪工程學院,江蘇 南京 210044)

0 引 言

早期的大量研究表明,電離層與低層大氣氣象活動是存在相互聯系的,低層大氣氣象活動會引起電離層結構和物理特性的變化,產生小尺度的電離層擾動[1].臺風便是一種典型且復雜的低層大氣氣象活動.1958年Bauer[2]首次發現颶風過境時電離層會產生擾動:當颶風向臺站靠近時,F2層臨界頻率foF2開始增加;當颶風最接近臺站時,地面氣壓降到最低,foF2增加到最大;同時,他假設性地提出了大氣輻散/輻合模式是颶風產生電離層擾動的機理.1978年美國NASA的Hung等[3]利用高頻多普勒探測儀觀測到龍卷風爆發期間電離層F層存在重力波并發生中尺度擾動,并以此認為重力波可能是產生電離層擾動的主要原因.臺灣的Huang等[4]研究發現高頻多普勒探測儀對由臺風激發的能夠對電離層產生擾動的聲重力波有效探測性較差.但是,肖冠賽等[5]通過北京大學高頻多普勒臺站,對24次臺風事件進行了分析,其中有明顯擾動的記錄高達22次.1982年,沈長壽[6]統計分析發現:臺風對電離層foF2的變化有顯著影響,且因臺風的遠近而不同,臺風登陸前后電離層foF2明顯下降,并假設性地提出了湍流層頂移動是中低層大氣活動與電離層耦合的可能機理.劉依謀等[7]采用一維電離層模型驗證了沈長壽認為的臺風期間湍流層頂的抬升可能是臺風影響電離層F2區的一種有效機理的觀點.余濤等[8]通過對登陸廈門的一次臺風分析認為:由于臺風登陸前后,強烈的海氣、陸氣相互作用可能影響到電離層高度,從而導致電離層foF2、Es、和擴展F等參量的變化.隨著GPS技術的不斷發展,越來越多的學者開始借助GPS-TEC來研究臺風與電離層異常的耦合問題.毛田等[9]利用GPS臺站資料,研究了臺風“麥莎”對電離層電子濃度總含量(TEC)的影響.Rice等[10]分析GPS觀測的電離層TEC發現:Melor臺風期間TEC顯著增強,很可能是由臺風引起的.Lin等[11]認為臺風期間的電離層異常是由垂直聲重力波引起的,且利用非線性主成分分析法可以準確探測出TEC異常在地圖中的位置.Kong等[12]提出臺風激發的聲重力波引起了電離層行進式擾動(TID),并給出了擾動的速度和方向.

關于臺風過程電離層擾動形態特征,不同學者的研究結果仍存在一些差異.這主要是由臺風路徑上方監測的電離層數據精度不高、可靠性低、不連續及缺失造成的[13].本文利用中國氣象局布設在廣州、福建等地的GPS氣象臺站提取區域電離層TEC,并結合IGS提供的全球電離層地理(GIM)數據,借助滑動四分位距法[14]建立TEC的背景值與異常探測區間.綜合分析超級臺風莫蘭蒂登陸對廈門地區電離層TEC的影響.

1 數據與方法

本文以1614號超級臺風“莫蘭蒂”為研究對象,“莫蘭蒂”于2016年9月10日14時,在西北太平洋生產,9月13日晚間風速達到極值(75 m/s),9月15日3時5分登陸廈門,登陸時中心最大風速52 m/s.9月16日凌晨,臺風消散.圖1中粗實線表示臺風“莫蘭蒂”經過的路徑.

圖2為9月1日—9月19日空間環境指數F10.7、Kp和Dst的變化圖,虛線范圍內為臺風生成到消散的時間段.圖中所示時間段內的F10.7指數均小于100,尤其在臺風期間,F10.7指數均小于90.雖然在9月1日—9月10日,臺風生成之前的這段時間,Kp指數大部分都大于3,且Dst指數也均小于-30.但是臺風生成期間,Kp指數均小于3,且Dst指數絕大部分都大于-30(只有9月15日有一段時間小于-30).通常認為,F10.7指數小于100,太陽活動水平較低;Kp指數小于等于3和Dst指數大于-30,地磁活動較平靜.因此在臺風“莫蘭蒂”生成到消散的這一時間段內,太陽和地磁活動都較為平靜.

本文所用數據為中國氣象局布設于廣東、福建、廣西、云南和貴州的GPS氣象站資料,站點分布如圖1所示.其中正三角標記為提取廈門地區TEC所用的GPS站點,倒三角為提取參照地區TEC所用的GPS站點.已有研究表明TEC的緯度效應明顯,經度效應較小[15].即同緯度地區,經度變化對TEC產生的影響較小.因此本文利用廈門區域和參照GPS站點分別提取TEC,圖1中實心圓標記為臺風登陸廈門點(24.5°N,118.3°E)與參考點(24.5°N,110.3°E)位置.臺風期間兩地區的TEC變化差異有可能就是臺風對廈門地區TEC的影響.

基于雙頻GPS數據,利用雙頻P碼之差消去衛星、接收機的鐘差和對流層延遲可得:

(1)

帶入頻率f可得接收機至衛星路徑上的TEC觀測方程為

(2)

采用電離層單層模型(SLM)可將接收機至衛星路徑上的TEC投影到穿刺點的天頂方向:

(3)

式中:R為地球半徑;h為電離層薄層高度,一般取450 km;α為衛星的地平高度角.

因此,提取TEC的關鍵就是偽距差的優化和接收機、衛星硬件延遲的準確獲取.本文采用載波平滑偽距[16]的方法來得到較高精度的偽距差;然后利用區域球諧函數[17]方法進行接收機、衛星硬件延遲估計.

bnmsinms) .

(4)

聯合式(3)和式(4)即可建立穿刺點的VTEC方程.因為每個測站每天有約20 000個觀測值,聯立方程組,基于最小二乘原理即可求出接收機、衛星硬件延遲和球諧函數的系數[18].從而準確提取區域內任意點的VTEC.

2 結果與分析

2.1 本文方法獲取的電離層VTEC與IGS提供的VTEC精度對比

利用圖1所示站點的GPS雙頻數據,采用本文方法提取廈門和參考點上空的VTEC時間序列,并與IGS提供的GIM數據內插得到的VTEC時間序列比較.如圖3所示,兩種方法得到的VTEC變化趨勢一致,其均方差分別為4.12 TECU、3.58 TECU.因此采用本文方法,通過區域建模得到的VTEC的精度較好,能夠用于臺風影響電離層TEC的研究.

2.2 臺風引起電離層擾動分析

為分析臺風“莫蘭蒂”對廈門區域電離層TEC的影響,利用本文提取TEC的方法,建立廈門與同緯度不同經度的參照地區,從臺風生成到消失的TEC時間序列如圖4所示,虛線為臺風登陸時刻.

圖4顯示,廈門地區和參照地區的TEC時序圖趨勢基本相同,這進一步說明了TEC的經度效應不明顯.在臺風登陸前,兩地的電離層TEC都較背景值有所上升.尤其是9月13日的上午(世界時),上升最明顯,這一時刻臺風正好開始達到風速極值.隨著臺風的不斷臨近,兩地的TEC都開始趨于背景值.在臺風登陸后的第二天,兩地的TEC又較背景值有所上升.圖4(c)為兩地TEC與背景值之差,從中可以看出,廈門地區的TEC變化在臺風登陸前要大于參照地區,而在臺風登陸后的第二天里,廈門地區的TEC變化又小于參照地區.之后兩地變化趨于相同.這正說明廈門地區TEC受臺風影響,先上升后降低.

利用IGS提供的GIM數據建立臺風路徑區域的電離層TEC與其背景值之差的分布如圖5所示，背景值的建立方法為滑動四分位距法.

臺風登陸前一天,即世界時9月13日6：00,福建、臺灣地區上空電離層TEC值相較于其背景值有15～20 TECU的正擾動,且越靠近臺風中心,擾動強度越大.這時臺風開始影響臺灣島,且臺風風速達到極值(75 m/s).同一天的12:00,TEC異常值的范圍在12～15 TECU.9月14日6:00,臺風登陸臺灣島,此時臺風影響區域上空的電離層TEC相較其背景值的正擾動為1 TECU左右.而12:00時的TEC相較其背景值出現了1 TECU左右的負擾動,這時臺風已經離開臺灣島.

世界時9月14日18:00臺風登陸福建廈門,此時福建、臺灣地區上空電離層TEC值并沒有出現較大的擾動,其與背景值的差值在1～2 TECU.

借助滑動四分位距法,以1.5倍四分位距為限差,對2016年9月13日全球電離層TEC異常進行探測.如圖6所示,6:00時臺風周邊區域上空的電離層TEC異常已經形成,10:00和12:00異常有所降低.16:00和18:00時,TEC異常降低到8～10 TECU,圖6中，當臺風區域上空產生電離層擾動時，其赤道共軛區域也產生了電離層擾動，但是持續時間較短，形成原因還需要進一步的深入研究.

綜上,臺風“莫蘭蒂”登陸臺灣島前一天,電離層TEC表現為較強的正擾動,隨著臺風的臨近,TEC逐漸恢復到背景值,之后又出現了小尺度的負擾動.

2.3 可能機制分析

當臺風臨近臺灣島時,臺灣地區和廈門地區上空電離層都出現了強烈的正擾動;當臺風登陸廈門時,電離層TEC并沒有像臺風臨近臺灣島那樣出現明顯異常.一方面是由于臺風先經過臺灣島后登陸廈門,所以廈門上空電離層TEC在臺風臨近臺灣島時就發生了異常擾動.另一方面,當臺風臨近臺灣島時,強烈的風暴遭遇臺灣中央山脈的阻攔,為聲重力波的產生創造了一個絕佳的條件.之后,聲重力波向上傳播,并影響到電離層高度.而廈門地區的地形條件沒有起到這樣的阻攔作用.

3 結束語

本文基于雙頻GPS數據提取的區域電離層TEC,結合IGS提供的GIM數據,共同分析了臺風“莫蘭蒂”對TEC的影響.得出結論:

1)臺風登陸廈門前一天,即世界時9月13日,廈門區域電離層TEC相較于其背景值和同緯度的參考區域都出現了明顯的正異常擾動.

2)臺風引起的區域電離層異常擾動,在其臨近臺灣島前就已經出現.9月13日的6：00異常值在15～20 TECU之間.隨著時間的推移,異常值逐漸降低,10小時后,異常消失;在世界范圍內,9月13日這一天,只有臺風影響區域及其赤道共軛區出現了電離層TEC的異常.這一天的太陽和地磁活動都比較平靜,所以此次電離層異常擾動應該是臺風引起的.考慮到臺灣島的特殊地形(中央山脈縱貫全島,長約320 km,有62座山峰高度在3 000 m以上),當臺風臨近時,強烈的對流風暴遭遇中央山脈的阻攔,由此產生的聲重力波向上傳播,造成了電離層的異常擾動.