• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    中國耕地景觀細碎度時空變化特征及其影響因素

    2023-01-13 01:07:26徐曉凡
    農(nóng)業(yè)工程學報 2022年16期
    關(guān)鍵詞:縣域耕地斑塊

    王 學,徐曉凡,2

    中國耕地景觀細碎度時空變化特征及其影響因素

    王 學1,徐曉凡1,2

    (1. 中國科學院地理科學與資源研究所,中國科學院陸地表層格局與模擬重點實驗室,北京 100101;2. 中國科學院大學,北京 100049)

    耕地細碎化是中國耕地利用效率提高和農(nóng)業(yè)現(xiàn)代化轉(zhuǎn)型的關(guān)鍵制約因素。2014年以來中央政府試圖通過鼓勵土地經(jīng)營權(quán)流轉(zhuǎn)以實現(xiàn)耕地細碎度降低和適度規(guī)模經(jīng)營耕地的目標。然而,少有研究開展中國全域尺度的耕地細碎化趨勢分析,土地經(jīng)營權(quán)流轉(zhuǎn)政策的相關(guān)實施效果也尚不明晰。該研究基于長時序土地利用數(shù)據(jù),輔以景觀格局指數(shù)、趨勢分析和動態(tài)面板模型等方法,對1990—2020年中國耕地景觀細碎度的時空變化特征進行了識別,并剖析了土地流轉(zhuǎn)政策實施對耕地景觀細碎化的影響。結(jié)果表明,1)中國耕地景觀細碎度大致呈“平原、盆地低、山地、高原高”的分布格局,東北平原、黃淮海平原、長江中下游平原以及四川盆地、準噶爾和塔里木盆地等地區(qū)耕地景觀細碎度低,而二三級地勢階梯過渡地帶以及黃土高原、云貴高原和東南丘陵等地區(qū)耕地景觀相對細碎;2)1990—2020年中國耕地景觀細碎度大致呈“持續(xù)細碎化-細碎化改善-波動細碎化”變化特征,其中2004—2013年耕地景觀整體呈細碎化趨勢,超過70%的縣域耕地斑塊密度(Patch Density,PD)/形狀指數(shù)(Landscape Shape Index,LSI)顯著增加、聚集度指數(shù)(Aggregation Index,AI)顯著減小,集中分布在黃淮海平原、長江中下游平原、四川盆地以及兩廣地區(qū),但2014—2020年70%的省份耕地景觀細碎度有所改善,以兩廣地區(qū)表現(xiàn)最為典型;3)土地流轉(zhuǎn)政策顯著降低了耕地景觀細碎度,此外,區(qū)域土地利用變化、非糧化耕地利用、坡度以及灌溉條件等因素也是影響中國耕地景觀細碎度變化的重要因素。該研究從全域尺度剖析了中國耕地景觀細碎度的時空變化特征并揭示土地流轉(zhuǎn)政策等因素對其的影響,研究結(jié)果可為區(qū)域耕地景觀細碎度變化以及不同土地利用/生態(tài)系統(tǒng)類型變化分析提供新思路。

    土地利用;空間格局;動態(tài)模型;耕地細碎化;景觀格局指數(shù);趨勢分析;土地流轉(zhuǎn)

    0 引 言

    耕地資源是保障國家糧食安全的基石,同時也在農(nóng)業(yè)現(xiàn)代化建設和鄉(xiāng)村振興過程中擔負著重要使命[1-3]。然而,中國耕地細碎化問題突出,一方面影響了耕地規(guī)模利用效應的發(fā)揮,成為制約農(nóng)業(yè)生產(chǎn)效率提高和農(nóng)業(yè)現(xiàn)代化轉(zhuǎn)型的關(guān)鍵因素[4-5];另一方面,在農(nóng)業(yè)生產(chǎn)要素價格不斷上漲和農(nóng)業(yè)勞動力非農(nóng)轉(zhuǎn)移趨勢下,細碎化耕地的經(jīng)營效益低下,嚴重降低了農(nóng)民的生產(chǎn)積極性并引發(fā)了較為普遍的耕地粗放利用、甚至撂荒現(xiàn)象,給國家糧食安全造成了一定威脅[6-7]。

    為了降低耕地細碎度、提高農(nóng)業(yè)生產(chǎn)率和優(yōu)化土地資源配置,國家在2014年印發(fā)了《關(guān)于引導農(nóng)村土地經(jīng)營權(quán)有序流轉(zhuǎn)發(fā)展農(nóng)業(yè)適度規(guī)模經(jīng)營的意見》,鼓勵土地經(jīng)營權(quán)有序流轉(zhuǎn),推進耕地適度規(guī)模經(jīng)營。此后,國家及地方政府頒發(fā)了多項政策法規(guī)文件,以保障農(nóng)村土地經(jīng)營權(quán)流轉(zhuǎn)的順利實施。截至2020年,全國流轉(zhuǎn)耕地面積占承包經(jīng)營耕地總面積的比例已達34%[8]。然而,需要指出的是,土地經(jīng)營權(quán)流轉(zhuǎn)與耕地經(jīng)營規(guī)模之間并不一定是線性相關(guān)關(guān)系,低水平、零散化的土地流轉(zhuǎn)甚至可能會在一定程度上加劇當?shù)氐母丶毸榛潭萚9-11]。因此,土地經(jīng)營權(quán)流轉(zhuǎn)相關(guān)政策實施后,中國耕地細碎度的變化情況并不十分明晰。

    近年來,耕地細碎化研究多圍繞其內(nèi)涵的兩個方面,也即權(quán)屬細碎化和景觀細碎化展開[3]。其中,耕地權(quán)屬細碎化主要從生產(chǎn)經(jīng)營主體這一微觀視角出發(fā),強調(diào)細碎化是產(chǎn)權(quán)主體無限細分的過程,相關(guān)研究多借助農(nóng)戶調(diào)查數(shù)據(jù),進行不同規(guī)模經(jīng)營主體的耕地利用效率對比以及細碎化對耕地利用方式的影響評估[12-15];耕地景觀細碎化則側(cè)重宏觀視角,認為細碎化是耕地地塊在形態(tài)上不斷分割、變小,且空間上不相連的過程,相關(guān)研究多借助土地利用數(shù)據(jù)和景觀格局指數(shù),進行景觀細碎化特征分析及影響因素識別[16-19]??傮w而言,受研究對象和數(shù)據(jù)可獲取性的制約,已有研究多局限在區(qū)域或局地尺度,且多為某一時點的結(jié)果,中國全域長時序耕地細碎度變化研究并不多見[20]。

    鑒于此,本文聚焦中國全域耕地景觀細碎度的變化特征,利用長時間序列土地利用數(shù)據(jù),輔以景觀格局指數(shù)、趨勢分析和計量模型,對1990—2020年中國耕地景觀細碎度的時空變化特征進行分析;同時結(jié)合區(qū)域土地利用情況、農(nóng)業(yè)生產(chǎn)管理情況、耕地資源稟賦特征、區(qū)域發(fā)展水平等,剖析土地流轉(zhuǎn)政策對耕地景觀細碎化的影響并識別耕地景觀細碎化的其他影響因素。該項研究不僅能夠揭示土地經(jīng)營權(quán)流轉(zhuǎn)政策實施前后中國全域耕地景觀細碎度的變化情況,為耕地適度規(guī)模經(jīng)營以及優(yōu)化利用等相關(guān)政策制定提供科學參考,還可為區(qū)域耕地景觀細碎度變化以及不同土地利用/生態(tài)系統(tǒng)類型變化分析提供新思路和新方法。

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

    1.1 數(shù)據(jù)來源

    耕地分布數(shù)據(jù)是進行耕地景觀細碎度分析的關(guān)鍵基礎數(shù)據(jù)。本文所用的耕地分布數(shù)據(jù)來源于中國1990—2020年長時間序列土地利用數(shù)據(jù)集(China Land Cover Dataset,CLCD),空間分辨率為30 m。該數(shù)據(jù)集是基于Google Earth Engine平臺和30余萬張Landsat影像、利用從中國土地利用/覆被數(shù)據(jù)集[21]以及Google Earth影像中提取的訓練樣本和隨機森林分類法,獲取的1990—2020年全國逐年土地利用數(shù)據(jù)(共31期),劃分了耕地、林地、灌木林地、草地、水域、冰雪、裸地、不透水面和濕地共9種土地利用類型(每種土地利用類型不再細分)[22]。對比發(fā)現(xiàn),CLCD數(shù)據(jù)集的整體準確性優(yōu)于MCD12Q1、ESACCI_LC、FROM_GLC和GlobeLand30等數(shù)據(jù)集,總體精度達到79.31%,其中各年份耕地的總體精度大都超過80%。

    除CLCD數(shù)據(jù)集,進行耕地景觀細碎化影響因素分析時還需要使用氣象數(shù)據(jù)、DEM數(shù)據(jù)和統(tǒng)計數(shù)據(jù)。其中,氣象數(shù)據(jù)包括2000—2020年全國逐年1 km×1 km分辨率年均溫和年降水量數(shù)據(jù),DEM數(shù)據(jù)空間分辨率為30 m × 30 m,兩者均來源于中科院資源環(huán)境科學與數(shù)據(jù)中心;統(tǒng)計數(shù)據(jù)主要包括2008—2020年各省份(港澳臺數(shù)據(jù)暫缺;經(jīng)計算,西藏基本無耕地,文中忽略未統(tǒng)計)耕地灌溉面積、糧食作物和農(nóng)作物播種面積、生態(tài)退耕面積、農(nóng)村人口數(shù)量、城鎮(zhèn)人口和總?cè)丝?、第一產(chǎn)業(yè)增加值和地區(qū)生產(chǎn)總值,來源于歷年的《中國統(tǒng)計年鑒》《中國農(nóng)村統(tǒng)計年鑒》《中國人口與就業(yè)統(tǒng)計年鑒》等。

    1.2 研究方法

    1.2.1 景觀格局指數(shù)計算

    耕地地塊在形態(tài)上不斷分割變小、在空間上互不相連和形狀各異可被視為耕地景觀細碎化的主要表征。已有研究通常采用斑塊數(shù)量(Number of Patches,NP)、斑塊面積(Total Class Area,CA)、斑塊密度(Patch Density,PD)、最大斑塊指數(shù)(Largest Patch Index,LPI)、形狀指數(shù)(Landscape Shape Index,LSI)、聚集度指數(shù)(Aggregation Index,AI)、分離度指數(shù)(Landscape Division Index,DIVISION)、凝聚度指數(shù)(Patch Cohesion Index,COHESION)等多個景觀格局指數(shù)進行耕地景觀細碎度評估[19-20];也有研究為降低指數(shù)間的冗余性,利用多個指數(shù)構(gòu)建綜合性指標,對區(qū)域或局地耕地景觀細碎度變化特征進行分析[23-24]。本文從耕地斑塊的易耕性和易規(guī)模利用視角出發(fā),結(jié)合耕地景觀細碎化內(nèi)涵以及已有研究成果,選用了表征耕地斑塊大小、形狀規(guī)則程度和空間連通性的3個關(guān)鍵指數(shù),PD、LSI和AI,以更直觀地揭示全國耕地景觀細碎度變化情況。具體來說,PD為單位面積的耕地斑塊數(shù)量,側(cè)重從區(qū)域內(nèi)耕地斑塊的相對大小和密集程度這一角度測量耕地景觀細碎度。PD越大,耕地斑塊面積越小,易耕性越差,耕地細碎度也越高。LSI為某一耕地斑塊形狀與相同面積的規(guī)則圖形(此處為正方形)之間的偏離程度,側(cè)重從耕地斑塊形狀規(guī)則程度這一視角衡量耕地細碎度。相同PD條件下,LSI越接近1,耕地斑塊形狀越接近易于耕作的正方形,否則耕地斑塊形狀越不規(guī)則,越難耕作,耕地景觀細碎度也越高。AI為耕地斑塊的空間連通程度,側(cè)重從耕地斑塊空間位置是否離散這一視角衡量耕地景觀細碎度。相同PD和LSI條件下,AI越大,耕地斑塊分布越集中連片,越容易實現(xiàn)規(guī)模利用;反之,耕地斑塊更趨分散分布,耕地景觀也更為細碎化??傮w而言,三個指標各有側(cè)重,又互為補充,能夠?qū)Ω鼐坝^細碎度進行較為全面地評估。

    具體計算式如下:

    式中PD、LSI和AI分別為縣域內(nèi)耕地的斑塊密度指數(shù)(個/km2)、斑塊形狀指數(shù)和聚集度指數(shù);N、CAE分別為縣域內(nèi)耕地的斑塊數(shù)量(個)、面積(km2)和斑塊邊界總長度(km);P為縣域內(nèi)耕地斑塊的周長(km)。一般而言,縣域的PD和LSI的值越高、AI的值越低,其耕地斑塊越細碎、形狀越不規(guī)則、空間連通性也越差,相應地耕地景觀細碎化程度越高;反之,則相反。

    考慮到本文耕地景觀細碎度指標的計算以縣域為基本單元,且涵蓋1990—2020年共計31個年份的數(shù)據(jù),常用的Fragstats軟件難以滿足如此巨大的計算任務,因此采用Python編程實現(xiàn)。

    1.2.2 趨勢分析

    自然地理學界通常采用趨勢分析獲取長時序的氣溫、降水、歸一化植被指數(shù)和凈初級生產(chǎn)力等自然地理要素的變化趨勢[25-27]。主要思路是通過一元線性回歸、非參數(shù)化趨勢度等方法進行上述地理要素的變化趨勢識別,同時借助T檢驗、Mann-Kendall統(tǒng)計檢驗等方法進行顯著性檢驗。本文將趨勢分析應用于縣域長時序耕地景觀細碎度指數(shù)的變化分析。具體來說,利用一元線性回歸模型對縣域耕地景觀細碎度指數(shù)進行擬合,通過斜率(Slope)綜合反映其變化趨勢;構(gòu)造統(tǒng)計量并計算其值,進而查閱分布表獲取可表征顯著性水平的值。相關(guān)計算式如下:

    式中Slope為縣域某個耕地景觀細碎度指數(shù)的變化趨勢;LF為縣域第年的某個耕地景觀細碎度指數(shù)值(PD、LSI或AI);為研究年份;e為樣本殘差;當Slope>0時,表示研究時段內(nèi)縣域該項耕地景觀細碎度指數(shù)呈增加趨勢;當Slope<0時,表示該時段內(nèi)縣域該項耕地景觀細碎度指數(shù)呈下降趨勢;根據(jù)值和分布表可獲取值,≤0.05表示顯著性變化;>0.05表示不顯著變化。

    1.2.3 動態(tài)面板模型

    在進行耕地景觀細碎度變化影響因素分析時,考慮到被解釋變量為長時序的耕地景觀細碎度指數(shù)數(shù)據(jù),且當年的耕地景觀細碎度指數(shù)與其上一年度情況具有較強的相關(guān)關(guān)系,也即被解釋變量受其滯后項的影響。因此,構(gòu)建計量經(jīng)濟模型時需將耕地景觀細碎度指數(shù)的滯后項作為關(guān)鍵影響因素進行考慮。此時,處理靜態(tài)面板數(shù)據(jù)時可以采用的混合OLS回歸和固定/隨機效應回歸方法往往難以解決被解釋變量滯后項的內(nèi)生性問題,需構(gòu)建動態(tài)面板模型,采用系統(tǒng)廣義矩估計(System GMM)方法對模型中的系數(shù)進行估計[28-29]。動態(tài)面板模型的基本設定如下:

    根據(jù)已有研究,自然和人為雙重因素共同導致了耕地景觀細碎化[3]。其中,自然因素包括耕地的地形、豐富度、區(qū)位等,人為因素包括區(qū)域土地利用情況、經(jīng)營主體的農(nóng)業(yè)生產(chǎn)和管理行為、地區(qū)經(jīng)濟發(fā)展水平以及政策因素等[3,17-20,24,30-31]。基于此,解釋變量選擇時考慮數(shù)據(jù)的可獲取性和代表性,本文以省域作為基本評估單元,除被解釋變量的滯后項外,重點考慮土地流轉(zhuǎn)政策對耕地景觀細碎化的影響;此外,還主要從區(qū)域土地利用情況、農(nóng)業(yè)生產(chǎn)管理情況、耕地資源稟賦特征、區(qū)域發(fā)展水平等方面選取可能影響耕地景觀細碎化的解釋變量。相關(guān)變量的定義和描述性統(tǒng)計分析結(jié)果見表1。

    表1 解釋變量的定義與描述性統(tǒng)計分析

    注:為避免內(nèi)生性,灌溉面積占比、第一產(chǎn)業(yè)增加值占比、降水、氣溫、高程、坡度、人均耕地面積均采用其一階滯后項;生態(tài)退耕主要為退耕還林(草);GDP取其自然對數(shù)納入模型中。

    Note: In order to avoid endogeneity, the proportion of irrigated area, the proportion of added value of primary industry, precipitation, temperature, elevation, slope and cultivated land per capita were all adopted first-order lag terms. Ecological conservation mainly includes conversion of cultivated land to forest (grass), and the natural logarithm of GDP is incorporated into the model.

    采用System GMM方法進行參數(shù)估計前,首先進行異方差、序列相關(guān)、截面相關(guān)等系列檢驗。檢驗結(jié)果顯示面板數(shù)據(jù)存在顯著的異方差、序列相關(guān)與截面相關(guān),因此需采用兩階段無偏穩(wěn)健性System GMM估計方法。此外,因被解釋變量滯后項存在內(nèi)生性問題,需采用工具變量進行參數(shù)估計。因此,進一步對模型模擬結(jié)果進行了工具變量設定合理性檢驗(Sargan檢驗和Hansen檢驗)和干擾項序列相關(guān)檢驗,前者用于驗證工具變量的選擇是否合理,后者則用于確定是否存在被解釋變量的二階序列相關(guān)。Sargan檢驗和Hansen檢驗的結(jié)果顯示各模型的工具變量設定均較為合理;干擾項序列相關(guān)檢驗的結(jié)果說明各模型僅存在一階序列相關(guān),因此僅將被解釋變量的滯后項納入模型即可。上述過程均借助Stata15相關(guān)命令實現(xiàn)。

    2 中國耕地景觀細碎度時空變化特征

    2.1 中國耕地景觀細碎度空間分布特征

    為揭示中國現(xiàn)狀耕地景觀細碎度的空間分布特征,圖1給出了縣域和省域尺度耕地景觀細碎度指數(shù)2018—2020年3年均值的空間分布圖。在縣域尺度(圖1a),PD低值區(qū)集中在東北平原、黃淮海平原、長江中下游平原、四川盆地、準噶爾盆地和塔里木盆地等地勢低平區(qū)域,PD高值區(qū)集中在第二級地勢階梯向第三級地勢階梯過渡地帶以及黃土高原、云貴高原和東南丘陵等地區(qū);LSI與PD的空間分布近似,但低值區(qū)有所收縮,更集中在東北平原、黃淮海平原、長江中下游平原和四川盆地,而高值區(qū)在二三級地勢階梯過渡地帶有所擴張;AI與PD大致呈反相位變化特征,即AI高值區(qū)與PD低值區(qū)在空間上基本吻合,主要分布在地勢低平的平原和盆地區(qū),而AI低值區(qū)與PD高值區(qū)在空間上也基本一致,主要分布在二三級地勢階梯過渡地帶和云貴高原、黃土高原以及東南丘陵區(qū)。因此,在縣域尺度,二三級地勢階梯過渡地帶以及黃土高原、云貴高原和東南丘陵等地區(qū)的耕地斑塊較細碎、形狀較不規(guī)則、空間連通性也較差,景觀細碎化程度較高;與之相比,地勢低平的平原盆地區(qū)耕地景觀細碎化程度則較高。

    注:該圖基于自然資源部標準地圖服務網(wǎng)站標準地圖制作,底圖邊界無修改。港澳臺數(shù)據(jù)暫缺。PD,斑塊密度;LSI,形狀指數(shù);AI,聚集度指數(shù)。下同。

    Note: This map is made based on the standard map download from the standard map service website of the Ministry of Natural Resources. The boundary of the base map is not modified. Data for Hong Kong, Macao and Taiwan are temporarily unavailable. PD, Patch Density; LSI, Landscape Shape Index; AI, Aggregation Index. Same as below.

    圖1 中國耕地景觀細碎度指數(shù)空間分布

    Fig.1 Spatial distribution of landscape fragmentation indices of cultivated land in China

    省域尺度上(圖1b),PD和AI仍呈明顯的反相位變化特征,河南、山東、安徽、河北、江蘇、四川、黑龍江、吉林、新疆、湖北、湖南等省區(qū)PD值較低而AI值較高,說明這些省份耕地斑塊規(guī)模較大且空間連通性高;貴州、甘肅、廣東、海南、云南、青海、北京、上海、福建等省份PD值較高而AI值較低,因此省域內(nèi)耕地斑塊較為細碎且連通性差。與PD和AI的省際差異略有不同,LSI值較高的省份除上述PD值較高和AI值較低的省份外,還包括湖南、重慶、內(nèi)蒙古、廣西等,說明這4個省區(qū)雖耕地斑塊規(guī)模較大且連通性較好,但形狀較不規(guī)則,青海、北京和上海等PD高值省市的LSI值則較低,說明盡管這3個省市耕地斑塊規(guī)模小且空間連通性差,但形狀較為規(guī)則。

    2.2 中國耕地景觀細碎度時序變化特征

    2.2.1 全國尺度

    1990—2020年中國耕地PD、LSI和AI均值整體均呈顯著下降趨勢,降幅分別為7.47個/(km2·10a)(=0.001)、1.13/10a(=0.026)和0.23/10a(=0.004)。具體來說,近30年中國耕地景觀細碎度變化可大致劃分為3個階段(圖2):1)持續(xù)細碎化階段(1990—2003年),PD和LSI均值顯著下降,降幅分別為3.52個/(km2·10a)和6.58/10a,AI均值顯著增加,增幅為0.67/10a;2)細碎化改善階段(2004—2013年),PD和LSI均值顯著上升,增幅分別為1.62個/(km2·10a)和2.59/10a,AI均值則顯著下降,降幅為1.03/10a;3)波動細碎化階段(2014 —2020年),PD和LSI的均值呈現(xiàn)波動下降趨勢,降幅分別為1.41個/(km2·10a)和7.60/10a,AI均值則波動增加,增幅為0.92/10a。

    圖2 1990—2020年中國耕地景觀細碎度指數(shù)平均值時序變化

    需要指出的是,兩個時間節(jié)點(2004年和2014年)分別為國家取消農(nóng)業(yè)稅和推行土地經(jīng)營權(quán)流轉(zhuǎn)政策的始點,可見政策因素可能對中國耕地景觀細碎度變化產(chǎn)生重要影響。考慮到本文重點考察土地流轉(zhuǎn)政策實施前后中國全域耕地景觀細碎度的變化情況,因此保持取消農(nóng)業(yè)稅這一政策背景不變,后文將主要從縣域和省域兩個尺度對2004—2020年耕地景觀細碎度指數(shù)的變化趨勢進行分析。

    2.2.2 省域尺度

    省域尺度上,2004—2013年絕大部分省份的耕地PD、LSI和AI顯著變化縣域占比超過50%,而2014—2020年各指數(shù)顯著變化縣域占比超過 50%的省份數(shù)量有所降低(圖3)。具體來說,2004—2013年,73.33%和80%省份的耕地PD和LSI整體呈增加趨勢,超90%省份的耕地AI整體呈下降趨勢;2014—2020年,耕地PD和LSI整體呈上升趨勢的省份占比降至不到1/3和1/5,AI整體呈下降趨勢的省份占比降至不足30%。因此,相較于2004 —2013年,2014—2020年90%省份的耕地景觀細碎度指數(shù)發(fā)生了部分或全部改變。其中,黑龍江、江西、遼寧、廣西、甘肅、廣東、海南和福建的耕地PD和LSI由整體增加轉(zhuǎn)為整體減少,而AI則由整體減少轉(zhuǎn)為整體增加,表征上述省份耕地景觀細碎度整體變化趨勢發(fā)生逆轉(zhuǎn),耕地斑塊細碎度、形狀規(guī)則性和空間連通性各方面均明顯改善。此外,僅安徽、寧夏、天津3個省級行政區(qū)在兩個時段的耕地景觀各指數(shù)變化趨勢保持一致。

    2.2.3 縣域尺度

    2004—2013年全國縣域耕地景觀細碎度指數(shù)變化圖顯示(圖4),全國大部分(>70%)縣域的耕地PD、LSI和AI呈現(xiàn)顯著變化趨勢(表2)。具體來說,PD和AI變化趨勢也呈明顯的反相位變化特征:PD顯著增加和AI顯著下降的縣域占比分別為50.13%和61.54%,主要分布在黃淮海平原、長江中下游平原和四川盆地以及兩廣地區(qū),說明該時段全國超 50%的縣域耕地斑塊更趨細碎且空間連通性變差;PD顯著減少和AI顯著增加的縣域占比分別為23.34%和15.86%,主要分布在東北平原、準噶爾盆地、塔里木盆地、二三級地勢過渡階梯地帶、黃土高原、云貴高原以及湖南、江西和福建地區(qū),說明這些區(qū)域耕地斑塊規(guī)模有所增加、連通性也有所改善(圖 4a、4b)。LSI變化趨勢的分布則有所不同:LSI顯著增加的縣域占比為57.39%,集中分布在黃淮海平原、長江中下游平原、四川盆地、塔里木盆地以及云貴高原和東南丘陵地區(qū),此外還散布在東北平原和準噶爾盆地區(qū),說明該時段全國超過50%縣級行政區(qū)的耕地形狀趨于不規(guī)則;LSI顯著減少的縣域占比為17.73%,集中分布在二三級地勢階梯過渡地帶、黃土高原和新疆西北部,說明上述地區(qū)耕地斑塊形狀更趨規(guī)整(圖4c)。

    圖3 各省份耕地景觀細碎度指數(shù)顯著變化縣域占比及整體趨勢

    表2 中國耕地景觀細碎度指數(shù)變化的縣域個數(shù)及比例

    總體來說,2004—2013年中國耕地景觀整體呈現(xiàn)細碎化趨勢,以黃淮海平原、長江中下游平原、四川盆地以及兩廣地區(qū)表現(xiàn)尤為典型,二三級地勢階梯過渡地帶、黃土高原和新疆西北部的細碎化程度則趨于改善。

    2014—2020年全國縣域耕地景觀細碎度指數(shù)變化情況顯示,PD、LSI和AI顯著變化的縣域占比分別為58.11%,56.24%和39.25%(表2)。具體來說,PD和AI仍大致呈反相位變化特征,PD顯著增加和AI顯著下降的縣域占比僅為20.91%和15.66%,均比2004—2013年的比例要低,空間上集中分布在黃淮海平原和四川盆地,這些地區(qū)LSI也大致呈增加趨勢,因此耕地景觀整體仍向細碎化方向發(fā)展。PD顯著減少和AI顯著增加的縣域占比分別為37.2%和23.59%,均比2004—2013年的比例要高,空間集中分布區(qū)除2004—2013時段的同方向變化區(qū)域外,還包括了兩廣地區(qū),這些地區(qū)LSI也呈顯著減少趨勢,因此耕地景觀細碎化程度有所改善。除上述區(qū)域外,全國一半左右縣域PD、LSI和AI變化趨勢不顯著,集中分布在長江中下游平原和西北地區(qū)。

    a. PD變化趨勢a. PD change trendb. LSI變化趨勢b. LSI change trendc. AI變化趨勢c. AI change trend

    d. PD顯著性d. PD significancee. LSI顯著性e. LSI significancef. AI顯著性f. AI significance

    a. PD變化趨勢a. PD change trendb. LSI變化趨勢b. LSI change trendc. AI變化趨勢c. AI change trend

    d. PD顯著性d. PD significancee. LSI顯著性e. LSI significancef. AI顯著性f. AI significance

    3 耕地景觀細碎度變化的影響因素

    考慮到土地流轉(zhuǎn)政策對耕地景觀細碎度變化的影響是本文研究的重點之一,在進行影響因素分析時,為盡可能降低國家取消農(nóng)業(yè)稅政策對耕地景觀細碎度影響的干擾,選擇2008—2020年作為研究時段,此時段國家取消農(nóng)業(yè)稅政策已實施5年,其對耕地景觀細碎化的年際影響趨于穩(wěn)定,耕地細碎度變化更有可能是土地流轉(zhuǎn)政策和其他因素共同作用的結(jié)果。將耕地景觀細碎度各指數(shù)分別作為被解釋變量構(gòu)建動態(tài)面板模型,利用System GMM方法對解釋變量的系數(shù)進行估計。最終得到的穩(wěn)健型估計結(jié)果見表3。需要指出的是,考慮一階滯后項的動態(tài)面板模型中,其他變量主要用于解釋耕地景觀細碎度指數(shù)的年際變化情況。

    1)土地流轉(zhuǎn)政策是耕地景觀細碎度變化的重要影響因素。根據(jù)表3,土地流轉(zhuǎn)政策對PD和LSI的影響均顯著為負,對AI的影響則顯著為正。具體來說,土地流轉(zhuǎn)政策實施后,PD和LSI整體減少了0.137個/km2和0.526,而AI整體增加了0.136,可見土地流轉(zhuǎn)政策實施后,土地確權(quán)與制度改革、新型農(nóng)業(yè)經(jīng)營主體培育、農(nóng)業(yè)社會化服務體系建設等推進土地經(jīng)營權(quán)流轉(zhuǎn)的各項舉措效果顯著,在降低耕地景觀細碎度,特別是在規(guī)整地塊、改善空間連通度等方面產(chǎn)生了積極效果,有效改善了中國的耕地細碎化狀況。

    表3 動態(tài)面板模型的穩(wěn)健型回歸結(jié)果

    注:*、**和***分別表示估計參數(shù)在10%、5%和1%的水平上顯著。

    Note: *, ** and *** indicate that the estimated parameters are significant at the 10%, 5% and 1% levels, respectively.

    2)區(qū)域土地利用變化直接且顯著影響耕地細碎度。具體來說,耕地擴張和建設用地占用耕地均對PD、LSI和AI產(chǎn)生了顯著影響,耕地每擴張1×104hm2,PD和LSI分別增加0.010個/km2和0.026,而AI減少0.006,而建設用地每占用1×104hm2耕地,PD和LSI分別增加0.067個/km2和0.687,AI則減少0.173。這與已有研究的結(jié)論是一致的,即耕地擴張和建設用地擴張均使得耕地斑塊更細碎、邊界更不規(guī)則、空間連通度更差,耕地景觀細碎度上升[20,31]。生態(tài)退耕主要對LSI和AI產(chǎn)生顯著影響,且生態(tài)退耕面積每增加1×104hm2,LSI減少0.111,AI增加0.016。這可能是因為生態(tài)退耕主要面向山區(qū)坡度大于25°的耕地,這些耕地本身細碎化程度較高,其退出耕作后反而使得區(qū)域耕地細碎度降低、空間連通度提升[20]。

    3)非糧化耕地利用行為也顯著促進了耕地細碎化。第一產(chǎn)業(yè)增加值占比對PD的影響是顯著負向的,說明農(nóng)業(yè)產(chǎn)值占比較高的地區(qū)更有可能降低耕地細碎度,實現(xiàn)耕地適度規(guī)模經(jīng)營。然而,非糧化率對PD和LSI產(chǎn)生了顯著正面影響,對AI產(chǎn)生了顯著負面影響,說明非糧化耕地利用行為使得耕地景觀更為細碎。非糧化率每提高1個百分點,PD和LSI分別增加0.002個/km2和0.007,AI則減少0.001。這應該與非糧化耕地利用行為能在單位面積土地上獲得高額利潤有關(guān),相應農(nóng)戶進行規(guī)模經(jīng)營的意愿因而較低;且非糧作物,特別是瓜果蔬菜等多采用大棚種植,對耕地景觀細碎度,包括斑塊大小、形狀和空間連通性產(chǎn)生負面影響。

    4)坡度對耕地景觀細碎度變化產(chǎn)生了一定影響。具體來說,坡度主要對PD產(chǎn)生顯著負向影響,說明耕地平均坡度越大的地區(qū),PD越有可能呈現(xiàn)減少趨勢。這可能與山區(qū)坡耕地撂荒和生態(tài)退耕政策的實施有關(guān)。已有研究表明,坡度越大且越細碎的地塊越容易被撂荒[32-34],加之這些地塊同時也是生態(tài)退耕政策需要退耕的目標,因此,隨著撂荒和退耕政策的實施,相關(guān)區(qū)域耕地的PD有所減小,耕地細碎度也有所降低。

    5)區(qū)域發(fā)展水平等因素對耕地景觀細碎度變化的影響并不顯著。GDP、城市化率、人均耕地面積、林業(yè)工程建設、灌溉面積占比、除坡度外的其他耕地資源稟賦特征等因素的系數(shù)均不顯著,說明其對耕地景觀細碎度的影響較為有限。需要指出的是,灌溉面積占比的各項系數(shù)接近顯著性臨界值(0.1),說明其對耕地細碎度仍產(chǎn)生了一定的影響。具體來說,灌溉面積占比每增加1個百分點,PD和LSI分別減少0.004個/km2和0.008,AI增加0.001,也即灌溉條件的改善可能會降低耕地景觀細碎度,對耕地適度規(guī)模經(jīng)營具有一定的促進作用。

    4 討 論

    2014年國家提出推行土地經(jīng)營權(quán)流轉(zhuǎn)和耕地適度規(guī)模經(jīng)營,以改變中國細碎化耕地的低效利用模式,推進農(nóng)業(yè)現(xiàn)代化進程,實現(xiàn)農(nóng)民增產(chǎn)增收。然而,很少學者針對長時序耕地細碎度變化開展研究,中國耕地細碎度時空變化特征以及土地經(jīng)營權(quán)流轉(zhuǎn)政策在改善耕地細碎度方面的效果均尚未可知。本文主要從宏觀尺度出發(fā),借助長時序中高分辨率土地利用數(shù)據(jù)集(CLCD),進行了中國全域耕地景觀細碎度變化的趨勢分析,并對土地流轉(zhuǎn)政策等因素對耕地景觀細碎度變化的可能影響進行了剖析。一方面使得人們對中國耕地景觀細碎度時空變化特征有了一定認識,另一方面初步揭示了土地流轉(zhuǎn)政策等因素對耕地景觀細碎度的可能影響,同時也為區(qū)域耕地景觀細碎度變化以及不同土地利用/生態(tài)系統(tǒng)類型變化分析提供了新思路。

    考慮到中國耕地景觀細碎度變化存在顯著的空間差異,特別是土地經(jīng)營權(quán)流轉(zhuǎn)政策實施之后(2014—2020年),盡管70%省份的耕地景觀細碎度整體改善,但在黃淮海平原和四川盆地仍呈現(xiàn)PD增加和AI減小趨勢,這些地區(qū)應作為重點區(qū)域,加強新型農(nóng)業(yè)經(jīng)營主體培育和農(nóng)業(yè)社會化服務體系建設,完善土地流轉(zhuǎn)市場,以提升該區(qū)域的政策實施效果[6,35-36]。在該區(qū)域以及全國層面上,還應避免耕地無序開發(fā)并加強耕地整治,嚴格管控耕地非農(nóng)占用,同時推進農(nóng)業(yè)良種改革和優(yōu)化農(nóng)業(yè)補貼政策以降低非糧化程度,并保護山區(qū)退耕還林(草)政策實施效果,以進一步降低耕地細碎度[20]。此外,以灌溉條件改善為代表的耕地質(zhì)量提升措施也能在一定程度上改善耕地細碎度,因此,建議加大全國高標準農(nóng)田建設的力度和廣度,進一步激發(fā)區(qū)域土地流轉(zhuǎn)潛力,推進耕地適度規(guī)模經(jīng)營。

    本文的不足之處,一是采用了CLCD數(shù)據(jù)集作為耕地景觀細碎度分析的關(guān)鍵基礎數(shù)據(jù),該數(shù)據(jù)集是遙感解譯數(shù)據(jù),主要表達的是地物光譜信息和覆蓋信息,可能存在同一地塊因種植作物不同或田埂存在而被劃分成不同斑塊的情況,也可能存在同一片梯田按丘劃塊等情況,能否準確表達耕地細碎度及其變化情況值得商榷。然而,考慮到該數(shù)據(jù)集是目前時空分辨率和總體精度均較高的土地利用數(shù)據(jù)集,利用其進行耕地景觀細碎度變化研究可被視為一次有益的嘗試,由此獲取的中國全域耕地景觀細碎度變化特征也可為相關(guān)研究提供一定的參考。今后,可開展與其他數(shù)據(jù)源研究結(jié)果的交互驗證工作,以更全面、精準地認識中國耕地細碎化特征。二是受數(shù)據(jù)可獲取性限制,本文以省域作為基本單元進行了影響因素評估,識別了土地流轉(zhuǎn)政策對該層面耕地景觀細碎度變化的影響以及其他影響因素,今后可將景觀細碎化與微觀農(nóng)業(yè)經(jīng)營主體相關(guān)研究結(jié)合,并構(gòu)建理論模型,從多層面、多視角對耕地細碎度變化進行解釋。

    5 結(jié) 論

    利用中國1990—2020年長時間序列土地利用數(shù)據(jù)集,輔以景觀格局指數(shù)、趨勢分析和動態(tài)面板模型等方法,剖析了中國耕地景觀細碎度指數(shù)的時空變化特征及其影響因素。本文得出的主要結(jié)論如下:

    1)中國耕地景觀細碎度大致呈“平原盆地低、山地高原高”的分布格局,斑塊密度(PD)/形狀指數(shù)(LSI)低值區(qū)和聚集度指數(shù)(AI)高值區(qū)主要分布在東北平原、黃淮海平原、長江中下游平原和四川盆地以及新疆的準噶爾和塔里木盆地等,PD/LSI高值區(qū)和AI低值區(qū)則集中在二三級地勢階梯過渡地帶以及黃土高原、云貴高原和東南丘陵等地區(qū)。

    2)1990—2020年中國耕地景觀細碎度呈現(xiàn)“持續(xù)細碎化-細碎化改善-波動細碎化”三階段變化特征。土地經(jīng)營權(quán)流轉(zhuǎn)政策實施后(2014—2020年),PD和LSI均值呈波動下降趨勢,降幅分別為1.41個/(km2·10 a)和7.60/10 a,AI均值則波動上升,增幅為0.92/10 a。

    3)2004—2013年中國耕地景觀整體趨于細碎化,PD、LSI顯著增加的縣域占比(50.13%、57.39%)高于顯著減少的縣域占比(23.34%和17.73%),AI顯著增加的縣域占比(15.86%)低于顯著減少的縣域占比(61.54%)。但2014—2020年90%省份的耕地景觀細碎度有所改善。

    4)中國耕地景觀細碎度變化受多種因素的共同影響。其中,土地流轉(zhuǎn)政策是影響耕地景觀細碎度變化的重要因素,政策實施后全國耕地景觀細碎度整體呈下降趨勢。此外,耕地景觀細碎度變化還受區(qū)域土地利用變化、非糧化耕地利用行為、坡度以及灌溉條件等因素的共同影響。

    [1] Zhou Y, Li X, Liu Y. Cultivated land protection and rational use in China[J]. Land Use Policy, 2021, 106: 105454.

    [2] Lai Z, Chen M, Liu T. Changes in and prospects for cultivated land use since the reform and opening up in China[J]. Land Use Policy, 2020, 97: 104781.

    [3] 杜國明,蓋兆雪,王洪彥. 中國耕地細碎化的理論解析與研究框架[J]. 地球科學與環(huán)境學報,2021,43(6):997-1008.

    Du Guoming, Gai Zhaoxue, Wang Hongyan. Theoretical explanation and research framework of cultivated-land fragmentation in China[J]. Journal of Earth Sciences and Environment, 2021, 43(6): 997-1008. (in Chinese with English abstract)

    [4] 盧華,胡浩. 土地細碎化增加農(nóng)業(yè)生產(chǎn)成本了嗎?來自江蘇省的微觀調(diào)查[J]. 經(jīng)濟評論,2015,195(5):129-140.

    [5] 秦立建,張妮妮,蔣中一. 土地細碎化、勞動力轉(zhuǎn)移與中國農(nóng)戶糧食生產(chǎn):基于安徽省的調(diào)查[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟,2011,(11):16-23.

    [6] 王學,李秀彬,辛良杰,等. 華北平原農(nóng)戶主動退耕冬小麥的影響因素模型分析及政策啟示[J]. 農(nóng)業(yè)工程學報,2018,34(9):248-257.

    Wang Xue, Li Xiubin, Xin Liangjie, et al. Driving factors on winter wheat abandonment by farmers in North China Plain and its policy implication[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(9): 248-257. (in Chinese with English abstract)

    [7] Yan J, Yang Z, Li Z, et al. Drivers of cropland abandonment in mountainous areas: A household decision model on farming scale in Southwest China[J]. Land Use Policy, 2016, 57: 459-469.

    [8] 農(nóng)業(yè)農(nóng)村部政策與改革司. 2020中國衣村政策與改革統(tǒng)計年報[M]. 北京:中國農(nóng)業(yè)出版社,2020.

    [9] Wang Y, Li X, Xin L, et al. Spatiotemporal changes in Chinese land circulation between 2003 and 2013[J]. Journal of Geographical Sciences, 2018, 28(6): 707-724.

    [10] Xie H, Lu H. Impact of land fragmentation and non-agricultural labor supply on circulation of agricultural land management rights[J]. Land Use Policy, 2017, 68: 355-364.

    [11] 王春超,李兆能. 農(nóng)村土地流轉(zhuǎn)中的困境:來自湖北的農(nóng)戶調(diào)查[J]. 華中師范大學學報(人文社會科學版),2008,47(4):51-56.

    [12] Cao Y, Zou J, Fang X, et al. Effect of land tenure fragmentation on the decision-making and scale of agricultural land transfer in China[J]. Land Use Policy, 2020, 99: 104996.

    [13] 王亞輝,李秀彬,辛良杰. 耕地地塊細碎程度及其對山區(qū)農(nóng)業(yè)生產(chǎn)成本的影響[J]. 自然資源學報,2019,34(12):2658-2672.

    Wang Yahui, Li Xiubin, Xin Liangjie. Characteristics of cropland fragmentation and its impact on agricultural production costs in mountainous areas[J]. Journal of Natural Resources, 2019, 34(12): 2658-2672. (in Chinese with English abstract)

    [14] Lu H, Xie H, He Y, et al. Assessing the impacts of land fragmentation and plot size on yields and costs: A translog production model and cost function approach[J]. Agricultural Systems, 2018, 161: 81-88.

    [15] 王嫚嫚,劉穎,蒯昊,等. 土地細碎化、耕地地力對糧食生產(chǎn)效率的影響:基于江漢平原354個水稻種植戶的研究[J].資源科學,2017,39(8):1488-1496.

    Wang Manman, Liu Ying, Kuai Hao, et al. The effects of land fragmentation and land quality on the technical efficiency of grain production based on 354 rice planters on the Jianghan Plain[J]. Resources Science, 2017, 39(8): 1488-1496. (in Chinese with English abstract)

    [16] Cheng L, Xia N, Jiang P, et al. Analysis of farmland fragmentation in China modernization demonstration zone since “Reform and Openness”: A case study of South Jiangsu Province[J]. Scientific Reports, 2015, 5: 11797.

    [17] Liu J, Jin X, Xu W, et al. Evolution of cultivated land fragmentation and its driving mechanism in rural development: A case study of Jiangsu Province[J]. Journal of Rural Studies, 2022, 91: 58-72.

    [18] Xu W, Jin X, Liu J, et al. Analysis of influencing factors of cultivated land fragmentation based on hierarchical linear model: A case study of Jiangsu Province, China[J]. Land Use Policy, 2021, 101: 105119.

    [19] 劉立文,段永紅,徐立帥,等. 黃土丘陵區(qū)耕地細碎化特征及影響因素研究:以山西省和順縣為例[J]. 中國農(nóng)學通報,2021,37(22):90-96.

    Liu Liwen, Duan Yonghong, Xu Lishuai, et al. Characteristics and influencing factors of cultivated land fragmentation in Loess Hilly Region: A case study of Heshun County in Shanxi Province[J]. Chinese Agricultural Science Bulletin, 2021, 37(22): 90-96. (in Chinese with English abstract)

    [20] Yu Q, Hu Q, van Vliet J, et al. GlobeLand30 shows little cropland area loss but greater fragmentation in China[J]. International Journal of Applied Earth Observation and Geoinformation, 2018, 66: 37-45.

    [21] 劉紀遠,寧佳,匡文慧,等. 2010—2015年中國土地利用變化的時空格局與新特征[J]. 地理學報,2018,73(5):789-802.

    Liu Jiyuan, Ning Jia, Kuang Wenhui, et al. Spatio-temporal patterns and characteristics of land-use change in China during 2010-2015[J]. Acta Geographica Sinica, 2018, 73(5): 789-802. (in Chinese with English abstract)

    [22] Yang J, Huang X. 30 m annual land cover and its dynamics in China from 1990 to 2019 [J]. Earth System Science Data, 2021: 13(8): 3907-3925.

    [23] Wei L, Luo Y, Wang M, et al. Essential fragmentation metrics for agricultural policies: Linking landscape pattern, ecosystem service and land use management in urbanizing China[J]. Agricultural Systems, 2020, 182: 102833.

    [24] Liang J, Pan S, Chen W, et al. Cultivated land fragmentation and its influencing factors detection: A case study in Huaihe River Basin, China[J]. International Journal of Environmental Research and Public Health, 2021, 19(1): 138.

    [25] 劉憲鋒,潘耀忠,朱秀芳,等. 2000-2014年秦巴山區(qū)植被覆蓋時空變化特征及其歸因[J]. 地理學報,2015,70(5):705-716.

    Liu Xianfeng, Pan Yaozhong, Zhu Xiufang, et al. Spatiotemporal variation of vegetation coverage in Qinling-Daba Mountains in relation to environmental factors[J]. Acta Geographica Sinica, 2015, 70(5): 705-716. (in Chinese with English abstract)

    [26] 金凱,王飛,韓劍橋,等. 1982—2015年中國氣候變化和人類活動對植被NDVI變化的影響[J]. 地理學報,2020,75(5):961-974.

    Jin Kai, Wang Fei, Han Jianqiao, et al. Contribution of climatic change and human activities to vegetation NDVI change over China during 1982-2015[J]. Acta Geographica Sinica, 2020, 75(5): 961-974. (in Chinese with English abstract)

    [27] 程小強,周兆葉,李旺平,等. 基于MODIS數(shù)據(jù)的中亞地區(qū)旱情監(jiān)測及影響因素分析[J]. 農(nóng)業(yè)工程學報,2022,38(10):128-137.

    Cheng Xiaoqiang, Zhou Zhaoye, Li Wangping, et al. Monitoring drought situation and analyzing influencing factors in Central Asia using MODIS data[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(10): 128-137. (in Chinese with English abstract)

    [28] 董洪超,蔣伏心. 交通基礎設施對中國區(qū)域市場一體化的影響研究:基于動態(tài)面板模型的實證分析[J]. 經(jīng)濟問題探索,2020(5):26-39.

    [29] 陳沖,吳煒聰. 人口老齡化對農(nóng)業(yè)技術(shù)進步的影響機制分析:基于DEA-Malmquist的技術(shù)評價與動態(tài)面板模型[J]. 中國農(nóng)業(yè)資源與區(qū)劃,2021,42(1):231-238.

    Chen Chong, Wu Weicong. Analysis on the influence mechanism of population aging on agricultural technological progress[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2021, 42(1): 231-238. (in Chinese with English abstract)

    [30] Liu Y, Feng Y, Zhao Z, et al. Socioeconomic drivers of forest loss and fragmentation: A comparison between different land use planning schemes and policy implications[J]. Land Use Policy, 2016, 54: 58-68.

    [31] Tan S, Heerink N, Qu F. Land fragmentation and its driving forces in China[J]. Land Use Policy, 2006, 23: 272-285.

    [32] Shi T, Li X, Xin L, et al. The spatial distribution of farmland abandonment and its influential factors at the township level: A case study in the mountainous area of China[J]. Land Use Policy, 2018, 70: 510-520.

    [33] Li Y, Zhou T, Jiang G, et al. Spatial pattern and mechanisms of farmland abandonment in Agricultural and Pastoral Areas of Qingzang Plateau[J]. Geography and Sustainability, 2021, 2(3): 139-150.

    [34] Li S, Li X. Global understanding of farmland abandonment: A review and prospects[J]. Journal of Geographical Sciences, 2017, 27(9): 1123-1150.

    [35] 王學,李秀彬,辛良杰,等. 華北地下水超采區(qū)冬小麥退耕的生態(tài)補償問題探討[J]. 地理學報,2016,71(5):829-839.

    Wang Xue, Li Xiubin, Xin Liangjie, et al. Ecological compensation for winter wheat abandonment in groundwater over-exploited areas in the North China Plain[J]. Acta Geographica Sinica, 2016, 71(5): 829-839. (in Chinese with English abstract)

    [36] 牛海鵬,孫壹鳴. 農(nóng)戶土地承包經(jīng)營權(quán)退出意愿影響因素及退出模式研究[J]. 農(nóng)業(yè)工程學報,2019,35(11):265-275.

    Niu Haipeng, Sun Yiming. Factors affecting farmers’ willingness and mode of farmland usufruct abandonment for rural households land contractual operation right [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(11): 265-275. (in Chinese with English abstract)

    Spatiotemporal characteristics and influencing factors of landscape fragmentation of cultivated land in China

    Wang Xue1, Xu Xiaofan1,2

    (1.,,,100101,;2.,100049,)

    Cultivated land fragmentation has been widely known as one of the key limiting factors for cultivated land use efficiency, as well as agricultural modernization and transformation in China. Since 2014, much effort has been made to encourage the transfer of land management rights for the less cultivated land fragmentation, even for the cultivated land at a moderate scale. Most studies have also been conducted at the local or regional scales. However, only a few studies are focused on the fine-grained changes in cultivated land fragmentation from a national perspective. It is still lacking in the implementation effect of land management right transfer policy on the cultivated land fragmentation. In this study, a systematic analysis was made of the temporal and spatial variation characteristics of China's cultivated land fragmentation from 1990 to 2020. The influencing factors were also identified near the implementation of the land transfer policy (2008-2020). The long-term land use data was selected to supplement the landscape pattern index, trend analysis, and dynamic panel model. The cultivated land fragmentation was roughly divided into ownership fragmentation and landscape fragmentation, in terms of connotation. The landscape fragmentation of cultivated land (CLF) was then used in this case. The results showed that: 1) China’s CLF presented a distribution pattern of low in the plain and basin areas, and high in the mountainous and plateau areas. Specifically, the CLF was relatively low in the Northeast Plain, the Huang-Huai-Hai Plain, the Middle and Lower Reaches of the Yangtze River Plain, as well as the Sichuan Basin, Junggar and Tarim Basin; while the cultivated land in the transition zone of the second and third-level topographic steps, as well as the Loess Plateau, Yunnan-Guizhou Plateau, and the Southeastern Hills, was relatively fragmented. 2) China’s CLF generally showed a changing trend of continuous-improved-fluctuating fragmentation from 1990 to 2020. Specifically, the overall cultivated land showed a trend of landscape fragmentation from 2004 to 2013. More than 70% of the counties showed a significant increase in the PD/LSI or a significant decrease in the AI, which were concentrated in the Huang-Huai-Hai Plain, the Middle and Lower Reaches of the Yangtze River Plain, the Sichuan Basin, and the Guangdong-Guangzhou region. However, 70% of the provinces improved the landscape fragmentation of cultivated land from 2014 to 2020. Spatially, the counties with the improved CLF (that is, the counties with the significantly reduced PD and LSI, but increased AI) were mainly located in the transition zone of the second and third-level topographic steps, the Loess Plateau, northwest Xinjiang, and the Guangdong-Guangzhou region. 3) The land transfer policy significantly reduced the CLF. In addition, some important factors were the regional land use change, non-grained utilization of cultivated land, slope, and irrigation conditions. The spatial-temporal variation characteristics of CLF in China on a nationwide scale revealed the impact of land transfer policies and other factors on the CLF. The finding can provide a new research paradigm for the changes in the regional CLF and different land use/ecosystem types.

    land use; spatial distribution; dynamic models; cultivate land fragmentation; landscape pattern index; trend analysis; land transfer

    10.11975/j.issn.1002-6819.2022.16.002

    S17

    A

    1002-6819(2022)-16-0011-10

    王學,徐曉凡. 中國耕地景觀細碎度時空變化特征及其影響因素[J]. 農(nóng)業(yè)工程學報,2022,38(16):11-20.doi:10.11975/j.issn.1002-6819.2022.16.002 http://www.tcsae.org

    Wang Xue, Xu Xiaofan. Spatiotemporal characteristics and influencing factors of landscape fragmentation of cultivated land in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(16): 11-20. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.16.002 http://www.tcsae.org

    2022-06-07

    2022-08-09

    國家自然科學基金項目(42171264)

    王學,博士,副研究員,研究方向為土地利用變化。Email:wangxue@igsnrr.ac.cn

    猜你喜歡
    縣域耕地斑塊
    自然資源部:加強黑土耕地保護
    我國將加快制定耕地保護法
    保護耕地
    北京測繪(2021年12期)2022-01-22 03:33:36
    頸動脈的斑塊逆轉(zhuǎn)看“軟硬”
    自我保健(2021年2期)2021-11-30 10:12:31
    新增200億元列入耕地地力保護補貼支出
    一篇文章了解頸動脈斑塊
    婦女之友(2021年9期)2021-09-26 14:29:36
    縣域消防專項規(guī)劃研究
    microRNA-146a與冠心病患者斑塊穩(wěn)定性的相關(guān)性
    有頸動脈斑塊未必要吃降脂藥
    百姓生活(2019年2期)2019-03-20 06:06:16
    山東縣域GDP排名出爐
    商周刊(2017年6期)2017-08-22 03:42:46
    秋霞伦理黄片| 国产在视频线精品| 18禁在线播放成人免费| 久久久久人妻精品一区果冻| 免费观看无遮挡的男女| 日本爱情动作片www.在线观看| 嫩草影院入口| videossex国产| 老司机影院成人| 高清毛片免费看| 国产不卡av网站在线观看| 人妻 亚洲 视频| 免费看不卡的av| 夫妻性生交免费视频一级片| 亚洲精品视频女| 视频在线观看一区二区三区| 日日撸夜夜添| 有码 亚洲区| 激情五月婷婷亚洲| 亚洲人与动物交配视频| 国产精品久久久久成人av| 国产 精品1| 纯流量卡能插随身wifi吗| 国产欧美日韩一区二区三区在线 | 看免费成人av毛片| 美女大奶头黄色视频| 日韩不卡一区二区三区视频在线| 五月玫瑰六月丁香| 五月玫瑰六月丁香| 免费久久久久久久精品成人欧美视频 | 狠狠精品人妻久久久久久综合| 久久久久久人妻| 一级爰片在线观看| 亚洲国产精品999| 少妇高潮的动态图| 久久国产精品大桥未久av| 蜜臀久久99精品久久宅男| 这个男人来自地球电影免费观看 | 天天操日日干夜夜撸| freevideosex欧美| 国产亚洲午夜精品一区二区久久| 秋霞在线观看毛片| 欧美+日韩+精品| 国产极品天堂在线| 国产精品女同一区二区软件| 国产一级毛片在线| 国产在线免费精品| 久久久久久久久大av| 99久久综合免费| 99热这里只有精品一区| 国产精品久久久久久久电影| 99久久中文字幕三级久久日本| 99久久综合免费| 国产深夜福利视频在线观看| 九草在线视频观看| 亚洲国产精品专区欧美| 国产精品久久久久久av不卡| 在线观看人妻少妇| 我的老师免费观看完整版| 五月玫瑰六月丁香| 亚洲精品第二区| 久久久精品免费免费高清| 久久99一区二区三区| 国产精品人妻久久久影院| 免费观看的影片在线观看| 日本91视频免费播放| 99国产综合亚洲精品| 日本vs欧美在线观看视频| 毛片一级片免费看久久久久| 伊人久久国产一区二区| 日本免费在线观看一区| 国产精品欧美亚洲77777| 人妻少妇偷人精品九色| h视频一区二区三区| 91久久精品电影网| 91久久精品国产一区二区三区| 97精品久久久久久久久久精品| 久久久精品免费免费高清| 亚洲精品亚洲一区二区| 青春草亚洲视频在线观看| 三级国产精品片| 日日摸夜夜添夜夜添av毛片| 3wmmmm亚洲av在线观看| 国产免费视频播放在线视频| 日本欧美视频一区| 国产成人精品无人区| 大片电影免费在线观看免费| 少妇熟女欧美另类| videosex国产| 国产欧美另类精品又又久久亚洲欧美| 成年人午夜在线观看视频| 亚洲精品久久成人aⅴ小说 | 亚洲精品日本国产第一区| 亚洲国产毛片av蜜桃av| 亚洲精品亚洲一区二区| 黄色怎么调成土黄色| 久久久久久久国产电影| 久久99热6这里只有精品| 一本大道久久a久久精品| 国产精品99久久99久久久不卡 | 人人妻人人爽人人添夜夜欢视频| 成人毛片a级毛片在线播放| 又大又黄又爽视频免费| 一级毛片 在线播放| 狂野欧美激情性xxxx在线观看| 欧美日韩av久久| 欧美 日韩 精品 国产| 蜜桃在线观看..| 九九在线视频观看精品| 超色免费av| 人人妻人人澡人人看| 蜜臀久久99精品久久宅男| 日产精品乱码卡一卡2卡三| 日本wwww免费看| 亚洲欧美日韩卡通动漫| 久久久久久久亚洲中文字幕| 大香蕉久久成人网| av国产久精品久网站免费入址| 精品亚洲成国产av| 亚洲av国产av综合av卡| 啦啦啦中文免费视频观看日本| 九九爱精品视频在线观看| 精品国产一区二区三区久久久樱花| 日韩三级伦理在线观看| 国产69精品久久久久777片| 伊人久久精品亚洲午夜| 亚洲精品乱久久久久久| 丰满饥渴人妻一区二区三| 国产一区二区三区综合在线观看 | 亚洲色图综合在线观看| 国产免费又黄又爽又色| 亚洲美女黄色视频免费看| 男女边吃奶边做爰视频| 五月天丁香电影| 赤兔流量卡办理| 国产 一区精品| 国产成人免费无遮挡视频| 免费黄频网站在线观看国产| 久久99蜜桃精品久久| 一区二区三区四区激情视频| 欧美97在线视频| 日产精品乱码卡一卡2卡三| 日本猛色少妇xxxxx猛交久久| 亚洲,欧美,日韩| 少妇被粗大的猛进出69影院 | 九色亚洲精品在线播放| 久久久久久人妻| 国产一区二区三区综合在线观看 | 亚洲国产精品专区欧美| 国模一区二区三区四区视频| 亚洲精品456在线播放app| 亚洲性久久影院| 亚洲av欧美aⅴ国产| 久久久久视频综合| 亚洲综合精品二区| 一级爰片在线观看| 成人毛片a级毛片在线播放| 一级毛片 在线播放| 18禁在线无遮挡免费观看视频| 丝袜在线中文字幕| 亚洲精品美女久久av网站| 日韩一本色道免费dvd| a级毛片免费高清观看在线播放| 日韩 亚洲 欧美在线| 免费少妇av软件| 你懂的网址亚洲精品在线观看| 国产国语露脸激情在线看| 日日摸夜夜添夜夜爱| 国产 精品1| 亚洲国产毛片av蜜桃av| 精品国产露脸久久av麻豆| 色婷婷av一区二区三区视频| 国产免费一区二区三区四区乱码| 久久久久久久大尺度免费视频| 免费看不卡的av| 母亲3免费完整高清在线观看 | 亚洲av综合色区一区| 91精品三级在线观看| 母亲3免费完整高清在线观看 | 亚洲国产欧美日韩在线播放| av专区在线播放| 在线天堂最新版资源| 日本wwww免费看| 在线观看美女被高潮喷水网站| 亚洲精品日本国产第一区| 亚洲av在线观看美女高潮| 夜夜爽夜夜爽视频| 五月开心婷婷网| 一区二区三区精品91| 日本猛色少妇xxxxx猛交久久| 婷婷色综合www| av在线app专区| 色94色欧美一区二区| 五月伊人婷婷丁香| 成人漫画全彩无遮挡| 最近的中文字幕免费完整| 久久久欧美国产精品| 亚洲欧洲国产日韩| 美女国产高潮福利片在线看| 人人妻人人爽人人添夜夜欢视频| 成人黄色视频免费在线看| 视频中文字幕在线观看| 国产乱来视频区| 18+在线观看网站| 制服诱惑二区| 国产av国产精品国产| 三上悠亚av全集在线观看| 免费大片黄手机在线观看| 亚洲欧美成人综合另类久久久| 欧美日韩国产mv在线观看视频| 久久久国产欧美日韩av| 大陆偷拍与自拍| 大码成人一级视频| 一级黄片播放器| 亚洲欧美成人精品一区二区| 国产精品偷伦视频观看了| 精品少妇久久久久久888优播| 亚洲成人一二三区av| 人妻制服诱惑在线中文字幕| 午夜av观看不卡| 极品人妻少妇av视频| 一二三四中文在线观看免费高清| 亚洲精品久久成人aⅴ小说 | 人妻系列 视频| 精品国产露脸久久av麻豆| 一级二级三级毛片免费看| 国产亚洲精品第一综合不卡 | 亚洲情色 制服丝袜| 91成人精品电影| 日产精品乱码卡一卡2卡三| 精品久久久精品久久久| 伊人久久精品亚洲午夜| 九色成人免费人妻av| 午夜精品国产一区二区电影| 街头女战士在线观看网站| 日本欧美国产在线视频| 少妇高潮的动态图| 99九九线精品视频在线观看视频| 日韩欧美一区视频在线观看| 欧美日韩精品成人综合77777| 亚洲欧美中文字幕日韩二区| 欧美最新免费一区二区三区| 久久97久久精品| 久久婷婷青草| av又黄又爽大尺度在线免费看| 欧美成人精品欧美一级黄| 午夜福利在线观看免费完整高清在| 精品少妇久久久久久888优播| 精品视频人人做人人爽| 各种免费的搞黄视频| 亚洲成人一二三区av| 午夜福利网站1000一区二区三区| 中文字幕av电影在线播放| 少妇的逼水好多| 亚洲一级一片aⅴ在线观看| 另类精品久久| 亚洲精品久久久久久婷婷小说| 婷婷色综合www| 日本黄色日本黄色录像| 国产精品女同一区二区软件| 久久精品熟女亚洲av麻豆精品| 搡女人真爽免费视频火全软件| 国产视频内射| 涩涩av久久男人的天堂| 午夜日本视频在线| 一本大道久久a久久精品| 成年人免费黄色播放视频| 国产精品三级大全| 一个人免费看片子| 日本vs欧美在线观看视频| 另类精品久久| 久久99精品国语久久久| 国产成人aa在线观看| 高清午夜精品一区二区三区| av国产久精品久网站免费入址| 日韩成人伦理影院| 久久ye,这里只有精品| 老女人水多毛片| 亚洲熟女精品中文字幕| 热re99久久精品国产66热6| 亚洲精品aⅴ在线观看| 91精品伊人久久大香线蕉| 秋霞伦理黄片| 少妇的逼好多水| 久久久久久久久久久丰满| 亚洲色图 男人天堂 中文字幕 | 在线亚洲精品国产二区图片欧美 | 国产爽快片一区二区三区| av在线观看视频网站免费| 亚洲av免费高清在线观看| 亚洲在久久综合| 亚洲欧洲日产国产| 精品视频人人做人人爽| 久久久国产一区二区| 亚洲人成网站在线播| 欧美xxⅹ黑人| 五月伊人婷婷丁香| 这个男人来自地球电影免费观看 | 自拍欧美九色日韩亚洲蝌蚪91| 下体分泌物呈黄色| 不卡视频在线观看欧美| 成人国语在线视频| av国产久精品久网站免费入址| 自线自在国产av| 国产av码专区亚洲av| 国产一区亚洲一区在线观看| 欧美日韩av久久| 91精品一卡2卡3卡4卡| 18在线观看网站| 亚洲av电影在线观看一区二区三区| 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 精品一区二区免费观看| 一本一本综合久久| 国产午夜精品一二区理论片| 午夜av观看不卡| 一区二区三区免费毛片| 国产精品麻豆人妻色哟哟久久| 边亲边吃奶的免费视频| 国产日韩一区二区三区精品不卡 | 国产高清不卡午夜福利| 91久久精品国产一区二区成人| 熟女av电影| 国产精品三级大全| 国产色婷婷99| 国产精品国产三级国产av玫瑰| 亚洲不卡免费看| 亚洲精品日韩av片在线观看| 成人二区视频| 哪个播放器可以免费观看大片| 国产乱来视频区| 日韩av在线免费看完整版不卡| 亚洲国产色片| 寂寞人妻少妇视频99o| 插阴视频在线观看视频| 人体艺术视频欧美日本| 久久午夜福利片| 性高湖久久久久久久久免费观看| 嫩草影院入口| videosex国产| 亚洲国产最新在线播放| 美女脱内裤让男人舔精品视频| 国产成人a∨麻豆精品| 国产成人精品福利久久| 亚洲成人手机| 九九在线视频观看精品| 人妻人人澡人人爽人人| 在线天堂最新版资源| 丰满少妇做爰视频| 肉色欧美久久久久久久蜜桃| av不卡在线播放| 狠狠婷婷综合久久久久久88av| 好男人视频免费观看在线| 亚洲中文av在线| 亚洲欧美日韩卡通动漫| 成人亚洲精品一区在线观看| 免费观看无遮挡的男女| 免费av不卡在线播放| 国产av国产精品国产| 久久鲁丝午夜福利片| 精品少妇久久久久久888优播| 高清视频免费观看一区二区| 色视频在线一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 制服丝袜香蕉在线| 精品一区二区三区视频在线| 亚洲av欧美aⅴ国产| 一本大道久久a久久精品| 不卡视频在线观看欧美| 亚洲欧洲日产国产| 99久久人妻综合| 亚洲精品乱码久久久久久按摩| 卡戴珊不雅视频在线播放| 国产精品麻豆人妻色哟哟久久| 这个男人来自地球电影免费观看 | 国产视频内射| 国产成人精品在线电影| 激情五月婷婷亚洲| 午夜精品国产一区二区电影| 国产av码专区亚洲av| 久久久精品94久久精品| 久久久久久久大尺度免费视频| 亚洲欧美精品自产自拍| 夫妻午夜视频| 少妇熟女欧美另类| 国产男女超爽视频在线观看| 99久久中文字幕三级久久日本| 人妻系列 视频| 人人妻人人澡人人爽人人夜夜| 一级a做视频免费观看| 丰满饥渴人妻一区二区三| 日本免费在线观看一区| 9色porny在线观看| 新久久久久国产一级毛片| 26uuu在线亚洲综合色| 自线自在国产av| 日韩制服骚丝袜av| 国产一区有黄有色的免费视频| 在线观看三级黄色| 秋霞伦理黄片| 免费大片18禁| 日韩熟女老妇一区二区性免费视频| 99久久综合免费| 美女脱内裤让男人舔精品视频| 我要看黄色一级片免费的| 精品熟女少妇av免费看| av女优亚洲男人天堂| 久久久久久久久久久久大奶| 少妇 在线观看| 欧美日韩成人在线一区二区| 嘟嘟电影网在线观看| 欧美人与性动交α欧美精品济南到 | 国产熟女午夜一区二区三区 | 亚洲国产欧美在线一区| 午夜免费观看性视频| 免费av中文字幕在线| 18+在线观看网站| 色视频在线一区二区三区| 九草在线视频观看| 一本久久精品| 久久久午夜欧美精品| 亚洲人与动物交配视频| 亚洲色图 男人天堂 中文字幕 | 9色porny在线观看| www.av在线官网国产| 王馨瑶露胸无遮挡在线观看| 寂寞人妻少妇视频99o| 日韩 亚洲 欧美在线| 91午夜精品亚洲一区二区三区| 少妇 在线观看| 亚洲美女黄色视频免费看| 中文字幕亚洲精品专区| 桃花免费在线播放| 精品一品国产午夜福利视频| 黄片播放在线免费| .国产精品久久| 2021少妇久久久久久久久久久| 亚洲av欧美aⅴ国产| 国产乱来视频区| 69精品国产乱码久久久| 欧美成人午夜免费资源| xxx大片免费视频| av专区在线播放| 丰满迷人的少妇在线观看| 亚洲精品美女久久av网站| 国产免费一区二区三区四区乱码| 午夜福利网站1000一区二区三区| 国产精品人妻久久久影院| 天天影视国产精品| 久久ye,这里只有精品| 国产成人freesex在线| 好男人视频免费观看在线| 国产成人精品久久久久久| 51国产日韩欧美| 亚洲精品日韩av片在线观看| videos熟女内射| 亚洲综合色网址| 少妇被粗大的猛进出69影院 | 精品久久久精品久久久| av.在线天堂| 婷婷色综合www| 全区人妻精品视频| 一级二级三级毛片免费看| 欧美日韩国产mv在线观看视频| 韩国高清视频一区二区三区| 亚洲美女黄色视频免费看| 国产视频内射| 亚洲精品自拍成人| 水蜜桃什么品种好| 男女啪啪激烈高潮av片| 欧美精品亚洲一区二区| 美女大奶头黄色视频| 老司机影院成人| 夜夜看夜夜爽夜夜摸| 久久韩国三级中文字幕| 欧美日本中文国产一区发布| 国产精品女同一区二区软件| 国产午夜精品一二区理论片| 国产男女内射视频| 欧美97在线视频| 各种免费的搞黄视频| 日日啪夜夜爽| 日韩不卡一区二区三区视频在线| 亚洲av不卡在线观看| av播播在线观看一区| 99热这里只有是精品在线观看| 18在线观看网站| 亚洲欧美日韩另类电影网站| 免费黄网站久久成人精品| 全区人妻精品视频| 婷婷色综合www| 亚洲人成77777在线视频| 日本爱情动作片www.在线观看| 十八禁网站网址无遮挡| 日本av手机在线免费观看| 国产不卡av网站在线观看| 丰满乱子伦码专区| 午夜视频国产福利| 激情五月婷婷亚洲| 日韩在线高清观看一区二区三区| 制服人妻中文乱码| 97超碰精品成人国产| 有码 亚洲区| 亚洲三级黄色毛片| 久久久久久久久大av| 亚洲人成网站在线播| 久久精品国产自在天天线| 特大巨黑吊av在线直播| 欧美人与性动交α欧美精品济南到 | 欧美丝袜亚洲另类| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 大话2 男鬼变身卡| 老司机影院成人| 91久久精品电影网| 日韩 亚洲 欧美在线| 亚洲精品乱码久久久久久按摩| 少妇 在线观看| 免费高清在线观看视频在线观看| 久久精品国产亚洲网站| 91国产中文字幕| 亚洲精品中文字幕在线视频| 亚洲精品色激情综合| 国产成人一区二区在线| 高清欧美精品videossex| 久久久午夜欧美精品| 伊人久久国产一区二区| 制服人妻中文乱码| 高清午夜精品一区二区三区| 欧美丝袜亚洲另类| 日韩三级伦理在线观看| 91精品三级在线观看| 国内精品宾馆在线| 哪个播放器可以免费观看大片| 国产一区有黄有色的免费视频| 久久综合国产亚洲精品| 成年美女黄网站色视频大全免费 | 桃花免费在线播放| 亚洲无线观看免费| 99久久精品国产国产毛片| 欧美精品一区二区免费开放| 国产日韩欧美在线精品| 18禁在线播放成人免费| 在线天堂最新版资源| 国产69精品久久久久777片| 国产成人免费观看mmmm| 尾随美女入室| 2021少妇久久久久久久久久久| 制服诱惑二区| 成人综合一区亚洲| 一二三四中文在线观看免费高清| 男女啪啪激烈高潮av片| 高清在线视频一区二区三区| 99九九线精品视频在线观看视频| 中文字幕人妻丝袜制服| 亚洲精品中文字幕在线视频| 一级片'在线观看视频| av有码第一页| 在线免费观看不下载黄p国产| 日韩不卡一区二区三区视频在线| 热99久久久久精品小说推荐| 国产在视频线精品| 日本色播在线视频| 亚洲三级黄色毛片| 国产精品久久久久久av不卡| 亚洲精品自拍成人| 日本91视频免费播放| 国产成人av激情在线播放 | 久久午夜福利片| av在线观看视频网站免费| tube8黄色片| 91在线精品国自产拍蜜月| 母亲3免费完整高清在线观看 | 精品久久久精品久久久| 欧美精品一区二区免费开放| 日韩视频在线欧美| 毛片一级片免费看久久久久| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av蜜桃| 精品久久蜜臀av无| 国产午夜精品一二区理论片| 久久久久久伊人网av| 中文字幕av电影在线播放| 性高湖久久久久久久久免费观看| 嫩草影院入口| 国产av精品麻豆| 中文字幕精品免费在线观看视频 | 国产精品 国内视频| 国产无遮挡羞羞视频在线观看| 久久久久精品久久久久真实原创| 精品亚洲成a人片在线观看| 色哟哟·www| 十八禁网站网址无遮挡| 校园人妻丝袜中文字幕| 在线精品无人区一区二区三| 亚洲色图综合在线观看| 最近中文字幕2019免费版| 老司机影院成人| a级毛片在线看网站| 成人二区视频| 亚洲综合色惰| 一区二区日韩欧美中文字幕 | 国产成人免费观看mmmm| 日韩精品免费视频一区二区三区 | 黑人欧美特级aaaaaa片| videosex国产| 在线 av 中文字幕| 久久99蜜桃精品久久| 内地一区二区视频在线| 啦啦啦在线观看免费高清www| 日日摸夜夜添夜夜添av毛片| 91精品国产国语对白视频| 欧美日韩一区二区视频在线观看视频在线| 九色亚洲精品在线播放| 欧美日韩视频精品一区|