張蔚然 劉黎平 吳翀

1 中國(guó)氣象科學(xué)研究院災(zāi)害天氣重點(diǎn)實(shí)驗(yàn)室, 北京 100081

2 陜西省氣象臺(tái), 西安 710014

1 引言

青藏高原地處我國(guó)西南部，是世界上海拔最高、范圍最大、地形最復(fù)雜的高原。獨(dú)特的地理環(huán)境趨勢(shì)使該地區(qū)產(chǎn)生了特殊的熱力和動(dòng)力作用，使得高原對(duì)流活動(dòng)頻發(fā)，對(duì)亞洲甚至全世界的氣候變化都有著重要的影響（Yanai et al., 1992; Webster et al.,1998; Yu et al., 2004），因此對(duì)青藏高原地區(qū)的相關(guān)研究十分必要。魏麗和鐘強(qiáng)（1997）利用國(guó)際衛(wèi)星云氣候計(jì)劃（ISCCP）資料和全球地面觀測(cè)云氣候資料，對(duì)青藏高原地區(qū)（20°～50°N，70°～110°E）云的水平與垂直分布進(jìn)行研究，發(fā)現(xiàn)夏季高原云量明顯多于冬季，且云狀存在明顯的地區(qū)性分布特征，高原主體以積雨云為主，北部高云出現(xiàn)較多，東南部層狀低云出現(xiàn)較多。

青藏高原大地形產(chǎn)生的特殊熱力和動(dòng)力作用，使該地可能會(huì)產(chǎn)生一些獨(dú)特的區(qū)域日降水特征（Kuo and Qian, 1981; Fujinami et al., 2005）。此前一些研究對(duì)青藏高原地區(qū)降水日變化特征進(jìn)行了分析。Yu et al.（2007）通過(guò)分析觀測(cè)臺(tái)站逐時(shí)降水資料，發(fā)現(xiàn)在青藏高原大部分地區(qū)降水高峰出現(xiàn)在午夜前后。通過(guò)對(duì)在高原主體那曲地區(qū)GAMETibet（GEWEX Asian Monsoon Experiment）期 間獲取的多普勒雷達(dá)資料、地面測(cè)站資料等的分析，Liu et al.（2002）發(fā)現(xiàn)那曲降水量在傍晚達(dá)到最大值，清晨降水量最小。Ueno et al.（2001）同樣表明那曲的小時(shí)降水強(qiáng)度在夜間大于白天。通過(guò)對(duì)1998～2007 年夏季風(fēng)季節(jié)（6～8 月）TRMM 降水雷達(dá)數(shù)據(jù)資料分析，Singh and Nakamura（2009）研究了高原中部降水特征，指出丘陵地區(qū)降水在傍晚最強(qiáng)，山谷和湖泊降水峰值則出現(xiàn)在夜間。對(duì)于青藏高原東緣地區(qū)，Wang et al.（2004）利用地球靜止氣象衛(wèi)星（GMS）1998～2001 年5～8 月的每小時(shí)紅外亮溫?cái)?shù)據(jù)了解到該地對(duì)流日變化在午后或傍晚達(dá)到峰值。Zhou et al.（2008）利用2000～2004 年雨量計(jì)數(shù)據(jù)也報(bào)告了青藏高原東部的降水有近午夜最大值。對(duì)于青藏高原南部，Barros et al.（2004）通過(guò)分析Meteosat-5 的紅外衛(wèi)星數(shù)據(jù)，發(fā)現(xiàn)夏季風(fēng)期間喜馬拉雅山南坡降水日循環(huán)明顯，當(dāng)?shù)貢r(shí)夜晚00:00～03:00 達(dá)到峰值（除特別說(shuō)明為當(dāng)?shù)貢r(shí)外，均為北京時(shí)，當(dāng)?shù)貢r(shí)=北京時(shí)-2 小時(shí)）。Bhatt and Nakamura（2005）利用TRMM 降水雷達(dá)資料研究發(fā)現(xiàn)在喜馬拉雅山南坡，季風(fēng)爆發(fā)前降水主要發(fā)生在下午，季風(fēng)爆發(fā)后午夜至清晨降水量最大。Chen et al.（2012）通過(guò)分析暖季青藏高原東南部色赤拉山山坡觀測(cè)站與常規(guī)山谷中觀測(cè)站的資料，探討了測(cè)站位置對(duì)臺(tái)站記錄與TRMM 衛(wèi)星數(shù)據(jù)間不同日降水量變化的影響。結(jié)果表明降水在山坡站表現(xiàn)為午后峰值與TRMM 數(shù)據(jù)所示相似，在山谷站則盛行夜間降水。Fu et al.（2018）通過(guò)分析TRMM 衛(wèi)星資料指出對(duì)于高原南坡區(qū)域，溫暖濕潤(rùn)的夏季風(fēng)帶來(lái)的充足水汽使該地區(qū)形成了獨(dú)特的云和降水。尤其在高原東南坡，夏季基本被云系所覆蓋（譚瑞婷等, 2018）。Chen et al.（2017）通過(guò)分析CloudSat 衛(wèi)星資料揭示了高原南坡不同地形下云團(tuán)的基本特征，指出云團(tuán)的性質(zhì)差異可能是強(qiáng)上升運(yùn)動(dòng)與喜馬拉雅地形相互作用的結(jié)果。Yu et al.（2020）利用CloudSat/CALIPSO 資料分析了高原喇叭口地區(qū)各種云類、云層、云水路徑以及垂直結(jié)構(gòu)的日夜差異特征并討論了造成這種差異的物理機(jī)制。另外，此前研究（Shimizu et al., 2001;Kurosaki and Kimura, 2002）表明青藏高原在亞洲季風(fēng)爆發(fā)后，高原上的云頂高度會(huì)顯著增高，云降水特征發(fā)生明顯變化。Uyeda et al.（2001）與馮錦明等（2002）使用GAME-Tibet 期間獲取的資料，發(fā)現(xiàn)在高原雨季來(lái)臨之后，那曲地區(qū)對(duì)流活動(dòng)明顯增強(qiáng)，雷達(dá)回波強(qiáng)度增大，回波頂高和回波面積增加。對(duì)流活動(dòng)在白天發(fā)展，對(duì)流強(qiáng)度在當(dāng)?shù)貢r(shí)下午14:00 達(dá)到最強(qiáng)，夜晚到次日凌晨逐漸消散。

上述研究多借助高原臺(tái)站和衛(wèi)星遙感資料進(jìn)行。但青藏高原地處偏遠(yuǎn)、地形復(fù)雜，觀測(cè)臺(tái)站極為稀疏，臺(tái)站資料代表性差。衛(wèi)星觀測(cè)在近地面處易受雜波干擾，且衛(wèi)星的時(shí)空分辨率相對(duì)較低。另外，此前僅在青藏高原主體那曲地區(qū)進(jìn)行過(guò)多普勒雷達(dá)觀測(cè)，目前對(duì)于墨脫云降水特征的研究基本是空白的。墨脫以高山峽谷地形為主，呈亞熱帶濕潤(rùn)氣候，其復(fù)雜的地形和特殊的氣候使該地云降水特征與高原主體那曲地區(qū)可能有顯著的差別?？紤]到墨脫地區(qū)交通不便，電力供應(yīng)差，地形遮擋等條件，第二次青藏高原綜合科學(xué)考察研究項(xiàng)目在墨脫設(shè)置了一部X 波段相控陣偏振雷達(dá)，在國(guó)內(nèi)首次采用最先進(jìn)的雙偏振相控陣?yán)走_(dá)，實(shí)現(xiàn)了對(duì)河谷地區(qū)云降水的雷達(dá)連續(xù)觀測(cè)。采用X 波段相控陣天氣雷達(dá)進(jìn)行觀測(cè)的目的主要是提高雷達(dá)的穩(wěn)定性和可靠性，實(shí)現(xiàn)無(wú)人值守的連續(xù)觀測(cè)。目前，我國(guó)已經(jīng)在大灣區(qū)建設(shè)了28 部X 波段相控陣天氣雷達(dá)。為了分析該雷達(dá)觀測(cè)數(shù)據(jù)的一致性，張蔚然等（2021）利用廣州市業(yè)務(wù)運(yùn)行的S 波段雙線偏振雷達(dá)作為“標(biāo)準(zhǔn)”，分析了X 波段雙偏振相控陣天氣雷達(dá)回波強(qiáng)度、差分反射率和差分相位移率的系統(tǒng)偏差及隨仰角的變化，初步驗(yàn)證了該型號(hào)雷達(dá)的數(shù)據(jù)可靠性。本文使用墨脫X 波段雙偏振相控陣?yán)走_(dá)獲取的2019 年11 月至2020 年10 月一整年的觀測(cè)數(shù)據(jù)資料，定量分析了墨脫地區(qū)云降水的回波強(qiáng)度、回波頂高等參數(shù)的月變化、日變化和高度變化，并與青藏高原主體那曲地區(qū)夏季季風(fēng)時(shí)期的觀測(cè)結(jié)果進(jìn)行對(duì)比，這為認(rèn)識(shí)河灣地區(qū)水汽輸送結(jié)構(gòu)特征，分析季風(fēng)、地形、山谷風(fēng)等對(duì)該地區(qū)降水的影響提供了可能，對(duì)亞洲水塔水資源變化評(píng)估具有重要的科學(xué)與應(yīng)用價(jià)值。

2 數(shù)據(jù)與方法

2.1 觀測(cè)站點(diǎn)、設(shè)備和數(shù)據(jù)說(shuō)明

青藏高原東南部是青藏高原水汽輸送過(guò)程的關(guān)鍵區(qū)，來(lái)自印度洋和孟加拉灣的大量水汽通過(guò)高原東南部的雅魯藏布大峽谷輸送到高原內(nèi)部（高登義等, 1985）。墨脫縣位于高原東南部，喜馬拉雅山南坡，該地平均海拔高度為1200 m，為典型河谷地形。該地區(qū)三面環(huán)山，地勢(shì)北高南低，海拔跨度范圍大，屬于亞熱帶濕潤(rùn)氣候，年平均氣溫為18°C，年平均降水量為2093 mm（陳萍和李波,2018）。強(qiáng)大的暖濕氣流使墨脫地區(qū)常年云霧繚繞、溫暖濕潤(rùn)，形成了一個(gè)較強(qiáng)的降水帶。圖1 為利用2019 年11 月至2020 年10 月ERA-5 整層水汽通量數(shù)據(jù)計(jì)算的平均水汽通量疊加在地形上的示意圖。

圖1 青藏高原地形（陰影，單位：m）與墨脫站地理位置（29°18′46″N，95°19′03″E，海拔高度為1305 m）以及2019 年11月至2020 年10 月平均整層水汽通量（箭頭，單位：kg m-1 s-1）Fig.1 Topography of Qinghai-Tibet Plateau (shaded, units: m),geographical location of Motuo (29°18′46″N，95°19′03″E, 1305 m ASL), and average whole layer water vapor flux from November 2019 to October 2020 (arrow, units: kg m-1 s-1)

為了了解高原東南部河谷地區(qū)的云降水特征，第二次青藏高原綜合科學(xué)考察研究項(xiàng)目，在墨脫站（29°18′46″N，95°19′03″E，海拔高度為1305 m）設(shè)置了1 部X 波段雙偏振相控陣?yán)走_(dá)（下文簡(jiǎn)稱X-PAR），該雷達(dá)由珠海納睿達(dá)公司研制，它采用全相參脈沖多普勒、雙線偏振和相控陣體制，主要技術(shù)參數(shù)見表1。該雷達(dá)的峰值功率為256 W，采用間隔1.8°共12 層的距離高度掃描（RHI），形成周期為92 s 的體掃數(shù)據(jù)，可獲取42 km 范圍內(nèi)的回波強(qiáng)度（ZH）、差分反射率因子（ZDR）、差分傳播相移率（KDP）、相關(guān)系數(shù)（ρhv）等數(shù)據(jù)，其距離分辨率為30 m。

表1 X-PAR 雷達(dá)參數(shù)Table 1 The parameters of X-PAR

已知X 波段電磁波的衰減問(wèn)題比起S 波段與C 波段較為嚴(yán)重。散射模擬計(jì)算指出：X 波段電磁波的單程衰減率（AH）分別是C、S 波段的7～8倍和10 倍以上（Park et al., 2005），因此衰減訂正工作十分必要。雷達(dá)反射率衰減訂正的本質(zhì)是估算AH，如公式（1）所示：

其中，ZHa為訂正前反射率因子值、ZHe為訂正前的反射率因子值，r為雷達(dá)探測(cè)目標(biāo)與雷達(dá)中心之間的距離。

對(duì)于X-PAR 的ZH觀測(cè)數(shù)據(jù)，本文使用Testud et al.（2000）提出“ZPHI 降水廓線訂正算法”進(jìn)行衰減訂正。ZPHI 訂正法存在一個(gè)約束條件：路徑總衰減（PIA）與整個(gè)分段區(qū)間（r0，r1）里的ΦDP的變化ΔΦDP是一致的。ZPHI 算法要進(jìn)行雨區(qū)分段AH的計(jì)算就是基于整個(gè)約束條件。通過(guò)設(shè)定n+1 個(gè)邊界將1 個(gè)徑向分成n段來(lái)進(jìn)行訂正。在每分段中，AH的計(jì)算如下：

式中：

其中，b為常數(shù)，范圍為0.757～0.804，γ在X 波段取0.32，ΔΦ為與區(qū)范圍內(nèi)的差分傳播相移的變化總量。

圖2 為X-PAR 可探測(cè)范圍內(nèi)的墨脫地形圖（每個(gè)圈為14 km），從圖中可以看出雷達(dá)處于峽谷之中。雷達(dá)周圍地形極其復(fù)雜，海拔高度跨度大，其北部為崗日嘎布山脈，西部與西北部屬于喜馬拉雅高山地段，為南迦巴瓦峰，其東部為米什米山脈。墨脫地處喜馬拉雅山東段南坡（西南季風(fēng)迎風(fēng)坡）以及外寬內(nèi)窄的山谷地形使該地降水類型主要以地形雨為主。暖濕夏季風(fēng)不斷到喜馬拉雅山南坡的陡峭地帶，有利于形成獨(dú)特的云系。本文使用墨脫XPAR 在2019 年11 月至2020 年10 月獲取的觀測(cè)基數(shù)據(jù)資料對(duì)該地區(qū)的云降水宏觀特征進(jìn)行統(tǒng)計(jì)分析，需指出的是本文進(jìn)行的是所有降水過(guò)程的統(tǒng)計(jì)分析，并未進(jìn)行天氣類型的分類分析。

圖2 觀測(cè)區(qū)地形示意圖Fig.2 Terrain of the observation region

2.2 數(shù)據(jù)處理方法

2.2.1 地物回波與雜散點(diǎn)回波的消除

由于墨脫地區(qū)地形復(fù)雜，雷達(dá)站周圍受到山脈的阻擋，雷達(dá)數(shù)據(jù)的地物回波勢(shì)必會(huì)影響墨脫地區(qū)云降水特征的統(tǒng)計(jì)結(jié)果，雖然X-PAR 使用高斯自適應(yīng)濾波器方式進(jìn)行了地物回波濾除，但數(shù)據(jù)中仍有地物回波存在。通過(guò)觀察長(zhǎng)時(shí)間的墨脫站XPAR 的數(shù)據(jù)平面顯示圖（下文簡(jiǎn)稱PPI）發(fā)現(xiàn)雷達(dá)參量在每層仰角上仍受到一些地物與雜散點(diǎn)回波的影響，其主要位于X-PAR 的東南方向與西北方向（圖3a-d）。在上述兩方向區(qū)域，地物雜散點(diǎn)回波的ZH與降水回波ZH近乎無(wú)變化（圖3a），但雙偏振參量有較為明顯的變化，其中，ρhv小于0.85（圖3c），ZDR小于-1 dB（圖3b），KDP大于5° km-1（圖3d）。根據(jù)上述特征，剔除地物回波。另外，為了除去X-PAR 回波中的雜散點(diǎn)，選擇某點(diǎn)數(shù)據(jù)與相鄰徑向、距離庫(kù)的數(shù)據(jù)組成的3×3 數(shù)據(jù)格點(diǎn)，若是有數(shù)據(jù)格點(diǎn)比例超過(guò)總格點(diǎn)的70%，則將其視為降水回波數(shù)據(jù)，否則為雜散點(diǎn)回波。圖3e-f 顯示利用上述特征有效地去除了X-PAR 紅圈內(nèi)地物與雜散點(diǎn)回波。

圖3 2020 年10 月4 日04:53:46 X-PAR 觀測(cè)第7 層仰角（10.7°）的（a）ZH，（b）ZDR，（c）ρhv，（d）KDP，除掉地物雜點(diǎn)回波的（e）ZH，（f）ZDR，（g）ρhv，（h）KDP。紅圈位置為地物和雜點(diǎn)Fig.3 Plan position indicator (PPI) of (a) ZH, (b) ZDR, (c) ρhv, and (d) KDP at the elevation of the 7th layer (10.7°) observed by X-PAR on October 4,2020, at 0453:46 BJT; (e) ZH, (f) ZDR, (g) ρhv, and (h) KDP are those whose ground clutter and spurious echo have been removed.The position of the red circle is the echo of ground clutter and spurious

2.2.2 雷達(dá)有效探測(cè)回波

墨脫地區(qū)地形復(fù)雜，海拔跨度大，在不同位置地形對(duì)雷達(dá)回波的遮擋不同。本文利用X-PAR 掃描模式，將雷達(dá)體掃球坐標(biāo)轉(zhuǎn)換為笛卡爾直角坐標(biāo)格點(diǎn)，結(jié)合地形高度圖給出了X-PAR 不同海拔高度的探測(cè)范圍（圖4），并統(tǒng)計(jì)了雷達(dá)在不同海拔高度可探測(cè)格點(diǎn)數(shù)的變化情況（圖4，格點(diǎn)大小為30m×30m）。圖5 顯示X-PAR 的可探測(cè)回波點(diǎn)數(shù)隨海拔高度升高先增加后減小，在2 km 以下高度，雷達(dá)探測(cè)到的回波格點(diǎn)十分有限，在5 km 左右雷達(dá)可探測(cè)回波格點(diǎn)數(shù)達(dá)到最大，隨后可探測(cè)點(diǎn)數(shù)隨高度升高而減少。上述現(xiàn)象是因?yàn)槔走_(dá)有效探測(cè)范圍除了受地形遮擋影響之外，還受雷達(dá)頂部的圓形盲區(qū)影響。在探測(cè)高度較低時(shí)，雷達(dá)探測(cè)主要受地形遮擋影響，隨著海拔高度的升高，雷達(dá)頂部盲區(qū)的影響越來(lái)越大。

圖4 部分高度層X-PAR 有效探測(cè)回波范圍（黑色為可探測(cè)到的區(qū)域）：（a）2 km；（b）3 km；（c）4 km；（d）5 km；（e）9 km；（f）10 km；（g）11 km；（h）12 kmFig.4 Effective detection echo range of X-PAR is in some altitude layers (black area represents the detectable area): (a) 2 km; (b) 3 km; (c) 4 km;(d) 5 km; (e) 9 km; (f) 10 km; (g) 11 km; (h) 12 km

通過(guò)圖5 可確定雷達(dá)觀測(cè)時(shí)間段內(nèi)、不同海拔高度理論可探測(cè)格點(diǎn)數(shù)，結(jié)合雷達(dá)實(shí)際觀測(cè)到的回波點(diǎn)數(shù)，就可以統(tǒng)計(jì)不同月份、不同小時(shí)、不同高度雷達(dá)探測(cè)回波的實(shí)際分布情況。

圖5 墨脫X-PAR 海拔高度可探測(cè)點(diǎn)數(shù)廓線圖（高度分辨率：100 m）Fig.5 Profile of detectable points of Motuo X-PAR with altitude(Height resolution: 100 m)

2.2.3 雷達(dá)觀測(cè)數(shù)據(jù)統(tǒng)計(jì)方法

為了探究西藏墨脫云降水特征，需要對(duì)XPAR 觀測(cè)到的數(shù)據(jù)進(jìn)行詳細(xì)的統(tǒng)計(jì)分析，其中包括對(duì)降水回波強(qiáng)度、回波頂高、回波面積、強(qiáng)回波發(fā)生頻率以及回波發(fā)生頻率分布的統(tǒng)計(jì)。這些雷達(dá)參數(shù)能夠反映對(duì)流系統(tǒng)的強(qiáng)弱、降水范圍的大小，下面分別說(shuō)明這些參數(shù)的計(jì)算方法。

對(duì)回波強(qiáng)度的統(tǒng)計(jì)實(shí)際是對(duì)雷達(dá)單個(gè)體掃的平均回波強(qiáng)度ZVT進(jìn)行統(tǒng)計(jì)，ZVT計(jì)算公式為

受墨脫地形的限制，當(dāng)海拔高度較低時(shí)，XPAR 可探測(cè)到的回波面積十分有限。另外，墨脫地區(qū)零度亮帶層高度較低，因此選擇某一高度層來(lái)統(tǒng)計(jì)回波區(qū)域面積較為困難。為了了解回波區(qū)域面積的客觀變化規(guī)律，選擇對(duì)雷達(dá)體掃的組合反射率面積（CR）進(jìn)行統(tǒng)計(jì)。CR 是指在一個(gè)雷達(dá)體掃過(guò)程中，將常定仰角方位角掃描中發(fā)現(xiàn)的最大反射率因子投影到笛卡爾格點(diǎn)上的產(chǎn)品（俞小鼎等,2007）。在資料分析過(guò)程中，首先將體掃資料通過(guò)雷達(dá)投影變換，轉(zhuǎn)換為網(wǎng)格間距30 m 的柵格資料，計(jì)算12 個(gè)仰角層的CR，最終形成網(wǎng)格間距為30 m的CR 柵格資料（肖艷姣和劉黎平, 2006）。那么雷達(dá)體掃CR 回波面積ACR為

其中，ΔA=30 m×30 m=0.0009 km2，為柵格面積，N為柵格數(shù)。

雷達(dá)強(qiáng)回波區(qū)域往往預(yù)示著強(qiáng)對(duì)流活動(dòng)的生成與發(fā)展，選取回波強(qiáng)度＆gt;30 dBZ的值表示強(qiáng)回波。對(duì)雷達(dá)強(qiáng)回波發(fā)生頻率的統(tǒng)計(jì)實(shí)際是對(duì)雷達(dá)單個(gè)體掃中＆gt;30 dBZ回波的發(fā)生頻率FZ30的統(tǒng)計(jì)，F(xiàn)Z30的計(jì)算公式如下：

其中，NZ30表示單個(gè)體掃中回波強(qiáng)度大于30 dBZ的點(diǎn)數(shù)，NZ則表示雷達(dá)體掃在墨脫地形下理論可探測(cè)回波點(diǎn)總數(shù)。

根據(jù)2.2.2 節(jié)對(duì)雷達(dá)有效探測(cè)回波的分析，可確定雷達(dá)在某時(shí)間段內(nèi)、某海拔高度理論探測(cè)數(shù)據(jù)格點(diǎn)數(shù)。據(jù)此統(tǒng)計(jì)不同月份、不同小時(shí)降水回波強(qiáng)度、頂高發(fā)生頻率的分布情況，以及雷達(dá)降水回波的垂直分布規(guī)律。下面給出降水回波（分辨率為1 dBZ）月變化（PMi-Zj）和日變化（PHi-Z j）的出現(xiàn)頻率分布計(jì)算方法：

其中，ETj表示回波頂高等于j（單位：km），NMi-ETj為第i月回波頂高等于j的點(diǎn)數(shù)，NHi-ETj為第i小時(shí)回波強(qiáng)度等于j的點(diǎn)數(shù)，NMi-ET為第i月中所有雷達(dá)體掃笛卡爾坐標(biāo)格點(diǎn)總數(shù)，NHi-ETj第i小時(shí)中所有雷達(dá)體掃笛卡爾直角坐標(biāo)格點(diǎn)的總和。

X-PAR 不同海拔高度對(duì)應(yīng)的降水回波（分辨率1 dBZ）出現(xiàn)頻率PHgi-Z j的計(jì)算方法為

其中，Hgi表示海拔高度i（單位：km），NHgi-Zj在高度i探測(cè)到的回波強(qiáng)度為j的點(diǎn)數(shù)，NHgi為理論上雷達(dá)在海拔高度i時(shí)能夠探測(cè)到的回波總數(shù)。

3 墨脫地區(qū)云降水宏觀特征統(tǒng)計(jì)分析

3.1 云降水宏觀特征的月變化

為了探究西藏墨脫云降水宏觀特征，首先對(duì)該地區(qū)云降水特征逐月變化情況進(jìn)行詳細(xì)的統(tǒng)計(jì)分析，其中包括對(duì)不同月份降水回波強(qiáng)度、頂高發(fā)生頻率分布的統(tǒng)計(jì)與降水回波強(qiáng)度、頂高、面積、強(qiáng)回波發(fā)生頻率逐月變化箱型圖的分析。

不同月份降水回波發(fā)生頻率分布（圖6a）顯示，2019 年11 月至2020 年3 月降水回波強(qiáng)度主要集中在11～31 dBZ，2020 年4 月至2020 年10月集中在5～33 dBZ。不同月份降水回波頂高頻率分布（圖6b）顯示，2019 年11 月至2020 年3 月頂高主要分布在1～5 km，最大頂高頻率對(duì)應(yīng)3 km處。2020 年4 月至2020 年10 月頂高分布在1～7 km，其中4～5 月與8～10 月最大頂高頻率對(duì)應(yīng)頂高分布在3～4 km，5～6 月最大頂高頻率對(duì)應(yīng)頂高為4～5 km。

根據(jù)圖6 可知，墨脫地區(qū)降水回波強(qiáng)度、頂高的分布范圍在2020 年4 月至2020 年10 月大于2019 年11 月至2020 年3 月。進(jìn)入4 月后，降水回波發(fā)生頻率突然增大，最大頂高頻率對(duì)應(yīng)高度也隨之升高，其中6 月達(dá)到最大，隨后開始減小。據(jù)此，認(rèn)為墨脫4～10 月降水較11～3 月多，且云降水的垂直發(fā)展更加旺盛，其中5～7 月在一年中降水最旺盛。

為了檢驗(yàn)圖6 的分析結(jié)果是否合理，對(duì)相控陣?yán)走_(dá)觀測(cè)時(shí)間段內(nèi)墨脫站的月降水量進(jìn)行統(tǒng)計(jì)，結(jié)果如圖7 所示。墨脫月降水量時(shí)間變化圖顯示其降水量在4～10 月遠(yuǎn)大于11～3 月，其中6、7 月為一年中墨脫降水量最大的月份，其隨時(shí)間變化趨勢(shì)與圖6 的分析結(jié)果較為一致。另外，圖7 中8 月降水量的減少同時(shí)解釋了圖6 中8 月降水回波發(fā)生頻率與回波頂高分布的突然減小。

圖6 不同月份（a）回波強(qiáng)度與（b）回波頂高發(fā)生頻率的分布Fig.6 Frequency distribution of (a) echo intensity and (b) echo top height in different months

圖7 墨脫站月降水量隨時(shí)間變化Fig.7 Time variation diagram of monthly precipitation at Motuo station

在完成對(duì)不同月份降水回波強(qiáng)度、頂高發(fā)生頻率分布的討論后，統(tǒng)計(jì)分析2019 年11 月至2020年10 月雷達(dá)降水回波數(shù)據(jù)強(qiáng)度、頂高、面積以及強(qiáng)回波發(fā)生頻率的逐月變化情況，如圖8 所示。箱型圖中最高最低兩點(diǎn)分別為最大值和最小值，盒子上下橫線分別為上四分位數(shù)點(diǎn)（75%）與下四分位數(shù)點(diǎn)（25%），盒子中間橫線為中位數(shù)點(diǎn)（50%）。

圖8a 中折線為月平均回波強(qiáng)度的連線，其顯示X-PAR 平均回波強(qiáng)度在11～3 月穩(wěn)定在26 dBZ左右，到4 月平均值有所下降，4～9 月平均值在22 dBZ附近波動(dòng)，10 月平均值上升至25 dBZ。圖8a中的上四分位數(shù)點(diǎn)與下四分位數(shù)點(diǎn)組成的箱子顯示11～3 月與10 月回波強(qiáng)度主要分布在20～31 dBZ，4～9 月上四位數(shù)點(diǎn)值下降至25 dBZ左右。圖中顯示降水回波強(qiáng)度箱型圖與平均值折線的逐月變化趨勢(shì)相同，從4 月開始，墨脫地區(qū)的降水回波強(qiáng)度整體有所減弱，該變化持續(xù)到9 月。

圖8b 中降水回波頂高箱型數(shù)據(jù)與平均值逐月變化趨勢(shì)相同。平均頂高與箱型數(shù)據(jù)均是從11 月開始波動(dòng)上升，6 月平均值達(dá)到最高的4 km，上四分位點(diǎn)與下四分位點(diǎn)也在6 月達(dá)到最大值，7 月后開始持續(xù)下降。圖8b 指出6 月墨脫降水垂直發(fā)展最為旺盛。

圖8 觀測(cè)數(shù)據(jù)箱型圖的逐月變化，箱型圖中最高最低兩點(diǎn)分別為最大值和最小值，盒子上下橫線分別為上四分位數(shù)點(diǎn)（75%）與下四分位數(shù)點(diǎn)（25%），盒子中間橫線為中位數(shù)點(diǎn)。（a）回波強(qiáng)度，圖中折線為月平均回波強(qiáng)度的連線；（b）回波頂高，折線為月平均回波頂高的連線；（c）回波面積，折線為月平均回波面積值連線；（d）強(qiáng)回波發(fā)生頻率，折線為每月強(qiáng)回波發(fā)生頻率平均值的連線Fig.8 Monthly variation of box chart of observational data; the highest and lowest points in the box diagram indicate the maximum and minimum values, upper and lower horizontal lines of the box represent the upper quartile points (75%) and lower quartile points (25%), respectively, and the middle line of the box represents the median point.(a) echo intensity, the line is the average of echo intensity per month in the figure; (b) echo top height, the line is the average of echo top height per month; (c) echo area, the line is the average of echo area per month; and (d) frequency of strong echo, the line is the average of the frequency of strong echo per month

圖8c 中折線為月平均回波面積連線。折線顯示平均面積從11 月開始波動(dòng)上升至6 月達(dá)到最大值754.1 km2，7 月后開始波動(dòng)下降。回波面積箱型圖（圖8c）中的數(shù)據(jù)上四分位數(shù)點(diǎn)從11 月開始波動(dòng)上升，至4 月達(dá)到最大值，隨后波動(dòng)下降，下四分位數(shù)點(diǎn)則是在6 月與7 月有明顯的增加，其中6 月為下四分位數(shù)點(diǎn)最大月份。圖8c 中平均面積與箱型圖逐月變化趨勢(shì)近似一致，4～7 月為一年中回波面積最大的四個(gè)月，即該觀測(cè)時(shí)間段內(nèi)降水范圍分布較廣，其中6 月降水分布最廣。

墨脫地區(qū)強(qiáng)回波所占比例逐月變化圖（圖8d）中平均值與箱型數(shù)據(jù)變化較為一致，箱型數(shù)據(jù)與平均值廓線均顯示強(qiáng)回波比例自3 月起開始增加，6月強(qiáng)回波所占比例達(dá)到最大，隨后開始減小。據(jù)此認(rèn)為6 月是一年中是強(qiáng)降水發(fā)生最多的月份。

以上分析表明墨脫地區(qū)云降水宏觀特征具有逐月變化規(guī)律。雷達(dá)探測(cè)到的降水回波頂高、面積以及強(qiáng)回波發(fā)生頻率均從11 月起逐步增大，在6 月達(dá)到最大值，隨后波動(dòng)下降。該現(xiàn)象表明在進(jìn)入4月后，降水頻次、對(duì)流性降水、降水范圍均在逐漸增大，其中以6 月最為顯著。但從4 月起降水回波強(qiáng)度值卻有所下降，這是由于4 月后墨脫降水大量增加，增加的以弱降水回波為主（圖6a）。圖8分析結(jié)果與圖6 一致。

3.2 云降水宏觀特征日變化

根據(jù)上述3.1 節(jié)對(duì)墨脫X-PAR 月降水回波強(qiáng)度、頂高發(fā)生頻率分布的統(tǒng)計(jì)與回波強(qiáng)度、頂高、面積與強(qiáng)回波發(fā)生頻率逐月變化特征的分析，發(fā)現(xiàn)該地云降水宏觀特征具有明顯的逐月變化規(guī)律。結(jié)合Zhou et al.（2017）提出的高原季風(fēng)指數(shù)（ZPMI），其中ZPMI 在4～10 月為正，在11～3 月為負(fù)，將一年分為雨季（4～10 月）與旱季（11～3 月）。隨后對(duì)旱季與雨季中云降水宏觀特征的日變化規(guī)律進(jìn)行討論。

圖9 給出了兩個(gè)季節(jié)降水回波強(qiáng)度、頂高發(fā)生頻率分布的日變化，其顯示雨季回波強(qiáng)度、頂高分布范圍大于旱季，雨季降水回波發(fā)生頻次明顯高于旱季。降水的發(fā)生頻次與超過(guò)4 km 的回波頂高有明顯的日變化，旱季降水回波主要發(fā)生在18:00～20:00 與00:00～02:00，雨季發(fā)生在23:00～09:00。在00:00～09:00 雨季降水發(fā)生頻次近乎是旱季的2 倍。

圖9 不同時(shí)間回波強(qiáng)度、頂高發(fā)生頻率分布：（a）旱季回波強(qiáng)度；（b）雨季回波強(qiáng)度；（c）旱季回波頂高；（d）雨季回波頂高Fig.9 Frequency distribution of echo intensity and echo top height at different time.Echo intensity in (a) dry season and (b) rainy season; echo top height in (c) the dry season and (d) the rainy season

為了定量分析這些參量的日變化，圖10 給出了回波強(qiáng)度、頂高、面積以及強(qiáng)回波發(fā)生頻率平均值的日變化，圖中空心菱形表示旱季，實(shí)心菱形表示雨季。其中雨季的平均回波頂高約比旱季高1 km（圖10b），回波面積比旱季大150 km2左右（圖10c），強(qiáng)回波所占比例始終大于旱季，說(shuō)明雨季降水垂直發(fā)展更加旺盛，強(qiáng)降水出現(xiàn)次數(shù)增加，降水范圍增大。圖10a 中雨季的平均回波強(qiáng)度約比旱季小1 dBZ，這是因?yàn)檫M(jìn)入雨季后，降水頻率較旱季顯著增加，但增加的降水主要以弱回波（15～25 dBZ）降水為主（圖9a），從而導(dǎo)致雨季降水回波強(qiáng)度平均值減小。兩季平均頂高日變化趨勢(shì)近似一致，均是在下午（旱季17:00，雨季15:00）達(dá)到最高值，夜晚存在有頂高次峰值（旱季00:00，雨季01:00），兩季的日最強(qiáng)對(duì)流均發(fā)生在下午。在旱季，平均回波面積日變化趨勢(shì)與頂高近似相同，面積自上午（10:00）開始增加，下午（19:00）達(dá)到最大，凌晨（02:00）開始減小。雨季的回波面積日變化趨勢(shì)不同于頂高，平均面積從夜晚（22:00）開始增加，在凌晨（03:00）達(dá)到最大。上述參量日變化情況結(jié)合圖9 對(duì)回波發(fā)生頻次的分析，得出旱季日降水主要出現(xiàn)在下午和上半夜。對(duì)于雨季，雖然雨季下午強(qiáng)對(duì)流活動(dòng)旺盛，但其下午的回波面積與回波發(fā)生頻率遠(yuǎn)小于夜晚，因此認(rèn)為雨季日降水峰值主要出現(xiàn)在下半夜。

圖10 觀測(cè)數(shù)據(jù)平均值的日變化：（a）回波強(qiáng)度；（b）回波頂高；（c）回波面積；（d）強(qiáng)回波發(fā)生頻率Fig.10 Daily variation of the average of the observed data: (a) Echo intensity; (b) echo top height; (c) echo area; (d) frequency of strong echo

3.3 云降水垂直分布特征

3.3.1 旱季與雨季云降水垂直分布特征

為了探究墨脫地區(qū)云降水垂直分布特征，對(duì)回波強(qiáng)度的垂直分布進(jìn)行了統(tǒng)計(jì)，結(jié)果如圖11 所示。旱季回波強(qiáng)度分布在10～30 dBZ，回波發(fā)生海拔高度集中在2～5 km。雨季回波強(qiáng)度范圍為5～33 dBZ，高度2～7 km。雨季回波頻率增加部分回波強(qiáng)度主要分布在10～30 dBZ，對(duì)應(yīng)高度在2～3 km。對(duì)比圖11a 與圖11b 發(fā)現(xiàn)，雨季的回波發(fā)生頻率、回波強(qiáng)度分布范圍以及回波高度范圍均大于旱季，即雨季相較旱季降水回波數(shù)量明顯增多，對(duì)流降水活動(dòng)發(fā)展更加旺盛。

圖11 觀測(cè)時(shí)間段內(nèi)回波強(qiáng)度的垂直分布：（a）旱季；（b）雨季Fig.11 Vertical distribution of echo intensity in observation period during (a) the dry season and (b) the rainy season

3.3.2 云降水垂直分布特征日變化

為了進(jìn)一步探究墨脫地區(qū)云降水垂直分布日變化特征，統(tǒng)計(jì)兩個(gè)季節(jié)的不同高度回波強(qiáng)度分布的小時(shí)變化，結(jié)果如圖12 與圖13 所示。圖12 旱季日變化顯示，降水發(fā)生頻率從14:00～15:00 開始增加（圖12h），18:00～19:00 頻率達(dá)到最大（圖12j），隨后開始減小，00:00～00:01 其值再次增大（圖12a），04:00～13:00（圖12c-g）頻率持續(xù)減小。頻率增大區(qū)域回波強(qiáng)度主要分布在15～25 dBZ，海拔高度分布在3～4.5 km。圖13 雨季日變化顯示頻率在18:00～19:00 開始增大（圖13j），00:00～07:00達(dá)到一天中的最大時(shí)段（圖13a-d），隨后頻率持續(xù)減小。頻率增大的區(qū)域回波強(qiáng)度分布在10～30 dBZ，海拔高度分布在2～3 km。

對(duì)比旱季與雨季回波強(qiáng)度垂直分布日變化特征發(fā)現(xiàn)（圖12 與圖13），雨季降水回波發(fā)生頻率整體在各個(gè)時(shí)刻均大于旱季。旱季回波發(fā)生頻率存在兩個(gè)日峰值，分別為當(dāng)?shù)貢r(shí)下午16:00～17:00（圖12j）與夜晚00:00～00:01（圖12b）。雨季，在當(dāng)?shù)貢r(shí)夜晚22:00～05:00（圖13a-d）頻率最大。對(duì)比兩季回波發(fā)生頻率增大區(qū)域，發(fā)現(xiàn)旱季3 km 以上回波頻次高，而雨季3 km 以下較高。據(jù)此，認(rèn)為旱季有兩個(gè)日降水峰值分別是下午與上半夜，雨季日降水峰值則發(fā)生在下半夜。雨季雷達(dá)降水回波數(shù)目遠(yuǎn)大于旱季。雨季降水更多，更旺盛。

圖12 旱季回波強(qiáng)度垂直分布的逐時(shí)變化Fig.12 Hourly variation of the vertical distribution of echo intensity in the dry season

圖13 雨季回波強(qiáng)度垂直分布的逐時(shí)變化Fig.13 Hourly variation of the vertical distribution of echo intensity in the rainy season

3.4 日變化結(jié)果分析與對(duì)比討論

通過(guò)定量分析墨脫X-PAR 觀測(cè)數(shù)據(jù)的日變化發(fā)現(xiàn)旱季日降水主要發(fā)生在午后與上半夜，雨季降水以下半夜降水為主。墨脫地區(qū)三面環(huán)山，呈高山峽谷地形，是雅魯藏布峽谷水汽通道的主體入口。受水汽通道的影響，該地區(qū)水汽充足，異常濕潤(rùn)，長(zhǎng)期云霧繚繞。考慮到墨脫地形與該地的氣候認(rèn)為墨脫旱季雨季夜晚降水多是由云頂輻射強(qiáng)迫冷卻（Bhatt and Nakamura, 2006）與夜間山風(fēng)的盛行（Chen et al., 2012）共同引起的。在季風(fēng)爆發(fā)后，來(lái)自印度洋的大量水汽涌入墨脫使墨脫雨季降水更加頻繁，對(duì)流活動(dòng)發(fā)展更加旺盛。

Chen et al.（2012）分析夏季季風(fēng)時(shí)期（6～8 月）藏東南山谷站雨量計(jì)數(shù)據(jù)，結(jié)果表明在27個(gè)臺(tái)站中的24 個(gè)臺(tái)站顯示夜間峰值主導(dǎo)了這些臺(tái)站的日降雨模式。平均降雨量日變化曲線顯示在當(dāng)?shù)貢r(shí)03:00 時(shí)達(dá)到日最大值。另外，這些臺(tái)站在下午17:00 存在有弱峰值。圖14 統(tǒng)計(jì)了墨脫X-PAR 6～8 月觀測(cè)數(shù)據(jù)的日變化情況。夏季季風(fēng)時(shí)期平均回波頂高自上午10:00 開始升高，下午15:00 達(dá)到最高平均頂高3.6 km，次頂高峰值出現(xiàn)在01:00（圖14a）。平均回波面積與頂高日變化趨勢(shì)不同，其自16:00 起增大，最大面積平均值770 km2出現(xiàn)在凌晨04:00（圖14b）。回波發(fā)生頻率在不同小時(shí)間的分布指出6～8 月回波數(shù)量主要出現(xiàn)在00:00～09:00（圖14c）。頂高發(fā)生頻率日變化趨勢(shì)分布與回波強(qiáng)度的發(fā)生頻率對(duì)應(yīng)關(guān)系較為一致。圖14 的分析中指出在夏季季風(fēng)時(shí)期回波面積在夜晚大，回波頂高高，回波發(fā)生頻率大，因而得出夜晚降水峰值的結(jié)論。該雷達(dá)回波面積最大值與回波數(shù)量峰值出現(xiàn)時(shí)間與Chen et al.（2012）指出的降雨夜間峰值時(shí)間較為一致。

圖14 夏季季風(fēng)期間（6～8 月）(a)平均回波頂高、(b)平均回波面積日變化曲線與不同時(shí)間(c)回波強(qiáng)度與(d)回波頂高的分布Fig.14 Diurnal variation of the average of (a) echo intensity and (b) echo area, and the distribution of (c) echo intensity and (d) echo top height at different time during the summer monsoon (June-August)

Uyeda et al.（2001）與馮錦明等（2002）通過(guò)分析GAME-Tibet 期間（6～8 月）在高原主體那曲地區(qū)獲取的多普勒雷達(dá)資料，指出那曲地區(qū)回波頂高與面積日變化趨勢(shì)一致，從當(dāng)?shù)貢r(shí)06:00 開始增大，當(dāng)?shù)貢r(shí)下午14:00 為頂高最高、面積最大的時(shí)刻，此時(shí)云中對(duì)流活動(dòng)最強(qiáng)，夜晚到次日凌晨逐漸消散，呈單峰變化。圖14 墨脫夏季季風(fēng)時(shí)期XPAR 觀測(cè)數(shù)據(jù)統(tǒng)計(jì)結(jié)果與那曲相比較發(fā)現(xiàn)，那曲地區(qū)的回波頂高高于墨脫，墨脫的頂高與面積日變化趨勢(shì)與那曲并不相同。墨脫頂高平均值日變化顯示一天中頂高存在兩個(gè)峰值分別為下午的最強(qiáng)峰值與夜晚的次峰值。墨脫回波面積雖然呈單峰變化但區(qū)別于那曲的下午峰值，峰值出現(xiàn)在凌晨。另外，兩地區(qū)日降水峰值也存在差異，那曲降水多發(fā)生在下午與上半夜，墨脫降水多發(fā)生在下半夜。兩地夜間降水形成的原因也不相同，那曲夜間降水主要是由于對(duì)流云開始消散轉(zhuǎn)向?qū)訝钤平邓ǔ５t和郭學(xué)良, 2016），而山谷風(fēng)與輻射冷卻是導(dǎo)致墨脫夜間降水的主要原因。本文對(duì)墨脫旱季日變化特征的研究中指出旱季回波頂高、回波面積均在下午達(dá)到最大值，夜晚至凌晨消散，降水多發(fā)生在下午及上半夜，其云降水特征與那曲夏季風(fēng)時(shí)期降水特征較為相似。

4 結(jié)論

為了探究藏東南墨脫地區(qū)的云降水宏觀特征，本文利用X-PAR 在該地區(qū)獲取的2019 年11 月至2020 年10 月一整年的觀測(cè)數(shù)據(jù)，對(duì)該地云降水宏觀特征的月變化、日變化、垂直變化進(jìn)行了統(tǒng)計(jì)并完成了初步討論，得到了以下結(jié)論：

（1）墨脫地區(qū)云降水特征具有逐月變化規(guī)律。X-PAR 降水回波頂高、面積、強(qiáng)回波發(fā)生頻率自4 月起開始增加，6 月達(dá)到最大值，隨后開始下降。一年中降水以6 月最為顯著，11、12 月降水最少、對(duì)流活動(dòng)少。4～9 月中增加的降水回波主要以弱回波（15～25 dBZ）為主，該現(xiàn)象導(dǎo)致4～9 月降水回波平均值略小于10～3 月。根據(jù)降水回波月變化分析討論結(jié)果同時(shí)結(jié)合高原季風(fēng)指數(shù)，將一年分為旱季（11～3 月）與雨季（4～10 月）。

（2）通過(guò)對(duì)墨脫地區(qū)云降水回波日變化特征的研究可知，墨脫雨季降水回波發(fā)生頻率、頂高、面積、強(qiáng)回波所占比例均大于旱季，且降水回波強(qiáng)度、頂高的分布范圍也大于旱季，表明雨季降水?dāng)?shù)量多于旱季，對(duì)流活動(dòng)頻繁。另旱季、雨季的回波頂高日變化趨勢(shì)相似，自上午10:00 起升高，分別在下午17:00 與15:00 達(dá)到最高，夜晚次峰值分別出現(xiàn)在00:00 與01:00。兩季回波面積日變化趨勢(shì)不同，旱季最大回波面積出現(xiàn)在19:00，面積次峰值在01:00，雨季面積自22:00 起增大，在3:00 左右達(dá)到最大。降水回波發(fā)生頻率日分布顯示旱季在18:00 與01:00 回波發(fā)生頻繁，雨季03:00 回波最頻繁。旱季降水回波頂高、面積、發(fā)生頻率在下午與夜晚均存在有較大值，認(rèn)為旱季降水主要發(fā)生在下午與上半夜。雨季雖然頂高在下午與夜晚較高，但面積與回波發(fā)生頻率在下午值較小，在凌晨有最大值，認(rèn)為雨季降水主要發(fā)生在下半夜。

（3）兩季節(jié)回波強(qiáng)度垂直分布顯示墨脫降水回波強(qiáng)度大小幾乎不超過(guò)30 dBZ，旱季降水回波主要發(fā)生在3 km 以上，雨季回波多發(fā)生在3 km以下。

（4）通過(guò)對(duì)比墨脫X-PAR 與那曲多普勒雷達(dá)在夏季季風(fēng)期間降水回波頂高、面積的日變化發(fā)現(xiàn)，夏季季風(fēng)期間那曲地區(qū)的回波頂高高于墨脫地區(qū)，兩地頂高與面積的日變化趨勢(shì)存在差異。那曲回波頂高與面積自當(dāng)?shù)貢r(shí)上午10:00 開始增大，下午14:00 達(dá)到最大，在夜晚至凌晨消散。而墨脫頂高日變化顯示頂高有兩個(gè)峰值，在當(dāng)?shù)貢r(shí)13:00 有最大頂高，23:00 為頂高次峰值。墨脫面積自當(dāng)?shù)貢r(shí)間19:00 起增大，凌晨02:00 達(dá)到最大回波面積。夏季季風(fēng)期間兩地日最強(qiáng)降水出現(xiàn)時(shí)間也不同，墨脫最強(qiáng)降水出現(xiàn)在下半夜，而那曲出現(xiàn)在下午與上半夜。那曲夏季季風(fēng)期間的云降水日變化特征近似墨脫旱季。