段麗君， 申紅艷， 余 迪， 馬有絢， 白文蓉， 李萬(wàn)志

（1.青海省氣候中心，青海西寧810001；2.青海省防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室，青海西寧810001）

0 引言

青藏高原（以下簡(jiǎn)稱“高原”）平均海拔4 000多米，在垂直方向上約占對(duì)流層厚度的四分之一。從喜馬拉雅山南側(cè)的常綠林到北部的荒漠草原、高山冰雪帶，地形地貌復(fù)雜，自然帶垂直差異顯著，植被分布不均勻、種類繁多。境內(nèi)湖泊面積和冰川儲(chǔ)量分別占中國(guó)總量的52%和80%，素有“世界屋脊”［1］和“亞洲水塔”之稱［2］。高原及周邊地區(qū)也是全球重要的冰川資源富集地［3］，冰川對(duì)于氣候變化有敏感的指示作用，冰凍圈對(duì)水資源短缺的中國(guó)西部地區(qū)的水文與水資源影響巨大［4］。

降水作為影響高原地區(qū)自然生態(tài)系統(tǒng)最活躍、最直接的氣候因子，對(duì)生態(tài)環(huán)境具有顯著影響，同時(shí)降水對(duì)我國(guó)西北地區(qū)及高原的徑流變化也存在一定影響［4］。水汽是形成降水的必要條件，水汽輸送與收支是影響區(qū)域水分平衡的直接因素［5-7］。夏季高原最強(qiáng)的熱力壓迫作用形成較強(qiáng)的低壓槽［8-11］。夏季亞洲季風(fēng)區(qū)為強(qiáng)大的水汽匯，南支西風(fēng)向東亞輸送水汽［12-13］。夏季高原上空為明顯的大氣水汽含量高中心，4—9月高原可降水量變化顯著，高原的增濕速度小于減濕速度，來(lái)自阿拉伯海-印度洋和南海的水汽在孟加拉灣匯合后所形成的偏南風(fēng)水汽，在高原南緣和東南緣分別存在經(jīng)向、緯向不同分量水汽流型，進(jìn)而影響高原地區(qū)的降水分布［14-16］。陳濤等［17］指出，青藏高原西部春季降水與其西南部輻合上升運(yùn)動(dòng)及阿拉伯海北部水汽含量變化存在聯(lián)系。南亞季風(fēng)對(duì)青藏高原東北緣地區(qū)降水量影響顯著，為顯著正相關(guān)［18］。荀學(xué)義等［19］研究指出，高原夏季降水的起訖與高原及北側(cè)氣壓系統(tǒng)、高原東南或南部水汽輸送條件息息相關(guān)。解承瑩等［20］提出“青藏高原南緣水汽輸送關(guān)鍵區(qū)”并討論關(guān)鍵區(qū)各邊界水汽收支與高原夏季降水分布的關(guān)系，發(fā)現(xiàn)各邊界水汽收支與印度熱低壓和南海夏季風(fēng)活動(dòng)關(guān)系密切。

高原旱季、雨季分明，降雨主要集中在5—9月，占全年總雨量80%以上，高原降水由東南向西北遞減，而且年際變化具有一定的多元性［21-22］。高原雨季也是全年中氣溫高、濕度大、風(fēng)速小的時(shí)段，是農(nóng)作物、牧草等生長(zhǎng)發(fā)育的最佳時(shí)期，因此雨季是一年中高原地區(qū)非常重要的階段。雨季持續(xù)時(shí)間和雨量強(qiáng)弱直接關(guān)系高原地區(qū)的旱澇異常，對(duì)高原自然生態(tài)系統(tǒng)、水資源，以及三江源區(qū)、祁連山區(qū)、環(huán)青海湖等國(guó)家重點(diǎn)生態(tài)功能區(qū)和絲綢之路經(jīng)濟(jì)帶的區(qū)域生態(tài)環(huán)境等有重要影響。徐國(guó)昌等［23］研究指出，高原是我國(guó)雨季最顯著的地區(qū)。章凝丹等［24］基于1951—1970年降水資料分析發(fā)現(xiàn)，高原雨季自東南至西北開始，結(jié)束正好與此相反，因此雨季具有西部短、東部長(zhǎng)的特點(diǎn)，且高原雨季年際變化很大?？妴埖龋?5］分析得出高原地區(qū)1961—2000年降水呈增加趨勢(shì)，在1978年由少雨期轉(zhuǎn)為多雨期，且高原南部降水量增加明顯，而北部變化較小；高原春季降水量年際變化較大，夏季降水量值較大但變化幅度較小，冬季變化與夏季相反［26］。

由以上研究可以看出，青藏高原及周邊地區(qū)的水汽輸送對(duì)降水起著重要作用。前人較多關(guān)注了高原雨季降水最主要的分布型水汽條件，而對(duì)其他分布型水汽條件研究偏少。本文擬在已有研究的基礎(chǔ)上，分析高原雨季降水在不同分布型時(shí)的水汽來(lái)源、輸送路徑和水汽收支特征，以及對(duì)高原主體水汽輸送在雨季不同階段有重要貢獻(xiàn)的因子，并討論雨季期間高原各邊界的水汽貢獻(xiàn)，對(duì)認(rèn)識(shí)高原雨季水資源的分布和利用有重要意義。

1 資料與方法

1.1 資料

本文所用到的降水資料為中國(guó)國(guó)家氣象信息中心整編的中國(guó)青藏高原地區(qū)地面氣象觀測(cè)站（基準(zhǔn)站、基本站、一般站）從建站至2017年109站逐日（20：00—20：00，北京時(shí)間）降水資料。其中，青海50站、西藏38站、四川17站、甘肅4站（圖1），站點(diǎn)名稱等信息見表1。青藏高原區(qū)域多數(shù)臺(tái)站始建于1958—1960年，1961年開始觀測(cè)，考慮到記錄的完整性和可比性，選取1961—2017年有較完整降水資料的站進(jìn)行分析。個(gè)別站1967年、1968年資料有缺失，部分站20世紀(jì)70年代初建站，針對(duì)這類缺測(cè)情況選取臨近站點(diǎn)利用線性回歸和比值訂正法進(jìn)行訂正。

圖1 青藏高原氣象站點(diǎn)分布Fig.1 The distribution of weather stations over the Tibetan Plateau（shaded：altitude，unit：m）

表1 青藏高原氣象站點(diǎn)（109站）信息表Table 1 Information of the 109 meteorological stations on the Tibetan Plateau

格點(diǎn)資料采用NCEP/NCAR（National Centers for Environmental Prediction/National Center for Atmospheric Research）提供的1961—2017年逐月再分析資料，垂直方向從1 000 hPa到300 hPa，包括緯向風(fēng)、經(jīng)向風(fēng)、位勢(shì)高度、相對(duì)濕度、地表氣壓等常規(guī)變量，格點(diǎn)資料的水平分辨率統(tǒng)一插值為2.5°×2.5°的 粗 網(wǎng) 格 數(shù) 據(jù)［27］。文 中 高 原 雨 季 指5—9月。

1.2 方法

1.2.1 水汽計(jì)算方法

（1）水汽通量及其散度[28-30]

水汽通量由下式進(jìn)行計(jì)算。

其中，包括緯向水汽通量Qλ和經(jīng)向水汽通量Qφ。

水汽通量散度由下式進(jìn)行計(jì)算。

式中：g為重力加速度；u為緯向風(fēng)；v為經(jīng)向風(fēng)；ps為下邊界氣壓；pu為上邊界氣壓；q為比濕；V為單位氣柱各層大氣的風(fēng)速矢量；λ為緯度；φ為經(jīng)度；a為地球半徑。

（2）水汽收支

各邊界積分的水汽輸送表示為

區(qū)域總的水汽收支為

式中：QW、QE、QS、QN分別為西邊界、東邊界、南邊界、北邊界的水汽收支；QT為總的凈水汽收支；λ1、λ2、φ1、φ2分別為各邊界對(duì)應(yīng)的緯度和經(jīng)度；x、y為所選取經(jīng)緯度內(nèi)的格點(diǎn)數(shù)；t為研究時(shí)段內(nèi)的年份數(shù)。

1.2.2 比值訂正法

兩個(gè)距離遠(yuǎn)的臺(tái)站，并非所有氣象要素的插值都趨近于常數(shù)，但它們所對(duì)應(yīng)的比值是相對(duì)穩(wěn)定的，即

式（6）的變化較小，趨近于常數(shù)，所表示的意義就是y增加的倍數(shù)和x增加的倍數(shù)逐年或逐月是一樣的［31］。

1.2.3 統(tǒng)計(jì)方法

文中主要采用經(jīng)驗(yàn)正交函數(shù)分解法（empirical orthogonal function，EOF）、去趨勢(shì)、合成分析及顯著性t檢驗(yàn)等常規(guī)數(shù)理統(tǒng)計(jì)方法［32-35］。

2 青藏高原雨季水汽輸送特征

2.1 雨季平均水汽變化

首先了解高原雨季降水的多年變化趨勢(shì)特征。為了能夠更為清晰的研究其是否具有顯著的年際變化特征，以33° N為界將高原分為南、北兩個(gè)區(qū)域，將大于等于33°N的站點(diǎn)（58站）劃分為北區(qū)，小于33°N（51站）的站點(diǎn)劃分為南區(qū)。從圖2可以看出，去除年代際變化后的高原南區(qū)、北區(qū)雨季降水的年際變化更明顯，以高原南區(qū)、北區(qū)雨季降水呈反位相變化為主要特征，其中以1962年、1967年、1990年、1995年南、北區(qū)域反位相變化最為明顯，即高原南區(qū)雨季降水最?。ù螅r(shí)，高原北區(qū)雨季降水為最大（?。５?，從圖2也可以看出，高原雨季降水也存在南、北區(qū)雨季降水一致偏多（少）的變化特征，尤其是在1997年之后，即高原北區(qū)雨季降水偏多（少）時(shí)，南區(qū)降水偏多（少）。

圖2 1961—2017年青藏高原南區(qū)（a）、北區(qū)（b）雨季降水去趨勢(shì)后區(qū)域平均降水量年際序列Fig.2 The inter-annual series of region mean precipitation after the detrending of rainy season in southern（a）and northern（b）Tibetan Plateau during 1961—2017

2.2 雨季降水異常的水汽條件

2.2.1 典型年水汽輸送特征

為了驗(yàn)證高原南、北區(qū)雨季降水的一致性是否為主要的分布型，對(duì)所選取的高原109個(gè)站點(diǎn)1961—2017年雨季降水量進(jìn)行EOF分解，得到其主要空間模態(tài)及各模態(tài)對(duì)應(yīng)的時(shí)間系數(shù)，前5個(gè)模態(tài)的累積方差貢獻(xiàn)率達(dá)54.8%（表2），經(jīng)North檢驗(yàn)，能反映出青藏高原雨季降水量的主要空間分布特征。前兩個(gè)模態(tài)能夠體現(xiàn)高原雨季降水典型空間分布型，對(duì)前兩個(gè)模態(tài)時(shí)間系數(shù)進(jìn)行標(biāo)準(zhǔn)化處理，并根據(jù)不同的時(shí)間系數(shù)標(biāo)準(zhǔn)化值定義各模態(tài)正、負(fù)異常年。

表2 青藏高原雨季降水EOF前5個(gè)模態(tài)的方差貢獻(xiàn)率及累積方差貢獻(xiàn)率Table 2 The variance and accumulate variance contribution of the top five modes of EOF for rainy season precipitation on the Tibetan Plateau

圖3 給出了EOF展開后第一、第二模態(tài)對(duì)應(yīng)空間型及其時(shí)間系數(shù)（PC），反映出高原雨季降水的兩種主要的空間分布型。從第一模態(tài)空間分布型［圖3（a）］可以看出，高原雨季降水北多南少，負(fù)值中心位于高原南緣，雅魯藏布江中下游地區(qū)，正值中心位于青海湖南部地區(qū)，方差貢獻(xiàn)為22.1%，這種分布與高原地形及西太平洋副熱帶高壓的位置有關(guān)，當(dāng)西太平洋副熱帶高壓北跳時(shí)，高原北部降水通常是增加的；第一時(shí)間系數(shù)標(biāo)準(zhǔn)化值（PC1）呈現(xiàn)明顯的年代際變化，且正（負(fù)）位相分別對(duì)應(yīng)高原雨季降水北多（少）南少（多）分布型［圖3（c）］。

圖3 1961—2017年青藏高原雨季降水異常EOF分解前二模態(tài)空間分布型（a、b）及其標(biāo)準(zhǔn)化時(shí)間系數(shù)序列（c、d）Fig.3 The spatial distribution patterns（a，b）and their normalized time series（c，d）of the two dominant EOF modes of rainy season in Tibetan Plateau during 1961—2017

第二模態(tài)空間分布型［圖3（b）］反映高原雨季降水干濕變化的一致性，即高原整體降水偏多或偏少，高值中心位于高原腹地，方差貢獻(xiàn)率為12.6%；第二時(shí)間系數(shù)標(biāo)準(zhǔn)化值（PC2）存在明顯的年際和年代際變化特征，且正（負(fù)）位相分別對(duì)應(yīng)高原雨季降水一致偏多（少）分布型，1997年以來(lái)，PC2以正值為主［圖3（d）］，青藏高原在1997年存在明顯的增暖突變［36］，表明氣候變暖背景下高原雨季降水具有一致偏多的響應(yīng)特征。

由以上研究可知，高原降水南北一致分布型可以作為高原雨季降水的主要特征。那么這種分布型對(duì)應(yīng)的水汽條件是怎樣的？因此，下面將著重分析高原雨季降水一致偏多（少）時(shí)的水汽輸送情況。

根據(jù)EOF第二特征向量時(shí)間系數(shù)［圖3（d）］，挑選出多雨年（大于4）6年（1974年、1985年、1989年、1998年、2003年、2004年），少雨年（小于-4）7年（1972年、1973年、1977年、1986年、1994年、1997年、2006年），進(jìn)一步研究多雨、少雨年的水汽輸送的演變及差異。

在對(duì)異常年進(jìn)行研究之前，首先了解氣候態(tài)（1981—2010年）高原雨季期間的水汽輸送形勢(shì)。受多種季風(fēng)等的影響，高原及附近水汽輸送的情況較為復(fù)雜［37-38］。從圖4可以看出，來(lái)自阿拉伯海的偏西風(fēng)水汽輸送在孟加拉灣附近分為3支水汽輸送氣流：一支向北輸送，通過高原南部河谷等有利地形有少量水汽進(jìn)入高原，一支在南海附近轉(zhuǎn)為偏南風(fēng)水汽輸送，一支受高原大地形的阻擋作用轉(zhuǎn)為偏西風(fēng)水汽輸送。

圖4 1981—2010年從地面至300 hPa垂直積分的青藏高原雨季水汽通量（矢量，單位：kg·m-1·s-2）及散度（陰影，單位：10-5 kg·m-2·s-2）分布Fig.4 The distribution of vertically integral water vapor flux（arrows，unit：kg·m-1·s-2）and divergence（shading，unit：10-5 kg·m-2·s-2）of rainy season from surface to 300 hPa in Tibetan Plateau during 1981—2010

圖5 和圖6分別給出了PC2正、負(fù)異常年的水汽通量及其散度場(chǎng)。從圖5可以看出，PC2正異常年，即高原雨季降水偏強(qiáng)年，高原主體呈現(xiàn)北部弱輻合南部輻散的分布形態(tài)，來(lái)自孟加拉灣的偏東南暖濕氣流在高原南部形成反氣旋性環(huán)流，在高原南緣存在水汽輻散中心，南海附近存在偏南風(fēng)水汽輸送，并在高原東部邊緣形成水汽輻合中心。圖6在高原主體呈現(xiàn)出與圖5相反的變化特征，即受來(lái)自阿拉伯海偏南水汽流在25° N附近轉(zhuǎn)為偏東水氣流，在高原南緣形成較強(qiáng)的水汽輻合中心，此時(shí)，高原北部及東部為水汽輻散，并在東部邊緣形成較強(qiáng)的水汽輻散中心。

圖5 PC2正異常年從地面至300 hPa垂直積分的青藏高原雨季水汽通量（矢量，單位：kg·m-1·s-2）及散度（陰影，單位：10-5 kg·m-2·s--2）分布Fig.5 The distribution of vertically integral water vapor flux（arrows，units：kg·m-1·s-2）and divergence（shading，units：10-5 kg·m-2·s-2）of rainy season from surface to 300 hPa in Tibetan Plateau in positive years of the PC2

圖6 PC2負(fù)異常年從地面至300 hPa垂直積分的青藏高原雨季水汽通量（矢量，單位：kg·m-1·s-2）及散度（陰影，單位：10-5 kg·m-2·s-2）分布Fig.6 The distribution of vertically integral water vapor flux（arrows，unit：kg·m-1·s-2）and divergence（shading，unit：10-5 kg·m-2·s-2）of rainy season from surface to 300 hPa in Tibetan Plateau in negative years of the PC2

因此，從以上分析可以看出，降水偏少年高原主體呈現(xiàn)北部水汽輻散南部水汽輻合的特征，降水偏多年，高原主體北部水汽輻合南部水汽輻散，兩者呈相反的變化。

2.2.2 水汽收支特征

雨季期間，高原主體各個(gè)邊界水汽收支的年際變化反映了該區(qū)域內(nèi)水汽輸送強(qiáng)度的變化。用“箱體”模型描述高原雨季水汽收支變化特征，各邊界的區(qū)域?yàn)椋簴|邊界（26°N、104°E至40°N、104°E），西 邊 界（26° N、74° E至40° N、74° E），南 邊 界（26° N、74° E至26° N、104° E），北邊界（40° N、74°E至40°N、104°E）。

圖7 為氣候態(tài)雨季期間高原主體的水汽收支情況，南邊界的水汽輸入最強(qiáng)，西邊界次之，東邊界存在較弱的水汽輸入，北邊界為較弱的水汽輸出。從圖8可以看出，各邊界水汽收支存在顯著的年際變化，北邊界以水汽輸出為主且輸出量較小，個(gè)別年份存在較弱的水汽輸入，該現(xiàn)象與中高緯度地區(qū)的增濕有一定關(guān)系［39］，東邊界以弱水汽輸入為主，西邊界和南邊界水汽輸入量較大，為高原的主要水汽來(lái)源，其中南邊界在1992年水汽輸送量達(dá)到最大，為27.3×106kg·s-1，西邊界在1982年水汽輸送量為最大25.4×106kg·s-1，西邊界與南邊界的水汽輸入在20世紀(jì)80年代至90年代中期呈反位相變化，之后呈大致相同的變化趨勢(shì)?？傮w來(lái)看，西邊界水汽輸入呈增加趨勢(shì)，南邊界水汽輸入變化波動(dòng)較大。結(jié)合以上分析可知，高原雨季期間高原主體的主要水汽來(lái)源為阿拉伯海、南海，次要來(lái)源為西風(fēng)帶。

圖7 1981—2010年氣候態(tài)高原主體各邊界雨季平均的水汽收支Fig.7 Climate mean water budget distribution of the four boundaries of rainy season in Tibetan Plateau during 1981—2010

圖8 1981—2017年高原主體各邊界雨季平均的水汽收支年變化Fig.8 The averaged water budget series over the four boundaries of rainy season in Tibetan Plateau during 1981—2017

3 結(jié)論與展望

通過以上分析，得到以下主要結(jié)論：

（1）1961—2017年高原南區(qū)、北區(qū)雨季降水存在顯著的年際變化，以高原南區(qū)、北區(qū)雨季降水呈反位相變化為主要特征，1997年之后也存在南、北區(qū)雨季降水一致偏多（少）的變化特征。

（2）高原雨季降水存在兩種主要的模態(tài)，第一模態(tài)為青藏高原北多（少）南少（多）型分布，第二模態(tài)為全區(qū)一致型分布。

（3）氣候態(tài)（1981—2010年）雨季平均高原及其鄰近地區(qū)上空的水汽輸送路徑為來(lái)自阿拉伯海的偏西風(fēng)水汽輸送在孟加拉灣附近分為3支水汽輸送氣流：一部分向北輸送，一部分在南海附近轉(zhuǎn)為偏南風(fēng)水汽輸送，一部分受高原大地形的阻擋作用轉(zhuǎn)為偏西風(fēng)水汽輸送。

（4）第二模態(tài)正異常年，高原主體呈現(xiàn)北部弱輻合南部輻散的分布形態(tài)。來(lái)自孟加拉灣的偏東南暖濕氣流在高原南部形成反氣旋性環(huán)流，在高原南緣存在水汽輻散中心，南海附近存在偏南風(fēng)水汽輸送，并在高原東部邊緣形成水汽輻合中心。第二模態(tài)負(fù)異常年，受來(lái)自阿拉伯海偏南水汽流在25° N附近轉(zhuǎn)為偏東水氣流，在高原南緣形成較強(qiáng)的水汽輻合中心。

（5）各邊界水汽收支存在顯著的年際變化，北邊界以水汽輸出為主且輸出量較小，東邊界以弱水汽輸入為主，西邊界和南邊界水汽輸入量較大且在20世紀(jì)80年代至90年代中期呈反位相變化，之后呈大致相同的變化趨勢(shì)。總體來(lái)看，西邊界水汽輸入呈增加趨勢(shì)，南邊界水汽輸入變化波動(dòng)較大。

由上述分析可以看出，本文研究的重點(diǎn)為青藏高原雨季降水第二模態(tài)（全區(qū)一致型）正負(fù)異常年的同期水汽條件及邊界水汽收支情況，文中僅研究了第二模態(tài)的同期水汽條件、已設(shè)定邊界的各邊界水汽收支情況等。值得考慮的是，第二模態(tài)對(duì)應(yīng)的前期及同期環(huán)流信號(hào)有哪些特征？環(huán)流的演變和影響機(jī)理是什么？邊界設(shè)置對(duì)其數(shù)值的是否存在影響？如果將四個(gè)邊界的范圍擴(kuò)大，南邊界的水汽輸入量大于文中的18.2×106kg·s-1，可能是因?yàn)榍嗖馗咴匦蔚挠绊懹兴鶞p弱，來(lái)自孟加拉灣、阿拉伯海等地區(qū)的水汽輸送量較大；西邊界的水汽輸入量略高于文中的13.4×106kg·s-1，說(shuō)明來(lái)自西風(fēng)帶的水汽在越接近高原時(shí)的水汽輸送量越小，但依然是高原雨季水汽的主要來(lái)源之一；東邊界和北邊界范圍擴(kuò)大后，水汽輸入量和輸出量與文中的數(shù)值相比變化很小。從大致的評(píng)估可以看出，邊界設(shè)置對(duì)南邊界、西邊界的影響較大，對(duì)北邊界、東邊界的影響較小。因此，下一步將著重研究第二模態(tài)對(duì)應(yīng)的環(huán)流特征的前期及同期信號(hào)，并診斷其演變和影響機(jī)理，深入評(píng)估邊界設(shè)置對(duì)水汽收支的影響。