田思婷,王毅勇,周超凡

(寶雞文理學(xué)院 地理與環(huán)境學(xué)院,陜西 寶雞 721013)

氣候變化已成為人類十分關(guān)注的問題,對生態(tài)系統(tǒng)[1]和社會經(jīng)濟(jì)[2]發(fā)展造成了嚴(yán)重的威脅。氣候變化引發(fā)的極端降水造成了洪澇災(zāi)害和干旱事件的頻發(fā),危害人類的生命健康[3]和財(cái)產(chǎn)安全。研究表明,未來中國區(qū)域之間的脆弱性差別很大,并且存在幾個明顯的高脆弱區(qū),主要分布在東北和西北地區(qū)(甘肅和寧夏)[4]。在我國西北干旱區(qū),極端高溫事件明顯增多,氣溫變化也存在季節(jié)和區(qū)域差異[5-7]。早期已經(jīng)有不少國外的學(xué)者開始研究氣候變化[8-9]。在極端氣溫研究方面,區(qū)域性變化是極端高溫變化的一個特征[10]。在全球陸地降水的研究中,不同尺度的降水空間變化所對應(yīng)的時間序列有很大的差別,并且具有顯著的年際特征[11]。目前,眾多學(xué)者采用不同分析方法對不同地區(qū)的降水和溫度時空變化趨勢進(jìn)行了研究[12-14]。無可置疑,水資源與氣溫已經(jīng)成為衡量全球氣候變化十分重要的因素。

渭河是黃河的最大支流,位于我國西北地區(qū)東部生態(tài)脆弱帶,自然災(zāi)害多發(fā),這也是限制這一區(qū)域的國民經(jīng)濟(jì)尤其是農(nóng)業(yè)經(jīng)濟(jì)以及農(nóng)業(yè)可持續(xù)發(fā)展的一個重要原因[15]。在氣候變暖的背景下,渭河流域降水和熱量的時空分布發(fā)生了相應(yīng)的變化,這對渭河流域生態(tài)系統(tǒng)的穩(wěn)定、人類的生活以及農(nóng)業(yè)的發(fā)展等都帶來了十分嚴(yán)重的影響。近年來,許多學(xué)者分別從氣溫及降水等角度分析了渭河流域極端氣候情況[16-18],但上述研究均是針對渭河流域整體在宏觀尺度上開展的,對渭河上游地區(qū)的研究較少。渭河上游地區(qū)大部分處于甘肅省,主要以丘陵山地為主,地貌和氣候類型復(fù)雜多樣[19];夏秋季受太平洋副熱帶高壓影響,高溫降水集中且多暴雨。有必要對渭河上游地區(qū)的氣候變化進(jìn)行較為精細(xì)化的研究,為渭河源區(qū)的生態(tài)保育與環(huán)境治理提供科學(xué)依據(jù)。

1 研究區(qū)概況

渭河全長818 km,控制流域面積約13.48×104km2。自陜西省寶雞峽以上河段為其上游河段,主要流經(jīng)甘肅省天水市、寧夏回族自治區(qū)固原市和陜西省寶雞市。渭河流域上游的地形主要以山地丘陵為主,處于我國半濕潤半干旱氣候過渡區(qū),降水集中在夏天,年平均降水量在500～800 mm之間,時間分配上極不均勻[16]。圖1為渭河上游地形圖和氣象站點(diǎn)分布圖。

圖1 渭河上游地形圖和氣象站點(diǎn)分布圖Fig. 1 Topographic map and distribution map of meteorological stations in the upper reaches of Weihe River

2 數(shù)據(jù)來源及方法

2.1 數(shù)據(jù)來源

本文數(shù)據(jù)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://data.cma.cn/)。選取渭河上游段較完整的22個氣象站點(diǎn)降水、氣溫?cái)?shù)據(jù),且訂正了不同原因造成的錯誤。數(shù)據(jù)系列包括1969-2017年期間的日降水量和日平均氣溫?cái)?shù)據(jù)。

2.2 研究方法

本文采用五年滑動平均法和線性回歸分析法分析降水和氣溫的年際變化特征;通過M-K(Mann-Kendall)突變檢驗(yàn)[20]對降水和氣溫進(jìn)行突變分析;通過小波分析法[21]對降水和氣溫進(jìn)行周期分析,通過經(jīng)驗(yàn)正交分解[22-24]和克里金插值法進(jìn)行空間分析,并繪制了渭河上游降水和氣溫變化的空間分布圖。

3 結(jié)果與分析

3.1 渭河上游降水量時空變化特征

3.1.1 渭河上游降水量時間變化特征

(1)年平均降水量的時間變化特征

運(yùn)用五年滑動平均和線性回歸方法分析渭河上游1969-2017年年均降水量的時間變化趨勢。如圖2所示,1969-2017年渭河上游年平均降水量為510.60 mm,2003年最大(686.71 mm),1997年最少(326.86 mm),二者相差359.85 mm。近49 a渭河上游地區(qū)的年均降水量呈微降趨勢,其氣候傾向率為-7.07 mm/10 a;從五年滑動平均來看,20世紀(jì)90年代存在相對較大下降-上升態(tài)勢,其他時間段降水量變化不明顯,年均降水量在1969-1994年期間呈上升趨勢,1995-2012年呈下降趨勢,UF和UB出現(xiàn)多個交點(diǎn),UF值均未超過信度值,表明全年降水變化不顯著,由Pettitt檢驗(yàn)可知年均降水量在1991年經(jīng)歷了由多變少的不顯著突變(p>0.05)。

圖2 1969-2017年渭河上游年均降水量變化及M-K檢驗(yàn)Fig. 2 Annual average precipitation change and M-K test in the upper Weihe River from 1969 to 2017

(2)四季降水量突變分析

對渭河上游地區(qū)近49 a四季降水量的時間序列進(jìn)行M-K突變檢驗(yàn)。如圖3所示,春季降水量整體呈下降趨勢(UF<0),UF值在1982年超出信度的下限,表明降水量顯著下降,由Pettitt檢驗(yàn)可知春季降水量在1992年經(jīng)歷了由多變少的不顯著突變(p>0.05);夏季降水量整體呈上升趨勢(UF>0),UF值在1980-1981年期間超出信度上限,表明該時段降水量顯著上升,由Pettitt檢驗(yàn)可知夏季降水量在1976年經(jīng)歷了由少變多的不顯著突變(p>0.05);秋季降水量整體呈現(xiàn)下降趨勢(UF<0),UF和UB出現(xiàn)多個交點(diǎn),存在多個突變點(diǎn),表明此期間震蕩劇烈,UF值在1999年超出信度的下限,表明降水量顯著下降,由Pettitt檢驗(yàn)可知秋季降水量在1986年經(jīng)歷了由多變少的不顯著突變(p>0.05);冬季降水量在1969-1990年總體呈現(xiàn)下降趨勢(UF<0),但是1991-2017年總體呈上升趨勢(UF>0),UF和UB出現(xiàn)多個交點(diǎn),存在多個突變值,UF值未超出信度值,表明此階段降水量變化不顯著,由Pettitt檢驗(yàn)可知冬季降水量在1988年經(jīng)歷了由少變多的不顯著突變(p>0.05)。

圖3 1969-2017年渭河上游四季降水量的M-K檢驗(yàn)Fig. 3 M-K test of seasonal precipitation in the upper reaches of Weihe River from 1969 to 2017

通過MATLAB R2018b小波周期分析方法得到渭河上游年際降水量小波方差及小波實(shí)部數(shù)據(jù),再利用Origin 2019b繪制小波實(shí)部圖。小波方差的峰值越大,則代表該尺度的周期性越明顯。由圖4可知,降水小波方差圖中有3個顯著峰值,以峰值的高低可以看出,第一、二、三主周期分別為55 a,35 a和15 a。55 a尺度對應(yīng)小波方差最大峰值,說明55 a尺度左右的周期性更強(qiáng),主要控制渭河上游地區(qū)降水的變化特征。對應(yīng)降水的小波實(shí)部圖可以看出,在55 a年時間尺度條件下它的周期是31 a;在35 a時間尺度條件下它的周期是23 a;在15 a時間尺度條件下它的周期是9～11 a。

圖4 降水小波方差及小波實(shí)部Fig. 4 Precipitation wavelet variance and wavelet real part

3.1.2 渭河上游降水量空間變化特征

利用ARCGIS軟件,采用普通克里金(Kriging)插值法進(jìn)行空間差值分析,生成年均和四季降水量的空間分布圖(圖5)?？梢钥闯?多年降水量整體呈現(xiàn)從東南到西北遞減的變化趨勢,太白站降雨最多,年均降水量為709.2 mm,降水最少的為通渭站,年均降水量只有385.7 mm,兩者相差323.5 mm,空間差異明顯;四季降水量的空間分布與年均降水相似,均呈現(xiàn)從東南到西北減少的趨勢。

圖5 1969-2017年年均和四季降水空間分布圖Fig. 5 Spatial distribution of annual and seasonal precipitation from 1969 to 2017

為了進(jìn)一步研究渭河上游降水量的空間分布特征,對渭河上游22個氣象站點(diǎn)1969-2017年的降水?dāng)?shù)據(jù)進(jìn)行經(jīng)驗(yàn)正交函數(shù)(Empirical Orthogonal Function, EOF)分析,并進(jìn)行North檢測,得出主要空間分布型。

1969-2017年渭河上游地區(qū)年均降水量前3個特征向量的特征值和貢獻(xiàn)率如表1所示,其前2個特征向量的累積貢獻(xiàn)率接近80%,說明這2個特征值可以很好地反映渭河上游地區(qū)49 a年均降水量的2種分布類型,如圖6所示。

表1 年降水前3個特征向量特征值及其方差貢獻(xiàn)率Tab. 1 The first three eigenvector eigenvalues and their variance contribution rates of annual precipitation

圖6 1969-2017年渭河上游年降水前2個特征向量Fig. 6 The first two eigenvectors of annual precipitation in the upper reaches of Weihe River from 1969 to 2017

第一特征向量的方差貢獻(xiàn)率為60.92%,是渭河上游的主要空間分布型。該模態(tài)特征向量的高值區(qū)出現(xiàn)在張家川站,低值出現(xiàn)在千陽站,且其高值中心在最西端,反映出該區(qū)域降水變化量大,特征向量值從東南向西北遞增,以東南為負(fù)值區(qū),以西北為正值區(qū)。該降水型存在明顯的東南或西北部的差異;第二特征向量的累積方差貢獻(xiàn)率為77.32%,該模態(tài)特征向量的高值區(qū)出現(xiàn)在天水站,低值出現(xiàn)在清水站,該降水型存在明顯的東西或南北差異。綜上可知,降水存在空間差異性。

3.2 渭河上游氣溫時空變化特征

3.2.1 渭河上游氣溫時間變化特征

(1)年平均氣溫的時間變化特征

運(yùn)用五年滑動平均和線性回歸方法分析渭河上游1969-2017年年均溫度的時間變化趨勢,如圖7所示,年均溫度整體呈上升趨勢,其氣候傾向率為0.344 ℃/10 a,年均溫度最高達(dá)9.9 ℃,最低為7.3 ℃,兩者相差2.6 ℃。UF和UB相交于1997年,存在1個突變點(diǎn),UF值在1999-2017年期間超出信度上限,表明該時間段溫度顯著上升,由Pettitt檢驗(yàn)可知在1997年發(fā)生顯著突變增溫(p<0.05)。

圖7 1969-2017年渭河上游年均溫度變化及M-K檢驗(yàn)Fig. 7 Annual average temperature changes and M-K test in the upper reaches of Weihe River from 1969 to 2017

(2)四季氣溫突變分析

對1969-2017年渭河上游地區(qū)四季氣溫的時間序列進(jìn)行M-K趨勢突變檢驗(yàn),結(jié)果如圖8所示。春季氣溫總體呈現(xiàn)先下降再上升趨勢,UF和UB相交于1999年,存在1個突變點(diǎn),UF值在2004-2017年期間超出信度的上限,表明該時段溫度顯著上升,由Pettitt檢驗(yàn)可知春季溫度在1997年發(fā)生顯著突變(p<0.05);夏季氣溫總體呈現(xiàn)先下降再上升趨勢,UF和UB分別相交于2005年、2007年和2008年,存在3個突變點(diǎn),UF值在1976-1994年期間超過信度的下限,表明該時間段溫度顯著下降,在2009-2017年期間超過信度的上限,表明該時間段溫度顯著上升,由Pettitt檢驗(yàn)可知夏季溫度在1994年發(fā)生顯著突變(p<0.05);秋季總體呈上升趨勢,UF和UB相交于1995年和2000年,存在2個突變點(diǎn),UF值在1996年和1998-2017年期間超過信度的上限,表明該時間段溫度顯著上升,由Pettitt檢驗(yàn)可知秋季溫度在1994年發(fā)生顯著突變(p<0.05);冬季UF和UB相交于1994年,存在1個突變點(diǎn),UF值在1998-2017年期間超出信度的上限,表明該時間段溫度顯著上升,由Pettitt檢驗(yàn)可知秋季溫度在1997年發(fā)生顯著突變(p<0.05)。

圖8 1969-2017年渭河上游四季溫度的M-K檢驗(yàn)Fig. 8 M-K test of seasonal temperature in the upper reaches of Weihe River from 1969 to 2017

從圖9可知,溫度小波方差圖中有1個顯著峰值,以峰值的高低可以看出,第一、二、三主周期分別為56 a,44 a和31 a。56 a尺度對應(yīng)小波方差最大峰值,說明56 a尺度左右的周期性更強(qiáng),主要控制渭河上游地區(qū)氣溫的變化特征。對應(yīng)溫度的小波實(shí)部圖可以看出,在56 a年時間尺度條件下它的周期是35 a,在44 a時間尺度條件下它的周期是27 a,在31 a時間尺度條件下它的周期是17～21 a。

圖9 溫度小波方差及小波實(shí)部Fig. 9 Temperature wavelet variance and wavelet real part

3.2.2 渭河上游氣溫空間變化特征

利用ArcMap10.5中的克里金插值得到渭河上游年際和四季氣溫的空間變化,如圖10所示?？梢钥闯?多年溫度整體呈現(xiàn)從西北到東南遞增的變化趨勢,眉縣站的年均溫度最高(13 ℃),華家?guī)X站最低(4 ℃),兩者相差9 ℃,空間存在差異性。四季的溫度空間分布與年溫度變化相似,均呈現(xiàn)從西北到東南遞增的趨勢。

為了進(jìn)一步研究渭河上游溫度的空間分布特征,對渭河上游22個氣象站點(diǎn)1969-2017年的溫度數(shù)據(jù)進(jìn)行經(jīng)驗(yàn)正交函數(shù)(EOF)分析,得出主要空間分布型。

1969-2017年渭河上游地區(qū)年均氣溫前3個特征向量的特征值和貢獻(xiàn)率如表2所示,其前2個特征向量的累積貢獻(xiàn)率達(dá)到91.91%,說明這2個特征值能夠很好地反映渭河上游地區(qū)近49 a氣溫的2種分布類型,如圖11所示。第一特征向量的方差貢獻(xiàn)率為87.09%,是渭河上游的主要空間分布型。該模態(tài)特征向量的高值區(qū)出現(xiàn)在甘谷站,低值出現(xiàn)在隆德站,以寶雞、天水、甘谷、通渭為界,東北部為負(fù)值,存在東北和西南部的差異。第二特征向量的累積方差貢獻(xiàn)率為91.91%,該模態(tài)特征向量的高值區(qū)出現(xiàn)在通渭站,低值出現(xiàn)在麟游站,以渭源、通渭、靜寧、西吉為界的西北部區(qū)域?yàn)檎?大部分區(qū)域?yàn)樨?fù)值,存在東西和南北差異性。綜上可知,溫度存在區(qū)域差異性。

表2 年均溫度前3個特征向量特征值及其方差貢獻(xiàn)率Tab. 2 The first three eigenvector eigenvalues and their variance contribution rates of annual average temperature

圖11 1969-2017年渭河上游年均溫度前2個特征向量Fig. 11 The first two eigenvectors of annual average temperature in the upper reaches of Weihe River from 1969 to 2017

4 結(jié)論

(1)從時間上看,1969-2017年渭河上游區(qū)域年均降水量無明顯變化趨勢,M-K突變分析表明,年均、春季和秋季降水量分別在1991年,1992年和1986年經(jīng)歷了由多變少的突變;夏季和冬季降水量分別在1976年和1988年經(jīng)歷了由少變多的突變。年均降水量存在55 a,35 a和15 a左右尺度的周期變化。

1969-2017年渭河上游區(qū)域年均氣溫整體呈現(xiàn)上升趨勢,M-K突變分析表明,夏季和秋季氣溫在1994年發(fā)生顯著突變,年均、春季和冬季氣溫在1997年發(fā)生顯著突變,均表現(xiàn)出增溫現(xiàn)象。年際氣溫存在56 a,44 a和31 a左右尺度的周期變化。

(2)從空間上看,年均降水量和年均氣溫在空間尺度上整體的變化均呈現(xiàn)從西北到東南遞增的規(guī)律。EOF分析表明,渭河上游年降水量和氣溫均存在明顯的空間差異性,第一特征向量降水型存在明顯的東南或西北部的差異,第二特征向量降水型存在明顯的東西或南北差異。第一特征向量氣溫型存在東北或西南部的差異;第二特征向量氣溫型存在東西或南北差異性。