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

    中國耕地多功能耦合協(xié)調時空演變及其驅動因素

    2023-06-12 05:23:04代亞強陳媛媛柯新利
    農(nóng)業(yè)工程學報 2023年7期
    關鍵詞:耕地耦合驅動

    張 玥,代亞強,陳媛媛,柯新利

    ?土地保障與生態(tài)安全?

    中國耕地多功能耦合協(xié)調時空演變及其驅動因素

    張 玥,代亞強,陳媛媛,柯新利※

    (華中農(nóng)業(yè)大學公共管理學院,武漢 430070)

    為探究耕地多功能耦合協(xié)調演變規(guī)律及其驅動因素,該研究依據(jù)2000-2018年中國30個省級行政區(qū)統(tǒng)計數(shù)據(jù),綜合評價耕地“生產(chǎn)-生活-生態(tài)”三維功能水平并采用核密度估計和可視化制圖刻畫耕地多功能時空特征,進而采用耦合協(xié)調度模型分析耕地多功能耦合協(xié)調關系,最終采用地理探測器識別耕地多功能耦合協(xié)調演變的驅動因素并解析其驅動機制。研究結果表明:1)耕地生產(chǎn)功能水平、生活功能水平和生態(tài)功能水平分別處于[0.030, 0.608]、[0.042, 0.672]和[0.058, 0.897]范圍,生產(chǎn)功能呈現(xiàn)上升-下降-上升波動變化,生活功能水平先上升后下降,生態(tài)功能水平則整體小幅提升,各項功能均具有極化特征。耕地生產(chǎn)功能、生活功能和生態(tài)功能空間格局分別呈現(xiàn)“東高西低”“北高南低”和“西高東低”分布特征。2)耕地多功能耦合協(xié)調水平處于[0.093, 0.554]范圍,整體穩(wěn)步提升但水平仍然較低,多數(shù)省份耕地多功能仍為失調狀態(tài)。耕地多功能耦合協(xié)調空間格局逐漸呈現(xiàn)集聚提升的演化特征,東北地區(qū)、京津冀地區(qū)和長三角地區(qū)耕地多功能耦合協(xié)調水平較高,西北地區(qū)則處于較低水平。3)耕地多功能耦合協(xié)調度演變受經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平影響顯著,其中,經(jīng)營主體條件和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調度演變的驅動力較強,農(nóng)業(yè)發(fā)展水平和工業(yè)發(fā)展水平次之。各驅動因素交互作用類型均為雙因子增強或非線性增強,因素交互作用正向強化了各驅動因素對耕地多功能耦合協(xié)調度演變的驅動力。研究結果為充實豐富耕地多功能研究提供了實證指導,同時為推進耕地資源差異化治理和耕地多功能互促提升提供了決策支撐。

    土地利用;耕地;多功能;耦合協(xié)調度;時空演變;驅動因素

    0 引 言

    耕地資源是人類賴以生存和發(fā)展的基礎自然資源和重要空間載體,其功能衍生與需求層次升級和社會經(jīng)濟革新密切相關[1]。伴隨中國由農(nóng)耕文明過渡為工業(yè)文明進而走向生態(tài)文明,人的需求也經(jīng)歷著“生存-享受-發(fā)展”的漸次演進,耕地多功能已然成為耕地利用的合理方向。從家庭聯(lián)產(chǎn)承包責任制改革著力解決人民溫飽問題,到改革開放推進工業(yè)化和城鎮(zhèn)化進程,再到鄉(xiāng)村振興和生態(tài)文明建設提升人類福祉、促進可持續(xù)發(fā)展,耕地利用承擔著促進經(jīng)濟發(fā)展、維護社會穩(wěn)定和保障生態(tài)安全等多重任務[2-3]。與此同時,耕地資源的稀缺性日益凸顯,糧食安全、生態(tài)保護與城鎮(zhèn)發(fā)展的用地矛盾和空間競爭愈發(fā)激烈,嚴重催化了耕地資源隱性流失,抑制了耕地生產(chǎn)潛能發(fā)揮,同時削弱了其社會貢獻價值并引發(fā)生態(tài)系統(tǒng)退化[4-5],有違耕地保護底線原則和永續(xù)利用科學理念。當前,傳統(tǒng)的、單一的耕地生產(chǎn)功能管理模式已難以適配多元化的耕地利用現(xiàn)實形勢,同時因其忽視耕地利用復合價值、抑制微觀主體保護動力而未能充分助力耕地保護實踐。因此,在中國新發(fā)展階段與新發(fā)展格局背景下,耕地多功能協(xié)調利用不僅是順應社會經(jīng)濟發(fā)展、滿足人民美好需求的必然選擇,同時是平衡多方利益關系、強化耕地保護實效的重要手段。對此,積極開展耕地多功能實證研究與機制探討,尤其是探清耕地多功能耦合協(xié)調規(guī)律表征與現(xiàn)實癥結,對提升耕地利用綜合效能、促進耕地管理創(chuàng)新轉型具有重要的科學意義和現(xiàn)實意義。

    長期以來,耕地多功能研究已逐漸形成“內涵界定-現(xiàn)狀評價-管理調控”的研究脈絡,并在宏觀-中觀-微觀不同研究尺度上加以表達。在內涵界定方面,多功能概念起源于農(nóng)業(yè)領域,并逐漸應用至土地利用、景觀管理、生態(tài)系統(tǒng)服務等領域[6-8]。耕地系統(tǒng)是自然、經(jīng)濟、社會和生態(tài)子系統(tǒng)高度耦合的復合系統(tǒng),耕地多功能是指在復雜因素交互作用下耕地提供產(chǎn)品和服務以滿足人類需求的各項能力[9-10]。依據(jù)耕地多功能內涵對其類型加以細分,耕地多功能不僅包括提供糧食、蔬菜等重要農(nóng)作物的生產(chǎn)能力,同時包括保障農(nóng)民生活、維持農(nóng)村穩(wěn)定等的生活功能和提供氣候調節(jié)、水源涵養(yǎng)、環(huán)境凈化等生態(tài)系統(tǒng)服務的生態(tài)能力,兼具阻隔城鎮(zhèn)空間、承載鄉(xiāng)土文化等其他衍生功能[11-12]。在現(xiàn)狀評價方面,眾多學者多基于耕地多功能內涵界定及類型劃分構建評價體系,采用各類綜合評價方法、三角模型、空間收斂性法等從多角度揭示耕地多功能時空演變特征及規(guī)律[13-16]。在管理調控方面,基于不同研究單元的耕地多功能評價結果,耕地多功能管理可綜合運用分區(qū)調控、重點保護、典型模式推廣等手段進行合理引導和科學調控[2,17]?,F(xiàn)有耕地多功能研究已初步構建從內涵認知到管理應用的研究體系,為本文提供了重要的理論基礎和實踐支撐。耕地利用系統(tǒng)是經(jīng)濟-社會-生態(tài)子系統(tǒng)構成的復合系統(tǒng),耕地多功能共存于同一區(qū)域實體且具有此消彼長的內在關聯(lián),因而耕地多功能研究需充分關注并著力厘清功能結構性、共生性和協(xié)同性特征表達及其演化邏輯。然而,現(xiàn)有耕地多功能研究多關注功能時空分異特征和權衡協(xié)同關系,并多在區(qū)域、省級、市級尺度開展實證研究,因而其未能從全局視角充分研判耕地多功能耦合協(xié)調總體水平及其階段性規(guī)律,而且研究結論和實踐建議具有明顯的地域性、典型性特征。耕地多功能耦合協(xié)調水平是衡量耕地子功能交互影響程度的重要指標,其時序更替性和空間特定性可反映耕地利用綜合效能的實現(xiàn)程度。中國不同地區(qū)自然資源稟賦、耕地利用模式、社會經(jīng)濟水平等各具特點,耕地利用主導功能及其與其他衍生功能關系不盡相同,針對國家層面省域尺度的宏觀研究有待開展。同時,鮮少研究定量識別耕地多功能耦合協(xié)調的驅動因素,尤其是未能全面考察耕地利用系統(tǒng)“內生-外緣”作用力以深入剖析耕地多功能演進的內在機理,難以為耕地多功能優(yōu)化管控和統(tǒng)籌協(xié)調提供充分的決策支撐。耕地利用系統(tǒng)是受社會經(jīng)濟環(huán)境影響的復雜巨系統(tǒng),要素間相互關聯(lián)、彼此制約的關系推動著各子功能演化重塑,影響著耕地多功能發(fā)揮程度和協(xié)調穩(wěn)態(tài)。因此,診斷耕地多功能耦合協(xié)調演變“內生-外緣”驅動因素,深入解析各因素作用機制,可為促進耕地功能協(xié)同融合提供科學依據(jù)。

    鑒于此,本文構建“生產(chǎn)-生活-生態(tài)”三維耕地多功能評價體系,采用2000-2018年中國30個省級行政區(qū)面板數(shù)據(jù)開展耕地多功能實證探究,刻畫并分析耕地多功能時序特征和空間格局,并采用地理探測器識別耕地多功能耦合協(xié)調度時空演變的驅動因素,以期為制定耕地優(yōu)化配置方案、實施耕地資源差異化治理提供科學依據(jù)和政策建議。

    1 研究方法與數(shù)據(jù)來源

    1.1 研究方法

    本研究構建“生產(chǎn)-生活-生態(tài)”三維耕地多功能評價模型,綜合評價2000-2018年中國30個省級行政區(qū)耕地多功能水平,采用核密度估計和可視化制圖刻畫耕地多功能時空演進特征,進而采用耦合協(xié)調度模型分析耕地多功能耦合協(xié)調關系及其演變規(guī)律,最終采用地理探測器識別耕地多功能耦合協(xié)調度演變的驅動因素(圖1)。

    1.1.1 耕地多功能評價模型

    耕地系統(tǒng)是自然-人工的復合系統(tǒng),具有經(jīng)濟、社會、生態(tài)等多重屬性,其功能因人類利用活動而存在。中國耕地利用和保護的重心在于保障糧食供給、維持社會穩(wěn)定和維護生態(tài)安全3個維度,分別對應耕地的生產(chǎn)、生活和生態(tài)3種功能。結合現(xiàn)有研究成果[2,17-18],依據(jù)系統(tǒng)性、科學性和數(shù)據(jù)可獲取性等原則,本文選取12項評價指標綜合評價耕地多功能,指標體系及具體釋義如表1所示。

    1)生產(chǎn)功能。生產(chǎn)功能是耕地最基本、最核心的功能,反映耕地生產(chǎn)活動提供農(nóng)產(chǎn)品的能力[18]。本文綜合考慮耕地利用的實際產(chǎn)出與耕作狀況,選取糧食單產(chǎn)水平、耕地地均產(chǎn)值、復種指數(shù)和土地墾殖率4項指標予以衡量。其中,糧食單產(chǎn)水平與耕地地均產(chǎn)值分別表征實物產(chǎn)量水平和經(jīng)濟產(chǎn)出水平以反映耕地利用的實際產(chǎn)出,復種指數(shù)和土地墾殖率則分別表征農(nóng)業(yè)生產(chǎn)活動強度和耕地開發(fā)更新程度以反映耕地利用的耕作狀況。

    2)生活功能。生活功能是耕地最主要的衍生功能之一,反映耕地保障糧食安全和農(nóng)民生活的能力[19]。為體現(xiàn)耕地利用滿足城鄉(xiāng)居民糧食需求和保障糧食穩(wěn)定自給的糧食安全保障功能,本文選取人均糧食保障率予以反映,具體參照國際慣例,以每年人均400 kg糧食為基準進行測算[2]。為體現(xiàn)耕地利用承載農(nóng)村剩余勞動力和滿足促進農(nóng)村社會經(jīng)濟發(fā)展需求的農(nóng)民生活保障功能,本文選取家庭農(nóng)業(yè)收入比重、農(nóng)業(yè)從業(yè)人員比重和人均農(nóng)業(yè)機械化水平3項指標予以反映。為與下文農(nóng)業(yè)發(fā)展水平驅動因素中的第一產(chǎn)業(yè)從業(yè)人員數(shù)進行區(qū)分,此處采用從業(yè)人員比重值衡量耕地生活功能的勞動人員承載能力。

    3)生態(tài)功能。生態(tài)功能反映耕地維持生物多樣性、提供生態(tài)系統(tǒng)服務和維護農(nóng)田生態(tài)健康的能力,不僅與耕地自身生態(tài)本底有關,同時受農(nóng)業(yè)生產(chǎn)方式的影響[20-21]。本文選取農(nóng)田生態(tài)系統(tǒng)多樣性和耕地占生態(tài)用地比例2項指標反映耕地的生態(tài)本底條件[4,18,22],其中,農(nóng)田生態(tài)系統(tǒng)多樣性通過各品種作物播種面積與農(nóng)作物播種面積之比加以測算。同時,本文選取人均耕地生態(tài)承載力和耕地利用化學負荷反映耕地的生態(tài)系統(tǒng)服務供給和農(nóng)田生態(tài)環(huán)境狀況。

    表1 耕地多功能評價指標體系

    依據(jù)耕地多功能評價指標體系,采用極差標準化法對原始數(shù)據(jù)進行標準化處理,采用熵值法確定指標權重,最后采用線性加權法計算得到各省耕地多功能水平。計算方法如下:

    式中為第個省份第項評價指標標準化數(shù)值;為第項評價指標權重;為第個省份的耕地多功能綜合分值,取值區(qū)間為[0,1],數(shù)值越大,則耕地多功能水平越高。

    1.1.2 核密度估計

    核密度估計(kernel density estimation,KDE)主要優(yōu)勢在于不依賴數(shù)據(jù)分布的先驗知識,不設置參數(shù)模型的基本假設,從數(shù)據(jù)本身出發(fā)用連續(xù)的密度函數(shù)描述變量階段分布特征和時序演進規(guī)律[23-24]。因此,借助Eviews8.0軟件,本文采用核密度估計方法刻畫耕地多功能階段性、動態(tài)性特征,具體計算方法可參考文獻[23]。

    1.1.3 耦合協(xié)調度模型

    耦合協(xié)調度常用以表征兩個及以上子系統(tǒng)在發(fā)展演化過程中的協(xié)調一致程度[25]。本文構建耕地多功能耦合協(xié)調度模型,用以定量評價耕地生產(chǎn)功能、生活功能和生態(tài)功能的交互作用。在耦合協(xié)調度模型中,耦合度表征子系統(tǒng)之間的相互作用強弱程度;協(xié)調度表征子系統(tǒng)相互作用關系中良性互動強弱程度[26-27];最終,耦合協(xié)調度表征子系統(tǒng)耦合交互的協(xié)同程度。耦合協(xié)調度具體計算方法及類型劃分依據(jù)可參考文獻[28-29]。

    1.1.4 地理探測器

    耕地多功能耦合協(xié)調度時空演變受內生驅動和外在驅動的雙重影響,識別耕地多功能耦合協(xié)調度時空演變的驅動因素對實現(xiàn)耕地充分利用和有效保護具有重要的參考意義。依據(jù)現(xiàn)有研究[4,18,30-32],本文從經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平4個方面考察各因素對耕地多功能耦合協(xié)調度時空演變的影響程度(表2)。依據(jù)代表性和數(shù)據(jù)可獲取性原則,在經(jīng)營主體條件方面,選取農(nóng)村居民人均純收入和農(nóng)村居民消費水平作為表征指標,反映耕地綜合利用理念的認知程度和功能實現(xiàn)程度的保障條件[33];在農(nóng)業(yè)發(fā)展水平方面,選取第一產(chǎn)業(yè)從業(yè)人員數(shù)和第一產(chǎn)業(yè)增加值占比作為表征指標,反映農(nóng)業(yè)發(fā)展集約化、現(xiàn)代化水平。為與上文耕地多功能評價指標體系中的農(nóng)業(yè)從業(yè)人員比重進行區(qū)分,此處采用第一產(chǎn)業(yè)從業(yè)人員數(shù)和增加值占比共同反映農(nóng)業(yè)發(fā)展水平驅動作用;在工業(yè)發(fā)展水平方面,選取第二產(chǎn)業(yè)增加值占比和規(guī)模以上工企業(yè)個數(shù)作為表征指標,反映工業(yè)發(fā)展現(xiàn)狀對耕地多功能發(fā)揮的經(jīng)濟支撐和需求引導;在城鎮(zhèn)化水平方面,選取城鎮(zhèn)人口占比和城鎮(zhèn)道路面積作為表征指標,反映城鎮(zhèn)化發(fā)展對耕地資源本底條件及其功能水平的外力改造和強度加載。

    表2 耕地多功能耦合協(xié)調度驅動因素指標選取

    地理探測器(geographical detector)是探析地理現(xiàn)象空間分異形成機理的重要方法,其優(yōu)勢在于不設置線性假設和條件限制,可客觀探測定性因子和定量因子的驅動作用及不同因子間的交互作用[34]。因此,本文首先運用自然斷點法對驅動因素進行離散化和類別化處理,進而采用地理探測器定量識別各驅動因素作用程度及因素交互作用。地理探測器具體計算過程和驅動因素交互作用類型判別依據(jù)可參考文獻[34]。

    1.2 研究區(qū)域與數(shù)據(jù)來源

    由于西藏、香港、澳門及臺灣數(shù)據(jù)缺失,本文選取中國30個省級行政區(qū)為研究單元,其中包括22個省、4個自治區(qū)、4個直轄市。本文使用的數(shù)據(jù)包括統(tǒng)計數(shù)據(jù)與行政區(qū)劃數(shù)據(jù)。其中,統(tǒng)計數(shù)據(jù)均來源于2001-2019年《中國統(tǒng)計年鑒》《中國農(nóng)村統(tǒng)計年鑒》《中國環(huán)境統(tǒng)計年鑒》、中國環(huán)境數(shù)據(jù)庫。因上述指標相鄰年份數(shù)據(jù)具有關聯(lián)性、較少發(fā)生異常變動,因此個別年份缺失數(shù)據(jù)采用趨勢外推或插值法進行近似獲取[35]。行政區(qū)劃數(shù)據(jù)來源于中科院地理科學與資源研究所資源環(huán)境科學與數(shù)據(jù)中心(http:// www. resdc.cn/)。

    2 結果與分析

    2.1 耕地多功能評價

    2.1.1 時序演變特征

    依據(jù)耕地各功能水平評價結果,繪制2000-2018年耕地生產(chǎn)、生活和生態(tài)功能水平核密度曲線圖(圖2)。1)從耕地生產(chǎn)功能核密度曲線來看,曲線重心經(jīng)歷“右-左-右”遷移軌跡,曲線始終為單一波峰且高度逐漸下降,表明研究期內中國耕地生產(chǎn)功能水平整體呈現(xiàn)上升-下降-上升的波動趨勢,各省生產(chǎn)功能處于極化狀態(tài)且差異有所擴大;2)從耕地生活功能核密度曲線來看,2000-2015年重心持續(xù)向右遷移,2015-2018年重心輕微向左遷移,曲線保持單一波峰且高度持續(xù)下降,表明研究期內中國耕地生活功能水平整體呈現(xiàn)先上升后下降的趨勢,各省生活功能呈現(xiàn)兩極分化狀態(tài)且差異逐漸擴大;3)從耕地生態(tài)功能核密度曲線來看,曲線重心輕微向右遷移,曲線單一波峰于2000-2015年持續(xù)下降,后于2015-2018年上升,表明中國耕地生態(tài)功能水平整體呈現(xiàn)小幅提升,省域間耕地生態(tài)功能水平差異呈現(xiàn)先擴大后縮小的特征。

    圖2 2000-2018年中國耕地多功能時序特征

    2.1.2 空間格局特征

    為直觀展現(xiàn)各功能空間分布,繪制耕地生產(chǎn)功能、生活功能和生態(tài)功能空間格局圖。從耕地生產(chǎn)功能來看,研究期內中國耕地生產(chǎn)功能水平處于[0.030, 0.608]范圍,高水平區(qū)集中位于黃淮海平原并逐漸向東北平原和長江中下游地區(qū)延伸,低水平區(qū)則分散位于西北地區(qū)和西南地區(qū),大體呈現(xiàn)“東高西低”的空間非均衡格局(圖3)。具體來看:1)2000年,耕地生產(chǎn)功能高水平區(qū)集中在黃淮海平原,得益于其平坦廣闊的地形和較為完善的農(nóng)業(yè)設施;2)2005-2015年,耕地生產(chǎn)功能高水平區(qū)集中在黃淮海平原和長江中下游地區(qū),四川盆地耕地生產(chǎn)功能有所下降而云貴地區(qū)、華南部分地區(qū)則有所提升;3)2018年,耕地生產(chǎn)功能高水平區(qū)集中位于東北平原、黃淮海平原和長江中下游地區(qū),低水平區(qū)位于西北地區(qū)和西南地區(qū),耕地生產(chǎn)功能受限于崎嶇地形、水熱條件較差、農(nóng)業(yè)設施等條件[36-37]。總體而言,研究期內耕地生產(chǎn)功能空間格局以黃淮海平原、長江中下游地區(qū)等重要農(nóng)業(yè)區(qū)為核心輻射延展,耕地生產(chǎn)功能整體水平有所提升。

    從耕地生活功能來看,研究期內中國耕地生活功能水平處于[0.042, 0.672]范圍,高水平區(qū)分散位于黃淮海平原、東北平原和新疆地區(qū),隨后延伸至江漢平原和長江中下游地區(qū),低水平區(qū)則位于京津地區(qū)、云貴地區(qū)和華南地區(qū),大體呈現(xiàn)“北高南低”的空間非均衡格局(圖4)。1)2000-2005年,耕地生活功能高水平省份為新疆、黑龍江和江蘇,新疆產(chǎn)業(yè)結構以第一產(chǎn)業(yè)為主,其余兩省為重要的糧食主產(chǎn)區(qū),耕地資源可為農(nóng)民生活提供較大保障;2)2010-2015年,耕地生活功能高水平區(qū)以上述三省為核心延伸至北方干旱半干旱區(qū)、東北平原、黃淮海地區(qū)和長江中下游北部地區(qū);3)2018年,耕地生活功能較高水平區(qū)向長江中下游南部地區(qū)拓展,該空間格局與耕地生產(chǎn)功能高水平區(qū)布局具有相似性。研究期內,耕地生活功能低水平區(qū)分散位于云貴地區(qū)、華南地區(qū)和西北三省。其中,華南地區(qū)經(jīng)濟發(fā)達、農(nóng)民生計選擇多樣,云貴地區(qū)和西北三省人口外遷嚴重,耕地生活保障作用較為薄弱。

    從耕地生態(tài)功能來看,研究期內中國耕地生態(tài)功能水平處于[0.058, 0.897]范圍,高水平區(qū)集中位于東北地區(qū)和西南地區(qū)并逐漸延伸至黃淮海平原和四川盆地及周邊地區(qū),低水平區(qū)則位于京津冀地區(qū)和長江中部地區(qū),大體呈現(xiàn)“西高東低”的空間格局(圖5)。具體來看:1)2000-2005年,耕地生態(tài)功能高水平區(qū)集中位于東北地區(qū)和西南地區(qū),其中,東北地區(qū)耕地資源豐富、規(guī)模經(jīng)營較為成熟,西南地區(qū)化肥、農(nóng)藥等農(nóng)資投放量較少。同時,內蒙古人均耕地資源較多且規(guī)模經(jīng)營范圍較廣,因而耕地生態(tài)功能穩(wěn)定保持高水平。耕地生態(tài)功能低水平區(qū)則位于京津冀地區(qū)、長三角地區(qū)和長江中部地區(qū),上述地區(qū)耕地資源約束趨緊且利用強度增大,因而其耕地生態(tài)系統(tǒng)愈加脆弱、生態(tài)功能逐漸退化;2)2010-2015年,耕地生態(tài)功能高水平區(qū)向黃淮海平原和四川盆地及周邊地區(qū)拓展,同時在生態(tài)文明理念倡導與經(jīng)濟發(fā)展基礎保障的雙重作用下,低水平區(qū)生態(tài)逐漸改善;3)2018年,耕地生態(tài)功能仍然保持“西高東低”的格局特征。

    注:西藏、香港、澳門、臺灣數(shù)據(jù)暫缺。下同。

    圖4 2000-2018年中國耕地生活功能空間格局

    2.2 耕地多功能耦合協(xié)調度評價

    從耕地多功能耦合協(xié)調時序演變來看,研究期內中國耕地多功能耦合協(xié)調水平處于[0.093, 0.554]范圍,整體呈現(xiàn)穩(wěn)步提升趨勢但水平仍然較低,多數(shù)省份耕地多功能為失調狀態(tài),表明中國耕地多功能正在由無序沖突轉為有序發(fā)展,但耕地多功能耦合協(xié)調水平仍存在較大提升空間。從耕地多功能耦合協(xié)調空間格局來看,研究期內中國耕地多功能耦合協(xié)調度逐漸呈現(xiàn)集聚提升的演化特征(圖6)。具體來看:1)2000年,除內蒙古、黑龍江和江蘇耕地多功能為輕度失調外,中度失調省份多位于京津冀地區(qū)、東南沿海地區(qū)和中部地區(qū),嚴重失調和極度失調省份位于西北地區(qū);2)2005-2015年,吉林耕地多功能提升至勉強協(xié)調類型,黑龍江和上海提升至瀕臨失調類型,輕度失調省份多位于四川盆地和長江中下游地區(qū);3)2018年,東北地區(qū)耕地多功能耦合協(xié)調水平較高,西北地區(qū)雖有所提升但耦合協(xié)調水平仍然較低。由此可見,耕地多功能退耦現(xiàn)象有所改善,多功能逐漸形成協(xié)同提升的良性演化關系。

    圖5 2000-2018年中國耕地生態(tài)功能空間格局

    圖6 2000-2018年中國耕地多功能耦合協(xié)調度空間格局

    耕地多功能是耕地利用的內在屬性和客觀需求,其實質是耕地系統(tǒng)與經(jīng)濟、社會和環(huán)境持續(xù)交換物質、能量與信息的過程[9]。耕地多功能是單項功能關聯(lián)反饋、互動融合而形成的有機整體,并非單項功能的獨立發(fā)揮與機械疊加。在耕地資源長期稀缺的前提下,耕地利用需求不斷衍生引致耕地功能競爭消長、損益平衡,宏觀政策制定執(zhí)行與微觀農(nóng)戶決策行動均折射為耕地多功能耦合協(xié)調水平的動態(tài)演化和區(qū)域異質分布[9]。京津冀地區(qū)、長三角地區(qū)等經(jīng)濟發(fā)達區(qū)農(nóng)業(yè)機械化、專業(yè)化程度較高,同時,其綠色的市場消費導向和先進的生態(tài)保護理念進一步促進耕地多功能有序分化[11,38]。東北地區(qū)具有優(yōu)越的資源稟賦條件和農(nóng)業(yè)生產(chǎn)基礎,加之其積極推廣采用保護性耕作技術和治理修復措施,因而該地區(qū)耕地多功能水平較高。西北地區(qū)耕地利用與改造難度較大,同時該地區(qū)技術推廣較為滯后、市場化程度提升緩慢,相對滯后的耕地利用認知程度和改造能力阻礙了耕地多功能優(yōu)化提升,仍需合理引導和積極促進。

    2.3 耕地多功能耦合協(xié)調度驅動因素

    2.3.1 單因子驅動探測

    在揭示中國耕地多功能耦合協(xié)調度演變規(guī)律的基礎上,采用地理探測器模型探析各驅動因素對耕地多功能耦合協(xié)調度的驅動力(表3)。各驅動因素均通過1%的顯著性水平檢驗,表明經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調度演變均具有驅動作用,但不同因素的驅動力不盡相同,經(jīng)營主體條件和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調度演變的驅動力較強,農(nóng)業(yè)發(fā)展水平和工業(yè)發(fā)展水平的驅動力次之。

    2.3.2 雙因子交互探測

    在識別耕地多功能耦合協(xié)調度演變驅動因素的基礎上,進一步探究各驅動因素交互作用對其影響(圖7)。具體如下:1)研究期內,各驅動因素交互作用類型均為雙因子增強或非線性增強,表明雙因素交互作用正向強化了各驅動因素對耕地多功能耦合協(xié)調度演變的驅動力。同時,雙因子增強類型由8組減至2組,非線性增強類型由20組增至26組,表明驅動因素交互作用漸趨增強、合力作用趨勢逐漸顯現(xiàn);2)農(nóng)村居民人均純收入(1)、農(nóng)村居民消費水平(2)、第一產(chǎn)業(yè)從業(yè)人員數(shù)(3)與其他驅動因素交互作用有所增強,第一產(chǎn)業(yè)增加值占比(4)、第二產(chǎn)業(yè)增加值占比(5)與其他驅動因素交互作用有所減弱,進一步驗證,經(jīng)營主體條件和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調度演變的驅動力較為主導,兩者與農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平共同作用引致耕地多功能耦合協(xié)調度呈現(xiàn)上述演變規(guī)律。因此,本文將從經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平四個方面具體闡述耕地多功能耦合協(xié)調度演變驅動機制。

    表3 耕地多功能耦合協(xié)調演變驅動因素探測結果

    圖7 耕地多功能耦合協(xié)調演變驅動因素交互作用探測結果

    2.3.3 驅動機制解析

    耕地系統(tǒng)是開放的、受人為干預明顯的復雜巨系統(tǒng),其多功能耦合協(xié)調度演變受社會經(jīng)濟諸多要素影響,即內生驅動和外生驅動綜合作用(圖8),具體驅動機制如下:

    1)經(jīng)營主體條件反映農(nóng)村居民對耕地綜合利用理念的認知程度和對耕地綜合利用開展的保障條件,是影響耕地功能發(fā)揮的基礎因素。一方面,農(nóng)村居民開始適應并追求多元化從業(yè)方式,農(nóng)村剩余勞動力加快析出并逐步獲取穩(wěn)定就業(yè)機會,以致農(nóng)村居民逐漸脫離對耕地資源的生存依賴,有助于激發(fā)耕地流轉行為,實現(xiàn)規(guī)?;?jīng)營,對耕地生產(chǎn)功能和生活功能提升具有促進作用[32,39];另一方面,農(nóng)村居民感知市場需求信息和消費升級導向的敏銳性增強,逐漸關注耕地生態(tài)功能提升和生態(tài)效益挖掘,進一步促進了耕地多功能相互融合與良性互饋。

    2)農(nóng)業(yè)發(fā)展水平反映農(nóng)業(yè)發(fā)展集約化、現(xiàn)代化水平,是影響耕地功能發(fā)揮的重要因素之一。伴隨著傳統(tǒng)農(nóng)業(yè)發(fā)展理念的逐步改觀,現(xiàn)代農(nóng)業(yè)不僅要求耕地資源發(fā)揮糧食供給作用,同時提倡耕地利用經(jīng)濟-社會-生態(tài)效益的統(tǒng)籌兼顧。農(nóng)業(yè)發(fā)展水平提高為耕地多功能利用創(chuàng)造了良好的基礎保障,促進多功能協(xié)調發(fā)展并正向反饋于農(nóng)業(yè)可持續(xù)發(fā)展。但是,農(nóng)業(yè)發(fā)展水平加速了農(nóng)業(yè)機械化推廣與水利化改造建設,雖有助于高效利用耕地,但一定程度上破壞了耕地本底生態(tài),加重了耕地利用負擔,進而影響生態(tài)功能的保持和發(fā)揮,以致耕地多功能難以有效協(xié)同。

    3)工業(yè)發(fā)展水平反映工業(yè)發(fā)展對耕地多功能顯化的支撐能力和引導作用,是影響耕地功能發(fā)揮的誘導因素。工業(yè)發(fā)展刺激消費需求增長和消費結構升級,誘發(fā)耕地利用投入結構、種植結構、經(jīng)營方式等發(fā)生改變[30],驅動傳統(tǒng)耕地利用方式轉型。另一方面,工業(yè)快速發(fā)展會攫取甚至破壞耕地資源,對耕地資源數(shù)量和質量均產(chǎn)生負面影響。為穩(wěn)定并提高有限耕地資源的經(jīng)濟產(chǎn)出,大量追加化肥、農(nóng)藥等化學生產(chǎn)要素是促進生產(chǎn)的直接手段之一。但是,高強度、高能耗的耕地利用方式抑制了耕地生態(tài)功能的發(fā)揮,從而影響耕地多功能全面協(xié)調發(fā)展。

    4)城鎮(zhèn)化水平反映城鎮(zhèn)化發(fā)展對耕地利用的改造方向和強度,是影響耕地功能發(fā)揮的關鍵因素。城鎮(zhèn)化發(fā)展逐漸破除城鄉(xiāng)發(fā)展壁壘,城鎮(zhèn)人口聚集增加和生活水平改善提升對耕地利用提出了更高品質的效益需求,推動了耕地多功能協(xié)調發(fā)展[4,30]。但同時,城市空間擴張迫使耕地向城鎮(zhèn)建設用地轉變,對耕地生產(chǎn)本底造成了負面影響。另外,城鎮(zhèn)化發(fā)展也會導致農(nóng)業(yè)勞動力質量降低、農(nóng)業(yè)發(fā)展活力喪失,雖耕地多重效益需求膨脹但多功能利用開展存在一定難度,因此城鎮(zhèn)化發(fā)展會對耕地多功能培育和發(fā)展產(chǎn)生負向影響。

    圖8 耕地多功能耦合協(xié)調演變驅動機制分析框架

    3 討 論

    耕地是自然-人工充分交互的復雜巨系統(tǒng),其多功能發(fā)揮具有明顯的階段分化和地域分異特征。因此,本文從時空二維視角綜合評價耕地多功能水平及其演變規(guī)律,同時進一步考察耕地多功能耦合協(xié)調關系,通過定量探測驅動因素和定性解構驅動機制相結合的方法深入探究耕地多功能耦合協(xié)調演變的作用機理,嘗試遵循“表征刻畫-動因識別-機理剖析”的研究邏輯豐富耕地多功能理性認識。耕地多功能為耕地利用提供了合理方向并為耕地價值顯化提供了具體路徑,耕地多功能管理應結合資源優(yōu)勢與發(fā)展環(huán)境以制定精細化、差異化的管控方案及保障措施。

    在“高效利用+高度保護”的可持續(xù)理念引導下,耕地多功能管理不僅應該針對性滿足區(qū)域社會經(jīng)濟增長和居民生活改善對農(nóng)業(yè)和耕地的適度要求,充分利用社會多元需求和現(xiàn)代化發(fā)展提供的有利機遇顯化自身價值,同時應強調耕地質量提升和生態(tài)系統(tǒng)穩(wěn)定以實現(xiàn)耕地永續(xù)利用。在耕地多功能利用發(fā)展較為活躍的地區(qū),政府和各經(jīng)營主體應充分利用地區(qū)資源優(yōu)勢和市場優(yōu)勢,通過釋放市場消費潛力、激發(fā)耕地多功能利用活力。在耕地資源有限但市場需求旺盛的背景下,該類地區(qū)應充分利用發(fā)展耕地多功能的先行條件,著力打造生態(tài)綠色農(nóng)業(yè)、現(xiàn)代設施農(nóng)業(yè)、典型特色農(nóng)業(yè)等,通過推動多樣化、高端化、綠色化的農(nóng)業(yè)產(chǎn)業(yè)融合,保障經(jīng)營主體收益、助力耕地多功能價值顯化。同時,這類地區(qū)應自下而上地推進耕地多功能管理制度完善和路徑創(chuàng)新,為耕地可持續(xù)利用提供良好的制度支撐。在耕地多功能耦合協(xié)調水平較低的地區(qū),決策者可以通過加強規(guī)劃布局引導、指標控制、經(jīng)濟補貼、監(jiān)測預警等多種方式實行耕地多功能管理,在不降低耕地商品性生產(chǎn)功能的前提下,激勵農(nóng)戶參與耕地衍生功能的培育和挖掘,發(fā)揮農(nóng)戶保護與利用耕地的主觀能動性,因勢利導、因地制宜地實現(xiàn)耕地利用轉型。

    4 結 論

    本文以2000-2018年中國30個省級行政區(qū)為研究單元,開展“生產(chǎn)-生活-生態(tài)”三維耕地多功能評價并采用核密度估計和可視化制圖刻畫耕地多功能時空演進特征,進而采用耦合協(xié)調度模型分析耕地多功能耦合協(xié)調關系及其演變規(guī)律,最終采用地理探測器識別耕地多功能耦合協(xié)調演變的驅動因素。研究主要結論如下:

    1)從時序演變特征來看,耕地生產(chǎn)功能水平呈現(xiàn)上升-下降-上升的波動變化,生活功能水平先上升后下降,生態(tài)功能水平則整體小幅提升,各功能均具有極化特征。從空間格局特征來看,耕地生產(chǎn)功能、生活功能和生態(tài)功能分別呈現(xiàn)“東高西低”“北高南低”和“西高東低”的空間非均衡格局。

    2)從時序演變來看,中國耕地多功能耦合協(xié)調水平整體穩(wěn)步提升但水平仍然較低,多數(shù)省份耕地多功能為失調狀態(tài),表明中國耕地多功能正在由無序沖突轉為有序發(fā)展,但耕地多功能耦合協(xié)調水平仍需進一步提升。從空間格局來看,中國耕地多功能耦合協(xié)調度逐漸呈現(xiàn)集聚提升的演化特征。

    3)耕地多功能耦合協(xié)調度演變受經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平影響顯著,但不同因素的驅動力不盡相同。各驅動因素交互作用類型均為雙因子增強或非線性增強,且非線性增強類型逐漸增多,表明耕地多功能耦合協(xié)調演變?yōu)閮壬寗雍屯庠隍寗与p重影響的結果。

    本文針對長時序、省級尺度耕地多功能開展實證研究,受限于統(tǒng)計數(shù)據(jù)和調研資料的可獲取性,并未于更細尺度、更長時序揭示耕地多功能時空規(guī)律及演化機理,也未能將耕地多功能評價結果與農(nóng)戶走訪調查情況進行對比修正。因此,未來研究將采用多源數(shù)據(jù)、多種技術手段全面考察并深入挖掘耕地多功能演化表征和驅動機制,同時借助專家訪談、實地調研等多種研究方法輔助修正研究結果,以期充分揭示中國耕地多功能結構性、趨勢性、系統(tǒng)性的發(fā)展規(guī)律和內在機理,為實現(xiàn)耕地資源優(yōu)化配置和功能協(xié)同提升提供決策支撐。

    [1] 鄒利林,李裕瑞,劉彥隨,等. 基于要素視角的耕地“三生”功能理論建構與實證研究[J]. 地理研究,2021,40(3):839-855.

    ZOU Lilin, LI Yurui, LIU Yansui, et al. Theory building and empirical research of production-living-ecological function of cultivated land based on the elements[J]. Geographical Research, 2021, 40(3): 839-855. (in Chinese with English abstract)

    [2] 熊昌盛,張永蕾,王雅娟,等. 中國耕地多功能評價及分區(qū)管控[J]. 中國土地科學,2021,35(10):104-114.

    XIONG Changsheng, ZHANG Yonglei, WANG Yajuan, et al. Multi-function evaluation and zoning control of cultivated land in China[J]. Chinese Land Science, 2021, 35(10): 104-114. (in Chinese with English abstract)

    [3] 朱從謀,李武艷,杜瑩瑩,等. 浙江省耕地多功能價值時空變化與權衡-協(xié)同關系[J]. 農(nóng)業(yè)工程學報,2020,36(14):263-272.

    ZHU Congmou, LI Wuyan, DU Yingying, et al. Spatial-temporal change, trade-off and synergy relationships of cropland multifunctional value in Zhejiang Province, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 263-272. (in Chinese with English abstract)

    [4] 張英男,龍花樓,戈大專,等. 黃淮海平原耕地功能演變的時空特征及其驅動機制[J]. 地理學報,2018,73(3):518-534.

    ZHANG Yingnan, LONG Hualou, GE Dazhuan, et al. Spatio-temporal characteristics and dynamic mechanism of farmland functions evolution in the Huang-Huai-Hai Plain[J]. Acta Geographica Sinica, 2018, 73(3): 518-534. (in Chinese with English abstract)

    [5] 安悅,譚雪蘭,李印齊,等. 洞庭湖地區(qū)耕地功能時空演變特征及影響因素研究[J]. 地理科學,2022,42(7):1272-1282.

    AN Yue, TAN Xuelan, LI Yinqi, et al.Spatio-temporal evolution characteristics and influencing factors of cultivated land functions in the Dongting Lake Aera[J]. Scientia Geographica Sinica, 2022, 42(7): 1272-1282. (in Chinese with English abstract)

    [6] 宋小青,歐陽竹. 耕地多功能內涵及其對耕地保護的啟示[J]. 地理科學進展,2012,31(7):859-868.

    SONG Xiaoqing, OUYANG Zhu. Connotation of multifunctional cultivated land and its implications for cultivated land protection[J]. Progress in Geography, 2012, 31(7): 859-868. (in Chinese with English abstract)

    [7] HELMING K, PEREZ-SOBA M, TABBUSH P. Sustainability Impact Assessment of Land Use Changes[M]. Berlin: Springer, 2008: 375-404.

    [8] MANDER U, HELMING K, WIGGERING H. Multifunctional Land Use: Meeting Future Demands for Landscape Goods and Services[M]. Berlin: Springer, 2007: 1-13.

    [9] 宋小青,歐陽竹. 中國耕地多功能管理的實踐路徑探討[J]. 自然資源學報,2012,27(4):540-551.

    SONG Xiaoqing, OUYANG Zhu. Route of multifunctional cultivated land management in China[J]. Journal of Natural Resources, 2012, 27(4): 540-551. (in Chinese with English abstract)

    [10] 宋小青,李心怡. 區(qū)域耕地利用功能轉型的理論解釋與實證[J]. 地理學報,2019,74(5):992-1010.

    SONG Xiaoqing, LI Xinyi. Theoretical explanation and case study of regional cultivated land use function transition[J]. Acta Geographica Sinica, 2019, 74(5): 992-1010. (in Chinese with English abstract)

    [11] 范業(yè)婷,金曉斌,項曉敏,等. 蘇南地區(qū)耕地多功能評價與空間特征分析[J]. 資源科學,2018,40(5):980-992.

    FAN Yeting, JIN Xiaobin, XIANG Xiaomin, et al. Evaluation and spatial characteristics of arable land multifunction in southern Jiangsu[J]. Resources Science, 2018, 40(5): 980-992. (in Chinese with English abstract)

    [12] WANG X J, WANG D Y, WU S Z, et al. Cultivated land multifunctionality in undeveloped peri-urban agriculture areas in China: Implications for sustainable land management[J]. Journal of Environmental Management, 2023, 325: 116500.

    [13] ZHANG S Y, HU W Y, LI M R, et al. Multiscale research on spatial supply-demand mismatches and synergic strategies of multifunctional cultivated land[J]. Journal of Environmental Management, 2021, 299: 113605.

    [14] JIANG G H, WANG M Z, QU Y B, et al. Towards cultivated land multifunction assessment in China: Applying the “influencing factors-functions-products-demands” integrated framework[J]. Land Use Policy, 2020, 99: 104982.

    [15] 馬才學,金瑩,柯新利,等. 基于全排列多邊形圖示法的湖北省耕地多功能強度與協(xié)調度典型模式探究[J]. 中國土地科學,2018,32(4):51-58.

    MA Caixue, JIN Ying, KE Xinli, et al. Study on the typical modes of farmland multi-functional intensity and coordination in Hubei province based on the entire-array-polygon indicator method[J]. Chinese Land Science, 2018, 32(4): 51-58. (in Chinese with English abstract)

    [16] 向敬偉,廖曉莉,宋小青,等. 中國耕地多功能的區(qū)域收斂性[J]. 資源科學,2019,41(11):1959-1971.

    XIANG Jingwei, LIAO Xiaoli, SONG Xiaoqing, et al. Regional convergence of cultivated land multifunctions in China[J]. Resources Science, 2019, 41(11): 1959-1971. (in Chinese with English abstract)

    [17] 方瑩,王靜,孔雪松,等. 耕地利用多功能權衡關系測度與分區(qū)優(yōu)化——以河南省為例[J]. 中國土地科學,2018,32(11):57-64.

    FANG Ying, WANG Jing, KONG Xuesong, et al. Trade-off relation measurement and zoning optimization of multi-functionality of cultivated land use: A case study of Henan province[J]. Chinese Land Science, 2018, 32(11): 57-64. (in Chinese with English abstract)

    [18] 宋小青,吳志峰,歐陽竹. 1949年以來中國耕地功能變化[J]. 地理學報,2014,69(4):435-447.

    SONG Xiaoqing, WU Zhifeng, OUYANG Zhu. Changes of cultivated land function in China since 1949[J]. Acta Geographica Sinica, 2014, 69(4): 435-447. (in Chinese with English abstract)

    [19] 胡偉艷,朱慶瑩,張安錄,等. 總量與結構視角耕地多功能對農(nóng)業(yè)經(jīng)濟增長的影響:以湖北省為例[J]. 中國土地科學,2018,32(5):62-70.

    HU Weiyan, ZHU Qingying, ZHANG Anlu, et al. Impacts of multifunctional farmland on agricultural economic growth at county-level from the perspective of quantity and structure: Evidence from Hubei province[J]. Chinese Land Science, 2018, 32(5): 62-70. (in Chinese with English abstract)

    [20] 許多藝,濮勵杰,黃思華,等. 江蘇省耕地多功能時空動態(tài)分析及對耕地數(shù)量變化響應研究[J]. 長江流域資源與環(huán)境,2022,31(3):575-587.

    XU Duoyi, PU Lijie, HUANG Sihua, et al. Spatial and temporal dynamic analysis of cultivated land multifunction in Jiangsu Province and its response to cultivated land change[J]. Resources and Environment in the Yangtze Basin, 2022, 31(3): 575-587. (in Chinese with English abstract)

    [21] 高星,宋昭穎,李晨曦,等. 城鄉(xiāng)梯度下的耕地多功能價值空間分異特征[J]. 農(nóng)業(yè)工程學報,2021,37(16):251-259.

    GAO Xing, SONG Zhaoying, LI Chenxi, et al. Spatial differentiation characteristics of cultivated land multifunctional value under urban-rural gradient[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(16): 251-259. (in Chinese with English abstract)

    [22] 盧新海,唐一峰,易家林,等. 基于空間計量模型的耕地利用轉型對農(nóng)業(yè)經(jīng)濟增長影響研究[J]. 中國土地科學,2019,33(6):53-61.

    LU Xinhai, TANG Yifeng, YI Jialin, et al. Study on the impact of cultivated land use transition on agricultural economic growth based on spatial econometric model[J]. Chinese Land Science, 2019, 33(6): 53-61. (in Chinese with English abstract)

    [23] 盧新海,楊喜,陳澤秀. 中國城市土地綠色利用效率測度及其時空演變特征[J]. 中國人口·資源與環(huán)境,2020,30(8):83-91.

    LU Xinhai, YANG Xi, CHEN Zexiu. Measurement and temporal-spatial evolution characteristics of urban land green use efficiency in China[J]. China Population, Resources and Environment, 2020, 30(8): 83-91. (in Chinese with English abstract)

    [24] LU X H, KUANG B, LI J. Regional difference decomposition and policy implications of China’s urban land use efficiency under the environmental restriction[J]. Habitat International, 2018, 77: 32-39.

    [25] 徐維祥,李露,周建平,等. 鄉(xiāng)村振興與新型城鎮(zhèn)化耦合協(xié)調的動態(tài)演進及其驅動機制[J]. 自然資源學報,2020,35(9):2044-2062.

    XU Weixiang, LI Lu, ZHOU Jianping, et al. The dynamic evolution and its driving mechanism of coordination of rural rejuvenation and new urbanization[J]. Journal of Natural Resources, 2020, 35(9): 2044-2062. (in Chinese with English abstract)

    [26] CAI J, LI X P, LIU L J, et al. Coupling and coordinated development of new urbanization and agro-ecological environment in China[J]. Science of the Total Environment, 2021, 776: 145837.

    [27] LI W J, WANG Y, XIE S Y, et al. Coupling coordination analysis and spatiotemporal heterogeneity between urbanization and ecosystem health in Chongqing municipality, China[J]. Science of the Total Environment, 2021, 791: 148311.

    [28] 朱慶瑩,陳銀蓉,胡偉艷,等. 中國土地集約利用與區(qū)域生態(tài)效率耦合協(xié)調度時空格局[J]. 農(nóng)業(yè)工程學報,2020,36(4):234-243.

    ZHU Qingying, CHEN Yingrong, HU Weiyan, et al. Spatiotemporal pattern of coupling coordination degree between land intensive use and regional ecological efficiency in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 234-243. (in Chinese with English abstract)

    [29] 蓋美,張福祥. 遼寧省區(qū)域碳排放-經(jīng)濟發(fā)展-環(huán)境保護耦合協(xié)調分析[J]. 地理科學,2018,38(5):764-772.

    GAI Mei, ZHANG Fuxiang. Regional carbon emissions, economic development and environmental protection coupling in Liaoning province[J]. Scientia Geographica Sinica, 2018, 38(5): 764-772. (in Chinese with English abstract)

    [30] 張文斌,張志斌,董建紅,等. 多尺度視角下耕地利用功能轉型及驅動力分析:以甘肅省為例[J]. 地理科學,2021,41(5):900-910.

    ZHANG Wenbin, ZHANG Zhibin, DONG Jianhong, et al. Transformation and driving forces of cultivated land utilization function from a multi-scale perspective in Gansu Province[J]. Scientia Geographica Sinica, 2021, 41(5): 900-910. (in Chinese with English abstract)

    [31] 衛(wèi)新東,林良國,羅平平,等. 耕地多功能耦合協(xié)調發(fā)展時空格局與驅動力分析[J]. 農(nóng)業(yè)工程學報,2022,38(4):260-269.

    WEI Xindong, LIN Liangguo, LUO Pingping, et al. Spatiotemporal pattern and driving force analysis of multi-functional coupling coordinated development of cultivated land[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(4): 260-269. (in Chinese with English abstract)

    [32] 李怡,方斌,李裕瑞,等. 城鎮(zhèn)化進程中耕地多功能權衡/協(xié)同關系演變及其驅動機制[J]. 農(nóng)業(yè)工程學報,2022,38(8):244-254.

    LI Yi, FANG Bin, LI Yurui, et al. Trade-off and synergy evolution of farmland functions and its dynamic mechanism in the process of urbanization[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(8): 244-254.

    [33] 張玥,代亞強,陳媛媛,等. 中國耕地低碳利用效率時空演變及其驅動因素[J]. 農(nóng)業(yè)工程學報,2022,38(8):234-243.

    ZHANG Yue, DAI Yaqiang, CHEN Yuanyuan, et al. Spatial-temporal evolution and driving factors of low-carbon use efficiency of cultivated land in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(8): 234-243. (in Chinese with English abstract)

    [34] 王勁峰,徐成東. 地理探測器:原理與展望[J]. 地理學報,2017,72(1):116-134.

    WANG Jinfeng, XU Chengdong. Geodetector: Principle and prospective[J]. Acta Geographica Sinica, 2017, 72(1): 116-134. (in Chinese with English abstract)

    [35] 匡兵,范翔宇,盧新海. 中國耕地利用綠色轉型效率的時空分異特征及其影響因素[J]. 農(nóng)業(yè)工程學報,2021,37(21):269-277.

    KUANG Bing, FAN Xiangyu, LU Xinhai. Spatial-temporal differentiation characteristics of the efficiency of green transformation of cultivated land use and its affecting factors in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(21): 269-277. (in Chinese with English abstract)

    [36] 陳正發(fā),史東梅,何偉,等. 1980-2015年云南坡耕地資源時空分布及演變特征分析[J]. 農(nóng)業(yè)工程學報,2019,35(15):256-265.

    CHEN Zhengfa, SHI Dongmei, HE Wei, et al. Spatio-temporal distribution and evolution characteristics of slope farmland resources in Yunnan from 1980 to 2015[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(15): 256-265. (in Chinese with English abstract)

    [37] 李升發(fā),李秀彬. 中國山區(qū)耕地利用邊際化表現(xiàn)及其機理[J]. 地理學報,2018,73(5):803-817.

    LI Shengfa, LI Xiubin. Economic characteristics and the mechanism of farmland marginalization in mountainous areas of China[J]. Acta Geographica Sinica, 2018, 73(5): 803-817. (in Chinese with English abstract)

    [38] 付慧,劉艷軍,孫宏日,等. 京津冀地區(qū)耕地利用轉型時空分異及驅動機制[J]. 地理科學進展,2020,39(12):1985-1998.

    FU Hui, LIU Yanjun, SUN Hongri, et al. Spatiotemporal characteristics and dynamic mechanism of cultivated land use transition in the Beijing-Tianjin-Hebei region[J]. Progress in Geography, 2020, 39(12): 1985-1998. (in Chinese with English abstract)

    [39] 秦彥杰,劉欣,趙艷霞,等. 1949年以來河北省耕地利用轉型特征及過程[J]. 農(nóng)業(yè)工程學報,2022,38(12):260-269.

    QIN Yanjie, LIU Xin, ZHAO Yanxia, et al. Characteristics and process of the cultivated land use transformation in Hebei Province since 1949[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(12): 260-269. (in Chinese with English abstract)

    Spatial-temporal evolution and driving factors of cultivated land multifunctional coupling coordination development in China

    ZHANG Yue, DAI Yaqiang, CHEN Yuanyuan, KE Xinli※

    (,,430070,)

    To explore the spatial-temporal evolution and driving factors of cultivated land multifunctional coupling coordination development, this paper applying the data of 30 provinces in China from 2000 to 2018 evaluated cultivated land production-living-ecological functions and revealed the spatial-temporal evolution of cultivated land multifunction with kernel density estimation and visualization mapping. Furthermore, this paper adopted coupling coordination degree model to analyze the cultivated land multifunctional coupling coordination development. Finally, this paper explored driving factors and driving mechanism of cultivated land multifunctional coupling coordination development by geographical detector. The results showed that: 1) from the perspective of temporal evolution, cultivated land production function level showed a fluctuation evolution trend of “rising-falling-rising”, the living function level showed an evolution trend of “rising-falling”, and the ecological function level increased slightly. All cultivated land functions had polarized characteristics. From the perspective of spatial pattern, the cultivated land production function level changed within [0.030, 0.608], showing a non-equilibrium spatial pattern of “high in the east, low in the west”. The living function level changed within [0.042, 0.672], showing a non-equilibrium spatial pattern of “high in the north, low in the south”. The ecological function level changed within [0.058, 0.897], showing a non-equilibrium spatial pattern of “high in the west, low in the east”. The areas with high level of production and living function were concentrated in the Northeast Plain, the Huang-Huai-Hai Plain and the Yangtze Region. The areas with low level of production function were located in the northwest and southwest regions, and the areas with low level of living function were located in the Beijing-Tianjin, Yunnan-Guizhou and southern regions in China. The area with high level of ecological function areas were located in the northeast and southwest regions, while the low-level areas were located in the Beijing-Tianjin-Hebei region, the Yangtze River Delta region and the central region of the Yangtze River. 2) The cultivated land multifunctional coupling coordination level changed within [0.093, 0.554], and the overall level had steadily improved but was still low. The cultivated land multifunctional coupling coordination levels of most provinces were still out of balance. The spatial pattern of cultivated land multifunctional coupling coordination level showed an evolution trend of agglomeration and promotion. In detail, the cultivated land multifunctional coupling coordination level of the northeast China, Beijing-Tianjin-Hebei region and the Yangtze River Delta were relatively high, while that was low in the northwestern part of China. 3) The evolution of cultivated land multifunctional coupling coordination level was significantly affected by the business entities conditions, agricultural development level, industrial development level and urbanization level. In short, the driving force of the business entities conditions and urbanization level were more significant than that of the agricultural development level and industrial development level. Moreover, the interaction type of two driving factors was two-factor enhancement or nonlinear enhancement, which indicated that the driving force of various factors had obvious divergence and interaction of factors positively strengthened the single-factor driving force, indicated that the divergence and interaction of various driving factors positively strengthened the single-factor driving force. In conclusion, the research results provide the guiding effect for deepening the study on cultivated land multifunction and decision support for promoting the differential governance of cultivated land resources and mutual promotion of cultivated land functions.

    land use; cultivated land; multifunction; coupling coordination degree; spatial-temporal evolution; driving factors

    2022-09-23

    2023-02-10

    國家自然科學基金面上項目(41971240);國家社科基金后期資助項目(19FGLB071);中央高?;究蒲袠I(yè)務費專項資金資助項目(2662022XCZX012)

    張玥,博士生,研究方向為土地資源管理。Email:yuez@webmail.hzau.edu.cn

    柯新利,博士,教授,博士生導師,研究方向為土地利用優(yōu)化。Email:kexl@mail.hzau.edu.cn

    10.11975/j.issn.1002-6819.202209185

    S28; F301.21

    A

    1002-6819(2023)-07-0244-12

    張玥,代亞強,陳媛媛,等. 中國耕地多功能耦合協(xié)調時空演變及其驅動因素[J]. 農(nóng)業(yè)工程學報,2023,39(7):244-255. doi:10.11975/j.issn.1002-6819.202209185 http://www.tcsae.org

    ZHANG Yue, DAI Yaqiang, CHEN Yuanyuan, et al. Spatial-temporal evolution and driving factors of cultivated land multifunctional coupling coordination development in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2023, 39(7): 244-255. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.202209185 http://www.tcsae.org

    猜你喜歡
    耕地耦合驅動
    自然資源部:加強黑土耕地保護
    我國將加快制定耕地保護法
    非Lipschitz條件下超前帶跳倒向耦合隨機微分方程的Wong-Zakai逼近
    基于模糊PI控制的驅動防滑仿真系統(tǒng)分析
    保護耕地
    北京測繪(2021年12期)2022-01-22 03:33:36
    新增200億元列入耕地地力保護補貼支出
    屈宏斌:未來五年,雙輪驅動,砥礪前行
    軌旁ATC系統(tǒng)門控柜接收/驅動板改造
    基于S3C6410的Wi-Fi驅動移植實現(xiàn)
    基于“殼-固”耦合方法模擬焊接裝配
    大型鑄鍛件(2015年5期)2015-12-16 11:43:20
    日本成人三级电影网站| 国产成人午夜福利电影在线观看| 级片在线观看| 久久久久久久久久久丰满| 中文字幕制服av| 午夜精品国产一区二区电影 | 少妇猛男粗大的猛烈进出视频 | 亚洲久久久久久中文字幕| 亚洲不卡免费看| 乱码一卡2卡4卡精品| 国产在线精品亚洲第一网站| 国产探花在线观看一区二区| 中文资源天堂在线| 九九久久精品国产亚洲av麻豆| 免费观看精品视频网站| 中文字幕免费在线视频6| 亚洲成人久久爱视频| 午夜激情欧美在线| 久久精品国产99精品国产亚洲性色| 亚洲人与动物交配视频| 能在线免费看毛片的网站| 成年女人永久免费观看视频| 欧美一区二区国产精品久久精品| 校园人妻丝袜中文字幕| 一区二区三区免费毛片| 亚洲国产色片| 高清在线视频一区二区三区 | 国产人妻一区二区三区在| 一级黄片播放器| 国产精品久久电影中文字幕| 国产高清三级在线| 国产综合懂色| 亚洲不卡免费看| 国产在线男女| 大香蕉久久网| 男的添女的下面高潮视频| eeuss影院久久| 高清日韩中文字幕在线| 男女下面进入的视频免费午夜| 男人舔女人下体高潮全视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 老女人水多毛片| 非洲黑人性xxxx精品又粗又长| 久久久久久国产a免费观看| 久久久久久久久久黄片| 国产成人91sexporn| 婷婷六月久久综合丁香| 黄片wwwwww| 国内揄拍国产精品人妻在线| 亚洲欧美日韩无卡精品| 亚洲精品乱码久久久久久按摩| 男人舔奶头视频| 国产极品天堂在线| av免费观看日本| 国产 一区精品| 日本黄色视频三级网站网址| 一个人看的www免费观看视频| 人妻制服诱惑在线中文字幕| 91狼人影院| 天堂av国产一区二区熟女人妻| 精品久久久久久久久久免费视频| 亚洲天堂国产精品一区在线| 婷婷色av中文字幕| 99在线视频只有这里精品首页| 91久久精品国产一区二区成人| 久久久精品欧美日韩精品| 深夜精品福利| 26uuu在线亚洲综合色| 国产一区二区三区av在线 | 九九热线精品视视频播放| 床上黄色一级片| 九色成人免费人妻av| 免费大片18禁| 最近手机中文字幕大全| 韩国av在线不卡| 国产精品久久久久久精品电影小说 | 能在线免费看毛片的网站| 国产片特级美女逼逼视频| 国产一级毛片在线| 免费观看精品视频网站| 国产黄片美女视频| 91午夜精品亚洲一区二区三区| 国产不卡一卡二| 色综合站精品国产| а√天堂www在线а√下载| 亚洲四区av| 嘟嘟电影网在线观看| 人妻系列 视频| 成熟少妇高潮喷水视频| 十八禁国产超污无遮挡网站| 午夜免费男女啪啪视频观看| 日韩国内少妇激情av| 伦理电影大哥的女人| 精品熟女少妇av免费看| 日韩精品青青久久久久久| 精品久久久久久久久亚洲| 热99在线观看视频| 精品久久国产蜜桃| 久久精品国产99精品国产亚洲性色| 人妻制服诱惑在线中文字幕| 亚洲18禁久久av| 亚洲国产色片| 在线国产一区二区在线| 日韩精品有码人妻一区| 午夜福利在线观看吧| 亚洲人成网站高清观看| 99国产极品粉嫩在线观看| 偷拍熟女少妇极品色| 性欧美人与动物交配| 日日摸夜夜添夜夜添av毛片| 亚洲国产高清在线一区二区三| 黄片无遮挡物在线观看| 男人舔奶头视频| 美女国产视频在线观看| 我要看日韩黄色一级片| 一级黄片播放器| 伊人久久精品亚洲午夜| 免费电影在线观看免费观看| 国产精品一二三区在线看| 日韩成人伦理影院| 丰满的人妻完整版| 午夜激情福利司机影院| 伊人久久精品亚洲午夜| 欧美精品国产亚洲| 亚洲一区高清亚洲精品| 国产老妇女一区| 国产一区亚洲一区在线观看| 成人二区视频| 午夜福利在线观看免费完整高清在 | 性欧美人与动物交配| 久久久久久久久久久免费av| 亚洲精品粉嫩美女一区| 91精品国产九色| 成年版毛片免费区| 精品久久久久久久久av| 中文字幕制服av| 三级男女做爰猛烈吃奶摸视频| 国产一区亚洲一区在线观看| 如何舔出高潮| 人体艺术视频欧美日本| 亚洲国产欧美人成| 欧美成人免费av一区二区三区| 国产白丝娇喘喷水9色精品| 国产精品,欧美在线| 人妻少妇偷人精品九色| 老熟妇乱子伦视频在线观看| 国产黄a三级三级三级人| 99热6这里只有精品| 搞女人的毛片| 亚洲av成人av| 亚洲欧美精品综合久久99| av福利片在线观看| 我要看日韩黄色一级片| 亚洲在线自拍视频| 亚洲乱码一区二区免费版| 69av精品久久久久久| 观看免费一级毛片| 久久欧美精品欧美久久欧美| а√天堂www在线а√下载| 国产又黄又爽又无遮挡在线| 午夜福利在线在线| 一本久久精品| 国产伦一二天堂av在线观看| 久久中文看片网| 可以在线观看的亚洲视频| 91狼人影院| www.色视频.com| 成人av在线播放网站| 欧美+日韩+精品| 男女做爰动态图高潮gif福利片| 麻豆国产av国片精品| 中文资源天堂在线| 亚洲精品色激情综合| 美女 人体艺术 gogo| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 欧美成人免费av一区二区三区| 亚洲第一区二区三区不卡| 噜噜噜噜噜久久久久久91| 国产91av在线免费观看| 一区二区三区高清视频在线| 91久久精品国产一区二区成人| 国产日韩欧美在线精品| 亚洲成人中文字幕在线播放| 国产精品久久久久久精品电影小说 | 婷婷亚洲欧美| 国产精品爽爽va在线观看网站| 国产91av在线免费观看| 久久欧美精品欧美久久欧美| 久久九九热精品免费| 国产伦在线观看视频一区| 久久这里只有精品中国| 亚洲一区高清亚洲精品| 全区人妻精品视频| 校园人妻丝袜中文字幕| 可以在线观看毛片的网站| 亚洲国产精品sss在线观看| 插阴视频在线观看视频| 亚洲在久久综合| 我的女老师完整版在线观看| 成熟少妇高潮喷水视频| 日本五十路高清| 天天躁日日操中文字幕| 午夜免费男女啪啪视频观看| 国产伦精品一区二区三区视频9| 亚洲国产色片| 熟女电影av网| 卡戴珊不雅视频在线播放| 日韩视频在线欧美| 五月伊人婷婷丁香| 中文字幕精品亚洲无线码一区| 亚洲色图av天堂| 国产成年人精品一区二区| 日本撒尿小便嘘嘘汇集6| 成人亚洲欧美一区二区av| 成人特级av手机在线观看| 99久久人妻综合| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 国产探花在线观看一区二区| 国产精品一及| 午夜福利成人在线免费观看| 精品人妻偷拍中文字幕| 一本久久中文字幕| 久久久久久久久久久免费av| 黄色视频,在线免费观看| 级片在线观看| avwww免费| 乱系列少妇在线播放| 老熟妇乱子伦视频在线观看| 免费看美女性在线毛片视频| 日韩一区二区视频免费看| 久久精品夜夜夜夜夜久久蜜豆| 一本一本综合久久| 国产老妇女一区| 欧美日韩乱码在线| 欧美高清成人免费视频www| 精品一区二区免费观看| 如何舔出高潮| 亚洲乱码一区二区免费版| 男人和女人高潮做爰伦理| 熟妇人妻久久中文字幕3abv| 久久久久久久午夜电影| 亚洲欧美日韩高清在线视频| 色吧在线观看| 成人亚洲精品av一区二区| 成人毛片60女人毛片免费| 日日干狠狠操夜夜爽| 免费人成在线观看视频色| 久久精品久久久久久噜噜老黄 | 小蜜桃在线观看免费完整版高清| 亚洲自偷自拍三级| 色尼玛亚洲综合影院| 日本与韩国留学比较| 午夜精品在线福利| 免费av不卡在线播放| 嫩草影院精品99| 午夜精品一区二区三区免费看| 精品久久久久久久久亚洲| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| .国产精品久久| 国产熟女欧美一区二区| 一级黄色大片毛片| 久久99蜜桃精品久久| 91在线精品国自产拍蜜月| 久久人妻av系列| 在线观看66精品国产| avwww免费| 国产精品嫩草影院av在线观看| 亚洲在线自拍视频| 欧美日韩一区二区视频在线观看视频在线 | 精品久久久久久久久久免费视频| 亚洲中文字幕一区二区三区有码在线看| av免费在线看不卡| 日韩欧美精品免费久久| 联通29元200g的流量卡| 日日啪夜夜撸| 日韩成人av中文字幕在线观看| 亚洲久久久久久中文字幕| 波野结衣二区三区在线| 日韩欧美一区二区三区在线观看| 午夜福利视频1000在线观看| 日本一本二区三区精品| 国产黄片美女视频| 亚洲精品日韩av片在线观看| 夜夜爽天天搞| 日本三级黄在线观看| 免费观看在线日韩| 亚洲综合色惰| 国产精品精品国产色婷婷| 亚洲av免费高清在线观看| 精品久久久久久久久久久久久| 中文亚洲av片在线观看爽| 变态另类丝袜制服| 成人二区视频| 女同久久另类99精品国产91| 国产成人精品婷婷| 秋霞在线观看毛片| 黄色日韩在线| 色综合站精品国产| 日韩欧美国产在线观看| 亚洲自拍偷在线| 久久久色成人| 亚洲av二区三区四区| 国国产精品蜜臀av免费| 一级黄片播放器| 午夜免费男女啪啪视频观看| 久久99精品国语久久久| 久久6这里有精品| 国产乱人偷精品视频| 婷婷亚洲欧美| 国产精品久久久久久精品电影| 日本欧美国产在线视频| 久久久久久久亚洲中文字幕| 国产一级毛片七仙女欲春2| 麻豆国产97在线/欧美| 国产片特级美女逼逼视频| 国产午夜精品一二区理论片| 最好的美女福利视频网| 偷拍熟女少妇极品色| 嘟嘟电影网在线观看| 色综合色国产| 亚洲一级一片aⅴ在线观看| 亚洲av免费高清在线观看| 国产伦精品一区二区三区四那| 男人舔奶头视频| 青春草亚洲视频在线观看| 欧美成人精品欧美一级黄| 夜夜看夜夜爽夜夜摸| 美女xxoo啪啪120秒动态图| 精品久久久久久久人妻蜜臀av| 哪里可以看免费的av片| 国产精品电影一区二区三区| 国产黄a三级三级三级人| 亚洲在线自拍视频| 国产国拍精品亚洲av在线观看| av在线老鸭窝| 国产精品免费一区二区三区在线| 高清毛片免费看| av免费观看日本| 男人狂女人下面高潮的视频| 国产精品无大码| 色视频www国产| 日本五十路高清| 国产中年淑女户外野战色| 成人毛片60女人毛片免费| 久久午夜亚洲精品久久| 国内精品美女久久久久久| 日本欧美国产在线视频| 成人毛片a级毛片在线播放| 国产片特级美女逼逼视频| 国产成人精品久久久久久| 午夜福利视频1000在线观看| 日本黄大片高清| 国内精品宾馆在线| 男的添女的下面高潮视频| 亚洲国产精品成人久久小说 | 国产老妇伦熟女老妇高清| 天堂中文最新版在线下载 | 精品少妇黑人巨大在线播放 | 国产精品伦人一区二区| 亚洲国产精品sss在线观看| 成人鲁丝片一二三区免费| 在线观看av片永久免费下载| 人妻系列 视频| 99国产极品粉嫩在线观看| 久久精品国产亚洲网站| 亚洲精品日韩av片在线观看| 亚洲久久久久久中文字幕| 亚洲欧美成人综合另类久久久 | 婷婷精品国产亚洲av| 久久6这里有精品| 色综合色国产| 亚洲高清免费不卡视频| 国产高清不卡午夜福利| 男插女下体视频免费在线播放| 小说图片视频综合网站| 午夜精品在线福利| 日韩av不卡免费在线播放| 三级毛片av免费| 亚洲激情五月婷婷啪啪| 国产日本99.免费观看| 欧美日本亚洲视频在线播放| 又爽又黄无遮挡网站| 97人妻精品一区二区三区麻豆| 又爽又黄无遮挡网站| 亚洲成人中文字幕在线播放| 久久精品国产亚洲av天美| 国产伦精品一区二区三区视频9| 日本三级黄在线观看| 美女国产视频在线观看| 国内精品久久久久精免费| 亚洲熟妇中文字幕五十中出| 少妇的逼水好多| 国产一区二区三区在线臀色熟女| 亚洲精华国产精华液的使用体验 | 国产黄片美女视频| 午夜福利视频1000在线观看| 午夜福利高清视频| a级毛色黄片| 少妇人妻一区二区三区视频| 国产久久久一区二区三区| 99热网站在线观看| 国内精品久久久久精免费| 久久久久久国产a免费观看| 免费av观看视频| 少妇裸体淫交视频免费看高清| 成人高潮视频无遮挡免费网站| 卡戴珊不雅视频在线播放| 在线a可以看的网站| 日韩高清综合在线| 我的老师免费观看完整版| 久久久国产成人免费| 欧美+亚洲+日韩+国产| 日本黄色片子视频| 亚洲欧美成人综合另类久久久 | 亚洲第一区二区三区不卡| 一本久久精品| 国产伦理片在线播放av一区 | 白带黄色成豆腐渣| 久久精品国产清高在天天线| 精品久久久久久久久久久久久| 久久这里只有精品中国| 又黄又爽又刺激的免费视频.| 国产高清视频在线观看网站| 一级毛片aaaaaa免费看小| 国产精品日韩av在线免费观看| 人人妻人人看人人澡| 日韩成人av中文字幕在线观看| 大香蕉久久网| 插阴视频在线观看视频| 国产乱人视频| 国产精品久久久久久久久免| 亚洲av.av天堂| 大又大粗又爽又黄少妇毛片口| 男女边吃奶边做爰视频| 国产探花极品一区二区| 欧美成人a在线观看| 成年女人看的毛片在线观看| 亚洲人与动物交配视频| 精华霜和精华液先用哪个| 在线播放国产精品三级| 九色成人免费人妻av| 日韩成人伦理影院| 性插视频无遮挡在线免费观看| 欧美三级亚洲精品| 三级毛片av免费| 国产 一区精品| 久久精品国产鲁丝片午夜精品| 亚洲国产欧美在线一区| 亚洲无线观看免费| 18禁在线无遮挡免费观看视频| 久久久久国产网址| av在线天堂中文字幕| 中文字幕av在线有码专区| 精品久久久久久久久av| 精品久久久久久久久久免费视频| 美女国产视频在线观看| 国产激情偷乱视频一区二区| 成人毛片60女人毛片免费| 欧美xxxx性猛交bbbb| 精品人妻视频免费看| 国产乱人偷精品视频| 变态另类丝袜制服| 免费看光身美女| 国产成人影院久久av| 久久久久久九九精品二区国产| 欧美3d第一页| 高清日韩中文字幕在线| 99久久精品一区二区三区| 偷拍熟女少妇极品色| 黄色日韩在线| 熟女人妻精品中文字幕| 少妇被粗大猛烈的视频| 国产高清三级在线| 欧美性感艳星| 欧美区成人在线视频| 亚洲国产精品成人久久小说 | 亚洲天堂国产精品一区在线| 色综合站精品国产| 国产精品.久久久| 久久九九热精品免费| 亚洲精品久久久久久婷婷小说 | 美女被艹到高潮喷水动态| 观看美女的网站| 国产熟女欧美一区二区| 99热全是精品| 精品国产三级普通话版| 国产色婷婷99| 日本av手机在线免费观看| 熟女电影av网| 夜夜夜夜夜久久久久| 久久久成人免费电影| 国产精品电影一区二区三区| 一区二区三区四区激情视频 | 久久精品久久久久久久性| 久久99热这里只有精品18| 日韩大尺度精品在线看网址| 免费观看人在逋| 插逼视频在线观看| 麻豆精品久久久久久蜜桃| avwww免费| 免费人成在线观看视频色| 一夜夜www| 午夜福利在线观看免费完整高清在 | 国产精品久久久久久久电影| 在线免费观看不下载黄p国产| 久久精品夜夜夜夜夜久久蜜豆| 联通29元200g的流量卡| av视频在线观看入口| 亚洲电影在线观看av| 一级黄片播放器| 日韩制服骚丝袜av| 日韩欧美在线乱码| 亚洲精品影视一区二区三区av| 久久久久久九九精品二区国产| 人妻久久中文字幕网| 又粗又硬又长又爽又黄的视频 | 国产精品爽爽va在线观看网站| 三级经典国产精品| 三级毛片av免费| 国产伦精品一区二区三区视频9| a级一级毛片免费在线观看| 91久久精品国产一区二区成人| 成人毛片a级毛片在线播放| 女的被弄到高潮叫床怎么办| 九色成人免费人妻av| 亚洲国产精品合色在线| 国产精品一区www在线观看| 亚洲精品成人久久久久久| 国内精品久久久久精免费| 国产黄片视频在线免费观看| 国产免费一级a男人的天堂| 晚上一个人看的免费电影| 五月玫瑰六月丁香| 久久久久国产网址| 国产精品美女特级片免费视频播放器| 日韩中字成人| 啦啦啦韩国在线观看视频| 亚洲av第一区精品v没综合| 天堂中文最新版在线下载 | 日本五十路高清| 中文字幕人妻熟人妻熟丝袜美| 成人特级av手机在线观看| av天堂中文字幕网| 一级黄片播放器| 高清日韩中文字幕在线| 国产精品99久久久久久久久| 老司机影院成人| 18禁黄网站禁片免费观看直播| 国产探花在线观看一区二区| 久久韩国三级中文字幕| 国产精品久久久久久亚洲av鲁大| 中文在线观看免费www的网站| 欧美一区二区亚洲| 色哟哟哟哟哟哟| 中出人妻视频一区二区| 欧美精品国产亚洲| 亚洲色图av天堂| 一级黄片播放器| 丰满人妻一区二区三区视频av| 久久久久久久久久成人| 亚洲四区av| 久久精品国产自在天天线| 成人午夜精彩视频在线观看| 日韩欧美精品免费久久| 欧美日韩在线观看h| 亚洲三级黄色毛片| 国产精品久久久久久av不卡| 波多野结衣高清无吗| av黄色大香蕉| 美女xxoo啪啪120秒动态图| 寂寞人妻少妇视频99o| 特大巨黑吊av在线直播| 一级毛片电影观看 | 国产三级中文精品| 最近最新中文字幕大全电影3| 日本av手机在线免费观看| 欧美一区二区精品小视频在线| 直男gayav资源| 青春草亚洲视频在线观看| 国产成人精品久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 床上黄色一级片| 成人亚洲欧美一区二区av| 干丝袜人妻中文字幕| 国产成人福利小说| 欧美日韩综合久久久久久| 国产单亲对白刺激| 久久精品影院6| av天堂中文字幕网| 男人狂女人下面高潮的视频| 精品久久国产蜜桃| 九九在线视频观看精品| 99久国产av精品| 亚洲自偷自拍三级| 成年版毛片免费区| 成年免费大片在线观看| 变态另类丝袜制服| 欧美区成人在线视频| 亚洲欧美精品综合久久99| a级毛片a级免费在线| 在线天堂最新版资源| 欧美极品一区二区三区四区| 哪里可以看免费的av片| 日本成人三级电影网站| 国产v大片淫在线免费观看| 久久久久久大精品| 不卡一级毛片| 国产精品永久免费网站| 青青草视频在线视频观看| 天天躁夜夜躁狠狠久久av| 真实男女啪啪啪动态图| 久久韩国三级中文字幕|