黨皓飛 毛 穎 張春桂 吳作航

(1.福建省氣象科學研究所，福建 福州 350008;2.福建省災害天氣重點實驗室，福建 福州 350008;3.福建省氣象災害防御技術中心，福建 福州 350008)

1 概述

暖區(qū)暴雨是指發(fā)生在地面鋒面南側(cè)暖區(qū)或華南地區(qū)無冷空氣影響時產(chǎn)生的暴雨[1-2]。暖區(qū)暴雨對流活動顯著，短時雨強大、突發(fā)性強，易引發(fā)災害，給群眾生命財產(chǎn)安全與經(jīng)濟建設帶來極大威脅[3-4]。由于暖區(qū)暴雨機理復雜且具有上述幾類特點，導致對暖區(qū)暴雨的預報與研究成為氣象工作中的一個重點及難點[5]。吳亞麗[6]基于再分析資料、雷達、地面降水資料進行了華南前汛期暖區(qū)暴雨的發(fā)展過程分析；肖紅茹等[7]將四川盆地暖區(qū)暴雨分為四類，并選取了典型個例，利用地面降水數(shù)據(jù)分析了暴雨成因與中尺度特征。在對暖區(qū)暴雨的研究中，及時了解其降水量的空間分布情況尤為重要[8]。傳統(tǒng)的地基雷達與地面雨量觀測站雖然是目前運用最普遍的觀測方法，但雷達的降水反演精度相對較低、觀測站分布不均、地形地勢影響等一些問題依舊存在[9]。

隨著遙感技術的發(fā)展，使得大范圍暖區(qū)暴雨降水觀測成為了可能[10]。氣象衛(wèi)星不受自然條件與地理條件限制，可實現(xiàn)高時空分辨率、大范圍的暖區(qū)暴雨監(jiān)測，與地面觀測數(shù)據(jù)實現(xiàn)互補[11]。目前常見的衛(wèi)星遙感反演降水的方法主要有可見光與紅外[12]、微波[13]、多傳感器組合[14]等。夏雙等[15]使用FY-2C紅外波段與人工神經(jīng)網(wǎng)絡方法進行了藏北高原典型區(qū)降水模擬，效果較優(yōu)；李萬彪等[16]利用阜陽天氣雷達雨量資料與微波成像資料，對TRMM/TMI的微波通道及其組合進行了動態(tài)聚類，構建了降水反演算法，并與日本NASDA降水產(chǎn)品進行了對比，證明了其反演方法的可靠性；張春桂等[17]結(jié)合風云二號、風云三號衛(wèi)星數(shù)據(jù)、天氣雷達與降水數(shù)據(jù)，運用逐步線性回歸方法進行了云團降水的模擬與預報。國內(nèi)外也通過衛(wèi)星數(shù)據(jù)制作了多種降水反演產(chǎn)品，如國外的TRMM、GPM、CMOPRH以及我國的FY衛(wèi)星系列降水反演產(chǎn)品[18]。隨著風云三號C星與風云三號D星(FY-3C/D)的升空并實現(xiàn)組網(wǎng)運行，大氣探測精度大大提高，也為暖區(qū)暴雨監(jiān)測提供了切實可靠的數(shù)據(jù)產(chǎn)品。

從目前相關研究文獻看，利用衛(wèi)星微波觀測資料開展暖區(qū)降水監(jiān)測的研究少有報道。本研究將利用我國自主發(fā)射的風云衛(wèi)星微波觀測資料，針對福建暖區(qū)暴雨的遙感監(jiān)測開展探索研究，以期得到相關初步結(jié)論，同時由于衛(wèi)星數(shù)據(jù)、降水反演過程與地理位置不同，不可避免與實際情況產(chǎn)生誤差[19]，故需要對FY-3C/D降水反演產(chǎn)品進行星地校驗評估，探尋FY-3C/D降水反演產(chǎn)品在福建暖區(qū)降水監(jiān)測中的實際應用效果。一些模式的模擬研究指出，不同模式對云宏、微觀物理特性的描述差異會直接導致模擬結(jié)果存在很大差異，從而使得云成為影響氣候變化預測準確性的重要因子。現(xiàn)有大量研究根據(jù)云的宏、微觀物理特性對云進行了分類。中國云圖分類標準和國際衛(wèi)星云氣候計劃(ISCCP)將云劃分為卷云、高層云、高積云等類型。研究結(jié)果表明，不同類型的云受其宏觀和微觀物理特性差異的影響，在輻射效應和降水類型等方面存在明顯差異。因此，研究不同云類型對于準確理解地氣系統(tǒng)輻射能量收支、降雨的分布，以及全球氣候變暖等問題都具有重要意義。

2 數(shù)據(jù)和研究方法

2.1 研究區(qū)域概況

福建省地處中國東南沿海，全省面積約12.4萬km2，包含9個設區(qū)市與1個綜合實驗區(qū)。省內(nèi)地勢從西北向東南降低，丘陵延綿，山地聳峙，有著“八山一水一分田”之稱。福建省東面沿海，海岸線綿長，為亞熱帶季風氣候，光照充足、雨量充沛，氣候濕潤，且處在熱帶東風帶與中緯度西風帶交替影響的過渡地帶，因此是中國暴雨高頻區(qū)之一。不僅如此，近些年通過觀測發(fā)現(xiàn)，暖區(qū)暴雨是福建省暴雨的重要組成部分，而暖區(qū)暴雨降水強度大、集中，無疑會造成滑坡、山洪、泥石流等重大災害，給人民的生命安全和財產(chǎn)安全產(chǎn)生巨大威脅。

2.2 實測樣本采集與預處理

實測樣本基于全國綜合氣象信息共享平臺(CIMISS)，并通過氣象數(shù)據(jù)統(tǒng)一服務接口(MUSIC)獲取。實測數(shù)據(jù)包括：2016—2019年間，福建省18場暖區(qū)暴雨過程中各測站點的名稱、站號、經(jīng)緯度、具體觀測時間，以及每小時觀測1次的過去1小時的累計雨量等。實測樣本使用前，將缺測以及不在合理范圍內(nèi)的降水數(shù)據(jù)予以剔除。

2.3 衛(wèi)星數(shù)據(jù)采集與預處理

風云三號(FY-3)氣象衛(wèi)星是我國第二代極軌氣象衛(wèi)星，相較于風云一號氣象衛(wèi)星，其功能與技術上皆有很大的提升。FY-3氣象衛(wèi)星不僅大幅提高了全球資料的獲取能力，同時也提升了地表和云區(qū)特征的遙感能力等。本研究采集的衛(wèi)星數(shù)據(jù)包括風云三號C星(FY-3C)大氣垂直探測系統(tǒng)(VASS)中大氣溫濕廓線(AVP)產(chǎn)品(分辨率星下點約17 km)、風云三號D星(FY-3D)MWTS/MWHS融合(TSHS)中的大氣溫濕廓線(AVP)產(chǎn)品(分辨率星下點約33 km)。具體提取的產(chǎn)品參數(shù)為：經(jīng)緯度以及微波濕度計中的亮溫資料。相關數(shù)據(jù)可從國家衛(wèi)星氣象中心下載(http://satellite.nsmc.org.cn/PortalSite/Data/Satellite.aspx)。注意到，微波濕度計中的亮溫資料共有15個通道，其中心頻率分別為89.0GHz、118.75±0.08GHz、118.75±0.2GHz、118.75±0.3GHz、118.75±0.8GHz、118.75±1.1GHz、118.75±2.5GHz、118.75±3.0GHz、118.75±5.0GHz、150.0GHz、183.31±1GHz、183.31±1.8GHz、183.31±3GHz、183.31±4.5GHz與183.31±7GHz，除89.0GHz與150.0GHz為V極化外，其余的波段皆為H極化?；谖⒉穸扔嫸嗤ǖ懒翜匦畔⑦b感反演，可在一定程度上提高暴雨、臺風等災害性強對流天氣的監(jiān)測和預警能力。本研究中，衛(wèi)星數(shù)據(jù)應用前，將經(jīng)緯度異常以及亮溫不在合理范圍內(nèi)的數(shù)據(jù)剔除。

2.4 衛(wèi)星數(shù)據(jù)與實測樣本的匹配

由于在數(shù)據(jù)的時空采樣方面，地面實測降水資料與衛(wèi)星遙感影像間存在一定差異，且為了探索FY-3C/D衛(wèi)星對暖區(qū)暴雨提前1小時的預報預測可行性，本研究在對FY-3C VASS與FY-3D TSHS亮溫數(shù)據(jù)與實測降水樣本間開展的星地匹配過程中進行了時空上的約束。具體做法為：衛(wèi)星影像的獲取時間比降水實測數(shù)據(jù)提早1個小時(如：過去1小時累計降水實測數(shù)據(jù)記錄的時間為2016年4月11日15時，則衛(wèi)星觀測的時間需為2016年4月11日14時)。為避免星地匹配數(shù)據(jù)點存在較大的地理位置差異，因此衛(wèi)星數(shù)據(jù)與地面實測降水站點的經(jīng)緯度差的平方和需滿足≤0.002的條件?；谝陨戏椒ǎ罱K共獲取1851組星地匹配數(shù)據(jù)。

2.5 降水反演算法

本研究采用多元線性回歸方法進行FY-3C/D的降水遙感反演。多元線性回歸是用線性擬合的方法把多個自變量x1，x2，...，xi和因變量y之間存在的線性關系表達出來，這種線性關系可表示為：

y=b0+b1x1+b2x2+...+bixi

(1)

式(1)中，b0為截距，b1、b2、...bi為回歸系數(shù)。本研究中，多元線性回歸算法反演暖區(qū)暴雨降水時，自變量設定為15個，分別對應著FY-3C VASS與FY-3D TSHS中微波濕度計15個通道的亮溫數(shù)據(jù)減去一個基數(shù)(本研究設為160)以后的數(shù)值，因變量則為實測降水數(shù)據(jù)。最后通過公式(1)擬合亮溫數(shù)據(jù)與實測降水數(shù)據(jù)間的關系，從而確定公式(1)中的回歸系數(shù)，用作之后的暖區(qū)暴雨降水值反演。在反演算法建立之前，本研究從1851組星地匹配數(shù)據(jù)抽取70%的數(shù)據(jù)(1296組)進行FY-3C/D遙感衛(wèi)星建模，剩余30%的數(shù)據(jù)(555組)開展所建立模型的精度驗證，以確保建模與驗證數(shù)據(jù)集中降水資料的數(shù)值范圍是均勻分布的。

2.6 模型精度評價

本研究采用的模型精度評價指標包括相關系數(shù)R，均方根誤差RMSE，以及平均相對誤差的絕對值MAPE。具體計算公式如下：

(2)

(3)

其中，N為樣本數(shù)，i為第i個樣本，x與y分別代表降雨實測值與模型反演值。

3 結(jié)果與分析

3.1 基于FY-3C/D遙感衛(wèi)星資料的暖區(qū)降水模型

通過多次試驗發(fā)現(xiàn)，將建模的降水數(shù)據(jù)按照一定的數(shù)值范圍進行分組，并分別對不同組別的數(shù)據(jù)進行多元線性回歸擬合，即依據(jù)降水數(shù)值的大小進行分段建模時，可以得到較好的衛(wèi)星反演結(jié)果。本研究中，建模降水數(shù)據(jù)從小到大分成3組，分別是0mm<降雨量≤10mm，10mm<降雨量≤25mm以及降雨量>25mm，得到的分段模型為：當0mm<降雨量≤10mm時，多元線性回歸擬合結(jié)果如公式(4)所示；當10mm<降雨量≤25mm時，多元線性回歸擬合結(jié)果如公式(5)所示；當降雨量>25mm時，多元線性回歸擬合結(jié)果如公式(6)所示。

y=2.9463+0.0527x1+0.083x2+0.1589x3-0.2332x4+0.1161x5-0.0123x6-0.1939x7+0.1601x8-0.1568x9-0.0636x10+0.0939x11-0.0208x12-0.0925x13-0.0212x14+0.2047x15

(4)

y=32.818-0.1451x1-0.5968x2+0.4669x3-0.0069x4-1.5924x5+1.8999x6-0.9451x7-1.0072x8+1.4428x9+0.024x10+1.4318x11-2.13x12+1.4182x13-0.5374x14+0.0029x15

(5)

y=2.0279x1-1.0359x2+2.2855x3+0.6423x5-3.3112x9-0.6189x10+0.9371x11-4.6602x14+5.3032x15

(6)

為定量計算所建立衛(wèi)星模型的反演精度，本研究基于2.6節(jié)計算出相應的模型精度評價指標，如表1所示。為了更加直觀地表示模型反演質(zhì)量，給出了如圖1所示的基于建模數(shù)據(jù)的算法反演結(jié)果與實測降水數(shù)據(jù)對比散點圖。從表1可知，降雨量>25mm這一量級內(nèi)，星地匹配的數(shù)據(jù)量最少，但反演精度最好；10mm<降水量≤25mm這一量級內(nèi)，星地匹配數(shù)據(jù)量與反演精度次之；而0mm<降水量≤10mm這一量級內(nèi)，星地匹配數(shù)據(jù)量最多，反演精度卻最差，MAPE值甚至高達330.06%。這說明在FY-3C/D衛(wèi)星建模中，其對于較小量級降水的反演結(jié)果較差，而對于較大量級降水的反演結(jié)果較好。從圖1也可看出，0mm<降水量≤10mm的反演結(jié)果與實測數(shù)據(jù)對比散點圖中出現(xiàn)了明顯的離散現(xiàn)象，而降雨量>25mm的反演結(jié)果與實測數(shù)據(jù)散點則較好地分布在1∶1線附近，這進一步表明在建模數(shù)據(jù)中，F(xiàn)Y-3C/D對于較大量級的降水遙感反演要優(yōu)于較小量級的降水遙感反演。

表1 基于建模數(shù)據(jù)的多元線性回歸擬合降水精度

(a)實測降水量級為：0mm<降雨量≤10mm (b)實測降水量級為：10mm<降雨量≤25mm

3.2 暖區(qū)降水模型的精度驗證

將除去建模外的剩余數(shù)據(jù)(555組)對反演算法進行驗證。表2為基于誤差驗證數(shù)據(jù)的多元線性回歸擬合降水精度結(jié)果，可見驗證數(shù)據(jù)的反演精度結(jié)果較差于建模數(shù)據(jù)的反演結(jié)果，但也出現(xiàn)了較大量級降水星地匹配組數(shù)小于較小量級降水，而較大量級降水反演精度質(zhì)量要優(yōu)于較小量級降水的情況。圖2 為基于驗證數(shù)據(jù)的算法反演結(jié)果與實測降水數(shù)據(jù)對比散點圖，可見0mm<降水量≤10mm量級的反演結(jié)果出現(xiàn)的離散現(xiàn)象比10mm<降水量≤25mm量級反演結(jié)果的離散情況明顯，而降雨量>25mm量級出現(xiàn)了一定的低估。

表2 基于誤差驗證數(shù)據(jù)的多元線性回歸擬合降水精度

(a)實測降水量級為：0mm<降雨量≤10mm (b)實測降水量級為：10mm<降雨量≤25mm

4 結(jié)論與展望

在本次模型的建立與精度驗證中，F(xiàn)Y-3C/D衛(wèi)星對于較小量級降水(0mm<降水量≤10mm)的反演結(jié)果較差，離散現(xiàn)象明顯。對于較大量級降水(降雨量>25mm)，盡管在驗證結(jié)果中出現(xiàn)了一定的低估現(xiàn)象，但總體反演結(jié)果較好。因此，基于FY-3C/D微波觀測資料的福建暖區(qū)降水遙感監(jiān)測模型對于較大量級的降水反演具有一定的指示意義。由于本文中采用的是多元線性回歸模型，因此模型的精度可能隨著時間與地域的變化發(fā)生較大改變。在未來研究中建議收集更多數(shù)據(jù)，并嘗試以物理模型等方法解決這一問題。