• <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
    国产极品粉嫩免费观看在线 | 熟女av电影| 欧美 亚洲 国产 日韩一| 亚洲国产成人一精品久久久| 国产在线免费精品| 少妇精品久久久久久久| 一级黄片播放器| 如何舔出高潮| 久久久久网色| 久久精品夜色国产| 亚洲国产成人一精品久久久| 国产精品麻豆人妻色哟哟久久| 国产淫语在线视频| 尾随美女入室| a级片在线免费高清观看视频| 日韩欧美一区视频在线观看| 人妻制服诱惑在线中文字幕| 国产精品不卡视频一区二区| 波野结衣二区三区在线| 十八禁网站网址无遮挡| 在线天堂最新版资源| 久久精品久久精品一区二区三区| 国产精品一区二区在线不卡| 九色成人免费人妻av| 欧美另类一区| 人人妻人人添人人爽欧美一区卜| 丝袜美足系列| 欧美日韩av久久| 丁香六月天网| 国产精品99久久久久久久久| 99精国产麻豆久久婷婷| 亚洲精品乱码久久久久久按摩| 日本午夜av视频| 在线观看三级黄色| 视频中文字幕在线观看| 日韩免费高清中文字幕av| 亚洲精品成人av观看孕妇| 中文字幕亚洲精品专区| 久久国产精品男人的天堂亚洲 | 亚洲第一av免费看| 精品酒店卫生间| 午夜福利,免费看| 大陆偷拍与自拍| 日韩av不卡免费在线播放| 国产极品天堂在线| 国产黄色免费在线视频| 国产一区二区三区av在线| 看非洲黑人一级黄片| 日韩三级伦理在线观看| 精品国产一区二区久久| 亚洲一区二区三区欧美精品| 大陆偷拍与自拍| 日本黄色日本黄色录像| 亚洲精品久久成人aⅴ小说 | 免费大片黄手机在线观看| 国产亚洲av片在线观看秒播厂| 久久久久久久大尺度免费视频| 99热全是精品| 亚洲国产精品成人久久小说| 国产成人a∨麻豆精品| 青青草视频在线视频观看| 一边亲一边摸免费视频| 日日撸夜夜添| 日韩中字成人| 亚洲精品日韩av片在线观看| 各种免费的搞黄视频| 亚洲情色 制服丝袜| 男男h啪啪无遮挡| 一本久久精品| 国产在视频线精品| 水蜜桃什么品种好| 国产一级毛片在线| 美女中出高潮动态图| 日韩精品免费视频一区二区三区 | 国产精品熟女久久久久浪| 日韩一区二区三区影片| 蜜桃国产av成人99| 日本91视频免费播放| 尾随美女入室| 天堂中文最新版在线下载| 国产精品偷伦视频观看了| 最近最新中文字幕免费大全7| 在线观看美女被高潮喷水网站| 一本久久精品| 亚洲成人一二三区av| 日韩视频在线欧美| 国产成人精品福利久久| 成人毛片a级毛片在线播放| 精品久久久久久久久亚洲| 亚洲精品国产色婷婷电影| 亚洲,一卡二卡三卡| 国产老妇伦熟女老妇高清| 纯流量卡能插随身wifi吗| 精品一区二区三卡| 亚洲国产精品999| 欧美日韩综合久久久久久| 亚洲av中文av极速乱| 久久精品国产自在天天线| 国产精品人妻久久久影院| 亚洲,一卡二卡三卡| 国产成人aa在线观看| 精品人妻一区二区三区麻豆| 人人妻人人爽人人添夜夜欢视频| 人体艺术视频欧美日本| 亚洲精品自拍成人| 69精品国产乱码久久久| 日韩av免费高清视频| 美女福利国产在线| 久久99热6这里只有精品| 一个人免费看片子| 女人久久www免费人成看片| 91午夜精品亚洲一区二区三区| 亚洲综合色惰| 人人澡人人妻人| 成年人午夜在线观看视频| 久久精品国产亚洲av天美| 欧美激情国产日韩精品一区| 亚洲精品自拍成人| √禁漫天堂资源中文www| 一区二区三区乱码不卡18| 999精品在线视频| 国产探花极品一区二区| 97超视频在线观看视频| 丰满饥渴人妻一区二区三| 国产成人精品福利久久| 又粗又硬又长又爽又黄的视频| 国产成人aa在线观看| 久久久久久久久久成人| 免费看av在线观看网站| 免费看av在线观看网站| 一本久久精品| 日日撸夜夜添| 国产成人精品婷婷| 亚洲婷婷狠狠爱综合网| 如何舔出高潮| 尾随美女入室| 久久免费观看电影| 色吧在线观看| 国产成人免费无遮挡视频| 99热全是精品| 男人添女人高潮全过程视频| 人体艺术视频欧美日本| 69精品国产乱码久久久| 免费av中文字幕在线| 国产在视频线精品| 丝瓜视频免费看黄片| 国产一区二区三区综合在线观看 | 久久久午夜欧美精品| 一级二级三级毛片免费看| 性色av一级| 岛国毛片在线播放| 国产男女内射视频| 99热这里只有精品一区| 男女啪啪激烈高潮av片| 人人妻人人添人人爽欧美一区卜| 乱码一卡2卡4卡精品| 欧美精品人与动牲交sv欧美| 毛片一级片免费看久久久久| 午夜视频国产福利| 又大又黄又爽视频免费| 国产精品久久久久久精品电影小说| 蜜桃在线观看..| 少妇丰满av| 色婷婷久久久亚洲欧美| 亚洲国产精品999| 亚洲怡红院男人天堂| 色婷婷久久久亚洲欧美| 久久久久国产网址| 精品久久国产蜜桃| 丝瓜视频免费看黄片| 中文天堂在线官网| 91精品国产国语对白视频| 在线观看一区二区三区激情| 制服诱惑二区| 五月天丁香电影| 亚洲无线观看免费| 亚洲国产欧美在线一区| 亚洲av不卡在线观看| 国产成人精品久久久久久| 青春草亚洲视频在线观看| 自线自在国产av| 国产成人精品无人区| 日韩成人av中文字幕在线观看| av黄色大香蕉| 黄片无遮挡物在线观看| 婷婷色综合www| 亚洲精华国产精华液的使用体验| 亚洲一级一片aⅴ在线观看| 免费黄色在线免费观看| 高清午夜精品一区二区三区| 欧美日韩视频高清一区二区三区二| 乱码一卡2卡4卡精品| 日产精品乱码卡一卡2卡三| 91久久精品国产一区二区成人| 国产日韩欧美在线精品| 最新中文字幕久久久久| 亚洲人与动物交配视频| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 91成人精品电影| 国产黄色免费在线视频| 老司机亚洲免费影院| 色5月婷婷丁香| 日日撸夜夜添| 久久久久久久大尺度免费视频| 亚洲成人av在线免费| 国产精品国产三级国产av玫瑰| 国产亚洲欧美精品永久| 91午夜精品亚洲一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| av卡一久久| 中文字幕亚洲精品专区| 久久久精品免费免费高清| 亚洲av日韩在线播放| 一边亲一边摸免费视频| 国产精品嫩草影院av在线观看| 成人综合一区亚洲| 永久免费av网站大全| 久久久国产一区二区| 丝袜脚勾引网站| 久热这里只有精品99| 亚洲精品自拍成人| 国产精品一区二区在线不卡| 性色av一级| 国产探花极品一区二区| 美女主播在线视频| 国产成人免费无遮挡视频| 日日啪夜夜爽| 边亲边吃奶的免费视频| av一本久久久久| 国产精品国产三级国产专区5o| 国产成人精品在线电影| 国产综合精华液| 一级爰片在线观看| 亚洲怡红院男人天堂| 韩国av在线不卡| 夜夜看夜夜爽夜夜摸| 国产精品女同一区二区软件| 午夜免费观看性视频| 九九久久精品国产亚洲av麻豆| 青春草国产在线视频| 你懂的网址亚洲精品在线观看| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 欧美 日韩 精品 国产| 91久久精品电影网| 国内精品宾馆在线| 国产熟女欧美一区二区| 午夜福利视频精品| 十八禁网站网址无遮挡| 99九九线精品视频在线观看视频| 80岁老熟妇乱子伦牲交| 自拍欧美九色日韩亚洲蝌蚪91| 晚上一个人看的免费电影| 久久久久久久久久久久大奶| 国产精品熟女久久久久浪| av有码第一页| 亚洲色图综合在线观看| 尾随美女入室| 国产高清国产精品国产三级| 精品一品国产午夜福利视频| 精品久久蜜臀av无| 午夜激情av网站| 国产色婷婷99| 99久久综合免费| 飞空精品影院首页| 久久这里有精品视频免费| 日韩av在线免费看完整版不卡| 中文字幕最新亚洲高清| 另类精品久久| 九草在线视频观看| 日本vs欧美在线观看视频| 一本大道久久a久久精品| 成年av动漫网址| 久久女婷五月综合色啪小说| 国产欧美另类精品又又久久亚洲欧美| 亚洲欧美日韩另类电影网站| 人人妻人人澡人人爽人人夜夜| 大片电影免费在线观看免费| 亚洲五月色婷婷综合| 中国美白少妇内射xxxbb| 国产欧美日韩综合在线一区二区| 亚洲人与动物交配视频| 22中文网久久字幕| 成年人午夜在线观看视频| 99精国产麻豆久久婷婷| 午夜91福利影院| 久久精品国产亚洲av涩爱| 男女高潮啪啪啪动态图| 日韩欧美一区视频在线观看| 各种免费的搞黄视频| 丝瓜视频免费看黄片| 国产精品一国产av| 国产精品国产三级国产专区5o| 成人毛片60女人毛片免费| 国模一区二区三区四区视频| 老司机影院毛片| 日本爱情动作片www.在线观看| 80岁老熟妇乱子伦牲交| 国国产精品蜜臀av免费| av专区在线播放| 黑人高潮一二区| 国产乱人偷精品视频| 亚洲天堂av无毛| 免费黄色在线免费观看| 欧美精品一区二区大全| 国产精品一区www在线观看| 久久人妻熟女aⅴ| 大片免费播放器 马上看| 国产成人freesex在线| 久久久久人妻精品一区果冻| 18+在线观看网站| 日韩中文字幕视频在线看片| 高清欧美精品videossex| 成人影院久久| 亚洲欧美中文字幕日韩二区| 国产成人一区二区在线| 成人黄色视频免费在线看| 亚洲婷婷狠狠爱综合网| 国产白丝娇喘喷水9色精品| 日韩成人av中文字幕在线观看| 亚洲精华国产精华液的使用体验| 中国国产av一级| 亚洲人成网站在线观看播放| 国产精品欧美亚洲77777| 亚洲av成人精品一区久久| 女人精品久久久久毛片| 丝袜喷水一区| 日韩中文字幕视频在线看片| 日韩视频在线欧美| 一级a做视频免费观看| 中文字幕久久专区| 国产精品人妻久久久影院| 午夜福利视频在线观看免费| 免费不卡的大黄色大毛片视频在线观看| 丝瓜视频免费看黄片| 精品国产乱码久久久久久小说| 婷婷成人精品国产| 欧美日韩精品成人综合77777| 超碰97精品在线观看| 欧美亚洲 丝袜 人妻 在线| 人妻少妇偷人精品九色| 天天躁夜夜躁狠狠久久av| 中文乱码字字幕精品一区二区三区| 国产av码专区亚洲av| 成人无遮挡网站| 国产精品一区二区在线观看99| 高清不卡的av网站| 国产成人精品在线电影| 人妻人人澡人人爽人人| 日本-黄色视频高清免费观看| 日韩免费高清中文字幕av| 插逼视频在线观看| 春色校园在线视频观看| 国产成人精品久久久久久| 国产亚洲精品久久久com| 狠狠婷婷综合久久久久久88av| 插阴视频在线观看视频| 熟女电影av网| 日韩亚洲欧美综合| 欧美精品一区二区大全| 一级黄片播放器| 2022亚洲国产成人精品| 热re99久久精品国产66热6| 午夜日本视频在线| 另类精品久久| 国产欧美日韩综合在线一区二区| 99久久人妻综合| 亚洲美女视频黄频| videos熟女内射| 七月丁香在线播放| 日韩视频在线欧美| 久久久精品94久久精品| 中国美白少妇内射xxxbb| 9色porny在线观看| 九草在线视频观看| av国产久精品久网站免费入址| 国产精品一区二区在线观看99| 精品卡一卡二卡四卡免费| 成人亚洲欧美一区二区av| 大片免费播放器 马上看| 黑丝袜美女国产一区| 王馨瑶露胸无遮挡在线观看| 99视频精品全部免费 在线| 午夜91福利影院| 插逼视频在线观看| 男人爽女人下面视频在线观看| 成人毛片a级毛片在线播放| 国产成人aa在线观看| 亚洲精品日韩在线中文字幕| 亚洲精品中文字幕在线视频| 国产一区二区在线观看日韩| av黄色大香蕉| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 久久综合国产亚洲精品| 成人免费观看视频高清| 亚洲欧美清纯卡通| 人成视频在线观看免费观看| 少妇人妻精品综合一区二区| 日本av免费视频播放| 女的被弄到高潮叫床怎么办| 亚洲综合色惰| 免费人成在线观看视频色| 高清av免费在线| 午夜免费观看性视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品久久午夜乱码| 精品亚洲成国产av| 男女啪啪激烈高潮av片| 国产在视频线精品| 免费大片黄手机在线观看| 亚洲精品国产色婷婷电影| 亚洲国产毛片av蜜桃av| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 最近中文字幕2019免费版| 精品少妇久久久久久888优播| 美女中出高潮动态图| 精品一区二区免费观看| 国产日韩欧美视频二区| 久久人妻熟女aⅴ| 精品少妇内射三级| 久久亚洲国产成人精品v| 精品久久国产蜜桃| 久久久久久久久久成人| 美女内射精品一级片tv| 亚洲第一区二区三区不卡| 久久久久网色| 国产精品偷伦视频观看了| 中文字幕亚洲精品专区| 亚洲精品乱久久久久久| 亚洲人成网站在线观看播放| av国产精品久久久久影院| 女性生殖器流出的白浆| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 999精品在线视频| 搡女人真爽免费视频火全软件| 精品亚洲成a人片在线观看| 又粗又硬又长又爽又黄的视频| 欧美变态另类bdsm刘玥| 国产av国产精品国产| 黄色配什么色好看| 在线观看国产h片| 欧美成人午夜免费资源| 狠狠精品人妻久久久久久综合| 欧美三级亚洲精品| 日本av手机在线免费观看| 人妻 亚洲 视频| 国产成人精品婷婷| 午夜福利视频在线观看免费| 久久精品久久久久久噜噜老黄| 日本黄色片子视频| 欧美精品人与动牲交sv欧美| 99久国产av精品国产电影| 国产高清不卡午夜福利| 欧美日韩综合久久久久久| 一级毛片我不卡| 午夜福利网站1000一区二区三区| 亚洲精品第二区| 欧美日韩精品成人综合77777| videos熟女内射| 多毛熟女@视频| 亚洲性久久影院| 美女内射精品一级片tv| 高清黄色对白视频在线免费看| 麻豆乱淫一区二区| 日本色播在线视频| 精品亚洲成a人片在线观看| 国产亚洲午夜精品一区二区久久| 亚洲精华国产精华液的使用体验| 中文字幕av电影在线播放| 国产熟女欧美一区二区| 亚洲人与动物交配视频| 日本午夜av视频| 久久久a久久爽久久v久久| av一本久久久久| 国产成人精品福利久久| av不卡在线播放| 欧美日韩成人在线一区二区| 欧美成人午夜免费资源| 国产成人精品婷婷| 赤兔流量卡办理| 伦理电影免费视频| 久久久精品免费免费高清| 欧美老熟妇乱子伦牲交| xxx大片免费视频| 乱码一卡2卡4卡精品| 少妇丰满av| 久久精品久久久久久久性| 中文天堂在线官网| 久久综合国产亚洲精品| 99久久精品一区二区三区| 一级片'在线观看视频| 亚洲国产最新在线播放| av一本久久久久| 久久狼人影院| 免费播放大片免费观看视频在线观看| 国产精品蜜桃在线观看| 亚洲精品,欧美精品| 老司机亚洲免费影院| 免费看光身美女| 欧美日韩视频精品一区| 男女国产视频网站| 久久精品国产亚洲av天美| 精品久久久久久久久亚洲| 亚洲欧洲日产国产| 亚洲五月色婷婷综合| 亚洲内射少妇av| 国产精品国产三级国产专区5o| 尾随美女入室| 精品久久久久久久久亚洲| 毛片一级片免费看久久久久| 天天影视国产精品| 欧美日韩在线观看h| 黑人欧美特级aaaaaa片| 久热久热在线精品观看| 成人亚洲精品一区在线观看| 亚洲欧洲精品一区二区精品久久久 | 啦啦啦中文免费视频观看日本| 久久久久久久久久久丰满| 大陆偷拍与自拍| 五月伊人婷婷丁香| 国产综合精华液| 丝瓜视频免费看黄片| 欧美人与善性xxx| 99久久精品一区二区三区| 久久久久久久久久成人| 精品国产乱码久久久久久小说| 日韩免费高清中文字幕av| 22中文网久久字幕| 亚洲国产欧美日韩在线播放| 菩萨蛮人人尽说江南好唐韦庄| 久久精品熟女亚洲av麻豆精品| 国产日韩一区二区三区精品不卡 | 国产精品一区二区在线不卡| 五月开心婷婷网| 国产精品嫩草影院av在线观看| 2018国产大陆天天弄谢| 777米奇影视久久| 五月天丁香电影| 色5月婷婷丁香| 黄色怎么调成土黄色| 国产av精品麻豆| 午夜激情福利司机影院| 亚洲欧洲国产日韩| av专区在线播放| 国产成人a∨麻豆精品| 人成视频在线观看免费观看| 999精品在线视频| www.色视频.com| 考比视频在线观看| 美女福利国产在线| 高清黄色对白视频在线免费看| 国产无遮挡羞羞视频在线观看| 日日撸夜夜添| 在线播放无遮挡| 丰满少妇做爰视频| 交换朋友夫妻互换小说| 18禁动态无遮挡网站| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲在久久综合| 热99国产精品久久久久久7| 欧美精品一区二区免费开放| 久久久国产欧美日韩av| 国产日韩一区二区三区精品不卡 | 亚洲国产av影院在线观看| 99久久精品一区二区三区| 精品少妇内射三级| 男女啪啪激烈高潮av片| 亚洲不卡免费看| 男女高潮啪啪啪动态图| .国产精品久久| 亚洲精品第二区| www.色视频.com| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美女主播在线视频| 秋霞伦理黄片| 丝袜美足系列| 久久99精品国语久久久| 久久久国产欧美日韩av| 国产欧美日韩综合在线一区二区| 国产视频内射| 两个人的视频大全免费| 日本wwww免费看| av在线播放精品| 国产深夜福利视频在线观看| 人妻夜夜爽99麻豆av| 99九九线精品视频在线观看视频| 国产精品蜜桃在线观看| 全区人妻精品视频| 丰满迷人的少妇在线观看| 日韩欧美一区视频在线观看| 草草在线视频免费看| 青春草国产在线视频| 高清黄色对白视频在线免费看| 性色avwww在线观看| 性高湖久久久久久久久免费观看| 9色porny在线观看| 蜜桃在线观看..| 国产精品欧美亚洲77777| 一区二区日韩欧美中文字幕 | 午夜福利,免费看| 丰满迷人的少妇在线观看| 午夜福利,免费看| 人人妻人人添人人爽欧美一区卜| a级片在线免费高清观看视频| 亚洲av成人精品一二三区| 国产精品麻豆人妻色哟哟久久| 国产成人av激情在线播放 | 午夜91福利影院| 国产一区二区三区av在线| 国产免费又黄又爽又色|