利用FY-3C衛(wèi)星GNSS掩星數(shù)據(jù)分析中國區(qū)域?qū)α鲗禹攨?shù)變化

郭佳賓 金雙根,2

1 南京信息工程大學(xué)遙感與測繪工程學(xué)院，南京市寧六路219號，210044 2 中國科學(xué)院上海天文臺，上海市南丹路80號, 200030

對流層是地球大氣層中密度最大的一層，占據(jù)了整個大氣層75%的質(zhì)量。大氣層中90%的水汽集中在對流層，且對流頻繁，因此，眾多天氣現(xiàn)象都發(fā)生在對流層中[1]。對流層頂是平流層和對流層的過渡區(qū)域，在這個區(qū)域?qū)α鲗优c平流層頻繁交換氣團(tuán)、水汽和能量。世界各國學(xué)者利用探空氣球[2]、MST雷達(dá)[3]以及數(shù)值天氣預(yù)報再分析(ERA)等手段與資料，針對對流層頂開展了廣泛研究，也取得一定成果。但傳統(tǒng)探空手段均無法提供較高垂直分辨率的對流層頂結(jié)構(gòu)數(shù)據(jù)，且成本高，已經(jīng)無法滿足更深層次的對流層頂結(jié)構(gòu)與氣候變化關(guān)系的研究需求。

GNSS無線電掩星技術(shù)具有較高的垂直探測分辨率以及全球覆蓋的特點，有效打破了探空氣球、MST雷達(dá)等觀測手段的局限，被廣泛應(yīng)用于地球大氣的三維立體探測。針對眾多無線電掩星產(chǎn)品，學(xué)者們開展了大量科學(xué)實驗。劉艷等[4]利用COSMIC/GPS掩星折射率資料研究得出全球海洋邊界層頂高度的季節(jié)變化、年際變化和日變化的氣候?qū)W特點，徐寄遙等[5]利用COSMIC數(shù)據(jù)分析全球?qū)α鲗禹敎囟群透叨鹊淖兓卣?，劉久偉等[6]利用COSMIC/GPS 掩星干溫及干壓資料探測對流層頂?shù)母叨龋珺orsche等[7]利用CHAMP掩星觀測資料分析熱帶對流層頂溫度和高度的分布及變化特征，Schmidt等[8]利用CHAMP掩星觀測資料分析對流層頂高度、溫度、氣壓等參數(shù)的空間變化及年際變化等。這些研究基本上都是利用國外的GNSS掩星數(shù)據(jù)，且數(shù)據(jù)量不足。

2013-09我國發(fā)射FY-3C衛(wèi)星，其搭載的GNOS載荷是國內(nèi)第一個星上GNSS無線電掩星探測儀，也是國際上首臺兼容北斗導(dǎo)航系統(tǒng)(BDS)和全球定位系統(tǒng)(GPS)的掩星探測儀[9]。FY-3C衛(wèi)星每天可以分別得到200次北斗掩星事件以及500次GPS掩星事件[10]，這為氣象研究提供了豐富的高精度、全方位、全天時觀測資料。為了驗證我國首個無線電掩星探測儀的數(shù)據(jù)價值，本文使用國家衛(wèi)星氣象中心(http://satellite.nsmc.org.cn)提供的2017-03-01～2018-02-28期間GNSS無線電掩星附加相位數(shù)據(jù)，通過ROPP軟件[11]反演得到大氣溫度廓線，并分析中國區(qū)域?qū)α鲗禹攨?shù)隨緯度的分布狀況以及季節(jié)性變化特征。

1 觀測數(shù)據(jù)與方法

1.1 觀測數(shù)據(jù)

無線電掩星觀測技術(shù)可以完成對大氣剖面高垂直分辨率的觀測，其探測高度可以從地面延伸到60 km高度處，覆蓋整個對流層與平流層。但是每一條掩星廓線的位置是由FY-3C衛(wèi)星和被掩GNSS衛(wèi)星的位置決定的，分布比較稀疏。圖1(a)展示了2017-03-01～03-07掩星事件在中國區(qū)域的分布狀況，可以看到，在整個中國區(qū)域只有較少的掩星數(shù)據(jù)。圖1(b)展示了2017年整個3月份的掩星事件在中國區(qū)域的分布狀況，此時整個中國區(qū)域覆蓋了大量的掩星事件，且分布較為均勻，這為我們使用掩星數(shù)據(jù)分析中國區(qū)域?qū)α鲗禹攨?shù)提供了思路。

圖1 中國區(qū)域掩星事件分布Fig.1 Distribution of occultation events in China

本文選取2017-03-01～2018-02-28共計148 349條Level-1級數(shù)據(jù)的附加相位數(shù)據(jù)，借助掩星處理軟件ROPP計算得出大氣溫度廓線。

1.2 中國區(qū)域時空特征劃分

我國地域遼闊，東西、南北跨度大，且具有四季分明的特點。因此，在空間上，使用2°×2°的空間分辨率將整個中國區(qū)域劃分為608個網(wǎng)格，然后統(tǒng)計落在每個網(wǎng)格內(nèi)的掩星廓線數(shù)量，計算出每個網(wǎng)格點的對流層頂高度及溫度的平均值；在時間上，考慮到GNOS掩星每天的掩星點落在中國區(qū)域較少的狀況，通過分析不同月份的對流層頂參數(shù)變化，進(jìn)而得到對流層頂參數(shù)的季節(jié)性變化情況。

1.3 掩星數(shù)據(jù)精度驗證與估計對流層頂高度

在使用FY-3C掩星數(shù)據(jù)計算對流層頂高度時，首先要對其進(jìn)行精度驗證，以保證結(jié)果的準(zhǔn)確性。本文分別使用探空氣球數(shù)據(jù)以及ERA5數(shù)據(jù)對掩星數(shù)據(jù)反演溫度剖面精度進(jìn)行驗證。為了保證實驗數(shù)據(jù)的普遍準(zhǔn)確度，選擇2018-09-01～09-30的FY-3C掩星數(shù)據(jù)進(jìn)行驗證。

在與探空氣球數(shù)據(jù)比較時，首先對掩星數(shù)據(jù)和探空氣球數(shù)據(jù)進(jìn)行時空匹配。匹配標(biāo)準(zhǔn)為：二者觀測的時間差值在2 h之內(nèi)，水平距離在150 km之內(nèi)?；谠摌?biāo)準(zhǔn)匹配后，共得到134對匹配數(shù)據(jù)。在134對數(shù)據(jù)中，由于探空氣球數(shù)據(jù)文件在記錄時存在較多的以數(shù)值9999或者-9999代替缺省值，故需要對存在缺省值的探空數(shù)據(jù)進(jìn)行剔除。剔除后，整個9月一共得到89組匹配數(shù)據(jù)。在與ERA5數(shù)據(jù)進(jìn)行比較時，為計算方便，選擇2018-09每一天中午12：00的ERA5溫度數(shù)據(jù)；然后對掩星數(shù)據(jù)和ERA5溫度數(shù)據(jù)進(jìn)行時間和空間匹配，匹配標(biāo)準(zhǔn)與探空氣球一致；基于該標(biāo)準(zhǔn)匹配后，整個9月一共得到30組匹配數(shù)據(jù)。

國際上對于對流層頂高度的計算，主要有最冷點對流層頂與溫度遞減率對流層頂兩種方法。其中，溫度最低點對流層頂法(cold point tropopause，CPT)主要是依據(jù)溫度垂直廓線上最冷點對應(yīng)的高度(該高度所對應(yīng)的溫度稱為對流層頂溫度)來確定[12]。而用大氣溫度遞減率法LRT(lapse rate tropopause)來確定對流層頂高度，根據(jù)WMO于1957年給出的標(biāo)準(zhǔn)，如果遞減率小于2 K/km，且在此高度向上2 km的范圍內(nèi)溫度遞減率平均不超過2 K/km時，則認(rèn)為該范圍的最低高度為對流層頂高度[13]。大氣溫度遞減率的定義為：

(1)

其中，T為溫度，h為高度，按照由低到高的方法計算遞減率。

當(dāng)使用探空資料作為對照組時，探空氣球使用CPT法計算對流層頂高度，F(xiàn)Y-3C GNOS資料分別使用CPT法與LRT法計算對流層頂高度，二者的對比結(jié)果如圖2所示。

圖2 不同方法計算的結(jié)果對比Fig.2 Comparison of calculation resultsby different methods

由圖2可見，兩種方法計算得出的對流層頂高度與使用探空氣球計算得出的對流層頂高度差別不大，其中GNOSCPT-RadiosondeCPT的偏差基本在±2 km以內(nèi)，GNOSLRT-RadiosondeCPT的偏差均在±2.5 km以內(nèi)。將上述數(shù)據(jù)進(jìn)行統(tǒng)計，結(jié)果見表1。

由表1可知，使用兩種方法計算得出的對流層頂高度與使用探空氣球數(shù)據(jù)計算結(jié)果的平均偏差均不超過1 km，均方根誤差不超過1.3 km?？梢钥闯?，使用FY-3C掩星數(shù)據(jù)計算對流層頂高度具備可行性。

表1 兩種計算方法的精度對比Tab.1 Comparison of precision of two calculation methods

當(dāng)使用ERA5資料作為對照組時，ERA5數(shù)據(jù)使用CPT方法計算對流層頂高度，F(xiàn)Y-3C GNOS資料分別使用CPT法與LRT法計算對流層頂高度，二者的對比結(jié)果如圖3所示。

圖3 不同方法計算結(jié)果對比Fig.3 Comparison of calculation results by different methods

由圖3可見，兩種方法計算得出的對流層頂高度與使用ERA5數(shù)據(jù)計算得出的對流層頂高度差別不大，其中GNOSCPT-ERA5CPT的偏差均在±3 km以內(nèi)，GNOSLRT-ERA5CPT的偏差均在±2 km以內(nèi)。將上述數(shù)據(jù)進(jìn)行統(tǒng)計，結(jié)果見表2。

表2 兩種計算方法的精度對比Tab.2 Comparison of precision of two calculation methods

由表2可知，使用兩種方法計算得出的對流層頂高度與使用ERA5數(shù)據(jù)計算的結(jié)果平均偏差均不超過0.7 km，均方根誤差不超過1.4 km?？梢钥闯?，使用FY-3C掩星數(shù)據(jù)計算對流層頂高度具備可行性。

在驗證FY-3C掩星數(shù)據(jù)精度后，下文的對流層頂高度計算使用了CPT方法。

2 結(jié)果與分析

2.1 對流層頂參數(shù)空間特征變化

圖4展示了掩星廓線數(shù)量隨緯度的變化情況，可以看出，在32°～34°N上數(shù)據(jù)量最少，為251條；在48°～50°N上數(shù)據(jù)量最多，為471條。

圖4 掩星事件數(shù)目隨緯度的變化特征Fig.4 Variation of occultation events with latitude

圖5展示了掩星事件數(shù)目隨經(jīng)緯度的變化情況，可以看出，在16～24°N和42°～56°N上的掩星廓線數(shù)目較多，在24°～42°N上的每個柵格里的掩星廓線數(shù)目較少。綜合圖4、圖5可以看出，F(xiàn)Y-3C掩星數(shù)據(jù)基本覆蓋了整個中國區(qū)域。

圖5 掩星事件數(shù)目分布Fig.5 Distribution of occultation events

圖6展示了研究時段內(nèi)對流層頂高度、溫度隨經(jīng)緯度的變化，可以看出，中國區(qū)域的對流層頂高度以及溫度具有明顯的緯度分布特征，即隨緯度變化較為明顯，隨經(jīng)度的變化較小。由圖6(a)可見，中國南部區(qū)域?qū)α鲗禹敻叨雀哂谥袊辈繀^(qū)域，中國西部區(qū)域?qū)α鲗禹敻叨雀哂谥袊鴸|部區(qū)域。在整個16°～54°N范圍內(nèi)，可以明顯看到對流層頂高度隨緯度分帶分布的現(xiàn)象。由圖6(b)可見，中國南部區(qū)域?qū)α鲗禹敎囟鹊陀诒辈繀^(qū)域，中國西部區(qū)域?qū)α鲗禹敎囟嚷缘陀跂|部區(qū)域，在整個16°～54°N范圍內(nèi)，也可以看到明顯的條帶分布現(xiàn)象。

圖6 對流層頂高度、溫度隨經(jīng)緯度的變化Fig.6 Variation of tropopause height, temperature with latitude and longitude

從上文可知，對流層頂參數(shù)隨緯度變化較為明顯。圖7展示了對流層頂高度、溫度隨緯度的變化。總體來看，隨著緯度的升高，對流層頂高度逐漸減小，對流層頂溫度逐漸增大，但是到42°N開始出現(xiàn)略微下降。由圖7(a)可見，在研究時段內(nèi)，對流層頂高度隨著緯度的升高總體上呈現(xiàn)出遞減的趨勢。在22°～24°N存在最大值17.29 km，在52°～54°N存在最小值10.45 km，16°～28°N的低緯度處減小速率較緩，在28°～40°N處的中緯度地區(qū)對流層頂高度驟減且減小速率較大，40°～56°N處減小速率開始趨緩。由圖7(b)可見，對流層頂溫度隨緯度的升高而呈現(xiàn)升高的趨勢，在16°～18°N處有最小值-81.62 ℃，在42°～44°N處有最大值-55.84 ℃，在16°～22°N區(qū)間內(nèi)增長速度較緩，在22°～40°N區(qū)間內(nèi)增長較為迅速，在40°～56°N區(qū)間內(nèi)增長放緩。

圖7 對流層頂高度、溫度隨緯度的變化Fig.7 Variation of tropopause height and temperature with latitude

2.2 對流層頂參數(shù)的季節(jié)性變化

上文研究表明，在空間特征上，中國區(qū)域內(nèi)對流層頂參數(shù)隨緯度變化較為明顯，在經(jīng)度上變化并不明顯。所以可以在保持緯度值不變的情況下，通過觀察對流層頂參數(shù)隨時間的變化，來研究對流層頂隨時間的變化特征。表3展示了25°N、35°N、45°N對流層頂高度與溫度的最大值與最小值以及最大值與最小值出現(xiàn)的月份差，由表可知，對流層頂高度最大值到最小值的轉(zhuǎn)變，二者相差了5～6個月；對流層頂溫度最大值到最小值的轉(zhuǎn)變，二者相差了7～8個月。因此可以確定，對流層頂參數(shù)具有較為明顯的隨時間變化的特征。

表3 對流層頂高度的最大值與最小值Tab.3 Maximum and minimum of tropopause height

表4 對流層頂溫度最大值與最小值Tab.4 Maximum and minimum of tropopause temperature

圖8(a)～(d)分別為研究時段內(nèi)春、夏、秋、冬季所對應(yīng)的對流層頂高度，可以看出，隨著春-夏-秋-冬季節(jié)的輪替，對流層頂高度出現(xiàn)升高-降低-降低的趨勢，這一結(jié)果與單映穎[14]等人對中國區(qū)域?qū)α鲗禹數(shù)难芯拷Y(jié)果一致。在4個季節(jié)時間段內(nèi)，均存在對流層頂高度急劇下降的緯度值，春季為28°N，夏季為34°N，秋季為35°N，冬季為30°N?？梢钥闯?，對流層頂高度最小值出現(xiàn)在冬季的東北地區(qū)。春季的最大值出現(xiàn)在46°～48°N、22°～76°E的柵格內(nèi)，最小值出現(xiàn)在22°～24°N、76°～78°E的柵格內(nèi)；夏季的最大值出現(xiàn)在50°～52°N、132°～134°E的柵格內(nèi)，最小值出現(xiàn)在16°～18°N、114°～116°E的柵格內(nèi)；秋季的最大值出現(xiàn)在52°～54°N、132°～134°E的柵格內(nèi)，最小值出現(xiàn)在22°～24°N、112°～114°E的柵格內(nèi)；冬季的最大值出現(xiàn)在48°～50°N、128°～130°E的柵格內(nèi)，最小值出現(xiàn)在24°～26°N、132°～134°E的柵格內(nèi)。

圖8 中國區(qū)域?qū)α鲗禹敻叨仍?個季節(jié)的分布Fig.8 Distribution of tropopause height over China in spring,summer, autumn and winter

圖9(a)～(d)分別為研究時段內(nèi)春、夏、秋、冬季所對應(yīng)的對流層頂溫度。在每一個分圖上都可以明顯看到3個條帶。每一個季節(jié)里，北方區(qū)域的對流層頂溫度都要高于南方。伴隨著春-夏-秋-冬季節(jié)的輪替，對流層頂溫度出現(xiàn)升高-降低-降低的趨勢。春季的最大溫度值-47.68 ℃出現(xiàn)在46°～48°N、98°～100°E的柵格內(nèi)，最小值-86.76 ℃出現(xiàn)在16°～18°N、92°～94°E的柵格內(nèi)；夏季的最大溫度值-42.81 ℃出現(xiàn)在40°～42°N、72°～74°E的柵格內(nèi)，最小值-86.12 ℃出現(xiàn)在16°～18°N、82°～84°E的柵格內(nèi)；秋季的最大溫度值-47.68 ℃出現(xiàn)在48°～50°N、74°～76°E的柵格內(nèi)，最小值-86.17 ℃出現(xiàn)在20°～22°N、110°～112°E的柵格內(nèi)；冬季的最大溫度值-49.19 ℃出現(xiàn)在42°～44°N、132°～134°E的柵格內(nèi)，最小值-87.68 ℃出現(xiàn)在18°～20°N、94°～96°E的柵格內(nèi)。

圖9 中國區(qū)域?qū)α鲗禹敎囟仍?個季節(jié)的分布Fig.9 Distribution of tropopause temperature in spring, summer, autumn and winter in China

在4個季節(jié)時間段內(nèi)，春季與秋季的溫度隨緯度變化特征較為相似，且東部地區(qū)溫度要高于西部地區(qū)；在夏季 ，西北地區(qū)對流層頂溫度高于東北地區(qū)，西南地區(qū)對流層頂溫度低于東南地區(qū)；冬季正好與夏季相反，東北地區(qū)對流層頂溫度高于西北地區(qū)，西南地區(qū)對流層頂溫度高于東南地區(qū)。

3 結(jié) 語

本文首次使用FY-3C掩星數(shù)據(jù)對2017-03-01～2018-02-28中國區(qū)域的對流層頂結(jié)構(gòu)參數(shù)進(jìn)行研究。首先利用ROPP軟件處理相位延遲數(shù)據(jù)，獲得溫度廓線，估計對流層頂參數(shù)；然后分別使用探空氣球數(shù)據(jù)和ERA5數(shù)據(jù)對反演結(jié)果進(jìn)行驗證；最后根據(jù)CPT方法計估計了中國區(qū)域?qū)α鲗禹數(shù)臏囟群透叨葏?shù)。通過分析統(tǒng)計，結(jié)論如下：

1)FY-3C掩星資料反演得出的對流層頂參數(shù)與ERA5和探空氣球數(shù)據(jù)得出的結(jié)果比較一致，表明FY-3C掩星數(shù)據(jù)能夠估計對流層頂?shù)淖兓?/p>

2)中國區(qū)域?qū)α鲗禹敻叨纫约皽囟入S緯度變化呈現(xiàn)出明顯的分帶現(xiàn)象。其中對流層頂高度隨緯度的升高而降低，即北方對流層頂高度低于南方，且在中緯度地區(qū)出現(xiàn)了對流層頂高度隨緯度降低速率急劇加快的現(xiàn)象；對流層頂溫度隨緯度的升高而升高，即北方對流層頂溫度高于南方，且在中緯度地區(qū)出現(xiàn)了對流層頂溫度增加速率急劇加快的現(xiàn)象。

3)對流層頂高度與溫度具有隨時間變化的特征。在研究時間段內(nèi)，對流層頂高度參數(shù)從最大值到最小值或最小值到最大值的過渡過程中，均有較大的時間間隔。通過將研究時間段分為春、夏、秋、冬4個季節(jié)進(jìn)行統(tǒng)計發(fā)現(xiàn)，對流層頂高度隨春-夏-秋-冬季節(jié)變化表現(xiàn)出升高-降低-降低的趨勢；對流層頂溫度表現(xiàn)為春、秋季節(jié)變化相似，冬、夏季節(jié)變化相反的現(xiàn)象。