• <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一区二区精品久久| a级片在线免费高清观看视频| 午夜成年电影在线免费观看| 亚洲 国产 在线| 男女做爰动态图高潮gif福利片 | 51午夜福利影视在线观看| 日日爽夜夜爽网站| 男女做爰动态图高潮gif福利片 | 国产亚洲精品一区二区www| 搡老岳熟女国产| 久久精品91无色码中文字幕| 9热在线视频观看99| 人妻丰满熟妇av一区二区三区| a在线观看视频网站| 亚洲av熟女| 五月开心婷婷网| 乱人伦中国视频| 韩国av一区二区三区四区| 国产97色在线日韩免费| 长腿黑丝高跟| 日韩精品青青久久久久久| 久久久久久免费高清国产稀缺| 精品久久久精品久久久| 欧美乱码精品一区二区三区| 欧美日韩瑟瑟在线播放| 天堂俺去俺来也www色官网| 好男人电影高清在线观看| 欧美日韩黄片免| 亚洲精品成人av观看孕妇| 免费av中文字幕在线| 神马国产精品三级电影在线观看 | 久久久久国产精品人妻aⅴ院| 国产麻豆69| 亚洲成国产人片在线观看| 欧美一级毛片孕妇| 国产麻豆69| 精品一区二区三区四区五区乱码| 欧美丝袜亚洲另类 | 激情视频va一区二区三区| 欧美精品亚洲一区二区| 亚洲中文字幕日韩| 国产亚洲av高清不卡| 老熟妇仑乱视频hdxx| 国产高清激情床上av| av视频免费观看在线观看| 女生性感内裤真人,穿戴方法视频| 在线看a的网站| 国产av又大| 黄片大片在线免费观看| 18禁黄网站禁片午夜丰满| av片东京热男人的天堂| 亚洲成a人片在线一区二区| 最新美女视频免费是黄的| 亚洲av成人一区二区三| 欧美日韩瑟瑟在线播放| 黄网站色视频无遮挡免费观看| 久久久久亚洲av毛片大全| 国产精品影院久久| 男女午夜视频在线观看| 久热这里只有精品99| 久久久精品欧美日韩精品| 午夜福利在线观看吧| 这个男人来自地球电影免费观看| 久久久久久大精品| 精品久久久久久电影网| 黑人操中国人逼视频| 天天躁夜夜躁狠狠躁躁| 欧美日韩亚洲综合一区二区三区_| aaaaa片日本免费| 精品高清国产在线一区| 日韩成人在线观看一区二区三区| 久久中文字幕一级| 欧美日韩亚洲国产一区二区在线观看| 亚洲五月色婷婷综合| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 怎么达到女性高潮| 精品人妻1区二区| 黑丝袜美女国产一区| 精品无人区乱码1区二区| 亚洲精品av麻豆狂野| 国产精品综合久久久久久久免费 | 亚洲精品粉嫩美女一区| 精品国产一区二区久久| 老熟妇乱子伦视频在线观看| 亚洲第一青青草原| 热99re8久久精品国产| netflix在线观看网站| 久久久久亚洲av毛片大全| 叶爱在线成人免费视频播放| 欧美日韩亚洲国产一区二区在线观看| 激情视频va一区二区三区| 99国产精品免费福利视频| 丝袜人妻中文字幕| 国产1区2区3区精品| 好男人电影高清在线观看| 日韩人妻精品一区2区三区| 亚洲一区高清亚洲精品| 国产精品免费视频内射| 久久精品影院6| 两性夫妻黄色片| xxxhd国产人妻xxx| 新久久久久国产一级毛片| 久久久精品欧美日韩精品| 黄色片一级片一级黄色片| 亚洲色图 男人天堂 中文字幕| 无遮挡黄片免费观看| 91成人精品电影| 男女下面进入的视频免费午夜 | 精品乱码久久久久久99久播| 91精品三级在线观看| 国产免费现黄频在线看| 欧美激情高清一区二区三区| 欧美日韩乱码在线| 亚洲精品美女久久av网站| 国产精品爽爽va在线观看网站 | 亚洲色图av天堂| 日韩欧美一区视频在线观看| 久久久久国内视频| 久久人人爽av亚洲精品天堂| av天堂在线播放| 中文字幕人妻丝袜制服| 欧美中文综合在线视频| 老司机靠b影院| 成年人免费黄色播放视频| 视频区图区小说| netflix在线观看网站| 搡老岳熟女国产| 久久天躁狠狠躁夜夜2o2o| 50天的宝宝边吃奶边哭怎么回事| 久久久久亚洲av毛片大全| 久久精品国产亚洲av高清一级| 亚洲avbb在线观看| 12—13女人毛片做爰片一| 亚洲av电影在线进入| 首页视频小说图片口味搜索| 国产视频一区二区在线看| 99精品欧美一区二区三区四区| 欧美日韩一级在线毛片| 成年版毛片免费区| 成人亚洲精品一区在线观看| 亚洲国产中文字幕在线视频| www.自偷自拍.com| 电影成人av| xxx96com| 亚洲第一av免费看| 99精品久久久久人妻精品| 一本综合久久免费| 亚洲av美国av| 国产单亲对白刺激| 动漫黄色视频在线观看| 国产亚洲精品久久久久5区| 国产精品久久电影中文字幕| 嫩草影视91久久| 黄频高清免费视频| 国产99久久九九免费精品| 大型黄色视频在线免费观看| videosex国产| svipshipincom国产片| 夜夜看夜夜爽夜夜摸 | 久久久久久亚洲精品国产蜜桃av| 欧美日韩av久久| 亚洲精品美女久久久久99蜜臀| 午夜老司机福利片| 亚洲专区国产一区二区| 免费高清在线观看日韩| 久久影院123| 纯流量卡能插随身wifi吗| 久久中文字幕人妻熟女| 久久午夜综合久久蜜桃| 精品久久久久久,| 久久精品91无色码中文字幕| 香蕉久久夜色| 亚洲男人天堂网一区| 欧美日韩视频精品一区| av片东京热男人的天堂| 搡老乐熟女国产| 国产高清videossex| 曰老女人黄片| 99国产极品粉嫩在线观看| 国产不卡一卡二| av在线播放免费不卡| 国产精品自产拍在线观看55亚洲| 电影成人av| 国产99久久九九免费精品| 欧美+亚洲+日韩+国产| 欧美激情高清一区二区三区| 亚洲精品久久午夜乱码| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品偷伦视频观看了| 美女福利国产在线| 亚洲第一青青草原| 成人三级做爰电影| 在线观看www视频免费| 日韩视频一区二区在线观看| 精品第一国产精品| 涩涩av久久男人的天堂| 中文字幕人妻熟女乱码| 亚洲专区国产一区二区| 性色av乱码一区二区三区2| 欧美大码av| 国产深夜福利视频在线观看| 一区二区三区精品91| 在线国产一区二区在线| 黄片小视频在线播放| 精品无人区乱码1区二区| 亚洲精品一卡2卡三卡4卡5卡| 女性生殖器流出的白浆| 国产精品一区二区免费欧美| 亚洲第一av免费看| 免费在线观看亚洲国产| 超碰97精品在线观看| 精品国产乱码久久久久久男人| 国产片内射在线| 伊人久久大香线蕉亚洲五| 精品福利永久在线观看| 欧美日韩国产mv在线观看视频| 国产蜜桃级精品一区二区三区| 久99久视频精品免费| 97人妻天天添夜夜摸| 久久久久久久午夜电影 | 长腿黑丝高跟| 国产高清国产精品国产三级| 一级a爱片免费观看的视频| 欧美精品一区二区免费开放| 亚洲专区国产一区二区| 欧美成人性av电影在线观看| 色播在线永久视频| 99国产精品免费福利视频| 91麻豆av在线| 国产av在哪里看| 日韩高清综合在线| 亚洲av五月六月丁香网| 人人澡人人妻人| 12—13女人毛片做爰片一| 日日干狠狠操夜夜爽| 亚洲精品久久午夜乱码| 国产精品免费一区二区三区在线| 在线av久久热| 日韩精品青青久久久久久| 在线观看午夜福利视频| 咕卡用的链子| 天堂√8在线中文| 亚洲第一欧美日韩一区二区三区| 嫩草影视91久久| 国产精品 国内视频| 夫妻午夜视频| 看黄色毛片网站| 少妇粗大呻吟视频| 国产精品香港三级国产av潘金莲| 视频区欧美日本亚洲| 99在线视频只有这里精品首页| 美女福利国产在线| 国产又爽黄色视频| 欧美色视频一区免费| 久久精品国产综合久久久| 亚洲全国av大片| 精品免费久久久久久久清纯| 在线观看舔阴道视频| 久久天堂一区二区三区四区| 精品欧美一区二区三区在线| 午夜影院日韩av| 后天国语完整版免费观看| 欧美性长视频在线观看| 黑丝袜美女国产一区| 黑人欧美特级aaaaaa片| 美女大奶头视频| 在线观看日韩欧美| 亚洲成人免费电影在线观看| 国产亚洲精品久久久久5区| 日日夜夜操网爽| 日本 av在线| 很黄的视频免费| 精品国产超薄肉色丝袜足j| 男女之事视频高清在线观看| 美国免费a级毛片| 19禁男女啪啪无遮挡网站| 久久伊人香网站| 涩涩av久久男人的天堂| 国产欧美日韩综合在线一区二区| 国产视频一区二区在线看| 亚洲欧美一区二区三区黑人| 精品久久久精品久久久| 成在线人永久免费视频| 欧美久久黑人一区二区| 欧美老熟妇乱子伦牲交| 国产精品 欧美亚洲| 亚洲精品在线观看二区| 国产精品美女特级片免费视频播放器 | 欧洲精品卡2卡3卡4卡5卡区| 天天影视国产精品| 欧美黄色淫秽网站| 国产av在哪里看| 久久久国产成人免费| 久热爱精品视频在线9| 免费女性裸体啪啪无遮挡网站| 怎么达到女性高潮| 99国产精品免费福利视频| 欧美另类亚洲清纯唯美| 亚洲欧美日韩无卡精品| 夜夜夜夜夜久久久久| 欧美亚洲日本最大视频资源| 99riav亚洲国产免费| www日本在线高清视频| 在线天堂中文资源库| 久久伊人香网站| 国产xxxxx性猛交| 视频区图区小说| 身体一侧抽搐| 超色免费av| 国产精品国产高清国产av| 国产一区二区三区在线臀色熟女 | 一夜夜www| 男男h啪啪无遮挡| 深夜精品福利| 搡老乐熟女国产| 国产主播在线观看一区二区| 熟女少妇亚洲综合色aaa.| 欧美日本亚洲视频在线播放| 很黄的视频免费| 亚洲精品国产区一区二| 两性夫妻黄色片| 国产亚洲av高清不卡| 亚洲av第一区精品v没综合| 久久国产亚洲av麻豆专区| 亚洲av美国av| 在线播放国产精品三级| 男女下面进入的视频免费午夜 | 大型av网站在线播放| 91九色精品人成在线观看| 久久久久久大精品| 怎么达到女性高潮| 长腿黑丝高跟| 亚洲精品av麻豆狂野| 80岁老熟妇乱子伦牲交| 老司机福利观看| 91成人精品电影| 18禁裸乳无遮挡免费网站照片 | 久久亚洲真实| 丁香六月欧美| 亚洲精品美女久久久久99蜜臀| 国产精品一区二区在线不卡| 午夜免费观看网址| 国产片内射在线| 国产蜜桃级精品一区二区三区| 十分钟在线观看高清视频www| 免费在线观看影片大全网站| 国产精品美女特级片免费视频播放器 | 国产精品免费视频内射| av超薄肉色丝袜交足视频| 国产精品乱码一区二三区的特点 | 国产成人一区二区三区免费视频网站| 在线观看舔阴道视频| 亚洲人成伊人成综合网2020| 最新美女视频免费是黄的| 国产亚洲精品久久久久久毛片| 亚洲狠狠婷婷综合久久图片| 中文字幕最新亚洲高清| 91成人精品电影| 亚洲五月天丁香| 热re99久久国产66热| 日本欧美视频一区| 欧美成人性av电影在线观看| 国产片内射在线| 久热这里只有精品99| 成熟少妇高潮喷水视频| 美女大奶头视频| 中出人妻视频一区二区| 国产又色又爽无遮挡免费看| www.精华液| 男女下面进入的视频免费午夜 | 亚洲av成人av| 久久国产亚洲av麻豆专区| 黄色视频不卡| 男男h啪啪无遮挡| 99在线人妻在线中文字幕| 欧美久久黑人一区二区| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 啦啦啦 在线观看视频| 黄色毛片三级朝国网站| 色综合站精品国产| 午夜福利影视在线免费观看| 欧美日韩亚洲综合一区二区三区_| 韩国av一区二区三区四区| 国产有黄有色有爽视频| 色婷婷久久久亚洲欧美| 亚洲一区二区三区欧美精品| 亚洲 欧美 日韩 在线 免费| 久9热在线精品视频| 国产成人欧美| x7x7x7水蜜桃| 精品久久久久久久久久免费视频 | 精品欧美一区二区三区在线| 精品人妻在线不人妻| 国产成人啪精品午夜网站| 国产高清激情床上av| 国产av一区在线观看免费| 一区二区日韩欧美中文字幕| 免费av中文字幕在线| 国产高清激情床上av| 校园春色视频在线观看| 国产高清激情床上av| 国产高清videossex| 欧美午夜高清在线| 淫妇啪啪啪对白视频| www.999成人在线观看| 黑人欧美特级aaaaaa片| 亚洲熟妇熟女久久| 搡老熟女国产l中国老女人| 亚洲九九香蕉| 欧美激情久久久久久爽电影 | e午夜精品久久久久久久| 国产三级在线视频| 精品久久蜜臀av无| 国产又色又爽无遮挡免费看| 欧美另类亚洲清纯唯美| 国产免费av片在线观看野外av| 中文字幕最新亚洲高清| 久久精品影院6| 国产1区2区3区精品| 淫秽高清视频在线观看| 午夜91福利影院| 亚洲国产精品一区二区三区在线| 亚洲av五月六月丁香网| 两性夫妻黄色片| 亚洲一区二区三区色噜噜 | 水蜜桃什么品种好| 在线观看午夜福利视频| 亚洲成人久久性| bbb黄色大片| 亚洲av成人不卡在线观看播放网| 性欧美人与动物交配| 国产成人系列免费观看| 免费久久久久久久精品成人欧美视频| 久久精品国产清高在天天线| 高潮久久久久久久久久久不卡| 国产欧美日韩精品亚洲av| 黄色毛片三级朝国网站| 嫩草影院精品99| 香蕉久久夜色| a级毛片在线看网站| 午夜激情av网站| 这个男人来自地球电影免费观看| 超碰97精品在线观看| 国产成人一区二区三区免费视频网站| 亚洲国产精品sss在线观看 | 午夜福利欧美成人| www.熟女人妻精品国产| 欧美黄色片欧美黄色片| 亚洲av成人一区二区三| 色精品久久人妻99蜜桃| 国产精品99久久99久久久不卡| 亚洲精华国产精华精| 亚洲成人免费电影在线观看| 操出白浆在线播放| 日本欧美视频一区| 国产精品一区二区免费欧美| 国产一区二区三区在线臀色熟女 | 欧美日韩视频精品一区| 亚洲少妇的诱惑av| 老司机福利观看| 久久中文看片网| 亚洲 欧美一区二区三区| 99国产综合亚洲精品| 亚洲国产精品合色在线| 午夜免费激情av| 一进一出抽搐gif免费好疼 | 国产精品国产av在线观看| 久久影院123| 天堂√8在线中文| 精品国产国语对白av| 丝袜美腿诱惑在线| 亚洲一区中文字幕在线| 亚洲色图 男人天堂 中文字幕| 精品电影一区二区在线| 久久草成人影院| 国产主播在线观看一区二区| 精品国内亚洲2022精品成人| 国产精品自产拍在线观看55亚洲| 亚洲色图av天堂| 桃红色精品国产亚洲av| 好男人电影高清在线观看| 国产成人欧美在线观看| 亚洲第一青青草原| 欧美日韩国产mv在线观看视频| 夜夜爽天天搞| 国产欧美日韩一区二区三区在线| www国产在线视频色| 宅男免费午夜| 免费一级毛片在线播放高清视频 | 精品午夜福利视频在线观看一区| 亚洲成人国产一区在线观看| 99国产精品99久久久久| 午夜亚洲福利在线播放| 欧美日韩av久久| 男人的好看免费观看在线视频 | 亚洲中文av在线| 亚洲国产精品一区二区三区在线| 电影成人av| 精品国产亚洲在线| 十八禁人妻一区二区| 黑人巨大精品欧美一区二区mp4| 夜夜爽天天搞| 变态另类成人亚洲欧美熟女 | a级毛片在线看网站| 久久香蕉精品热| 成人黄色视频免费在线看| 超碰97精品在线观看| 久99久视频精品免费| 乱人伦中国视频| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 国产精品秋霞免费鲁丝片| 国产av在哪里看| 精品日产1卡2卡| 99re在线观看精品视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲片人在线观看| 中文字幕最新亚洲高清| 免费高清视频大片| 看片在线看免费视频| 亚洲av电影在线进入| 欧美国产精品va在线观看不卡| 高清黄色对白视频在线免费看| 欧美日韩瑟瑟在线播放| 在线观看午夜福利视频| 热re99久久精品国产66热6| 中文字幕人妻熟女乱码| 日韩欧美国产一区二区入口| 国产精品亚洲av一区麻豆| 久久国产精品男人的天堂亚洲| 99久久综合精品五月天人人| 巨乳人妻的诱惑在线观看| 国产成人欧美在线观看| 美女大奶头视频| 国产精品久久电影中文字幕| 99热只有精品国产| 免费在线观看日本一区| 免费高清在线观看日韩| 亚洲一区二区三区不卡视频| 一夜夜www| 欧美不卡视频在线免费观看 | 国产精品亚洲av一区麻豆| 国产野战对白在线观看| 老司机午夜十八禁免费视频| 久久久久亚洲av毛片大全| 日本撒尿小便嘘嘘汇集6| 日韩欧美国产一区二区入口| 免费不卡黄色视频| 女人高潮潮喷娇喘18禁视频| aaaaa片日本免费| 成人三级黄色视频| 国产精品亚洲av一区麻豆| 国产又爽黄色视频| 午夜久久久在线观看| 男人舔女人下体高潮全视频| av网站在线播放免费| 精品国产美女av久久久久小说| 91老司机精品| 日韩欧美三级三区| svipshipincom国产片| 他把我摸到了高潮在线观看| 精品久久久久久久久久免费视频 | 国产精品美女特级片免费视频播放器 | 亚洲自拍偷在线| 伦理电影免费视频| 91av网站免费观看| 超碰97精品在线观看| 精品日产1卡2卡| 亚洲熟妇熟女久久| 亚洲av五月六月丁香网| 免费在线观看亚洲国产| 国产精品久久久人人做人人爽| 在线观看舔阴道视频| 这个男人来自地球电影免费观看| 无遮挡黄片免费观看| 黄色怎么调成土黄色| av网站免费在线观看视频| 午夜日韩欧美国产| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲欧美98| 亚洲精品一二三| 久久精品亚洲熟妇少妇任你| 成年人黄色毛片网站| av欧美777| 99国产精品99久久久久| 久久久久九九精品影院| 欧美av亚洲av综合av国产av| 免费在线观看日本一区| 久久精品91无色码中文字幕| 免费看a级黄色片| 国产在线精品亚洲第一网站| 18禁美女被吸乳视频| 欧美日韩亚洲高清精品| 午夜成年电影在线免费观看| 免费高清在线观看日韩| 亚洲一区二区三区色噜噜 | 国产亚洲精品第一综合不卡| 中出人妻视频一区二区| 91麻豆精品激情在线观看国产 | 真人做人爱边吃奶动态| 亚洲av成人av| 国产1区2区3区精品| 午夜福利在线免费观看网站| 午夜福利影视在线免费观看| 亚洲欧美一区二区三区黑人| 日韩中文字幕欧美一区二区| 大码成人一级视频| 欧美日本中文国产一区发布| 国产精品98久久久久久宅男小说| 亚洲成a人片在线一区二区| 天堂动漫精品| 激情在线观看视频在线高清| 18禁黄网站禁片午夜丰满|