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

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

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

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

1.1 研究區(qū)概況

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

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

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

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

1.3 分析方法

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

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

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

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

式中:Ki為研究時間范圍內(nèi)第i 種單一土地利用類型的動態(tài)度;Ua為研究初期a 時期第i 種土地利用類型的數(shù)量;Ub為研究末期b 時期第i 種土地利用類型的數(shù)量;T 為研究時段，當(dāng)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］。其表達(dá)式為:

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

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

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

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

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

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

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

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

2 土地利用變化特征

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

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

(1)2000—2010年，林地面積增加2 719.15 km2，所占比例增加1.3 個百分點(diǎn);草地面積增加69.44 km2，所占比例增加0.03 個百分點(diǎn);濕地面積增加594.08 km2，所占比例增加0.28 個百分點(diǎn);耕地面積減少明顯，減少面積達(dá)5 840.27 km2，所占比例下降2.8 個百分點(diǎn);人工表面用地面積迅速擴(kuò)大，增加了2 337.03 km2，所占比例增加1.12 個百分點(diǎn);其他用地面積增大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)濟(jì)區(qū)2005—2010年土地利用變化速度大于2000—2005年，呈現(xiàn)加快變化趨勢。

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

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

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

耕地主要轉(zhuǎn)化為人工表面用地和林地。10年間共有5 907.15 km2耕地轉(zhuǎn)出為其他幾種土地利用類型，其中，3.58%的耕地轉(zhuǎn)化為林地，轉(zhuǎn)化面積達(dá)3 705.78 km2;1.81%的耕地轉(zhuǎn)化為人工表面，轉(zhuǎn)化面積達(dá)1 870.61 km2;0.24%的耕地轉(zhuǎn)化為濕地，轉(zhuǎn)化面積為246.90 km2。而轉(zhuǎn)出為林地和人工表面的部分分別占所有轉(zhuǎn)出耕地總面積的62.73%和31.67%，合計(jì)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%，合計(jì)98.81%。從空間上看，轉(zhuǎn)入的人工表面用地主要分布于重慶市和成都市(圖3)。新轉(zhuǎn)入人工表面用地占總轉(zhuǎn)入人工表面用地面積比例排名前26 位的區(qū)縣均屬于重慶市或成都市，合計(jì)達(dá)75.41%，其中，前10 位分別為重慶主城區(qū)、涪陵區(qū)、成都市主城、長壽區(qū)、雙流縣、郫縣、萬州區(qū)、豐都縣、新都區(qū)、江津區(qū)，合計(jì)達(dá)52.08%。

表2 成渝經(jīng)濟(jì)區(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)濟(jì)區(qū)成都部分以地級市為單位、重慶部分以區(qū)縣為單位進(jìn)行統(tǒng)計(jì)得出，2000—2010年新增林地占各轄區(qū)行政面積比例超過3.0%的市縣依次為榮昌縣＞眉山市＞綦江區(qū)＞自貢市＞長壽區(qū)＞永川區(qū)＞重慶主城區(qū)＞成都市＞大足區(qū)＞云陽縣，主要分布在岷江上游、沱江上游、嘉陵江上游和長江干流重慶段(圖4)。

圖4 2000—2010年成渝經(jīng)濟(jì)區(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)濟(jì)區(qū)土地利用變化驅(qū)動力分析中，應(yīng)重點(diǎn)分析土地利用變化的社會驅(qū)動因素［25-27］。相關(guān)研究表明，總?cè)丝?、非農(nóng)業(yè)人口、國內(nèi)生產(chǎn)總值、固定資產(chǎn)投資額、人均國內(nèi)生產(chǎn)總值、第二產(chǎn)業(yè)比例、公路通車?yán)锍膛c土地利用變化密切相關(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)完成偏最小二乘回歸分析，計(jì)算結(jié)果(表4)包括模型效應(yīng)和因變量權(quán)數(shù)、模型效應(yīng)負(fù)荷量、標(biāo)準(zhǔn)回歸系數(shù)和回歸方程。根據(jù)交叉驗(yàn)證法≥0.097 5，選取在自變量和因變量中提取的第一成分t1和u1構(gòu)建偏最小二乘回歸模型。

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

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

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

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

3.2 驅(qū)動力分析

3.2.1 人口壓力驅(qū)動

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

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

圖5 2000—2010年成渝經(jīng)濟(jì)區(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)濟(jì)發(fā)展驅(qū)動

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

2000—2010年，成渝經(jīng)濟(jì)區(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)比例的快速上升實(shí)際反映出成渝經(jīng)濟(jì)區(qū)工業(yè)化快速發(fā)展的態(tài)勢，產(chǎn)業(yè)工業(yè)園區(qū)的入駐占用大量耕地及林草地，“雙核兩帶”的產(chǎn)業(yè)布局促進(jìn)了區(qū)內(nèi)土地利用系統(tǒng)快速轉(zhuǎn)化，導(dǎo)致成都和重慶兩大都市區(qū)、成德綿城市經(jīng)濟(jì)帶和沿長江城市經(jīng)濟(jì)帶人工表面用地的迅速擴(kuò)張。

3.3.3 國家政策驅(qū)動

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

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

4 結(jié)論

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

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

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

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

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

［3］劉紀(jì)遠(yuǎn)，張?jiān)鱿?，俆新良，等.21 世紀(jì)初中國土地利用變化的空間格局與驅(qū)動分析［J］. 地理學(xué)報(bào)，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］莊大方，劉紀(jì)遠(yuǎn). 中國土地利用程度的區(qū)域分異模型研究［J］.自然資源學(xué)報(bào)，1997，12(2):105-111.

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

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

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

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

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

［18］朱會義，李秀彬.關(guān)于區(qū)域土地利用變化指數(shù)模型方法的討論［J］.地理學(xué)報(bào)，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é)報(bào)，2004，19(6):804-810.

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

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

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

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

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

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

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