石紅艷，楊 寧，潘華盛

(1.雙鴨山市氣象局，黑龍江 雙鴨山 155100；2．黑龍江省氣象服務中心，哈爾濱 150030；3.黑龍江省氣象臺，哈爾濱 150030)

0 引 言

黑龍江省1998 年嫩江、松花江流域由于夏季超強降水影響發(fā)生了超百年一遇的特大洪水, 國民經(jīng)濟直接損失達296.1 億元人民幣。2013 年全省降水異常偏多，全省夏季平均降水量為456.9 mm，比常年偏多33%，為1961 年以來第1位，導致自1998 年之后的又一次大洪水。冷春香等[1]采用EOF經(jīng)驗正交函數(shù)分析和聚類相結(jié)合的方法，對中國主汛期降水進行了8 類雨型的空間分布分析，對我國的降水預報起到了很好的作用。對降水異常的成因分析，大多采用環(huán)流特征[2-8]、ENSO的影響[9-12]和統(tǒng)計及數(shù)值模擬[13-16]等技術(shù)。本文基于降水趨勢分析、EOF展開、環(huán)流成因分析等方法，揭示出黑龍江省降水存在16 a多、少轉(zhuǎn)換周期,2012年開始進入多水期。

1 資料與方法

資料：①國家氣候中心氣候系統(tǒng)監(jiān)測預測評估網(wǎng)站提供西太平洋副熱帶高壓資料；②NCEP/NCAR 再分析資料包括北半球1961—2017年夏季500 hPa高度經(jīng)緯網(wǎng)格點5°×5°；③黑龍江省氣象信息中心提供全省62站夏季平均降水量。

方法： 自然正交函數(shù)展開(EOF)及相關(guān)系數(shù)計算[17]。

圖1 1961—2017年夏季降水變化 Fig.1 Summer precipitation variance from 1961 to 2017

2 黑龍江省夏季降水變化特點

2.1 降水年際年代際變化

1961—2017年黑龍江省夏季降水呈稍降趨勢，見圖1。由多項式擬合劃分3個長時段：1967—1980年14 a少水，平均降水301.1 mm；1981—1998年18 a多水,平均降水362.7 mm；1999—2012年14 a少水,平均312. 2 mm。2013—2017年多水,平均降水349.7 mm，已進入多水時期。降水最多年份出現(xiàn)在2013年456.9 mm,最少年份出現(xiàn)在2007年201.1 mm。多水段雨量與第1和第3段少水段雨量之差分別為61.6 mm、50.5 mm。另外1961—1966年，也為多水時段，但它屬于前期多水期的一部分。由降水等級劃分(表1)取降水量特多特少量級統(tǒng)計，在多水年時段共有29 a，其中有25 a異常偏多占多水時段的86%；而少水時段共28 a，出現(xiàn)異常偏少26 a占少水時段93%；多、少水時段的劃分接近t=0.10顯著性檢驗。多少水變化周期在33 a左右，平均半周期為16 a左右，而且多、少水年數(shù)基本持平。從年代際變化來看，上世紀80—90年代多水，70年代和本世紀初10年代少水。

表1 夏季降水等級標準

2.2 夏季降水分型

2.2.1 夏季降水EOF 展開分型

將1961—2017年全省62站夏季平均降水，進行EOF展開取前4個特征向量，它們的方差見表2。第一特征向量方差最大為0.40,第二特征向量方差為0.11，其它特征向量越來越小。全省62站降水展開的第一特征向量值、第二特征向量值、第三特征向量值、第四特征向量值分布見圖2。第一特征向量分布型式即全省降水分布一致多(或少)見圖2(a)，其降水中心位于齊齊哈爾市的克山和拜泉，哈爾濱市的阿城、五常、巴彥、木蘭、延壽、依蘭，佳木斯市和鶴崗市。這種降水形勢主要受來自南來氣旋和西移蒙古氣旋共同影響。第二特征向量分布型式即西多東少型(或西少東多型)見圖2(b)，多水主要分布在黑河以南包括齊齊哈爾、大慶、綏化所轄區(qū)，其中齊齊哈爾市的甘南、龍江、克山、拜泉降水最大，降水最少區(qū)主要分布在哈爾濱市的通河、依蘭，佳木斯市的富錦，牡丹江市的林口，七臺河市的勃利，雞西市和雙鴨山市的寶清及鶴崗市的蘿北等地。主要受來自西部的蒙古氣旋和河套氣旋以及來自西部的鄂霍次克海阻塞高壓共同影響。第三特征向量分布型式為南多北少加東部少型見圖2(c)，其主要降水在哈爾濱市的阿城、雙城、五常、尚志和齊齊哈爾市的泰來，大慶市的杜蒙、林甸。少水主要分布在大興安嶺和黑河市及轄區(qū)，降水中心在黑河市及孫吳、遜克。主要受到來自西南的河套氣旋和西北的貝加爾湖高壓的影響。第四特征向量分布型式為中多東西部少型見圖2(d)，其主要降水中心分布在伊春市，綏化市的慶安。少水主要分布在東部的完達山脈臨近區(qū)域，包括牡丹江市東部各縣(市)，雞西轄區(qū)各縣(市)還有綏芬河市。西部少水區(qū)域主要包括齊齊哈爾市全區(qū)及大慶西部的杜蒙和林甸；少水中心在齊齊哈爾市的龍江、甘南和泰來。主要受到來自南部的華北、渤海氣旋及臺風的影響。

表2 黑龍江省夏季降水EOF展開前四個特征向量方差及累計方差

圖2 1961—2017年夏季降水EOF展開前四個特征向量分布型Fig.2 Distribution types of the first four eigenvector EOF expansion on the summer precipitation from 1961 to 2017

2.2.2 EOF展開特征向量時間系數(shù)與降水變化關(guān)系

2.2.2.1 夏季降水與EOF1時間系數(shù)變化關(guān)系

全省夏季降水與EOF1時間系數(shù)變化見圖3，兩條曲線升降變化基本一致，相關(guān)系數(shù)為0.64，遠超0.01信度檢驗。降水距平符號與EOF1時間系數(shù)“+、-”變化一致率為0.74。如取降水距平百分率≥0.15或≤0.15，則與時間系數(shù)相同符號達0.83，僅有1982、1991和2009年符號關(guān)系相反，表明變化關(guān)系不在EOF1時間系數(shù)變化之內(nèi)，而是在其它降水形式之內(nèi)，如1982年夏季降水應在EOF3形勢控制之下，1991年應在EOF4形式控制之下。全省降水EOF1特征向量值分布與時間系數(shù)變化，可以反映全省夏季降水的變化情景。

2.2.2.2 夏季降水與EOF2時間系數(shù)變化關(guān)系

選取圖2(b)降水西多東少型代表站西部拜泉和東部勃利降水作與EOF2時間變化關(guān)系見圖4。1961—2017年拜泉降水與EOF2時間變化呈反位相變化關(guān)系，相關(guān)系數(shù)為-0.28通過0.05信度檢驗；勃利降水與EOF2時間系數(shù)呈正位相變化關(guān)系，相關(guān)系數(shù)為0.49通過0.001信度檢驗，東部降水與EOF2時間系數(shù)變化關(guān)系要好于西部。經(jīng)統(tǒng)計拜泉降水取“+、-”距平符號對應時間系數(shù)正確率僅為50%，取降水距平百分率≥0.20或≤0.20對應時間系數(shù)正確率達60%，取降水距平百分率≥0.40或≤0.40對應時間系數(shù)正確率達83%，隨降水的增多與減少與時間系數(shù)的反相變化的正確率在增加。勃利降水同樣取“+、-”距平符號對應時間系數(shù)正確率達70%, 取降水距平百分率≥0.20或≤0.20對應時間系數(shù)正確率達70%，取降水距平百分率≥0.40或≤0.40對應時間系數(shù)正確率達87%。不論西部或東部降水越大(小)，東西部降水反相變化大，則EOF2型分布越明顯。

圖3 黑龍江省夏季降水與EOF1時間系數(shù)變化關(guān)系 Fig.3 Relationship between the Heilongjiang summer precipitation and EOF1 time coefficient change

圖4 拜泉、勃利夏季降水與EOF2時間系數(shù)變化關(guān)系Fig.4 Relationship between the Baiquan and Boli summer precipitation and EOF2 time coefficient change

3 大氣環(huán)流變化對黑龍江省主要降水形勢的影響

3.1 大氣環(huán)流變化對全省降水EOF1型的影響

3.1.1 EOF1時間系數(shù)與500 hPa高度場相關(guān)分析

全省降水EOF1型占總方差40.39%，EOF1型時間系數(shù)與全省平均降水相關(guān)系數(shù)為0.64，因此EOF1時間系數(shù)變化類似全省降水變化。1961—2017年6—8月EOF1降水時間系數(shù)與北半球高空500 hPa高度場相關(guān)計算分析見圖5 。在東亞區(qū)500 hPa高度場中高緯地區(qū)“+、-”相關(guān)區(qū)分布明顯，其主要“-”相關(guān)區(qū)位于100～130°E，45～50°N相關(guān)系數(shù)為-0.30通過0.01信度檢驗。中心位于110～125°E，45～50N°相關(guān)系數(shù)為-0.40以上通過0.001信度檢驗相關(guān)很高。此區(qū)正是位于貝加爾湖東南部，蒙古國以東內(nèi)蒙東北及黑龍江全省?！?”相關(guān)則表明系數(shù)越低，則EOF1時間系數(shù)值越大，黑龍江省降水越多的反位相變化。同時也表明內(nèi)蒙東北及黑龍江受高空東亞低槽長時間控制，地面經(jīng)常受蒙古低壓、河套低壓及南來系統(tǒng)的影響造成多雨。位于西伯利亞地區(qū)80～100°E,65～75°N為“+”相關(guān)區(qū)，相關(guān)系數(shù)為0.30以上通過0.01信度檢驗，表明系數(shù)越高則EOF1時間系數(shù)值越大，黑龍江省降水越大的正位相變化。另外在我國東海東部海域，即西太平洋副高區(qū)域也維持正相關(guān)較高值，而且位置偏北。也即西伯利亞的泰米爾地區(qū)長時間維持阻塞高壓，我國東北區(qū)受低槽控制則黑龍江省降水偏多。從整個東亞中高緯地區(qū)高度場形勢來看，呈北高南低形勢，黑龍江省降水則盛行EOF1分布形勢，全省降水偏多或特多。

圖5 全省降水EOF1時間系數(shù)與500 hPa高度場相關(guān)系數(shù)分布(等值線為相關(guān)系數(shù)分布，陰影為通過顯著性水平α=0.05的t檢驗)Fig.5 Distribution of the relationship coefficient between provincial precipitation EOF1 time coefficient and 500 hPa height fields

3.1.2 EOF1型下多水年高空500 hPa環(huán)流形式分析

選取EOF1時間系數(shù)大又是黑龍江省夏季降水特多的1981、1985、1991、1994和2013年做高空500 hPa高度場平均合成圖6 。在中高緯地區(qū)盛行兩脊兩槽型，我國東北區(qū)與哈薩克斯坦上空為低槽控制，鄂霍茨克海和貝加爾湖為高脊控制。冷空氣沿貝加爾湖高脊前部進入東北與西南方向和副高外緣東南方向相疊加的暖濕空氣相遇交綏，前部又受到穩(wěn)定的鄂海高壓阻擋，西太平洋副高脊線偏北，西伸脊點位置偏東保證有充足的水汽供應，造成黑龍江省長時間降水。

圖6 多水年1981、1985、1991、1994、2013年夏季500 hPa平均高度場合成Fig.6 Average composition of the 500 hPa height on 1981、1985、1991、1994、2013 of much water years

圖7 全省降水EOF2時間系數(shù)與500 hPa高度場相關(guān)系數(shù)分布(等值線為相關(guān)系數(shù)分布，陰影為通過顯著性水平α=0.05的t檢驗)Fig.7 Distribution of the relationship coefficient between provincial precipitation EOF2 time coefficient and 500 hPa height fields

3.2 大氣環(huán)流變化對全省降水EOF2型的影響

3.2.1 EOF2時間系數(shù)與500 hPa高度場相關(guān)分析

全省降水EOF2型占總方差11.16%，與夏季降水相關(guān)系數(shù)為0.12,這種形式的降水在全省占第二位，主要影響的是區(qū)域降水。同樣做1961—2017年6—8月EOF2型降水時間系數(shù)與北半球高空500 hPa高度場相關(guān)計算分析見圖7 。在東亞區(qū)500 hPa高度場中高緯地區(qū)主要盛行一脊一槽型。位于鄂霍茨克海及周邊地區(qū)長期維持阻塞高壓(簡稱鄂高)中心140～150°E,50～60°N,相關(guān)系數(shù)為0.30以上。另在80°E以西，75°N以北喀拉海區(qū)域也維持一個通過信度檢驗的正相關(guān)區(qū)，即受高壓控制。而在我國東北，華北東部及內(nèi)蒙東北為低槽所控制，雖然“-”相關(guān)區(qū)未通過顯著性檢驗，但對于全省降水EOF2型分布即西多東少型分布是顯著的。另外此相關(guān)形式所見西太平洋副熱帶高壓區(qū)脊線是偏北偏西的。

3.2.2 EOF2型下多水年高空500 hPa環(huán)流形式分析

選取EOF2時間系數(shù)大又是黑龍江省降水特多的1966、1969、1984、1998和2003年做高空500 hPa高度場平均合成見圖8，其基本形式相似圖7相關(guān)系數(shù)分布。但不同的是在中高緯地區(qū)基本呈北脊南槽型，65°N及以北地區(qū)為高壓帶，貝加爾湖以東中國東北華北為高空低槽控制。此種形式對短期天氣過程來說，常有切斷低壓或東北冷渦較長期盛行維持造成多雨。而西太平洋副高脊線大部偏南西伸脊點偏西，大量的暖濕氣體由副高西部邊緣隨強大的西南氣流進入黑龍江，與貝加爾湖東南下的冷空氣相遇造成強降水。這種形式基本符合1969和1998年夏季黑龍江省西部發(fā)生強降水造成嫩江-松花江大洪水。

圖8 多水年的西多東少型1966、1969、1984、1998、2003年夏季500 hPa平均高度場合成Fig.8 Average composition of the summer 500 hPa height for west more east little type much water years in 1966、1969、1984、1998、2003

圖9 西太平洋副熱帶高壓脊線與夏季降水變化關(guān)系Fig.9 Relationship between the western pacific subtropical high ridge position and summer precipitation change

3.2.3 夏季西太平洋副熱帶高壓對黑龍江省夏季降水的影響

夏季西太平洋副熱帶高壓(以下簡稱副高)變化對我國東部地區(qū)降水影響是很大的[7]，黑龍江降水也同樣受它的影響。經(jīng)計算夏季副高面積、強度、脊線、西伸與黑龍江降水的相關(guān)系數(shù)分別為0.03、0.02、0.2、0.0,僅有脊線位置變化與降水相關(guān)一項剛好達到0.10信度檢驗標準。由副高脊線與降水關(guān)系(圖9)分析，降水距平符號“+、-”檢驗正確率達60%，降水偏多年距平百分率≥0.11共16年，其中脊線偏北有11年正確率達69%；降水特多年距平百分率≥0.15共11年，其中1961、1963、1981、1984、1985、1991、1994、2009、2013年脊線偏北正確率達82%；如降水距平百分率≥0.20有1981、1985、1991、1994、2009、2013年脊線偏北正確率達100%。降水偏多年份脊線大部偏北且西伸脊點大多偏東；1998年降水特多發(fā)生百年不遇的嫩江-松花江特大洪水，脊線偏南脊點且異常偏西；與它類似1969年也是降水偏多發(fā)生了嫩江-松花江大洪水，同樣是脊線偏南脊點且異常偏西，表明若脊線偏南脊點偏西降水中心在黑龍江省西部，盛行全省降水西多東少型。另外降水特少大部脊線偏南，如其中個別脊線偏北的，但也是絕大多數(shù)西伸脊點偏西。

4 結(jié) 論

黑龍江省夏季經(jīng)常受到來自西風帶系統(tǒng)和南來系統(tǒng)影響且水汽充足，降水量較大且規(guī)律性明顯具有多、少15～16 a轉(zhuǎn)換時段；1961—2017年降水趨勢變化不大。采用EOF方法對黑龍江省降水進行展開，分析了全省降水EOF1和EOF2型與高空500 hPa環(huán)流影響系統(tǒng)的關(guān)系，指出了位于西北部泰米爾高壓發(fā)展和在它東南部的我國東北和內(nèi)蒙為低槽控制的一脊一槽形勢和東北部鄂霍茨克海高壓控制下的一槽一脊形勢，加之西太平洋副高脊線偏北共同影響造成降水異常偏多。