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

    基于InVEST模型的岷江流域土地利用變化對生境質(zhì)量的影響研究

    2022-09-17 03:07:50趙慶建吳曉珍
    生態(tài)科學 2022年6期
    關鍵詞:岷江生境土地利用

    趙慶建, 吳曉珍

    基于InVEST模型的岷江流域土地利用變化對生境質(zhì)量的影響研究

    趙慶建*, 吳曉珍

    南京林業(yè)大學經(jīng)濟管理學院, 南京 210037

    土地利用及生境質(zhì)量變化過程對維護自然系統(tǒng)平衡, 區(qū)域社會經(jīng)濟可持續(xù)發(fā)展具有重大意義。借助InVEST模型生境質(zhì)量模塊, 以岷江流域作為研究區(qū), 在分析岷江流域土地利用變化的基礎上對其2000、2008、2013、2017年的生境質(zhì)量變化特征進行分析。結果表明: (1)在2000—2017年, 岷江流域的建設用地呈現(xiàn)不斷增長的趨勢, 而耕地、林地和草地呈現(xiàn)縮減的趨勢; 耕地、林地和建設用地之間的轉換較為明顯, 建設用地轉入最大的來源主要是耕地。(2)從時間尺度上看, 岷江流域生境質(zhì)量屬于較差和良好等級的占比總體上分別上升71.34%和2.41%; 中等和優(yōu)等等級的占比總體上下降了4.28%和2.61%。(3)從空間尺度上看, 岷江流域上游地區(qū)的生境質(zhì)量總體上高于中下游地區(qū); 建設用地大面積擴張, 侵占了大量的耕地和林地, 從而導致了中下游平原區(qū)域的整體生境質(zhì)量下降, 其中處于較差和中等等級的區(qū)域顯著擴大。

    岷江流域; 土地利用變化; 生境質(zhì)量; InVEST模型

    0 前言

    生境質(zhì)量是與人類密切相關的自然資源質(zhì)量, 能夠反映人類生存環(huán)境的優(yōu)劣程度, 被視為區(qū)域生物多樣性和生態(tài)服務水平的重要表征及反映, 是區(qū)域生態(tài)安全保障和人類福祉提升的關鍵環(huán)節(jié)[1-2]。土地利用通過影響生境斑塊之間的物質(zhì)流、能量流循環(huán)過程, 對區(qū)域生境質(zhì)量分布格局和功能產(chǎn)生重大影響[3]。近年來, 人口增長和社會經(jīng)濟快速發(fā)展, 使得對土地資源的需求日益增長, 自然資源被過度開發(fā)。土地利用方式及強度的變化使生態(tài)系統(tǒng)受到嚴重破壞, 進而導致自然資源的價值不斷下降[4-5]。深入分析土地利用變化對區(qū)域生境質(zhì)量的影響是實現(xiàn)資源可持續(xù)發(fā)展、制定區(qū)域生態(tài)保護政策的基礎依據(jù)。

    目前基于土地利用變化對生境質(zhì)量的影響研究主要有兩種方法。第一種是基于景觀格局的指標體系方法。第二種是基于生態(tài)評估的定量研究[6-9], 常見的評價模型有HIS模型、SoLVES模型、InVEST模型等, 其中InVEST模型中的生境質(zhì)量模塊能夠通過評估生境類型、植被類型的范圍和退化程度來反應生境質(zhì)量及生物多樣性, 并結合生態(tài)系統(tǒng)內(nèi)的土地利用植被信息和生物多樣性威脅因素來得到生境質(zhì)量的分布信息[10-11]。目前國內(nèi)學者重點運用該模型對我國部分具有重要生態(tài)價值的區(qū)域如自然保護區(qū)[12-13]、流域[14-15]、城市[16-19]的生境質(zhì)量進行評價, 如谷雙喜[20]等運用InVEST模型分析了皖江流域生境質(zhì)量時空變化特征, 發(fā)現(xiàn)該流域的土地利用類型轉移主要為耕地轉向建設用地且整體的生境質(zhì)量下降明顯; 謝余初[21]等以甘肅白龍江流域為例, 結合GIS和InVEST模型選取區(qū)域生境質(zhì)量、植物凈初級生產(chǎn)力和景觀狀態(tài)指數(shù)為評價指標, 在柵格像元尺度上開展白龍江流域景觀生物多樣性時空變化特征分析。

    岷江流域作為長江經(jīng)濟帶上游綠色發(fā)展的生態(tài)屏障, 區(qū)域內(nèi)土地利用類型豐富, 林地、草地所占比重較大, 對流域內(nèi)生境質(zhì)量、水源涵養(yǎng)等生態(tài)功能具有至關重要的作用。同時, 四川省作為西部地區(qū)重要的經(jīng)濟中心, 近年來, 其工業(yè)化、城鎮(zhèn)化、農(nóng)業(yè)現(xiàn)代化都得到了發(fā)展, 岷江流域的土地利用格局也發(fā)生了不同程度的變化。在2000—2017年間, 岷江流域土地利用總體上以草地、林地和耕地為主, 建設用地總量迅速擴張, 主要以占用耕地和草地實現(xiàn)不同程度的土地利用變化。各地類面積的變化以及相互間的轉換, 使得生境狀況發(fā)生顯著變化[16]。因此, 借助相關軟件以岷江流域為研究對象進行定量分析十分必要, 此外, 岷江流域上中下游地勢復雜, 科學研究進展緩慢, 基礎資料缺失, 利用對數(shù)據(jù)要求較小、適用性強的InVEST模型對區(qū)域生境質(zhì)量進行深入研究具有重要的現(xiàn)實意義。本文在分析岷江流域的土地利用變化的基礎上, 利用InVEST模型對其2000、2008、2013和2017年的生境質(zhì)量時空變化特征進行深入分析從而為岷江流域在上中下游開展生態(tài)保護, 推進生態(tài)示范區(qū)建設的過程中提供重要的理論支撐。

    1 研究區(qū)概況

    岷江流域全長約1279千米, 面積約135881平方公里, 地處北緯28°38′—33°10′, 東經(jīng)102°35′—104°51′。岷江可分為上中下游, 以都江堰市以上為上游, 都江堰市至樂山市大渡河匯入處為中游, 樂山至宜賓長江匯合處為下游段[22]。岷江上游地形復雜, 峽谷連綿, 地廣人稀; 中下游地勢較低, 多為丘陵平原, 人口和耕地較多。隨著人類活動范圍的不斷擴大, 流域土地利用和地表覆蓋發(fā)生了顯著變化, 導致支流河道斷流、水污染嚴重、水土流失、濕地退化等生態(tài)環(huán)境問題頻繁發(fā)生, 流域植被面積與生物多樣性都在急劇減少, 生境質(zhì)量不斷退化, 嚴重威脅到區(qū)域的可持續(xù)發(fā)展。

    2 研究數(shù)據(jù)與方法

    2.1 研究數(shù)據(jù)

    岷江流域各縣市的遙感影像數(shù)據(jù)來自地理數(shù)據(jù)空間云(http://www.gscloud.cn/search), 選取2000、2008、2013和2017年分辨率為30 m的Landsat TM遙感影像。根據(jù)《土地利用現(xiàn)狀分類標準》(GB T21010—2007), 將土地利用體系分類合并為耕地、林地、草地、水域、建設用地和未利用土地6類(圖1)。

    2.2 土地利用轉移

    土地利用轉移矩陣是對不同時期地類變化的定量描述。通過計算不同時期同種地類的面積或者百分比差值得到土地變化度, 從而反映出特定區(qū)域特定時間各個地類的面積數(shù)據(jù), 同時反映出各個地類初期的轉出面積以及末期轉入面積。具體數(shù)學表達形式為:

    圖1 2000—2017年岷江流域土地利用空間分布

    Figure 1 Spatial distribution of land use in the Minjiang River Basin in 2000-2017

    式中:S表示研究時段初期類土地到研究末期轉為類土地的面積,為土地利用類型。此外, 轉移矩陣的每一行總和表示研究初期該土地利用類型的面積總數(shù), 每個行值表示該土地類型的轉移去向和大小; 每一類的總和表示研究末期該土地類型的面積總數(shù), 每個列值表示該土地類型的所有轉入類型及大小。

    通過ArcGIS對不同時期岷江流域土地利用現(xiàn)狀圖進行融合疊置, 計算面積后導出屬性表, 從而完成轉移矩陣的數(shù)據(jù)提取。

    2.3 生境質(zhì)量模型

    InVEST模型是由美國斯坦福大學主導開發(fā)的一種可以量化多種生態(tài)服務系統(tǒng)的評估模型, 其中生境質(zhì)量模塊是利用土地利用分類圖, 基于威脅源的影響距離及空間權重, 生境適宜性以及對威脅源的敏感性來對研究區(qū)的生境質(zhì)量進行空間量化評估的, 具有引入數(shù)據(jù)量少, 輸出數(shù)據(jù)量大的特點[23]。

    2.3.1 生境退化度評價原理

    生境退化度代表了威脅因子對生境結構的影響程度, 其計算基于以下假設, 即生態(tài)系統(tǒng)中的某種土地利用類型對威脅因子的敏感度越高, 則該地類的退化度越大, 其大小與生境中各地類以及威脅因子的數(shù)量等因素密切相關。土地利用類型為時柵格的生境退化度D的計算公式如下:

    式中,為威脅因子;是威脅柵格圖的柵格數(shù);Y是威脅柵格圖中的一組威脅柵格;w為威脅因子權重, 取值范圍為0—1;r表征柵格是否為威脅柵格(0或1);i為威脅柵格的威脅因子值對區(qū)域內(nèi)柵格的威脅水平;β為柵格的可達性水平, 取值范圍0—1;K是土地利用類型為威脅因子的敏感度, 取值范圍0—1。i由以下公式計算可得:

    2.3.2 生境質(zhì)量評價原理

    生境質(zhì)量是生態(tài)環(huán)境提供給生物個體和種群生存所需要的環(huán)境水平, 是一個連續(xù)變量, 數(shù)值范圍由低到高, 生境質(zhì)量越高的斑塊生態(tài)結構和功能越穩(wěn)定。人類對土地利用的方式和強度決定了生境質(zhì)量的高低, 土地利用強度越大, 生境質(zhì)量下降越明顯。

    在InVEST生境質(zhì)量模型中, 生境質(zhì)量基于生境退化度進行計算, 分值隨生境退化分值的增大而減小, 計算公式如下:

    式中,Q表示土地利用類型為時柵格的生境質(zhì)量;H表示土地利用類型的生境適宜性;D表示土地利用類型為時柵格的生境退化程度;和為比例因子,常數(shù)為半飽和常數(shù), 通常等于值, 可使它取D最大值的一半,定義為常數(shù)2.5。

    本研究將人類活動最為集中、對生境質(zhì)量直接影響比較大的耕地、城鎮(zhèn)用地、農(nóng)村居民點和交通運輸用地4種土地利用類型定義為生境威脅因子, 參照InVEST模型使用指南, 在模型中應用最大脅迫因子對生境影響范圍(MAX-DIST), 脅迫因子權重(0—1)、脅迫因子對生境影響遞減指數(shù)(0—1)、生境對每種脅迫因子敏感度指數(shù)(0—1)。表1中參數(shù)來源于模型推薦參考值, 表2參數(shù)參考模型推薦的參考值及相關文獻[24-26]。

    3 結果與分析

    3.1 2000—2017年岷江流域土地利用變化情況

    在2000—2008年, 岷江流域的土地利用變化主要體現(xiàn)在耕地向建設用地轉變, 其中草地、建設用地、林地、水域以及未利用土地有所增長, 建設用地增長最為迅速, 增長了189.13 km2, 約增長18%, 而耕地面積下降了約243 km2, 約下降1.2%。同時, 一方面岷江流域的總人口由1330.9萬人增長至1390.5萬人, 人口增長率約為4.5%, 對住宅、工業(yè)和其他設施用地的需求越來越大, 城鎮(zhèn)村范圍不斷向外擴展; 另一方面, 隨著城市化的快速發(fā)展, 岷江流域內(nèi)的中心城區(qū)和各縣GDP總值從787.5億元增長至1793.3億元, 各地區(qū)的基礎設施和交通網(wǎng)絡的完善占用了大量耕地, 導致建設用地面積持續(xù)增加, 耕地面積大幅減少。

    表1 威脅因子屬性表

    表2 生境適宜度及其對威脅因子的敏感度

    在2008—2013年, 岷江流域的土地利用變化主要體現(xiàn)在草地和耕地向建設用地轉變, 其中草地和耕地大幅度減少, 分別減少了712 km2和340 km2, 下降幅度達到4.7%和1.8%; 建設用地、林地、水域和未利用土地不斷增長, 建設用地約增長443 km2, 增長幅度達37%。與此同時, 流域內(nèi)總人口由1346.8萬人上升至1381萬人, GDP總值由1790.3億元增長至4409.6億元, 人口的增加對土地不斷提出新的索取, 導致土地利用類型的變化顯著, 這與彭文甫等人的研究結果也相一致[27], 該期間經(jīng)濟快速發(fā)展帶動了土地城鎮(zhèn)化速度加快。

    在2013—2017年,岷江流域的土地利用變化表現(xiàn)為草地和耕地向建設用地和水域轉變, 其中草地、耕地、林地面積和未利用土地分別減少了0.06%、1.4%、0.07%和5.6%; 建設用地和水域分別增加了14.2%和8.3%。其中建設用地增長面積中由耕地轉變的約占94%; 水域增長面積中由耕地轉變的約占60.5%。相比于2008—2013年, 流域內(nèi)人口上升至1398萬人, GDP總值也達6421.1億元, 但耕地變化呈現(xiàn)小幅度下降趨勢。在此期間, 四川省制定“一帶六片八廊”的省域生態(tài)安全格局, 嚴格保護耕地, 防止城鎮(zhèn)建設連綿, 同時加強生態(tài)型經(jīng)濟發(fā)展, 協(xié)調(diào)社會經(jīng)濟發(fā)展與保護資源的關系, 逐步改善了土地生態(tài)狀況。

    總體來看, 岷江流域在2000—2017年間建設用地面積不斷擴增, 約增長了858 km2, 耕地和草地面積呈下降趨勢, 分別減少了830 km2和701 km2, 約減少了4.7%和4.3%。流域內(nèi)人地關系緊張的基本格局沒有改變, 大量耕地、草地被開墾利用, 亂砍亂伐的不合理土地利用現(xiàn)象依然嚴重, 水源涵養(yǎng)、凈化空氣等自然資源價值同步呈現(xiàn)下降趨勢, 但隨著公眾環(huán)保意識的增強以及國家退耕還林工程的啟動, 在2008—2017年耕地面積的減小幅度呈現(xiàn)下降趨勢。

    3.2 2000—2017年岷江流域生境質(zhì)量時空變化

    InVEST模型采用生境質(zhì)量指數(shù)和生境退化指數(shù)來反映研究區(qū)域的生境質(zhì)量總體狀況。其中生境質(zhì)量指數(shù)是指環(huán)境提供給適合物種和種群持續(xù)生存和繁衍所需條件和資源能力的大小, 模型中以柵格圖層上0—1的連續(xù)變化值呈現(xiàn), 值越靠近1, 生境質(zhì)量越好。通常土地利用強度的增加會引起威脅源地的增加和強度的增大, 使其附近的生境質(zhì)量退化。將準備好的數(shù)據(jù)輸入InVEST模型, 運行生境質(zhì)量評價模塊, 得到2000、2008、2013和2017年的生境質(zhì)量指數(shù)和生境退化程度。結合岷江流域的實際生境質(zhì)量情況進行等級劃分, 區(qū)分為較差、中等、良好、優(yōu)等4個等級, 對應的生境質(zhì)量指數(shù)范圍分別為0—0.2、0.2—0.5、0.5—0.8和0.8—1, 統(tǒng)計2000—2017年各個生境質(zhì)量等級的柵格數(shù)所占比例(表4)。

    表3 2000—2017年岷江流域土地利用轉移矩陣

    在時間上分析, 在2000—2017年, 生境質(zhì)量較差的流域面積一直在增加, 在2017年較差流域面積相比于2000年約增長了71.34%; 生境質(zhì)量處于中等水平的流域面積呈現(xiàn)下降趨勢, 分別占當年柵格數(shù)的35.97%、34.71%、35.35%和35.02%; 生境質(zhì)量處于良好和優(yōu)等水平的流域面積呈現(xiàn)上升、下降再上升的變化趨勢。2000—2017生境質(zhì)量變化趨勢與該階段人口增長變化、經(jīng)濟發(fā)展路徑特點相吻合, 在2000—2013年我國處于高速發(fā)展階段, 對土地資源進行過度開發(fā), 之后轉向高質(zhì)量發(fā)展, 堅持綠色發(fā)展道路, 大力實施天保、退耕還林等生態(tài)保護和建設工程, 有力促進了岷江流域生境質(zhì)量的進一步改善。

    在空間上分析, 從圖2可以看出在2000—2017年, 岷江流域整體生境質(zhì)量從岷江上游到岷江下游逐漸減弱。岷江流域上游受到地貌高度、坡度、水分、土壤等條件的限制, 林地的覆蓋面積大, 地形破碎程度比較高, 土地利用類型主要是草地和林地, 不利于進行生產(chǎn)生活開發(fā), 適合居住或耕種的土地所占的比例很小, 受到人類活動影響較小, 因此生境質(zhì)量指數(shù)相對比較高。從縣級行政區(qū)劃分來看, 松潘縣、黑水縣、茂縣、理縣、汶川縣、都江堰市、蘆山縣等地區(qū)生境質(zhì)量指數(shù)相對較高, 處于0.7以上; 岷江下游的地區(qū)大多數(shù)為平原和沿江流域, 地形破碎程度比較小, 受人為擾動影響較大, 同時隨著城市化進程的加快以及建設用地需求的不斷增加, 耕地、草地面積呈現(xiàn)不斷減少的趨勢, 其中崇州市、郫縣、邛崍市、蒲江縣、彭山縣、眉山縣、仁壽縣、樂山市、青神縣等地區(qū)生境質(zhì)量指數(shù)相對處于中等水平, 生境質(zhì)量較差的區(qū)域主要呈現(xiàn)散點狀分布, 主要集中在岷江流域東南區(qū)域內(nèi)的大中小城市。從土地利用及其轉移變化來看, 作為生境威脅源的建設用地迅速增加, 耕地、林地等生境適宜性較高的土地利用類型在不斷減少, 從而環(huán)境的破壞程度比較嚴重, 對岷江流域的生境質(zhì)量造成巨大威脅, 如仁壽縣、榮縣、井研縣等區(qū)域生境質(zhì)量指數(shù)相對較低, 處于0.3以下。

    表4 2000—2017年岷江流域生境質(zhì)量變化占比

    圖2 2000—2017年岷江流域生境質(zhì)量指數(shù)

    Figure 2 Habitat quality index of Minjiang river basin in 2000-2017

    3.3 2000—2017年岷江流域生境退化度時空變化

    對于生境退化指數(shù), 由于InVEST模型的輸出結果中直接表達的結果值具有連續(xù)性, 不易進行定性分析和空間分布規(guī)律的解讀, 因此利用ArcGIS中的自然斷點分類方法, 基于數(shù)據(jù)的分布情況采用類內(nèi)差異最小、類間差異最大的原則進行重分類和空間統(tǒng)計。將岷江流域生境退化度劃分為四個等級, 分別為基本不變, 輕度退化、中度退化和高度退化四類, 生境退化度指數(shù)的取值范圍分別對應0—0.008, 0.008—0.04, 0.04—0.07, 0.07以上。生境退化指數(shù)的高低反映了該柵格受到威脅因子的影響程度, 生境退化指數(shù)越高表明受影響越大, 退化程度越嚴重。

    在時間上分析, 從圖3可以看出在2000—2017年, 岷江流域生境質(zhì)量屬于輕度退化的變化較為明顯, 從16.38%增長至17.63%; 生境質(zhì)量基本不變和中度退化的流域面積呈先上升后下降的趨勢, 其中處于基本不變的流域面積在2013年高達43.37%; 生境質(zhì)量高度退化的流域在2000—2013年呈現(xiàn)下降的趨勢且在2017年略微增長, 區(qū)域占比分別約為17.56%、17.49%、17.16%和17.20%。

    在空間上分析, 從圖4可以看出在2000—2017年, 岷江流域生境退化度較高的地方主要分布在岷江流域的東南部, 如仁壽縣、榮縣、井研縣、樂山市、新津縣、屏山縣、沐川縣等耕地、建設用地分布廣、城市較為集中的地區(qū); 而松潘縣、黑水縣、茂縣、理縣、汶川縣、洪雅縣等地區(qū)的生境退化度幾乎為0, 這些地區(qū)海拔相對較高, 人類活動地區(qū)范圍比較小且土地利用類型主要以林地、草地為主, 生態(tài)恢復系統(tǒng)較其他區(qū)域相比較高; 生境質(zhì)量屬于中度退化的行政區(qū)主要有屏山縣、沐川縣、宜賓縣等, 這些地區(qū)的主要土地利用類型主要以耕地和水域為主, 同時從土地利用變化轉移矩陣中可以看出主要土地利用變化以耕地向草地或草地向水域轉變?yōu)橹? 受到人類活動的影響相對較小。

    4 討論

    生境質(zhì)量能夠表征區(qū)域內(nèi)生態(tài)系統(tǒng)的健康狀況, 是區(qū)域生態(tài)安全保障和人類福祉提升的關鍵。近些年, 國內(nèi)外一些學者利用生態(tài)指標法、AHP法和綜合指標法等方法對不同區(qū)域的生境質(zhì)量進行深入分析。王瓊以清河為例, 采用熵權法和主成分分析賦權法, 構建河流生境質(zhì)量評價指數(shù)[28]; 楊建強采用生態(tài)環(huán)境綜合指數(shù)法分析了萊州灣的生境質(zhì)量[29]; 吳佳俁從“社會—生產(chǎn)—生態(tài)”視角構建土地利用多功能評價指標體系, 分區(qū)研究四川省的土地利用多功能性[30]。上述研究結果準確, 但在時效性、空間化等方面有待提高。InVEST模型融合了遙感數(shù)據(jù)、數(shù)學高程模型以及氣候插值數(shù)據(jù)等對地觀測數(shù)據(jù), 克服了“以點帶面”的缺點, 具有運行參數(shù)相對較少、基礎數(shù)據(jù)獲取較易、輸出數(shù)據(jù)量大、評價結果定量化和空間可視化等特點, 使缺少實測數(shù)據(jù)資料的大中尺度區(qū)域生境質(zhì)量的定量化、可視化和精細化分析和評估成為了可能, 能夠較為直觀描述區(qū)域生境質(zhì)量的空間異質(zhì)性[31]。

    圖3 2000—2017年岷江流域生境質(zhì)量退化變化

    Figure 3 The habitat quality degradation changes in the Minjiang river basin in 2000-2017

    圖4 2000—2017年岷江流域年生境退化度空間分布圖

    Figure 4 Spatial distribution of habitat degradation in the Minjiang river basin in 2000-2017

    本研究基于InVEST模型的評估結果, 直觀準確地揭示了岷江流域在2000—2017年間生境質(zhì)量時空演變規(guī)律: 在時間上, 生境質(zhì)量較差面積呈現(xiàn)上升趨勢, 中等水平面積呈現(xiàn)下降趨勢。在此期間土地利用面積變化顯著, 建設用地面積快速增長, 耕地、林地、草地面積不斷下降, 進而導致生境質(zhì)量不斷退化。在空間上, 岷江流域上游地區(qū)的主要植被類型為草地和林地, 能夠為整個岷江流域提供涵養(yǎng)水源和凈化水質(zhì)的作用, 人為干擾比較小, 生境質(zhì)量處于優(yōu)等狀態(tài)。岷江中下游以耕地為主要植被類型, 人口密度大, 耕作、施肥、生活污水和工業(yè)廢水對岷江流域的生境質(zhì)量造成很大的破壞, 部分地區(qū)甚至陷入通過開閘放水也難以解決耕作生活用水的困境。因此, 為維護岷江流域生物多樣性, 必須科學合法管理開發(fā)利用土地, 加強對草地和耕地的建設保護, 保證生態(tài)用地的適宜面積并注重土地結構組成和空間配置, 提高生態(tài)用地的整體質(zhì)量, 從而提升流域內(nèi)的整體生境質(zhì)量。

    但InVEST模型在生境質(zhì)量評價中也具有局限性。一是模型中部分參數(shù)如土地利用類型對各個脅迫因子的敏感度是參照文獻以及InVEST模型用戶指南確定的, 而生境質(zhì)量的時空變化受到多種因素的影響, 分析結果存在一定的誤差, 只限于對岷江流域總體上的生境質(zhì)量時空變化進行分析。二是模型中土地利用類型上的所有威脅都是簡單相加, 而多種威脅的總影響總是遠遠大于各種影響的算術和[32]。因此, 未來如何科學綜合流域內(nèi)多種生態(tài)服務, 深入分析不同生態(tài)系統(tǒng)服務之間的復雜關系, 充分結合當?shù)貙崪y從而獲得更加精準的計算服務量是研究生境質(zhì)量的重點方向, 以便能更好地揭示岷江流域生境質(zhì)量的發(fā)展演變規(guī)律。

    5 結論

    本文采用InVEST模型的生境質(zhì)量評估模塊對岷江流域2000、2008、2013和2017年四個時期的生境質(zhì)量變化進行空間異質(zhì)性分析, 得出以下結論:

    (1)在2000—2017年, 岷江流域的建設用地呈現(xiàn)不斷增長的趨勢, 耕地、林地和草地呈現(xiàn)縮減的趨勢。其中耕地、草地、林地和建設用地之間的轉換較為明顯且建設用地轉入最大的來源主要是耕地。

    (2)從時間上來看, 岷江流域生境質(zhì)量處于較差水平呈增長趨勢, 中等水平呈下降趨勢; 處于良好和優(yōu)等水平的流域面積呈上升, 后小幅度減少再大幅度增加的趨勢。岷江流域的生境質(zhì)量與土地利用之間存在很顯著的相關性, 耕地、林地和草地的生境質(zhì)量高, 而建設用地的生境質(zhì)量較低且耕地和草地面積下降是生境質(zhì)量下降的主要原因。

    (3)在空間上來看, 岷江流域上游地區(qū)生境質(zhì)量總體普遍高于中下部地區(qū), 生境質(zhì)量較差的區(qū)域大多集中在耕地、建設用地聚集的人類活動區(qū)域。

    [1] 劉園, 周勇, 杜越天. 基于InVEST模型的長江中游經(jīng)濟帶生境質(zhì)量的時空分異特征及其地形梯度效應[J]. 長江流域資源與環(huán)境, 2019, 28(10): 2429–2440.

    [2] Jason B F, ERAN H, WILLIAM D, et al. Stream Physical Characteristics Impact Habitat Quality for Pacific Salmon in Two Temperate Coastal Watersheds[J]. PLoSONE, 2015, 10(7): 1–16.

    [3] 吳健生, 曹祺文, 石淑芹, 等. 基于土地利用變化的京津冀生境質(zhì)量時空演變[J]. 應用生態(tài)學報, 2015, 26(11): 3457–3466.

    [4] 張暉, 吳霜, 張燕媛, 等. 流域生態(tài)補償政策對受償?shù)貐^(qū)經(jīng)濟增長的影響研究——以安徽省黃山市為例[J]. 長江流域資源與環(huán)境, 2019, 28(12): 2848–2856.

    [5] 李敏納, 蔡舒, 覃成林. 中美西部土地資源開發(fā)異同比較及啟示[J]. 西部論壇, 2016, 2(4): 47–57.

    [6] 王燕, 高吉喜, 金宇, 等. 基于2005—2015年土地利用變化和InVEST模型的內(nèi)蒙古巴林右旗農(nóng)牧交錯帶生境質(zhì)量研究[J]. 生態(tài)與農(nóng)村環(huán)境學報, 2020, 36(5): 654– 662.

    [7] ZHONG Lina, WANG Jun. Evaluation on effect of land consolidation on habitat quality based on InVEST model[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2017, 33(1): 250–255.

    [8] MENESES B M, REIS R, VALE M J, et al. Land use and land cover changes in Zêzere watershed (Portugal) - water quality implications[J]. Science of the Total Environment, 2015, 527(9): 439–447.

    [9] 焦云騰, 程文仕. 基于InVEST模型的土地整治對生境質(zhì)量的影響研究[J]. 國土與自然資源研究, 2019, 42(2): 39–42.

    [10] 黃從紅, 楊軍, 張文娟. 生態(tài)系統(tǒng)服務功能評估模型研究進展[J]. 生態(tài)學雜志, 2013, 32(12): 3360–3367.

    [11] CERIAN GIBBES, DAVID G. HAVLICK, JOSEPH R. ROBB. Land use and land cover in a transitioning militarized landscape[J]. Journal of Land Use Science, 2017, 12(23): 182–196.

    [12] 顧羊羊, 黃賢峰, 鄒長新, 等. 沅江源自然保護區(qū)生境質(zhì)量變化遙感監(jiān)測[J]. 生態(tài)與農(nóng)村環(huán)境學報, 2019, 35(6): 764–772.

    [13] 朱燕, 李怡然, 李雪梅. 基于InVEST模型的昌黎黃金海岸國家級自然保護區(qū)生境質(zhì)量評價[J]. 環(huán)境與可持續(xù)發(fā)展, 2019, 44(6): 156–160.

    [14] 任涵, 張靜靜, 朱文博, 等. 太行山淇河流域土地利用變化對生境的影響[J]. 地理科學進展, 2018, 37(12): 1693– 1704.

    [15] 劉園, 周勇, 杜越天. 基于InVEST模型的長江中游經(jīng)濟帶生境質(zhì)量的時空分異特征及其地形梯度效應[J]. 長江流域資源與環(huán)境, 2019, 28(10): 2429–2440.

    [16] 楊志鵬, 許嘉巍, 馮興華, 等. 基于InVEST模型的東北地區(qū)土地利用變化對生境的影響研究[J]. 生態(tài)科學, 2018, 37(6): 139–147.

    [17] 張文靜, 孫小銀, 單瑞峰. 基于InVEST模型研究山東半島沿海地區(qū)土地利用變化及其對生境質(zhì)量的影響[J]. 環(huán)境生態(tài)學, 2019, 1(5): 15–23.

    [18] ZHAO Qingjian, WEN Zuomin, CHEN Shulin, et al. Quantifying Land Use/Land Cover and Landscape Pattern Changes and Impacts on Ecosystem Services[J]. International journal of environmental research and public health, 2019, 17(1): 119–126.

    [19] 陳書林, 虞宋杰. 長三角地區(qū)1985~2015年土地利用變化及其驅動力分析[J]. 林業(yè)經(jīng)濟, 2019, 41(11): 49–54.

    [20] 谷雙喜, 楊海燕, 李衡均, 等. 基于土地利用的皖江流域生境質(zhì)量時空變化特征分析[J]. 黑龍江工程學院學報, 2019, 33(2): 1–8.

    [21] 謝余初, 鞏杰, 張素欣, 等. 基于遙感和InVEST模型的白龍江流域景觀生物多樣性時空格局研究[J]. 地理科學, 2018, 38(6): 979–986.

    [22] 杜華明. 岷江流域徑流演變及其對氣候變化的響應[J]. 宜賓學院學報, 2014, 14(12): 36–40.

    [23] 黃從紅, 楊軍, 張文娟. 生態(tài)系統(tǒng)服務功能評估模型研究進展[J]. 生態(tài)學雜志, 2013, 32(12): 3360–3367.

    [24] 陳珊珊, 臧淑英, 孫麗. 基于InVEST模型的土地利用變化對生境質(zhì)量的影響研究——以松嫩平原為例[J]. 中國農(nóng)學通報, 2020, 36(3): 74–80.

    [25] 陳妍, 喬飛, 江磊. 基于InVEST模型的土地利用格局變化對區(qū)域尺度生境質(zhì)量的影響研究——以北京為例[J]. 北京大學學報(自然科學版), 2016, 52(3): 553–562.

    [26] CROOKS K R, BURDETT C L, THEOBALD D M, et al. Quantification of habitat fragmentation reveals extinction risk in terrestrial mammals[J]. Proceedings of the National Academy of Sciences, 2017, 114(29): 7635–7640.

    [27] 彭文甫, 周介銘, 楊存建, 等. 基于土地利用變化的四川省生態(tài)系統(tǒng)服務價值研究[J]. 長江流域資源與環(huán)境, 2014, 23(7): 1053–1062.

    [28] 王瓊, 盧聰, 李法云, 等. 基于主成分分析和熵權法的河流生境質(zhì)量評價方法——以清河為例[J]. 生態(tài)科學, 2017, 36(4): 185–193.

    [29] 楊建強, 朱永貴, 宋文鵬, 等. 基于生境質(zhì)量和生態(tài)響應的萊州灣生態(tài)環(huán)境質(zhì)量評價[J]. 生態(tài)學報, 2014, 34(1): 105–114.

    [30] 吳佳俁. 四川省縣域土地利用多功能性評價及分區(qū)研究[D]. 四川師范大學, 2018, 12(3): 71–141.

    [31] 李國慶, 劉長成, 劉玉國, 等. 物種分布模型理論研究進展[J]. 生態(tài)學報, 2013, 33(16): 4827–4835.

    [32] BAGSTAD K J, SEMMENS D J, WAAGE S, et al. A comparative assessment of decision-support tools for ecosystem services quantification and valuation[J]. Ecosystem services 2013, 5: 27–39.

    Research on the impact of land use change on habitat quality in Minjiang river basin based on InVEST Model

    ZHAO Qingjian*, WU Xiaozhen

    College of Economics and Management, Nanjing Forestry University, Nanjing 210037, China

    The process of land use and habitat quality change is of great significance to maintaining the balance of natural systems and the sustainable development of regional socio-economics. Taking the Minjiang River as the research area, the InVEST Model was used to evaluate the spatial and temporal changes in habitat quality based on the land use change. The results show that: (1) The construction land in the Minjiang River Basin showed an increasing trend, while the cultivated land, woodland and grassland showed a shrinking trend in 2000-2017; the conversion between cultivated land, woodland and construction land was more obvious, and the construction land changed. The largest source of construction land transfer was mainly cultivated land. (2) In terms of time scale, the proportion of poor and good habitat quality increased by 71.34% and 2.41% respectively; the proportion of medium and excellent grades dropped by 4.28% and 2.61% respectively. (3) In terms of spatial scale, the habitat quality of upper reaches of the Minjiang River Basin was generally higher than that of the middle and lower reaches. The large-scale expansion of construction land had encroached on a large amount of arable land and woodland, which led to a decline in the overall habitat quality of the middle and lower reaches of the plain area. Among them, the poor and middle-level areas were significantly expanded.

    Minjiang river basin; land use change; habitat quality; InVEST Model

    10.14108/j.cnki.1008-8873.2022.06.001

    Q148

    A

    1008-8873(2022)06-001-10

    2020-09-10;

    2020-11-25

    國家自然科學基金項目(71373125); 江蘇高校哲學社會科學研究重大項目(2020SJZDA073)

    吳曉珍(1996—), 女, 江蘇南通人, 碩士研究生, 主要從事生態(tài)經(jīng)濟研究, E-mail: 1026742937@qq.com

    通信作者:趙慶建, 男, 博士, 教授, 主要從事生態(tài)經(jīng)濟?復雜社會生態(tài)系統(tǒng)研究, E-mail: zhao5190@126.com

    趙慶建, 吳曉珍. 基于InVEST模型的岷江流域土地利用變化對生境質(zhì)量的影響研究[J]. 生態(tài)科學, 2022, 41(6): 1–10.

    ZHAO Qingjian, WU Xiaozhen. Research on the impact of land use change on habitat quality in Minjiang river basin based on InVEST Model[J]. Ecological Science, 2022, 41(6): 1–10.

    猜你喜歡
    岷江生境土地利用
    岷江
    棗樹適應干旱生境研究進展
    落葉果樹(2021年6期)2021-02-12 01:29:10
    岷江行(外一首)
    岷峨詩稿(2019年4期)2019-04-20 09:02:06
    瘋狂的蟲草,瘋狂的松茸和瘋狂的岷江柏*——專訪作家阿來
    阿來研究(2018年1期)2018-07-10 09:39:08
    土地利用生態(tài)系統(tǒng)服務研究進展及啟示
    岷江同大渡河相會樂山
    廉政瞭望(2016年13期)2016-08-11 11:22:02
    金秀瑤族長鼓舞傳承的文化生境
    體育科技(2016年2期)2016-02-28 17:06:07
    濱??h土地利用挖潛方向在哪里
    不同生境中入侵植物喜旱蓮子草與本地種接骨草的光合特性
    論低碳經(jīng)濟與轉變土地利用方式
    久久精品久久久久久久性| 一级av片app| 国产黄色视频一区二区在线观看 | 精品久久久噜噜| 亚洲人成网站高清观看| 欧美成人午夜免费资源| 亚洲精品乱码久久久v下载方式| 亚洲久久久久久中文字幕| 精品久久久久久成人av| 国内精品一区二区在线观看| 特级一级黄色大片| 高清视频免费观看一区二区 | 欧美+日韩+精品| 麻豆成人av视频| 国产成人一区二区在线| 亚洲无线观看免费| 特级一级黄色大片| h日本视频在线播放| av天堂中文字幕网| 高清午夜精品一区二区三区| 久久精品91蜜桃| ponron亚洲| 亚洲国产高清在线一区二区三| 免费观看在线日韩| 国产免费一级a男人的天堂| 国内少妇人妻偷人精品xxx网站| 免费观看在线日韩| 日韩av不卡免费在线播放| 亚洲精品,欧美精品| 国产不卡一卡二| 黑人高潮一二区| 乱人视频在线观看| 欧美zozozo另类| 欧美成人一区二区免费高清观看| 亚洲欧美中文字幕日韩二区| 热99re8久久精品国产| 久久久久久久国产电影| 国产综合懂色| 久久精品熟女亚洲av麻豆精品 | 丝袜喷水一区| 91精品一卡2卡3卡4卡| 亚洲成人av在线免费| 日韩精品青青久久久久久| 久久这里只有精品中国| 日韩大片免费观看网站 | 国产精品美女特级片免费视频播放器| 99热这里只有是精品在线观看| 午夜福利成人在线免费观看| 成人一区二区视频在线观看| 亚洲精品自拍成人| 变态另类丝袜制服| 在线观看美女被高潮喷水网站| 在线观看美女被高潮喷水网站| 一级毛片久久久久久久久女| 亚洲国产精品专区欧美| 99热网站在线观看| 久久精品夜夜夜夜夜久久蜜豆| 你懂的网址亚洲精品在线观看 | 日韩人妻高清精品专区| 久久久久久久久久久免费av| 国产精品人妻久久久久久| 好男人在线观看高清免费视频| 久久99热6这里只有精品| 狂野欧美白嫩少妇大欣赏| 男插女下体视频免费在线播放| 日日摸夜夜添夜夜爱| 我要看日韩黄色一级片| 亚洲精品色激情综合| 日韩欧美精品v在线| av播播在线观看一区| 亚洲精品日韩av片在线观看| 寂寞人妻少妇视频99o| 日本免费在线观看一区| 亚洲精品日韩在线中文字幕| 神马国产精品三级电影在线观看| 日本五十路高清| 亚洲自拍偷在线| 亚洲成人精品中文字幕电影| 国产一区二区三区av在线| 97人妻精品一区二区三区麻豆| 午夜a级毛片| 精品久久久久久电影网 | 在线免费观看不下载黄p国产| 三级国产精品片| 99在线视频只有这里精品首页| 能在线免费观看的黄片| 最近2019中文字幕mv第一页| 国产在视频线精品| 亚洲国产精品成人综合色| 一个人观看的视频www高清免费观看| 最近手机中文字幕大全| 国产亚洲精品久久久com| 国模一区二区三区四区视频| 男插女下体视频免费在线播放| 2021少妇久久久久久久久久久| 久久久久网色| 在线免费观看不下载黄p国产| 一级黄片播放器| 成人二区视频| 直男gayav资源| 又爽又黄a免费视频| 国产综合懂色| 免费观看精品视频网站| 天堂网av新在线| 黄色欧美视频在线观看| 成人综合一区亚洲| 午夜亚洲福利在线播放| 综合色丁香网| 日日干狠狠操夜夜爽| 搡老妇女老女人老熟妇| 亚洲五月天丁香| 久久国内精品自在自线图片| 1024手机看黄色片| 亚洲av不卡在线观看| 蜜桃亚洲精品一区二区三区| 精品99又大又爽又粗少妇毛片| 精品久久久久久久久久久久久| 精品人妻一区二区三区麻豆| 一级爰片在线观看| 一级毛片aaaaaa免费看小| 免费av观看视频| 国产精品一区二区三区四区免费观看| www.av在线官网国产| 日韩,欧美,国产一区二区三区 | 国产国拍精品亚洲av在线观看| 日本猛色少妇xxxxx猛交久久| 久久久久久大精品| 中国国产av一级| 久久精品久久精品一区二区三区| 欧美成人一区二区免费高清观看| 亚洲四区av| 最近的中文字幕免费完整| 国产单亲对白刺激| 九九爱精品视频在线观看| 日韩欧美国产在线观看| 久久精品熟女亚洲av麻豆精品 | 国产麻豆成人av免费视频| 国产老妇伦熟女老妇高清| 最近2019中文字幕mv第一页| 麻豆精品久久久久久蜜桃| 特大巨黑吊av在线直播| 亚洲不卡免费看| 国产精品永久免费网站| 欧美潮喷喷水| 免费av观看视频| 波野结衣二区三区在线| av在线蜜桃| 哪个播放器可以免费观看大片| 国产精品久久久久久精品电影| 日韩欧美三级三区| 麻豆国产97在线/欧美| 中文字幕熟女人妻在线| 国产成人a∨麻豆精品| 久久久久网色| 如何舔出高潮| 久久久a久久爽久久v久久| 国产欧美日韩精品一区二区| 亚洲成人av在线免费| 欧美精品一区二区大全| 日本与韩国留学比较| 国产成人91sexporn| 国产亚洲午夜精品一区二区久久 | 日韩欧美国产在线观看| 精品久久久久久久人妻蜜臀av| 岛国毛片在线播放| 欧美精品一区二区大全| 久久精品夜夜夜夜夜久久蜜豆| av在线天堂中文字幕| 久久午夜福利片| 久久精品夜色国产| 亚洲美女搞黄在线观看| 成人欧美大片| 黄色一级大片看看| 国产精品久久久久久久电影| 天天躁夜夜躁狠狠久久av| 午夜a级毛片| 国产一区二区在线av高清观看| 亚洲aⅴ乱码一区二区在线播放| 春色校园在线视频观看| 久久精品久久久久久久性| 国产亚洲av片在线观看秒播厂 | 一区二区三区乱码不卡18| 亚洲成人久久爱视频| 欧美另类亚洲清纯唯美| 国产成人精品一,二区| 国产探花极品一区二区| 国产熟女欧美一区二区| 2021天堂中文幕一二区在线观| 九九久久精品国产亚洲av麻豆| 久久久久久久午夜电影| 一二三四中文在线观看免费高清| 国产亚洲av嫩草精品影院| 三级国产精品片| 在线播放国产精品三级| 97超碰精品成人国产| 白带黄色成豆腐渣| 成人毛片a级毛片在线播放| 欧美日韩精品成人综合77777| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品女同一区二区软件| 18禁裸乳无遮挡免费网站照片| 免费在线观看成人毛片| 老师上课跳d突然被开到最大视频| 国产在线男女| 成人亚洲欧美一区二区av| 国产精品福利在线免费观看| 国产精品av视频在线免费观看| 久久久久久大精品| 少妇丰满av| 免费观看人在逋| 干丝袜人妻中文字幕| 国产午夜精品论理片| 亚洲真实伦在线观看| 亚洲av福利一区| 成人午夜高清在线视频| 一个人观看的视频www高清免费观看| 禁无遮挡网站| 热99re8久久精品国产| 国产亚洲av片在线观看秒播厂 | 水蜜桃什么品种好| 蜜桃亚洲精品一区二区三区| 九九爱精品视频在线观看| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| 国产精品.久久久| 日本免费一区二区三区高清不卡| 亚洲欧美日韩无卡精品| 在线观看一区二区三区| 色哟哟·www| 国模一区二区三区四区视频| 亚洲成人久久爱视频| 日本三级黄在线观看| 久久久欧美国产精品| .国产精品久久| 国产三级在线视频| 爱豆传媒免费全集在线观看| 免费大片18禁| 免费观看的影片在线观看| 久久人人爽人人片av| 在线播放国产精品三级| 国产精品不卡视频一区二区| 亚洲成人av在线免费| 国产精品国产三级专区第一集| 免费av毛片视频| 欧美一区二区亚洲| 免费观看的影片在线观看| 三级国产精品片| 国产又黄又爽又无遮挡在线| 欧美人与善性xxx| 亚洲精品aⅴ在线观看| 久久99蜜桃精品久久| 国产伦精品一区二区三区视频9| or卡值多少钱| 丰满少妇做爰视频| a级毛色黄片| 久久人妻av系列| 黄色配什么色好看| 日韩av在线免费看完整版不卡| 色吧在线观看| 久久综合国产亚洲精品| 亚洲四区av| 九九爱精品视频在线观看| 亚洲性久久影院| 22中文网久久字幕| 99久久中文字幕三级久久日本| 一二三四中文在线观看免费高清| 一级黄片播放器| 中文字幕久久专区| 成人高潮视频无遮挡免费网站| 日韩精品青青久久久久久| 国产黄a三级三级三级人| 搡老妇女老女人老熟妇| 国产真实乱freesex| av黄色大香蕉| 亚洲国产精品合色在线| 人人妻人人澡欧美一区二区| 国产高清国产精品国产三级 | 你懂的网址亚洲精品在线观看 | 国产黄色视频一区二区在线观看 | 国产av码专区亚洲av| 久久国产乱子免费精品| 国产视频首页在线观看| 简卡轻食公司| 国产免费又黄又爽又色| 小蜜桃在线观看免费完整版高清| 欧美三级亚洲精品| АⅤ资源中文在线天堂| 禁无遮挡网站| 级片在线观看| 日韩,欧美,国产一区二区三区 | 少妇猛男粗大的猛烈进出视频 | a级一级毛片免费在线观看| 最近中文字幕2019免费版| 深夜a级毛片| 免费大片18禁| 欧美xxxx性猛交bbbb| 午夜久久久久精精品| 亚洲成色77777| 国产一区二区在线观看日韩| 中文字幕熟女人妻在线| 久99久视频精品免费| av在线老鸭窝| 人体艺术视频欧美日本| 久久精品夜夜夜夜夜久久蜜豆| 欧美一区二区亚洲| 国产乱来视频区| 日韩一区二区视频免费看| 夫妻性生交免费视频一级片| 免费观看a级毛片全部| 欧美性感艳星| 中文字幕精品亚洲无线码一区| 成人欧美大片| 亚洲色图av天堂| 日本与韩国留学比较| 欧美性猛交黑人性爽| 亚洲欧美日韩高清专用| 国产精品一二三区在线看| 国产欧美日韩精品一区二区| 亚洲中文字幕日韩| 亚洲综合色惰| av.在线天堂| 久久久久精品久久久久真实原创| av在线播放精品| 日韩 亚洲 欧美在线| 十八禁国产超污无遮挡网站| 国产av码专区亚洲av| 草草在线视频免费看| 免费搜索国产男女视频| 精品99又大又爽又粗少妇毛片| 岛国毛片在线播放| 亚洲成人中文字幕在线播放| 日本免费a在线| 久久精品夜色国产| 成人综合一区亚洲| 中文亚洲av片在线观看爽| 26uuu在线亚洲综合色| 最近中文字幕高清免费大全6| a级毛色黄片| 成人午夜精彩视频在线观看| 身体一侧抽搐| 五月玫瑰六月丁香| 亚洲va在线va天堂va国产| 欧美成人精品欧美一级黄| 少妇高潮的动态图| 成人特级av手机在线观看| 日本三级黄在线观看| 久久久午夜欧美精品| 亚洲国产最新在线播放| 国产精品乱码一区二三区的特点| 一边亲一边摸免费视频| 精品人妻一区二区三区麻豆| 国产一区二区三区av在线| 欧美一区二区亚洲| 亚洲av电影不卡..在线观看| 日本wwww免费看| 国产精品美女特级片免费视频播放器| 国产精品.久久久| 久久婷婷人人爽人人干人人爱| 欧美成人一区二区免费高清观看| 婷婷色av中文字幕| 亚洲最大成人中文| 我的女老师完整版在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲成人av在线免费| 国产高清国产精品国产三级 | 99久久九九国产精品国产免费| 亚洲伊人久久精品综合 | 69人妻影院| 人妻系列 视频| 欧美潮喷喷水| 久久久久久九九精品二区国产| 精华霜和精华液先用哪个| 日本猛色少妇xxxxx猛交久久| 亚洲av免费高清在线观看| 观看免费一级毛片| 全区人妻精品视频| 国产成年人精品一区二区| 男插女下体视频免费在线播放| 99热这里只有精品一区| 国产精品乱码一区二三区的特点| av免费观看日本| av视频在线观看入口| 成年免费大片在线观看| 国产伦在线观看视频一区| 深夜a级毛片| 中文字幕熟女人妻在线| 欧美xxxx性猛交bbbb| 看黄色毛片网站| 婷婷色麻豆天堂久久 | 国产在视频线精品| 国产精品伦人一区二区| 最近最新中文字幕大全电影3| 国产黄片美女视频| 亚洲国产精品专区欧美| 精品久久久久久久末码| 内地一区二区视频在线| 老司机影院毛片| 精品人妻视频免费看| 天天躁日日操中文字幕| av线在线观看网站| 小说图片视频综合网站| www.色视频.com| 日韩av在线大香蕉| 狂野欧美白嫩少妇大欣赏| 日韩视频在线欧美| 中文天堂在线官网| 综合色av麻豆| 成年免费大片在线观看| 久久久久国产网址| 色综合亚洲欧美另类图片| 禁无遮挡网站| 天堂网av新在线| 欧美一区二区精品小视频在线| or卡值多少钱| 久热久热在线精品观看| 精品国内亚洲2022精品成人| ponron亚洲| 亚洲熟妇中文字幕五十中出| 久99久视频精品免费| 中国美白少妇内射xxxbb| 欧美日韩在线观看h| 丝袜美腿在线中文| 国产v大片淫在线免费观看| 国产探花极品一区二区| 2022亚洲国产成人精品| 亚洲美女视频黄频| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久久免| 精品久久久久久久人妻蜜臀av| 黄片无遮挡物在线观看| 亚洲av熟女| 国产大屁股一区二区在线视频| 国产免费一级a男人的天堂| 久久久精品94久久精品| 男人狂女人下面高潮的视频| 99视频精品全部免费 在线| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美国产在线观看| 2022亚洲国产成人精品| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂 | 午夜福利视频1000在线观看| 午夜免费男女啪啪视频观看| 日本午夜av视频| 亚洲aⅴ乱码一区二区在线播放| 观看免费一级毛片| 汤姆久久久久久久影院中文字幕 | 国产老妇伦熟女老妇高清| 精品国内亚洲2022精品成人| 免费看a级黄色片| 91在线精品国自产拍蜜月| 久久精品国产亚洲av天美| 成年女人看的毛片在线观看| 高清av免费在线| 日韩一区二区三区影片| h日本视频在线播放| 国产欧美另类精品又又久久亚洲欧美| 1024手机看黄色片| 99久国产av精品| 国产高清不卡午夜福利| 精品久久久久久久久av| 久热久热在线精品观看| 久久99蜜桃精品久久| 99热6这里只有精品| 两个人的视频大全免费| 2021少妇久久久久久久久久久| 一级黄色大片毛片| 一二三四中文在线观看免费高清| 亚洲成av人片在线播放无| 亚洲国产欧洲综合997久久,| 久久人妻av系列| 亚洲av免费在线观看| 天堂网av新在线| 亚洲精品乱久久久久久| 精品99又大又爽又粗少妇毛片| 男人舔女人下体高潮全视频| 欧美成人精品欧美一级黄| 2022亚洲国产成人精品| 99热全是精品| 久久欧美精品欧美久久欧美| 熟女人妻精品中文字幕| 日韩精品青青久久久久久| 能在线免费观看的黄片| 国产爱豆传媒在线观看| 热99re8久久精品国产| 国产精品1区2区在线观看.| 97在线视频观看| 成年av动漫网址| 日本免费在线观看一区| 噜噜噜噜噜久久久久久91| 国产精品久久久久久精品电影小说 | 国产在线一区二区三区精 | 男女啪啪激烈高潮av片| 亚洲精品色激情综合| 国产人妻一区二区三区在| 一个人看视频在线观看www免费| 日本免费一区二区三区高清不卡| 一级二级三级毛片免费看| 久久久久久久久久久丰满| 黄色配什么色好看| 有码 亚洲区| 高清毛片免费看| 欧美精品一区二区大全| 99久久无色码亚洲精品果冻| 国产成人午夜福利电影在线观看| 黄片wwwwww| 国产真实乱freesex| 国产白丝娇喘喷水9色精品| 亚洲内射少妇av| 久久久国产成人精品二区| 精品国产一区二区三区久久久樱花 | 欧美极品一区二区三区四区| 亚洲成色77777| 亚洲欧洲日产国产| 成年免费大片在线观看| 小说图片视频综合网站| 校园人妻丝袜中文字幕| av线在线观看网站| 国产极品天堂在线| 美女内射精品一级片tv| 免费一级毛片在线播放高清视频| 男女下面进入的视频免费午夜| 亚洲四区av| 日韩国内少妇激情av| 国产精品,欧美在线| 51国产日韩欧美| 国产精品国产高清国产av| 精品久久久久久久久av| 成年av动漫网址| 久久久久久久久久久免费av| 国产视频首页在线观看| 国产黄色视频一区二区在线观看 | 久久精品国产自在天天线| 在现免费观看毛片| 亚洲欧美日韩东京热| 国产黄a三级三级三级人| 亚洲国产成人一精品久久久| 国产白丝娇喘喷水9色精品| 激情 狠狠 欧美| 男人狂女人下面高潮的视频| 成人毛片60女人毛片免费| 卡戴珊不雅视频在线播放| 久久久久精品久久久久真实原创| 黄片wwwwww| 熟女人妻精品中文字幕| 在线观看av片永久免费下载| 51国产日韩欧美| 日韩精品青青久久久久久| 少妇熟女欧美另类| 日韩av不卡免费在线播放| 午夜激情欧美在线| av在线观看视频网站免费| 国产精品伦人一区二区| 亚洲怡红院男人天堂| 在线播放国产精品三级| 夫妻性生交免费视频一级片| 变态另类丝袜制服| 午夜激情福利司机影院| 亚洲真实伦在线观看| 精品午夜福利在线看| 欧美极品一区二区三区四区| 国产精品无大码| 久久6这里有精品| 国产精品麻豆人妻色哟哟久久 | 日本猛色少妇xxxxx猛交久久| 日日干狠狠操夜夜爽| 国产精品不卡视频一区二区| 亚州av有码| 国产久久久一区二区三区| 国产精品一区www在线观看| 一个人看视频在线观看www免费| 色综合亚洲欧美另类图片| 久久久精品欧美日韩精品| 国产精品永久免费网站| 日韩欧美在线乱码| 国产精品一区二区性色av| 黑人高潮一二区| 亚洲综合精品二区| 中文字幕久久专区| 婷婷色综合大香蕉| 青春草国产在线视频| 只有这里有精品99| 成人av在线播放网站| 美女脱内裤让男人舔精品视频| 亚洲欧洲日产国产| 中文字幕av在线有码专区| 国产男人的电影天堂91| 内地一区二区视频在线| 日韩在线高清观看一区二区三区| 2022亚洲国产成人精品| 欧美丝袜亚洲另类| 女人十人毛片免费观看3o分钟| 国产毛片a区久久久久| 亚洲精品国产av成人精品| 美女脱内裤让男人舔精品视频| av女优亚洲男人天堂| 国产午夜精品久久久久久一区二区三区| 精品99又大又爽又粗少妇毛片| 婷婷色综合大香蕉| 久久精品国产鲁丝片午夜精品| 亚洲欧洲国产日韩| 狂野欧美白嫩少妇大欣赏| 精品人妻视频免费看| 最近最新中文字幕大全电影3| 成年女人永久免费观看视频| 日韩在线高清观看一区二区三区| 成人午夜高清在线视频| 高清午夜精品一区二区三区| 水蜜桃什么品种好| 亚洲国产精品sss在线观看| 亚洲国产精品国产精品| АⅤ资源中文在线天堂| 尾随美女入室| 亚洲高清免费不卡视频| 永久免费av网站大全|