糧食安全視角下長三角耕地利用格局變化及預測分析

曹宗龍 劉錦洋

（1江蘇蘇地仁合土地房地產資產評估測繪造價咨詢有限公司，江蘇南京 210000；2安徽師范大學地理與旅游學院，安徽蕪湖 241000）

“十三五”規劃提出，堅持最嚴格的耕地保護制度，堅守耕地紅線，實施“藏糧于地、藏糧于技”戰略[1]。耕地是保障國家糧食安全的基石，一方面要保證耕地的數量，另一方面要提高耕地的質量[2]。中國的耕地變化一直以來備受關注[3]，國內外學者基于遙感影像、“3S”技術等[4]，將耕地保護與利用與景觀格局、生態脆弱性、生物多樣性、耕地破碎性、耕地生態安全[5]等領域相交互，通過構建預測模型、驅動力分析[6]等，形成了相對完善的評價體系[7]。作為六大世界級城市群之一的長三角，不僅是參與國際競爭的重要平臺，也是重要的糧食基地，但當前對于長三角耕地變化的研究仍相對較少。本文從面積、形態、類型、格局、預測對長三角耕地進行研究，為長三角糧食安全、耕地保護等相關政策制定提供支撐。

1 研究區概況與數據來源

1.1 研究區概況

長江三角洲城市群（簡稱長三角城市群），位于長江入海之前的沖積平原。長三角城市群面積僅占全國總面積2.3%，卻擁有2.25億人口，并貢獻了全國25%左右的GDP。

長三角城市群是“一帶一路”與長江經濟帶的重要交匯地帶，在中國國家現代化建設大局和開放格局中具有舉足輕重的戰略地位。

長三角作為六大世界級城市群之一，是中國參與國際競爭的重要平臺、經濟社會發展的重要引擎。長三角城市群要建設面向全球、輻射亞太、引領全國的世界級城市群，建成最具經濟活力的資源配置中心、具有全球影響力的科技創新高地、全球重要的現代服務業和先進制造業中心、亞太地區重要國際門戶、全國新一輪改革開放排頭兵、美麗中國建設示范區。

1.2 數據來源

該研究采用GLOBELAND 30數據，對長三角2000、2010、2020年土地利用數據進行拼接、裁剪、融合等相關處理（圖1）。

圖1 長三角2000、2010、2020年土地利用

2 方法與模型

2.1 耕地利用動態度

耕地動態度可表征某個時段研究區或格網內耕地變化幅度與速度[8]。計算公式：

式中，Sa、Sb分別為研究期初與期末的耕地面積（km2）；T為研究時間,當T設定為年時，K為研究時段的耕地面積動態度。

2.2 景觀格局指數法

利用景觀分析軟件Fragstats 3.3軟件從景觀水平對土地利用景觀格局進行分析，得出不同時期的土地利用景觀格局和變化趨勢[9]。選用6種具有代表性的景觀指數：景觀面積（TA）、拼塊所占景觀面積的比例（PLAND）、最大斑塊所占景觀面積的比例（LPI）、平均分維數（FRAC_MN）、分離度（SHLIT）、聚合度（AI）。

（1）拼塊類型面積（TA），單位：hm2，范圍：TA＞0。

描述：TA等于一個景觀的總面積，除以10 000后轉化為hm2。

（2）拼塊所占景觀面積的比例（PLAND），單位：%，范圍：0＜PLAND≤100。

描述：PLAND等于1個街區類型的總面積占總景觀面積的百分比。當其值趨近于0時，表明該地塊類型在景觀中非常罕見。當其值為100時，表示整個景觀僅由1種地塊構成。

（3）最大斑塊所占景觀面積的比例（LPI），單位：百分比，范圍：0＜LPI≤100。

描述：LPI等于斑塊類型中最大斑塊所占景觀總面積的百分比。

（4）平均分維數（FRAC_MN），單位：無，范圍：1.0≤FRAC_MN≤2.0。

公式描述：D=2ln(P/4)/ln(A)，D表示分維數，P為斑塊周長，A為斑塊面積。D值越小，表明斑塊形狀越簡單，D值的理論范圍為1.0～2.0，1.0代表形狀最簡單的正方形斑塊，2.0表示等面積下周邊最復雜的斑塊。

（5）分離度（SHLIT），單位：%，范圍：0＜SHLIT≤100。

描述：SPLIT等于景觀面積的平方除以所有斑塊面積的平方之和，其中這些斑塊屬于相關斑塊類型。

（6）聚合度（AI），單位：%，范圍：0＜AI≤100。

描述：AI等于分裂類似鄰接的所有類型，包括相關的塊類相似的鄰接的最大數量可能包括相關的地塊類型，乘以之和的比例值的這種類型的景觀組成，最后乘以100%。

3 結果與分析

3.1 長三角耕地利用變化分析

由表1可知，長三角2000年耕地總面積為19 764 321 hm2，2010年為19 131 997 hm2，2020年為17 571 185 hm2。從2000年至2010年長三角耕地減少632 324 hm2，而自2010年到2020年長三角耕地減少1 560 812 hm2，由此可見。2000—2020年，長三角耕地數量持續減少，尤其2010—2020年，長三角耕地數量大幅度減少。

結合表1、圖2、圖3可知，2000年安徽省耕地總面積為8 326 144 hm2，2010年為8 201 182 hm2，2000—2010年安徽省的耕地動態度為-1.50%，而2020年安徽省的耕地數量為7 952 145 hm2，2010—2020年安徽省的耕地動態度為-3.04%。安徽省作為長三角中耕地面積的最大的省份，對糧食安全起到重要的作用。整體上來看，安徽省耕地數量雖然持續下降，但下降幅度變化不大，耕地利用動態度呈低強度減少。耕地在空間上的變化主要體現為城市擴張，耕地減少（圖4）。

圖4 安徽省3期耕地分布

表1 長三角耕地數量 單位：hm2

圖2 2000—2010年耕地利用動態度

圖3 2010—2020年耕地利用動態度

江蘇省2000年耕地總面積為7 511 408 hm2，2010年江蘇省耕總面積為7 237 584 hm2，2000—2010年江蘇省的耕地動態度為-3.65%，而2020年江蘇省的耕地數量為6 482 350 hm2，2010—2020年江蘇省的耕地動態度為-10.43%。整體上看，江蘇省耕地數量不斷下降（圖5），且2010—2020年耕地數量減少較大，耕地利用動態度呈中強度減少。

圖5 江蘇省3期耕地分布

浙江省2000年耕地總面積為3 465 255 hm2，2010年浙江省耕總面積為3 309 669 hm2，2000—2010年浙江省的耕地動態度為-4.49%，而2020年浙江省的耕地數量為2 811 009 hm2，2010—2020年浙江省的耕地動態度為-15.07%。整體上看，浙江省耕地數量不斷下降（圖6），耕地利用動態度呈不斷擴大趨勢，2010—2020年耕地數量減少十分明顯，耕地利用動態度呈高強度減少。

圖6 浙江省3期耕地分布

上海市2000年耕地總面積為461 514 hm2，2010年上海市耕地總面積為383 562 hm2，2020年上海市耕地總面積為325 681 hm2。上海市耕地數量較少，主要承載經濟功能。整體上看，上海市耕地數量減少（圖7），但耕地利用動態度呈微強度增加。

圖7 上海市3期耕地分布

由圖8可知，耕地的減少形態主要有以下方面：以上海市為首的耕地炸裂式減少；以長江經濟帶臨岸城市，沿長江呈帶狀式減少；以快速發展的合肥市、寧波市為首的耕地星狀式減少。這可以看出耕地的減少形態主要與經濟的發展息息相關。2000—2020年，耕地的減少主要集中在長江沿岸、沿海城市以及掌握優質資源快速發展的城市，發展導致城市擴張，耕地減少。

圖8 長三角耕地變化

3.2 長三角耕地轉移矩陣分析

2000—2010年，長三角耕地未發生變化的面積為18 007 381 hm2，耕地主要向建設用地、林地、水域方面轉移。同時，林地、水域、建設用地向耕地轉移的數量也較大。通過轉入轉出的數據進行對比，林地、草地、建設用地由耕地轉出數量大于轉入耕地數量。水域、未利用地轉入數量大于轉出數量。整體上看，林地、草地、未利用地，轉入和轉出量基本趨于平衡。水域轉入耕地數量較大于轉出數量，而耕地轉入建設用地數量明顯大于建設用地轉入數量，差額為665 187 hm2，轉移數量顯著。由此可見，2000—2010年，主要是由耕地向建設用地方面轉移，具體見表2。

表2 2000—2010年耕地轉入轉出 單位：hm2

2010—2020年長三角耕地未發生變化的面積為16 297 373 hm2，耕地主要向建設用地、水域、林地、草地方面轉移。建設用地、林地、水域向耕地轉移的數量也較大。通過轉入轉出的數據進行對比，林地、草地、水域、建設用地，由耕地轉出數量大于轉入耕地數量。僅未利用地轉入耕地數量大于轉出數量，且數量較少，這也體現了耕地整體向其他地類轉移，意味著耕地數量大幅度減少。整體上看，草地、未利用地，轉入和轉出量基本趨于平衡。林地、水域、建設用地由耕地轉出數量均較明顯大于轉入耕地數量。尤其是建設用地，差額為1 306 391 hm2，較2000—2010年的凈轉出量也明顯提升。由此可見，2000—2020年，主要是由耕地向建設用地方面轉移。尤其是2010—2020年，轉移數量十分明顯，具體見表3。

表3 2010—2020年耕地轉入轉出單位：hm2

3.3 長三角耕地景觀格局分析

基于GLOBELAND 30中2000、2010、2020年3期土地利用數據，將其放到Fragstats 3.3中class層級進行處理，選用6種具有代表性的景觀指數：景觀面積（TA）、拼塊所占景觀面積的比例（PLAND）、最大斑塊所占景觀面積的比例（LPI）、平均分維數（FRAC_MN）、分離度（SHLIT）、聚合度（AI）[5]，得到2000、2010、2020年關于耕地的相關景觀指數。

耕地景觀規模主要是研究景觀面積（TA）、拼塊所占景觀面積的比例（PLAND）、最大斑塊所占景觀面積的比例（LPI）。從表4可以看出，2000—2020年長三角耕地景觀面積（TA）、拼塊所占景觀面積的比例（PLAND）都在不斷下降，耕地景觀數量減少，同時所占比例減少。最大斑塊所占景觀面積的比例（LPI）也在不斷下降。這表現為長三角耕地的規模經營正在不斷削弱，集約利用耕地、保護利用耕地有待加強。

耕地景觀分離度主要是研究分離度（SHLIT）、聚合度（AI）。從表4可以看出，2000—2020年長三角耕地聚合度（AI）不斷下降，而分離度（SHLIT）正不斷上升。這表現為長三角耕地的細碎化程度不斷上升，耕地正在被拆分為更細小的地塊，耕地處于不斷離散狀態。

景觀形狀指數主要研究平均分維數（FRAC_MN）。從表4可以看出，2000—2020年長三角耕地平均分維數（FRAC_MN）不斷下降，逐漸向1.0趨近，這表明長三角耕地的形狀特征變得更加簡單，耕地斑塊向正方形斑塊趨近。

表4 耕地景觀格局變化

3.4 長三角耕地預測分析

根據2000、2010、2020年長三角耕地分布數據，對長三角3期耕地進行質心提取，2000年長三角的耕地質心在南京市的建鄴區，位于長江以南。2010年長三角耕地重心轉移至南京市浦口區中東部邊緣，并由長江以南轉向長江以北。2020年繼續向西北部轉移，轉移至南京市浦口區中部。根據2000、2010、2020年3期數據進行二階擬合分析，得到2030年耕地重心的經度為118.463 631°，緯度為32.121 696°（圖9）。

圖9 長三角耕地預測

預測2030年長三角耕地重心可能轉移至安徽省滁州市南譙區（圖10）。整體上看，長三角耕地重心一直向西北方向轉移，由長江以南轉移至長江以北。至2030年長三角耕地重心可能由江蘇省轉移至安徽省，由南京市轉移至滁州市。東南方向耕地減少較快，導致整個耕地格局不斷向西北方向進行轉移，并跨過長江。

圖10 長三角耕地重心轉移

4 結論與建議

4.1 結論

長三角耕地呈北方高、南方低的格局，與長三角地形分布格局基本一致；主要是由耕地轉向建設用地，耕地利用動態度呈先增加后減少的變化趨勢。耕地減少形態主要分為炸裂式、帶狀式、星狀式。耕地規模經營減弱，碎片化程度不斷上升，但耕地形狀變得更加簡單、規整。耕地重心由長江以南轉向長江以北，從南京建鄴區轉向浦口區，未來可能會由江蘇省南京市轉向安徽省滁州市。

4.2 建議

長三角耕地數量不斷減少，今后應更加嚴格征地政策，平衡好城市擴張、經濟發展與糧食安全的關系。加強耕地的規模經營，使耕地更加集約利用，減少經營成本。未來的耕地重心可能會由江蘇省南京市轉向安徽省滁州市，耕地重心向西北方移動，安徽省和江蘇省作為糧食大省，加強耕地保護尤其重要。