昆山市耕地總量動(dòng)態(tài)平衡邊界研究

程紹鉑 ，楊桂山 ，吳建楠 ，胡 艷

（1.安徽大學(xué)經(jīng)濟(jì)學(xué)院，安徽合肥 230601；2.中國(guó)科學(xué)院南京地理與湖泊研究所，江蘇南京210008；3.安徽大學(xué)長(zhǎng)三角經(jīng)濟(jì)與社會(huì)發(fā)展研究中心，安徽合肥230601；4.中山大學(xué)廣東決策科學(xué)研究院，廣東廣州510000）

如何在現(xiàn)有的國(guó)情條件下和日益開(kāi)放的國(guó)際環(huán)境下，在實(shí)現(xiàn)工業(yè)化和城市化過(guò)程中保持農(nóng)業(yè)可持續(xù)發(fā)展保護(hù)有限的耕地資源，是中國(guó)經(jīng)濟(jì)和社會(huì)可持續(xù)發(fā)展的全局性戰(zhàn)略問(wèn)題[1]。中國(guó)政府在歷經(jīng)1978—1985年耕地保護(hù)制度恢復(fù)期、1986—1997年的耕地保護(hù)制度發(fā)展期和1998至今的以國(guó)土資源部成立為標(biāo)志的耕地保護(hù)制度完善期[2]，先后出臺(tái)了一系列耕地保護(hù)政策，在相當(dāng)大的程度上實(shí)現(xiàn)了耕地保護(hù)目標(biāo)，也為社會(huì)經(jīng)濟(jì)尤其是農(nóng)業(yè)生產(chǎn)的可持續(xù)發(fā)展提供了必要的保障[3]。耕地總量動(dòng)態(tài)平衡政策作為耕地保護(hù)重要政策之一，在一定程度上保證了耕地總量的平衡，有助于遏制耕地加速流失。但是，由于目前中國(guó)正處于城市化發(fā)展和經(jīng)濟(jì)發(fā)展的快速時(shí)期，對(duì)建設(shè)用地的需求不斷增長(zhǎng)，給耕地總量動(dòng)態(tài)平衡政策的實(shí)施帶來(lái)了巨大的挑戰(zhàn)。由于國(guó)情和政治體制的不同，國(guó)外對(duì)耕地總量動(dòng)態(tài)平衡邊界問(wèn)題的研究較少。國(guó)內(nèi)由黃賢金等率先開(kāi)展了耕地總量動(dòng)態(tài)平衡邊界研究[1,4]，但主要集中于區(qū)域尺度。實(shí)際上，在更小的空間尺度上研究耕地總量動(dòng)態(tài)平衡邊界問(wèn)題，更有助于加強(qiáng)對(duì)耕地總量動(dòng)態(tài)平衡實(shí)現(xiàn)規(guī)律的認(rèn)識(shí)，有助于地方政府對(duì)不同土地經(jīng)濟(jì)發(fā)展階段的把握，從而制定出適合于當(dāng)?shù)匕l(fā)展情況的耕地保護(hù)對(duì)策。本文以昆山市為研究對(duì)象，研究在現(xiàn)有耕地保護(hù)政策實(shí)施力度保持不變趨勢(shì)下的耕地總量動(dòng)態(tài)平衡邊界，有助于從實(shí)際情況出發(fā)，衡量耕地總量動(dòng)態(tài)平衡政策的實(shí)施效果，增強(qiáng)政策的可操作性，高效地實(shí)現(xiàn)耕地保護(hù)的目標(biāo)。同時(shí)，本研究對(duì)中國(guó)其他正處于經(jīng)濟(jì)高速發(fā)展面臨耕地保護(hù)問(wèn)題的縣級(jí)市均有借鑒意義。

1 研究區(qū)概況

昆山市地處江蘇省東南部，土地總面積928 km2，其中水域面積占29.92%。昆山市經(jīng)濟(jì)發(fā)達(dá)，據(jù)昆山市政府網(wǎng)站公布數(shù)據(jù)顯示，2010年全市地區(qū)生產(chǎn)總值達(dá)到2100.28億元，成為國(guó)內(nèi)首個(gè)GDP總量超過(guò)2000億元的縣級(jí)市，同時(shí)按常住人口計(jì)算人均地區(qū)生產(chǎn)總值達(dá)到14.45萬(wàn)元。改革開(kāi)放以來(lái)，昆山市工業(yè)化和城市化進(jìn)程迅速，成為以開(kāi)放型經(jīng)濟(jì)為主導(dǎo)，三大產(chǎn)業(yè)協(xié)調(diào)發(fā)展的沿海新興城市之一。在國(guó)家統(tǒng)計(jì)局公布的全國(guó)百?gòu)?qiáng)縣（市）排名中，昆山市位居全國(guó)百?gòu)?qiáng)縣之首。昆山市高速的經(jīng)濟(jì)發(fā)展引起了快速的土地利用變化，成為一個(gè)典型的縣市級(jí)尺度的研究區(qū)域，為揭示經(jīng)濟(jì)高速發(fā)展和城市化過(guò)程中的耕地變化提供了一個(gè)難得的案例[4]。分析經(jīng)濟(jì)社會(huì)快速發(fā)展的昆山市耕地總量動(dòng)態(tài)平衡邊界，對(duì)其他城市的耕地平衡政策的實(shí)施具有啟示性意義。

2 昆山市耕地變化態(tài)勢(shì)分析與影響因素分析

根據(jù)昆山市統(tǒng)計(jì)年鑒資料，昆山市1985—2009年間耕地?cái)?shù)量數(shù)量變化如圖1所示。昆山市耕地總量1985—2009年間呈減少趨勢(shì)，25年間總共減少了68.2%，2009年耕地面積僅相當(dāng)于1985年耕地面積的32.8%。除1999年耕地面積略有上升外，其他年份耕地面積均不斷下降，耕地面積的減少是一個(gè)長(zhǎng)期積累的過(guò)程。

圖1 昆山市耕地?cái)?shù)量變化趨勢(shì)圖Fig.1 The trend of cultivated land change

耕地面積變化受多種因素影響，昆山市自然因素對(duì)耕地面積的影響較小，耕地面積變化主要受政策因素和社會(huì)經(jīng)濟(jì)因素影響。根據(jù)昆山市實(shí)際情況以及資料的可獲得性和可靠性，并且參考已有學(xué)者的研究成果[5-8]，選用總戶籍人口X1、非農(nóng)業(yè)人口X2、地區(qū)生產(chǎn)總值X3、人均地區(qū)生產(chǎn)總值X4、第二產(chǎn)業(yè)生產(chǎn)總值X5、糧食總產(chǎn)量X6、社會(huì)固定資產(chǎn)投資額X7、城市化率X8等8種指標(biāo)來(lái)反映社會(huì)經(jīng)濟(jì)變化對(duì)耕地面積的影響，數(shù)據(jù)均來(lái)自于昆山市統(tǒng)計(jì)年鑒；選用耕地保護(hù)政策實(shí)施強(qiáng)度X9反映耕地面積對(duì)保護(hù)政策實(shí)施力度的反饋。為了量化政策因素對(duì)耕地面積的影響，按照已有專家的評(píng)價(jià)方法[9]并根據(jù)昆山市的實(shí)際情況，將不同時(shí)期耕地保護(hù)政策實(shí)施的強(qiáng)度分別定義：將1985年前耕地保護(hù)政策強(qiáng)度賦值為“-1”，此時(shí)為弱耕地保護(hù)狀態(tài)；將1986—1994年賦值為“-3”，這一時(shí)期開(kāi)始實(shí)施《土地管理法》，可以定義為較強(qiáng)耕地保護(hù)狀態(tài)；將1995—2003年賦值為“-5”，開(kāi)始實(shí)施基本農(nóng)田保護(hù)區(qū)制度，定義為強(qiáng)耕地資源保護(hù)狀態(tài)；將2003年后定義為“-8”，雖然自1997年開(kāi)始實(shí)施土地用途管制制度，但結(jié)合昆山市具體情況以及考慮制度實(shí)施的時(shí)間滯后問(wèn)題，將2003年后定義為最嚴(yán)格的耕地資源保護(hù)狀態(tài)。

從昆山市社會(huì)經(jīng)濟(jì)因素和政策因素驅(qū)動(dòng)變量與耕地?cái)?shù)量的相關(guān)分析可以看出，各驅(qū)動(dòng)變量與耕地?cái)?shù)量均呈極顯著相關(guān)，說(shuō)明所選指標(biāo)的合理性，這些指標(biāo)能夠反映耕地變化的驅(qū)動(dòng)。其中，總戶籍人口、非農(nóng)業(yè)人口、地區(qū)生產(chǎn)總值、人均地區(qū)生產(chǎn)總值、第二產(chǎn)業(yè)生產(chǎn)總值、社會(huì)固定資產(chǎn)投資額和城市化率與耕地?cái)?shù)量呈負(fù)相關(guān)，即這些變量增加會(huì)導(dǎo)致耕地?cái)?shù)量的減少。糧食總產(chǎn)量和耕地保護(hù)政策實(shí)施強(qiáng)度與耕地?cái)?shù)量呈正相關(guān)，即這些變量值越大，耕地面積越大。從耕地保護(hù)政策實(shí)施強(qiáng)度與耕地?cái)?shù)量呈現(xiàn)極顯著相關(guān)，且相關(guān)系數(shù)達(dá)到0.937，可以看出，耕地保護(hù)政策實(shí)施力度的增強(qiáng)能夠有效遏制耕地?cái)?shù)量的減少。

表1 昆山市耕地?cái)?shù)量與驅(qū)動(dòng)變量的相關(guān)分析結(jié)果Tab.1 Correlation result of cultivated land and driving factors of Kunshan

以耕地?cái)?shù)量為因變量，表1中所列的驅(qū)動(dòng)變量為自變量，對(duì)耕地?cái)?shù)量進(jìn)行線性回歸擬合。得到昆山市耕地?cái)?shù)量變化模型為：

模型的復(fù)相關(guān)系數(shù)為0.997，調(diào)整后的復(fù)相關(guān)系數(shù)的平方為0.991，表示自變量總體對(duì)因變量的解釋程度為99.1%，也就是說(shuō)，所選的自變量幾乎完全解釋了耕地?cái)?shù)量的變化。模型的顯著性檢驗(yàn)中F值為302.874，通過(guò)檢驗(yàn)。實(shí)際上，由于所有的變量都是時(shí)間序列變量，本身隨時(shí)間的變化具有一定的規(guī)律，使得變量間相關(guān)程度增加[10]。因此，抽取其中少數(shù)幾個(gè)自變量的組合對(duì)耕地?cái)?shù)量進(jìn)行多元線性回歸，仍然能得到較好的結(jié)果：

模型的復(fù)相關(guān)系數(shù)為0.987，調(diào)整后的復(fù)相關(guān)系數(shù)的平方為0.970，這表示自變量總體對(duì)因變量的解釋程度為97.0%，表明地區(qū)生產(chǎn)總值X3、第二產(chǎn)業(yè)生產(chǎn)總值X5和耕地保護(hù)政策實(shí)施強(qiáng)度X9這3個(gè)自變量基本上能夠完全解釋耕地?cái)?shù)量的變化。模型的顯著性檢驗(yàn)中F值為227.229，通過(guò)檢驗(yàn)。總體來(lái)看，模型的擬合效果較好。也就是說(shuō)，這3個(gè)自變量能夠反映與耕地?cái)?shù)量變化的伴隨關(guān)系，用地區(qū)生產(chǎn)總值、第二產(chǎn)業(yè)生產(chǎn)總值和政策保護(hù)實(shí)施強(qiáng)度來(lái)表示多因素綜合作用下驅(qū)動(dòng)耕地?cái)?shù)量變化的機(jī)制，是比較符合實(shí)際的，用此式替代式1是可行的。

這種驅(qū)動(dòng)力模型雖然解釋程度較高，能夠通過(guò)所選的自變量基本上完全解釋耕地面積的變化，但這是在不考慮耕地變化的邊際收益傾向隨時(shí)間變化的前提下得到的，實(shí)際上隨著經(jīng)濟(jì)的不斷發(fā)展與技術(shù)的進(jìn)步，耕地變化的邊際收益隨時(shí)間不斷減少，呈現(xiàn)出逐漸遞減的趨勢(shì)。

以式2中的經(jīng)濟(jì)發(fā)展驅(qū)動(dòng)變量地區(qū)生產(chǎn)總值X3和第二產(chǎn)業(yè)生產(chǎn)總值X5來(lái)考察耕地面積變化和經(jīng)濟(jì)發(fā)展之間的數(shù)量轉(zhuǎn)換關(guān)系。以“Δ耕地/Δ地區(qū)生產(chǎn)總值”和“Δ耕地/Δ第二產(chǎn)業(yè)生產(chǎn)總值”作為經(jīng)濟(jì)發(fā)展對(duì)耕地資源占用的邊際傾向分析指標(biāo)，考察這兩個(gè)變量隨時(shí)間的變化趨勢(shì)，通過(guò)這一分析可以看出昆山市地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值每增加一單位，耕地?cái)?shù)量的變化情況。

由昆山市耕地變化的邊際占用傾向趨勢(shì)分析圖2和圖3可以看出，雖然由于昆山市耕地變化的邊際收益傾向出現(xiàn)階段性的波動(dòng)，但是總體上遞減的趨勢(shì)還是主要的，增加一個(gè)單位地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值所需占用的耕地資源總體上呈現(xiàn)出減少的趨勢(shì)，表明昆山市經(jīng)濟(jì)發(fā)展正逐漸向節(jié)約和集約利用土地方向發(fā)展。

3 昆山市耕地資源邊際占用傾向的經(jīng)濟(jì)邊界分析

昆山市地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值的增加對(duì)耕地資源的邊際占用傾向是遞減的，在維持目前相對(duì)穩(wěn)定高速的經(jīng)濟(jì)發(fā)展趨勢(shì)與最為嚴(yán)格的耕地保護(hù)政策實(shí)施強(qiáng)度不變的前提下，可以認(rèn)為當(dāng)邊際數(shù)值接近于零的時(shí)候，耕地資源將不再被占用。此時(shí)的耕地面積為耕地資源占用的經(jīng)濟(jì)邊界。

3.1 耕地資源的地區(qū)生產(chǎn)總值邊際占用傾向的邊界分析

以耕地?cái)?shù)量為因變量Y，以Δ耕地/Δ地區(qū)生產(chǎn)總值為自變量X建立邊界分析函數(shù)：

式3中，決定系數(shù)R=0.843，說(shuō)明方程的擬合效果不錯(cuò)。當(dāng)X=0時(shí)的得到昆山市耕地面積的地區(qū)生產(chǎn)總值發(fā)展邊界，此時(shí)耕地面積為250843.7 hm2。也就是說(shuō)，在不考慮其他因素的影響下，當(dāng)耕地面積減少到250843.7 hm2之后，耕地面積就不再隨地區(qū)生產(chǎn)總值的增加而減少。這一面積相當(dāng)于2006年昆山市耕地面積的83.2%，相當(dāng)于在2006年的基礎(chǔ)上耕地再減少50727.27 hm2。

以時(shí)間為因變量T，Δ耕地/Δ地區(qū)生產(chǎn)總值為自變量X建立邊界分析函數(shù)：

式4中，決定系數(shù)R=0.826。當(dāng)X=0時(shí)可以等到昆山市耕地面積的地區(qū)生產(chǎn)總值邊際占用停滯時(shí)間為2011.8，也就是會(huì)說(shuō)在只考慮地區(qū)生產(chǎn)總值對(duì)耕地占用的情況下，在2011—2012年耕地資源將停止被占用。

圖2 昆山市地區(qū)生產(chǎn)總值的耕地資源邊際占用傾向趨勢(shì)圖Fig.2 The trend of the marginal effect on cultivated land occupation due to GDP growth

圖3 昆山市第二產(chǎn)業(yè)生產(chǎn)總值的耕地資源邊際占用傾向趨勢(shì)圖Fig.3 The trend of the marginal effect on cultivated occupation due to Gop growth in the secondary sector

3.2 耕地資源的第二產(chǎn)業(yè)生產(chǎn)總值邊際占用傾向的邊界分析

以耕地面積為因變量Y，以Δ耕地/Δ第二產(chǎn)業(yè)生產(chǎn)總值為自變量X建立邊界分析函數(shù)，得到：

式5中，決定系數(shù)R=0.786。當(dāng)X=0時(shí)，得到昆山市耕地面積的第二產(chǎn)業(yè)生產(chǎn)總值發(fā)展邊界，此時(shí)的耕地面積為276933.2 hm2。即在不考慮其他因素的影響下，當(dāng)耕地面積減少到276933.2 hm2之后，耕地面積就不再隨第二產(chǎn)業(yè)生產(chǎn)總值的增加而減少。這一面積相當(dāng)于昆山市2006年耕地面積的91.8%，比昆山市地區(qū)生產(chǎn)總值邊界的結(jié)果要小，說(shuō)明相比于地區(qū)生產(chǎn)總值昆山市第二產(chǎn)業(yè)生產(chǎn)總值更早達(dá)到發(fā)展邊界，工業(yè)用地更早的實(shí)現(xiàn)土地資本的存量調(diào)劑，零增量占用耕地資源的發(fā)展模式。

以時(shí)間為自變量T，Δ耕地/Δ第二產(chǎn)業(yè)生產(chǎn)總值為自變量X建立邊界分析函數(shù)：

式6中，決定系數(shù)R=0.843，擬合效果較好。當(dāng)X=0時(shí)可以得到昆山市耕地面積的第二產(chǎn)業(yè)生產(chǎn)總值邊際占用停滯時(shí)間為2010.9。也就是會(huì)說(shuō)在只考慮第二產(chǎn)業(yè)生產(chǎn)總值對(duì)耕地占用的情況下，在2010—2011年間耕地資源可能會(huì)停止被占用。這一停滯時(shí)間同樣比由地區(qū)生產(chǎn)總值停滯邊界得到的停滯時(shí)間要早，進(jìn)一步證實(shí)了相對(duì)于地區(qū)生產(chǎn)總值來(lái)說(shuō)，昆山市第二產(chǎn)業(yè)生產(chǎn)總值更早達(dá)到發(fā)展邊界。

3.3 耕地資源占用態(tài)勢(shì)分析

通過(guò)建立耕地?cái)?shù)量隨時(shí)間變化的回歸模型來(lái)驗(yàn)證前面對(duì)耕地占用停滯邊界的分析。以年份為因變量T，耕地面積為自變量X建立函數(shù)得到：

式7中，決定系數(shù)R=0.961，擬合效果很好。令耕地面積為地區(qū)生產(chǎn)總值占用邊界值即X=250843.7時(shí)得到T=2012.93，這個(gè)值與由式4得到的只考慮地區(qū)生產(chǎn)總值對(duì)耕地占用的情況下耕地會(huì)停止被占用的時(shí)間接近。可以認(rèn)為式3—4與本式得到了互為印證，可以認(rèn)為昆山市耕地?cái)?shù)量的地區(qū)生產(chǎn)總值邊界值為250843.7 hm2，當(dāng)耕地?cái)?shù)量減少到此邊界值之后，耕地?cái)?shù)量不再隨地區(qū)生產(chǎn)總值的增加而減少，其達(dá)到時(shí)間大致為2012年左右。

令耕地面積為第二產(chǎn)業(yè)生產(chǎn)總值占用邊界值即X=276933.2時(shí)得到T=2011.2，這個(gè)值與由式6得到的只考慮第二產(chǎn)業(yè)生產(chǎn)總值對(duì)耕地占用的情況下耕地會(huì)停止被占用的時(shí)間非常接近。可以認(rèn)為式5—6與式7得到了相互印證，因此基于本分析可以認(rèn)為昆山市耕地?cái)?shù)量的第二產(chǎn)業(yè)生產(chǎn)總值邊界值為276933.2 hm2，當(dāng)耕地?cái)?shù)量減少到此邊界值之后，耕地?cái)?shù)量不再隨地區(qū)生產(chǎn)總值的增加而減少，其達(dá)到時(shí)間大致為2011年。

4 結(jié)論

單從經(jīng)濟(jì)效益的角度來(lái)看，耕地的經(jīng)濟(jì)產(chǎn)出效益遠(yuǎn)小于建設(shè)用地經(jīng)濟(jì)效益，所以耕地被占用轉(zhuǎn)變?yōu)榻ㄔO(shè)用地是符合經(jīng)濟(jì)理性的，屬于自發(fā)、主動(dòng)的經(jīng)濟(jì)行為。但隨著耕地占用成本的不斷增加、科學(xué)技術(shù)的進(jìn)步、區(qū)域產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整以及政府耕地保護(hù)政策實(shí)施力度的加強(qiáng)，耕地占用的邊際收益傾向正在不斷減小，昆山市經(jīng)濟(jì)發(fā)展正逐漸向節(jié)約和集約利用土地方向發(fā)展。從地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值來(lái)看，昆山市經(jīng)濟(jì)發(fā)展將由大量占用土地資本向相對(duì)高效占用土地資本轉(zhuǎn)變，且最終將實(shí)現(xiàn)土地資本存量調(diào)劑，經(jīng)濟(jì)發(fā)展零占用耕地資源的發(fā)展模式。單從地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值發(fā)展角度來(lái)說(shuō)，在維持目前的最為嚴(yán)格的耕地保護(hù)政策實(shí)施強(qiáng)度不變的情況下，昆山市耕地?cái)?shù)量將在24637.8—250843.7 hm2之間達(dá)到邊界值。此后地區(qū)生產(chǎn)總值和第二產(chǎn)業(yè)生產(chǎn)總值的發(fā)展將不再依賴于耕地占用，達(dá)到邊界值的時(shí)間大致為2011—2012年之間，且第二產(chǎn)業(yè)生產(chǎn)總值的耕地占用比地區(qū)生產(chǎn)總值的耕地占用更早達(dá)到邊界值。屆時(shí)，昆山市耕地總量動(dòng)態(tài)平衡有望得到更為高效的實(shí)現(xiàn)。

當(dāng)然，這只是在不考慮耕地省內(nèi)異地占補(bǔ)平衡以及假設(shè)耕地保護(hù)政策實(shí)施強(qiáng)度維持不變的情況下得出的，事實(shí)上耕地占用的實(shí)際情況更為復(fù)雜，并且當(dāng)經(jīng)濟(jì)社會(huì)進(jìn)入一個(gè)新的發(fā)展時(shí)期后，這一預(yù)測(cè)結(jié)果將發(fā)生相應(yīng)的變化。從昆山市近25年來(lái)耕地面積變化與經(jīng)濟(jì)高速增長(zhǎng)的相互關(guān)系中可以看出，與經(jīng)濟(jì)增長(zhǎng)條件下環(huán)境的庫(kù)茲涅茲倒U形曲線相似，耕地面積的減少也存在著同樣的動(dòng)態(tài)平衡點(diǎn)。該平衡點(diǎn)即為曲線的拐點(diǎn)，動(dòng)態(tài)平衡點(diǎn)達(dá)到的時(shí)間除受到社會(huì)經(jīng)濟(jì)發(fā)展變化的影響外，與耕地保護(hù)政策的實(shí)施力度也有很大關(guān)系。同時(shí)耕地保護(hù)政策實(shí)施力度也會(huì)對(duì)社會(huì)經(jīng)濟(jì)發(fā)展產(chǎn)生影響，如何建立它們之間的反饋關(guān)系，合理調(diào)整耕地保護(hù)政策，使其能盡快達(dá)到耕地總量的動(dòng)態(tài)平衡點(diǎn)，是今后需要進(jìn)一步研究的問(wèn)題。

（

）：

［1］ 翟文俠，黃賢金.我國(guó)耕地保護(hù)政策運(yùn)行效果分析［J］.中國(guó)土地科學(xué)，2003，17（2）：8-13.

［2］ 張效軍，歐名豪，李景剛.我國(guó)耕地保護(hù)制度變遷及其績(jī)效分析［J］.社會(huì)科學(xué)，2007，（8）：13-20.

［3］ 徐夢(mèng)潔，陳慧中，趙紫玉，等.耕地總量動(dòng)態(tài)平衡評(píng)價(jià)模型初探［J］.資源科學(xué)，2009，（12）：2137-2142.

［4］ 黃賢金，濮勵(lì)杰，彭補(bǔ)拙.城市土地利用變化及其響應(yīng)：模型構(gòu)建與實(shí)證研究［M］.北京：科學(xué)出版社，2008.

［5］ 楊桂山.長(zhǎng)江三角洲近50年耕地?cái)?shù)量變化的過(guò)程及驅(qū)動(dòng)機(jī)制研究［J］.自然資源學(xué)報(bào)，2001，16（2）.

［6］ Pu Li-jie,Yang Gui-shan,Xu Meng-jie.Characteristics and Process of Land Use Changes in the Yangtze River Delta,China［J］.Pedosphere,2001,11（3）:193-198.

［7］ 濮勵(lì)杰，黃賢金，彭補(bǔ)拙.耕地總量動(dòng)態(tài)平衡政策體系研究［J］.南京大學(xué)學(xué)報(bào)（自然科學(xué)），2001，37（6）：757-765.

［8］ 丁洪建，吳次芳，梁留科.耕地保護(hù)理念的創(chuàng)新研究［J］.中國(guó)土地科學(xué)，2002，16（4）：14-19.

［9］ 黃賢金，濮勵(lì)杰，周峰，等.長(zhǎng)江三角洲地區(qū)耕地總量動(dòng)態(tài)平衡政策目標(biāo)實(shí)現(xiàn)的可能性分析［J］.自然資源學(xué)報(bào)，2002，17（6）：670-676.

［10］姚鑫，楊桂山，萬(wàn)榮榮.昆山市 22年來(lái)耕地流失研究［J］.長(zhǎng)江流域資源與環(huán)境，2010，19（9）：1063-1068.