1981—2013華北平原氣候時空變化及其對植被覆蓋度的影響

阿 多,趙文吉,*,宮兆寧,張 敏,范云豹

1 首都師范大學(xué)資源環(huán)境與旅游學(xué)院, 北京 100048 2 三維信息獲取與應(yīng)用教育部重點實驗室, 北京 100048 3 資源環(huán)境與地理信息系統(tǒng)北京市重點實驗室, 北京 100048

1981—2013華北平原氣候時空變化及其對植被覆蓋度的影響

阿 多1,2,3,趙文吉1,2,3,*,宮兆寧1,2,3,張 敏1,2,3,范云豹1,2,3

1 首都師范大學(xué)資源環(huán)境與旅游學(xué)院, 北京 100048 2 三維信息獲取與應(yīng)用教育部重點實驗室, 北京 100048 3 資源環(huán)境與地理信息系統(tǒng)北京市重點實驗室, 北京 100048

基于1981—2013年華北平原氣象數(shù)據(jù),對華北平原近30a的氣候時空變化趨勢與突變情況進(jìn)行分析,并結(jié)合GIMMS-NDVI(1981—2006年)和MODIS-NDVI(2000—2013年)遙感數(shù)據(jù),探討氣候突變影響下,華北平原1981—2013年植被覆蓋度的空間分布和變化特征。結(jié)合生態(tài)學(xué)分區(qū),從不同時空尺度出發(fā),分析華北平原不同生態(tài)分區(qū)內(nèi)氣候因子與植被覆蓋度在年代際、年際和月變化的相關(guān)關(guān)系。結(jié)果發(fā)現(xiàn)：(1)1981—2013年,華北平原氣溫整體呈現(xiàn)顯著上升趨勢(0.20℃/10a,P<0.01),春季氣溫的升高、>15℃等溫線控制范圍的擴大和年均溫0℃等值線在華北平原的消失,是區(qū)域平均氣溫升高的誘因。華北平原降水整體呈現(xiàn)顯著減少趨勢(-1.75mm/10a,P>0.05),其中秋季降水量減少過快,400—600mm降水等值線控制范圍的擴大、600—800mm和800—1000mm降水等值線的范圍的縮小,共同造成區(qū)域降水量的減少。四季氣候傾向率的特征變化敏感區(qū)域主要位于北緯35°—39°之間。1991—1994年為華北平原氣候的突變時期。(2)華北平原植被覆蓋度總體呈現(xiàn)上升的趨勢,呈增加趨勢的面積占總面積的55%。人類活動不僅加速了區(qū)域植被覆蓋度的降低,也加劇了降低速率的變快。(3)總體上,研究區(qū)月尺度植被覆蓋度與氣候因子的相關(guān)性高于年尺度的值。植被覆蓋度與年降水量的偏相關(guān)性高于其與年均溫的偏相關(guān)性。年均溫對農(nóng)業(yè)生態(tài)區(qū)和森林生態(tài)區(qū)的植被覆蓋度的影響更大,草原生態(tài)區(qū)的植被覆蓋度對年降水量的依賴性更強。在月尺度上降水量對植被的影響具有時滯效應(yīng)。氣溫對草原生態(tài)區(qū)具有時滯效應(yīng),降水量對農(nóng)業(yè)生態(tài)區(qū)具有時滯效應(yīng)。(4)華北平原干熱化的氣候突變降低了植被覆蓋度的增加趨勢。從植被覆蓋度出發(fā),草原生態(tài)區(qū)對氣候突變的響應(yīng)最明顯；從變化速率角度出發(fā),農(nóng)田生態(tài)區(qū)的響應(yīng)最明顯。就整體而言,人類活動的影響力還在持續(xù)增強,且呈現(xiàn)出在退化區(qū)的作用力高于改善區(qū)的趨勢。氣候突變后,出現(xiàn)了人類活動在植被覆蓋度的改善區(qū)的相對作用力高于退化區(qū)的變化。

氣候變化；植被覆蓋度；趨勢分析；殘差分析；偏相關(guān)性

依據(jù)IPCC的研究報告,從1880到2012年,全球陸地和海洋表面的平均溫度升高了0.85℃[1]。在過去的50a里,我國的地表溫度升高了1.1℃,升高的速率約為0.22/10a,略高于全球水平[2]。氣候和陸氣間相互作用的變化,通過氣溫的升高,降水機制的改變以及季節(jié)周期性的波動,對植物的生長和分布產(chǎn)生重要影響[3]。因此,氣候變化引起的歸一化植被指數(shù)NDVI在一定范圍內(nèi)來回波動[4]。20a以上的長期研究揭示生物和自然系統(tǒng)變化與區(qū)域氣溫降水變化相關(guān)[5]。因此區(qū)域氣候特征變化和區(qū)域生態(tài)系統(tǒng)變化的內(nèi)在聯(lián)系,以及兩者間的響應(yīng)機制成為了研究熱點。

植被作為生態(tài)系統(tǒng)的重要組成部分,對諸如大氣系統(tǒng)和土壤系統(tǒng)等生態(tài)系統(tǒng)之間的能量傳輸有重要的影響作用[6]。與植被覆蓋度有正相關(guān)系數(shù)的植被指數(shù)NDVI是指示植被密度和長勢的重要指示器[7]。處在自然環(huán)境和人類活動共同影響下的植被覆蓋度,在全球氣候和碳循環(huán)系統(tǒng)中扮演重要的角色。所以,植被的變化情況揭示了自然環(huán)境和人類活動的相關(guān)影響,以及兩者之間的反饋影響。Myneny[8]和Tucker[9]研究發(fā)現(xiàn)長時間序列的植被覆蓋度呈現(xiàn)增加的趨勢,特別是在北半球中緯度地區(qū)。他們把這一趨勢歸因于全球氣候的變暖。一些研究發(fā)現(xiàn)不同地區(qū)的植被覆蓋度的增減趨勢是不一致的。Park[10]和Piao[11]研究發(fā)現(xiàn)在歐亞大陸增減趨勢是一致的。Wang[12]和Zhang[13]研究發(fā)現(xiàn)在北美洲,西北部和東北部的增減趨勢呈現(xiàn)相反趨勢。已有研究都表明,植被覆蓋度與氣候因子有內(nèi)在聯(lián)系,尤其是溫度因子和降水量因子[14-16]。研究表明在中國的西北地區(qū),植被覆蓋度的變化趨勢與降水量的變化趨勢一致,植被指數(shù)與降水量之間有顯著的正相關(guān)關(guān)系,與氣溫的相關(guān)性較低。依據(jù)以上研究發(fā)現(xiàn),植被覆蓋度的變化呈現(xiàn)明顯的區(qū)域差異性,與氣候因子關(guān)系的空間異質(zhì)性也較強。同時,氣候因子在一定的時間尺度內(nèi)會有突變現(xiàn)象發(fā)生[17]。氣候突變勢必影響動植物生長和生態(tài)結(jié)構(gòu)變化。因此,探討各生態(tài)系統(tǒng)內(nèi)植被覆蓋度對氣候變化的響應(yīng),對保護區(qū)域生物多樣性與生態(tài)安全具有重要意義。

基于植被覆蓋度變化與氣候因子之間的密切關(guān)系,以及植被覆蓋度時空變化的空間異質(zhì)性。本文選擇華北平原為研究區(qū),探討華北平原植被覆蓋度的空間分布特征和變化趨勢。并從氣候要素入手,結(jié)合生態(tài)系統(tǒng)分區(qū),從不同時空尺度出發(fā),分析了農(nóng)田生態(tài)區(qū)、草原生態(tài)區(qū)和森林生態(tài)區(qū)內(nèi)植被覆蓋度的年代際變化、年際變化和月變化對降水量因子和溫度因子的時空響應(yīng),并探討1981—2013年植被覆蓋度變化對氣候突變的響應(yīng)特征。

1 研究區(qū)地理概況及生態(tài)區(qū)劃分

華北平原位于30°00′—40°24′N、112°48′—122°45′E之間,包括河北、河南、山東、安徽、江蘇、北京及天津等5省2市的大部或部分地區(qū),總面積約40×104km2。研究區(qū)北起燕山,西沿太行山、伏牛山,南抵淮河干流及蘇北灌溉總渠,東臨黃海和渤海,除山東半島為丘陵外,大部分地區(qū)地勢平坦,平原面積約為33×104km2,耕地面積約占總面積的72%。研究區(qū)屬于溫帶大陸季風(fēng)氣候區(qū),光熱資源豐富,雨熱同季,降水量大且多集中于6—9月。冬小麥、夏玉米為主要種植農(nóng)作物,種植制度為一年兩熟[18-19]。

圖1 華北平原地理位置和氣象站點分布圖Fig.1 The distribution of the meteorological station

圖2 研究區(qū)生態(tài)區(qū)分區(qū)圖Fig.2 The distribution of the ecological area

參考Wu等[20]的研究成果,將研究區(qū)分為8個地理分區(qū)(圖2)。從南到北按照編碼順序依次為：淮南與長江中下游自然區(qū)(1)、漢中盆地自然區(qū)(2)、華北平原自然區(qū)(3)、魯中山地丘陵自然區(qū)(4)、遼東膠東山地丘陵區(qū)(5)、華北山地丘陵自然區(qū)(6)、大興安嶺南部自然區(qū)(7)和內(nèi)蒙古高原平原東部自然區(qū)(8)。為探討不同生態(tài)系統(tǒng)分區(qū)植被覆蓋度的變化特征,將已有自然分區(qū)按照生態(tài)系統(tǒng)類型合并為3個大類分區(qū)：森林生態(tài)系統(tǒng)區(qū)、草原生態(tài)系統(tǒng)區(qū)和農(nóng)業(yè)生態(tài)系統(tǒng)區(qū)。結(jié)果如圖2所示。草原生態(tài)區(qū)主要分布于河北省的西北部,包括張北地區(qū)和木蘭圍場。該生態(tài)區(qū)地處我國溫帶半濕潤地區(qū)與半干旱地區(qū)的交織帶,主要以溫性草甸草原、山地草甸和溫性草原為主[21]。森林生態(tài)區(qū)包括燕山山脈、太行山-伏牛山山脈和山東半島的丘陵山地,植被類型從山麓至山頂,大體分布落葉林、闊葉針葉混交林、針葉林和高山灌木草叢,主要以松林、側(cè)柏林、櫟林為主[22]。農(nóng)田生態(tài)區(qū)分布區(qū)域為除了草原生態(tài)區(qū)和森林生態(tài)區(qū)以外的區(qū)域。華北平原的農(nóng)業(yè)以一年兩熟為主,農(nóng)作物中冬小麥、玉米占絕對優(yōu)勢。

2 數(shù)據(jù)準(zhǔn)備與研究方法

2.1 數(shù)據(jù)準(zhǔn)備與一致性檢驗

美國馬里蘭大學(xué)提供的GIMMS-NDVI數(shù)據(jù)集包含1981—2006年半月最大合成數(shù)據(jù),空間分辨率為8km×8km,數(shù)據(jù)經(jīng)過輻射校正和幾何校正、大氣校正、除云等處理,GIMMS數(shù)據(jù)集的質(zhì)量優(yōu)于其他NDVI數(shù)據(jù)集[23]。GIMMS數(shù)據(jù)集是目前最長時間序列的NDVI數(shù)據(jù),與其他高分辨率數(shù)據(jù)集之間具有較好的相關(guān)性,已廣泛應(yīng)用在全球和區(qū)域尺度植被動力學(xué)研究中[24-25]。NASA提供的MODIS-NDVI數(shù)據(jù)集包含2000—2013年華北平原MOD13Q1級產(chǎn)品,時間分辨率為16d,空間分辨率為250m×250m。應(yīng)用MODIS Reprojection Tools (MRT)軟件對數(shù)據(jù)進(jìn)行投影轉(zhuǎn)換、數(shù)據(jù)拼接等處理。

在獲取植被覆蓋度的方法中,最大值合成法MVC(Maximum Value Composites)使用較為廣泛,它是基于NDVI旬?dāng)?shù)據(jù)合成獲得月度和年度NDVI數(shù)據(jù)的方法。具體公式如下[26]：

NDVIi=Max(NDVIij)

式中，NDVIi指第i個月或者第i年的NDVI, NDVIij指第i月的第j個15天的NDVI數(shù)據(jù)或第i年的第j月的NDVI 數(shù)據(jù)。

由于GIMMS和MODIS數(shù)據(jù)采用不同的傳感器,因此需要對兩者進(jìn)行一致性檢驗[27-28]。本文采用的GIMMS-NDVI數(shù)據(jù)的時間序列為1981—2006,MODIS-NDVI數(shù)據(jù)的時間序列為2000—2013年,依據(jù)2000—2006年的重疊數(shù)據(jù),對年最大NDVI數(shù)據(jù)進(jìn)行相關(guān)性分析,得到兩者相關(guān)系數(shù)為0.801,在0.05置信水平上顯著相關(guān),說明兩種數(shù)據(jù)的年度數(shù)據(jù)在全區(qū)域尺度上具有顯著一致性,二者的線性回歸方程為：NDVIGIMMS=0.5886×NDVIMODIS+0.0965(r2=0.6417,n=7,P<0.05)。采用同樣方法對兩者的月尺度數(shù)據(jù)進(jìn)行一致性檢驗,兩者月尺度數(shù)據(jù)的相關(guān)系數(shù)為0.957,在0.01置信水平上通過檢驗,線性回歸方程為：NDVIGIMMS=0.781NDVIMODIS-0.0469(r2=0.9161,n=7,P<0.01)。通過以上回歸方程,結(jié)合2007—2013年MODIS-NDVI數(shù)據(jù)插補2007—2013年度和月度GIMMS-NDVI數(shù)據(jù),將實驗數(shù)據(jù)的時間跨度延伸為1981—2013年。

2.2 突變檢測

氣候突變是普遍存在于氣候系統(tǒng)中的重要現(xiàn)象。氣候突變的普適性定義：氣候從一種穩(wěn)定態(tài)(或穩(wěn)定持續(xù)的變化趨勢)跳躍式地轉(zhuǎn)變到另一種穩(wěn)定態(tài)(或穩(wěn)定持續(xù)的變化趨勢)的現(xiàn)象,表現(xiàn)為氣候在時空上從一個統(tǒng)計特性到另一個統(tǒng)計特性的急劇變化。在氣候突變的檢測方法中,Mann-Kendall非參數(shù)檢測法的檢測范圍寬、人為性少,定量化程度高,不需要樣本遵從一定的分布,也不受少數(shù)異常值的干擾[29-30]。

Mann-Kendall檢測的基本原理如下：氣候序列設(shè)為平穩(wěn)序列(i=1,2,3,…,n),表示第i個樣本(1≤j≤i)的累計數(shù),定義統(tǒng)計量:

在原序列的隨機獨立等的假設(shè)下,均值、方差分別為：

var(dk)=k(k-1)(2k+5)/72

dk標(biāo)準(zhǔn)化后,

滿足標(biāo)準(zhǔn)正態(tài)分布,其概率a=prob(|U|>|UFk|)可以通過計算或者查表獲得。給出顯著性水平a(取a=0.05,=1.96),若|UF|>Ua則表明該氣候序列呈現(xiàn)顯著的變化趨勢。

將氣候序列(i=1,2,3,…,n)的逆序列,按照以上過程計算一次。同時使|UBk|=-UFk,k=n、n-1,...,1,UB1=0。若|UF|>0,此時序列呈上升趨勢；若|UF|<0序列呈下降趨勢。通過信度檢驗可知其是否有變化趨勢。當(dāng)曲線超過置信水平線時,表明序列的上升(下降)趨勢顯著.若UF和UB兩條曲線出現(xiàn)交點,且交點位于置信區(qū)間內(nèi),該交點對應(yīng)的時刻即序列突變開始的時刻。

2.3 植被覆蓋度變化趨勢

基于像元的植被覆蓋度的變化趨勢,可以通過一元線性回歸分析得出,計算公式為[31]：

式中,θslope為變化趨勢的斜率,n為監(jiān)測的年數(shù),Ci為第i年的年最大植被覆蓋度。植被覆蓋年際間變化的顯著性可以通過年時間序列和植被覆蓋度的相關(guān)關(guān)系獲得。正值代表植被覆蓋度上升,負(fù)值代表植被覆蓋度下降。變化趨勢的顯著性檢驗采用F檢驗。統(tǒng)計量計算公式為：

2.4 氣候因子與植被覆蓋度的相關(guān)關(guān)系

地理系統(tǒng)是一種多要素的復(fù)雜巨系統(tǒng),其中一個要素的變化必然影響到其他要素的變化。在多要素系統(tǒng)中,研究一個要素對另外一個要素的相關(guān)性時,暫不考慮其他要素的影響,則稱為偏相關(guān)[32]。偏相關(guān)系數(shù)的計算公式為:

式中,rxy*z為變量z固定后變量x與變量y的偏相關(guān)系數(shù)； rxy為變量x與變量y的相關(guān)系數(shù)；rxz為變量y與變量z的相關(guān)系數(shù)；ryz為變量x與變量z的相關(guān)系數(shù)。

其中,x、y和z之間的相關(guān)系數(shù)計算公式為：

式中,rxy為x與y之間的相關(guān)系數(shù),xij、yij分別為第i年第j月的研究變量值；xij、yij分別為變量在各時間尺度的均值。

偏相關(guān)系數(shù)的顯著性檢驗,一般采用t檢驗法。其統(tǒng)計量計算公式為：

式中,r12,34,...,m為偏相關(guān)系數(shù)；n為樣本數(shù)；m為自變量個數(shù)。顯著水平的臨界值ta可以通過查詢t分布表獲得,若t>ta,則表明偏相關(guān)顯著；若t

2.5 殘差分析

氣候變化對植被覆蓋度變化產(chǎn)生重要影響,降水因子和氣溫因子是氣候變化影響中的重要因子。通過降水因子和氣溫因子擬合出NDVI的預(yù)測值,即為氣候變化的影響量。計算其與NDVI真實值之間的差,可作為人為因素的影響,該方法即為殘差分析,已經(jīng)得到廣泛應(yīng)用[33-35]。

ε=NDVIreal-NDVIpre

式中,NDVIreal為NDVI真實值,NDVIpre為NDVI預(yù)測值,ε>0,說明人類活動產(chǎn)生正面影響;ε<0,說明人類活動產(chǎn)生負(fù)面影響;ε=0,說明人類活動影響比較微弱。

3 結(jié)果與討論

3.1 華北平原氣候變化特征與突變分析

1981—2013年華北平原氣溫和降水量的年際和四季變化趨勢如圖3所示,表1為統(tǒng)計的氣溫和降水量的年際和四季的氣候傾向率。從圖3中可見,華北平原年均溫呈現(xiàn)整體上升的趨勢,其氣候傾向率約為0.26℃/10a(P<0.01)。氣溫的季節(jié)性變化整體呈現(xiàn)升高趨勢,只是幅度略有不同,春季(0.34℃/10a)增溫最快,冬季(0.09℃/10a)最慢。春季均溫的氣候傾向率大于年均溫,說明春季增溫對年均溫的增長貢獻(xiàn)大于其他季節(jié)。該研究結(jié)果與趙少華的研究結(jié)果一致[36]。降水量呈現(xiàn)減少的趨勢,氣候傾向率約為-18.15mm/10a(P>0.05)。與氣溫的季節(jié)性整體升高趨勢不同,降水量季節(jié)性差異更加明顯。降水量在秋季減少最快(-10.70 mm/10a),夏季減少最慢(-0.35 mm/10a),冬季呈現(xiàn)略微增長的趨勢,氣候傾向率為1.4mm/10a。雖然冬季的降水呈增加趨勢,但增加幅度小,遠(yuǎn)不及其他季節(jié)降水量的減少幅度。從氣溫與降水量的變化趨勢可見,華北平原的氣候變化有向干熱化發(fā)展的趨勢。

表1 1981—2013年氣候因子氣候傾向率統(tǒng)計結(jié)果

圖3 華北平原氣溫與降水量數(shù)據(jù)的變化趨勢Fig.3 The change trend of temperature and precipitation

圖4為華北平原年降水量和年均溫的氣候傾向率的空間分布圖。1981—2013年間,華北平原82個氣象站點的年均溫氣候傾向率大部分呈現(xiàn)增溫變化趨勢,增溫趨勢為0.32℃/10a。增溫趨勢較高區(qū)域位于河北省西部太行山地區(qū)以及山東省的威海(0.64℃/10a)和青島(0.67℃/10a)。河北承德、山東東營站點年均溫的氣候傾向率呈現(xiàn)降低趨勢,氣候傾向率分別為-0.06℃/10a和-0.28℃/10a。河北省的氣候傾向率均值高于山東省,山東省的值高于河南省。說明在全球增暖背景下,華北平原氣溫的空間變化具有一定的緯度地帶性,即：氣溫的氣候傾向率隨緯度升高而升高。春季、夏季和秋季的氣候傾向率分布特征從南到北呈現(xiàn)高—低—高的趨勢,低值區(qū)域分布也相似,主要呈帶狀沿河南、山東和河北三省交界地區(qū)向東直到渤海一帶分布。冬季從南到北呈現(xiàn)逐漸增大的趨勢。華北平原年降水量的氣候傾向率為-1.75℃/10a,整體呈減少趨勢,減少幅度最快的站點為東營。年降水量增加的區(qū)域主要位于南陽—鄭州—菏澤一線的南側(cè)。春、秋兩季呈現(xiàn)南部減少,北部增加的趨勢,夏、冬兩季呈現(xiàn)南部增加北部減少的趨勢。季節(jié)氣候傾向率增減特征變化敏感區(qū)域主要位于35°—39°N之間。

圖4 華北平原年降水量氣候傾向率(a)與年均溫氣候傾向率(b)空間分布Fig.4 The spatial distribution of climate tendency rate of precipitation (a) and temperature (b)

分別提取各年代際的氣候因子的面積要素,探尋氣溫和降水量因子各量級在不同年代際的空間分布變化。氣溫因子按照<0、0—5、5—10、10—15℃、>15℃五類統(tǒng)計,降水量因子按照<400、400—600、600—800、800—1000、1000—1500、>1500mm六類統(tǒng)計。

從圖5中可以發(fā)現(xiàn),1960—2010年,華北平原大部分地區(qū)都被10—15℃范圍的溫度控制(多年均值面積占比78.52%),其次是5—10℃的控制范圍(15.49%),>15℃的控制范圍排在第三位(5.08%),<0℃和0—5℃的平均面積占比均未超過1%。各量級區(qū)間面積占比在不同年代際的變化各不相同。1980s是各溫度范圍控制面積變化的特征年代,>15℃的面積范圍在1990s超過了5—10℃的面積范圍?？傮w上,>15℃溫度的控制范圍在變大,其它等溫線控制的范圍,從1960s到2010s面積都在減小。結(jié)合研究區(qū)年均溫的整體升高趨勢,可以得出結(jié)論：一方面,全球變暖背景下的華北平原各站點的年均溫大部分呈升高趨勢；另一方面,>15℃等溫線控制范圍的擴張,也是年均溫升高的因素之一。從圖5可見,華北平原降水量400—600mm等值線(面積占比33.97%)和600—800mm(面積占比35.06%)等值線控制范圍的平均面積相當(dāng),其次是降水量1000—1500mm等值線控制范圍(面積占比20.15%)。<400mm等值線控制面積和1000—1500mm的控制面積的平均值相當(dāng),分別為4.51%和6.87%。研究區(qū)降水量>1500mm的面積占比多年來一直低于1%,但也呈現(xiàn)略微下降趨勢?？傮w上,全區(qū)大部分地區(qū)的降水量范圍為400—600mm和600—800mm,兩種類型斑塊面積平均占比之和達(dá)到79%。400—600mm等值線控制的面積整體呈現(xiàn)增加趨勢,其余等值線的控制范圍都呈現(xiàn)下降或者維持不變的態(tài)勢。1970s—1980s之間,400—600mm降水等值線控制范圍超過了600—800mm控制范圍,成為最大面積斑塊。1980s是各降水量等值線控制面積變化的特征年代。

圖5 華北平原氣溫與降水量各統(tǒng)計量變化趨勢圖Fig.5 The change trend chart of temperature and the precipitation in Huabei Plain

華北平原的氣候環(huán)境在過去的30多年內(nèi)經(jīng)歷了突變,尤其是氣溫的突變比較明顯(圖6)。M-K檢測結(jié)果顯示：1981—2013年間,UF和UB曲線在1994年相交,且交點在置信區(qū)間內(nèi)。累積距平統(tǒng)計結(jié)果呈“V”字形,轉(zhuǎn)折點出現(xiàn)在1993年,表明華北平原氣溫經(jīng)歷了從1981—1993年之間的降溫過程和1994—2009年之間的增溫過程,2009—2013年又有小幅回落。結(jié)合兩種檢測方法結(jié)果分析得出：1981—2013年,華北平原氣溫由冷到暖的轉(zhuǎn)變時間出現(xiàn)在1993—1994年間。華北平原年降水量的M-K檢測結(jié)果顯示(圖7),UF和UB曲線在置信水平內(nèi)有5個交點。依據(jù)累積距平的統(tǒng)計結(jié)果分析,1982—1991年間,降水量累積距平呈現(xiàn)增多的變化過程；1991—2013年間,累積距平變?yōu)闇p少趨勢,降水量累積距平的最大值出現(xiàn)在1991年。結(jié)合兩種檢測結(jié)果分析得出：在1981—2013年間,華北平原降水由多到少的轉(zhuǎn)變時間發(fā)生在1991年前后。結(jié)合年均溫和降水量的突變研究結(jié)果得出：1991—1994年間是氣候因子活躍變化時期,該時期為突變特征時期。

圖6 年均溫M-K檢測與累積距平圖Fig.6 M-K test and Cumulative anomaly of temperature

圖7 年均降水量M-K檢測與累積距平圖Fig.7 M-K test and Cumulative anomaly of precipitation

3.2 植被覆蓋度空間分布與變化趨勢特征

圖8是華北平原1981—2013年植被覆蓋度空間分布格局。研究區(qū)年平均植被覆蓋度為0.6471。城市群、較大的水面和海岸線地區(qū)的植被覆蓋度相對稀疏,甚至沒有植被覆蓋。植被覆蓋度較低和中等水平植被覆蓋度區(qū)域,主要分布在河北省西北部半干旱的農(nóng)牧交織帶以及河南省和山東省的低山丘陵地區(qū)。植被覆蓋度較高的區(qū)域主要集中在華北平原中部,包括河北平原、河南省的中東部、山東省西部以及安徽北部地區(qū)。植被覆蓋度最高區(qū)域集中分布在燕山和秦嶺東延的伏牛山山脈中的幾乎不受人類影響的區(qū)域。草原生態(tài)區(qū)多年平均植被覆蓋度約為0.5504。森林生態(tài)區(qū)平均植被覆蓋度約為0.6526。農(nóng)田生態(tài)區(qū)平均植被覆蓋度約為0.6520,水平與森林生態(tài)區(qū)差別不大。

從月尺度植被覆蓋度變化趨勢圖(圖9)中可以看出,在植被生長季,草原生態(tài)區(qū)和森林生態(tài)區(qū)的植被覆蓋度在植物生長期內(nèi)表現(xiàn)出先增加后下降的趨勢。兩者出現(xiàn)峰值均在8月份,分別為0.541和0.549。草原生態(tài)區(qū)的植被覆蓋度整體低于森林生態(tài)區(qū)。農(nóng)業(yè)生態(tài)區(qū)的植被覆蓋度在植物生長季的變化相對復(fù)雜,6月份植被覆蓋度出現(xiàn)明顯的降低,隨后持續(xù)增高,在8月出現(xiàn)峰值0.59。華北平原農(nóng)作物以冬小麥為主, 6月是小麥?zhǔn)斋@的季,大面積裸露的農(nóng)田造成了該生態(tài)區(qū)植被覆蓋度的突然降低。4—5月的峰值期,農(nóng)業(yè)生態(tài)區(qū)的植被類型主要是以小麥為主。8月的峰值期主要以玉米、高粱以及經(jīng)濟作物為主,且生態(tài)區(qū)內(nèi)栽培果樹等都達(dá)到生長季峰值。植被類型的不同,是造成8月份峰值高于4—5月峰值的主要原因。

圖8 植被覆蓋度空間分布Fig.8 Spatial patterns of vegetation coverage

圖9 植被生長期內(nèi)各生態(tài)區(qū)覆蓋度變化趨勢Fig.9 Variation trend of coverage of ecotope

圖10為基于像元的植被覆蓋度變化趨勢空間分布圖。華北平原植被覆蓋度呈現(xiàn)增加趨勢,增加速率為0.00039/10a,呈現(xiàn)增加趨勢的區(qū)域占總面積的55%,主要位于河北省西北部山區(qū)、山東省西南部、河南省東部以及安徽省北部地區(qū)。呈現(xiàn)減少趨勢的區(qū)域主要集中于人類活動相對頻繁的地區(qū),包括河南省中部和西部,山東省北部,以及人口密集京津唐地區(qū)。在各生態(tài)區(qū)的趨勢變化中,森林生態(tài)區(qū)的增加趨勢最快,約為0.000034/10a；農(nóng)業(yè)生態(tài)區(qū)的增加速率最慢,約為0.00024/10a,草原生態(tài)區(qū)增加的速率約為0.0029/10a。

3.3 氣候因子與植被覆蓋度的關(guān)系

3.3.1 不同年代際的水熱條件對植被覆蓋度的影響

植被生長受到氣象要素、土地利用和CO2的等因素的多重影響,其中氣溫和降水是植被生長的關(guān)鍵性因素[37],水熱條件是決定植被空間分布及其變化的主要非生物因素[23]。同時,地表植被覆蓋度的變化也影響局部氣候[38]。為探討不同水熱組合與植被覆蓋度之間的相互作用關(guān)系,選擇1980s、1990s和2010s 3個時期的相關(guān)數(shù)據(jù),分別計算溫度累積距平、降水量距平百分率和植被覆蓋度的變化趨勢,探尋不同水熱組合影響下的植被覆蓋度變化情況。1980s、1990s和2010s 3個時期的氣溫和降水相較于多年均值的盈虧情況如表2所示：1980s年處于氣溫偏低、降雨偏多的水熱狀況,1990s處于氣溫偏高,降雨偏多的水熱狀況,2010s處于氣溫偏高降雨偏少的水熱狀況。

圖10 植被覆蓋度變化趨勢空間分布Fig.10 Dynamic of vegetation coverage

圖11 1990s研究區(qū)氣候因子偏差分析Fig.11 Deviation analysis of climatic factors in 1990s

時間尺度Timescale溫度距平百分比Temperaturdeparturepercentage降水量距平百分比Precipitationdeparturepercentage1980s1990s2010s1980s1990s2010s春季Spring-0.5034-0.1510.791-0.021460.101-0.119夏季Summer-0.4150.0720.371-0.0150.014-0.003秋季Autumn-0.3480.0830.2790.140-0.105-0.010冬季W(wǎng)inter-0.6760.4070.191-0.179-0.0100.219年Year-0.5120.0960.4490.00020.008-0.011

在1980s,華北平原植被覆蓋度呈現(xiàn)整體增加的趨勢,速率約為0.03/10a。農(nóng)田生態(tài)區(qū)(0.08/10a)和森林生態(tài)區(qū)(0.04/10a)呈增加趨勢,草原生態(tài)區(qū)呈現(xiàn)減少趨勢(-0.03/10a)。該年代際的溫低水多的水熱條件對整體以及農(nóng)田和森林生態(tài)區(qū)的植被覆蓋度呈現(xiàn)增加的影響,對草原生態(tài)區(qū)則呈現(xiàn)減少的影響。

在1990s,華北平原植被覆蓋度的呈現(xiàn)整體減少趨勢,速率約為-0.02/10a。農(nóng)業(yè)生態(tài)區(qū)(-0.05/10a)和森林生態(tài)區(qū)(-0.03/10a)呈減少趨勢,農(nóng)業(yè)生態(tài)區(qū)的減少速率最快,草原生態(tài)區(qū)呈現(xiàn)增加的趨勢(0.02/10a)。結(jié)合該時期的年均溫和年降水量的距平百分比可以發(fā)現(xiàn),雨熱不同步抑制了植被的生長(圖11)。草原生態(tài)區(qū)的植被覆蓋度呈現(xiàn)增加趨勢,主要得益于5℃等溫線和400mm降水等值線的北移。以400mm降水等值線為例,多年均位于草原生態(tài)區(qū)的邊緣,而該時段卻幾乎覆蓋了整個草原生態(tài)區(qū)。相較于整個研究區(qū)的雨熱不均, 平原草原區(qū)降水的增多結(jié)合適宜的溫度,導(dǎo)致植被覆蓋度呈現(xiàn)增加趨勢。

在2010s,華北平原植被覆蓋度呈現(xiàn)整體增加趨勢,速率約為0.02/10a。農(nóng)田生態(tài)區(qū)(0.02/10a)、森林生態(tài)區(qū)(0.07/10a)和草原生態(tài)區(qū)(0.13/10a)均呈現(xiàn)增加趨勢。研究區(qū)的干熱化趨勢已經(jīng)呈現(xiàn),區(qū)域內(nèi)植被覆蓋度呈現(xiàn)整體增加的趨勢。雖然降水量整體有所降低,但是與多年均值的空間分布相比,沒有呈現(xiàn)明顯差異。說明區(qū)域內(nèi)各生態(tài)區(qū)的氣溫和降水量達(dá)到了較好的水熱組合條件,

綜合1981—2010的華北平原不同年代際的植被覆蓋度對不同水熱組合的分析結(jié)果得出：在基于年代際的時間尺度研究中,氣溫的增加,對華北平原農(nóng)業(yè)生態(tài)區(qū)的影響體現(xiàn)為植被覆蓋度的穩(wěn)定增加,但是速率在不斷降低。降水量的減少對農(nóng)業(yè)生態(tài)區(qū)植被覆蓋度的影響相對較小。草原生態(tài)區(qū)植被覆蓋度的變化趨勢在一定程度上與氣溫的變化趨勢一致。當(dāng)氣溫偏高降雨偏多時,森林生態(tài)區(qū)的植被覆蓋度呈減少趨勢,其余的水熱組合方式植被覆蓋度均呈現(xiàn)增加趨勢。

3.3.2 氣候突變前后植被覆蓋度變化特征

在3.1節(jié)得出結(jié)論1991—1994年為研究時段華北平原氣候因子突變的特征時段。因要對氣候突變前后植被覆蓋度對氣候因子響應(yīng)情況進(jìn)行比對,因此暫將華北平原的氣候突變年選在1994年,此時降水量因子和氣溫因子均已經(jīng)歷了突變變化。

將1981—1994年和1995—2013年兩個時間段的遙感影像分別計算,得出基于像元的植被覆蓋度分布和變化速率統(tǒng)計結(jié)果圖(圖12)。華北平原年均植被覆蓋度在氣候突變前后分別為0.662和0.663。結(jié)合各生態(tài)區(qū)年均植被覆蓋度的統(tǒng)計結(jié)果圖發(fā)現(xiàn),草原生態(tài)系統(tǒng)區(qū)在氣候突變后由0.544增加到0.590,顯著增幅約為8%。農(nóng)田生態(tài)系統(tǒng)區(qū)由0.668增加到0.669。森林生態(tài)系統(tǒng)區(qū)由0.668減少到0.663。說明植被覆蓋度在研究區(qū)內(nèi)部,空間異質(zhì)性較強,氣候突變對草原生態(tài)區(qū)的植被覆蓋度的影響比其他生態(tài)區(qū)顯著。氣候突變前后植被覆蓋度的平均增加趨勢為0.018/10a和0.006/10a,增加幅度降低了0.012/10a。在各生態(tài)區(qū)中,降低幅度最快的是農(nóng)田生態(tài)區(qū)(降低了0.016/10a),其次是森林生態(tài)區(qū)(降低了0.005/10a),草原生態(tài)區(qū)變化最小(降低了0.001/10a)。在氣候突變前,增加速率最快的是農(nóng)田生態(tài)區(qū),氣候突變后草原生態(tài)區(qū)植被覆蓋度增加速率最快。森林生態(tài)區(qū)的速率相對較慢。氣候突變對農(nóng)田生態(tài)區(qū)的植被覆蓋度變化速率影響最大,其次是森林生態(tài)區(qū),草原生態(tài)區(qū)幾乎不受氣候突變的影響?？梢?從植被覆蓋度高低的角度出發(fā),氣候突變對草原生態(tài)區(qū)影響最大；從變化速率角度出發(fā),對農(nóng)田生態(tài)區(qū)的影響最大。

圖12 氣候突變前后各生態(tài)區(qū)植被覆蓋度變化情況統(tǒng)計結(jié)果Fig.12 The results of climate change effects on vegetation coverage of ecological zones

3.3.3 植被覆蓋度與氣候因子的相關(guān)性分析

圖13為年際植被覆蓋度與年均溫和年降水量的偏相關(guān)性空間分布結(jié)果。如圖顯示：在1981—2013年之間,華北平原年際植被覆蓋度與年降水量的整體平均偏相關(guān)系數(shù)為0.018,年植被覆蓋度與年溫度的整體平均偏相關(guān)系數(shù)為0.006,其中呈正相關(guān)的像元分別占總像元面積的59%和51%。總體來說,植被覆蓋度與降水量的偏相關(guān)性,高于其與年均溫的偏相關(guān)性。兩者在空間分布上都表現(xiàn)出明顯的空間異質(zhì)性。

農(nóng)業(yè)生態(tài)區(qū)的植被覆蓋度與年均溫和年降水量的平均相關(guān)系數(shù)為0.01和-0.001,與之對應(yīng)的呈正相關(guān)的區(qū)域面積占分區(qū)面積的52%和48%,表明在農(nóng)業(yè)生態(tài)區(qū)氣溫是影響區(qū)域內(nèi)植被長勢的主導(dǎo)因素。從圖13b可以看出,農(nóng)業(yè)生態(tài)區(qū)與氣溫呈正相關(guān)的區(qū)域主要分布在河北平原、河南東南部、安徽北部,呈現(xiàn)負(fù)相關(guān)的區(qū)域主要集中于山東半島以及山東省西北部地區(qū),且相對集中。華北農(nóng)業(yè)主要以種植冬小麥、玉米、棉花為主,冬小麥生長周期長,且在主要生長期多抽取地下水灌溉,導(dǎo)致區(qū)域植被覆蓋度與降水量相關(guān)性較弱,甚至呈負(fù)相關(guān)關(guān)系。森林生態(tài)區(qū)的植被覆蓋度與年均溫和年降水量的平均相關(guān)系數(shù)為0.03和-0.001,與之對應(yīng)的呈正相關(guān)區(qū)域面積占比為53%和49%,表明氣溫是影響森林生態(tài)區(qū)植被覆蓋度的主導(dǎo)因子。森林生態(tài)區(qū)植被覆蓋度與降水呈現(xiàn)負(fù)相關(guān),原因在于降雨意味著輻射減少和溫度的降低,從而影響植被生長[13,27]。森林生態(tài)區(qū)植被覆蓋度與氣溫因子的相關(guān)性有區(qū)域差異。北部太行山、燕山山脈的植被生長更多受到低溫的影響,因此與氣溫呈正相關(guān)的區(qū)域面積更大,約為55%。山東半島丘陵地區(qū)的森林植被覆蓋度與氣溫呈現(xiàn)正相關(guān)的區(qū)域占比為35%。山東半島丘陵地區(qū)相較于太行山、燕山山脈更靠南,植被生長所需的熱量資源更豐富。熱量資源較為充足的地區(qū),增溫對植被的影響會變?nèi)鮗39]。太行山脈地區(qū)的降水量與植被覆蓋度相關(guān)性高于山東半島丘陵地區(qū)的相關(guān)性,燕山地區(qū)的相關(guān)性最低。除了各地區(qū)的氣候、地形、土壤類型不同以外,太行山山脈的植被類型以落葉灌和草甸為主,山東半島丘陵地區(qū)植被以溫帶針葉林和闊葉林為主,燕山山脈的植被類型以灌叢和溫帶闊葉林為主。不同的植被類型,對降水的響應(yīng)機制不同,增強了相關(guān)系數(shù)的空間異質(zhì)性。草原生態(tài)區(qū)植被覆蓋度與年均溫和年降水量的平均相關(guān)系數(shù)為-0.015和0.18,與之對應(yīng)的呈正相關(guān)的區(qū)域面積分別占分區(qū)面積的39%和92%.說明降水量是影響該區(qū)域植被生長的主要氣候因子。該結(jié)論與張戈麗和穆少杰[28,31]的結(jié)論一致。該區(qū)域內(nèi)的草原生態(tài)區(qū)位于內(nèi)蒙古高原,地處我國水資源相對緊缺地區(qū),因此降水對該地區(qū)的影響要高于氣溫的影響。

圖13 植被覆蓋度與降水量(a)和溫度(b)的相關(guān)系數(shù)Fig.13 Coefficient of vegetation and temperature (a) and precipitation (b)

3.4 氣候?qū)χ脖桓采w度的滯后性研究

以1982—2006年植被覆蓋度和氣象數(shù)據(jù)為基礎(chǔ),分別計算生長季3月—11月的月植被覆蓋度和月降水量、月均溫的偏相關(guān)系數(shù),以探討月植被覆蓋度與前一月以及當(dāng)月氣候因子的相關(guān)性,進(jìn)而得出植被覆蓋度對氣候因子的滯后性變化。

華北平原當(dāng)月植被覆蓋度與當(dāng)月的月均溫和月降水量的平均相關(guān)系數(shù)為0.14和0.13(圖14),與之對應(yīng)的呈正相關(guān)的區(qū)域面積分別占總面積的88%和94%,其中通過P<0.05檢驗的面積占比分別為27.45%和31.15%。華北平原月植被覆蓋度與前一月的月均溫和月降水量的平均相關(guān)系數(shù)為0.008和0.227(圖15),與之對應(yīng)的呈正相關(guān)的區(qū)域面積分別占總面積的47%和85%,通過P<0.05檢驗的面積占比分別為8.01%和54.65%。總體而言,植被覆蓋度與前一個月的月均溫的相關(guān)性低于當(dāng)月均溫的相關(guān)性,與前一個月的月降水量的相關(guān)性高于當(dāng)月降水量的相關(guān)性。說明在月尺度上,降水量對整體植被覆蓋度有時滯效應(yīng)。植被覆蓋度與前一月氣象因子關(guān)系的空間分布具有顯著緯度地帶性：月均溫與植被覆蓋度的相關(guān)系數(shù)呈現(xiàn)隨緯度升高而降低的總體趨勢,月降水量與植被覆蓋度的相關(guān)系數(shù)呈現(xiàn)隨緯度增高而升高的趨勢?？傮w月尺度的相關(guān)性要遠(yuǎn)高于年際水平。

農(nóng)業(yè)生態(tài)區(qū)植被覆蓋度與前一月降水量的平均相關(guān)系數(shù)為0.3,高于其與當(dāng)月降水量的相關(guān)系數(shù)0.131。其與前一月的月均溫的平均相關(guān)性為-0.08,低于其與當(dāng)月均溫的相關(guān)性0.109,這表明,在月水平上農(nóng)業(yè)生態(tài)區(qū)植被生長對降水量的響應(yīng)存在時滯效應(yīng),與月均溫則沒有明顯的時滯效應(yīng)。氣溫隨季節(jié)變化是一種自然規(guī)律,在植物生長季,氣溫是植被某些生長期的關(guān)鍵誘因,植被遵循氣溫的變化進(jìn)入不同的生長期,因而在月尺度上,氣溫對植物生長不存在時滯效應(yīng)。森林生態(tài)區(qū)植被覆蓋度與前一月的月均溫的平均相關(guān)系數(shù)為0.13,低于其與當(dāng)月的數(shù)值0.205,其與前一月降水量的相關(guān)系數(shù)為0.12,與當(dāng)月降水量的數(shù)值0.11基本相同,整體而言不存在時滯效應(yīng)。這與穆少杰[31]的研究結(jié)果一致。其中當(dāng)月的月均溫與月植被覆蓋度的相關(guān)系數(shù)呈正相關(guān)的面積占分區(qū)面積的93%,范圍分布在0—0.35之間。草原生態(tài)區(qū)植被覆蓋度與前一月的月均溫的平均相關(guān)系數(shù)為0.24,高于其與當(dāng)月平均相關(guān)性0.046,其與前一月降水量的相關(guān)系數(shù)為0.07,低于其與當(dāng)月平均相關(guān)性系數(shù)0.26。這表明,草原生態(tài)分區(qū)對氣溫的響應(yīng)具有時滯效應(yīng)。在月水平上,相關(guān)性結(jié)果表現(xiàn)出一定的經(jīng)度地帶性：從西到東,植被生長與降雨的相關(guān)性遞減,同時與溫度的相關(guān)性則表現(xiàn)出遞增趨勢。

圖14 植被覆蓋度與當(dāng)月溫度(a)和降水量(b)的相關(guān)系數(shù)Fig.14 The coefficient of vegetation and temperature (a) and precipitation (b) in month scale

圖15 植被覆蓋度與前一月溫度(a)和降水量(b)的相關(guān)系數(shù)Fig.15 The coefficient of vegetation and temperature (a) and precipitation (b) of the month bfore

3.5 人類活動對植被覆蓋度影響力的趨勢分析

圖16 人類活動影響力變化趨勢圖Fig.16 Dynamic of vegetation coverage

為探討人類活動對植被覆蓋度的影響,本文用殘差分析計算人類活動對植被覆蓋度的影響,并在此基礎(chǔ)上計算人類活動影響的變化趨勢。圖16為影響植被覆蓋度變化的人類活動影響力的變化趨勢空間分布圖。統(tǒng)計結(jié)果顯示,人類活動影響力的整體平均變化趨勢為0.004/10a,說明人類活動影響力呈現(xiàn)繼續(xù)增強的趨勢。圖中正值表示人類活動影響力繼續(xù)增強,負(fù)值表示影響力減弱。從圖中可見,河北北部和河南中東部地區(qū)人類活動影響力在持續(xù)下降,其中京津唐地區(qū)減弱速度最快。以北京市為例,不斷涌入的人口加快了城市化進(jìn)程,中心城區(qū)的不透水面建設(shè)接近飽和,導(dǎo)致人類對植被覆蓋度變化的影響力的變?nèi)?。山東省西南部、安徽和江蘇北部地區(qū)人類活動影響力呈現(xiàn)持續(xù)增強的態(tài)勢。從生態(tài)學(xué)角度出發(fā),草原生態(tài)區(qū)的人類影響力增加速率最快(0.005/10a),農(nóng)業(yè)生態(tài)區(qū)的影響力增加最慢(0.003/10a)。農(nóng)業(yè)生態(tài)區(qū)中,果樹和糧食混種區(qū)的人類影響增加的趨勢最快(0.005/10a),糧食和經(jīng)濟農(nóng)作物種植區(qū)次之(0.004/10a),果樹種植區(qū)最慢(0.003/10a)。同時,結(jié)合植被覆蓋度變化趨勢圖,統(tǒng)計研究區(qū)內(nèi)人類活動的整體作用力為-5%(其中,植被覆蓋度呈增加趨勢的區(qū)域的相對作用力為正值,植被覆蓋度變化趨勢呈退化趨勢區(qū)域的相對作用力為負(fù)值,兩者之差即為整體相對作用力),說明人類活動將繼續(xù)以負(fù)面影響的作用力為主。氣候突變前植被改善區(qū)的人類活動作用力低于植被退化區(qū)的相對作用力；氣候突變后植被改善區(qū)的人類活動作用力高于退化區(qū)。農(nóng)業(yè)生態(tài)區(qū)的作用力一直呈現(xiàn)負(fù)面作用,改善區(qū)的相對作用力減弱了4%。草原生態(tài)區(qū)的作用力呈現(xiàn)改善力增強的趨勢,幅度達(dá)到12%。

4 結(jié)論與討論

4.1 結(jié)論

本文基于GIMMS-NDVI和MODIS-NDVI遙感數(shù)據(jù),結(jié)合1981—2013年氣象數(shù)據(jù),利用殘差分析、趨勢分析、相關(guān)性分析,揭示了華北平原近30年氣候和植被覆蓋度的空間變化特征,并對氣候變化對各生態(tài)區(qū)內(nèi)植被覆蓋度的影響進(jìn)行了分析,結(jié)論如下：

(1)1981—2013年間,華北平原氣溫整體呈現(xiàn)顯著上升趨勢(0.20℃/10a,P<0.01),一方面全球氣候變暖背景下大部分站點氣溫整體升高,另一方面春季氣溫的升高、>15℃等溫線控制范圍的擴大和年均溫0℃等值線在華北平原的消失,也是區(qū)域平均氣溫升高的誘因。其中,15℃等值線的年代際間南北波動范圍達(dá)到2個緯度, 有大于200km的北移。華北平原降水整體呈現(xiàn)顯著減少趨勢(-1.75mm/10a,P>0.05),其中秋季降水量減少過快,加上400—600mm降水等值線控制范圍的擴大,壓縮了600—800mm和800—1000mm降水等值線的范圍,這些因素都造成了區(qū)域降水量的減少。四季氣候傾向率的特征變化敏感區(qū)域主要位于35°—39°N之間。1991—1994年為華北平原氣候的突變時期。

(2)華北平原植被覆蓋度空間分布具有顯著的地區(qū)差異性,整體上呈現(xiàn)上升趨勢,呈現(xiàn)增加趨勢的面積占比為55%。草原生態(tài)區(qū)的年均植被覆蓋度較低,森林生態(tài)區(qū)和農(nóng)業(yè)生態(tài)區(qū)相對較高。從植被覆蓋度和植被覆蓋度變化趨勢的空間分布來看,人類活動不僅加速了區(qū)域植被覆蓋度的降低,也加劇了降低速率的變快。在基于年代際的時間尺度研究中,氣溫的增加,對華北平原農(nóng)業(yè)生態(tài)區(qū)的影響體現(xiàn)為植被覆蓋度的穩(wěn)定增加,但是速率在不斷降低；草原生態(tài)區(qū)植被覆蓋度的變化趨勢在一定程度上與氣溫的變化趨勢一直。當(dāng)氣溫偏高降雨偏多時；森林生態(tài)區(qū)的植被覆蓋度呈減少趨勢,其余的水熱組合方式植被覆蓋度均呈現(xiàn)增加趨勢。

(3)總體上,研究區(qū)月尺度植被覆蓋度與氣候因子的相關(guān)性高于年尺度上的值。植被覆蓋度與年降水量的偏相關(guān)性高于其與年均溫的偏相關(guān)性,但是存在明顯的區(qū)域差異。年均溫對農(nóng)業(yè)生態(tài)區(qū)和森林生態(tài)區(qū)的植被覆蓋度的影響更大,草原生態(tài)區(qū)的植被覆蓋度對年降水量的依賴性更大。在月尺度上年降水量對植被的影響具有時滯效應(yīng)。氣溫對草原生態(tài)區(qū)具有時滯效應(yīng),降雨對農(nóng)業(yè)生態(tài)區(qū)具有時滯效應(yīng)。

(4)華北平原干熱化的氣候突變降低了植被覆蓋度的增加趨勢。在氣候突變前,增加速率最快的是農(nóng)田生態(tài)區(qū)；氣候突變后草原生態(tài)區(qū)植被覆蓋度增加速率最快。從植被覆蓋度高低的角度出發(fā),氣候突變對草原生態(tài)區(qū)影響最大；從變化速率角度出發(fā),對農(nóng)田生態(tài)區(qū)的影響最大。雖然氣候突變后,人類活動在植被覆蓋度的改善區(qū)的相對作用力高于退化區(qū),但是就整體而言,人類活動的影響力還在持續(xù)增強,且呈現(xiàn)出在退化區(qū)的作用力高于改善區(qū)的趨勢。

4.2 討論

本文從生態(tài)學(xué)角度出發(fā),對不同生態(tài)區(qū)植被覆蓋度的分布與變化趨勢進(jìn)行了分析討論,總結(jié)了基于生態(tài)學(xué)分區(qū)的植被覆蓋的分布和變化規(guī)律,探討了植被覆蓋度與氣候因子的關(guān)系。相較于傳統(tǒng)統(tǒng)計分析結(jié)果,相關(guān)系數(shù)有一些不同。與把整個研究區(qū)域作為統(tǒng)計單元的研究方法相比,基于像元的小樣本統(tǒng)計能更好的體現(xiàn)NDVI和氣候因子的空間差異性。在統(tǒng)計年限有限的情況下,顯著的空間差異可能會影響小樣本的統(tǒng)計分析結(jié)果,掩飾了氣候因子與NDVI的相關(guān)性。另外,植被生長季的有效降水跟NDVI的相關(guān)性可能更高,此次研究沒有排除非生長季的無關(guān)水影響。

水熱組合是決定植被生長的重要非生物因素。雖然本研究利用不同年代際的水熱組合對植被覆蓋度的響應(yīng)進(jìn)行了粗略探討,但是水熱組合的形式多樣,且不同植物種類對水熱組合的脅迫響應(yīng)機制也各不相同。適宜的水熱組合加速植物的生長,不適宜的水熱組合將對植被生長產(chǎn)生抑制作用。因此,在全球氣候不斷變暖,厄爾尼諾現(xiàn)象頻發(fā)的大背景下,研究植被對不同水熱組合的響應(yīng)機理,探討不同植被類型氣候變化的響應(yīng)和分布,可以做為下一步研究內(nèi)容。

除了氣候因子,人類活動也是影響植被覆蓋度變化的重要驅(qū)動因素,雖然本文對人類活動的整體影響做了簡單分析,但是人類活動對植被覆蓋度的影響機理錯綜復(fù)雜。特別是進(jìn)入21世紀(jì),人類活動強度和范圍逐漸擴大,包括人口增長、城鎮(zhèn)化發(fā)展、能源消耗結(jié)構(gòu)調(diào)整、退耕還林、退耕還牧、農(nóng)業(yè)結(jié)構(gòu)的調(diào)整等,都會對區(qū)域植被覆蓋度產(chǎn)生嚴(yán)重影響。因此,如何細(xì)化人類活動對植被覆蓋度的影響,定位重要的人類活動因子,消除影響植被變化的不利因素,將對植被覆蓋度的變化產(chǎn)生積極影響,對于區(qū)域生態(tài)安全也有積極意義。

[1] IPCC. Climate change 2013: the physical science basis // Stocker T F, Qin D, Plattner G K, Tignor M, Allen S K, Boschung J, Nauels A, Xia Y, Bex V, Midgley P M, eds. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom, New York, NY, USA: Cambridge University Press, 2013: 5-5.

[2] 任國玉, 郭軍, 徐銘志, 初子瑩, 張莉, 鄒旭凱, 李慶祥, 孫小寧. 近50年中國地面氣候變化基本特征. 氣象學(xué)報, 2005, 63(6): 942-956.

[3] Zhao M, Running S W. Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science, 2010, 329(5994): 940-943.

[4] Piao S L, Mohammat A, Fang J Y, Cai Q, Feng J M. NDVI-based increase in growth of temperate grasslands and its responses to climate changes in China. Global Environmental Change, 2006, 16(4): 340-348.

[5] 劉光生, 王根緒, 胡宏昌, 任東興, 孫向陽. 長江黃河源區(qū)近45年氣候變化特征分析. 資源科學(xué), 2010, 32(8): 1486-1492.

[6] Wang Q, Zhang Q P, Zhou W. Grassland coverage changes and analysis of the driving forces in Maqu county. Physics Procedia, 2012, 33: 1292-1297.

[7] 孫紅雨, 王長耀, 牛錚, 布和敖斯?fàn)? 李兵. 中國地表植被覆蓋度變化及其與氣候因子關(guān)系──基于NOAA時間序列數(shù)據(jù)分析. 遙感學(xué)報, 1998, 2(3): 204-210.

[8] Myneni R B, Keeling C D, Tucker C J, Asrar G, Nemani R R. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature, 1997, 386(6626): 698-702.

[9] Tucker C J, Slayback D A, Pinzon J E, Los S O, Myneni R B, Taylor M G. Higher northern latitude normalized difference vegetation index and growing season trends from 1982 to 1999. International Journal of Biometeorology, 2001, 45(4): 184-190.

[10] Park H S, Sohn B J. Recent trends in changes of vegetation over East Asia coupled with temperature and rainfall variations. Journal of Geophysical Research, 2010, 115(D14), doi: 10.1029/2009JD012752.

[11] Piao S L, Wang X H, Ciais P, Zhu B, Wang T, Liu J. Changes in satellite-derived vegetation growth trend in temperate and boreal Eurasia from 1982 to 2006. Global Change Biology, 2011, 17(10): 3228-3239.

[12] Wang X H, Piao S L, Ciais P, Li J S, Friedlingstein P, Koven C, Chen A P. Spring temperature change and its implication in the change of vegetation growth in North America from 1982 to 2006. Proceedings of the National Academy of Sciences of the United States of America, 2011, 108(4):1240-1245.

[13] 張學(xué)珍, 戴君虎, 葛全勝. 1982-2006 年中國東部春季植被變化的區(qū)域差異. 地理學(xué)報, 2012, 67(1): 53-61.

[14] Jing X, Yao W Q, Wang J H, Song X Y. A study on the relationship between dynamic change of vegetation coverage and precipitation in Beijing′s mountainous areas during the last 20 years. Mathematical and Computer Modelling, 2011, 54(3/4): 1079-1085.

[15] 畢曉麗, 王輝, 葛劍平. 植被歸一化指數(shù)(NDVI)及氣候因子相關(guān)起伏型時間序列變化分析. 應(yīng)用生態(tài)學(xué)報, 2005, 16(2): 284-288.

[16] 孫睿, 劉昌明, 朱啟疆. 黃河流域植被覆蓋度動態(tài)變化與降水的關(guān)系. 地理學(xué)報, 2001, 56(6): 667-672.

[17] 賀偉, 布仁倉, 熊在平, 胡遠(yuǎn)滿. 1960-2005年東北地區(qū)氣溫和降水變化趨勢. 生態(tài)學(xué)報, 2013, 33(2): 519-531.

[18] 劉園, 王穎, 楊曉光. 華北平原參考作物蒸散量變化特征及氣候影響因素. 生態(tài)學(xué)報, 2010, 30(4): 923-932.

[19] 莫興國, 劉蘇峽, 林忠輝, 邱建秀. 華北平原蒸散和GPP格局及其對氣候波動的響應(yīng). 地理學(xué)報, 2011, 66(5): 589-598.

[20] Wu S H, Yang Q Y, Zheng D. Delineation of eco-geographic regional system of China. Journal of Geographical Sciences, 2003, 13(3): 309-315.

[21] 徐敏云, 李運起, 王堃, 曹玉鳳, 于海良, 李雪楓, 李連樹, 景福軍, 李佳祥, 謝帆. 河北省草地資源分布及植被特征動態(tài). 草業(yè)學(xué)報, 2009, 18(6): 1-11.

[22] Hou X Y, Wu T, Yu L J, Qian S. Characteristics of multi-temporal scale variation of vegetation coverage in the Circum Bohai Bay Region, 1999-2009. Acta Ecologica Sinica, 2012, 32(6): 297-304.

[23] 陳效逑, 王恒. 1982-2003年內(nèi)蒙古植被帶和植被覆蓋度的時空變化. 地理學(xué)報, 2009, 64(1): 84-94.

[24] 趙玉萍, 張憲洲, 王景升, 沈振西. 1982年至2003年藏北高原草地生態(tài)系統(tǒng)NDVI與氣候因子的相關(guān)分析. 資源科學(xué), 2009, 31(11): 1988-1998.

[25] Tucker C J, Pinzon J E, Brown M E, Slaybak D A, Pak E W, Mahoney R, Vermote E F. An extended AVHRR 8 km NDVI dataset, compatible with MODIS and SPOT vegetation NDVI data. International Journal of Remote Sensing, 2005, 26(20): 4485-4498.

[26] 張月叢, 趙志強, 李雙成, 孟憲峰. 基于SPOT NDVI的華北北部地表植被覆蓋變化趨勢. 地理研究, 2008, 27(4): 745-754.

[27] 信忠保, 許炯心, 鄭偉. 氣候變化和人類活動對黃土高原植被覆蓋變化的影響. 中國科學(xué) D輯: 地球科學(xué), 2007, 37(11): 1504-1514.

[28] 張戈麗, 徐興良, 周才平, 張宏斌, 歐陽華. 近30年來呼倫貝爾地區(qū)草地植被變化對氣候變化的響應(yīng). 地理學(xué)報, 2011, 66(1): 47-58.

[29] 符淙斌, 王強. 氣候突變的定義和檢測方法. 大氣科學(xué), 1992, 16(4): 482-493.

[30] 王力群, 李崇銀, 譚言科, 嚴(yán)冬. 冬季蒙古高壓的時空變化和突變特征研究. 高原氣象, 2010, 29(1): 23-32.

[31] 穆少杰, 李建龍, 陳奕兆, 剛成誠, 周偉, 居為民. 2001-2010年內(nèi)蒙古植被覆蓋度時空變化特征. 地理學(xué)報, 2012, 67(9): 1255-1268.

[32] 徐建華. 現(xiàn)代地理學(xué)中的數(shù)學(xué)方法(第二版). 北京: 高等教育出版社, 2002: 43-47.

[33] 黃森旺, 李曉松, 吳炳方, 裴亮. 近25年三北防護林工程區(qū)土地退化及驅(qū)動力分析. 地理學(xué)報, 2012, 67(5): 589-598.

[34] 王靜, 郭妮, 蔡迪花, 鄧振鏞. 瑪曲縣草地退牧還草工程效果評價. 生態(tài)學(xué)報, 2009, 29(3): 1276-1284.

[35] 許端陽, 康相武, 劉志麗, 莊大方, 潘劍君. 氣候變化和人類活動在鄂爾多斯地區(qū)沙漠化過程中的相對作用研究. 中國科學(xué) D輯: 地球科學(xué), 2009, 39(4): 516-528.

[36] 趙少華, 楊永輝, 邱國玉, 樊彤, 胡玉昆. 河北平原34年來氣候變化趨勢分析. 資源科學(xué), 2007, 29(4): 109-113.

[37] 李小燕. 黃土高原植被對水熱狀況的響應(yīng)研究. 地理科學(xué), 2013, 33(7): 865-872.

[38] 張翀, 任志遠(yuǎn), 袁鑫. 西北地區(qū)NDVI對水熱條件年內(nèi)變化的響應(yīng)及其空間特征. 資源科學(xué), 2011, 33(12): 2356-2361.

[39] 徐興奎, 林超暉, 薛峰. 氣象因子與地表植被生長相關(guān)性分析. 生態(tài)學(xué)報, 2003, 23(2): 221-230.

Temporal analysis of climate change and its relationship with vegetation cover on the north china plain from 1981 to 2013

A Duo1,2,3, ZHAO Wenji1,2,3,*, GONG Zhaoning1,2,3, ZHANG Min1,2,3, FAN Yunbao1,2,3

1CollegeofResourcesEnvironment&Tourism,CapitalNormalUniversity,Beijing100048,China2KeyLaboratoryof3DInformationAcquisitionandApplicationofMinistry,Beijing100048,China3KeyLaboratoryofResourcesEnvironmentandGISofBeijingMunicipal,Beijing100048,China

Global climate change has resulted in significant vegetation changes in the past fifty years. The North China Plain, the most important grain production base of China, is experiencing prominent warming and drying. Vegetation cover, which is used to monitor vegetation change, is affected by climate change, predominantly temperature and precipitation. In the present study, GIMMS (Glaobal Inventory Modelling and Mapping Studies)NDVI (Normalized Difference Vegetation Index), MODIS (Moderate-resolution Imaging Spectrometer) NDVI, and climate data, during 1981—2013, were used to investigate the spatial distribution and changes of vegetation. The relationship between climate and vegetation on the North China Plain were also analyzed on different spatial (agricultural, forest, and grassland) and temporal (yearly, decadal, and monthly) scales. (1) It was found that the temperature increased slightly (0.20℃/10a,P< 0.01). This may be due to the disappearance of the 0℃ isotherm, the rise of spring temperatures, and the increase of the isotherm (>15℃) area. At the same time, precipitation decreased significantly (-1.75mm/10a,P> 0.05). The expanding of the 400—600 mm precipitation contour area, the narrowing of ranges of 600—800 and 800—1000 mm, and excessive reduction of autumn precipitation indicated less precipitation. The sensitive areas of characteristic climatic change were mainly located between the latitudes 35—39°. The climate mutation occurred during 1991—1994. (2) Vegetation cover increased slightly in 55% of the total study area, with a rate of change of 0.00039/10a. Human activities could accelerate the changes of the vegetation cover and affect the rate of these changes. (3) Overall, the correlation between the vegetation cover and climatic factors is higher for a monthly than yearly scale. The correlation analysis between vegetation cover and climate change showed that annual vegetation cover was positively correlated with precipitation in grassland biomes; but it showed a more positive correlation with temperature in the agricultural and forest biomes. In addition, vegetation cover was temporally sensitive to precipitation. (4) The vegetation cover increased before and after the climate variations, but the rate of increase slowed down. Moreover, the grassland biome had an obvious response to the climatic variations, but the agricultural biome showed a significant response. The effect of human activity in regions of degradation was higher than that in environmentally improved areas. However, after the climate abruptly changing, the effect of human activity in environmentally improvement areas was higher than that in the degraded regions, and the influence of human activity will continue in the future.

climatic change; vegetation coverage; trend analysis; residual analysis; partial correlation

國家國際科技合作專項資助項目(2014DFA21620)

2015-07-30;

日期:2016-06-13

10.5846/stxb201507301600

* 通訊作者Corresponding author.E-mail: zhwenji1215@163.com

阿多,趙文吉,宮兆寧,張敏,范云豹.1981—2013華北平原氣候時空變化及其對植被覆蓋度的影響.生態(tài)學(xué)報,2017,37(2):576-592.

A D, Zhao W J, Gong Z N, Zhang M, Fan Y B.Temporal analysis of climate change and its relationship with vegetation cover on the north china plain from 1981 to 2013.Acta Ecologica Sinica,2017,37(2):576-592.