孫 琳, 唐國平, 竇乙峰, 劉永林, 吳裕珍

(中山大學(xué) 水資源與環(huán)境系, 廣東 廣州 510275)

人類社會(huì)經(jīng)濟(jì)活動(dòng)和自然環(huán)境變化是造成地表土地利用/覆被類型發(fā)生變化的主要原因[1]，這種變化不僅改變陸地表層物質(zhì)循環(huán)和生命過程，也影響區(qū)域生態(tài)環(huán)境的質(zhì)量和資源的可持續(xù)利用[1-2]。鑒于土地利用/覆被類型變化在生態(tài)環(huán)境保護(hù)和資源可持續(xù)利用方面的重要性，探討其時(shí)空變化特征[3-4]、成因以及機(jī)理[5]已成為全球環(huán)境變化研究的前沿與熱點(diǎn)。

中國“十三五”規(guī)劃支持珠江三角洲地區(qū)建設(shè)開放創(chuàng)新轉(zhuǎn)型升級(jí)新高地，即建設(shè)粵港澳大灣區(qū)。東江流域是珠三角和香港的重要飲用水源地，也是珠三角可持續(xù)發(fā)展戰(zhàn)略中重要的一環(huán)。流域社會(huì)經(jīng)濟(jì)的發(fā)展和自然環(huán)境的變化不僅改變區(qū)域土地利用/覆被類型的結(jié)構(gòu)，也影響流域生態(tài)環(huán)境的質(zhì)量和資源的可持續(xù)利用[6]。本文旨在利用土地利用/覆被類型遙感影像數(shù)據(jù)，結(jié)合土地利用評(píng)價(jià)指標(biāo)和灰色關(guān)聯(lián)分析，量化東江流域2001—2013年土地利用/覆被類型變化的時(shí)空特征和成因，區(qū)分人類活動(dòng)和自然環(huán)境變化對(duì)流域土地利用/覆被類型變化的貢獻(xiàn)率，進(jìn)而為優(yōu)化流域土地利用方式、保護(hù)流域生態(tài)環(huán)境和支持粵港澳大灣區(qū)的建設(shè)提供有益的信息和科技支撐。

1 研究資料與方法

1.1 研究區(qū)概況

東江為珠江水系主要干流之一，發(fā)源于江西省尋烏縣椏髻缽，自東北向西南流入廣東省境內(nèi)，是廣東目前綜合開發(fā)利用程度最高的水系[7]。流域面積35 340 km2，干流全長562 km[8]。東江上游位于江西省龍川縣以上的流域范圍，多為低山丘陵地貌，平均縱比降為2.21‰。龍川縣合河壩至觀音閣流域?yàn)闁|江中游，以丘陵和平原地貌為主，平均縱比降為0.31‰。觀音閣至入?？诹饔?yàn)闁|江下游，以臺(tái)地和沖積平原地貌為主，平均縱比降為0.17‰[9]。

東江流域(113°52′—115°52′E，22°38′—25°14′N)屬亞熱帶季風(fēng)性氣候，雨熱同期，降水主要集中于夏秋兩季，年平均降雨量1 750 mm，年平均氣溫為20.4 ℃[10]。流域內(nèi)大部分區(qū)域植被狀態(tài)良好，以常綠闊葉林、草被和人工營造的針葉林為主[11]，森林覆蓋度為68.77%，但流域內(nèi)土壤侵蝕等生態(tài)環(huán)境問題嚴(yán)重，年平均侵蝕總量高達(dá)1.62×109t，土壤侵蝕總面積為2 457.04 km2，占全流域面積的7%[12]。東江流域土壤類型主要為地帶性土壤，赤紅壤主要分布在海拔<300 m的低丘和山腳，紅壤多分布在300～500 m的低山丘陵，黃壤分布在500～900 m的中低山區(qū)[13]。東江流域經(jīng)濟(jì)結(jié)構(gòu)和人口分布空間差異明顯，上中游地區(qū)人口密度較小，以農(nóng)耕活動(dòng)為主，高新技術(shù)產(chǎn)業(yè)薄弱。下游地區(qū)人口密度大，經(jīng)濟(jì)發(fā)展迅速，廣州、深圳和東莞等地成為廣東省經(jīng)濟(jì)最為發(fā)達(dá)的地區(qū)。

1.2 數(shù)據(jù)來源和處理

本文選用2001—2013年MODIS土地覆被類型遙感影像數(shù)據(jù)(MCD12Q1)馬里蘭大學(xué)(UMD格式)分類方案，數(shù)據(jù)空間分辨率為500 m，質(zhì)量控制(QC)數(shù)據(jù)范圍為0～97，說明像元已經(jīng)進(jìn)行反演，質(zhì)量優(yōu)。本文使用人機(jī)交互的方式提取東江流域土地利用/覆被類型遙感數(shù)據(jù)，參照中國科學(xué)院土地利用遙感監(jiān)測(cè)分類體系[14]以及實(shí)際樣本資料，將東江流域土地利用/覆被類型劃分為水域、常綠闊葉林、混合林、灌叢、草地、耕地和城鎮(zhèn)用地7類，并統(tǒng)一為WGS84坐標(biāo)系。選取245個(gè)典型特征點(diǎn)野外調(diào)查數(shù)據(jù)，結(jié)合廣東省、江西省歷年土地利用專題圖件對(duì)影像解譯結(jié)果進(jìn)行驗(yàn)證，解譯精度高達(dá)84.5%。

1.3 土地利用/覆被類型變化評(píng)價(jià)指標(biāo)及方法

本文采用土地利用類型變化量、變化率以及土地利用綜合指數(shù)3個(gè)指標(biāo)來定量分析2001—2013年東江流域土地利用/覆被類型變化的時(shí)空特征。

土地利用/覆被類型類型變化量(UT)表示某土地利用/覆被類型類型相對(duì)研究初期面積數(shù)量變化的大小，土地利用/覆被類型類型變化率(KT)表示研究末期土地利用/覆被類型類型面積相對(duì)研究初期浮動(dòng)的程度。計(jì)算公式為：

UT=Ub-Ua

(1)

(2)

式中：T——研究歷時(shí)，通常取年為單位；UT——時(shí)間序列T內(nèi)某種土地利用/覆被類型的面積變化量(km2)；KT——時(shí)間序列T內(nèi)某種土地利用/覆被類型的面積變化率(%)；Ua，Ub——研究期初與研究期末某種土地利用/覆被類型的面積數(shù)量(km2)。

土地利用綜合指數(shù)(La)表示研究區(qū)土地利用實(shí)際開發(fā)程度的高低，計(jì)算方法如下[14]：

(3)

式中：La——區(qū)域土地利用綜合指數(shù)；n——土地利用/覆被類型總數(shù)；Ai——第i種土地利用/覆被類型對(duì)應(yīng)的分級(jí)指數(shù)(見表1)；Ci——第i類土地利用/覆被類型占研究區(qū)總面積的百分比。

表1 土地利用程度分級(jí)賦值

此外，本文利用轉(zhuǎn)移矩陣來量化各土地利用/覆被類型間相互轉(zhuǎn)化的數(shù)量特征以及轉(zhuǎn)化趨勢(shì)，其基本形式如下[15]：

(4)

式中：Sij——研究期內(nèi)第i類土地利用/覆被類型轉(zhuǎn)化成第j類土地利用/覆被類型的面積(km2)；n——區(qū)域內(nèi)土地利用/覆被類型的總數(shù)。

1.4 分析土地利用/覆被類型變化的驅(qū)動(dòng)力

基于查閱廣東省、江西省各市縣統(tǒng)計(jì)年鑒和閱讀文獻(xiàn)獲得的2001—2013年東江流域社會(huì)經(jīng)濟(jì)等相關(guān)數(shù)據(jù)，本文定義了4類驅(qū)動(dòng)因子： ①人口(常住人口數(shù))； ②經(jīng)濟(jì)發(fā)展水平(地方生產(chǎn)總值、社會(huì)消費(fèi)品零售總額)和結(jié)構(gòu)(農(nóng)業(yè)總產(chǎn)值、林業(yè)總產(chǎn)值、規(guī)模以上工業(yè)總產(chǎn)值)； ③城市化建設(shè)程度(農(nóng)作物播種面積、施工建筑面積)； ④技術(shù)革新(農(nóng)業(yè)機(jī)械總動(dòng)力)?；谏鲜?類因子，本文采用灰色關(guān)聯(lián)分析法量化各類驅(qū)動(dòng)因子與各土地利用/覆被類型面積變化的關(guān)聯(lián)度。

灰色關(guān)聯(lián)分析的基本步驟如下： ①建立參考數(shù)列和比較數(shù)列。參考數(shù)列為各土地利用/覆被類型的面積數(shù)據(jù)，比較數(shù)列為經(jīng)濟(jì)指標(biāo)(驅(qū)動(dòng)因子)； ②利用均值化變換方法(公式5)對(duì)參考數(shù)列和比較數(shù)列進(jìn)行歸一化處理； ③計(jì)算關(guān)聯(lián)系數(shù)。首先計(jì)算歸一化處理后比較數(shù)列與參考數(shù)列差的絕對(duì)值，將差列絕對(duì)值的最大值(Δmax)和最小值(Δmin)代入公式(6)求得關(guān)聯(lián)系數(shù)； ④計(jì)算和排列關(guān)聯(lián)度。同一驅(qū)動(dòng)因子不同年份對(duì)應(yīng)的關(guān)聯(lián)系數(shù)平均值即為關(guān)聯(lián)度。關(guān)聯(lián)度大的驅(qū)動(dòng)因子即可看成是影響對(duì)應(yīng)土地類型變化的主要因素，反之亦然。

(5)

式中：yi(k)——第i級(jí)參考(比較)數(shù)列的第k項(xiàng)驅(qū)動(dòng)因子的歸一化數(shù)據(jù)；xi(k)——第i級(jí)參考(比較)數(shù)列的第k項(xiàng)驅(qū)動(dòng)因子原數(shù)據(jù)；m——每項(xiàng)驅(qū)動(dòng)因子對(duì)應(yīng)的數(shù)值個(gè)數(shù)。

(6)

式中：Li(k)——第i級(jí)比較數(shù)列第k項(xiàng)驅(qū)動(dòng)因子對(duì)應(yīng)的關(guān)聯(lián)系數(shù)； Δmax，Δmin——差列絕對(duì)值的最大值和最小值； Δi(k)——第i級(jí)比較數(shù)列第k項(xiàng)驅(qū)動(dòng)因子的歸一化處理后數(shù)值。

2 結(jié)果與分析

2.1 東江流域土地利用/覆被類型變化的時(shí)空模式

圖1展示了2001—2013年東江流域各土地利用/覆被類型面積的變化趨勢(shì)。結(jié)果表明，常綠闊葉林和城鎮(zhèn)用地呈現(xiàn)顯著遞增趨勢(shì)(p<0.02；圖1a,1f)，灌叢和草地呈現(xiàn)顯著遞減趨勢(shì)(p<0.01；圖1c,1d)?；旌狭趾透匾渤尸F(xiàn)遞減趨勢(shì)，但該趨勢(shì)統(tǒng)計(jì)上不顯著(p>0.10；圖1b,1e)。

流域內(nèi)常綠闊葉林在流域內(nèi)分布最為廣泛，并且隨時(shí)間推移，其優(yōu)勢(shì)地位日益凸顯。研究初期，常綠闊葉林主要分布在東江中游的東源市、新豐縣和連平縣，總面積占全流域的33.7%；隨時(shí)間推移，其分布范圍以新豐江水庫為中心呈放射狀增長；到研究末期，常綠闊葉林在全流域均有分布，覆被面積高達(dá)53.7%，以下游和中游面積增加最明顯?；旌狭盅芯砍跗诩s占全流域的12.4%，2001—2004年面積急劇下降，2004年后在3個(gè)流域變化相對(duì)平緩，至研究末期約占總流域面積4.6%，主要分布在東江源區(qū)的安遠(yuǎn)縣、尋烏縣以及中游的連平縣、和平縣(圖2b)。灌叢與草地的變化類似，草地面積2006年出現(xiàn)短暫的峰值(約占全流域的34.2%)，總體逐年遞減，至2013年下降到19.6%，草地面積遞減主要發(fā)生在東江中游和下游(圖2d)。耕地面積在2005年達(dá)到峰值，約占全流域面積的20.9%，但2012年減少至13.4%，主要發(fā)生在人口密集、地勢(shì)平坦的下游流域(圖2e)。城鎮(zhèn)用地面積下游>中游>上游，狀態(tài)比較穩(wěn)定。研究末期在東江下游發(fā)生顯著增長(圖2f)，隨著時(shí)間推移，這種空間的差異性將呈增大趨勢(shì)。

圖1 2001-2013年東江流域各土地利用/覆被類型面積的變化趨勢(shì)

2.2 流域各土地利用/覆被類型變化的程度分析

如表2所示，不同土地利用/覆被類型面積變化量可以看出，常綠闊葉林與草地、混合林呈此消彼長關(guān)系，上游面積變化量表現(xiàn)為：常綠闊葉林>混合林>草地，中下游面積變化量為：常綠闊葉林>草地>混合林。

流域內(nèi)灌叢和耕地面積也發(fā)生較大變化，耕地面積變化主要發(fā)生在下游(詳見圖2e)，灌叢面積變化主要在中下游。城鎮(zhèn)用地面積變化相對(duì)穩(wěn)定，僅下游略有增長。水域面積在中游呈增長趨勢(shì)，但在上游和下游遞減。

表2 東江流域2001-2013年土地利用/覆被類型變化情況

從變化率上看，常綠闊葉林在全流域變化率為59.4%，是增長速度最快的類型，增長速率為：上游>下游>中游。灌叢為流域內(nèi)面積減少最劇烈的類型，變化率高達(dá)-63.6%，上、中游近90%的灌叢消失，下游50.7%的灌叢轉(zhuǎn)化成了其他土地類型。與灌叢變化相似，混合林、草地在上中下游同樣呈現(xiàn)明顯的遞減趨勢(shì)(圖2c,2d)，以下游變化最為突出，混合林、草地分別減少75.0%和49.1%。耕地在全流域面積增長了6.8%，變化速率是下游>上游>中游。

流域土地利用程度綜合指數(shù)為224.9，維持在中游水平，說明流域開發(fā)程度位于發(fā)展期，仍具有充裕的上升空間，隨著時(shí)間的推移，受珠江三角洲經(jīng)濟(jì)群的輻射影響加強(qiáng)，流域土地開發(fā)利用如城市建設(shè)活動(dòng)將更加頻繁。

圖2 2001-2013年東江上中下游各土地利用/覆被類型面積的變化趨勢(shì)

2.3 各土地利用/覆被類型之間的轉(zhuǎn)化關(guān)系

如表3所示，研究期間流域內(nèi)常綠闊葉林面積劇增了59.4%(約6 661.0 km2)，其增加面積主要來源于草地(36.1%)和混合林(18.6%)?；旌狭洲D(zhuǎn)化量是上游>中游>下游，比例依次為43.1%，32.0%和27.0%，草地轉(zhuǎn)化量則是下游>中游>上游，比例依次為52.7%，60.1%和64.0%?；旌狭趾筒莸剞D(zhuǎn)出面積除了大部分轉(zhuǎn)化成常綠闊葉林外，二者之間的相互轉(zhuǎn)化以及轉(zhuǎn)化成耕地的比重也較大。耕地面積的增加主要來自于草地，但由于存在相互轉(zhuǎn)化，僅517 km2草地完全轉(zhuǎn)化為耕地。城鎮(zhèn)用地在上、中游向其他類型用地轉(zhuǎn)化不明顯，僅下游區(qū)域面積增長2.2%(28.8 km2)，其中50.4%來自于灌叢，48.7%來源于耕地。灌叢的轉(zhuǎn)化相對(duì)平均，面積減少中的20.7%轉(zhuǎn)化為常綠闊葉林，15.5%轉(zhuǎn)化為耕地，14.2%轉(zhuǎn)化為草地，10.3%轉(zhuǎn)化為混合林。

表3 2001-2013年東江流域土地利用/覆被類型轉(zhuǎn)移矩陣 km2

2.4 東江流域土地利用/覆被類型變化的驅(qū)動(dòng)力

圖3分別展示了灰色關(guān)聯(lián)分析在上中下游地區(qū)分析的結(jié)果。如圖3所示，東江流域2001—2013年上中下游常綠闊葉林變化主導(dǎo)驅(qū)動(dòng)力差異明顯。上游，常綠闊葉林變化與社會(huì)消費(fèi)品零售總額(0.93)、林業(yè)總產(chǎn)值(0.95)以及農(nóng)業(yè)機(jī)械化總動(dòng)力(0.95)等因子高度關(guān)聯(lián)(圖3a)，經(jīng)濟(jì)發(fā)展和技術(shù)革新為主要驅(qū)動(dòng)力。

注：人口驅(qū)動(dòng)力包括PRP(常住人口)； 經(jīng)濟(jì)驅(qū)動(dòng)力包括GDP(地方性生產(chǎn)總值)、SCG(社會(huì)消費(fèi)品零售總額)、APV(農(nóng)業(yè)總產(chǎn)值)、FPV(林業(yè)總產(chǎn)值)、GPV(規(guī)模以上工業(yè)總產(chǎn)值)； 城市化發(fā)展驅(qū)動(dòng)力包括CPA(農(nóng)作物播種面積)、BCA(建筑施工面積)； 技術(shù)驅(qū)動(dòng)力包括PAM(農(nóng)業(yè)機(jī)械化總動(dòng)力)。

中下游常綠闊葉林變化與各類驅(qū)動(dòng)因子均聯(lián)系密切，各驅(qū)動(dòng)力相輔相成，共同影響常綠闊葉林的變化。由于城鎮(zhèn)用地僅在下游發(fā)生改變，故僅對(duì)下游城鎮(zhèn)用地變化進(jìn)行灰色關(guān)聯(lián)分析。下游與城鎮(zhèn)用地關(guān)聯(lián)度較高的因子為農(nóng)作物播種面積、畜牧業(yè)總產(chǎn)值以及常住人口，關(guān)聯(lián)度依次為0.97，0.89，0.81(圖3f)，表明城市化建設(shè)、技術(shù)革新以及人口因素共同促進(jìn)城鎮(zhèn)用地的轉(zhuǎn)化。除常綠闊葉林和城鎮(zhèn)用地外，其他土地類型在3個(gè)地區(qū)內(nèi)受不同人為驅(qū)動(dòng)力的影響程度基本一致，與常住人口、農(nóng)作物播種面積以及建筑施工面積呈現(xiàn)顯著關(guān)聯(lián)的特征(圖3b—3e)，關(guān)聯(lián)系數(shù)平均為0.83，0.85，0.74?？傮w來說，按關(guān)聯(lián)度大小排序，影響東江流域土地利用/覆被類型變化的人為驅(qū)動(dòng)力依次為城市化建設(shè)、人口增長、技術(shù)革新以及經(jīng)濟(jì)發(fā)展。

3 討 論

3.1 東江流域土地利用/覆被類型變化的時(shí)空差異和原因

東江流域土地利用/覆被類型變化的時(shí)空特征與流域自然環(huán)境條件密切相關(guān)[7]。東江上游90%的面積為山地和丘陵，降水多且集中于夏秋兩季，易引發(fā)水土流失、山體滑坡等自然災(zāi)害。地形和自然災(zāi)害的易發(fā)性限制了上游的土地開發(fā)利用。中下游鄰近珠江三角洲，地勢(shì)平緩，土地具有較高的可塑性，利于開墾改造，城市化用地呈明顯增加趨勢(shì)。自然環(huán)境條件的空間分異特征與上游至下游耕地、城鎮(zhèn)用地面積變化的空間分布特征相吻合，如耕地面積比例(10.0%，10.2%和26.3%)與城鎮(zhèn)用地面積比例(0.1%，0.6%和8.6%)的分布為：下游>中游>上游，表明自然環(huán)境條件與東江流域土地利用/覆被類型空間分異特征顯著關(guān)聯(lián)。

另一方面，東江流域土地利用/覆被類型變化的時(shí)空特征也反映了人口增加、經(jīng)濟(jì)增長、區(qū)域規(guī)劃政策等社會(huì)因素的影響[7]。上中下游人口和經(jīng)濟(jì)發(fā)展程度下游>中游>上游。人口和經(jīng)濟(jì)發(fā)展越快，耕地向城鎮(zhèn)用地、草地向耕地轉(zhuǎn)化程度越顯著，與下游城鎮(zhèn)用地顯著增加，上游耕地減少明顯趨勢(shì)相一致，說明人口與經(jīng)濟(jì)發(fā)展水平的差異是造成城鎮(zhèn)用地、草地等在不同流域存在明顯空間差異的原因。

此外，政府決策的制定與實(shí)施也是造成東江流域土地利用/覆被類型變化具有時(shí)空異質(zhì)性的重要原因。研究前期，上游流域經(jīng)濟(jì)發(fā)展水平相對(duì)落后，地方財(cái)政收入主要來源于礦產(chǎn)資源開采及林業(yè)，對(duì)林地草地破壞嚴(yán)重。為保護(hù)東江源區(qū)生態(tài)環(huán)境，2003年江西省通過《關(guān)于加強(qiáng)東江源區(qū)生態(tài)環(huán)境保護(hù)和建設(shè)的決定》，要求加大森林資源管理和培育力度[16]。為改變消耗環(huán)境成本的經(jīng)濟(jì)模式，贛州政府提出打造“臍橙基地”的產(chǎn)業(yè)定位。從2005年起，上游全面停止采伐天然闊葉樹作為商品材。這種政策的實(shí)施與上游常綠闊葉林自2006年后面積逐年遞增直接相關(guān)(圖1a)。東江中游在研究初期積極響應(yīng)國家“十五”計(jì)劃，對(duì)農(nóng)業(yè)結(jié)構(gòu)進(jìn)行戰(zhàn)略性調(diào)整，提倡大力發(fā)展高效生態(tài)循環(huán)農(nóng)業(yè)和機(jī)械現(xiàn)代農(nóng)業(yè)。2005年，東江中游遭遇百年一遇洪澇災(zāi)害，政府加強(qiáng)森林資源保護(hù)，推進(jìn)東江、韓江水源林及珠江防護(hù)林的建設(shè)，最終導(dǎo)致林地面積2006年后逐年遞增，耕地面積2005—2011年逐年遞減(圖2e)。東江下游靠近珠江三角洲，自20世紀(jì)90年代改革開放以來，廣州的快速發(fā)展和深圳特區(qū)的迅速崛起帶動(dòng)惠州、東莞等周邊城市的經(jīng)濟(jì)發(fā)展，推動(dòng)下游城市化建設(shè)，這與下游城鎮(zhèn)用地面積在21世紀(jì)初明顯增長相呼應(yīng)(圖2f)。

3.2 東江流域土地利用/覆被類型變化的社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)力

3.2.1 人口增長與城市化發(fā)展 東江流域鄰近經(jīng)濟(jì)發(fā)達(dá)的珠江三角洲，其開放的市場(chǎng)環(huán)境、多元的文化交流、龐大的人才需求吸引大量的外來人口，常住人口由2001年1 109.1萬人增長至2013年2 214.6萬人，城鎮(zhèn)用地面積由2001年的1 293.8 km2增長至2013年的1 322.5 km2。人口數(shù)量的增加和城市化進(jìn)程的推進(jìn)，一方面導(dǎo)致糧食和畜牧產(chǎn)品需求壓力增大，為應(yīng)對(duì)不斷增長的需求，增大耕地面積，人們開荒、毀林、毀草造田，草地和灌叢面積持續(xù)遞減；另一方面，為滿足日益增長的住房和娛樂服務(wù)需求，非農(nóng)業(yè)用地對(duì)土地的需求急劇增長，大量耕地轉(zhuǎn)化為城鎮(zhèn)用地，加快了農(nóng)業(yè)用地的轉(zhuǎn)化[17]。

3.2.2 技術(shù)革新與社會(huì)經(jīng)濟(jì)進(jìn)步 隨著經(jīng)濟(jì)發(fā)展、農(nóng)業(yè)機(jī)械化水平提高以及城鎮(zhèn)用地的擴(kuò)張[16]，技術(shù)驅(qū)動(dòng)力和經(jīng)濟(jì)驅(qū)動(dòng)力的影響比重逐漸加大。東江流域農(nóng)作物播種面積由2001年的6.91×105km2減少至2013年的6.78×105km2。為彌補(bǔ)播種面積減少帶來的農(nóng)作物產(chǎn)量下降，東江流域促進(jìn)農(nóng)民合作社的發(fā)展和現(xiàn)代化農(nóng)業(yè)科學(xué)技術(shù)的利用，加強(qiáng)流域土地集約程度，提高耕地單產(chǎn)水平，實(shí)現(xiàn)以較少的耕地面積滿足居民的日常需求。東江流域2013年農(nóng)業(yè)總產(chǎn)值較2001年增加了96.4億元，農(nóng)業(yè)機(jī)械化總動(dòng)力高達(dá)2.81×1010kWh。在基本生活需求得以滿足的前提下，依靠市場(chǎng)調(diào)控、政府引導(dǎo)，東江流域農(nóng)牧業(yè)結(jié)構(gòu)開始優(yōu)化升級(jí)。到2013年，第一、第二產(chǎn)業(yè)在地方性經(jīng)濟(jì)中的比重分別減少了6.6%和5.3%，第三產(chǎn)業(yè)卻增長了11.8%。這種產(chǎn)業(yè)結(jié)構(gòu)的內(nèi)部調(diào)整，促進(jìn)各流域土地利用/覆被類型開始向其他方向轉(zhuǎn)化，如農(nóng)林用地轉(zhuǎn)為非農(nóng)用地，這也部分解釋了常綠闊葉林持續(xù)增長而耕地逐漸遞減的現(xiàn)象[7]。

3.2.3 政府決策 21世紀(jì)初，珠江三角洲已經(jīng)發(fā)展成全國市場(chǎng)化程度最高，市場(chǎng)體系最完備的地區(qū)。東江流域受市場(chǎng)經(jīng)濟(jì)體制的影響，現(xiàn)代化進(jìn)程不斷加快，土地資源逐漸進(jìn)入市場(chǎng)支配，資產(chǎn)價(jià)值日益凸顯。在流域相關(guān)土地開發(fā)政策及發(fā)展規(guī)劃的引導(dǎo)下，土地向其他高收益類型發(fā)生轉(zhuǎn)化[18]。為統(tǒng)籌東江流域的開發(fā)建設(shè)工程及生態(tài)環(huán)境保護(hù)工作，2007年廣東正式成立東江流域管理局，負(fù)責(zé)流域綜合規(guī)劃、工程管理和水土保持等政策的執(zhí)行與監(jiān)督。同時(shí)，廣東省積極響應(yīng)國家主體功能區(qū)規(guī)劃政策，遵循不同國土空間的自然屬性，著力構(gòu)建“五大戰(zhàn)略格局”。其中，東江下游廣州市、深圳市、惠州市和東莞市實(shí)行嚴(yán)格的建設(shè)用地增量控制[19]，城市發(fā)展由外延擴(kuò)張式轉(zhuǎn)為以節(jié)約集約為主，積極探索城市土地立體化利用模式；中、上游韶關(guān)市、河源市及梅州市為北部生態(tài)經(jīng)濟(jì)區(qū)，嚴(yán)格控制開發(fā)力度，限制城鎮(zhèn)發(fā)展用地和農(nóng)村居民點(diǎn)用地，對(duì)成片耕地進(jìn)行標(biāo)準(zhǔn)化基本農(nóng)田保護(hù)區(qū)建設(shè)。這也是為什么常綠闊葉林面積2007年后在上中下游都持續(xù)增長，但灌叢、草地和耕地面積卻持續(xù)減少，尤其在中、下游變化更明顯。此外，2003年江西省通過《關(guān)于加強(qiáng)東江源區(qū)生態(tài)環(huán)境保護(hù)和建設(shè)的決定》并成立東江源國家級(jí)生態(tài)功能保護(hù)區(qū)，對(duì)源區(qū)進(jìn)行統(tǒng)籌規(guī)劃，先后實(shí)施退耕還林、珠江防護(hù)林工程等一系列項(xiàng)目。這些生態(tài)環(huán)境保護(hù)和重大建設(shè)工程的實(shí)施明顯提高了東江流域森林覆蓋的面積，改善了流域生態(tài)環(huán)境和水資源狀態(tài)。

4 結(jié) 論

(1) 2001—2013年東江流域土地利用/覆被類型的分布格局變化顯著。土地利用/覆被類型以常綠闊葉林、草地和耕地為主，但隨著時(shí)間推移，常綠闊葉林的主導(dǎo)地位日益凸顯，并以新豐江水庫為中心呈放射狀增長。

(2) 研究期內(nèi)不同土地利用/覆被類型變化差異明顯。水域、耕地和城鎮(zhèn)用地變化輕微，灌叢、混合林和常綠闊葉林面積變化劇烈，變化率均超過50%。盡管如此，由于灌叢、混合林、常綠闊葉林和草地4種類型之間相互轉(zhuǎn)化頻繁，它們的總面積變化不大。

(3) 土地利用/覆被類型變化的區(qū)域差異明顯。上游耕地面積減少顯著，主要向常綠闊葉林和草地轉(zhuǎn)化；相反，中、下游流域總體表現(xiàn)混合林和草地減少明顯，耕地面積增加，并且伴隨人口增長和城市化水平的提高，耕地面積向城鎮(zhèn)用地轉(zhuǎn)化幅度逐漸增強(qiáng)。

(4) 人口增長、城市化建設(shè)、技術(shù)革新和政府決策是造成研究期間內(nèi)東江流域土地利用/覆被類型發(fā)生變化的主要驅(qū)動(dòng)力。但是，各類驅(qū)動(dòng)力并不是孤立地發(fā)揮作用，而是相輔相成共同驅(qū)使流域內(nèi)土地利用/覆被類型發(fā)生變化。

總體來說，最近10 a多東江流域的土地利用/覆被類型變化有利于流域生態(tài)環(huán)境和水資源的保護(hù)，進(jìn)而有助于粵港澳大灣區(qū)的建設(shè)。