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

    基于生物多樣性的樟江流域自然保護地空間優(yōu)化

    2022-03-04 12:50:26傅微甄廣潤
    生態(tài)科學(xué) 2022年1期
    關(guān)鍵詞:保護地棲息地代價

    傅微, 甄廣潤

    基于生物多樣性的樟江流域自然保護地空間優(yōu)化

    傅微,*, 甄廣潤

    北京建筑大學(xué)建筑與城市規(guī)劃學(xué)院, 北京 100044

    優(yōu)化生態(tài)保護空間格局是我國生態(tài)文明建設(shè)的重要內(nèi)容, 也是國際保護生態(tài)學(xué)研究的前沿方向。評估資源稟賦、識別并確定優(yōu)先次序和預(yù)測保護區(qū)對保護和管理瀕危物種至關(guān)重要。樟江流域內(nèi)已占11.9%流域面積的保護區(qū)仍出現(xiàn)77.9%流域面積的生境顯著退化, 發(fā)展與保護亟待合理平衡。Zonation4GUI軟件以互補性非窮盡式算法為基礎(chǔ), 在優(yōu)化已有保護網(wǎng)絡(luò)并表征稀有性和獨特性方面結(jié)果較為理想, 以此展開樟江流域自然保護地空間整合優(yōu)化路徑研究。先人工預(yù)處理物種棲息地分布格局、物種保護權(quán)重、規(guī)劃單元, 并迭代選擇輸入性生境退化條件圖層、保護代價圖層, 結(jié)合人工干擾情況提出三種不同情景下的優(yōu)化方案。結(jié)果表明: 保護優(yōu)先區(qū)對代理物種棲息地有良好的覆蓋效果, 生態(tài)完整性優(yōu)先情景和最小社會經(jīng)濟和土地資源代價的保護優(yōu)化情景均呈現(xiàn)擴建小于8%的一級優(yōu)先區(qū)域, 所有保護對象均受到保護; 加入27%左右的所有優(yōu)先區(qū)域, 兩種優(yōu)化方案覆蓋所有保護對象棲息地的比例提升70.0%和59.0%。雖然生態(tài)完整性優(yōu)先的優(yōu)化情景優(yōu)于最小社會經(jīng)濟和土地資源代價下的優(yōu)化情景, 但保護代價更高。人為干擾強烈的搶救優(yōu)先的優(yōu)化情景中, 90.7%的保護對象棲息地覆蓋范圍有明顯提升。從技術(shù)和操作層面, 為重構(gòu)中國自然保護地體系整合優(yōu)化方案提供路徑和方法。

    系統(tǒng)保護規(guī)劃; 保護空缺; 喀斯特地區(qū); 流域整體保護; 保護目標

    0 前言

    生物多樣性為人類的福祉和長期生存提供生態(tài)系統(tǒng)服務(wù)和商品, 保持生態(tài)系統(tǒng)穩(wěn)定性, 維持生態(tài)平衡[1-2], 是全球生態(tài)系統(tǒng)的不可分割的要素。生物多樣性在農(nóng)業(yè)、林業(yè)、漁業(yè)和旅游業(yè)等推動經(jīng)濟發(fā)展中發(fā)揮著重要作用。人口和人均消費的持續(xù)增長導(dǎo)致了不可持續(xù)的開發(fā)利用生物多樣性, 物種的過度開發(fā)、外來入侵物種的擴張、氣候變化、森林退化和獨特棲息地的破壞等驅(qū)動力嚴重影響生物多樣性[3]。隨著生物多樣性喪失危機的不斷加劇, 保護生物多樣性的意識和政策機制正在區(qū)域、國家和國際層面增強[4], 旨在最小化物種退化和生境破壞。為了增加保護工作的可行性, 并考慮到世界上大多數(shù)地區(qū)的生物多樣性清單都不完整, 通常使用物種多樣性代理確定支持多個物種和棲息地的保護區(qū)系統(tǒng)保護規(guī)劃[5-6]。在過去的三十年里, 許多研究采用多樣性格局代理(如物種豐富度)作為最小覆蓋集合表征研究區(qū)內(nèi)幾乎所有的物種[7], 選址潛在的自然保護區(qū)。根據(jù)物種豐富度選擇保護地是最傳統(tǒng)的方法之一(案例可見[8-9])。然而, 當(dāng)保護規(guī)劃的目標旨在優(yōu)化已有保護網(wǎng)絡(luò)并表征關(guān)鍵區(qū)域的最大物種數(shù)量時, 基于物種豐富度并非保護優(yōu)先級篩選的良好指標, 因其無法反映互補性(Complementarity), 即一種基于場地的保護價值衡量方法, 確保被選擇納入保護區(qū)網(wǎng)絡(luò)的地點補充已有保護空間; 同時無法說明特有性和稀缺性, 從而出現(xiàn)受保護物種的數(shù)量減少現(xiàn)象[10]?;パa算法可解釋物種稀有性, 選擇的地點是相互補充的, 每個地點不一定是物種最豐富的, 但每個地點在物種組成上都不同, 可產(chǎn)生一種包含所有物種的地點的解決方案, 在有限的范圍內(nèi)給予所有物種更高的優(yōu)先權(quán)[11]。

    樟江流域呈扇形展布, 地跨貴州省黔南布依族苗族自治州荔波、三都縣, 屬于珠江流域。主河道樟江發(fā)源于貴州省荔波縣佳榮鄉(xiāng)月亮山, 由東北向西南流, 流域面積1673.9 km2。樟江流域是中國生物多樣性熱點地區(qū)之一, 景觀類型發(fā)育具有典型性、稀有性、脆弱性以及多樣性, 表現(xiàn)在世界范圍內(nèi)喀斯特地貌上保存完好的區(qū)域, 是長江、珠江上游重要的生態(tài)屏障。其豐富的、獨特的生物多樣性為流域上游、下游提供生態(tài)系統(tǒng)服務(wù)和商品, 但由于其脆弱性和資源密集性, 易受到包括過度開發(fā)在內(nèi)的各種驅(qū)動力變化因素的影響。雖然已建立了5個自然保護地, 但這些保護地提供的保護并不總是有效的[12]。土地開發(fā)和環(huán)境變化導(dǎo)致棲息地生境減少, 可能影響物種的分布, 經(jīng)濟發(fā)展和生物多樣性保護的競爭需求也存在著激烈矛盾沖突[13]。

    在此背景下, 以樟江流域為研究區(qū)域旨在關(guān)注并回答以下問題: 近20年已有自然保護地的保護效益是否足夠?人類開發(fā)利用的干擾程度及與之帶來的生境退化趨勢如何?采用互補算法最大化保證瀕危物種數(shù)量和最優(yōu)化補充保護空間; 篩選考慮不同保護策略的保護空缺格局和保護效果。

    1 整合優(yōu)化路徑與方法

    1.1 自然稟賦研究

    1.1.1 對象選擇

    依據(jù)外業(yè)調(diào)查、專家知識, 選取研究區(qū)內(nèi)重點保護對象54種, 其中植物38種, 動物16種, 如表1所示。植物物種相關(guān)信息及分布來自中國科學(xué)院植物研究所中國珍稀瀕危植物信息系統(tǒng)(http:// www.iplant.cn/rep/), 動物物種相關(guān)信息及分布來源自中國動物主題庫(http://www.zoology.csdb.cn/), 并于2017年9月至2018年1月和2018年12月開展外業(yè)工作予以核實。

    表1 重點保護物種列表

    1.1.2 空間分布

    空間分布采用統(tǒng)計物種分布模型(Statistical Species Distribution Models), 該模型結(jié)合了物種名錄、自然生態(tài)系統(tǒng)與資源環(huán)境數(shù)據(jù), 提供棲息地適宜性的連續(xù)圖層預(yù)測[14]。首先查找物種分布范圍, 再依據(jù)棲息地類型對分布范圍進行修正: 植物物種的棲息地類型通過海拔高程、土壤類型、植物類型制圖疊加生成; 動物物種棲息地類型則通過植被類型、海拔高程及土地利用類型制圖疊加獲取, 并在ArcGIS支持下疊加生成物種適宜空間分布圖。此外, 提取每種保護對象的適宜棲息地面積, 作為保護目標分析指標。植被類型分布來源于《1:1 000 000 中國植被圖集》矢量化數(shù)據(jù); 土壤類型空間分布數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http:// www.resdc.cn); 數(shù)字高程模型SRTM_DEM數(shù)據(jù)(30 m分辨率), 來源于http://srtm.csi.cgiar.org/。

    1.1.3 保護資源權(quán)重

    依據(jù)自然保護地中典型性、稀有性、自然性、多樣性、科學(xué)價值、美學(xué)價值(觀賞價值)、社會和經(jīng)濟價值等指標, 整合保護對象權(quán)重。如公式(1)所示, 選用國家保護等級、地方特有性、物種瀕危程度、經(jīng)濟價值和潛在棲息地面積五個指數(shù)制定權(quán)重方案。

    =(endangered+protection+endemic+economic+area)/5 (1)

    式(1)中是物種的保護權(quán)重;endangered按物種的瀕危程度打分;protection按國家保護等級,endemic特有性打分;economic按照植物的經(jīng)濟用途打分, 很多植物具有藥用、色素、野菜、觀賞、油脂、纖維、蜜源、香料、野果等經(jīng)濟價值。area按物種潛在棲息地面積打分, 最高分5分。動物不考慮經(jīng)濟價值, 去掉economic后四項取平均值[15]。

    1.1.4 規(guī)劃單元

    為突出流域生態(tài)系統(tǒng)結(jié)構(gòu)的完整性以及流域上下游之間的連續(xù)性, 規(guī)劃單元采用ArcGIS Hydrology工具構(gòu)建1421個集水區(qū), 因為集水區(qū)單元具有更大的自然相似性, 符合物種的自然分布特征, 可維持保護對象生境和地貌單元的完整性。

    1.2 保護現(xiàn)狀梳理

    流域內(nèi)已有樟江國家級重點風(fēng)景名勝區(qū)、國家級茂蘭喀斯特森林自然保護區(qū)和3個縣級自然保護區(qū)(含核心區(qū)、緩沖區(qū)及實驗區(qū)), 分布數(shù)據(jù)來源于黔南州林業(yè)局根據(jù)實地調(diào)查記錄經(jīng)緯度信息在ArcGIS10.2中空間成圖, 作為去除掩膜圖層, 為完善保護優(yōu)化格局提供基礎(chǔ)。

    1.2.1 流域自然保護現(xiàn)狀條件

    流域自然保護現(xiàn)狀條件采用植被凈初級生產(chǎn)力2000年–2015年年平均值趨勢變化作為指標表征生境特征退化情況。植被凈初級生產(chǎn)力為生態(tài)系統(tǒng)在特定時期內(nèi)的生物總量, 表征生態(tài)質(zhì)量, 是評價生態(tài)系統(tǒng)功能特征的重要參量。數(shù)據(jù)基于CASA (Carnegie Ames Stanford Approach)模型計算, 采用逐像元一元線性回歸分析方法, 闡釋2000—2015年研究區(qū)植被凈初級生產(chǎn)力的時間序列變化趨勢, 基于斜率和Pearson相關(guān)系數(shù)判定退化程度。此外, 結(jié)合《樟江流域綜合規(guī)劃》資料, 評估河流物理特征、水文水資源特征、水質(zhì)狀況、水生生物、社會服務(wù)功能情況, 繪制流域自然保護現(xiàn)狀景觀條件圖。生境退化條件情況中條件值可以在0和1之間變化, 值為1表示原始狀態(tài), 生境適宜性或特征發(fā)生沒有退化; 0的值表示完全退化的狀態(tài), 任何負值(丟失的數(shù)據(jù))都將被視為零。

    1.2.2 人為干擾度量與保護代價成本

    基于生物多樣性的保護優(yōu)化格局完整性的實現(xiàn)受制于人類活動空間的干擾, 為維護完整性而改變?nèi)藶楦蓴_空間所需保護代價成本隨人為影響程度增加而增加。因此參照相關(guān)研究[16], 選取公路、鐵路、城鎮(zhèn)、農(nóng)村居民點、水壩作為計算規(guī)劃單元保護代價的因子, 將每個因子的度量標準化后, 賦予不同的權(quán)重(表2)。相加計算得到每個規(guī)劃單元的保護代價指數(shù)。將保護代價指數(shù)用保護區(qū)覆蓋面積比例進行調(diào)整, 得到最終保護代價成本分布圖(cost layer)。道路、居民點數(shù)據(jù)、土地利用數(shù)據(jù)均來源于黔南州國土部門2017年土地利用數(shù)據(jù); 大壩分布來源于《樟江流域綜合規(guī)劃》, 計算如公式(2-3)所示。

    表2 各影響因子的度量和權(quán)重

    C′=C(1–0.5) (3)

    式(2)中:V為每個保護單元內(nèi)因子的度量值;W為因子的權(quán)重;為每個保護單元內(nèi)因子的個數(shù); C為每個保護單元調(diào)整前的代價值; 式(3)中C′為每個保護單元經(jīng)過保護區(qū)調(diào)整后的代價值;為每個保護單元內(nèi)被保護區(qū)覆蓋的面積比。

    1.3 結(jié)構(gòu)優(yōu)化分析

    1.3.1 保護價值評估

    保護價值評估基于空間保護優(yōu)先級劃分和大尺度空間保護規(guī)劃的前沿軟件Zonation4GUI, 基于互補的算法對單元進行分層排序[11]。首先假設(shè)所有單元都可能是重要的, 并在每次迭代后逐步刪除最不重要的單元。通過最小化保護價值的邊際損失, 最大程度降低生物多樣性損失, 并考慮到保護對象及分布的典型性、稀有性、多樣性、高價值等特征的連通性需求和優(yōu)先級。最終保證整體連通性為焦點找出最能配合現(xiàn)有保護地的范圍, 即最有保護價值的地區(qū)[17]。

    1.3.2 分析保護空缺

    本研究區(qū)域原有自然保護地(Protected Area, 簡稱PA)占總面積的11.9%, 因此, 選取樟江流域景觀總面積的30%作為目標優(yōu)先區(qū), 將最高保護價值區(qū)域的10%劃分為一級優(yōu)先區(qū), 10%—20%劃分為二級優(yōu)先區(qū), 20%—30%劃分為三級優(yōu)先區(qū), 最終產(chǎn)生不同級別的保護優(yōu)先格局。

    1.4 優(yōu)化整合方案

    分別將前述優(yōu)化路徑中保護對象空間分布圖層、物種名錄、保護權(quán)重、保護規(guī)劃單元圖層、已有自然保護地圖層、自然保護現(xiàn)狀條件圖層、流域范圍圖層代入Zonation4GUI軟件, 選擇性代入保護代價圖層作為不同保護策略情景展開分析。采用了保持格局連通性的核心區(qū)移除規(guī)則(Core Area Zonation, 簡稱CAZ) 和邊緣移除規(guī)則(Rule of Edge Removal)[18], 分層優(yōu)先逐步得到一個高連通性景觀結(jié)構(gòu)的嵌套序列。進一步研究保護空缺區(qū)域與周邊自然保護地的關(guān)系, 將生態(tài)系統(tǒng)類型相同、主要保護對象和保護目標一致的保護空缺區(qū)域補充并入周邊自然保護地, 設(shè)立新的自然保護地。

    分層優(yōu)化過程中, 引入三種不同保護策略情景下的保護空間優(yōu)化方案, (1)生態(tài)完整性優(yōu)先的優(yōu)化方案, 該方案不考慮保護代價迭代, 將高保護價值規(guī)劃單元建立優(yōu)先保護格局, 最大程度降低生物多樣性損失; (2)最小的社會經(jīng)濟和土地資源代價下優(yōu)化方案, 該方案考慮保護代價, 并作為成本層參與迭代分析, 平衡利益相關(guān)者和建立較為可行的保護區(qū); (3)搶救性保護的優(yōu)化方案, 選取人為干擾強烈的分布區(qū)疊加生態(tài)完整性高保護價值區(qū)域作為優(yōu)先保護格局。

    2 結(jié)果與分析

    2.1 流域保護現(xiàn)狀堪憂

    當(dāng)現(xiàn)狀自然保護區(qū)、風(fēng)景名勝區(qū)和世界遺產(chǎn)地等已有自然保護地總面積為199.1 km2的情形下, 樟江流域生態(tài)退化現(xiàn)狀如圖1所示, 生境顯著退化面積置信度超過99%的退化面積占流域總面積的46.8%, 置信度超過95%的顯著退化面積占流域總面積的67.0%, 置信度超過90%的顯著退化面積占總面積的77.9%。可見, 樟江流域近16年時間內(nèi)生態(tài)退化明顯且嚴峻, 水源涵養(yǎng)能力下降, 流域的生態(tài)環(huán)境質(zhì)量呈逐年下降趨勢, 生態(tài)脆弱性加劇。根據(jù)樟江荔波水文站的資料, 生態(tài)流量干流部分退化25%, 從河流連通性、河道蜿蜒度、河道渠化程度、河岸穩(wěn)定性、河床穩(wěn)定性、河流形態(tài)指標表明河流物理結(jié)構(gòu)干流退化了31.2%, 上游水源區(qū)退化了7.3%, 導(dǎo)致水體污染, 水量減少等現(xiàn)象。

    根據(jù)適宜棲息地分布格局, 動物棲息地分布面積居高, 其中穿山甲、普通鵟、小鴉鵑、褐翅鴉鵑、小靈貓和金貓分布較廣, 輻花苣苔、中國蕨、白花兜蘭、單座苣苔、狹葉含笑等分布少, 輻花苣苔分布適宜范圍為899 m2。已有保護區(qū)內(nèi)不同類別物種棲息地分布情況發(fā)現(xiàn), 半楓荷、單性木蘭、輻花苣苔、福建柏、紅豆杉、黃杉、楠木、異性玉葉金花、中國蕨尚未覆蓋, 其他79.6%的覆蓋物種棲息地面積比例均較低, 平均比例為7.1%, 表明現(xiàn)有保護地空間保護情況有較大不足。

    保護代價成本分析結(jié)果顯示(圖2), 樟江流域改變需要的保護代價高的規(guī)劃單元較集中分布在樟江流域中線地帶——荔波縣城樟江干流、上游居民點集中區(qū)為軸線, 以及分布于三荔水庫大壩、的馬水庫大壩集水區(qū)區(qū)域。

    2.2 生態(tài)完整性優(yōu)先的優(yōu)化方案

    如圖3所示, 高保護價值且存在保護空缺的區(qū)域集中分布在樟江主河道兩岸亞熱帶常綠闊葉林帶、樟江源頭區(qū)域拉易河、爽姑河匯入水昔河的交匯處、水令河上游、爽姑河上游、拉灘河上游, 茂蘭鎮(zhèn)瑤麓喀斯特典型生態(tài)系統(tǒng)亞熱帶常綠闊葉林區(qū)域。這些地區(qū)物種類型分布豐富, 具有較好的優(yōu)先級。一、二、三級優(yōu)先區(qū)共補充保護空缺面積453.6 km2(表3), 統(tǒng)計每個物種的棲息地在保護優(yōu)先區(qū)內(nèi)分布情況, 共計45個物種的棲息地面積比例超過60%, 約占全部物種的83.3%, 實現(xiàn)對森林生態(tài)環(huán)境、生物多樣性、珍稀野生動植物和水源涵養(yǎng)林地的保護。其中, 一級優(yōu)先區(qū)范圍內(nèi)新增面積比例為7.4%, 即可實現(xiàn)涵蓋100%研究選取保護對象的部分或者全部適宜棲息地, 其中金貓、獼猴、普通鵟等鳥類和哺乳類動物分布空間比例最高, 對如輻花苣苔等此類適宜棲息地面積極小的物種也實現(xiàn)其棲息地的完全保護。二級優(yōu)先區(qū)范圍內(nèi)新增面積比9.7%, 進一步實現(xiàn)了單座苣苔、紅豆杉、楠木、四藥門花、狹葉含笑、中國蕨保護對象棲息地生境的完全保護。三級優(yōu)先區(qū)范圍內(nèi)不足總面積的0.0033%為原有保護區(qū)范圍, 繼續(xù)實現(xiàn)了福建柏、藍翅八色鶇、單性木花、翠柏、白花兜蘭棲息地生境的完全保護。

    圖 1 生態(tài)退化顯著性分布

    Figure 1 Distribution of condition layer

    圖 2 保護代價成本分布

    Figure 2 Distribution of cost layer

    2.3 最小社會經(jīng)濟和土地資源代價的優(yōu)化方案

    最小的社會經(jīng)濟和土地資源代價下最大程度開展保護方案, 高保護價值區(qū)域集中分布在樟江主河道兩岸亞熱帶常綠闊葉林地帶、水便河中下游、樟江源頭區(qū)域水昔河上游拉灘河集水區(qū)、爽姑河匯入水昔河的交匯處以及水令河上游(圖4)。通過互補算法Zonation軟件分析所得優(yōu)化空間, 雖然相對于動物的適宜性棲息地分布, 部分植被保護對象表現(xiàn)出可適應(yīng)的棲息地面積非常有限, 但在納入一級優(yōu)先區(qū)后, 已覆蓋全部38種植物和16種動物, 新增優(yōu)先區(qū)面積占流域總面積的8.0%。植物物種棲息地所占比例最少的4個對象為輻花苣苔、狹葉含笑、紅豆杉、單座苣苔、中國蕨, 面積分別為75 m2、257 m2、275 m2、276 m2、1235 m2。補充納入二級優(yōu)先區(qū)內(nèi)也同樣全覆蓋保護對象, 保護空缺占流域總面積9.5%, 輻花苣苔適宜棲息地得以完全保護。補充納入三級優(yōu)先區(qū)作為輔助補充提高物種可棲息地的面積比例, 保護空缺占比9.7%, 生物多樣性代理保護對象中, 單座苣苔、輻花苣苔、福建柏、紅豆杉、中國蕨的棲息地完全覆蓋。一、二、三級優(yōu)先區(qū)共新增優(yōu)化空間454.7 km2, 共計32個物種的棲息地面積比例超過50%、22個物種的棲息地面積超過60%, 該保護方案可進一步實現(xiàn)環(huán)境保護與經(jīng)濟發(fā)展同步開展。

    圖 3 生態(tài)完整性優(yōu)先保護方案

    Figure 3 Landscape rankings for scenarios1

    圖 4 最小的社會經(jīng)濟和土地資源代價下保護方案

    Figure 4 Landscape rankings for scenarios2

    2.4 搶救性保護的優(yōu)化方案

    搶救優(yōu)先的保護方案中(圖5和表3), 一、二、三級優(yōu)先區(qū)共新增優(yōu)化空間126.6 km2, 約占流域總面積的7.6%, 保護的物種棲息地除翠柏、紅豆杉、單座苣苔、楠木5類植物外, 均有涉及??臻g優(yōu)先區(qū)主要分布于沿樟江主河道兩岸地帶, 這一水陸交接地帶從棲息地分析結(jié)果發(fā)現(xiàn)存在超過半數(shù)代表性動植物物種和特有種, 表明了這一地帶存在搶救性保護空缺, 研究結(jié)果也與政府報告中提到主河道兩岸緩沖區(qū)保護效率不高相吻合。此外, 一級優(yōu)先區(qū)中, 保護空缺面積為29.2 km2, 主要分布于水賈河上游集水區(qū); 二級優(yōu)先區(qū)保護空缺面積39.3 km2, 主要分布于水各河和水令河上游集水區(qū); 三級優(yōu)先區(qū)保護空缺面積58.0 km2, 集中于水昔河和水便河上游山谷集水區(qū)地帶, 也是人為干擾明顯區(qū)域。

    圖 5 搶救性優(yōu)先保護方案

    Figure 5 Landscape rankings for scenarios3

    表3 三種預(yù)景方案下分級保護優(yōu)先區(qū)的分布比例與面積(單位: km2)

    3 討論

    研究揭示了至少兩個值得探討的方面: 首先, 用于保護規(guī)劃的資源通常是有限的, 這就要求采取盡可能高效的保護行動, 以最大限度地保護生物多樣性[7]。在確定物種的空間保護優(yōu)先權(quán)及地點保護價值時, 相對于物種豐富度的窮盡式算法, 互補性算法是非窮盡式的, 即并非極致追求物種最豐富的地點, 而是綜合物種組成, 集中代表所有的物種[11]。對比近年來采用物種豐富度算法的研究發(fā)現(xiàn), 探究已有保護區(qū)與算法模擬保護區(qū)展開額外重疊分析[19], 但從實際建設(shè)可行性層面, 互補性算法可在已有保護區(qū)的范圍內(nèi)優(yōu)化新增, 而非調(diào)整原有保護區(qū)微分, 顯著提高了景觀空間優(yōu)先化的有效性。互補算法Zonation的核心區(qū)域分區(qū)(CAZ), 也對物種稀有性和獨特性給予了較好的解析, 因為其優(yōu)先考慮對物種重要的核心區(qū)域, 測試了所有可能的地點選擇場景, 在有限的范圍內(nèi)給予那些物種更高的優(yōu)先權(quán)。近三十年來, 保護生物學(xué)家已經(jīng)認識到互補性比物種豐富度更有效; 然而, 它尚未成為生物多樣性保護研究的主流方法[20-21], Astudillo?Scalia等[22]學(xué)者猜測是因為這個概念還沒有被很好地理解, 特別是基于典型區(qū)域應(yīng)用實踐層面。

    其次, 通過比較保護代價層介入與否的系統(tǒng)保護規(guī)劃的有效性, 選擇作為保護優(yōu)先級的空間分布有明顯差別。取決于人類干擾對生態(tài)完整性造成的影響, 未納入保護代價的保護目標的實現(xiàn)情況在所需保護面積小于納入后。因此, 人為干擾加劇條件下, 未來所需擴大保護地的范圍也有待相應(yīng)提高。人為干擾強烈的搶救性優(yōu)先保護區(qū)初步建議開展生境恢復(fù)試點, 改善其生態(tài)功能。緩解自然保護地與社區(qū)之間在自然資源利用上的矛盾, 社區(qū)共管[23]與生態(tài)旅游[24]在這些地帶的重要體現(xiàn), 補充收入來源, 開展環(huán)境容量評估。尋找保護與周邊用地、產(chǎn)業(yè)發(fā)展的平衡點, 較大程度改善矛盾沖突的核心地帶。確定優(yōu)先保護區(qū)域是生物多樣性保護的重要第一步, 但關(guān)鍵是根據(jù)這些建議采取行動[25], 樟江流域在自然保護地空間范疇與管理范疇需協(xié)調(diào)一致的努力。

    4 結(jié)論

    近16年樟江流域生境退化置信度超過95%的面積占流域總面積的67.0%, 流域內(nèi)生物多樣性資源交易頻繁, 許多物種棲息地分布逐漸破碎, 特別是水系干流和上游水源涵養(yǎng)地帶仍未得到保護。已有保護區(qū)覆蓋代理保護對象棲息地面積比例不足8%, 基于生物多樣性的樟江流域自然保護地空間優(yōu)化迫在眉睫, 是保護地體系建設(shè)、管理和決策的重要依據(jù)。采用Zonation核心區(qū)域分區(qū)(CAZ)優(yōu)化路徑, 對棲息地分布小、物種豐富度不足的物種保護互補性較好?;谌N不同情景下的優(yōu)化結(jié)果表明, 保護優(yōu)先區(qū)對代理物種棲息地有良好的覆蓋效果, 生態(tài)完整性優(yōu)先情景和最小社會經(jīng)濟和土地資源代價的保護優(yōu)化情景均呈現(xiàn): 在現(xiàn)有格局中加入小于8%的1級優(yōu)先區(qū)域, 所有代理物種均受到保護。加入27%左右的所有優(yōu)先區(qū)域, 兩種優(yōu)化方案覆蓋所有保護對象棲息地的比例提升70.0%和59.0%。雖然生態(tài)完整性優(yōu)先的優(yōu)化情景優(yōu)于最小社會經(jīng)濟和土地資源代價下的優(yōu)化情景, 但保護代價更高。人為干擾強烈的搶救優(yōu)先的優(yōu)化情景中, 90.7%的保護對象棲息地覆蓋范圍有明顯提升。三個優(yōu)化情景中優(yōu)先區(qū)的分布均可確定樟江主河道兩岸、爽姑河上游、水令河上游等重點保護空缺地區(qū), 亟待新建自然保護地, 強化保護管理措施。樟江流域生物多樣性保護空間格局優(yōu)化研究對互補算法保護規(guī)劃具有典型意義, 對我國乃至世界上這一獨特的喀斯特資源及其所提供的生態(tài)系統(tǒng)服務(wù)具有重要保護作用和科學(xué)支撐。

    [1] RANDS M R, ADAMS W M, BENNUN L, et al. Biodiversity Conservation: Challenges Beyond 2010[J]. Science, 2010, 329(5997): 1298–1303.

    [2] BRUGIERE D, SCHOLTE P. Biodiversity gap analysis of the protected area system in poorly–documented Chad[J]. Journal for Nature Conservation, 2013, 21(5): 286–293.

    [3] BUTCHART S H, WALPOLE M, COLLEN B, et al. Global Biodiversity: Indicators of Recent Declines[J]. Science, 2010, 328(5982): 1164–1168.

    [4] KULLBERG P, MOILANEN A. How do recent spatial biodiversity analyses support the convention on biological diversity in the expansion of the global conservation area network[J]. Natureza & Conservacao, 2014, 12(1): 3–10.

    [5] RODRIGUES A S, BROOKS T M. Shortcuts for Biodiversity Conservation Planning: The Effectiveness of Surrogates[J]. Annual Review of Ecology, Evolution, and Systematics, 2007, 38(1): 713–737.

    [6] DE ALBUQUERQUE F S, BEIER P. Using abiotic variables to predict importance of sites for species representation[J]. Conservation Biology, 2015, 29(5): 1390–1400.

    [7] Moilanen A, Wilson K A, Possingham H p. Spatial conservation prioritization: quantitative methods & computational tools[M]. Oxford University Press, New York, 2009.

    [8] ROBERTS C M, MCCLEAN C J, VERON J, et al. Marine biodiversity hotspots and conservation priorities for tropical reefs[J]. Science, 2002, 295(5558): 1280–1284.

    [9] RAMIREZ F, AFAN I, DAVIS L S, et al. Climate impacts on global hot spots of marine biodiversity[J]. Science Advances, 2017, 3(2).

    [10] VEACH V, MININ E D, POUZOLS F M, et al. Species richness as criterion for global conservation area placement leads to large losses in coverage of biodiversity[J]. Diversity and Distributions, 2017, 23(7): 715–726.

    [11] MOILANEN A, POUZOLS F, MELLER L, et al. Zonation: spatial conservation planning and software v. 4 [EB/OL]. http: //cbig. it. helsinki. fi, 2014.

    [12] WATSON J E, DUDLEY N, SEGAN D B, et al. 2014. The performance and potential of protected areas. Nature [J], 515: 67–73.

    [13] MA Z, CHEN Y, MELVILLE D S, et al. Changes in area and number of nature reserves in China[J]. Conservation Biology, 2019, 33(5): 1066–1075.

    [14] ELITH J, FRANKLIN J 2013. Species distribution modeling [M], Encyclopedia of Biodiversity: Second Edition. Elsevier Inc. : 692–705.

    [15] 馬琳, 李俊清. 基于系統(tǒng)保護規(guī)劃的長白山闊葉紅松林保護網(wǎng)絡(luò)優(yōu)化研究[J]. 生態(tài)學(xué)報, 2019, 39(22): 8547– 8555.

    [16] 梁晨, 李曉文, 崔保山, 等. 中國濱海濕地優(yōu)先保護格局構(gòu)建[J]. 濕地科學(xué), 2015, 13(6): 660–666.

    [17] LEHTOM?KI J, KUOKKANEN P, HANSKI I, et al. Planning of forest conservation areas using high–resolution GIS data and software for spatial conservation prioritization [J]. Forest Ecology and Management, 2009, 258: 2439–2449.

    [18] MOILANEN A, ANDERSON B J, EIGENBROD F, et al. Balancing alternative land uses in conservation prioritization [J]. Ecological Applications, 2011, 21: 1419–1426.

    [19] 史雪威, 張路, 張晶晶, 等. 西南地區(qū)生物多樣性保護優(yōu)先格局評估[J]. 生態(tài)學(xué)雜志, 2018, 37(12): 3721–3728.

    [20] ALBUQUERQUE F, BEIER P. Rarity–Weighted Richness: A Simple and Reliable Alternative to Integer Programming and Heuristic Algorithms for Minimum Set and Maximum Coverage Problems in Conservation Planning[J]. PloS one, 2015, 10: 1–7.

    [21] ALBUQUERQUE F, BEIER P. Global patterns and environmental correlates of high–priority conservation areas for vertebrates[J]. Journal of Biogeography, 2015, 42: 1397–1405.

    [22] ASTUDILLO–SCALIA Y, ALBUQUERQUE F. Why should we reconsider using species richness in spatial conservation prioritization[J]. Biodiversity and Conserva-tion, 2020: 1–13.

    [23] 徐建英, 陳利頂, 呂一河, 等. 保護區(qū)與社區(qū)關(guān)系協(xié)調(diào): 方法和實踐經(jīng)驗[J]. 生態(tài)學(xué)雜志, 2005, 24(1): 102–107.

    [24] 劉洋, 呂一河, 陳利頂,等. 自然保護區(qū)生態(tài)旅游影響評價: 進展與啟示[J]. 自然資源學(xué)報, 2005,20(5): 771–779.

    [25] GAME E T, KAREIVA P, POSSINGHAM H P, et al. Six common mistakes in conservation priority setting[J]. Conservation Biology, 2013, 27(3): 480–485.

    Biodiversity-based spatial optimization for natural conservation areas of Zhangjiang River basin

    FU Wei*, ZHEN Guangrun

    School of Architecture and Urban Planning, Beijing University of Civil Engineering and Architecture, Beijing 100044, China

    It assumes an essential role to optimize the spatial pattern of ecological conservation within the scope of ecological civilization construction in China, which is also the frontier orientation of international research regarding ecological conservation. Assessing resource endowments, identifying and prioritizing and predicting conservation areas are critical to the protection and management of endangered species. Although conservation areas account for 11.9% of the Zhangjiang River watershed, significant habitat degradation has occurred across 77.9% of the watershed, highlighting an urgent need to balance development and protection. Based on a complementary non-exhaustive algorithm, the Zonation4GU software can produce ideal outcomes in optimizing existing conservative networks and characterizing rarity and uniqueness. This paper presents a study of the route towards spatial integration and optimization in natural conservations in the Zhangjiang River watershed based on the software. First, the distributional patterns of habitats of species, weights of species protection and planning units were manually preprocessed; the layer of conditions for imported habitat degradation and the layer of protection cost were iteratively selected; finally, optimization solutions to three different scenarios were proposed in combination of the artificial interference. The results show that the area prioritizing protection has an ideal coverage on habitats of proxy species, optimization solutions to scenario prioritizing ecological integrity and that prioritizing the minimum social, economic and land-resource costs both present as tire-1 priority areas with an construction expansion rate smaller than 8% and all protected targets receiving protection; after around 27% of the all priority areas were incorporated, the two optimization solutions exhibited a 70.0% and 59.0% rise, respectively, in covering the habitats of all protected targets. Although the optimization scenario prioritizing ecological integrity produces better performances than that prioritizing the minimum social, economic and land-resource costs, it entails a higher cost of protection. With respect to the optimization scenario prioritizing intense artificial interference such as urgent conservation, 90.7% of the inhabits for the protected targets have seen a significant coverage increase. From the technical and operative perspective, this study provides a viable route towards creating an integration and optimization solution to China’s natural conservation system.

    systematic conservation planning; conservation gap; Karst region; integrated watershed protection; conservation targets

    10.14108/j.cnki.1008-8873.2022.01.006

    傅微, 甄廣潤. 基于生物多樣性的樟江流域自然保護地空間優(yōu)化[J]. 生態(tài)科學(xué), 2022, 41(1): 50–58.

    FU Wei, ZHEN Guangrun, SUN Zhe, et al. Biodiversity-based spatial optimization for natural conservation areas of Zhangjiang River basin[J]. Ecological Science, 2022, 41(1): 50–58.

    K903

    A

    1008-8873(2022)01-050-09

    2020-05-21;

    2020-07-10基金項目:國家自然科學(xué)基金項目(41901220) ;北京市優(yōu)秀人才培養(yǎng)資助青年骨干個人項目(24269220006)

    傅微(1988—), 女, 湖南長沙人, 博士, 講師, 主要從事景觀生態(tài)規(guī)劃研究, E-mail: fuwei0807@163.com

    傅微

    猜你喜歡
    保護地棲息地代價
    四川大熊貓棲息地
    全國自然保護地大檢查6月啟動
    綠色中國(2018年11期)2018-11-06 07:16:30
    BEAN SCENES
    抵達棲息地
    廈門航空(2018年4期)2018-04-25 10:49:27
    愛的代價
    海峽姐妹(2017年12期)2018-01-31 02:12:22
    保護地蔬菜白粉虱的發(fā)生與防治
    代價
    赤峰市保護地黃瓜主要病害的無公害治理
    保護地栽培及存在問題
    成熟的代價
    亚洲av电影不卡..在线观看| 日本一本二区三区精品| 2021天堂中文幕一二区在线观 | 欧美人与性动交α欧美精品济南到| 嫩草影视91久久| 一级毛片精品| 国产熟女午夜一区二区三区| 哪里可以看免费的av片| 亚洲av片天天在线观看| 国产一卡二卡三卡精品| 最新美女视频免费是黄的| 欧美av亚洲av综合av国产av| 非洲黑人性xxxx精品又粗又长| 一区二区三区高清视频在线| 久久久国产成人免费| 成人av一区二区三区在线看| 男人的好看免费观看在线视频 | 美女高潮到喷水免费观看| 国产一区二区三区视频了| 国产av又大| 国产亚洲欧美在线一区二区| 动漫黄色视频在线观看| 国内精品久久久久精免费| 一级a爱视频在线免费观看| 又紧又爽又黄一区二区| 国产成人精品无人区| 激情在线观看视频在线高清| 色播亚洲综合网| 亚洲精品一卡2卡三卡4卡5卡| 麻豆av在线久日| 他把我摸到了高潮在线观看| 精品久久久久久久久久久久久 | 亚洲成人国产一区在线观看| 精品久久久久久久人妻蜜臀av| 国产成年人精品一区二区| 日韩大尺度精品在线看网址| 精品乱码久久久久久99久播| 夜夜躁狠狠躁天天躁| 久久精品影院6| 成人精品一区二区免费| 天天添夜夜摸| 亚洲最大成人中文| 国产精品久久视频播放| 中文亚洲av片在线观看爽| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 两人在一起打扑克的视频| 欧美在线一区亚洲| 国产亚洲精品第一综合不卡| 亚洲va日本ⅴa欧美va伊人久久| 十八禁人妻一区二区| 少妇 在线观看| 久久九九热精品免费| 日韩精品青青久久久久久| 欧美中文日本在线观看视频| 国产av一区二区精品久久| 免费观看精品视频网站| 国产精品自产拍在线观看55亚洲| 2021天堂中文幕一二区在线观 | 国产熟女午夜一区二区三区| 亚洲电影在线观看av| 很黄的视频免费| 亚洲国产看品久久| 淫妇啪啪啪对白视频| 国产一区二区三区在线臀色熟女| 亚洲av成人av| 中文在线观看免费www的网站 | 欧美日本亚洲视频在线播放| 色尼玛亚洲综合影院| 久久精品91蜜桃| 99久久无色码亚洲精品果冻| 777久久人妻少妇嫩草av网站| 青草久久国产| 波多野结衣高清无吗| 亚洲第一电影网av| 精品欧美一区二区三区在线| 99在线人妻在线中文字幕| 中文字幕av电影在线播放| 国内毛片毛片毛片毛片毛片| 欧美精品亚洲一区二区| 亚洲人成电影免费在线| 高清毛片免费观看视频网站| 久久久国产精品麻豆| 国产一区在线观看成人免费| 午夜视频精品福利| 欧美日韩精品网址| 欧美色视频一区免费| www日本黄色视频网| 日日爽夜夜爽网站| a级毛片在线看网站| 男人舔女人的私密视频| 亚洲精品国产一区二区精华液| 国产黄a三级三级三级人| 在线观看免费日韩欧美大片| tocl精华| 国产成人av教育| 特大巨黑吊av在线直播 | 18禁黄网站禁片免费观看直播| av欧美777| www日本黄色视频网| 精品高清国产在线一区| 欧美乱色亚洲激情| 一本大道久久a久久精品| 九色国产91popny在线| 国内久久婷婷六月综合欲色啪| 校园春色视频在线观看| 91老司机精品| 成人18禁在线播放| svipshipincom国产片| 久久 成人 亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 1024手机看黄色片| 亚洲av成人av| 亚洲片人在线观看| 久久精品aⅴ一区二区三区四区| 精品不卡国产一区二区三区| 成人国产综合亚洲| 黄色成人免费大全| 精品不卡国产一区二区三区| 中文亚洲av片在线观看爽| 女性被躁到高潮视频| 国产在线精品亚洲第一网站| 亚洲精品av麻豆狂野| 激情在线观看视频在线高清| 精品国产亚洲在线| 久久国产乱子伦精品免费另类| 亚洲男人天堂网一区| 成人免费观看视频高清| 精品久久久久久久末码| 亚洲av日韩精品久久久久久密| 每晚都被弄得嗷嗷叫到高潮| 亚洲中文字幕一区二区三区有码在线看 | 搞女人的毛片| 热re99久久国产66热| 亚洲专区字幕在线| 999久久久精品免费观看国产| 一区福利在线观看| 极品教师在线免费播放| 午夜免费激情av| 久99久视频精品免费| 男人的好看免费观看在线视频 | 亚洲一码二码三码区别大吗| 99国产精品一区二区蜜桃av| 久久九九热精品免费| 亚洲第一欧美日韩一区二区三区| 黄色a级毛片大全视频| xxx96com| 又黄又粗又硬又大视频| 欧美+亚洲+日韩+国产| 久久中文看片网| 丝袜美腿诱惑在线| 国产精品 国内视频| 最近最新中文字幕大全免费视频| 亚洲精品久久国产高清桃花| 一进一出抽搐动态| 国产aⅴ精品一区二区三区波| 美女国产高潮福利片在线看| 国产av在哪里看| 亚洲精品美女久久av网站| 国产成人精品无人区| 亚洲第一欧美日韩一区二区三区| 精品卡一卡二卡四卡免费| 亚洲五月色婷婷综合| 无遮挡黄片免费观看| 精品久久久久久久毛片微露脸| av视频在线观看入口| 两个人视频免费观看高清| 日韩三级视频一区二区三区| 免费人成视频x8x8入口观看| 久久中文看片网| 91av网站免费观看| 国产国语露脸激情在线看| 欧美 亚洲 国产 日韩一| 亚洲成人久久性| 两个人免费观看高清视频| 91av网站免费观看| 亚洲 欧美一区二区三区| 91麻豆精品激情在线观看国产| 日本 欧美在线| 怎么达到女性高潮| 国产精品亚洲一级av第二区| 亚洲国产精品sss在线观看| 久久精品91无色码中文字幕| 日本在线视频免费播放| 久久久久久人人人人人| 亚洲精品国产精品久久久不卡| 桃色一区二区三区在线观看| 中文资源天堂在线| 一区二区日韩欧美中文字幕| 中亚洲国语对白在线视频| 欧美黄色淫秽网站| 女同久久另类99精品国产91| 97人妻精品一区二区三区麻豆 | 十分钟在线观看高清视频www| 成人精品一区二区免费| 一区二区日韩欧美中文字幕| 18美女黄网站色大片免费观看| 最近最新中文字幕大全电影3 | 国产精品综合久久久久久久免费| 一二三四社区在线视频社区8| 美国免费a级毛片| 村上凉子中文字幕在线| 1024香蕉在线观看| 中文字幕精品亚洲无线码一区 | 波多野结衣av一区二区av| 久久草成人影院| 欧美国产精品va在线观看不卡| 97人妻精品一区二区三区麻豆 | 久久 成人 亚洲| 国产不卡一卡二| 国产成人精品久久二区二区91| 久99久视频精品免费| 女性被躁到高潮视频| www日本黄色视频网| 1024手机看黄色片| av天堂在线播放| 色综合站精品国产| 日本在线视频免费播放| 国产亚洲精品久久久久久毛片| 九九热线精品视视频播放| 少妇熟女欧美另类| 寂寞人妻少妇视频99o| 色吧在线观看| 国产麻豆成人av免费视频| 黄色日韩在线| 精品一区二区免费观看| 欧美色欧美亚洲另类二区| 国产伦在线观看视频一区| 国产成人91sexporn| 国产欧美日韩一区二区精品| 日韩欧美三级三区| 秋霞在线观看毛片| 真人做人爱边吃奶动态| 一a级毛片在线观看| 午夜福利在线观看吧| 国产成人a区在线观看| 国产精品综合久久久久久久免费| 最近手机中文字幕大全| 久久久久久久久中文| av免费在线看不卡| 少妇人妻精品综合一区二区 | 国产中年淑女户外野战色| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 少妇猛男粗大的猛烈进出视频 | 2021天堂中文幕一二区在线观| 精品免费久久久久久久清纯| 亚洲图色成人| 人人妻人人澡欧美一区二区| 精品免费久久久久久久清纯| 搡老熟女国产l中国老女人| 丝袜美腿在线中文| 欧美一级a爱片免费观看看| 亚洲欧美日韩无卡精品| 亚洲国产精品成人综合色| 男人狂女人下面高潮的视频| 成人特级黄色片久久久久久久| 嫩草影院新地址| 别揉我奶头 嗯啊视频| 久久久午夜欧美精品| 男女视频在线观看网站免费| 校园春色视频在线观看| 国产av在哪里看| 最近中文字幕高清免费大全6| 美女被艹到高潮喷水动态| 亚洲美女黄片视频| 国产精品爽爽va在线观看网站| 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 91av网一区二区| 午夜福利在线观看吧| 搞女人的毛片| 99国产极品粉嫩在线观看| 天堂动漫精品| 最近最新中文字幕大全电影3| av国产免费在线观看| 一进一出抽搐动态| 成人毛片a级毛片在线播放| 欧美成人a在线观看| 亚洲欧美成人精品一区二区| 成人特级av手机在线观看| 国产高清三级在线| 久久精品国产亚洲av香蕉五月| 麻豆乱淫一区二区| 日本在线视频免费播放| 国产真实伦视频高清在线观看| 亚洲国产色片| 久久精品91蜜桃| 少妇猛男粗大的猛烈进出视频 | 国内精品美女久久久久久| 久久韩国三级中文字幕| 少妇高潮的动态图| 日韩欧美免费精品| 亚洲成a人片在线一区二区| 午夜免费男女啪啪视频观看 | 国产黄片美女视频| 色综合站精品国产| 亚洲高清免费不卡视频| 最新在线观看一区二区三区| 少妇熟女欧美另类| 又黄又爽又刺激的免费视频.| 免费av不卡在线播放| 老熟妇仑乱视频hdxx| 国产欧美日韩一区二区精品| 久久久国产成人精品二区| 亚洲va在线va天堂va国产| 久久久久精品国产欧美久久久| 淫妇啪啪啪对白视频| 久久鲁丝午夜福利片| 久久精品国产自在天天线| 欧美成人a在线观看| 国产精品乱码一区二三区的特点| 欧美日韩在线观看h| 熟女人妻精品中文字幕| 国产一级毛片七仙女欲春2| 亚洲精品国产成人久久av| 精品一区二区三区av网在线观看| 3wmmmm亚洲av在线观看| 久久久久九九精品影院| 在线观看66精品国产| 精品人妻偷拍中文字幕| av免费在线看不卡| 国产人妻一区二区三区在| 尾随美女入室| 我要看日韩黄色一级片| 午夜精品一区二区三区免费看| ponron亚洲| 午夜福利视频1000在线观看| 日本熟妇午夜| 床上黄色一级片| 久久久久久久午夜电影| 久久午夜福利片| 精品福利观看| 久久久久久国产a免费观看| 一卡2卡三卡四卡精品乱码亚洲| 日本与韩国留学比较| 亚洲精品影视一区二区三区av| 小蜜桃在线观看免费完整版高清| 97超视频在线观看视频| 国内久久婷婷六月综合欲色啪| 国产精品福利在线免费观看| 国产人妻一区二区三区在| 两个人的视频大全免费| 嫩草影院入口| 亚洲国产日韩欧美精品在线观看| 亚洲精品粉嫩美女一区| 午夜精品国产一区二区电影 | 欧美一区二区亚洲| 欧美高清成人免费视频www| 99热精品在线国产| 午夜a级毛片| 国产黄色小视频在线观看| 国产乱人偷精品视频| 寂寞人妻少妇视频99o| 日韩欧美免费精品| 亚洲av中文字字幕乱码综合| 久久精品国产清高在天天线| 成年女人毛片免费观看观看9| 国产精品爽爽va在线观看网站| 国产91av在线免费观看| 日韩欧美 国产精品| 一个人免费在线观看电影| 久久精品国产鲁丝片午夜精品| 我的老师免费观看完整版| 1024手机看黄色片| 99视频精品全部免费 在线| 亚洲内射少妇av| 日日摸夜夜添夜夜添小说| 欧美bdsm另类| a级毛色黄片| 人人妻,人人澡人人爽秒播| 禁无遮挡网站| 成年版毛片免费区| av卡一久久| 亚洲av中文字字幕乱码综合| 国产精品久久久久久精品电影| 亚洲图色成人| 精品久久久久久久久亚洲| 最近2019中文字幕mv第一页| 波多野结衣高清无吗| 最近手机中文字幕大全| 人人妻人人澡人人爽人人夜夜 | 在线观看一区二区三区| 亚洲精品456在线播放app| 高清午夜精品一区二区三区 | 色综合站精品国产| 国产欧美日韩精品一区二区| 99久久中文字幕三级久久日本| 久久久欧美国产精品| 欧美激情国产日韩精品一区| 美女 人体艺术 gogo| 午夜亚洲福利在线播放| 一本久久中文字幕| 亚洲第一电影网av| 又黄又爽又刺激的免费视频.| av天堂中文字幕网| 少妇人妻一区二区三区视频| 男女做爰动态图高潮gif福利片| 久久精品夜夜夜夜夜久久蜜豆| 免费大片18禁| 俺也久久电影网| 欧美成人免费av一区二区三区| 亚洲,欧美,日韩| av国产免费在线观看| 国产精品日韩av在线免费观看| 亚洲最大成人手机在线| 国产又黄又爽又无遮挡在线| 国产精品一区二区性色av| 免费看美女性在线毛片视频| 国内精品久久久久精免费| av在线亚洲专区| 国产不卡一卡二| 成人性生交大片免费视频hd| 国产精品亚洲一级av第二区| 欧美成人a在线观看| 国产视频一区二区在线看| 久久久精品欧美日韩精品| 大型黄色视频在线免费观看| 一本久久中文字幕| 永久网站在线| 国产高清有码在线观看视频| 国产探花极品一区二区| videossex国产| 不卡视频在线观看欧美| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 国产精品国产三级国产av玫瑰| 在线看三级毛片| 精品一区二区三区视频在线观看免费| 日日摸夜夜添夜夜添小说| 99久久中文字幕三级久久日本| 人妻久久中文字幕网| 婷婷精品国产亚洲av在线| 亚洲精品国产av成人精品 | 丰满的人妻完整版| 在现免费观看毛片| 亚洲第一区二区三区不卡| 成人无遮挡网站| 波多野结衣巨乳人妻| 国产私拍福利视频在线观看| 久久精品国产亚洲av天美| 深夜精品福利| 又爽又黄a免费视频| 国产视频一区二区在线看| 国产亚洲精品久久久久久毛片| 欧美性猛交╳xxx乱大交人| 国产精品人妻久久久久久| 国产男靠女视频免费网站| 日韩欧美免费精品| 精品国内亚洲2022精品成人| 男人舔奶头视频| 午夜福利成人在线免费观看| 国产精品一区二区性色av| 亚洲国产欧美人成| 极品教师在线视频| 嫩草影视91久久| 国产精品久久久久久av不卡| 国产精品久久久久久久久免| 大型黄色视频在线免费观看| 嫩草影院新地址| 亚洲在线观看片| 一进一出好大好爽视频| 亚洲va在线va天堂va国产| 人人妻人人看人人澡| 99久久精品热视频| 欧美成人精品欧美一级黄| 亚洲成人av在线免费| 成人特级av手机在线观看| 国产日本99.免费观看| 国产熟女欧美一区二区| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| 波野结衣二区三区在线| 欧美性感艳星| 最近最新中文字幕大全电影3| 97人妻精品一区二区三区麻豆| 色视频www国产| 国产精品久久久久久亚洲av鲁大| 久久久久国内视频| 亚洲第一区二区三区不卡| 一级黄片播放器| 久久久久久九九精品二区国产| 夜夜夜夜夜久久久久| 国产精品一区二区性色av| 色视频www国产| 久久人人精品亚洲av| 久久久久久久久久久丰满| 午夜a级毛片| 啦啦啦观看免费观看视频高清| 国产男人的电影天堂91| 91在线观看av| 亚洲精品国产成人久久av| 亚洲av五月六月丁香网| 99riav亚洲国产免费| 男人舔奶头视频| 久久人人爽人人片av| 国产一区二区在线观看日韩| 免费黄网站久久成人精品| 天堂√8在线中文| 国产乱人偷精品视频| 欧美bdsm另类| 成人性生交大片免费视频hd| 国内久久婷婷六月综合欲色啪| 欧美色欧美亚洲另类二区| 伦精品一区二区三区| 永久网站在线| 色哟哟·www| 欧美xxxx性猛交bbbb| 大香蕉久久网| 91在线精品国自产拍蜜月| 精品一区二区三区av网在线观看| .国产精品久久| 狂野欧美激情性xxxx在线观看| 国产精品一区二区免费欧美| 中文字幕av成人在线电影| 91久久精品国产一区二区成人| 国产精品野战在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久精品夜色国产| 真人做人爱边吃奶动态| 日韩制服骚丝袜av| 在现免费观看毛片| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜添小说| 色吧在线观看| 国产精品爽爽va在线观看网站| 免费黄网站久久成人精品| 久久精品91蜜桃| 欧美xxxx黑人xx丫x性爽| 女人被狂操c到高潮| 真人做人爱边吃奶动态| 国产蜜桃级精品一区二区三区| 99在线视频只有这里精品首页| 国产精品1区2区在线观看.| 久久精品国产亚洲av香蕉五月| 国产精品乱码一区二三区的特点| 黑人高潮一二区| 欧美色欧美亚洲另类二区| 精品一区二区三区人妻视频| 天天躁日日操中文字幕| 亚洲国产日韩欧美精品在线观看| 日本色播在线视频| 欧美成人精品欧美一级黄| 日韩av不卡免费在线播放| 欧美zozozo另类| 欧美日韩综合久久久久久| 亚洲精品国产av成人精品 | 欧美三级亚洲精品| 国产69精品久久久久777片| 3wmmmm亚洲av在线观看| 六月丁香七月| av天堂在线播放| 国模一区二区三区四区视频| 精品一区二区三区人妻视频| 成人亚洲精品av一区二区| 精品一区二区三区视频在线观看免费| 日韩亚洲欧美综合| 久久久午夜欧美精品| 老女人水多毛片| 午夜a级毛片| 久久99热6这里只有精品| 天天躁夜夜躁狠狠久久av| 中文字幕免费在线视频6| 我的女老师完整版在线观看| 国内精品久久久久精免费| 丝袜美腿在线中文| 日本黄色视频三级网站网址| 亚洲四区av| 免费看光身美女| 18禁裸乳无遮挡免费网站照片| 国产精品乱码一区二三区的特点| 九九热线精品视视频播放| 国产欧美日韩精品一区二区| 91麻豆精品激情在线观看国产| 高清午夜精品一区二区三区 | 国产亚洲精品av在线| 久久精品国产亚洲av香蕉五月| 日本黄色片子视频| 日韩一本色道免费dvd| 亚洲七黄色美女视频| 成人一区二区视频在线观看| 美女黄网站色视频| a级毛色黄片| 日韩人妻高清精品专区| 亚洲aⅴ乱码一区二区在线播放| 色尼玛亚洲综合影院| 日本黄大片高清| 99riav亚洲国产免费| 大又大粗又爽又黄少妇毛片口| 久久99热这里只有精品18| 国产精品一区二区免费欧美| 国产一区二区三区av在线 | 亚洲美女黄片视频| 国产中年淑女户外野战色| 最好的美女福利视频网| 免费观看精品视频网站| 男女边吃奶边做爰视频| 国产精品1区2区在线观看.| 亚洲高清免费不卡视频| 欧美日韩在线观看h| 99久久成人亚洲精品观看| 国产精品一区二区三区四区免费观看 | 成人永久免费在线观看视频| 免费看美女性在线毛片视频| 欧美日韩综合久久久久久| а√天堂www在线а√下载| 亚洲专区国产一区二区| 乱系列少妇在线播放| 身体一侧抽搐| 精品久久久久久成人av| 高清午夜精品一区二区三区 | 久久久久久久久久久丰满|