北京平原區糧菜林三地類時空轉變特征

奉 婷，張鳳榮，聶 鑫，謝 臻，汪 晗

?

北京平原區糧菜林三地類時空轉變特征

奉 婷1,2，張鳳榮2※，聶 鑫1，謝 臻2，汪 晗1

（1. 廣西大學公共管理學院，南寧 530004； 2. 中國農業大學資源與環境學院，北京 100193）

該文利用遙感與地理信息系統技術研究北京平原區農用地內部糧田、菜地和林地三大類農地利用類型的變化。結果表明：2004年至2014年這10 a間，糧田轉變為林地與菜地是北京平原區農地利用變化主流。糧田轉變為林地面積是969.97km2，其土地利用變化重要性指數高達72.11%，在各類農地利用變化類型中居首。而糧田轉為菜地則主要是糧田轉變為設施菜地，10 a間糧田轉變為設施菜地面積是163.53km2，其土地利用變化重要性指數為12.16%，在各類農地利用變化類型中居次；同時，10 a間北京平原區農地利用變化存在空間分異特征。主城近郊區與遠郊平原區農地利用變化主要是糧田轉為林地，兩大分區中該變化類型土地利用變化重要性指數分別為84.18%、76.21%，明顯高于其他農地利用變化類型土地利用變化重要性指數。在城市新區，糧田則分別轉為林地與設施菜地；10 a間北京平原區農地利用空間分布均具有圈層特征，由近郊向遠郊，農地內部結構都呈現林地為主類型，到糧田、菜地與林地混合利用類型，再到糧田與林地組合利用類型的空間格局。而10 a間隨著主城區的外擴，上述農地利用圈層特征出現了由近郊向遠郊外移、且林地與蔬菜面積占比增加的現象。研究認為：農地內部糧、菜、林3地類的轉變和農地利用圈層外移是農戶追求經濟效益最大化與政府為改善生態環境所致；而改善生態不應以犧牲糧田為代價，應充分利用糧田的生態系統服務功能以實現生態綠化目的。

土地利用；遙感；地理信息系統；農地利用變化；時空轉變特征；北京平原區

0 引 言

自從1993年，“國際地圈與生物圈計劃”（international geosphere-biosphere program）和“全球變化人文因素計劃”（international human dimensions programme on global environmental change）兩大國際組織共同啟動“土地利用/土地覆被變化（land-use and land-cover change，LUCC）科學研究計劃”并將其作為全球變化研究核心內容后，土地利用/土地覆被變化成為了全球氣候變化與全球環境變化研究重要內容[1-4]。目前國際上以人類-環境耦合系統為核心的土地利用/覆被動態過程監測與模擬成為該研究領域熱點問題[5-6]。自20世紀90年代中后期受國際科學研究計劃影響，國內LUCC研究也活躍開來，主要包括土地利用變化時空特征，土地利用變化生態環境效應以及土地利用變化驅動機制研究內容等[7-11]。相較于宏觀尺度研究，在國家提出糧食安全保障背景影響下，國內中微觀尺度關于農地利用變化相關研究多圍繞耕地保護與利用展開，涉及到耕地質量、景觀格局以及集約度等方面變化情況[12-15]。另外，近些年來中國城鎮化進程高速發展，北京、上海等大都市區均面臨耕地資源保護難的問題，在這些區域農地利用諸多研究針對基本農田劃定、高標準基本農田建設以及耕地多功能性等內容展開了探討[16-18]。這些探索不僅基于中國國情豐富了土地利用/覆被變化研究內涵，同時還對農業發展與農地可持續利用具有理論指導價值。但目前農用地、耕地變化方面研究多從土地利用類型角度出發，鮮有將土地利用類型與其生產經營方式作為整體展開分析[19]，本研究在土地利用類型基礎上結合生產經營方式對農地利用變化進行探索，以期尋求突破。

受城市經濟發展影響，加之特殊的政治與文化地位，近10 a首都北京的土地利用/覆被發生顯著變化[20-21]；特別是平原區土地利用變化劇烈。除了同其他城市一樣存在著城市化造成建設用地的不斷擴張蠶食耕地外，另一個值得注意的是，北京市為建設宜居城市而進行的平原造林工程[22]。平原區土地利用率高，幾乎沒有未利用地；平原造林綠化也勢必與糧爭地。北京市龐大的人口對蔬菜需求巨大[23-24]，在城市擴張、平原綠化造林過程中，菜地又發生了怎樣的變化。因此，研究北京平原區的農地內部糧、菜、林三地類時空轉變特征，這對于揭示平原區大城市城市化和生態建設過程中的農地變化規律具有科學意義，以期為平原大城市土地利用優化與可持續發展提供參考。

1 研究區概況

北京市主要地貌類型有山地、丘陵和臺地、平原三種類型(圖1)。由于地形平緩、土壤肥沃以及水源豐富等自然稟賦條件，平原區已成為全市各類生產生活活動最主要場所。北京作為環渤海經濟圈的核心區，工業化與城市化程度很高。2013年全市生產總值19 500.6億元，一、二、三產業產值比為0.83∶22.32∶76.85。全市常住人口為2 114.8萬人。據《北京統計年鑒2014》數據顯示：2013年北京市農作物播種面積中糧食與蔬菜播種面積占比分別為65.70%、28.51%，近10 a間糧食與蔬菜一直為北京市最主要農作物，對耕地集中分布區——平原區而言更是如此。同時，據近10 a《北京農村年鑒》數據顯示：北京設施農業以蔬菜、瓜果種植為主，苗圃種植與牲畜養殖較少。另外，同期北京平原區造林工程實施對農地利用影響也不可小覷。受城鎮化和生態建設的影響，北京平原區目前不僅面臨耕地保護巨大壓力，其農地內部結構也必然發生變化。

2 數據來源與研究方法

2.1 數據來源

本研究的平原區是指高程100 m以下及坡度3°～5°以下的，北京市內包括了山前臺地和沖積平原區內的區域。但考慮到空間的連片性，不包括延慶盆地的平原區。土地利用變化采用了2004年Landsat TM 5影像和2014年Landsat OLI影像。根據北京市平原區主要農作物類型與農地利用特點，遙感數據分辨率與本研究需要，通過外業調查并結合農地利用方式將區域農地利用劃分為糧田、設施菜地、露天菜地、林地4大類型（表1）。其中，考慮到果園多分布在山區，平原區園地類型劃分對本研究意義不大，并且平原區造林在遙感影像上紋理特征與園地相近、不易區分，故將果園與有林地、未成林地、苗圃等統一劃分為林地；依據農作物類型不同劃分為糧田與菜地；又考慮到設施菜地與露天菜地在集約度、經濟收益和生態效應上的差異，將菜地又劃分為設施菜地和露天菜地。另外，據《2012年北京市土地利用變更調查》數據得到：平原區草地面積僅占該區域總面積0.02%，且近10 a內變化不顯著，故不將其納入本研究范疇。最終，本文將基于糧田、設施菜地、露天菜地、林地4大類型對北京平原區糧、菜、林3大主要農地利用類型時空轉變特征展開研究與討論。借助ENVI5.1與ArcGIS10.2軟件平臺，運用人機交互式解譯方法，得到北京平原區2004年與2014年2期的農地利用圖斑數據庫和屬性數據庫。

表1 北京平原區農地利用分類體系

在實地采集GPS樣本點基礎上結合高分辨率遙感影像對分類結果進行精度檢驗：其中2014年遙感影像，從采集GPS樣本點中抽取具有代表性的土地覆被驗證點125個，并在Google Earth軟件同時期高精度影像上補充采集175個，共形成300個驗證點進行檢驗；而2004年影像驗證點則主要利用Google Earth高分辨率影像，再結合北京市國土局提供的2004年土地利用數據，同樣挑選300個驗證點進行檢驗。分類精度檢驗結果是：2004年與2014年分別為87.23%和85.11%，滿足遙感土地利用解譯精度要求。

2.2 研究方法

土地利用變化特征包括土地利用數量規模、內部結構與空間格局變化等方面，參照姜廣輝等[25]研究將平原區劃分為主城區（東城和西城）、主城近郊區（朝陽、豐臺、石景山與海淀）、城市新區（通州、順義、大興全域和昌平與房山的平原區）以及遠郊平原區（平谷、門頭溝、懷柔與密云各區的平原部分），對各分區農地利用變化展開測度。

2.2.1 土地利用變化幅度

應用公式（1）研究區域土地利用總體變化態勢[26]

式中U、U分別為研究初與研究末期某一種土地利用類型數量，是土地利用變化幅度。

2.2.2 土地利用變化重要性指數

通過測算土地利用變化重要性指數（C），得到土地利用變化主要類型。公式如下[27]

（3）

式中C是第種變化類型的土地利用變化重要性指數，取值在0～1；A指第類土地變化面積，km2；是區域各類土地變化面積和，km2。當C越大，則說明第類土地變化為主導變化類型。對于土地變化種類分析，將C值從大到小降序排列，統計C值累計和大于80%的土地變化類型，由此得出區域主要土地變化類型；采用同樣的方法統計C值大于1%以及其累計和大于90%的農地變化類型作為農地內部調整主要類型。

2.2.3 洛倫茲曲線

洛倫茲曲線可用于分析空間單元的分布，也可用于對比2個空間單元分布狀況，是研究離散區域分布的重要方法。土地利用變化研究中可利用此方法分析區域土地利用結構集中化趨勢：當各土地利用類型面積越接近時，曲線越接近均勻分布線；而反之當區域土地利用集中于某一種或幾種類型時，曲線則越偏離均勻分布曲線[28-29]。本研究將用此曲線分析北京不同分區各土地利用類型空間分布集中和離散狀態。

2.2.4 組合類型

威弗—托馬斯（Weaver-Thomas）組合系數法能夠反映農地不同利用方式的結構特征與其組合類型[30]。該方法原理是利用土地假設分布（假設相對面積百分比）來逐步逼近實際分布（實際相對面積百分比），由此得到最接近實際分布土地利用組合類型。該方法思路是：先將各種農地利用類型面積比例由大到小排列。然后假設土地只為一種類型，則將這類型假設分布為100%，而其他類型假設分布為0；如果是分配給2種類型，則這2種類型假設分布各為50%，其他類型假設分布為0；以此類推，如果區域土地均勻分配給糧田、設施菜地、露天菜地和林地4種農地利用類型，則分布假設為25%。在此基礎上，測算每種假設分布與實際分布之差的平方和，即是組合系數。最后選擇組合系數最小的假設分布類型作為該區域農地利用組合類型。

3 結果與分析

3.1 糧田、設施菜地、露天菜地與林地數量變化

2004年與2014年北京平原區農地利用類型分布情況如圖2所示。借助ArcGIS平臺空間分析與統計功能，依據公式（1）對研究區2期農地利用圖形數據進行面積統計與分析，得到不同時期北京平原區各類農地利用面積及變化幅度（表2）。結果顯示：林地、設施菜地與露天菜地用地面積都有不同程度增加，林地凈增面積高達924.50 km2，由2004年864.34 km2增加至2014年的17 88.83 km2，增幅為106.96%；設施菜地凈增面積居次，由2004年119.56 km2增至250.17 km2，10 a間共增加了130.61 km2，增幅為109.24%；露天菜地由2004年23.16 km2增至2014年的49.37 km2，增幅為113.17%；相較于前三者類型的增加，糧田減少態勢明顯，平原區10 a間糧田共減少了1 652.89 km2，由2004年2 790.83 km2減至1 137.94 km2，減幅為59.23%。

3.2 糧田、設施菜地、露天菜地、林地空間轉化

3.2.1 平原區農地主要變化類型

依據公式（2）和（3）得到不同變化類型的土地利用變化重要性指數后，將其從大到小降序排列；將C值大于1%以及其累積之和大于90%的農地變化類型提取出來，篩選得到研究區2004―2014年農地利用內部結構調整情況（表2）。

分析得到，研究區農地變化面積為1 345.22 km2，變化面積占研究區總面積21.48%。在發生變化的12類農地內部調整類型中，糧田轉為林地與糧田轉為設施菜地這2種類型土地變化C值累積之和已高達84.28%，而其中糧田轉為林地變化類型C值高達72.11%，可見該區域糧田轉變為林地是農地內部結構調整主要趨勢。另一方面，糧田轉為設施菜地變化類型C值為12.16%，糧食作物種植呈下降趨勢，蔬菜比例上升；并且在各類主要農地利用內部調整過程中，蔬菜種植多以日光溫室與大棚等生產方式為主。

表2 2004—2014年農地利用結構調整情況

注：C為土地利用變化重要性指數。

Note: Cis land use change intensity index.

3.2.2 平原區農地變化類型空間分異

土地利用變化會因區位差異而出現空間分異。通過主城近郊區、城市新區與遠郊平原區2004年、2014年農地利用類型分布圖空間疊加，根據式（2）與（3）分析得到各分區土地利用變化重要性指數，提取C值累積之和大于90%的類型，最終得到不同分區2004－2014年農地內部結構調整空間分異情況（圖3）。

結果顯示：主城近郊區與遠郊平原區農地利用變化主要是糧田轉為林地，兩大分區中該變化類型土地利用變化重要性指數分別為84.18%、76.21%，明顯高于其他農地利用變化類型土地利用變化重要性指數。農地內部結構調整中，主城近郊區、城市新區與遠郊平原區呈現出變化類型由少變多再變少趨勢：其中三大分區均表現為糧田轉變為林地類型為首的特征。有所不同的是主城近郊區不僅大部分糧田轉為林地，還有部分設施菜地也轉為林地。由此可見，緊鄰主城區的近郊區農地結構調整趨勢是由集約度較高的農地利用類型向集約度更低的類型轉變。而遠郊平原區還有部分林地轉為了糧田、糧田又轉為設施菜地，可見該區域農地結構調整態勢為集約度較低的農地利用類型向集約度更高類型轉變。與這兩大分區有所不同的是，城市新區農地利用內部結構調整類型則更多樣化與均衡化。除了糧田、露天菜地轉為林地之外，還有很大一部分林地轉為糧田，少量糧田、林地轉為設施菜地，該區域整體上農地利用內部結構調整最為頻繁。作為北京過去的“糧倉”，大興、通州與順義等城市新區如今仍是農業生產活動最活躍區域。

3.3 糧田、設施菜地、露天菜地、林地集散狀態

根據洛倫茲曲線繪制方法，將各分區內糧田、設施菜地、露天菜地、林地占總面積百分比從大到小順序排列，然后計算其累積百分比，據此繪制北京平原區2004與2014年平原區及各分區洛倫茲曲線（圖4）。從洛倫茲曲線能夠清晰了解糧田、設施菜地、露天菜地、林地的集中化程度與空間差異。

a. 主城近郊區

a. Suburban area

b. 城市新區

b. New urban districts

從2期農地利用洛倫茲曲線變化來看，2014年平原區洛倫茲曲線比2004年曲線更靠近均勻曲線，整體而言，平原區農地利用發生著由糧田為主向多種農地利用類型共存變化。2004年糧田、林地、設施菜地與露天菜地分別占研究區農地總面積比例為73.48%、22.76%、3.15%與0.61%；2014年各比例分別變化為35.27%、55.45%、7.75%與1.53%。數據顯示：平原區農地利用2004年主要以糧食作物種植為主，因此農地集中化程度較高；發展至今農地利用呈現出林地與糧田為主，設施菜地增長態勢，農地利用多樣化水平提高。另外，不同分區農地利用集散狀態也有所差異：其中主城近郊區與城市新區變化最為明顯，前者在2004年是各分區中農地利用集中化水平最低區域，并且明顯低于平原區平均水平，發展至2014年其農地利用洛倫茲曲線最偏離平均水平，成為平原區農地利用集中化程度最高區域；后者則由2004年農地利用集中化水平最高區域變成2014年農地利用離散度最高區域。數據顯示：主城近郊區2004年糧田、林地占該區域農地總面積比例分別為54.14%與39.95%，到2014年該區域林地占農地總面積比例高達86.33%，而糧田比例降低至10.22%，林地躍升為主城近郊區最主要農地利用類型是該區域農地利用集中度劇增主因。而城市新區糧田、林地與設施菜地占該區域農地總面積比例由2004年78.94%、17.12%與3.28%分別為2014年的38.64%、50.14%與9.23%，可見該區域農地利用內部結構調整最活躍，大量糧田不僅轉為林地，同時部分糧田也轉為設施菜地，由過去糧田為主農地利用結構向林地、糧田與設施菜地多元農地利用均衡發展過渡。相較于前兩者，雖然遠郊平原區10 a間農地利用集中化程度變化并不顯著，但從統計數據看，這是糧田與林地占區域農地總面積比例互相置換的結果：2004年該區域糧田、林地占其農地總面積比例分別是61.93%與36.34%，而2014年糧田、林地比例分別為33.42%與62.16%；10 a間出現這種變化特征一方面源于該區域處于北京山前平原過渡帶，特有的自然地理條件決定其在研究初期就已形成“糧食—林地”農業生產結構，另一方面也是源于平原造林政策實施的結果。

3.4 糧田、設施菜地、露天菜地、林地組合類型變化

采用威弗—托馬斯組合系數法，以鄉鎮（街道）為單元對2004與2014年農地利用組合類型空間分布情況進行分析。

從圖5結果來看，雖然10 a間農地利用組合類型更多樣化，同時各種組合類型也發生了不同程度的空間變化，但2期農地利用組合類型空間分布都具有圈層特征：2004年，臨近主城區的近郊區街道、鄉鎮多以林地種植為主，呈現出環繞著中心城區的林地圈層；隨著中心城區向主城近郊區外圍延伸，開始出現糧田—設施菜地—林地混合型圈層，但其空間分布并不廣泛，以豐臺、朝陽與海淀居多；再向遠郊延伸，城市新區多為糧田—林地與糧田為主的組合類型，并且這些區域中，即使是糧田—林地組合類型也多以糧田為主、林地為輔；2004年糧田—林地與糧田為主的組合類型空間分布最廣泛，城市新區主要為這兩大組合類型，其圈層特征也最為明顯。最后延伸到遠郊北部鄰近山地區域開始出現林地為主的組合類型，多分布在平谷、懷柔與密云平原區北部。整體而言，2004年農地利用多以糧食作物種植為主，至于蔬菜種植、特別是設施菜地種植比例較低；由于這種農地利用類型的“單一化”，令其農地利用圈層特性在空間上并非是勻稱分布。

2014年與2004年農地利用空間分布特征相比，同具有圈層性：靠近主城區依舊分布的是林地為主的農地利用圈層；并且該圈層與2004年相較在空間上明顯發生擴張，幾乎已覆蓋整個主城近郊區。另外，緊挨著林地圈層的依然是糧田—設施菜地—林地混合型圈層，并且多分布在城市新區，以大興北部為主，部分位于通州中部與北部、房山東部，以及少量零星分布在順義與昌平；同時，遠郊平原區中的平谷與密云也開始出現該組合類型。而這一圈層隨著林地為主圈層在空間上的擴張，與2004年相比空間上也有所向遠郊遷移，并且空間分布較過去更加廣泛。再繼續向遠郊延伸，出現的依舊是糧田—林地圈層，主要分布在城市新區外圍以及遠郊平原區；依據糧田與林地占研究單元農地總面積比例大小差異發現，該圈層與2004年相比最大不同在于，2014年以糧田為主類型已不再有，而2004年糧田與林地組合類型中更傾向于以糧田為主，而如今該組合類型中則出現林地為主趨勢，并且林地為主、糧田為輔的組合類型在空間分布上比糧田為主、林地為輔或者糧田、林地均分的組合類型更臨近城區與近郊區，糧田為主、林地為輔或者糧田、林地均分組合類型則幾乎分布在房山西南部、大興與通州南部、平谷北部及密云等平原最遠郊區域。

4 討 論

盡管10 a間農地利用各圈層由中心城區向郊區在不斷外遷，但農地利用空間格局整體上還是保持著近郊蔬菜種植為主、遠郊糧田與林地種植為主特征；而這種隨著中心城區向農村腹地遷移，由集約度較高的設施菜地生產方式逐步向集約度較低的糧食、林地種植模式轉變的空間分異特征，即城鎮化發展下近城區地租“水漲船高”，為平衡農業生產過程中成本收益關系，農戶會根據地租變化適時調整種植作物與生產模式。農地利用以蔬菜—糧田—林地這種圈層結構關系并向遠郊區域“外移”的過程是對城鎮建設擴張下地租上漲壓力的響應。

然而農業區位論無法解釋的是：緊鄰主城區的近郊區農地利用中是以林地為主而并非菜地，這是因為大都市區政府為改善生態環境進行平原造林的結果。但導致平原區非常珍貴的耕地資源大量減少，使得國家耕地保護目標不僅受到建設占用的壓力，同時還要遭受生態建設的擠占。大都市區即要實現耕地保護目標又實現生態建設目標，則應考慮充分利用耕地的生產與生態雙重功能，而由于菜地的水肥藥投入高，顯然耕地種植糧食作物能夠更好發揮生態功能[31]。

5 結 論

1）城鎮化與平原綠化工程推動下，北京平原區農地利用內部各類型間都發生大面積流轉。在該區域農地內部各類型變化中，糧田、設施菜地與林地是變化顯著的農地利用類型，大量糧田轉變為林地與設施菜地，糧田轉變為林地面積是969.97 km2，其土地利用變化重要性指數高達72.11%，在各類農地利用變化類型中居首。而糧田轉為菜地則主要是糧田轉變為設施菜地，10 a間糧田轉變為設施菜地面積是163.53 km2，其土地利用變化重要性指數為12.16%。農地利用內部結構調整呈現空間分異特征：距離中心城區最近的主城近郊區與距離中心城區最遠的遠郊平原區均表現出糧田轉變林地為主；地域上居于前面兩大分區之間的城市新區農地利用內部變化類型最多，除了糧田與林地之間互相轉變外，露地與設施菜地跟糧田、林地間的轉變也多發生在該區域內。由此可見，城市新區中大興、順義與通州不僅在過去是北京糧食主產區，至今依然是農業產業結構調整與生產活動最活躍區域。從農地利用變化態勢來看，北京平原造林工程是導致該區域糧田轉為林地的主因；與此同時，糧田轉變為菜地，而菜地中設施菜地面積增幅又遠大于露天菜地。

2）農地利用類型空間集散狀態變化表現為：主城近郊區在造林工程實施作用下大量糧田轉變為林地，林地成為該區域最主要農地類型，農地利用集中度大幅升高。與此同時，城市新區也受到平原造林影響，大量糧田轉為林地；另外還有部分糧田轉變為設施菜地，農地利用由低集約種植向高集約種植轉型，導致該區域農地利用類型在空間分布上更趨多元化與均衡化。相對前面兩者，遠郊平原區農地利用類型依舊主要為糧田與林地，而這是10 a間糧田與林地面積在該區所占空間比例互相置換下所出現的結果。

3）近10 a北京平原區農地利用組合類型空間分布都具有圈層特征，由中心城區向遠郊延伸圈層依次為林地圈層、糧田—設施菜地—林地混合圈層以及糧田—林地組合圈層。兩期有所不同的是：緊鄰中心城區林地圈層發展至今空間分布更廣，幾乎已覆蓋整個主城近郊區；而臨近該圈層所出現設施菜地—糧田—林地混合圈層其所涉空間范圍也在變廣，多分布在大興與通州中北部，呈現出由近郊區向遠郊遷移態勢。而距離中心城區最遠的糧田—林地圈層2004－2014年內有所變化的是，過去該圈層糧田為主，林地為輔，而如今該圈層兩者比例則發生置換，呈現林地為主、糧田為輔特征。兩期農地利用圈層結構表明平原區由中心城區向遠郊輻射，農地利用集約度具有“低—高—低”變化態勢這一空間分異特征。

[1] Mooney H A, Duraiappah A, Larigauderie A. Evolution of natural and social science interactions in global change research programs.[J]. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(Supp.1): 3665－3672.

[2] Sterling S M, Ducharne A, Polcher J. The impact of global land-cover change on the terrestrial water cycle[J]. Nature Climate Change, 2013, 3(4): 385－390.

[3] 李秀彬. 全球環境變化研究的核心領域——土地利用/土地覆被變化的國際研究動向[J]. 地理學報，1996，51(6)：553－558.

Li Xiubin. A review of international researches on land use/land cover change[J]. Acta Geographica Sinica, 1996, 51(6): 553—558. (in Chinese with English abstract)

[4] Ojima D, Moran E, Mcconnell W, et al. Global land project: Science plan and implementation strategy[J]. Environmental Policy Collection, 2005.

[5] IGBP Secretariat. GLP(2005) Science plan and Implementation Strategy[R]. IGBP Report No.53/IHDP Report No.19, 2005, Stockholm: 64pp.

[6] Turner II 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.

[7] 劉紀遠，匡文慧，張增祥，等. 20世紀80年代末以來中國土地利用變化的基本特征與空間格局[J]. 地理學報，2014，69(1)：3－14.

Liu Jiyuan, Kuang Wenhui, Zhang Zengxiang, et al. Spatiotemporal characteristics, patterns and causes of land use changes in China since the late 1980s[J]. Journal of Geographical Sciences, 2014, 69(1): 3－14. (in Chinese with English abstract)

[8] 劉成武，李秀彬. 1980－2002年中國農地利用變化的時序特征[J]. 農業工程學報，2006，22(4)：194－198.

Liu Chengwu, Li Xiubin. Time-sequence characteristics of the annual changes of the agricultural land use in China during the period 1980－2002[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2006, 22(4): 194－198. (in Chinese with English abstract)

[9] 楊依天，鄭度，張雪芹，等. 1980－2010年和田綠洲土地利用變化空間耦合及其環境效應[J]. 地理學報，2013，68(6)：813－824.

Yang Yitian, Zheng Du, Zhang Xueqin, et al. The spatial coupling of land use changes and its environmental effects on Hotan oasis during 1980－2010[J]. Acta Geographica Sinica, 2013, 68(6): 813—824. (in Chinese with English abstract)

[10] 彭建，王仰麟，張源，等. 滇西北生態脆弱區土地利用變化及其生態效應——以云南省永勝縣為例[J]. 地理學報，2004，59(4)：629—638.

Peng Jian, Wang Yanglin, Zhang Yuan, et al. Land use change and its ecological effect in the ecotone of northwest of Yunnan province, China: A case study of Yongsheng county[J]. Acta Geographica Sinica, 2004, 59(4): 629－638. (in Chinese with English abstract)

[11] 劉紀遠，張增祥，徐新良，等．21世紀初中國土地利用變化的空間格局與驅動力分析[J]. 地理學報，2009，64(12)：1411－1420.

Liu Jiyuan, Zhang Zengxiang, Xu Xinliang, et al. Spatial patterns and driving forces of land use change in China in the early 21st century[J]. Acta Geographica Sinica, 2009, 64(12): 1411－1420. (in Chinese with English abstract)

[12] 孔祥斌，劉靈偉，秦靜. 基于農戶土地利用行為的北京大興區耕地質量評價[J]. 地理學報，2008，63(8)：856－868.

Kong Xiangbin, Liu Lingwei, Qin Jing. Arable land evaluation based on the household land use behavior in Daxing district of Beijing[J]. Acta Geographica Sinica, 2008, 63(8): 856－868. (in Chinese with English abstract)

[13] 姜廣輝，趙婷婷，段增強，等．北京山區耕地質量變化及未來趨勢模擬[J]. 農業工程學報，2010，26(10)：304－311.

Jiang Guanghui, Zhao Tingting, Duan Zengqiang, et al. Cultivated land quality change and its future trend modeling in Beijing mountainous area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(10): 304－311. (in Chinese with English abstract)

[14] 岳德鵬，王計平，劉永兵，等. 京郊西北地區農地利用與景觀格局時空特征[J]. 農業工程學報，2008，24(6)：89－95.

Yue Depeng, Wang Jiping, Liu Yongbing, et al. Spatial and temporal characteristics of agricultural land-use and landscape pattern changes in northwest of Beijing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(6): 89－95. (in Chinese with English abstract)

[15] 崔麗，許月卿. 河北省農用地利用集約度時空變異分析[J]. 地理科學進展，2007，26(2)：116－125.

Cui Li, Xu Yueqing. Spatial-temporal variance of agricultural land use intensity in Hebei Province[J]. Progress in Geography, 2007, 26(2): 116－125. (in Chinese with English abstract)

[16] 奉婷，張鳳榮，李燦，等. 基于耕地質量綜合評價的縣域基本農田空間布局[J]. 農業工程學報，2014，30(1)：200－210.

Feng Ting, Zhang Fengrong, Li Can, et al. Spatial distribution of prime farmland based on cultivated land quality comprehensive evaluation at county scale[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(1): 200－210. (in Chinese with English abstract)

[17] 王新盼，姜廣輝，張瑞娟，等. 高標準基本農田建設區域劃定方法[J]. 農業工程學報，2013，29(10)：241－250.

Wang Xinpan, Jiang Guanghui, Zhang Ruijuan, et al. Zoning approach of suitable areas for high quality capital farmland construction[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(10): 241—250. (in Chinese with English abstract)

[18] 楊雪，談明洪. 北京市耕地功能空間差異及其演變[J]. 地理研究，2014，33(6)：1106－1118.

Yang Xue, Tan Minghong. Spatial differences and evolution of arable land functions in Beijing[J]. Geographical Research, 2014, 33(6): 1106－1118. (in Chinese with English abstract)

[19] 張鳳榮，孔祥斌，徐艷. 開展農地利用方式變化規律研究探討土地可持續利用模式[J]. 中國農業科技導報，2002，4(3)：18－22.

Zhang Fengrong, Kong Xiangbin, Xu Yan. Studied on the change law of farmland management patterns to develop sustainable land management patterns[J]. Review of China Agricultural Science & Technology, 2002, 4(3): 18－22. (in Chinese with English abstract)

[20] 孫強，蔡運龍，王樂. 北京耕地流失的時空特征與驅動機制[J]. 資源科學，2007，29(4)：158－163.

Sun Qiang, Cai Yunlong, Wang Le. Temporal-spatial processes of farmland conversion and its driving mechanisms in Beijing[J]. Resources Science, 2007, 29(4): 158－163. (in Chinese with English abstract)

[21] 趙媛媛，何春陽，龔立萍，等. 北京城市擴展過程中耕地自然生產功能損失研究[J]. 中國土地科學，2009，23(7)：71－78.

Zhao Yuanyuan, He Chunyang, Gong Liping, et al. Production capacity loss of farmland during urbanization in Beijing, China[J]. China Land Science, 2009, 23(7): 71－78. (in Chinese with English abstract)

[22] 造林105萬畝植樹5 400萬株北京市平原百萬畝造林工程全面完成[J]. 綠化與生活，2016(1)：4－5.

[23] 趙霞，穆月英，潘鳳杰，等. 北京市自產蔬菜供需平衡分析——基于批發市場層面的初步測算[J]. 中國蔬菜，2011(21)：12－17.

[24] 北京市農林科學院蔬菜研究中心. 北京蔬菜產業發展基本特征和重大需求[J]. 蔬菜，2010(7)：1－3.

[25] 姜廣輝，付晶，譚雪晶，等. 北京國土空間結構與未來空間秩序研究——基于主體功能區劃框架[J]. 中國人口·資源與環境，2011，21(1)：20－27.

Jiang Guanghui, Fu Jing, Tan Xuejing, et al. Spatial structure and its future order of national land-use according to major function oriented zoning in Beijing[J]. China Population Resources and Environment, 2011, 21(1): 20－27. (in Chinese with English abstract)

[26] 韓會然，楊成鳳，宋金平. 北京市土地利用變化特征及驅動機制[J]. 經濟地理，2015，35(5)：148－154，197.

Han Huiran, Yang Chengfeng, Song Jinping. The spatial- temporal characteristic of land use change in Beijing and its driving mechanism[J]. Economic Geography, 2015, 35(5): 148－154, 197. (in Chinese with English abstract)

[27] 羅婭，楊勝天，劉曉燕，等. 黃河河口鎮—潼關區間1998-2010年土地利用變化特征[J]. 地理學報，2014，69(1)：42－53.

Luo Ya, Yang Shengtian, Liu Xiaoyan, et al. Land use change in the reach from Hekouzhen to Tongguan of the Yellow River during 1998－2010[J]. Acta Geographica Sinica, 2014, 69(1): 42—53. (in Chinese with English abstract)

[28] 扈傳榮，姜棟，唐旭，等. 基于洛倫茲曲線的全國城市土地利用現狀抽樣分析[J]. 中國土地科學，2009，23(12)：44－50. Hu Chuanrong, Jiang Dong, Tang Xu, et al. Sampling analysis of national urban land use status based on Lorentz Curves[J]. China Land Science, 2009, 23(12): 44－50. (in Chinese with English abstract)

[29] 陳其春，呂成文，李壁成，等. 縣級尺度土地利用結構特征定量分析[J]. 農業工程學報，2009，25(1)：223－231.

Chen Qichun, Lü Chengwen, Li Bicheng, et al. Quantitative analysis of land use structure characteristics at county scale[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2009, 25(1): 223－231. (in Chinese with English abstract)

[30] 李燦，張鳳榮，姜廣輝，等. 京郊衛星城區域農村居民點土地利用特征分析[J]. 農業工程學報，2013，29(19)：233—243.

Li Can, Zhang Fengrong, Jiang Guanghui, et al. Analysis on utilization characteristics of rural residential land use in satellite city area of Beijing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(19): 233－243. (in Chinese with English abstract)

[31] 張鳳榮，趙華甫，黃大全，等. 基于宜居城市建設的北京市農業產業空間布局[J]. 資源科學，2008，30(2)：162—168.

Zhang Fengrong, Zhao Huafu, Huang Daquan, et al. Spatial planning for agricultural production for making the Beijing city as a good living place[J]. Resources Science, 2008, 30(2): 162－168. (in Chinese with English abstract)

Spatial-temporal change characteristics of cropland, vegetable land and forest land in Beijing plain region

Feng Ting1,2, Zhang Fengrong2※, Nie Xin1, Xie Zhen2, Wang Han1

(1.530004;2.,100193,)

This research used the techniques of the remote sensing (RS) and geographic information system (GIS) to obtain the data of agricultural land use classifications by analyzing and interpreting the multi-spectrum remote sensing images, which were extracted from Landsat TM image in 2004 and Landsat OLI image in 2014. According to the situation of agricultural land use in Beijing plain, the agricultural land use were classified as the grain crop fields, vegetables cultivated in greenhouse, vegetables cultivated in the open fields and forest land. Then it analyzed spatial-temporal change characteristics of the agricultural land use in Beijing plain over the past decade (2004 to 2014) by combination of the quantitative geography model and geographic information system (GIS). The results showed that the grain crop fields were changed into forest land and greenhouse vegetable fields, which was the major trend. The area of grain crop fields changed into forest land was 969.97 km2whose land use change intensity index was 72.11%, at the top of all kinds of land use change intensity indexes. The area of grain crop fields changed into greenhouse vegetable fields was 163.53 km2during the past 10 years, whose land use change intensity index was 12.16%, and ranked at the second among all kinds of land use change intensity indexes. Meanwhile the agricultural land use change had spatial characteristic in Beijing plain over the past 10 years. The leading type of agricultural land use changes was that grain crop fields turned into forest land in the suburban area and the distant suburb plain of Beijing, where land use change intensity indexes of grain crop fields into forest land were 84.18% and 76.21%, respectively and higher than the other agricultural land use change intensity indexes obviously. The overwhelming urbanization was encroaching upon a lot of grain crop fields, which leaded the agricultural land use to concentrate on the forest land in the suburban districts. At the same time grain crop fields changed separately to forest land and greenhouse vegetable fields in the new urban area of Beijing, which prompted to the equalization of agricultural land use structure. The research area kept the specific spatial pattern of agricultural land use that the forest land was the leading classification in close proximity to the city center, the mixed classifications including grain crop fields, greenhouse vegetable fields and forest land were main structure in suburban area of Beijing plain and the combination type consisted of grain crop fields and forest land mainly appeared in the distant suburbs of Beijing plain. The spatial structure of the agricultural land use showed as a specific spatial pattern including 3 layer circles and had experienced some changes over the last decade. Firstly, the whole layer circles of agricultural land use had migrated outwards towards the suburbs with the rapid expansion of the central city over the time. Secondly, the spatial distributions of forest land in close proximity to city center and the vegetables cultivated in greenhouse in suburban areas of Beijing plain both had expanded outward due to the area of forest land and vegetables cultivated in greenhouse increased. Although grain crop fields and forest land were the main agricultural land use types in the distant suburbs of Beijing plain in 2004, whose proportion of grain crop fields and forest land had exchanged in 2014. Forest land had become the main type of agricultural land use in the distant suburb of Beijing plain instead of grain crop fields over the time. The study identified that the spatial structure and changes of 3 layer circles of the agricultural land use objectively reflect the different land rent, which drove the peasant household in different locations of suburbs to make different choices about the cultivation, operation and management of agricultural products. Specifically, the peasant households who lived in urban and suburban areas preferred vegetables to grain crops and preferred cultivating vegetables in greenhouse to cultivating vegetables in the open fields, because they hope to get more economic income by changing agricultural products and ways of operation and management. On the other hand, forest land which mainly provided ecological service function was closer to the city center than the grain crop fields which mainly provided production function, which meant the agricultural function had changed from production function mainly to productive and ecological functions. The research results provide a reference for the optimization of spatial distribution of agricultural land use in Beijing plain.

land use; remote sensing; geographic information system (GIS); agricultural land use change; temporal and spatial variation characteristic; Beijing plain region

10.11975/j.issn.1002-6819.2017.06.033

F301.24

A

1002-6819(2017)-06-0257-08

2016-07-13

2017-02-14

國家自然科學基金項目（71363005；71403063）和國家社會科學基金項目（13CGL109）。

奉婷，女，廣西桂林人，博士，講師，主要研究方向為土地利用規劃與可持續利用。廣西 廣西大學公共管理學院，530004。Email：dreamland_21@163.com

張鳳榮，男，河北滄縣人，博士，教授。主要研究方向為土地評價、規劃與可持續利用。北京 中國農業大學資源與環境學院，100193。Email：frzhang@cau.edu.cn