1956-2017年?yáng)|北地區(qū)氣溫和降水的時(shí)空變化特征

吳金華, 盛芝露, 杜加強(qiáng), 張楊成思, 張 靜,2

(1.中國(guó)環(huán)境科學(xué)研究院, 北京 100012; 2.蘭州大學(xué) 生命科學(xué)學(xué)院, 蘭州 730030)

氣候變化是全球變化研究的核心問(wèn)題和重要內(nèi)容，同時(shí)也是當(dāng)前人類社會(huì)面臨的最嚴(yán)重的環(huán)境問(wèn)題之一。聯(lián)合國(guó)政府間氣候變化專門委員會(huì)(Intergovernmental Panel on Climate Change，IPCC)最新發(fā)布的第五次評(píng)估報(bào)告[1-2]指出，近百年來(lái)全球氣候變暖毋庸置疑，1880—2012年，全球地表平均溫度大約升高了0.85 ℃，1951—2012年期間全球平均地表溫度的升溫速率為0.12 ℃/10 a，期間陸地比海洋增溫快，高緯度地區(qū)增溫比中緯度地區(qū)大，氣候突變事件發(fā)生的頻率和強(qiáng)度伴隨著全球氣溫升高可能會(huì)有所增強(qiáng)。在此背景下，不少研究表明中國(guó)的氣候變暖趨勢(shì)與全球基本一致[3-6]，中國(guó)1951—2018年增溫速率為0.24℃/10 a，明顯高于同期全球或北半球同期平均水平，中國(guó)北方地區(qū)增溫速率明顯大于南方，其中東北地區(qū)是增溫最顯著的地區(qū)之一[4]；中國(guó)極端強(qiáng)降水事件呈增多趨勢(shì)，東北地區(qū)降水量年際波動(dòng)增大[6]。

不少學(xué)者曾針對(duì)東北地區(qū)的氣候變化展開(kāi)研究[7-10]，這些研究對(duì)于正確認(rèn)識(shí)東北地區(qū)氣候變化特征有重要意義。已有研究表明，東北地區(qū)氣候在過(guò)去數(shù)幾十年表現(xiàn)為氣溫明顯升高，降水有所減少。然而綜合分析發(fā)現(xiàn)以往研究對(duì)于氣候突變多采用單一檢測(cè)方法，缺少對(duì)檢測(cè)結(jié)果的驗(yàn)證；其次是將氣候要素時(shí)間序列作為整體研究，缺少了對(duì)序列內(nèi)部階段性變化特征的挖掘；此外，以往研究所選取的氣象站點(diǎn)數(shù)和研究起止年份不盡一致，終止年份多為2000年前后，未覆蓋近二十年的氣候變化。為此，本文選取東北地區(qū)作為研究區(qū)，基于1956—2017年的月平均氣溫與月降水量數(shù)據(jù)，運(yùn)用Kriging插值及一系列時(shí)空分析方法，旨在解決以下3個(gè)問(wèn)題：(1) 現(xiàn)有研究得出的東北地區(qū)氣溫、降水的年際變化趨勢(shì)在2000年后是否仍然持續(xù)；(2) 東北地區(qū)氣溫、降水是否存在突變及階段變化特征；(3) 東北地區(qū)氣溫、降水變化的空間分布及階段變化顯著的區(qū)域。

1 材料與方法

1.1 研究區(qū)概況

東北地區(qū)位于我國(guó)大陸的東北部，疆土遼闊，是我國(guó)重要的地理、文化和經(jīng)濟(jì)分區(qū)，地理位置在115°31′—133°5′E，38°43′—53°34′N，包括了黑龍江省、吉林省、遼寧省的全部區(qū)域以及內(nèi)蒙古自治區(qū)東部的呼倫貝爾市、興安盟、通遼市、赤峰市。區(qū)域內(nèi)地形多樣，南臨渤、黃二海，東側(cè)和北側(cè)是鴨綠江、圖們江、烏蘇里江和黑龍江等流域低地，西側(cè)是呼倫貝爾高原及內(nèi)蒙古高原，內(nèi)側(cè)是由大興安嶺、小興安嶺和長(zhǎng)白山系形成的馬蹄形環(huán)繞的山地丘陵，中心部分則為三江平原、松嫩平原和遼河平原所組成的東北平原。東北地區(qū)自南向北跨越了中溫帶和寒溫帶，氣候類型為溫帶季風(fēng)氣候以及大興安嶺以西的溫帶大陸性氣候，冬季受大陸氣團(tuán)影響，寒冷干燥，盛行西北風(fēng)，夏季受海洋氣團(tuán)影響，暖熱多雨。自東南而西北，年降水量自1 000 mm降至300 mm以下，從濕潤(rùn)區(qū)、半濕潤(rùn)區(qū)過(guò)渡到半干旱區(qū)，其中夏季降雨量占了全年降雨量的絕大部分。

1.2 數(shù)據(jù)來(lái)源及處理

氣溫和降水是表征區(qū)域氣候特征的兩個(gè)最基本的氣候要素。本研究中采用中國(guó)氣象數(shù)據(jù)網(wǎng)所提供的東北地區(qū)月平均氣溫和月降水量的氣象站點(diǎn)觀測(cè)數(shù)據(jù)，剔除數(shù)據(jù)缺失嚴(yán)重的站點(diǎn)，選用區(qū)域內(nèi)共計(jì)102個(gè)站點(diǎn)數(shù)據(jù)。這些站點(diǎn)在區(qū)域內(nèi)均勻分布，可以體現(xiàn)區(qū)域的整體情況(圖1)。采用空間插值方法來(lái)推算站點(diǎn)以外區(qū)域的氣候要素情況，用插值形成的柵格數(shù)據(jù)來(lái)描述氣候要素的空間分布，能有效地提高氣象資料的利用水平。地統(tǒng)計(jì)學(xué)的Kriging插值方法在我國(guó)降雨量和氣溫研究上有很好的擬合優(yōu)度[11-13]，因此，本文也采用該插值方法生成1956—2017年期間的月平均氣溫和月降水量柵格數(shù)據(jù)，空間分辨率為1 km。

圖1 東北地區(qū)地理概況及氣象站點(diǎn)分布

1.3 研究方法

1.3.1 趨勢(shì)分析法 采用趨勢(shì)分析法擬合氣候要素的變化趨勢(shì)，通過(guò)計(jì)算氣候傾向率和累積距平值[14]定量分析氣候要素的變化特征。

1.3.2 氣候變率 氣候變率可以表示氣候變化幅度的大小[15]，并且氣候變率與災(zāi)害性天氣的發(fā)生頻率有關(guān)。絕對(duì)變率Va的計(jì)算公式為：

1.3.3 突變檢測(cè) 氣候突變是普遍存在于氣候系統(tǒng)的一種重要現(xiàn)象，是氣候從一種穩(wěn)定態(tài)跳躍式地轉(zhuǎn)變到另一種穩(wěn)定態(tài)的現(xiàn)象，表現(xiàn)為氣候在時(shí)空上從一個(gè)統(tǒng)計(jì)特性到另一個(gè)統(tǒng)計(jì)特性的急劇變化[16]。檢驗(yàn)突變情況的方法有很多種，本研究采用Mann-Kendall趨勢(shì)檢驗(yàn)法、滑動(dòng)t-檢驗(yàn)兩種檢驗(yàn)方法對(duì)氣候要素時(shí)間序列的突變情況進(jìn)行分析，以達(dá)到互相比較、互相檢驗(yàn)的效果，兩種檢驗(yàn)方法均設(shè)置0.05顯著性水平的信度線。

Mann-Kendall法是由Mann和Kendall提出的，最初用于檢測(cè)序列的一種變化趨勢(shì)，后被發(fā)展用來(lái)檢測(cè)氣候突變。該方法是世界氣象組織(WMO)推薦的一種在世界范圍內(nèi)廣泛使用的非參數(shù)檢驗(yàn)法，其優(yōu)點(diǎn)是不需要樣本遵從一定的分布，且不受少數(shù)異常值的干擾，計(jì)算簡(jiǎn)單，被許多學(xué)者用來(lái)分析氣溫、降水量、徑流和泥沙等要素時(shí)間序列的突變分析[17-19]。

滑動(dòng)t-檢驗(yàn)用來(lái)檢驗(yàn)兩組樣本平均值的差異是否顯著來(lái)檢驗(yàn)突變[20-21]，為此將連續(xù)序列分成兩段子序列，如果兩段子序列的均值差異超過(guò)了一定的顯著性水平，可以認(rèn)為均值發(fā)生了質(zhì)變，有突變發(fā)生。具體應(yīng)用中，為避免任意選擇子序列長(zhǎng)度造成突變點(diǎn)的漂移，可以反復(fù)調(diào)整子序列，提高計(jì)算結(jié)果的可靠性，在本研究中子序列長(zhǎng)度分別取5 a，10 a和15 a。

1.3.4 Fisher最優(yōu)分割法 Fisher最優(yōu)分割法是一種對(duì)有序樣本進(jìn)行聚類分析的方法，可用于對(duì)時(shí)間序列進(jìn)行分段以確定不同時(shí)段的特征[22-23]，較多應(yīng)用于地震分期[24-25]和汛期分期[26-27]。該方法的原理是在對(duì)有序樣本進(jìn)行分割時(shí)，總希望段內(nèi)差異愈小，一般采用離差平方和來(lái)表示段內(nèi)數(shù)據(jù)的變化程度，在分割法中簡(jiǎn)稱為變差，變差愈小則表明各段數(shù)據(jù)愈接近。當(dāng)分段數(shù)確定，用分割后的各段變差之和，即誤差函數(shù)作為比較不同分割的優(yōu)劣標(biāo)準(zhǔn)，以誤差函數(shù)值最小為最優(yōu)。

2 結(jié)果與分析

2.1 氣溫和降水年際變化趨勢(shì)特征

由圖2可知，1956—2017年?yáng)|北地區(qū)年平均氣溫在波動(dòng)中顯著上升，傾向率為0.327℃/10 a(p<0.001)。62 a的年平均氣溫為2.97℃，年平均溫度最高值為4.75℃，出現(xiàn)在2007年，年平均溫度最低值為0.91℃，出現(xiàn)在1969年，最高值與最低值間相差3.84℃。由圖3可知，年均溫累積距平值呈現(xiàn)先下降再上升的過(guò)程，即年均溫距平值在1987年前幾乎持續(xù)為負(fù)，1987年后距平值持續(xù)為正，表明1987年前后分別是氣溫相對(duì)較低和較高的兩個(gè)時(shí)期。

圖2 東北地區(qū)年平均氣溫年際變化趨勢(shì)

圖3 東北地區(qū)年平均氣溫累積距平值

由圖4可知，1956—2017年?yáng)|北地區(qū)的降水量波動(dòng)性較大，呈小幅減少趨勢(shì)，傾向率為-3.634 mm/10 a(p>0.05)。62 a平均降水量為499.52 mm，2013年年降水量最高，為649.02 mm，2001年年降水量最低，為377.66 mm，最高值與最低值間相差271.36 mm。由圖5可知，年降水量累積距平值有多個(gè)升降過(guò)程，降水量累積距平峰值出現(xiàn)在1964年和1998年，谷值出現(xiàn)在1982年和2009—2011年期間。峰值往往表明對(duì)應(yīng)年份前面一個(gè)時(shí)期降水較多而往后降水較少，谷值則相反。

圖4 東北地區(qū)年降水量年際變化趨勢(shì)

圖5 東北地區(qū)年降水量累積距平值

2.2 氣溫和降水序列突變檢測(cè)

年平均氣溫Mann-Kendall趨勢(shì)檢驗(yàn)法和滑動(dòng)t檢驗(yàn)突變檢測(cè)的結(jié)果見(jiàn)圖6。UF曲線整體位于0值以上，說(shuō)明氣溫時(shí)間序列呈上升趨勢(shì)。這種上升趨勢(shì)在1959年和1989年以后超過(guò)了0.05的置信水平，達(dá)到了顯著上升趨勢(shì)。UF曲線與UB曲線有一個(gè)交點(diǎn)，且交點(diǎn)在置信區(qū)間內(nèi)，說(shuō)明氣溫突變可能發(fā)生在1986—1987年。為了檢驗(yàn)結(jié)果的可靠性，本研究選用了5 a，10 a，15 a滑動(dòng)步長(zhǎng)進(jìn)行滑動(dòng)t-檢驗(yàn)，從滑動(dòng)檢驗(yàn)的結(jié)果來(lái)看，1987年附近的確出現(xiàn)了一次增溫的突變過(guò)程，從5 a滑動(dòng)檢驗(yàn)的圖中可以看出在2009年氣溫在顯著升高的過(guò)程還有一個(gè)變冷的突變。

圖6 東北地區(qū)年平均氣溫Mann-Kendall檢驗(yàn)和滑動(dòng)t檢驗(yàn)

東北地區(qū)降水量的Mann-Kendall趨勢(shì)檢驗(yàn)的結(jié)果見(jiàn)圖7，UF曲線整體位于0值以下，說(shuō)明區(qū)域內(nèi)的降水量在時(shí)間段內(nèi)呈下降趨勢(shì)，這種下降趨勢(shì)在1967—1985年期間超出了0.05的置信水平，達(dá)到了顯著下降的趨勢(shì)。UF曲線與UB曲線有4個(gè)交點(diǎn)，分別是：1957年、1959年、1961年、2016年，突變檢測(cè)效果不理想，結(jié)合滑動(dòng)t-檢驗(yàn)的結(jié)果來(lái)看，1964年左右存在一個(gè)的降水量減少的突變過(guò)程；在1983年降水量是一個(gè)增加的突變過(guò)程；之后降水量下降趨緩，在1998—1999年期間，降水量又有一個(gè)減少的突變。

圖7 東北地區(qū)年降水量Mann-Kendall檢驗(yàn)和滑動(dòng)t檢驗(yàn)

2.3 氣溫和降水階段性變化特征

結(jié)合東北地區(qū)年平均氣溫變化的特點(diǎn)，運(yùn)用Fisher最優(yōu)分割法將氣溫時(shí)間序列分割為：1956—1957年、1958—1987年、1988—2013年、2014—2017年4個(gè)時(shí)間段，表1為東北地區(qū)4個(gè)時(shí)間段年平均氣溫及其絕對(duì)變率和相對(duì)變率。各個(gè)時(shí)段的年平均氣溫呈上升趨勢(shì)，從1956—1957年的1.26℃升高到2014—2017年的3.96℃；從氣候變率上來(lái)看，1958—1987年年均溫的絕對(duì)變率和相對(duì)變率較大，氣候變率大代表著更大的氣候要素波動(dòng)，往往與異常天氣的頻率及強(qiáng)度有關(guān)。

表1 1956-2017年?yáng)|北地區(qū)不同時(shí)段年平均氣溫及其變率

結(jié)合東北地區(qū)年降水量變化的特點(diǎn)，運(yùn)用Fisher最優(yōu)分割法將降水量時(shí)間序列分割為：1956—1964年、1965—1982年、1983—1998年、1999—2009年、2010—2017年5個(gè)時(shí)間段，分割結(jié)果基本與累積距平和突變檢測(cè)的分析結(jié)果基本一致，表2為東北地區(qū)5個(gè)時(shí)段年降水量及其絕對(duì)變率和相對(duì)變率。從中可以看出，1956—1964年是年平均降水量最多的時(shí)期，降水量為543.26 mm，降水量最少的時(shí)期是1999—2009年，為444 mm。62 a間，降水量經(jīng)歷了一個(gè)周期波動(dòng)過(guò)程。2010—2017年年降水量的絕對(duì)變率和相對(duì)變率較大，代表這一時(shí)間段內(nèi)的降水情況不穩(wěn)定。降水的多變常常預(yù)示旱澇災(zāi)害發(fā)生的頻率有所增加。

查證該時(shí)期記錄旱澇情況的文獻(xiàn)資料，2011年、2012年?yáng)|北地區(qū)是明顯的干旱區(qū)之一[28-29]，2013年?yáng)|北地區(qū)的洪澇災(zāi)害嚴(yán)重，區(qū)域降水量較常年偏多10%～20%，黑龍江、松花江流域汛期平均降水量較常年偏多3～4成[30]，2014年以來(lái)雖然存在一些季節(jié)性的洪澇情況，但東北地區(qū)整體降水量偏少，干旱缺水情況在部分地區(qū)一直存在[31-34]。

表2 1956-2017年?yáng)|北地區(qū)不同時(shí)段年降水量及其變率

2.4 氣溫和降水變化的空間分布

計(jì)算1956—2017年?yáng)|北地區(qū)年平均氣溫的空間分布見(jiàn)圖8A所示，可以看出氣溫的緯度地帶性分布特征明顯，由南到北呈現(xiàn)遞減的趨勢(shì)，氣溫較高的區(qū)域?yàn)椴澈?、黃海海灣地區(qū)，年均溫在10℃左右；氣溫較低的為大興安嶺北部，年均溫低于-4℃；東北平原的氣溫相較于整個(gè)區(qū)域比較適中，年平均氣溫整體大于2 ℃。計(jì)算62 a區(qū)域的氣溫傾向率，結(jié)果見(jiàn)圖8B所示，可以發(fā)現(xiàn)整個(gè)區(qū)域的傾向率均為正值，說(shuō)明1956—2017年區(qū)域整體的氣溫都呈上升趨勢(shì)。氣溫傾向率較高的區(qū)域在小興安嶺區(qū)域，其次是呼倫貝爾高原、內(nèi)蒙古高原以及東北平原的部分地區(qū)，氣溫傾向率較低的區(qū)域主要在遼河流域。

基于Fisher最優(yōu)分割法分割年平均氣溫時(shí)間序列的結(jié)果，分別計(jì)算1956—1957年、1958—1987年、1988—2013年、2014—2017年4個(gè)時(shí)間段(下文將分別用氣溫時(shí)間分段1，2，3，4進(jìn)行替代)年平均氣溫的平均值，將相鄰時(shí)間段的平均氣溫柵格相減，得到時(shí)間段間的氣溫變化差值。分析結(jié)果可以得出：時(shí)間分段2相較于分段1，除大興安嶺北部區(qū)域氣溫下降之外區(qū)域整體氣溫上升，其中氣溫增加最多的區(qū)域是在呼倫貝爾高原區(qū)域，內(nèi)陸的增溫幅度要高于沿海地區(qū)；時(shí)間分段3相較于分段2，區(qū)域氣溫也是有所增加，增溫幅度較大的區(qū)域變成了小興安嶺地區(qū)，其次才是呼倫貝爾高原，增溫較小的區(qū)域集中在遼河流域以及區(qū)域東北沿海地區(qū)；時(shí)間分段4相較于分段3，大興安嶺和遼河平原的氣溫有所上升，而小興安嶺和三江平原部分區(qū)域氣溫有所下降；時(shí)段間的增溫幅度有所減緩。

圖8 1956-2017年年平均氣溫均值和傾向率的空間分布

計(jì)算1956—2017年?yáng)|北地區(qū)年降水量的空間分布見(jiàn)附圖12A所示。從圖中可以看出，區(qū)域內(nèi)降水量最多的區(qū)域在長(zhǎng)白山區(qū)域南部臨黃海的區(qū)域，平均年降水量超過(guò)800 mm；降水量最少的區(qū)域是在呼倫貝爾高原，平均年降水量少于350 mm，其次便是大興安嶺山脈和內(nèi)蒙古高原，平均年降水量不超過(guò)500 mm，在少雨的地區(qū)降水量呈現(xiàn)由東向西逐漸遞減的趨勢(shì)。附圖12B是區(qū)域62 a降水量?jī)A向率的空間分布，結(jié)果表明，區(qū)域內(nèi)降水量增長(zhǎng)和下降趨勢(shì)同時(shí)存在，大興安嶺北部呈現(xiàn)增長(zhǎng)趨勢(shì)，其余大部分區(qū)域的降水量都呈下降趨勢(shì)，包括遼東半島、遼河平原、三江平原、內(nèi)蒙古高原東部以及呼倫貝爾高原。

由Fisher最優(yōu)分割法分割年降水量時(shí)間序列的結(jié)果，分別計(jì)算1956—1964年、1965—1982年、1983—1998年、1999—2009年、2010—2017年5個(gè)時(shí)間段(下文中將分別用降水時(shí)間分段1，2，3，4，5進(jìn)行替代)年降水量的平均值，將相鄰時(shí)間段的平均降水量柵格相減，得到時(shí)段間的降水量變化差值。結(jié)果表明，區(qū)域整體的降水量均值經(jīng)歷了一個(gè)降低—升高—降低—升高的周期波動(dòng)過(guò)程，在此過(guò)程中不同區(qū)域呈現(xiàn)出不同的變化特征。降水時(shí)間分段2與分段1相比，區(qū)域整體降水量降低，降幅較大的區(qū)域是長(zhǎng)白山南部和三江平原，降幅較小的區(qū)域是大興安嶺山脈；降水時(shí)間分段3相較分段2，降水量有所升高，主要的升高區(qū)域在松嫩平原北部靠近大興安嶺的區(qū)域；降水時(shí)間分段4相較于分段3，降水量降低，東北地區(qū)西南部，包括大興安嶺、呼倫貝爾高原、內(nèi)蒙古高原以及遼河平原，降水量減少明顯；降水時(shí)間分段5與分段4相比，降水量整體上升，其中大興安嶺南部以及長(zhǎng)白山脈西側(cè)和南側(cè)的降水量有較大提升。由降水量的變化過(guò)程可看出，2010年后，降水量進(jìn)入相對(duì)豐水期。

3 討 論

在以往的氣候變化研究中，大多數(shù)學(xué)者計(jì)算區(qū)域平均溫度或降水量多用區(qū)域內(nèi)站點(diǎn)觀測(cè)數(shù)據(jù)的均值來(lái)表示[7-9]，用算數(shù)平均作為區(qū)域氣象要素平均情況的前提是各氣象站點(diǎn)的分布是獨(dú)立的，但實(shí)際情況是站點(diǎn)間氣象要素值存在空間自相關(guān)，故將觀測(cè)數(shù)據(jù)直接用于周邊地區(qū)均值的推算不是最優(yōu)的方法。另外一些學(xué)者提出城市化的影響是造成氣候變化分析中結(jié)果不確定性的重要因素，在研究氣溫和降水變化特征時(shí)應(yīng)該排除熱島效應(yīng)，否則會(huì)出現(xiàn)虛假的增暖和降水格局[35]，這也是一些研究計(jì)算出的氣候傾向率偏大的原因。如何區(qū)分長(zhǎng)期氣候變化中的自然氣候變換項(xiàng)和人類活動(dòng)作用項(xiàng)是一項(xiàng)關(guān)鍵的問(wèn)題，有學(xué)者針對(duì)氣候變化的城市化效應(yīng)進(jìn)行了定量評(píng)價(jià)[36-37]，但剔除氣候變化中的人類活動(dòng)作用目前還沒(méi)有一個(gè)公認(rèn)有效的辦法。本研究對(duì)氣象觀測(cè)資料進(jìn)行Kriging插值，考慮氣象要素的空間自相關(guān)，計(jì)算得到的區(qū)域均值更加合理，另外插值避免將城市氣象站點(diǎn)觀測(cè)值直接用作周邊地區(qū)均值，減小了城市化效應(yīng)對(duì)定量分析氣候變化情況的影響。

本研究運(yùn)用Fisher最優(yōu)分割法對(duì)長(zhǎng)時(shí)間序列進(jìn)行分段處理和分析，有效地彌補(bǔ)了將時(shí)間序列作為整體分析的缺陷，即在對(duì)完整時(shí)間序列分析變化趨勢(shì)時(shí)，會(huì)丟失序列內(nèi)部相對(duì)短期的波動(dòng)情況以及階段變化信息。本研究將Fisher最優(yōu)分割法與變率分析相結(jié)合，檢測(cè)到2010—2017年為降水相對(duì)波動(dòng)較大的時(shí)期，旱澇發(fā)生頻繁；將Fisher最優(yōu)分割法與空間分析方法相結(jié)合，分析得出不同地區(qū)的階段升溫情況以及時(shí)段間降水量呈現(xiàn)降低—升高—降低—升高的周期波動(dòng)過(guò)程。Fisher最優(yōu)分割與常用的時(shí)空分析方法結(jié)合對(duì)于研究一個(gè)長(zhǎng)時(shí)間序列中的局部特征提供一個(gè)新的思路。

東北地區(qū)耕地資源集中，是國(guó)家重要的“糧倉(cāng)”之一，氣候變暖趨勢(shì)和降水格局的改變正逐步改變區(qū)域農(nóng)業(yè)種植結(jié)構(gòu)。典型的例子就是在東北松嫩—三江平原，中晚熟水稻品種種植面積大幅擴(kuò)張，顯著提高了糧食產(chǎn)量[38]。隨著區(qū)域水土資源開(kāi)發(fā)強(qiáng)度不斷增大，社會(huì)經(jīng)濟(jì)迅速發(fā)展，水資源供需矛盾尤為突出，這對(duì)水資源利用效率和配置情況提出了較高的要求[39]。氣候變化還對(duì)東北地區(qū)的生態(tài)環(huán)境產(chǎn)生了重大影響。東北地區(qū)是我國(guó)濕地資源的主要分布區(qū)之一，氣候變暖以及水分條件的不穩(wěn)定導(dǎo)致了濕地生態(tài)系統(tǒng)的嚴(yán)重退化[40]，同時(shí)氣候變化也是東北平原西部地區(qū)的土地荒漠化的重要影響因素[41]。因此，正確認(rèn)識(shí)氣候變化特征對(duì)于東北地區(qū)生產(chǎn)生活的開(kāi)展和生態(tài)環(huán)境保護(hù)具有十分重要的意義。從本研究的結(jié)果來(lái)看，過(guò)去幾十年內(nèi)氣溫顯著升高以及降水的波動(dòng)減少使得東北地區(qū)的氣候呈現(xiàn)暖干化的趨勢(shì)，這一結(jié)論與其他學(xué)者的研究結(jié)果相一致[7-8]。從整個(gè)時(shí)間序列的氣溫變化情況來(lái)看，氣溫上升的趨勢(shì)有所減緩；另一方面，雖然東北地區(qū)的降水呈現(xiàn)減少的趨勢(shì)，但相比21世紀(jì)初的十年間，當(dāng)前的降水量處在相對(duì)豐水期，同時(shí)水分條件不穩(wěn)定，極端降水情況發(fā)生頻繁，需要引起足夠重視。建議通過(guò)水利工程措施和水源涵養(yǎng)、水分調(diào)節(jié)等生態(tài)服務(wù)功能提升措施，提高區(qū)域水資源利用效率，保證社會(huì)經(jīng)濟(jì)活動(dòng)正常進(jìn)行。

4 結(jié) 論

(1) 1956—2017年?yáng)|北地區(qū)氣溫顯著上升，平均氣溫傾向率為0.327 ℃/10 a(p<0.001)，年均溫的最高值與最低值相差3.84 ℃，1987年前后分別是氣溫相對(duì)較低和較高的兩個(gè)時(shí)期。1956—2017年?yáng)|北地區(qū)降水量的波動(dòng)性較大，總體呈減少趨勢(shì)，平均傾向率為-3.634 mm/10 a(p>0.05)，最高年份與最低年份的降水量相差271.36 mm，降水量的變化包含有多個(gè)升降過(guò)程；

(2) 氣溫在1987年有一個(gè)增溫的明顯突變，1989年后氣溫增長(zhǎng)的趨勢(shì)達(dá)到顯著；氣溫時(shí)間序列分割結(jié)果表明氣溫的上升過(guò)程具有階段性，其中1958—1987年氣溫波動(dòng)較大；降水量在1964年、1983年、1999年附近降水量都有一個(gè)突變的過(guò)程；降水時(shí)間序列分割結(jié)果表明降水量呈現(xiàn)周期波動(dòng)過(guò)程，周期變化的時(shí)間點(diǎn)與突變時(shí)間基本一致，其中2010—2017年的降水量變率最高，降水情況較不穩(wěn)定。

(3) 東北地區(qū)氣溫由南至北氣溫逐漸降低；1956—2017年區(qū)域整體普遍升溫，但不同區(qū)域的增溫速率存在明顯差異，其中增溫較快是小興安嶺地區(qū)；通過(guò)計(jì)算不同時(shí)段間的氣溫差值發(fā)現(xiàn)區(qū)域內(nèi)增溫幅度較大的區(qū)域先是小興安嶺和呼倫貝爾高原，再是大興安嶺和遼河平原。多年平均降水量的空間分布表現(xiàn)為由東南向西北遞減，1956—2017年除大興安嶺北部外的絕大多數(shù)區(qū)域降水量呈下降趨勢(shì)；降水量呈現(xiàn)降低—升高—降低—升高的周期變化過(guò)程，時(shí)段間的降水量變化存在區(qū)域分異；2010年后，降水量進(jìn)入相對(duì)豐水期。