淮北平原農業區耕地面積時空格局演變與驅動力分析

鐘博星　何太蓉　張喆　鄔玉琴

摘要：【目的】探索安徽淮北平原農業區耕地面積的時空格局演變及其驅動力，為保障區域耕地資源可持續利用與管理提供科學依據。【方法】根據2000～2013年相關統計資料，采用耕地動態度、耕地指數及耕地重心模型對淮北平原農業區的耕地面積時空格局演變進行研究，并運用主成分分析法對其耕地面積時空格局演變的驅動力進行分析。【結果】（1）2000～2013年淮北平原農業區總耕地面積及其大部分縣（區）耕地面積小幅減少，其余部分市轄區耕地面積有一定增加；（2）淮北平原農業區的耕地在空間上呈西多東少的分布格局，2000～2013年耕地重心總體上從西北向東南方向轉移，但遷移幅度不明顯，并未影響研究區西多東少的分布格局；（3）主成分分析結果表明淮北平原農業區耕地面積時空格局演變的主要驅動力為社會經濟與城鎮化、農業生產能力和災變。【建議】加強現有耕地資源的保護力度，合理開發后備土地資源；依靠科技提升農業生產能力，完善土地管理政策，保障農戶切身利益；控制人口數量，樹立耕地資源危機意識。

關鍵詞： 耕地面積；時空格局；驅動力；主成分分析；淮北平原農業區

中圖分類號： F323.211 文獻標志碼：A 文章編號：2095-1191（2017）01-0174-07

Abstract：【Objective】The present study explored spatial-temporal pattern evolution and its driving forces of agricultural areas in Huaibei plain， in order to provide scientific basis for sustainable utilization and management of cultivated land resources. 【Method】Based on the relevant statistical data from 2000 to 2013， this paper analyzed spatial-temporal pattern evolution of cultivated land acreage of agricultural areas in Huaibei plain， Anhui province. The main indexes included dynamic degree of cultivated land， cultivate land index and gravity center model of cultivated land. Then the driving forces in the spatial-temporal pattern evolution of cultivated land area were studied by principal component analysis.【Result】（1）Total cultivated land area in Huaibei plain and most counties（district） decreased to a small extent during 2000-2013， but cultivated land area in the rest parts increased. （2）The space distribution pattern of cultivated land area in the west was more than that in the east. The gravity center of cultivated land shifted from northwest to southeast during study period. But for the whole study area， the migration range was not large. And it did not affect the whole distribution pattern.（3）Principal component analysis showed that main driving forces were social economy and urbanization， capability of agricultural production and natural disasters. 【Suggestion】Following suggestions are proposed： strengthening the protection for existing cultivated land resources， rational development of reserve land resources； relying on science and technology to improve agricultural production capacity， improving land management policies to ensure the vital interests of farmers； controlling population and establishing the consciousness of cultivated land resources crisis.

Key words： cultivated land area； spatial-temporal pattern evolution； driving force； principal component analysis； agricultural areas in Huaibei plain

0 引言

【研究意義】耕地是人類賴以生存的物質基礎，是國家糧食安全和社會穩定的重要支撐和保障（王文剛等，2012）。自改革開放以來，我國處于城鎮化和工業化快速發展階段，國家經濟建設占用了一定數量的耕地（陳紅等，2010），且由于人口數量大，我國人均耕地面積遠低于世界平均水平。因此，分析耕地的時空格局演變及其驅動力，對合理利用耕地資源，實現區域耕地動態平衡和耕地資源的可持續利用（曹銀貴等，2007），以及協調好耕地保護與經濟發展、生態環境間的關系有著重要意義（梁海鷗，2012）。【前人研究進展】國外學者對耕地面積時空變化的研究主要依托對土地利用土地覆被變化（LUCC）的研究，如Barbier等（2010）以越南森林轉型為例，指出影響土地利用變化的驅動因素主要包括內生社會生態系統和外生社會經濟系統兩個方面，內生社會生態系統是穩定森林覆蓋率的主導因素，外生社會經濟系統是土地轉型的根本動力；Lambin等（2010）探討了土地利用與森林覆被轉型的關系，認為二者是不可分割的總體模式；Ferreira等（2015）基于航空相片和農業人口普查數據探索巴西大西洋甘蔗用地的動態變化，認為強制性政策因素在所有影響土地變化的因素中最顯著。國內學者則運用不同方法對一些地區進行了相應研究，如鄧楚雄等（2013）結合GIS分析了長株潭城市群地區耕地數量的時空差異，并對研究區耕地數量變化的驅動因子進行定量分析；王國剛等（2014）基于耕地集約利用模型分析了1996～2008年中國省域的耕地集約利用時空變化及其影響因素；楊春艷等（2015）以四期遙感資料為數據源，對1980～2010年西藏耕地面積時空變化特征進行研究。總之，這些相關研究分別從全國（王國剛等，2014；張榮天和焦華富，2015）或區域（鄧楚雄等，2013；楊春艷等，2015）等不同時空尺度對耕地的時空格局演變及其驅動力進行了有益探討，為今后的耕地保護提供了參考依據。【本研究切入點】淮北平原農業區是安徽省最大的糧、棉、油產區，也是我國重要的農業商品糧基地和農業開發治理的重點區域之一，對保障地區糧食安全和社會經濟發展起著關鍵作用（洪舒蔓等，2014）。近年來，該區社會經濟快速發展，人口激增，對土地的需求日益增多，人地矛盾突出，但目前有關該區域耕地面積時空演變及其驅動因素的研究鮮見報道。【擬解決的關鍵問題】以淮北平原農業區各縣（區）為研究樣本區域，揭示其2000～ 2013年耕地面積的時空格局演變過程及其特征，并對該區域耕地面積動態變化的驅動力加以研究，為其耕地資源可持續利用及社會經濟可持續發展提供科學依據。

1 研究區概況與研究方法

1. 1 研究區概況

淮北平原農業區位于安徽省淮河干流以北至沙潁河以南的廣大地區，東、北部與蘇、魯接壤，西、西北部與豫毗鄰，南臨淮河。該區地處暖溫帶的南緣，年平均氣溫14～15 ℃，晝夜溫差大，日照充足，雨熱同季，無霜期200～220 d，地下水資源豐富，全區以砂姜黑土和潮土為主，屬沖積、洪積平原，地勢平坦，微地形起伏，非常適于農業的綜合開發與利用。該區總面積約4×104 km2；2013年末總常住人口3152.48萬人，其行政區有淮北、亳州、宿州、蚌埠、阜陽和淮南6市，為研究方便，每個地級市所轄區合并為一個單元，即6個市轄區（圖1）。

1. 2 數據來源及研究方法

本研究數據來自《安徽省統計年鑒（2001～2014）》和研究區各市的統計信息網。選用的指標包括耕地動態度、耕地指數和耕地重心模型。運用Arcgis 9.3呈現耕地的空間差異及耕地重心的空間移動軌跡，應用SPSS 17.0進行主成分分析。

1. 2. 1 耕地動態度 耕地動態度用于表達研究區一定時段內耕地的數量變化情況（王秀蘭和包玉海，1999）。其表達式為：

K=××100 （1）

式中，K為耕地動態度，LA（i， t1）和LA（i， t2）分別為研究期初和期末時耕地的數量，T為時間。K的絕對值越大，說明變化幅度越大；當K>0，表明耕地數量增加，反之，耕地數量減少。

1. 2. 2 耕地指數 耕地指數是一種反映耕地利用空間格局的有效方法。其表達式為：

Ri（%）=（Fi /Si）×100（2）

式中，Ri為耕地指數，Fi為第i縣（區）的耕地面積，Si為第i縣（區）的行政區面積。

1. 2. 3 耕地重心模型 重心模型能夠可視化地表達出某要素的空間集聚或偏移程度。重心模型的坐標一般用經緯度表示，原始的重心坐標通常選取省、市（區、縣）政府所在地的經緯度。本研究使用的耕地重心模型表達式為：

Xt=（Cit×Xi）/Cit（3）

Yt=（Cit×Yi）/Cit（4）

式中，Xt和Yt分別代表t年耕地重心坐標，Cti是第i縣（區）的耕地面積，Xi和Yi是第i縣（區）政府所在地的經緯坐標。

2 淮北平原農業區耕地面積時空格局演變

2. 1 耕地面積時間格局演變

2. 1. 1 耕地總面積時間格局演變 由圖2可看出，2000～2013年淮北平原農業區的耕地總面積呈階段性遞減趨勢，耕地總面積從2000年的210.98萬ha減少到2013年的209.65萬ha，減少了13296.00 ha。淮北平原農業區耕地總面積變化大致可分為3個階段。第一階段：2000～2003年，此階段耕地總面積呈急劇減少態勢，由2000年的210.98萬ha減少到2003年的207.79萬ha，凈減少了31953.00 ha，其中2002年末耕地總面積變化幅度呈近似垂直減少趨勢，其原因是2000年之后隨著皖北經濟的快速發展，基建項目增多，造成耕地面積大幅減少。第二階段：2003～2008年，此階段耕地總面積呈明顯的波動狀態，一是2004年安徽省行政區劃調整把原本不在該行政區的屬于合肥市長豐縣的孔店鄉、三和鄉、曹庵鎮、史院鄉、孤堆鄉、楊公鎮和孫廟鄉等7個鄉（鎮）劃入淮南市管轄，因此2004年該行政區耕地面積顯示增加；二是相對于2004年，2005年耕地面積減少較多的原因是建設用地的大增；三是2005～2007年耕地面積遞增，其原因可能是此間政府推出土地復墾政策和措施的影響；四是耕地面積在2008年又迅速減少，其原因是安徽省實施“861”工程，興建交通水利設施占用了一定量耕地。第三階段：2008～2013年，此階段為相對穩定階段，每年維持約209.00萬ha的耕地保有量，主要原因是生態保護愈來愈被普遍關注，生態退耕政策得以有效貫徹實施。

2. 1. 2 各縣（區）耕地面積時間格局演變 用耕地動態度分析淮北平原農業區各縣（區）在2000～2013年耕地變化情況。根據公式（1）得到研究區23個縣（區）14年間的耕地動態度。由圖3可知，2000年以來，淮北平原農業區各縣（區）耕地面積及其變化速度存在較大差異。整個淮北平原農業區總耕地動態度為-0.05%，大部分縣（區）與其相似，說明整個農業區及其大部分縣（區）的耕地面積小幅減少；碭山、阜南、濉溪、固鎮、五河和懷遠等6縣的耕地動態度負值較大，說明其耕地面積減少幅度較大；蚌埠、淮北和淮南3個市轄區為偏大正值，說明其耕地面積增加，其原因可能是行政區劃調整，尤其是周邊農村調整到市轄區。

2. 2 耕地面積空間格局演變

2. 2. 1 耕地面積的空間分布格局 本研究采用耕地指數來表達耕地的空間分布格局。根據公式（2）計算出淮北平原農業區各縣（區）的耕地指數；運用Arcgis的自然斷裂法將該區域耕地空間分布格局劃分為3個部分，即稀疏區Ri∈[0，0.5]，一般區Ri∈（0.5，0.6]，密集區Ri∈（0.6，0.7]；最后生成淮北平原農業區2000和2013年的耕地面積空間格局（圖4）。由圖4可知，淮北平原農業區的耕地密集區主要集中在西部，其中太和、臨泉及渦陽3縣分布著可觀的耕地資源，而東部的耕地資源顯然不如西部，除濉溪、靈璧、固鎮和懷遠4縣分布有足量的耕地資源外，其余的淮北、宿州、蚌埠和淮南4個市轄區及鳳臺、五河、泗縣、蕭縣和碭山5縣均為耕地稀疏區。因此，淮北平原農業區耕地面積總體上呈西多東少的空間格局。

2. 2. 2 耕地重心變化 耕地重心代表了區域耕地分布在空間上的集中性，能夠較好地反映出耕地的空間形態變化（王曉文等，2014）。根據公式（3）和（4）可計算得到淮北平原農業區2000～2013年各年的耕地坐標。由圖5可看出，2000～2013年淮北平原農業區的耕地重心總體上向東南移動，從2000年起，耕地重心經向移動0.0082°E，即0.913 km，緯向移動0.0144°N，即1.603 km，14年間向東南方向總共移動1.778 km，但各年份耕地重心移動距離長短不一，移動方向不盡相同。其中，2003～2004年耕地重心遷移幅度最明顯，向東南移動1.425 km，其主要原因為行政區劃的調整增加了研究區的耕地面積；2004～2007年耕地重心出現回調現象，總體向西北方向遷移，即經向從116.5452°E到116.5419°E，回調了0.0033°，其原因主要是該時期研究區耕地總面積處于明顯波動的階段，導致耕地重心出現回調現象。值得注意的是，研究區在14年期間的耕地重心雖有一定幅度的遷移，但總體上遷移幅度較小，對耕地面積西多東少的分布格局影響不明顯，且耕地重心的移動現象與耕地總面積時間格局的演變階段基本吻合。

3 淮北平原農業區耕地面積時空格局演變驅動力分析

3. 1 指標選取與分析方法

影響耕地面積時空格局演變的因子較多，本研究根據主成分分析（徐秋艷，2011）的基本要求和淮北平原農業區各縣（區）現有資料，依據科學性、系統性、代表性和可獲取性等原則，選取了12個指標，分別為X1地區生產總值（億元）、X2城鎮化率（%）、X3房屋竣工建筑面積（ha）、X4公路里程（km）、X5受災面積（千ha）、X6復種指數（%）、X7農業化肥施用強度（kg/ha）、X8固定資產投資（億元）、X9三次產業比重（%）、X10農業總產值（億元）、X11農業機械總動力（萬kW）、X12農村用電量（萬kW），其中，城鎮化率采用人口指標計算。

利用主成分分析法對2000～2013年淮北平原農業區耕地面積時空格局演變的驅動力進行研究。首先得出前3個主成分因子的累計貢獻率，其累計貢獻率達98.065；然后采用方差最大旋轉法得出3個主成分因子載荷矩陣，將X1、X2、X3、X4、X8和X9歸為第一主成分，X6、X7、X11和X12歸為第二主成分，X5歸為第三主成分，同時，亦將這三大類分別納入社會經濟與城鎮化因素、農業生產能力和災變因素。

3. 2 驅動力分析

3. 2. 1 社會經濟與城鎮化因素 經濟發展是土地利用結構及其演變最根本的動力（陳朝，2007）。淮北平原農業區因其突出的農業地位，決定著農業生產是其社會經濟發展的前提和基礎，也是影響其耕地面積時空變化的決定性因素。社會經濟對耕地面積時空變化的驅動主要表現在農業結構調整和經濟結構優化兩個方面（表1，圖6）。由表1可知，淮北平原農業區種植業占大農業的比重從2000年的63.3%下降到2013年的56.6%，而同期牧業產值比重從24.7%增長到30.7%，農林牧漁服務業產值比重也略有上升，說明該區域農產品市場正趨于多樣化，以種植業為主的小農業比重逐漸減小。

優化經濟結構是淮北平原農業區實現農業現代化的有效途徑。在農業產值得以保障的前提下，各縣（區）大力發展第二、三產業，實現地區經濟的快速發展。如圖6所示，淮北平原農業區的地區總產值從2000年的923.24億元增長到2013年的5398.81億元，14年來增長了近5倍。由圖6可看出，第一產業在國民經濟中的比重逐年減少，第三產業發展較平緩，而第二產業發展較快，由2000年的34%增長到2013年的51%。李艷（2010）研究表明，此階段工礦用地規模和塌陷區面積快速增加，致使耕地面積大量流失；2000～2013年該地區固定資產投資額占地區生產總值比率始終小于全國，說明該區資金短缺、政府招商引資力度不足，投資對耕地面積時空演變的驅動作用較弱，與其他大城市有所不同（郇紅艷和孫君，2012）。

城鎮化是耕地面積時空演變的關鍵驅動因素。2000～2013年該地區房屋竣工建筑面積從612.97 ha增長到2406.53 ha，增長了近3倍；公路里程由2000年的24247 km增加到2013年53251 km，增加了29004 km。城鎮化加速城市擴張、舊城改造、公路交通設施建設等，推動城市房地產業的快速發展，同時大量農村勞動力脫離農耕生產，部分農村人口在城市安家落戶導致農村耕地閑置，而農村宅基地尚未轉化為農用地，出現城市農村“兩頭”占的局面，必然帶來耕地的占用問題。

3. 2. 2 農業生產能力因素 復種指數、農業化肥施用強度、農業機械總動力和農村用電量是構成農業生產能力的主要因素。復種指數能有效反映耕地的利用程度和集約化水平。據統計，2000～2013年淮北平原農業區復種指數從192%提高到210%，說明耕作方式不斷向集約化方向轉變；農業化肥施用強度、農業機械總動力和農村用電量分別由555 kg/ha、1514.86和172014萬kW增加到758 kg/ha、3265.90和486848萬kW，分別增長37%、116%和183%；糧食總產量由1234.29萬t增加到1870.32萬t，增加了636.03萬t。因此，農業化肥施用強度加大、農業機械總動力提高、農村用電量增加均在一定程度上提高糧食單產，實現以較少的耕地面積確保等量甚至增量的糧食總產量，確保了糧食安全（劉彥彤等，2011）。這也是該區耕地非農化的前提。

3. 2. 3 災變因素 淮北平原農業區地處我國秦淮線上，特殊的地理位置使其氣候具有明顯不穩定性，是我國氣候變化的敏感區之一。夏季旱澇災害頻繁，冬季低溫、寒潮災害多發，給當地的農業生產帶來嚴重挑戰。同時，因全球變暖病蟲害增多，農藥施用量增加，導致農業面源污染不斷惡化，災變與經濟因素間相互作用對耕地數量產生影響，使農民對農業生產的積極性大大降低，減少對耕地的耕作需求，出現耕地撂荒棄耕，可能引起耕地面積的減少。

4 建議

（1）加強現有耕地資源的保護力度，合理開發后備土地資源（趙永華等，2103）。農業生產、保障糧食安全是淮北平原農業區的工作重心，提升農業生產能力、提高農業產量是耕地非農化的前提，因此，建立健全耕地保護機制，依法管理耕地，減少耕地保護與未來建設的沖突。（2）依靠科技，提升農業生產能力；完善政策，保證農戶切身利益。一是要加強科學技術應用，建立“3S”耕地資源動態監測和災害預警體系；二是不斷增加農業科技投入，開發先進的農業生產技術和設備，研究并推廣耕地災害防治技術；三是完善耕地保護政策，依法管理土地，保障農產品價格，提高農民收入。（3）控制人口數量，樹立耕地資源危機意識。人口快速增長帶來耕地壓力的增大，淮北平原農業區除煤炭資源及部分非金屬資源外，大多數資源人均占有量相對較少，因此，亟需增強淮北平原農業區耕地保護意識。

參考文獻：

曹銀貴，王靜，程燁，劉愛霞，劉慶，吳劍. 2007. 三峽庫區開縣30年耕地變化分析[J]. 中國土地科學，21（1）：43-47.

Cao Y G，Wang J，Cheng Y，Liu A X，Liu Q，Wu J. 2007. Change of cultivated land area in Kaixian county in the past 30 years[J]. China Land Science，21（1）：43-47.

陳朝. 2007. 湖北耕地利用變化與驅動力分析[D]. 武漢：華中師范大學.

Chen C. 2007. An analysis on the changes and their driving forces of cultivated land use in Hubei[D]. Wuhan：Central China Normal University.

陳紅，吳世新，馮雪力. 2010. 新疆耕地時空變化特征[J]. 地理科學進展，29（3）：312-318.

Chen H，Wu S X，Feng X L. 2010. Research of changes in cultivated land in Xinjiang based on RS and GIS[J]. Progress in Geography，29（3）： 312-318.

鄧楚雄，李曉青，向云波，張曄，劉利科. 2013. 長株潭城市群地區耕地數量時空變化及其驅動力分析[J]. 經濟地理，33（6）：142-147.

Deng C X，Li X Q，Xiang Y B，Zhang Y，Liu L K. 2013. The spatio-temporal change and driving forces of cultivated land quantity in Chang-Zhu-Tan urban agglomeration[J]. Economic Geography，33（6）：142-147.

洪舒蔓，郝晉珉，周寧，陳麗，呂振宇. 2014. 黃淮海平原耕地變化及對糧食生產格局變化的影響[J]. 農業工程學報，30（21）：268-277.

Hong S M，Hao J M，Zhou N，Chen L，Lü Z Y. 2014. Change of cultivated land and its impact on grain production pattern in Huang-Huai-Hai plain[J]. Transactions of the Chinese Society of Agricultural Engineering，30（21）：268-277.

郇紅艷，孫君. 2012. 中部地區耕地非農化及其驅動因子的灰色關聯分析——以安徽省阜陽市為例[J]. 水土保持通報，32（1）：82-88.

Huan H Y，Sun J. 2012. Grey correlation analysis of cultivated land conversion and its driving factors in central China——A case study of Fuyang city， Anhui province[J]. Bulletin of Soil and Water Conservation，32（1）：82-88.

李艷. 2010. 安徽省耕地資源變化預測與保護研究[D]. 合肥：安徽農業大學.

Li Y. 2010. Research on prediction and protection of cultivated land change in Anhui province[D]. Hefei： Anhui Agricultural University.

梁海鷗. 2012. 巴彥縣耕地資源利用變化時空特征研究[D]. 哈爾濱：東北農業大學.

Liang H O. 2012. Research on spatial-temporal characteristics of cultivated land use change in Bayan[D]. Harbin：Northeast Agricultural University.

劉彥彤，張延軍，趙玲. 2011. 長春市耕地動態變化及其驅動力分析[J]. 地理科學，31（7）：868-873.

Liu Y T，Zhang Y J，Zhao L. 2011. Dynamic changes of cultiva-

ted land in Changchun city and their driving forces[J]. Scientia Geographica Sinica，31（7）：868-873.

王國剛，劉彥隨，陳秧分. 2014. 中國省域耕地集約利用態勢與驅動力分析[J]. 地理學報，69（7）：907-915.

Wang G G，Liu Y S，Chen Y F. 2014. Dynamic trends and dri-

ving forces of land-use intensification in the cultivated land of China[J]. Acta Geographica Sinica，69（7）：907-915.

王文剛，龐笑笑，宋玉祥，李汝資. 2012. 中國建設用地變化的空間分異特征[J]. 地域研究與開發，31（1）： 110-115.

Wang W G，Pang X X，Song Y X，Li R Z. 2012. The spatial differential features of construction land changes in China[J]. Areal Research and Development，31（1）：110-115.

王秀蘭，包玉海. 1999. 土地利用動態變化研究方法探討[J].地理科學進展，18（1）：81-87.

Wang X L，Bao Y H. 1999. Study on the methods of land use dynamic change research[J]. Progress in Geography，18（1）：81-87.

王曉文，田秉濤，李娟，尹春. 2014. 1985-2010年福建省耕地重心時空變化研究[J]. 亞熱帶資源與環境學報，9（2）：1-9.

Wang X W，Tian B T，Li J，Yin C. 2014. Temporal and spatial variations of cultivated land gravity centers in Fujian province during 1985 and 2010[J]. Journal of Subtropical Resources and Environment，9（2）：1-9.

徐秋艷. 2011. SPSS統計分析方法及應用實驗教程[M]. 北京：中國水利水電出版社：200-213.

Xu Q Y. 2011. SPSS Statistical Analysis Method and Application Experiment Course[M]. Beijing：China Water & Power Press： 200-213.

楊春艷，沈渭壽，王濤. 2015. 近30年西藏耕地面積時空變化特征[J]. 農業工程學報，31（1）：264-271.

Yang C Y，Shen W S，Wang T. 2015. Spatial-temporal characteristics of cultivated land in Tibet in recent 30 years[J]. Transactions of the Chinese Society of Agricultural Engineering，31（1）： 264-271.

張榮天，焦華富. 2015. 中國省際耕地利用效率時空格局分異與機制分析[J]. 農業工程學報，31（2）：277-287.

Zhang R T，Jiao H F. 2015. Spatial-temporal pattern differentiation and its mechanism analysis of using efficiency for provincial cultivated land in China[J]. Transactions of the Chinese Society of Agricultural Engineering，31（2）：277-287. 趙永華，劉曉靜，奧勇. 2013. 陜西省耕地資源變化及耕地壓力指數分析與預測[J]. 農業工程學報，29（11）：217-223.

Zhao Y H，Liu X J，Ao Y. 2013. Analysis of cultivated land change， pressure index and its prediction in Shaanxi province[J]. Transactions of the Chinese Society of Agricultural Engineering，29（11）：217-223.

Barbier E B，Burgess J C，Grainger A. 2010. The forest transition： Towards a more comprehensive theoretical framework[J]. Land Use Policy，27（2）：98-107.

Ferreira M P， Alves D S， ShimabukuroY E. 2015. Forest dynamics and land-use transitions in the Brazilian Atlantic Forest： the case of sugarcane expansion[J]. Region Environmental Change，15（2）：365-377.

Lambin E F，Meyfroidt P，Rudel T K，Schneider L，Uriarte M.2010. Land use transitions： Socio-ecological feedback versus socio-economic change[J]. Land Use Policy，27（2）：108-118.

（責任編輯 鄧慧靈）