包福祥

內蒙古橫跨中國東北、華北、西北三大地區(qū)，地域狹長。內蒙古秋季降水量少而且分布不勻，隨著全球氣候變化異?，F(xiàn)象頻繁，內蒙古秋季降水也受到不同程度的影響。秋季降水異常特征對內蒙古馬鈴薯、玉米、大豆等一些主要農(nóng)作物的成熟、收割及晾曬等產(chǎn)生影響。中國西北地區(qū)秋季降水日數(shù)在這1971～2010年間的變化趨勢更為顯著，降水量及降水日數(shù)整體呈現(xiàn)出增加的趨勢[1]。1961～2010年中國降水呈現(xiàn)減少趨勢，且存在區(qū)域差異，降水量減少幅度最大的是黃淮地區(qū)[2]。內蒙古夏季全區(qū)性異常多和異常少雨年主要集中出現(xiàn)在上世紀60年代和90年代以后，1961～2012年內蒙古全區(qū)少雨的站數(shù)百分率有顯著的增加變化趨勢[3]。

本文采用客觀方法對內蒙古秋季降水變化進行分區(qū)，并分區(qū)域分析各區(qū)域降水量變化特征及氣候突變特征。

1 資料與方法

本文分析的秋季降水資料選用內蒙古101個國家級地面氣象站1961～2017年秋季降水觀測資料。使用的經(jīng)驗正交函數(shù) （EOF）、旋轉經(jīng)驗正交函數(shù)（REOF）分解、相關分析、線性傾向估計、滑動平均、Mann-Kendall檢驗及滑動T檢驗等方法都是氣候統(tǒng)計診斷中應用較為普遍的方法。

2 結果與分析

2.1 用REOF方法對內蒙古秋季降水進行分區(qū)

在進行REOF分解之前，首先判斷旋轉經(jīng)驗正交函數(shù)的個數(shù)。對內蒙古秋季降水量標準化場進行EOF分解，得到前三個特征向量的累計方差貢獻達61.4%，第四特征向量的方差貢獻僅為4.6%，初步確定選取前3個正交函數(shù)作為旋轉經(jīng)驗正交函數(shù) （表1）。通過特征值對數(shù)曲線變化得到，在第三個特征值之后的直線斜率明顯變小。用North特征值誤差范圍判斷，取第五特征值時與第六特征值之差（18.8）小于North特征值誤差范圍（35.1），根據(jù)該方法判斷，從第五個特征開始經(jīng)驗正交函數(shù)已無有價值的信號。綜合上述三種方法選取前3個經(jīng)驗正交函數(shù)進行REOF分解。結果表明，第1旋轉載荷向量的高值區(qū)位于河套地區(qū)，等值線向東北逐漸減小的空間分布，全區(qū)域具有較好的一致性特征（圖1a）；第2旋轉載荷向量表現(xiàn)出“+-+”的分布特征，高值區(qū)位于赤峰市和通遼市（圖1b）；第3旋轉載荷向量表現(xiàn)出東西相反的分布特征，高值區(qū)位于呼倫貝爾市地區(qū)（圖1c）。選取載荷絕對值≥0.45作為區(qū)劃標準進行分區(qū)，同時結合歷年秋季降水分布特征，確定分區(qū)的邊界位置，將1961～2017年內蒙古秋季降水自西向東分成中西部、東南部和東北部3個區(qū)域，分區(qū)結果見圖1d（下文為敘述方便把3個區(qū)域分別稱為1、2、3區(qū)）。

表1 內蒙古秋季降水標準化場EOF分解前四個特征向量的方差貢獻和累計方差貢獻

圖1 內蒙古秋季降水REOF的第一、二、三特征量分布圖（a-c）及內蒙古秋季降水分區(qū)圖（d）

圖2 1961～2017年3個區(qū)域秋季降水量距平11年滑動平均曲線

圖31961 ～2017年1區(qū)秋季降水量Mann-Kendall檢驗

2.2 各區(qū)域降水量變化特征分析

2.2.1 各區(qū)域間降水相關性分析。為了解3個區(qū)域降水變化特征的相似性，按1961～1990年、1971～2000年、1981～2010年和1988～2017年4個時段分別分析各區(qū)域平均降水量的相關性特征，分析發(fā)現(xiàn)，1961～1990年、1971～2000年和 1981～2010年 3個時段各區(qū)域降水量相關系數(shù)均未通過0.1的顯著性水平檢驗，說明在這3個時段內各區(qū)域的秋季降水量年際變化不存在顯著的相關性特征；1988～2017年1區(qū)與3區(qū)降水量相關系數(shù)為0.048，未通過0.1的顯著性水平檢驗，而1區(qū)與2區(qū)、2區(qū)與3區(qū)降水量的相關系數(shù)分別為0.446和0.431，均大于0.349，通過0.05的顯著性水平檢驗，說明在1988～2017年1區(qū)與2區(qū)、2區(qū)與3區(qū)秋季降水量年際變化存在顯著的相關性（表2）。

表2 4個時段3個區(qū)域平均降水量相關分析結果

2.2.2 各區(qū)域秋季降水趨勢變化特征。利用線性傾向估計、相關分析和滑動平均方法分析了1961～2017年各區(qū)域降水量變化特征。分析得到，1～3區(qū)秋季降水量變化的氣候傾向率分別為 3.15、1.30、1.46mm·10a-1，相關系數(shù)分別為0.234、0.087、0.103。1區(qū)秋季降水量線性增多的趨勢較2、3區(qū)顯著。1961～2017年內蒙古秋季降水量距平11年滑動平均曲線表明，2區(qū)與3區(qū)的降水量趨勢變化較為相似，而1區(qū)與2、3區(qū)變化特征差異較大。1區(qū)秋季降水量在60年代至70年代初有增加的趨勢，70年代初至80年代末呈現(xiàn)明顯減少的趨勢，80年代末至21世紀10年代為增加的趨勢；2、3區(qū)降水量在60年代至70年代初有增加的趨勢，70年代初至80年代初略有減少趨勢，80年代初至90年代初有較為明顯增加趨勢，90年代至21世紀初期有減少的趨勢之后開始增加（圖2）。

2.2.3 氣候突變特征分析。為進一步分析各區(qū)域秋季降水量的突變特征，采用Mann-Kendall法對各區(qū)域降水量序列進行突變檢驗。結果顯示，1區(qū)降水量在20世紀70年初至21世紀初有顯著的減少趨勢 （圖3，虛直線表示0.05顯著性水平檢驗值±1.96）；3區(qū)降水量在20世紀60年代至70年代初有顯著的增加趨勢，在21世紀初減少的趨勢較為顯著（圖4）。

利用滑動T檢驗法對各區(qū)域1961～2017年降水量序列進行突變檢驗。結果顯示：1區(qū)降水沒有通過顯著性水平的突變點；2區(qū)降水量5年滑動和6年滑動T檢驗均在1967年前后有一次顯著的由少轉多過程 （圖5）；3區(qū)降水量12年滑動在1995年前后有一次顯著的由多轉少過程（圖6）。

圖4 1961～2017年3區(qū)秋季降水量Mann-Kendall檢驗

圖5 1961～2017年2區(qū)秋季降水量5年滑動和6年滑動T檢驗

圖6 1961～2017年3區(qū)秋季降水量12年滑動和16年滑動T檢驗

利用1961～2017年平均值計算1區(qū)秋季降水量序列的累計距平曲線，結果表明，在1978～2006年曲線呈下降的趨勢，2006年以后呈上升趨勢 （圖7），2006年是曲線的拐點，表明1區(qū)秋季降水量在2006年以后開始由少轉多。

綜合以上3種氣候突變檢驗方法對內蒙古3個區(qū)域秋季量的檢驗結果表明，內蒙古秋季降水沒有明顯的突變現(xiàn)象；中西部降水在1978～2006年呈減少的趨勢，2006年以后呈增多的趨勢；東部地區(qū)在20世紀60年代至70年代初有顯著的增加趨勢；東部偏南地區(qū)在1967年前后有一次顯著的由少轉多過程；東部偏北地區(qū)在1995年前后有一次顯著的由多轉少過程，從1999年前后至2008年前后呈顯著減少的趨勢。

圖7 1961～2017年內蒙古1區(qū)秋季降水量累計距平曲線

2.3 結論

2.3.1 利用REOF分解結果，結合歷年秋季降水分布特征，將1961～2017年內蒙古秋季降水分成中西部、東南部和東北部3個區(qū)域。

2.3.2 分析了 1961～1990 年、1971～2000 年、1981～2010年和1988～2017年4個時段內3區(qū)域降水相關性，只有1區(qū)與2區(qū)、2區(qū)與3區(qū)在1988～2017年秋季降水的相關系數(shù)通過0.05的顯著性水平檢驗。

2.3.3 3個區(qū)域的秋季降水變化趨勢差異較大，各區(qū)域降水量的突變特征也不一致。內蒙古秋季降水沒有明顯的突變現(xiàn)象，各區(qū)域均有顯著的降水趨勢轉折過程。