• <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麻豆| 欧美+日韩+精品| 性欧美人与动物交配| 狠狠狠狠99中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 日日摸夜夜添夜夜爱| av国产免费在线观看| 色综合色国产| 国产高清三级在线| 麻豆成人午夜福利视频| 在线观看美女被高潮喷水网站| 美女大奶头视频| 欧美日韩综合久久久久久| 黄色欧美视频在线观看| a级毛色黄片| 亚洲婷婷狠狠爱综合网| 天堂√8在线中文| 最好的美女福利视频网| 身体一侧抽搐| 婷婷亚洲欧美| 国产伦精品一区二区三区四那| 一级av片app| 久久午夜福利片| 精品久久久久久久久av| 亚洲乱码一区二区免费版| 可以在线观看的亚洲视频| 久久久久久九九精品二区国产| 内射极品少妇av片p| 亚洲欧美日韩高清在线视频| 99热全是精品| 一级毛片我不卡| 亚洲精品国产成人久久av| 日韩精品有码人妻一区| 亚洲国产高清在线一区二区三| 国产白丝娇喘喷水9色精品| 小说图片视频综合网站| 不卡一级毛片| 午夜福利高清视频| 大型黄色视频在线免费观看| 久久99精品国语久久久| 黄片wwwwww| 久久中文看片网| 91麻豆精品激情在线观看国产| 色综合亚洲欧美另类图片| 黄片wwwwww| 精品人妻偷拍中文字幕| 亚洲色图av天堂| 亚洲,欧美,日韩| 国产大屁股一区二区在线视频| 中文字幕免费在线视频6| 一本一本综合久久| 性插视频无遮挡在线免费观看| 99久久精品热视频| 悠悠久久av| 亚洲国产精品合色在线| 亚洲不卡免费看| 亚洲欧美日韩高清专用| 99九九线精品视频在线观看视频| 亚洲av免费在线观看| 色噜噜av男人的天堂激情| 久久久久久九九精品二区国产| 免费av不卡在线播放| 91狼人影院| 国产在视频线在精品| 国产精品蜜桃在线观看 | 成人性生交大片免费视频hd| 免费搜索国产男女视频| 少妇的逼好多水| 久久精品国产自在天天线| 久久国产乱子免费精品| 日韩精品青青久久久久久| 99久久九九国产精品国产免费| 尾随美女入室| 国产精品野战在线观看| 国产精品久久久久久久久免| 麻豆成人午夜福利视频| av卡一久久| 狂野欧美激情性xxxx在线观看| 亚洲av成人av| 欧美日本亚洲视频在线播放| 黄色一级大片看看| 国产av麻豆久久久久久久| 级片在线观看| 国产乱人偷精品视频| 美女cb高潮喷水在线观看| 深爱激情五月婷婷| 麻豆国产97在线/欧美| 中国国产av一级| 免费人成视频x8x8入口观看| 亚洲欧美清纯卡通| 午夜老司机福利剧场| 好男人视频免费观看在线| 亚洲欧美日韩高清专用| 国产精品久久久久久久久免| 欧美另类亚洲清纯唯美| 免费看日本二区| 免费看av在线观看网站| 欧美xxxx性猛交bbbb| 男人和女人高潮做爰伦理| 成人午夜高清在线视频| 国产成年人精品一区二区| 免费观看a级毛片全部| 成年女人看的毛片在线观看| 欧美日本视频| 真实男女啪啪啪动态图| 成人美女网站在线观看视频| 国产精品永久免费网站| 高清在线视频一区二区三区 | 黄片wwwwww| 我的女老师完整版在线观看| 色吧在线观看| 色综合站精品国产| 校园人妻丝袜中文字幕| 女人十人毛片免费观看3o分钟| 一区二区三区四区激情视频 | 国产色婷婷99| 成人一区二区视频在线观看| 国产精品一区www在线观看| 91在线精品国自产拍蜜月| 亚洲最大成人av| 久久久久久大精品| 有码 亚洲区| 亚洲一级一片aⅴ在线观看| 午夜a级毛片| 国产精品无大码| 蜜臀久久99精品久久宅男| 久久热精品热| 不卡一级毛片| 亚洲欧洲国产日韩| 国产成人福利小说| 精品午夜福利在线看| av免费在线看不卡| 久久精品国产亚洲av香蕉五月| 男人舔女人下体高潮全视频| 精品人妻视频免费看| 午夜福利在线在线| 国产乱人视频| 成人毛片a级毛片在线播放| 亚洲欧美精品自产自拍| 国产精品久久久久久久电影| 国产老妇女一区| 在线观看66精品国产| www.av在线官网国产| 日本熟妇午夜| 熟妇人妻久久中文字幕3abv| 久久久久免费精品人妻一区二区| 欧美变态另类bdsm刘玥| 国产伦精品一区二区三区四那| 精品国内亚洲2022精品成人| 老女人水多毛片| 中文字幕久久专区| 国产色婷婷99| 亚洲欧美日韩高清在线视频| 又黄又爽又刺激的免费视频.| 成人美女网站在线观看视频| 成人美女网站在线观看视频| 人人妻人人看人人澡| 午夜免费男女啪啪视频观看| 赤兔流量卡办理| 国产成人一区二区在线| 亚洲精华国产精华液的使用体验 | 亚洲精品久久久久久婷婷小说 | 国产一级毛片在线| 亚洲一区高清亚洲精品| 国产美女午夜福利| 色视频www国产| 欧美另类亚洲清纯唯美| 99热这里只有精品一区| 成人鲁丝片一二三区免费| 性欧美人与动物交配| 成人二区视频| 三级毛片av免费| 性欧美人与动物交配| 国产精品伦人一区二区| 成人鲁丝片一二三区免费| 中文字幕免费在线视频6| 直男gayav资源| 直男gayav资源| 精品99又大又爽又粗少妇毛片| 欧美极品一区二区三区四区| 国产大屁股一区二区在线视频| 99热这里只有是精品50| 免费人成在线观看视频色| 欧美人与善性xxx| 国产成人freesex在线| 成人欧美大片| 中文字幕人妻熟人妻熟丝袜美| 日本三级黄在线观看| 亚洲欧美日韩高清在线视频| 青春草视频在线免费观看| 免费一级毛片在线播放高清视频| 亚洲丝袜综合中文字幕| 欧美一级a爱片免费观看看| 人人妻人人看人人澡| 九草在线视频观看| 国内精品一区二区在线观看| 国产精品1区2区在线观看.| 老熟妇乱子伦视频在线观看| 国产成人freesex在线| 欧美+日韩+精品| 欧美+日韩+精品| 大香蕉久久网| 免费av不卡在线播放| 黄色一级大片看看| 欧美+亚洲+日韩+国产| av免费观看日本| 亚洲精品色激情综合| 观看美女的网站| 插阴视频在线观看视频| 一进一出抽搐动态| 亚洲在线观看片| 日本熟妇午夜| 精品熟女少妇av免费看| 国产三级在线视频| 黄片wwwwww| 久久久久久久久久黄片| 久久精品国产鲁丝片午夜精品| 免费黄网站久久成人精品| 国产又黄又爽又无遮挡在线| 国内精品久久久久精免费| 国产亚洲精品av在线| 色综合色国产| 亚洲国产高清在线一区二区三| 亚洲激情五月婷婷啪啪| 精品久久久久久久久亚洲| 久久久久久伊人网av| 国产成人福利小说| 国产高清有码在线观看视频| 99久久精品国产国产毛片| 又爽又黄a免费视频| 成人午夜精彩视频在线观看| 看黄色毛片网站| 国产免费男女视频| h日本视频在线播放| 久久精品91蜜桃| 午夜福利在线在线| 久久久久久伊人网av| 国产一区二区亚洲精品在线观看| 日日干狠狠操夜夜爽| 国产成人影院久久av| 99热全是精品| 国产av一区在线观看免费| 欧美最黄视频在线播放免费| 伦精品一区二区三区| 婷婷精品国产亚洲av| ponron亚洲| 色哟哟·www| 免费搜索国产男女视频| 51国产日韩欧美| 蜜桃久久精品国产亚洲av| 啦啦啦观看免费观看视频高清| 特级一级黄色大片| 99视频精品全部免费 在线| 麻豆久久精品国产亚洲av| 美女大奶头视频| 91麻豆精品激情在线观看国产| 久久精品综合一区二区三区| 校园人妻丝袜中文字幕| 久久久久九九精品影院| av女优亚洲男人天堂| 久久草成人影院| 国产又黄又爽又无遮挡在线| 中文资源天堂在线| 婷婷精品国产亚洲av| 日产精品乱码卡一卡2卡三| 国产精品久久久久久精品电影| 国产男人的电影天堂91| а√天堂www在线а√下载| 国产日本99.免费观看| 麻豆av噜噜一区二区三区| 最近最新中文字幕大全电影3| 国内精品一区二区在线观看| 亚洲无线观看免费| 久久精品国产亚洲av涩爱 | 日韩精品青青久久久久久| 欧美日韩精品成人综合77777| 国产av在哪里看| 亚洲欧美成人综合另类久久久 | 高清午夜精品一区二区三区 | 深爱激情五月婷婷| 国产精品精品国产色婷婷| 欧美zozozo另类| 久99久视频精品免费| 国产v大片淫在线免费观看| 亚洲国产日韩欧美精品在线观看| 亚洲欧美精品自产自拍| 亚洲精品日韩av片在线观看| 久久久精品94久久精品| 99久久精品国产国产毛片| 日本三级黄在线观看| 久久鲁丝午夜福利片| 黄片wwwwww| 两个人视频免费观看高清| 非洲黑人性xxxx精品又粗又长| 久久久精品大字幕| 精品不卡国产一区二区三区| 亚洲精品成人久久久久久| 欧美一区二区国产精品久久精品| 好男人视频免费观看在线| 欧美丝袜亚洲另类| 最近最新中文字幕大全电影3| 不卡视频在线观看欧美| 日韩一本色道免费dvd| 久久精品国产清高在天天线| 69人妻影院| 午夜激情福利司机影院| 最近最新中文字幕大全电影3| 日韩,欧美,国产一区二区三区 | 日本熟妇午夜| 91午夜精品亚洲一区二区三区| av免费在线看不卡| 69人妻影院| 99热6这里只有精品| 成人美女网站在线观看视频| 网址你懂的国产日韩在线| 神马国产精品三级电影在线观看| 爱豆传媒免费全集在线观看| 国产高清有码在线观看视频| 一级黄片播放器| 青春草视频在线免费观看| 夫妻性生交免费视频一级片| 国产 一区精品| 青春草亚洲视频在线观看| 直男gayav资源| 91午夜精品亚洲一区二区三区| 亚洲七黄色美女视频| 国产久久久一区二区三区| 美女脱内裤让男人舔精品视频 | 18禁在线无遮挡免费观看视频| 联通29元200g的流量卡| 26uuu在线亚洲综合色| 亚洲四区av| 欧美+日韩+精品| 亚洲乱码一区二区免费版| 亚洲一级一片aⅴ在线观看| 女人十人毛片免费观看3o分钟| 啦啦啦啦在线视频资源| 波多野结衣巨乳人妻| 久久久午夜欧美精品| 狂野欧美激情性xxxx在线观看| 国产真实乱freesex| 欧美激情国产日韩精品一区| 久久精品久久久久久噜噜老黄 | 91麻豆精品激情在线观看国产| 亚洲七黄色美女视频| 国产中年淑女户外野战色| 成人一区二区视频在线观看| 国产精品一区二区性色av| 日韩中字成人| 天堂√8在线中文| 99久国产av精品国产电影| av黄色大香蕉| 你懂的网址亚洲精品在线观看 | 久久亚洲精品不卡| 天天躁日日操中文字幕| 又爽又黄a免费视频| 九草在线视频观看| 成人三级黄色视频| 男女边吃奶边做爰视频| 亚洲国产精品成人久久小说 | www日本黄色视频网| 亚洲av一区综合| 乱系列少妇在线播放| 欧美最黄视频在线播放免费| 午夜福利高清视频| 深爱激情五月婷婷| 一级av片app| 黄色配什么色好看| 少妇被粗大猛烈的视频| 午夜福利在线观看吧| 亚洲天堂国产精品一区在线| 白带黄色成豆腐渣| 久久久成人免费电影| 97热精品久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 一级黄色大片毛片| 国产不卡一卡二| 日韩成人av中文字幕在线观看| 久久99热6这里只有精品| 欧美变态另类bdsm刘玥| 男女啪啪激烈高潮av片| 精品久久久久久久久av| 国产精品嫩草影院av在线观看| 国产色爽女视频免费观看| 夫妻性生交免费视频一级片| 国国产精品蜜臀av免费| 人妻少妇偷人精品九色| 久久久久久久久中文| 只有这里有精品99| 久久综合国产亚洲精品| 亚洲一区高清亚洲精品| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清激情床上av| 18禁在线无遮挡免费观看视频| 亚洲精品乱码久久久v下载方式| 波野结衣二区三区在线| 九九热线精品视视频播放| 禁无遮挡网站| 国产成人a区在线观看| 日韩国内少妇激情av| 热99在线观看视频| 午夜老司机福利剧场| 亚洲av成人精品一区久久| 欧美一区二区精品小视频在线| 欧美性猛交╳xxx乱大交人| 中出人妻视频一区二区| 伦理电影大哥的女人| 淫秽高清视频在线观看| 成人鲁丝片一二三区免费| 免费人成在线观看视频色| 亚洲色图av天堂| 免费无遮挡裸体视频| 亚洲成人久久性| 国产av在哪里看| 精品午夜福利在线看| 久久久久久久午夜电影| 欧美zozozo另类| 乱系列少妇在线播放| 直男gayav资源| 人妻制服诱惑在线中文字幕| 免费观看在线日韩| 国产成人91sexporn| 嫩草影院精品99| 熟女人妻精品中文字幕| 亚洲三级黄色毛片| 欧美日韩一区二区视频在线观看视频在线 | 日本爱情动作片www.在线观看| 男插女下体视频免费在线播放| 99热这里只有是精品在线观看| 国产综合懂色| 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美 国产精品| 日韩三级伦理在线观看| 一级黄片播放器| www.av在线官网国产| 亚洲av.av天堂| 狂野欧美白嫩少妇大欣赏| 黄色配什么色好看| 一边亲一边摸免费视频| 性欧美人与动物交配| 一级黄色大片毛片| 国国产精品蜜臀av免费| 免费看日本二区| 免费观看在线日韩| 国产视频首页在线观看| 国内久久婷婷六月综合欲色啪| 欧美精品一区二区大全| 观看美女的网站| 国产在视频线在精品| 91在线精品国自产拍蜜月| 国产精品永久免费网站| 国内少妇人妻偷人精品xxx网站| 黑人高潮一二区| 麻豆精品久久久久久蜜桃| 欧美最新免费一区二区三区| 国产成人福利小说| 成人亚洲精品av一区二区| 国产午夜福利久久久久久| or卡值多少钱| 久久精品国产鲁丝片午夜精品| ponron亚洲| av视频在线观看入口| 亚洲欧美日韩高清在线视频| 国产精品无大码| 最近手机中文字幕大全| 在线播放无遮挡| 久久婷婷人人爽人人干人人爱| 国产精品一二三区在线看| 天天躁夜夜躁狠狠久久av| 久久久久久久久久成人| 亚洲成人久久性| 狠狠狠狠99中文字幕| 成人综合一区亚洲| 日本色播在线视频| 一本一本综合久久| 中文资源天堂在线| 久久中文看片网| 欧美3d第一页| 亚洲中文字幕一区二区三区有码在线看| 偷拍熟女少妇极品色| 国产成人精品一,二区 | 性色avwww在线观看| 免费观看精品视频网站| 亚洲精品日韩av片在线观看| 国产片特级美女逼逼视频| 亚洲国产精品合色在线| 国产成人精品婷婷| 精品久久久久久成人av| 在线观看av片永久免费下载| 亚洲精品影视一区二区三区av| 国产成人午夜福利电影在线观看| 国内揄拍国产精品人妻在线| 亚洲国产精品成人久久小说 | 国国产精品蜜臀av免费| 特大巨黑吊av在线直播| av.在线天堂| 国产精品女同一区二区软件| 精品久久久久久久久久免费视频| 又粗又爽又猛毛片免费看| 男女视频在线观看网站免费| 波多野结衣巨乳人妻| 久久久成人免费电影| 99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 国产黄片视频在线免费观看| 黄色欧美视频在线观看| 久久人人精品亚洲av| 日本免费a在线| 欧美xxxx性猛交bbbb| 国产精品电影一区二区三区| 成人毛片60女人毛片免费| 国产精品一区二区三区四区久久| 久久久久九九精品影院| 真实男女啪啪啪动态图| 三级国产精品欧美在线观看| 一边亲一边摸免费视频| 国产日本99.免费观看| 熟妇人妻久久中文字幕3abv| 亚洲精品日韩av片在线观看| av视频在线观看入口| 日本一本二区三区精品| 男的添女的下面高潮视频| 久久久久久大精品| 久久久国产成人精品二区| 伦理电影大哥的女人| 色综合站精品国产| 日韩高清综合在线| 精品久久久噜噜| 国产三级中文精品| 综合色丁香网| 男人狂女人下面高潮的视频| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 国产成人freesex在线| 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 亚洲综合色惰| 亚洲人成网站在线播| 欧美日韩在线观看h| 在线免费观看不下载黄p国产| 国内精品美女久久久久久| 18禁裸乳无遮挡免费网站照片| 九草在线视频观看| 亚洲欧美日韩卡通动漫| 不卡视频在线观看欧美| 一级二级三级毛片免费看| 午夜视频国产福利| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 久久精品国产清高在天天线| 蜜桃久久精品国产亚洲av| www日本黄色视频网| 女人被狂操c到高潮| av天堂在线播放| 最近的中文字幕免费完整| 99久久久亚洲精品蜜臀av| 亚洲四区av| 国产极品精品免费视频能看的| 六月丁香七月| 黄色日韩在线| 中国国产av一级| 人人妻人人看人人澡| 国产精品人妻久久久影院| 中文在线观看免费www的网站| 日韩精品有码人妻一区| 成熟少妇高潮喷水视频| 欧美日韩精品成人综合77777| 又爽又黄无遮挡网站| 久久99精品国语久久久| 久久久久性生活片| 此物有八面人人有两片| 麻豆乱淫一区二区| 国产 一区精品| 99热这里只有精品一区| 久久这里只有精品中国| 国产在线精品亚洲第一网站| 精品国产三级普通话版| 国产综合懂色| 久久99热这里只有精品18| 国产精品永久免费网站| 国产蜜桃级精品一区二区三区| 日韩欧美国产在线观看| 久久精品国产清高在天天线| 精品人妻一区二区三区麻豆| 美女国产视频在线观看| 精品久久久久久久久久久久久| 九九久久精品国产亚洲av麻豆| 内地一区二区视频在线| 国产伦精品一区二区三区视频9| 亚洲图色成人| 国产蜜桃级精品一区二区三区| 国产精品永久免费网站| av免费在线看不卡| 日本一本二区三区精品| 亚洲自偷自拍三级| 国产精品电影一区二区三区| 成年女人永久免费观看视频| 97超碰精品成人国产| 美女xxoo啪啪120秒动态图| 国产极品精品免费视频能看的| 日本欧美国产在线视频| 久久久午夜欧美精品| 中出人妻视频一区二区| 国产精品久久久久久精品电影| 亚洲丝袜综合中文字幕| 欧美一区二区亚洲| 国产精品乱码一区二三区的特点| 一区二区三区四区激情视频 | 久久中文看片网| 亚洲三级黄色毛片| 综合色av麻豆| 国产精品乱码一区二三区的特点| 午夜福利在线观看吧|