吳坤，王文杰 ，劉軍會，張哲，胡圣武

1.湖南科技大學(xué)建筑與城鄉(xiāng)規(guī)劃學(xué)院，湖南 湘潭 411201

2.中國環(huán)境科學(xué)研究院環(huán)境信息研究所，北京 100012

3.河南理工大學(xué)測繪與國土信息工程學(xué)院，河南 焦作 454000

土地利用變化是人類活動作用于生態(tài)環(huán)境所呈現(xiàn)出的最顯著形式［1］，作為土地科學(xué)計劃研究的重要課題，已經(jīng)成為全球環(huán)境變化與可持續(xù)發(fā)展研究的重要內(nèi)容［2-4］。對土地利用過程、格局、驅(qū)動力的研究被認為是揭示土地利用時空變化規(guī)律的有效方法［3-6］，國內(nèi)外學(xué)者關(guān)于土地利用動態(tài)變化及其驅(qū)動力分析進行了大量研究，以B. L. Turner 等［6-9］為代表的北美學(xué)者通常從宏觀角度出發(fā)，采用對比案例的研究方法將土地利用變化與土地過程管理進行分析和建模，定性研究全球土地利用變化與環(huán)境變化的關(guān)系;歐洲學(xué)者［10］多將土地利用作為整體系統(tǒng)，對不同時空土地利用變化仿真，以福利分析為出發(fā)點，研究社會經(jīng)濟驅(qū)動因素;日本學(xué)者［11］主要將經(jīng)濟學(xué)模型和數(shù)量模型結(jié)合，采用相關(guān)分析方法選取說明變量，采用KSIM 等模擬方法對土地利用變化及政策環(huán)境進行模擬，定量研究區(qū)域土地利用變化過程;國內(nèi)學(xué)者則大多基于RS 和GIS 技術(shù)，以人地矛盾突出、經(jīng)濟發(fā)展迅速的區(qū)域為研究區(qū)，采用土地利用程度變化模型［12］、驅(qū)動力診斷模型［13］、土地利用結(jié)構(gòu)預(yù)測模型等模型與數(shù)理統(tǒng)計分析方法相結(jié)合的技術(shù)手段研究土地利用變化及其驅(qū)動因素。

全國主體功能區(qū)規(guī)劃［14］將成渝經(jīng)濟區(qū)列入國家層面的重點開發(fā)區(qū)域，近年來社會經(jīng)濟發(fā)展十分迅速，土地格局變化劇烈。目前對成渝經(jīng)濟區(qū)土地利用的研究主要有以下幾個方面:土地利用與生態(tài)服務(wù)價值動態(tài)分析［15］、城市土地集約利用評價及時空特征分析［16］等。而對于西部大開發(fā)計劃實施后的第一個10年，針對成渝經(jīng)濟區(qū)尺度土地利用變化特征及驅(qū)動力分析方面的研究還較少。筆者利用遙感技術(shù)(RS)和地理信息系統(tǒng)(GIS)結(jié)合的技術(shù)手段，分析2000—2010年成渝經(jīng)濟區(qū)土地利用數(shù)量變化特點，識別土地利用類型轉(zhuǎn)移規(guī)律及空間分布特征，并構(gòu)建偏最小二乘回歸模型對成渝經(jīng)濟區(qū)土地利用變化進行驅(qū)動力分析，以期為該區(qū)域人-地關(guān)系協(xié)調(diào)發(fā)展及土地利用科學(xué)管理提供理論依據(jù)。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

成渝經(jīng)濟區(qū)是我國典型的流域經(jīng)濟區(qū)，地處四川盆地，北接陜甘、南連云貴、西通青藏、東鄰湘鄂。地理位置上包括成都和重慶兩大城市圈，其中四川省含成都、德陽、綿陽、南充、達州、廣安、遂寧、資陽、內(nèi)江、眉山、自貢、樂山、宜賓、瀘州、雅安等15 個地市;重慶市含中心城區(qū)(9 個區(qū))、璧山、永川、大足、合川、江津、綦江、南川、長壽、涪陵、萬州、榮昌、潼南、銅梁、墊江、豐都、石柱、忠縣、梁平、云陽、開縣等29 個區(qū)縣(圖1)。研究區(qū)總面積20.86 ×104km2，其中，四川省面積為15.73 ×104km2，重慶市面積為5.13 ×104km2。成渝經(jīng)濟區(qū)生態(tài)系統(tǒng)復(fù)雜，自然資源豐富，是長江上游生物多樣性的寶庫，成渝經(jīng)濟區(qū)重要的地理位置和生態(tài)功能決定了其在長江流域生態(tài)屏障建設(shè)中的重要地位。

圖1 成渝經(jīng)濟區(qū)地理位置Fig.1 Location of Chengdu-Chongqing Economic Zone

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

土地利用數(shù)據(jù)來源于全國生態(tài)環(huán)境10年變化(2000—2010)遙感調(diào)查與評估專項:長江流域生態(tài)環(huán)境10年變化調(diào)查，一級為IPCC(政府間氣候變化專門委員會)土地覆被類型，二級基于碳收支的土地覆被類型。利用ArcGIS10.1 空間分析模塊，對土地利用的空間變化進行統(tǒng)計，得到研究區(qū)土地利用的轉(zhuǎn)移矩陣，用以分析不同土地利用類型間的轉(zhuǎn)換規(guī)律。分析還用到了2000年、2005年、2010年重慶市統(tǒng)計年鑒和四川省統(tǒng)計年鑒的社會經(jīng)濟統(tǒng)計數(shù)據(jù)。

1.3 分析方法

1.3.1 土地利用變化特征分析方法

土地利用的變化情況可以通過單一土地利用類型動態(tài)度、轉(zhuǎn)出率和轉(zhuǎn)入率來表示。

(1)土地利用類型動態(tài)度

單一土地利用類型動態(tài)度可以定量表達研究區(qū)一定時間范圍內(nèi)的土地利用類型數(shù)量的變化情況［17-18］，綜合土地利用動態(tài)度描述土地利用的變化速度，其表達式分別為:

式中:Ki為研究時間范圍內(nèi)第i 種單一土地利用類型的動態(tài)度;Ua為研究初期a 時期第i 種土地利用類型的數(shù)量;Ub為研究末期b 時期第i 種土地利用類型的數(shù)量;T 為研究時段，當T 單位為年時，Ki即為i 種土地類型的年變化率;Lc為綜合土地利用動態(tài)度，即土地利用年變化率;ΔLUi→j為T 時間段內(nèi)第i 種土地類型轉(zhuǎn)化其他利用類型的總面積;ΔLUi為監(jiān)測起始時間第i 種土地類型的面積。

(2)轉(zhuǎn)出率(α)與轉(zhuǎn)入率(β)

式中:αij為a 時期第i 種土地利用類型轉(zhuǎn)化為b 時期第j 種土地利用類型的比例;βij為b 時期第j 種土地利用類型從a 時期第i 種土地利用類型轉(zhuǎn)化而來的比例［19］。

(3)土地利用轉(zhuǎn)移矩陣

土地利用轉(zhuǎn)移矩陣是定量研究土地系統(tǒng)內(nèi)各種土地利用類型間相互轉(zhuǎn)化主要方法，能夠反映出土地利用變化的結(jié)構(gòu)特征以及不同類型土地的轉(zhuǎn)移方向［20］。其表達式為:

式中:Sij指a 時期第i 種土地利用類型轉(zhuǎn)化為b 時期第j 種土地利用類型的面積;n 為土地利用類型的數(shù)量。

1.3.2 驅(qū)動力定量分析方法

偏最小二乘回歸通過自變量主成分的提取，很好地概括自變量系統(tǒng)，在解釋因變量的同時排除系統(tǒng)存在的噪聲干擾，消除變量間的多重共線性，并且能在樣本容量小于變量個數(shù)的情況下進行回歸建模，適合用于土地利用變化驅(qū)動分析［21-22］。其基本算法如下［22-23］:

假設(shè)自變量(X)和因變量(Y)經(jīng)過數(shù)據(jù)標準化后分別得到數(shù)據(jù)矩陣E0和F0，記t1和u1分別是從E0和F0中提取的第一主成分，則有:

其中:w1為E0的第一主軸，且有‖w1‖=1

其中:v1為F0的第一主軸，且有‖v1‖=1

根據(jù)主成分分析原理，要求t1和u1的標準差趨于最大，由典型相關(guān)分析原理，t1和u1的協(xié)方差要達到最大，綜合起來即:

數(shù)學(xué)表述為求解優(yōu)化問題，即在‖w1‖ =1 和‖v1‖=1 的約束條件下求的最大值。求得w1和v1后即可獲得第一主成分t1和u1，再分別實施E0對t1的回歸和F0對u1的回歸，如果回歸方程滿足精度要求，則算法終止;否則，利用殘余信息繼續(xù)下一輪成分提取，直至滿足精度為止。如最終提取r 個主成分T=(t1，…，tr)，則偏最小二乘回歸首先建立Y 與T 的回歸，然后再表達成Y 與原始自變量X 的回歸方程。

2 土地利用變化特征

2.1 土地利用數(shù)量變化

根據(jù)2000年、2005年、2010年3 個時期遙感影像解譯數(shù)據(jù)，成渝經(jīng)濟區(qū)主要土地利用類型為耕地和林地(圖2)，對比分析3 個時期的土地利用數(shù)據(jù)(表1)，發(fā)現(xiàn)土地利用變化的特征如下:

(1)2000—2010年，林地面積增加2 719.15 km2，所占比例增加1.3 個百分點;草地面積增加69.44 km2，所占比例增加0.03 個百分點;濕地面積增加594.08 km2，所占比例增加0.28 個百分點;耕地面積減少明顯，減少面積達5 840.27 km2，所占比例下降2.8 個百分點;人工表面用地面積迅速擴大，增加了2 337.03 km2，所占比例增加1.12 個百分點;其他用地面積增大120.56 km2。

(2)2000—2010年土地利用動態(tài)度由高到低依次為人工表面＞濕地＞其他用地＞耕地＞林地＞草地。人工表面土地利用類型變化最為劇烈，2000—2010年動態(tài)度為8.35%;2000—2005年和2005—2010年耕地變化動態(tài)度均為-0.57%，耕地總面積呈現(xiàn)出以穩(wěn)定速度持續(xù)降低態(tài)勢。2000—2005年研究區(qū)綜合土地利用動態(tài)度為0.170 2%，2005—2010年綜合土地利用動態(tài)度為0.170 7%，表明成渝經(jīng)濟區(qū)2005—2010年土地利用變化速度大于2000—2005年，呈現(xiàn)加快變化趨勢。

(3)從人工表面內(nèi)部結(jié)構(gòu)分析，居民地、工業(yè)用地、交通用地和采礦用地面積均大幅度提升。其中，2000—2010年工業(yè)用地變化尤為劇烈，動態(tài)度達到302.19%，說明成渝經(jīng)濟區(qū)工業(yè)產(chǎn)業(yè)發(fā)展正處于十分活躍的態(tài)勢。從時間上看，2005—2010年工業(yè)用地和采礦用地發(fā)展速度較2000—2005年加快，而居住地和交通用地動態(tài)度的變化則相對變緩，說明成渝經(jīng)濟區(qū)人工表面用地中用于工業(yè)用地發(fā)展的投入較大。

圖2 2000年、2005年和2010年成渝經(jīng)濟區(qū)土地利用Fig.2 Land use map of Chengdu-Chongqing Economic Zone in 2000，2005，2010

表1 2000—2010年成渝經(jīng)濟區(qū)土地利用數(shù)量變化Table 1 Land use change from 2000 to 2010 in Chengdu-Chongqing Economic Zone

2.2 土地利用類型轉(zhuǎn)移

2000—2010年土地利用轉(zhuǎn)移矩陣表明，成渝經(jīng)濟區(qū)內(nèi)各土地利用類型之間轉(zhuǎn)換頻繁，其中，大面積和大比例的土地類型轉(zhuǎn)換主要發(fā)生在耕地、人工表面、林地和濕地之間(表2)。

耕地主要轉(zhuǎn)化為人工表面用地和林地。10年間共有5 907.15 km2耕地轉(zhuǎn)出為其他幾種土地利用類型，其中，3.58%的耕地轉(zhuǎn)化為林地，轉(zhuǎn)化面積達3 705.78 km2;1.81%的耕地轉(zhuǎn)化為人工表面，轉(zhuǎn)化面積達1 870.61 km2;0.24%的耕地轉(zhuǎn)化為濕地，轉(zhuǎn)化面積為246.90 km2。而轉(zhuǎn)出為林地和人工表面的部分分別占所有轉(zhuǎn)出耕地總面積的62.73%和31.67%，合計94.40%。

人工表面主要由耕地和林地轉(zhuǎn)化而來。10年間共有2 353.25 km2其他土地利用類型用地轉(zhuǎn)化為人工表面，其中，由耕地轉(zhuǎn)化來的人工表面面積為1 870.61 km2，由林地轉(zhuǎn)化來的人工表面面積為454.55 km2，分別占所有轉(zhuǎn)入人工表面總面積的79.49%和19.32%，合計98.81%。從空間上看，轉(zhuǎn)入的人工表面用地主要分布于重慶市和成都市(圖3)。新轉(zhuǎn)入人工表面用地占總轉(zhuǎn)入人工表面用地面積比例排名前26 位的區(qū)縣均屬于重慶市或成都市，合計達75.41%，其中，前10 位分別為重慶主城區(qū)、涪陵區(qū)、成都市主城、長壽區(qū)、雙流縣、郫縣、萬州區(qū)、豐都縣、新都區(qū)、江津區(qū)，合計達52.08%。

表2 成渝經(jīng)濟區(qū)2000—2010年土地利用轉(zhuǎn)移情況Table 2 Land use conversion matrix in research area from 2000 to 2010 in Chengdu-Chongqing Economic Zone

圖3 2000—2010年成都市、重慶市新增人工表面用地分布Fig.3 The distribution of new structures from 2000 to 2010 in Chengdu and Chongqing

(3)林地主要由耕地轉(zhuǎn)化而來。10年間共有 3 751.07 km2其他土地利用類型用地轉(zhuǎn)化為林地，由耕地轉(zhuǎn)化而來的面積則占98.79%。從空間上看，成渝經(jīng)濟區(qū)成都部分以地級市為單位、重慶部分以區(qū)縣為單位進行統(tǒng)計得出，2000—2010年新增林地占各轄區(qū)行政面積比例超過3.0%的市縣依次為榮昌縣＞眉山市＞綦江區(qū)＞自貢市＞長壽區(qū)＞永川區(qū)＞重慶主城區(qū)＞成都市＞大足區(qū)＞云陽縣，主要分布在岷江上游、沱江上游、嘉陵江上游和長江干流重慶段(圖4)。

圖4 2000—2010年成渝經(jīng)濟區(qū)新增林地分布及占行政區(qū)面積比例數(shù)量分布Fig.4 The distribution and proportion of new forest from 2000 to 2010 in Chengdu-Chongqing Economic Zone

3 土地利用變化驅(qū)動分析

3.1 偏最小二乘回歸模型構(gòu)建

土地利用變化的驅(qū)動因素主要來自于自然和社會兩方面的驅(qū)動，而短期內(nèi)區(qū)域土地利用的轉(zhuǎn)變主要是由人類活動造成的［24］。因此，在2000—2010年成渝經(jīng)濟區(qū)土地利用變化驅(qū)動力分析中，應(yīng)重點分析土地利用變化的社會驅(qū)動因素［25-27］。相關(guān)研究表明，總?cè)丝凇⒎寝r(nóng)業(yè)人口、國內(nèi)生產(chǎn)總值、固定資產(chǎn)投資額、人均國內(nèi)生產(chǎn)總值、第二產(chǎn)業(yè)比例、公路通車里程與土地利用變化密切相關(guān)［24］。由于各驅(qū)動因素之間存在不同程度的相關(guān)性，容易產(chǎn)生多重共線性而影響多元線性回歸建模的精度，而偏最小二乘回歸可以很好地解決這一問題。

由各驅(qū)動因素組成自變量數(shù)據(jù)X(表3)，由林地面積(y1)、草地面積(y2)、濕地面積(y3)、耕地面積(y4)、人工表面用地面積(y5)和其他用地面積(y6)組成因變量數(shù)據(jù)Y(表1)。

使用DPS 數(shù)據(jù)處理系統(tǒng)完成偏最小二乘回歸分析，計算結(jié)果(表4)包括模型效應(yīng)和因變量權(quán)數(shù)、模型效應(yīng)負荷量、標準回歸系數(shù)和回歸方程。根據(jù)交叉驗證法≥0.097 5，選取在自變量和因變量中提取的第一成分t1和u1構(gòu)建偏最小二乘回歸模型。

表3 自變量數(shù)據(jù)Table 3 The independent variables data

表4 模型效應(yīng)和因變量權(quán)數(shù)及模型效應(yīng)負荷量Table 4 Dependent variable weights and model effect loadings

回歸方程結(jié)果如下:

通過比較因變量權(quán)數(shù)和模型效應(yīng)負荷量，得出各自變量對因變量作用的強弱關(guān)系。模型表明，在選擇的土地利用變化各指標中，各自變量對成渝經(jīng)濟區(qū)土地利用變化的驅(qū)動強度不同，由強到弱排序為總?cè)丝?非農(nóng)人口＞國內(nèi)生產(chǎn)總值＞人均國內(nèi)生產(chǎn)總值＞固定資產(chǎn)投資額＞第二產(chǎn)業(yè)比例＞公路通車里程。可以看出，研究區(qū)內(nèi)人口壓力為土地利用變化的最主要驅(qū)動因素，國內(nèi)生產(chǎn)總值、人均國內(nèi)生產(chǎn)總值和固定資產(chǎn)投資額為次強驅(qū)動因素，再次為第二產(chǎn)業(yè)比例，而公路通車里程影響程度相對較弱。

3.2 驅(qū)動力分析

3.2.1 人口壓力驅(qū)動

“有機疏散”理論［5］認為，城鎮(zhèn)人口的快速增長將會引起城鎮(zhèn)人口和就業(yè)崗位的分布朝著可供合理發(fā)展并離開中心區(qū)域的方向發(fā)展，隨著人口的增長、結(jié)構(gòu)的變化和勞動力的聚集，必將導(dǎo)致土地利用規(guī)模的擴張和用地結(jié)構(gòu)的改變。人口增長的越快，人地矛盾就越顯突出，對人工表面用地需求也相應(yīng)增加，而土地總面積的有限性又決定必須通過土地系統(tǒng)中各類型土地互相轉(zhuǎn)化來滿足人類生產(chǎn)和生活的需求。

從圖5 和表5 可以看出，2000—2010年成渝經(jīng)濟區(qū)總?cè)丝谠黾?59.13 萬人，增幅為6.83%;非農(nóng)業(yè)人口增加1 152.25 萬人，增幅達到59.89%。其中，成都市總?cè)丝谠龇?3.40%，非農(nóng)業(yè)人口增幅高達88.18%;重慶市總?cè)丝谠龇鶠?.73%，非農(nóng)業(yè)人口增幅達71.58%。從農(nóng)業(yè)和非農(nóng)業(yè)人口比例來看，成渝經(jīng)濟區(qū)非農(nóng)業(yè)人口比例增加9.9%，重慶市和成都市非農(nóng)業(yè)人口比例分別增加13.3% 和22.5%，導(dǎo)致重慶和成都人工表面用地規(guī)模的快速擴增，特別是成都市由2010年的農(nóng)業(yè)人口占主導(dǎo)發(fā)展為2010年非農(nóng)業(yè)人口占主導(dǎo)，說明成渝經(jīng)濟區(qū)的人口結(jié)構(gòu)正朝著非農(nóng)業(yè)人口方向演變，城鎮(zhèn)化進程十分明顯。

圖5 2000—2010年成渝經(jīng)濟區(qū)新增非農(nóng)人口占總?cè)丝诒壤植糉ig.5 The distribution and proportion of new urban population from 2000 to 2010 in Chengdu-Chongqing Economic Zone

表5 2000—2010年重慶市、成都市人口變化情況Table 5 Population change from 2000 to 2010 in Chongqing and Chengdu

3.2.2 經(jīng)濟發(fā)展驅(qū)動

成渝經(jīng)濟區(qū)國內(nèi)生產(chǎn)總值由2000年的4 943.548億元增至2010年的22 220.228 億元，增幅高達349.48%;人均國內(nèi)生產(chǎn)總值由2000年的5 125 元 增 至 2010年 的 24 311 元，增 幅 為374.36%;全社會固定資產(chǎn)投資額由2000年的1 356.145億元增至2010年的17 424.482 億元，增幅高達1 184.85%;第二產(chǎn)業(yè)比例由2000年的43.98%增至50.31%。國內(nèi)生產(chǎn)總值和固定資產(chǎn)投資額的大幅度上漲，表明成渝經(jīng)濟區(qū)經(jīng)濟的迅猛發(fā)展，固定資產(chǎn)投資額增加的直接表現(xiàn)是人工表面用地的增加，尤其是基礎(chǔ)設(shè)施建設(shè)的增加和房地產(chǎn)行業(yè)的迅速發(fā)展，耕地則成為人工表面用地的最主要來源，耕地面積的減少引發(fā)農(nóng)業(yè)勞動力向城市轉(zhuǎn)移，從而加速成渝經(jīng)濟區(qū)城鎮(zhèn)化的進程。

2000—2010年，成渝經(jīng)濟區(qū)三產(chǎn)比例從18.39 ∶43.98 ∶37.63 轉(zhuǎn)變?yōu)?1.65 ∶50.31 ∶38.05，呈現(xiàn)出二產(chǎn)占主導(dǎo)的“二三一”型產(chǎn)業(yè)結(jié)構(gòu)模式，且第一產(chǎn)業(yè)比例下降、第二產(chǎn)業(yè)比例上升趨勢明顯，反映了工業(yè)和建筑業(yè)的迅速發(fā)展趨勢，而第二產(chǎn)業(yè)主要構(gòu)成為工業(yè)，因此二產(chǎn)比例的快速上升實際反映出成渝經(jīng)濟區(qū)工業(yè)化快速發(fā)展的態(tài)勢，產(chǎn)業(yè)工業(yè)園區(qū)的入駐占用大量耕地及林草地，“雙核兩帶”的產(chǎn)業(yè)布局促進了區(qū)內(nèi)土地利用系統(tǒng)快速轉(zhuǎn)化，導(dǎo)致成都和重慶兩大都市區(qū)、成德綿城市經(jīng)濟帶和沿長江城市經(jīng)濟帶人工表面用地的迅速擴張。

3.3.3 國家政策驅(qū)動

國家和地方的政策制度因素對區(qū)域土地利用方式及其變化起強制性作用［28-29］，土地利用結(jié)構(gòu)及其變化是特定時期經(jīng)濟體制及政策水平下的產(chǎn)物，各級政府經(jīng)濟政策、土地政策、環(huán)保政策、農(nóng)業(yè)政策、水利政策等宏觀調(diào)控將直接影響區(qū)域土地利用的狀態(tài)。2000年，成渝經(jīng)濟區(qū)被確定為西部大開發(fā)的三大重點區(qū)域之一，2006年，成都和重慶被納入全國統(tǒng)籌城鄉(xiāng)綜合配套改革試驗區(qū)，一系列政府宏觀政策支持推進了成渝經(jīng)濟區(qū)的快速發(fā)展和城鎮(zhèn)化進程，加速了居住地、工業(yè)用地、交通用地等人工表面的形成，原有城市不斷擴張，新興工業(yè)城市不斷興起，裝備制造業(yè)、電子信息、生物醫(yī)藥、汽摩、清潔能源、食品飲料業(yè)、國防科技工業(yè)等優(yōu)勢產(chǎn)業(yè)迅速發(fā)展，不同土地利用類型的轉(zhuǎn)化更加頻繁，逐漸形成以成都和重慶兩大城市為核心的土地利用空間格局。

1999年起，成渝經(jīng)濟區(qū)內(nèi)各市縣陸續(xù)啟動退耕還林工程，2000年以三峽庫區(qū)為界的長江上游地區(qū)開始實施國務(wù)院批準的天然林資源保護工程［3］，受西部開發(fā)“生態(tài)退耕”政策的影響，研究區(qū)內(nèi)有3 705.8 km2的耕地轉(zhuǎn)化為林地，退耕還林還草政策成效十分明顯。2000—2010年，長江重慶段下游開縣、云陽縣、萬州區(qū)、忠縣、墊江縣、長壽區(qū)、豐都縣等區(qū)縣濕地面積共增加322.19 km2，占整個成渝經(jīng)濟區(qū)10年來濕地增加總面積的53.21%，主要原因則為三峽大壩的竣工，三峽水位的抬高使庫區(qū)林地、耕地等土地利用類型轉(zhuǎn)化為濕地。

4 結(jié)論

(1)2000—2010年，成渝經(jīng)濟區(qū)土地利用格局發(fā)生了明顯變化，主要體現(xiàn)在:耕地面積減少，林地、草地、濕地、人工表面和其他用地面積增加，其中，人工表面增加面積最大且增長速度最快。從土地利用轉(zhuǎn)移變化看，耕地主要轉(zhuǎn)換為林地和人工表面，人工表面的增加主要來自于耕地、林地和草地。

(2)偏最小二乘回歸模型表明，成渝經(jīng)濟區(qū)人口壓力驅(qū)動，特別是非農(nóng)業(yè)人口的增加是研究區(qū)土地利用變化的最主要驅(qū)動因素;國內(nèi)生產(chǎn)總值、人均國內(nèi)生產(chǎn)總值和固定資產(chǎn)投資額等經(jīng)濟指標是影響土地利用變化的次強驅(qū)動因素。另一方面，國家宏觀經(jīng)濟政策的調(diào)控帶動區(qū)內(nèi)經(jīng)濟和產(chǎn)業(yè)的快速發(fā)展，促使著人工表面用地、耕地、林地、濕地等土地利用類型之間的快速轉(zhuǎn)化，是成渝經(jīng)濟區(qū)土地利用變化的又一重要驅(qū)動因素。

(3)通過RS 與GIS 結(jié)合的技術(shù)手段，研究成渝經(jīng)濟區(qū)土地利用變化特征及其驅(qū)動力，有利于解決快速城市化地區(qū)土地利用與社會發(fā)展之間的矛盾，協(xié)調(diào)區(qū)域人—地關(guān)系的協(xié)調(diào)發(fā)展，為人地系統(tǒng)的合理調(diào)控提供數(shù)據(jù)、資料支持，對區(qū)域國土資源安全與可持續(xù)發(fā)展具有重要的科學(xué)意義。土地利用強度和結(jié)構(gòu)的巨大變化會引起一系列區(qū)域性生態(tài)環(huán)境問題，成渝經(jīng)濟區(qū)土地利用變化對生態(tài)環(huán)境變化的影響將是今后課題研究的重要方向。

［1］黃金良，李青生，洪華生，等. 九龍江流域土地利用/景觀格局:水質(zhì)的初步關(guān)聯(lián)分析［J］.環(huán)境科學(xué)，2011，32(1):64-72.

［2］劉紀遠，劉明亮，莊大方，等. 中國近期土地利用變化的空間格局分析［J］.中國科學(xué):D 輯，2002，32(12):1031-1040.

［3］劉紀遠，張增祥，俆新良，等.21 世紀初中國土地利用變化的空間格局與驅(qū)動分析［J］. 地理學(xué)報，2009，64 (12):1411-1420.

［4］TURNER B L，LAMBIN E F，REENBERG A.The emergence of land change science for global environmental change and sustainability［J］.PNAS，2007，104(52):20666-20671.

［5］雷征，董捷.武漢城市圈土地利用結(jié)構(gòu)演變規(guī)律及驅(qū)動因素分析［J］.農(nóng)業(yè)現(xiàn)代化研究，2010，31(2):147-151.

［6］TURNER B L，MOSS R H，SKOLE D L，et al. Relation of land use and global land cover change:IGBP Report No.24 and HDP Report No. 5［R］. Stockholm:Royal Swedish Academy of Sciences，2003.

［7］TURNER B L，MEYER W B.Land use and land cover in global environmental change:considerations for study［J］. Journal of International Social Science，1991，130:69-79.

［8］SKOLE D，TUCKER C. Tropical deforestation and habitat fragmentation in the Amazon:satellite data from 1978-1988［J］.Science，1993，260:1905-1910.

［9］OJIMA D S，GALVIN K A，TURNER B L Ⅱ.The global impact of land-use change［J］.BioScience，1994，44(5):300-304.

［10］FISCHER G，SUN L X. Model based analysis of future land-use development in China ［J］. Agriculture， Ecosystems ＆Environment，2001，85(1/2/3):163-176.

［11］OTSUBO K.Land use for global environmental conservation (LU/GEC):final report of the LU/GEC first phase:1995-1997［R］.Tokyo: National Institute for Environmental Studies，Environmental Agency of Japan，1999.

［12］莊大方，劉紀遠. 中國土地利用程度的區(qū)域分異模型研究［J］.自然資源學(xué)報，1997，12(2):105-111.

［13］張明.榆林地區(qū)土地利用與土地覆蓋變化研究［D］. 北京:中國科學(xué)院，1998.

［14］國務(wù)院辦公廳.全國主體功能區(qū)規(guī)劃［R］.北京:國務(wù)院辦公廳，2010.

［15］袁興中，肖紅艷，顏文濤，等. 成渝經(jīng)濟區(qū)土地利用與生態(tài)服務(wù)價值動態(tài)分析［J］.生態(tài)學(xué)雜志，2012，31(1):180-186.

［16］邱磊，廖和平，龍雨嬌，等. 成渝經(jīng)濟區(qū)的城市土地集約利用評價及時空特征分析［J］. 西南師范大學(xué)學(xué)報:自然科學(xué)版，2010，35(1):208-213.

［17］王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討［J］. 地理科學(xué)進展，1999，18(1):81-87.

［18］朱會義，李秀彬.關(guān)于區(qū)域土地利用變化指數(shù)模型方法的討論［J］.地理學(xué)報，2003，58(5):643-650.

［19］高吉喜，呂世海，劉軍會，等. 中國生態(tài)交錯帶［M］. 北京:中國環(huán)境科學(xué)出版社，2009:260-263.

［20］任斐鵬，江源，熊興，等.東江流域近20年土地利用變化的時空差異特征分析［J］.資源科學(xué)，2011，33(1):143-152.

［21］張旸，周成虎，戴錦芳，等. 偏最小二乘回歸方法在土地利用結(jié)構(gòu)研究中的應(yīng)用［J］. 自然資源學(xué)報，2004，19(6):804-810.

［22］劉賢趙，張安定，李嘉竹.地理學(xué)數(shù)學(xué)方法［M］.北京:科學(xué)出版社，2009:69-75.

［23］王惠文.最小二乘回歸方法及其應(yīng)用［M］.北京:國防工業(yè)出版社，1999:150-157.

［24］張亞利，游揚聲，劉星.重慶市沙坪壩區(qū)土地利用變化的偏最小二乘驅(qū)動機制分析［J］. 地理與地理信息科學(xué)，2013，29(6):89-94.

［25］譚少華，倪紹祥.區(qū)域土地利用變化驅(qū)動力的成因分析［J］.地理與地理信息科學(xué)，2005，21(3):47-50.

［26］陳美球，劉序，陳文波，等.鄱陽湖地區(qū)1985—2000年土地利用格局變化及其社會經(jīng)濟驅(qū)動力研究:Ⅱ. 土地利用格局變化驅(qū)動力分析［J］.安徽農(nóng)業(yè)大學(xué)學(xué)報，2006，33(1):13-129.

［27］吳明發(fā)，歐名豪，廖榮浩，等. 經(jīng)濟發(fā)達地區(qū)土地利用變化及其驅(qū)動力分析:以廣東省為例［J］. 水土保持研究，2012，19(1):179-183.

［29］婁和震，楊勝天，周秋文.延河流域2000—2010年土地利用/覆蓋變化及驅(qū)動力分析［J］. 干旱區(qū)資源與環(huán)境，2014，28(4):15-21. ?