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

    小相嶺山系大熊貓廊道規(guī)劃

    2016-05-10 02:33:19楊志松齊敦武楊旭煜古曉東強(qiáng)1西華師范大學(xué)大熊貓自然保護(hù)與文化研究中心南充700中國(guó)科學(xué)院成都生物研究所成都1001四川雪寶頂國(guó)家級(jí)自然保護(hù)區(qū)綿陽(yáng)0四川大學(xué)生命科學(xué)學(xué)院自然生物資源與生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室成都100四川省瀕危野生動(dòng)物保護(hù)生物學(xué)省部共建實(shí)驗(yàn)室成都大熊貓繁育研究基地成都1008四川省野生動(dòng)物資源調(diào)查保護(hù)管理站成都10081
    生態(tài)學(xué)報(bào) 2016年4期

    青 菁,胥 池,楊 彪,楊志松,齊敦武,楊旭煜,古曉東,戴 強(qiáng)1西華師范大學(xué)大熊貓自然保護(hù)與文化研究中心,南充 700中國(guó)科學(xué)院成都生物研究所,成都 1001四川雪寶頂國(guó)家級(jí)自然保護(hù)區(qū),綿陽(yáng) 0四川大學(xué)生命科學(xué)學(xué)院自然生物資源與生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,成都 100四川省瀕危野生動(dòng)物保護(hù)生物學(xué)省部共建實(shí)驗(yàn)室,成都大熊貓繁育研究基地,成都 1008四川省野生動(dòng)物資源調(diào)查保護(hù)管理站,成都 10081

    ?

    小相嶺山系大熊貓廊道規(guī)劃

    青菁1,2,胥池3,楊彪4,楊志松1,*,齊敦武5,楊旭煜6,古曉東6,戴強(qiáng)2
    1西華師范大學(xué)大熊貓自然保護(hù)與文化研究中心,南充637002
    2中國(guó)科學(xué)院成都生物研究所,成都610041
    3四川雪寶頂國(guó)家級(jí)自然保護(hù)區(qū),綿陽(yáng)622550
    4四川大學(xué)生命科學(xué)學(xué)院自然生物資源與生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室,成都610064
    5四川省瀕危野生動(dòng)物保護(hù)生物學(xué)省部共建實(shí)驗(yàn)室,成都大熊貓繁育研究基地,成都610086
    6四川省野生動(dòng)物資源調(diào)查保護(hù)管理站,成都610081

    摘要:廊道能有效地增加棲息地斑塊間的大熊貓種群交流,降低種群滅絕的風(fēng)險(xiǎn)。小相嶺山系是現(xiàn)存大熊貓種群數(shù)量最少的山系之一。通過(guò)MAXENT(Maximum Entropy)分析,發(fā)現(xiàn)小相嶺山系存在7個(gè)相對(duì)隔離的棲息地斑塊,棲息地破碎化嚴(yán)重。為了增加小相嶺山系大熊貓種群交流,利用最小代價(jià)路徑原理和電流理論,從保護(hù)管理角度出發(fā),在該區(qū)域劃定了9條連接大熊貓棲息地斑塊的廊道。其中紫馬-麻麻地廊道連接該區(qū)域面積最大、大熊貓痕跡點(diǎn)最多的兩個(gè)棲息地斑塊,生態(tài)意義重大,建議優(yōu)先建設(shè)該廊道。此外,研究劃定的多個(gè)廊道彼此鄰近,甚至存在一定程度的重疊。將這些極其鄰近或者重疊的廊道合并,并提出“廊道群”的概念。廊道群是兩個(gè)或兩個(gè)以上極其相鄰的廊道構(gòu)成的網(wǎng)狀結(jié)構(gòu),廊道群作為景觀上的網(wǎng)狀結(jié)構(gòu),能同時(shí)連接多個(gè)棲息地斑塊,更有效地增加多個(gè)棲息地斑塊間的種群交流。由于農(nóng)田的擠壓,部分廊道難以達(dá)到大熊貓廊道的最低寬度需求(2230 m),制約了其生態(tài)效應(yīng)的發(fā)揮,建議在以后的保護(hù)工作中探索農(nóng)田區(qū)域的棲息地恢復(fù)方式。

    關(guān)鍵詞:大熊貓;小相嶺;最小代價(jià)路徑;電流理論;廊道群

    青菁,胥池,楊彪,楊志松,齊敦武,楊旭煜,古曉東,戴強(qiáng).小相嶺山系大熊貓廊道規(guī)劃.生態(tài)學(xué)報(bào),2016,36(4):1125-1133.

    Qing J,Xu C,Yang B,Yang Z S,Qi D W,Yang X Y,Gu X D,Dai Q.Corridor design for the giant panda in the Xiaoxiangling Mountains.Acta Ecologica Sinica,2016,36(4):1125-1133.

    由于人類活動(dòng)干擾,野生動(dòng)物棲息地呈現(xiàn)破碎化,種群被隔離為孤立的小種群,彼此間缺乏基因交流[1]。通過(guò)構(gòu)建野生動(dòng)物廊道則可以增加棲息地之間的連接性,加強(qiáng)種群間的交流,有助于緩解棲息地破碎化帶來(lái)的負(fù)面影響[2]。野生動(dòng)物廊道(Wildlife Corridor)提出于1974年[3],廊道是連接野生動(dòng)物棲息地的線性景觀要素[4],或是拉伸的面狀景觀要素[5]。大量研究表明,野生動(dòng)物廊道能有效增加種群間的基因交流,防止近交衰退,降低種群的滅絕風(fēng)險(xiǎn)[6-7]。目前,常用的廊道規(guī)劃的方法主要有4個(gè):專家意見模型(the expert model)[8]、最小代價(jià)路徑分析(least-cost path analysis,LCP analysis)[9-10]、電流分析(circuit analysis)[11]、基于個(gè)體的模型(individual-based models)[12]等。

    大熊貓(Ailuropoda melanoleuca)是我國(guó)所特有的珍稀瀕危動(dòng)物,現(xiàn)僅分布于秦嶺、岷山、邛崍山、涼山、大相嶺及小相嶺六大山系。小相嶺山系位于中國(guó)大熊貓分布的西南端,是現(xiàn)有大熊貓種群數(shù)量最少的山系之一。由于墾荒、砍伐、放牧等人為干擾,小相嶺山系大熊貓棲息地急劇縮小且破碎化嚴(yán)重。其中108國(guó)道從南向北將小相嶺大熊貓棲息地分隔為東西兩塊。近250年來(lái),小相嶺大熊貓種群下降了約60倍[13]。朱磊等[14]認(rèn)為如果缺乏基因交流或者外來(lái)補(bǔ)充個(gè)體,這樣的小種群在80a內(nèi)會(huì)趨于滅絕。

    為加強(qiáng)小相嶺大熊貓種群保護(hù),增進(jìn)棲息地間種群交流,本文基于最小代價(jià)路徑原理和電流理論,在評(píng)估小相嶺山系的棲息地適宜度,確定棲息地斑塊基礎(chǔ)上,對(duì)這些斑塊間的廊道進(jìn)行規(guī)劃。

    1 研究區(qū)域

    小相嶺山系位于大涼山之西,東起甘洛河和越西河,西至雅礱江和九龍河,北抵大渡河,南到喜德縣紅莫山北面的太皇山。南北長(zhǎng)約147 km,東西寬約107 km,面積約為11500 km2,山體呈南北走向。小相嶺山系處于四川東部臺(tái)地與西部地槽過(guò)渡帶,因東部古陸地的長(zhǎng)期隆起,屬于青藏高原的組成部分。地貌以中高山為主,兼有部分低山和河谷階地。土壤屬四川盆地濕潤(rùn)亞熱帶森林土壤區(qū)。以亞熱帶季風(fēng)為基帶的山地氣候,年均溫約11.7—14.4℃,年降水量800—1250 mm。植被帶譜大致為:海拔2000 m以下為旱生河谷灌叢; 2000—2400 m為山地闊葉林;2400—3200 m或4000 m為山地暗針葉林或亞高山暗針葉林[15]。

    2 研究方法

    2.1方法和原理

    本研究使用Maximum Entropy(MAXENT)模型評(píng)估研究區(qū)域大熊貓棲息地質(zhì)量,得到棲息地適宜度指數(shù)(habitat suitability index,HSI),并確定大熊貓棲息地斑塊。MAXENT模型是根據(jù)已知的物種分布點(diǎn)和棲息地生態(tài)環(huán)境數(shù)據(jù)來(lái)評(píng)估棲息地質(zhì)量和預(yù)測(cè)物種的可能分布范圍[16-18]?;跅⒌刭|(zhì)量確定阻力系數(shù),構(gòu)建阻力圖層。阻力系數(shù)是表征野生動(dòng)物穿過(guò)特定環(huán)境的可能或者意愿的指數(shù)。在穿過(guò)特定環(huán)境過(guò)程中,動(dòng)物個(gè)體能量和時(shí)間消耗少,死亡率低,則意味著該環(huán)境阻力低,動(dòng)物從該區(qū)域通過(guò)可能大[19]。棲息地適宜度指數(shù)可以近似作為阻力系數(shù),即假設(shè)棲息地質(zhì)量越好,動(dòng)物通過(guò)的阻力越低,動(dòng)物越可能從該區(qū)域通過(guò)[20]。本文主要基于最小代價(jià)路徑原理(LCP model)和電流理論模型(Circuit theory model)進(jìn)行廊道規(guī)劃。LCP原理假設(shè)個(gè)體優(yōu)先選擇阻力小(一般假設(shè)為棲息地質(zhì)量好)的區(qū)域進(jìn)行遷移活動(dòng),以此計(jì)算代價(jià)距離,以距離代價(jià)最小(Least-cost distance,LCD)的區(qū)域作為廊道[21]。電流理論是將阻力圖層模擬成平面導(dǎo)體,即阻力越大,電流越小,以此來(lái)模擬個(gè)體移動(dòng)的密集程度,以電流主要通過(guò)的路徑作為大熊貓通過(guò)的熱點(diǎn)區(qū)域[11]。

    2.2數(shù)據(jù)來(lái)源

    大熊貓分布點(diǎn)數(shù)據(jù)來(lái)源于全國(guó)第四次大熊貓調(diào)查及研究區(qū)域長(zhǎng)期野外監(jiān)測(cè)數(shù)據(jù)。根據(jù)大熊貓最小家域面積3.9 km2[22],估算大熊貓的最小活動(dòng)半徑為1125 m,為避免空間自相關(guān)影響模型預(yù)測(cè),剔除距離小于1125 m的冗余點(diǎn),最終103個(gè)大熊貓分布點(diǎn)進(jìn)入分析運(yùn)算。

    大量研究表明,大熊貓對(duì)棲息地的坡型、坡向、坡位等地形因素以及植被因素等生態(tài)環(huán)境因子具有選擇性[23-24],并且水源距離、可食竹的分布以及人類干擾狀況也會(huì)影響大熊貓的分布[25-26],故選取以下生態(tài)環(huán)境因子數(shù)據(jù)評(píng)估其棲息地適宜性:1)地形數(shù)據(jù):包括海拔、坡度、坡向、坡位、坡型、太陽(yáng)輻射指數(shù),基于來(lái)自中國(guó)科學(xué)院科學(xué)數(shù)據(jù)庫(kù)分辨率為30 m×30 m的數(shù)字高程地圖(DEM)數(shù)據(jù)計(jì)算得到;2)植被數(shù)據(jù):依據(jù)landsat TM遙感圖片,利用全國(guó)森林二類調(diào)查數(shù)據(jù)及全國(guó)第四次大熊貓調(diào)查植被樣方數(shù)據(jù)對(duì)其進(jìn)行監(jiān)督分類和校正;3)水體數(shù)據(jù):包括到大河流的距離和到小河的流距離;4)大熊貓可食竹分布數(shù)據(jù):來(lái)源于全國(guó)第四次大熊貓調(diào)查;5)人類干擾數(shù)據(jù):選取到居民用地的距離、到主要公路的距離、到耕地的距離、居民點(diǎn)密度等指標(biāo)。

    2.3數(shù)據(jù)分析

    2.3.1評(píng)價(jià)棲息地質(zhì)量與劃定棲息地斑塊

    本研究利用大熊貓分布點(diǎn)和棲息地生態(tài)環(huán)境因子,通過(guò)MAXENT 3.3軟件包進(jìn)行MAXENT模型分析,采用交叉驗(yàn)證(cross validation)法重復(fù)計(jì)算10次,以10次計(jì)算結(jié)果的均值作為HSI。約登指數(shù)(Yourdon's index)是靈敏度和特異度之和減去1,指數(shù)越大,其真實(shí)性越高,故選用最大約登指數(shù)為閾值判定核心棲息地范圍[27]。

    2.3.2構(gòu)建阻力圖層

    本研究基于棲息地適宜度指數(shù),并結(jié)合專家意見[28],以1—1000作為阻力值構(gòu)建阻力圖層。根據(jù)第四次大熊貓調(diào)查數(shù)據(jù),四川省各山系99%大熊貓痕跡點(diǎn)分布于1600—3800 m,因此我們認(rèn)為此區(qū)間以外的海拔范圍大熊貓通過(guò)概率低,可以排除在廊道規(guī)劃之外。研究區(qū)域存在的大型水庫(kù),大熊貓難以直接通過(guò),也排除在廊道規(guī)劃之外。此外,雖然大熊貓可能通過(guò)農(nóng)田區(qū)域,但是從保護(hù)管理角度看,難以在農(nóng)田區(qū)域開展廊道恢復(fù),因此農(nóng)田區(qū)域也不納入廊道規(guī)劃。

    2.3.3規(guī)劃廊道

    基于構(gòu)建的阻力圖層和棲息地斑塊,采用野生動(dòng)物棲息地連接性分析工具箱Linkage Mapper 1.0.8[29]進(jìn)行廊道設(shè)計(jì)與規(guī)劃。首先通過(guò)LCP法分析確定了最小代價(jià)路徑,得到路徑代價(jià)柵格圖,依據(jù)大熊貓最小家域3.9 km2[22],設(shè)置廊道的最小寬度為2300 m,優(yōu)先選擇代價(jià)距離最小的區(qū)域作為廊道的范圍。

    在此基礎(chǔ)上,通過(guò)障礙分析,計(jì)算改善指數(shù)識(shí)別棲息地亟需得到優(yōu)先恢復(fù)的關(guān)鍵區(qū)域。改善指數(shù)高意味著該區(qū)域的棲息地若得到恢復(fù),可以最大程度地提高個(gè)體在所連接斑塊之間的移動(dòng)可能[30]。最后,基于電流理論模型進(jìn)行流量分析計(jì)算大熊貓穿過(guò)廊道的流通密度。對(duì)廊道中大熊貓可能通過(guò)的熱點(diǎn)區(qū)域進(jìn)行識(shí)別,評(píng)估廊道的效用。

    3 研究結(jié)果

    通過(guò)MAXENT模型進(jìn)行棲息地評(píng)估,得到HSI。受試者工作特征曲線(ROC曲線)分析表明,訓(xùn)練集和驗(yàn)證集的AUC(AUC)值分別為0.9939和0.9890,表明模型預(yù)測(cè)結(jié)果優(yōu)秀[31-32]。

    基于HSI,以最大約登指數(shù)(0.1961)為閾值判定棲息地斑塊。小相嶺大熊貓棲息地可以分為7個(gè)斑塊,即麻麻地、紫馬、魯壩、高橋、石灰窯、大火地、灣壩(圖1)。其中麻麻地和紫馬的大熊貓痕跡點(diǎn)分布較為密集,而魯壩沒有記錄到大熊貓活動(dòng)痕跡。

    圖1 小相嶺山系大熊貓廊道最小代價(jià)路徑分析Fig.1Least cost path analysis of giant pandas corridors in the Xiaoxiangling Mountains

    3.1廊道范圍

    通過(guò)LCP分析在整個(gè)山系總共規(guī)劃了9條廊道,即麻麻地—紫馬(以下簡(jiǎn)稱:麻紫)、麻麻地—魯壩(麻魯)、紫馬—魯壩(紫魯)、紫馬—高橋(紫高)、紫馬—石灰窯(紫石)、石灰窯—高橋(石高)、石灰窯—大火地(石大)、石灰窯—魯壩(石魯)、大火地—魯壩(大魯)(圖1)。其中,紫高廊道、紫石廊道和石高廊道彼此鄰近,麻紫廊道和紫魯廊道也極其鄰近,此外,石大廊道、石魯廊道和大魯廊道相互毗鄰甚至相交。

    LCP分析顯示,有兩片適宜區(qū)域可以規(guī)劃為連接麻麻地和紫馬斑塊的廊道。其中,北側(cè)區(qū)域(圖1中B廊道)穿過(guò)108國(guó)道處農(nóng)田較多(圖2),嚴(yán)重?cái)D壓廊道。鑒于大面積改造農(nóng)田難度較大,基于管理考慮,將該區(qū)域作為潛在廊道,而選用南側(cè)區(qū)域(圖1中A廊道)作為麻紫廊道。

    麻紫廊道和麻魯廊道主要連接了108國(guó)道東西兩側(cè)的大熊貓棲息地(圖1),其余7條是西側(cè)各大熊貓局部種群間交流的重要通道。受北側(cè)的高速公路和周邊農(nóng)田擠壓(圖2),連接麻麻地和魯壩的適宜區(qū)域形成一個(gè)瓶頸區(qū)域,最狹窄處只有大約40 m,除非該區(qū)域環(huán)境得到改善,否則麻魯廊道很難有效發(fā)揮其廊道作用。

    彎壩斑塊與其他斑塊隔離嚴(yán)重,僅通過(guò)紫馬斑塊連接到其他斑塊,與其余幾個(gè)斑塊間的間接通道被大范圍海拔4000 m以上的山脊阻隔(圖2)。并且彎壩-紫馬斑連接通道的距離很長(zhǎng),最小代價(jià)路徑大約為36 km(表1),可能不適于大熊貓通過(guò),因此也將其規(guī)劃為潛在廊道,不作進(jìn)一步分析。

    圖2 小相嶺山系廊道周圍農(nóng)田分布圖Fig.2The farmland around the giant pandas corridors in the Xiaoxiangling Mountains

    表1 小相嶺山系大熊貓廊道面積以及最小代價(jià)路徑長(zhǎng)度Table1 The area and the least cost path length of giant pandas corridors in the Xiaoxiangling Mountains

    石高廊道和紫石廊道有較大范圍的重疊。如果將其合并可以達(dá)到同時(shí)連接紫馬、高橋和石灰窯斑塊的效果。石大廊道和石魯廊道,石魯廊道和大魯廊道局部也存在一定的重疊現(xiàn)象,合并也可以同時(shí)連接石灰窯、大火地和魯壩斑塊。

    3.2障礙分析

    障礙分析結(jié)果確定了9條廊道內(nèi)部亟需優(yōu)先恢復(fù)的關(guān)鍵區(qū)域(圖3)。紫高廊道的優(yōu)先恢復(fù)區(qū)域主要集中在廊道的西北區(qū)域,石高和紫石廊道的優(yōu)先恢復(fù)區(qū)域,則分別集中在廊道的西南和東北區(qū)域,而麻紫等其它6條廊道的優(yōu)先恢復(fù)區(qū)域主要集中在廊道中部。若這些區(qū)域的棲息地得到了恢復(fù),則能使整個(gè)廊道的阻力得到最大程度的降低,極大的促進(jìn)棲息地斑塊間的連接性,增進(jìn)大熊貓種群的交流。

    3.3流量分析

    流量分析結(jié)果顯示,小相嶺山系大熊貓通過(guò)廊道的密集區(qū)域主要集中在廊道的狹窄區(qū)域(圖4)。由于受到人類干擾和地形條件的限制,大熊貓只能選擇這些狹窄區(qū)域通過(guò),若能改善這些區(qū)域周邊的棲息地質(zhì)量,則能擴(kuò)寬廊道范圍,降低大熊貓通行的阻力。其中,紫石廊道和高石廊道的西側(cè)受到高海拔的擠壓,形成瓶頸區(qū)域。而紫石廊道和石魯廊道由于受到水庫(kù)的影響,同樣在廊道內(nèi)部形成了瓶頸區(qū)域。此外,在農(nóng)田侵蝕下,麻魯廊道也出現(xiàn)了瓶頸區(qū),這也顯示如果能將該區(qū)域的農(nóng)田恢復(fù)自然狀態(tài),將極大地促進(jìn)大熊貓棲息地斑塊之間的連通。

    圖3 小相嶺山系大熊貓廊道障礙分析Fig.3Barrier detection analysis of giant pandas corridors in the Xiaoxiangling Mountains

    圖4 小相嶺山系大熊貓廊道流量分析Fig.4Current analysis of giant pandas corridors in the Xiaoxiangling Mountains

    3.4廊道區(qū)域內(nèi)土地覆被類型

    廊道區(qū)域內(nèi)的土地利用與土地覆被類型包括常綠闊葉林、落葉闊葉林、暖性針葉林、溫性針葉林、寒溫性針葉林、溫性針闊葉混交林、竹林、常綠灌叢、落葉灌叢、灌草叢、草甸、栽培森林植被、水域、公路及荒廢的建設(shè)用地(表2)。其中廊道范圍內(nèi)通過(guò)的公路主要有高速公路G5(雅安—西昌段)、國(guó)道108(石棉—冕寧段)以及栗子坪到冶勒的鄉(xiāng)村道。

    表2 小相嶺山系大熊貓廊道土地利用與土地覆被類型Table2 Land use and vegetation coverage of giant pandas corridors in the Xiaoxiangling Mountains

    9條廊道的總面積為11835.58 hm2,與大熊貓生境密切相關(guān)的針葉林、闊葉林(本研究區(qū)域包括常綠闊葉林類的青岡林和落葉闊葉林)和針闊混交林[33-34]廣泛分布于廊道內(nèi),總面積達(dá)6312.69 hm2,占廊道總面積的53.34%。其中紫魯廊道范圍內(nèi)這幾種植被類型所占面積比例最大,達(dá)85.87%,其次為麻紫和紫高廊道,比例分別為69.32%、68.88%。廊道范圍內(nèi)的竹林基本上全為箭竹林,其中分布比例最大的廊道是紫石廊道,占廊道面積的28.00%。

    4 討論

    野生動(dòng)物廊道有助于被斑塊化棲息地分隔的動(dòng)物種群完成遷移、擴(kuò)散、繁殖等關(guān)鍵生活史過(guò)程,增加種群間的基因和物種交流,降低種群滅絕風(fēng)險(xiǎn)[35]。本文基于MAXENT模型計(jì)算的棲息地適宜度指數(shù)構(gòu)建阻力圖層,結(jié)合LCP分析法、障礙分析法、流量分析法綜合規(guī)劃廊道,得到較為全面的分析結(jié)果。

    目前,關(guān)于大熊貓廊道的研究較少。Li等[36]基于最小代價(jià)路徑模型,分析了汶川大地震后臥龍自然保護(hù)區(qū)大熊貓種群分布格局及其之間的連接通道;而Xu等[37]提出了邛崍山系大熊貓保護(hù)方案,認(rèn)為該區(qū)域亟需建設(shè)大熊貓廊道。此外,Qi等[38]也曾用LCP法對(duì)小相嶺山系大熊貓棲息地進(jìn)行了線性廊道規(guī)劃,而本研究則進(jìn)一步依據(jù)代價(jià)距離梯度規(guī)劃了帶狀的廊道區(qū)域(拉伸的面狀景觀要素)。與沒有寬度的線性廊道相比,帶狀的廊道生態(tài)價(jià)值更大,在保護(hù)管理工作中也更容易操作。

    此外,本文還提出了網(wǎng)狀的廊道群規(guī)劃。研究結(jié)果顯示部分廊道彼此鄰近甚至有一定程度的重疊,形成了“廊道群”。廊道群是兩個(gè)或兩個(gè)以上極其相鄰的廊道構(gòu)成的網(wǎng)狀結(jié)構(gòu),這樣的網(wǎng)狀廊道結(jié)構(gòu),不僅能夠連接兩個(gè)相鄰的棲息地斑塊,而且能夠促進(jìn)多個(gè)棲息地斑塊間大熊貓種群的交流,實(shí)現(xiàn)對(duì)種群和棲息地更有效的保護(hù)。本文中紫高廊道、紫石廊道和石高廊道彼此鄰近,形成了廊道群;石大廊道、石魯廊道和大魯廊道相互毗鄰甚至相交,同時(shí)連接石灰窯斑塊、大火地斑塊和魯壩斑塊的大熊貓棲息地,實(shí)現(xiàn)了棲息地斑塊之間的兩兩互連。

    Xu等[37]認(rèn)為廊道長(zhǎng)度應(yīng)該小于大熊貓家域范圍(約2 km),否則會(huì)影響大熊貓遷移的效率。但是潘文石等[39]利用無(wú)線電跟蹤兩個(gè)亞成體大熊貓的遷移時(shí)發(fā)現(xiàn),它們分別移動(dòng)了34 km和24 km。此外,從種群遺傳學(xué)角度看,即使很少量個(gè)體的交流,也可以極大地增加被分隔的局部種群的遺傳多樣性,對(duì)種群穩(wěn)定有著非常重要的意義。因此,本文沒有采用這一限制條件。

    麻麻地—紫馬廊道連接了兩塊面積最大,同時(shí)也是大熊貓痕跡點(diǎn)最多的兩個(gè)棲息地斑塊:麻麻地和紫馬斑塊。這兩個(gè)棲息地斑塊對(duì)整個(gè)小相嶺區(qū)域大熊貓種群穩(wěn)定有著極為重要的意義。但是這兩個(gè)斑塊間有高速公路G5(雅安—西昌段)、國(guó)道108(石棉—冕寧段)以及栗子坪到冶勒的鄉(xiāng)村道通過(guò),溝谷地帶農(nóng)田較多,人類活動(dòng)頻繁,植被也不利于大熊貓遷移。因此,建議優(yōu)先加強(qiáng)這一廊道的建設(shè)。

    灣壩斑塊位于九龍縣,野外調(diào)查記錄到極少大熊貓活動(dòng)痕跡。由于大范圍的高海拔山脊,與鄰近的高橋、石灰窯、大火地等斑塊形成隔離,即使大熊貓偶爾能夠跨越這些4000 m以上的山脊,但兩側(cè)的種群交流也必然受到很大影響。若避開高海拔區(qū)域,彎壩斑塊只可能與紫馬斑塊連通,但是該通道長(zhǎng)達(dá)36.19 km,大熊貓利用可能性也很小。因此,灣壩斑塊實(shí)際上是一個(gè)孤島,建議加強(qiáng)該區(qū)域種群的保護(hù)和管理。

    全國(guó)第四次大熊貓調(diào)查并未在魯壩斑塊發(fā)現(xiàn)大熊貓活動(dòng)痕跡,但該區(qū)域植被較好,竹類分布豐富,通過(guò)廊道連通,該區(qū)域有望恢復(fù)成為大熊貓棲息地。此外,從空間布局上看,魯壩斑塊的存在可以極大地改善麻麻地、紫馬與其他幾個(gè)斑塊的大熊貓種群連接性,對(duì)當(dāng)?shù)卮笮茇埛N群的穩(wěn)定有重要價(jià)值。因此,建議通過(guò)廊道建設(shè)和進(jìn)一步的棲息地保護(hù)恢復(fù)該區(qū)域的種群。

    本研究規(guī)劃的廊道區(qū)域周邊有大量農(nóng)田分布,嚴(yán)重?cái)D壓大熊貓廊道,形成狹窄的瓶頸。如果能對(duì)農(nóng)田區(qū)域進(jìn)行退耕恢復(fù),則能有效地提高小相嶺山系大熊貓種群的連通性。但是農(nóng)田區(qū)域的退耕恢復(fù)的管理成本很高,這也極大地制約了小相嶺山系大熊貓棲息地的恢復(fù)工作。因此,在今后的工作中,需要從保護(hù)管理角度探索農(nóng)田區(qū)域改造的方式,以促進(jìn)小相嶺山系大熊貓種群和棲息地更有效地保護(hù)。

    參考文獻(xiàn)(References):

    [1]程宏毅,鮑毅新,葛寶明,鄭祥,胡一中.棲息地片斷化對(duì)動(dòng)物種群間基因流的影響及其測(cè)定方法.生態(tài)學(xué)雜志,2006,25(7):863-868.

    [2]Crooks K R,Sanjayan M.Connectivity Conservation.Cambridge:Cambridge University Press,2006.

    [3]Willis E O.Populations and local extinctions of birds on barro colorado island,panama.Ecological Monographs,1974,44(2):153-169.

    [4]Spellerberg I F,Gaywood M J.Linear features:linear habitats and wildlife corridors.Peterborough:English Nature,1993.

    [5]Brookes C J.A parameterized region-growing programme for site allocation on raster suitability maps.International Journal of Geographical Information Science,1997,11(4):375-396.

    [6]Wilson E O,Willis E O.Applied biogeography//Cody M L,Diamond J M,eds.Ecology and Evolution of Communities.Cambridge:Harvard University Press,1975:522-534.

    [7]武正軍,李義明.生境破碎化對(duì)動(dòng)物種群存活的影響.生態(tài)學(xué)報(bào),2003,23(11):2424-2435.

    [8]Sawyer S C,Epps C W,Brashares J S.Placing linkages among fragmented habitats:do least-cost models reflect how animals use landscapes?.Journal of Applied Ecology,2011,48(3):668-678.

    [9]Knaapen J P,Scheffer M,Harms B.Estimating habitat isolation in landscape planning.Landscape and Urban Planning,1992,23(1):1-16.

    [10]李紀(jì)宏,劉雪華.基于最小費(fèi)用距離模型的自然保護(hù)區(qū)功能分區(qū).自然資源學(xué)報(bào),2006,21(2):217-224.

    [11]Mcrae B H,Dickson B G,Keitt T H,Shah V B.Using circuit theory to model connectivity in ecology,evolution,and conservation.Ecology,2008,89(10):2712-2724.

    [12]Tracey J A.Individual-based modeling as a tool for conserving connectivity//Crooks K R,Sanjayan M,eds.Connectivity Conservation.Cambridge:Cambridge University Press,2006:343-368.

    [13]Zhu L F,Zhan X J,Wu H,Zhang S N,Meng T,Bruford M W,Wei F W.Conservation implications of drastic reductions in the smallest and most isolated populations of giant pandas.Conservation Biology,2010,24(5):1299-1306.

    [14]朱磊,吳攀文,張洪,胡錦矗.小相嶺山系大熊貓種群生存力分析.西華師范大學(xué)學(xué)報(bào):自然科學(xué)版,2008,29(2):112-116.

    [15]馮文和,李光漢.拯救大熊貓.成都:四川科學(xué)技術(shù)出版社,2000.

    [16]Jaynes E T.Information theory and statistical mechanics.Physical Review,1957,106(4):620-630.

    [17]Phillips S J,Anderson R P,Schapire R E.Maximum entropy modeling of species geographic distributions.Ecological Modelling,2006,190(3/ 4):231-259.

    [18]Hernandez P A,F(xiàn)ranke I,Herzog S K,Pacheco V,Paniagua L,Quintana H L,Soto A,Swenson J J,Tovar C,Valqui T H,Vargas J,Young B E.Predicting species distributions in poorly-studied landscapes.Biodiversity and Conservation,2008,17(6):1353-1366.

    [19]Zeller K A,McGarigal K,Whiteley A R.Estimating landscape resistance to movement:a review.Landscape Ecology,2012,27(6):777-797.

    [20]Spear S F,Balkenhol N,F(xiàn)ortin M J,McRae B H,Scribner K.Use of resistance surfaces for landscape genetic studies:considerations for parameterization and analysis.Molecular Ecology,2010,19(17):3576-3591.

    [21]Adriaensen F,Chardon J P,De Blust G,Swinnen E,Villalba S,Gulinck H,Matthysen E.The application of‘least-cost’modelling as a functional landscape model.Landscape and Urban Planning,2003,64(4):233-247.

    [22]胡錦矗,夏勒,潘文石,朱靖.臥龍的大熊貓.成都:四川科學(xué)技術(shù)出版社,1985.

    [23]魏輔文,馮祚建,王祖望.相嶺山系大熊貓和小熊貓對(duì)生境的選擇.動(dòng)物學(xué)報(bào),1999,45(1):57-71.

    [24]冉江洪,劉少英,王鴻加,孫治宇,曾宗永,劉世昌.放牧對(duì)冶勒自然保護(hù)區(qū)大熊貓生境的影響.獸類學(xué)報(bào),2003,23(4):288-294.

    [25]張澤鈞,胡錦矗,吳華.邛崍山系大熊貓和小熊貓生境選擇的比較.獸類學(xué)報(bào),2002,22(3):161-168.

    [26]王學(xué)志,徐衛(wèi)華,歐陽(yáng)志云,劉建國(guó),肖燚,陳佑平,趙聯(lián)軍,黃俊忠.生態(tài)位因子分析在大熊貓(Ailuropoda melanoleuca)生境評(píng)價(jià)中的應(yīng)用.生態(tài)學(xué)報(bào),2008,28(2):821-828.

    [27]Jiménez-Valverde A,Lobo J M.Threshold criteria for conversion of probability of species presence to either-or presence-absence.Acta Oecologica,2007,31(3):361-369.

    [28]Beier P,Majka D R,Spencer W D.Forks in the road:choices in procedures for designing wildland linkages.Conservation Biology,2008,22(4): 836-851.

    [29]Carroll C,Mcrae B H,Brookes A.Use of linkage mapping and centrality analysis across habitat gradients to conserve connectivity of gray wolf populations in Western North America.Conservation Biology,2012,26(1):78-87.

    [30]McRae B H,Hall S A,Beier P,Theobald D M.Where to restore ecological connectivity?Detecting barriers and quantifying restoration benefits.PLoS One,2012,7(12):e52604.

    [31]Swets J A.Measuring the accuracy of diagnostic systems.Science,1988,240(4857):1285-1293.

    [32]Araújo M B,Pearson R G,Thuiller W,Erhard M.Validation of species–climate impact models under climate change.Global Change Biology,2005,11(9):1504-1513.

    [33]魏輔文,周昂,胡錦矗,王維,楊光.馬邊大風(fēng)頂自然保護(hù)區(qū)大熊貓對(duì)生境的選擇.獸類學(xué)報(bào),1996,16(4):241-245.

    [34]徐衛(wèi)華,歐陽(yáng)志云,蔣澤銀,鄭華,劉建國(guó).大相嶺山系大熊貓生境評(píng)價(jià)與保護(hù)對(duì)策研究.生物多樣性,2006,14(3):223-231.

    [35]呂海燕,李政海,李建東,宋國(guó)寶.廊道研究進(jìn)展與主要研究方法.安徽農(nóng)業(yè)科學(xué),2007,35(15):4484.

    [36]Li H L,Li D H,Li T,Qiao Q,Yang J,Zhang H M.Application of least-cost path model to identify a giant panda dispersal corridor network after the Wenchuan earthquake-case study of Wolong Nature Reserve in China.Ecological Modelling,2010,221(6):944-952.

    [37]Xu W H,Ouyang Z Y,Vi?a A,Zheng H,Liu J G,Xiao Y.Designing a conservation plan for protecting the habitat for giant pandas in the Qionglai Mountain range,China.Diversity and Distributions,2006,12(5):610-619.

    [38]Qi D W,Hu Y B,Gu X D,Yang X Y,Yang G,Wei F W.Quantifying landscape linkages among giant panda subpopulations in regional scale conservation.Integrative Zoology,2012,7(2):165-174.

    [39]潘文石,呂植,朱小健,王大軍,王昊,龍玉,付達(dá)莉,周欣.繼續(xù)生存的機(jī)會(huì).北京:北京大學(xué)出版社,2001.

    Corridor design for the giant panda in the Xiaoxiangling Mountains

    QING Jing1,2,XU Chi3,YANG Biao4,YANG Zhisong1,*,QI Dunwu5,YANG Xuyu6,GU Xiaodong6,DAI Qiang2
    1 Giant Panda Conservation and Culture Research Center,China West Normal University,Nanchong 637002,China
    2 Chengdu Institute of Biology,the Chinese Academy of Sciences,Chengdu 610041,China
    3 Sichuan Xue Baoding Nature Reserve,Mianyang 622550,China
    4 Key Laboratory of Bio-Resources and Eco-Environment of Ministry Education,College of Life Science,Sichuan University,Chengdu 610064,China
    5 Sichuan Key Laboratory of Conservation Biology for Endangered Wildlife,Chengdu Research Base of Giant Panda Breeding,Chengdu 610086,China
    6 Sichuan Station of Wild life survey and Management,Chengdu 610081,China

    Abstract:The giant panda(Ailuropoda melanoleuca)is regarded as one of the most imperiled mammals in the world.The quality and quantity of giant panda habitat have been continually and dramatically degraded due to increasing human activity and to natural disasters.In order to conserve the remaining panda populations effectively,there is a critical need to establish or restore corridors between remaining islands of habitat.This is considered essential for reducing the risk of local populationbook=1126,ebook=238extinctions.The population of giant pandas in the Xiaoxiangling Mountains is the smallest,and the habitat there is seriously fragmented.To connect the existing habitat patches and protect its local population of giant pandas,we used quantitative habitat assessment methods to assess the current structure of the Xiaoxiangling Mountains habitat and to design wildlife corridors to increase connectivity.The habitat suitability index in the study area was assessed using Maximum Entropy analysis,with“area under curve”of 0.9939 and 0.9890 for training and test sets,respectively.The habitat patches were identified following the threshold defined by the maximum Yourdon's index.There are seven isolated habitat patches in the Xiaoxiangling Mountains(i.e.,Zima,Mamadi,Gaoqiao,Shihuiyao,Luba,Dahuodi,and Wanba).Based on the leastcost path and circuit theory models,nine corridors between those habitat patches were designed.Among these corridors,the Zima-Mamadi corridor connects the two largest habitat patches,which are also where most of the giant panda activities have been recorded.Given its ecological significance,we propose that construction of Zima-Mamadi corridor should be considered top priority.In addition,some of the corridors adjoined or even overlapped.Therefore,we suggest merging these corridors into a“corridor group”.A“corridor group”is defined as a network constituted by two or more closely adjacent corridors,which could interconnect several habitat patches simultaneously and effectively increase the exchange between populations in those habitat patches.Unfortunately,because of the presence of farmlands,some sections of the corridors formed extremely narrow channels.Some of these channels are narrower than the minimum efficient width of giant panda corridor(2230 m),and so their intended function could be severely impacted.However,given the complexities of land ownership,determining a way to construct corridors among farmlands should be a focus of future conservation work.

    Key Words:giant panda;xiaoxiangling;least-cost path;circuit theory;corridor group

    *通訊作者

    Corresponding author.E-mail:yangzhisong@126.com

    收稿日期:2014-06-13;網(wǎng)絡(luò)出版日期:2015-07-09

    基金項(xiàng)目:國(guó)家林業(yè)局“邛崍山與大相嶺山系連接地帶構(gòu)成大熊貓走廊帶關(guān)鍵因素研究”(林護(hù)發(fā)(2011)285號(hào));四川省教育廳重點(diǎn)項(xiàng)目“泥巴山大熊貓走廊帶關(guān)鍵構(gòu)成因素”(10ZA027);香港海洋公園保育基金支持項(xiàng)目“小相嶺大熊貓棲息地走廊帶規(guī)劃與保護(hù)管理策略制定”(GP061213)

    DOI:10.5846/stxb201406131225

    少妇的丰满在线观看| 久久九九热精品免费| 日韩三级视频一区二区三区| 中文亚洲av片在线观看爽| 国产成人av激情在线播放| 久久午夜综合久久蜜桃| 高清在线国产一区| 99久久精品国产亚洲精品| 亚洲人成伊人成综合网2020| 青草久久国产| 精品熟女少妇八av免费久了| 最近在线观看免费完整版| 国产精品美女特级片免费视频播放器 | 欧美日韩瑟瑟在线播放| 亚洲男人的天堂狠狠| 国产精品 国内视频| 999精品在线视频| 脱女人内裤的视频| 91成年电影在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲第一电影网av| 免费在线观看成人毛片| 国产欧美日韩一区二区精品| 国产三级黄色录像| 国产成人系列免费观看| 天堂√8在线中文| 亚洲人成77777在线视频| 国产精品日韩av在线免费观看| 国产午夜精品久久久久久| 啦啦啦 在线观看视频| 深夜精品福利| 国产欧美日韩精品亚洲av| 国产精品久久久久久精品电影 | 51午夜福利影视在线观看| 国产免费av片在线观看野外av| 精品国产超薄肉色丝袜足j| 人人澡人人妻人| 狂野欧美激情性xxxx| 亚洲av五月六月丁香网| 精品久久久久久久久久久久久 | 国产精品亚洲av一区麻豆| 制服人妻中文乱码| www.熟女人妻精品国产| 可以在线观看毛片的网站| 久久精品亚洲精品国产色婷小说| 午夜两性在线视频| 人人妻人人澡欧美一区二区| 美女大奶头视频| 长腿黑丝高跟| 视频区欧美日本亚洲| 午夜亚洲福利在线播放| 一本一本综合久久| 亚洲成国产人片在线观看| 午夜激情福利司机影院| 一区二区三区精品91| 少妇裸体淫交视频免费看高清 | 窝窝影院91人妻| 99久久综合精品五月天人人| 在线观看免费日韩欧美大片| 中国美女看黄片| 日韩欧美免费精品| 99re在线观看精品视频| 亚洲熟女毛片儿| 老熟妇乱子伦视频在线观看| 欧美乱码精品一区二区三区| 国产精品一区二区免费欧美| 欧美成人免费av一区二区三区| 一进一出抽搐动态| 18禁裸乳无遮挡免费网站照片 | 十八禁人妻一区二区| 久久中文字幕一级| av有码第一页| 国产欧美日韩一区二区三| 一级a爱片免费观看的视频| 欧美又色又爽又黄视频| 久久久久久九九精品二区国产 | 搡老熟女国产l中国老女人| 日日摸夜夜添夜夜添小说| 观看免费一级毛片| 精品久久久久久久久久久久久 | 亚洲 欧美一区二区三区| 日韩大尺度精品在线看网址| √禁漫天堂资源中文www| 法律面前人人平等表现在哪些方面| 在线观看免费午夜福利视频| 一本一本综合久久| 男女做爰动态图高潮gif福利片| 日韩 欧美 亚洲 中文字幕| 美女高潮喷水抽搐中文字幕| 久久久久久久精品吃奶| 国产成人欧美在线观看| 成人18禁在线播放| 国产成+人综合+亚洲专区| 国产黄片美女视频| 18禁观看日本| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看视频国产中文字幕亚洲| 变态另类成人亚洲欧美熟女| 一区二区三区高清视频在线| 欧美在线黄色| 美女午夜性视频免费| 在线观看一区二区三区| 亚洲av美国av| 久久人妻福利社区极品人妻图片| 女性生殖器流出的白浆| 人成视频在线观看免费观看| 国产精品久久视频播放| 后天国语完整版免费观看| 香蕉国产在线看| 亚洲成a人片在线一区二区| 首页视频小说图片口味搜索| 欧美日韩亚洲国产一区二区在线观看| 可以在线观看的亚洲视频| 18禁黄网站禁片午夜丰满| 国产麻豆成人av免费视频| 亚洲欧美激情综合另类| 免费在线观看亚洲国产| 午夜福利成人在线免费观看| 最新在线观看一区二区三区| 国产日本99.免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕精品亚洲无线码一区 | 69av精品久久久久久| 国产精品精品国产色婷婷| 欧美激情极品国产一区二区三区| av有码第一页| 性色av乱码一区二区三区2| 亚洲精品中文字幕一二三四区| 免费观看人在逋| 一边摸一边抽搐一进一小说| 亚洲五月婷婷丁香| 好男人在线观看高清免费视频 | 女人爽到高潮嗷嗷叫在线视频| 黄色a级毛片大全视频| 波多野结衣巨乳人妻| 波多野结衣巨乳人妻| 亚洲国产日韩欧美精品在线观看 | 美国免费a级毛片| 欧美乱色亚洲激情| 国产激情偷乱视频一区二区| 制服人妻中文乱码| 成人亚洲精品av一区二区| 欧美另类亚洲清纯唯美| 亚洲欧美精品综合久久99| 国产精品亚洲av一区麻豆| 中文字幕精品亚洲无线码一区 | 十分钟在线观看高清视频www| 亚洲aⅴ乱码一区二区在线播放 | 在线免费观看的www视频| 亚洲精品美女久久av网站| 亚洲一码二码三码区别大吗| 国产极品粉嫩免费观看在线| 午夜成年电影在线免费观看| 校园春色视频在线观看| 国产激情久久老熟女| 精品国内亚洲2022精品成人| av电影中文网址| netflix在线观看网站| 亚洲国产日韩欧美精品在线观看 | 变态另类成人亚洲欧美熟女| 成人特级黄色片久久久久久久| 丁香欧美五月| 久久精品国产清高在天天线| 亚洲无线在线观看| 国产精品久久视频播放| 国产欧美日韩一区二区精品| 狂野欧美激情性xxxx| 色综合婷婷激情| 国产乱人伦免费视频| 亚洲精品国产区一区二| 麻豆成人av在线观看| 日韩精品中文字幕看吧| 成人精品一区二区免费| 久久 成人 亚洲| 日日夜夜操网爽| 国产成+人综合+亚洲专区| 又紧又爽又黄一区二区| 人人妻人人澡欧美一区二区| 免费av毛片视频| 天天添夜夜摸| 一边摸一边抽搐一进一小说| 99久久无色码亚洲精品果冻| 两个人视频免费观看高清| 成人午夜高清在线视频 | 亚洲色图av天堂| 中文字幕人成人乱码亚洲影| 黄色视频,在线免费观看| 亚洲国产精品成人综合色| 两人在一起打扑克的视频| 最近最新中文字幕大全电影3 | 婷婷精品国产亚洲av在线| 日韩欧美一区视频在线观看| 国内毛片毛片毛片毛片毛片| 1024手机看黄色片| 欧美中文综合在线视频| 级片在线观看| 脱女人内裤的视频| 老熟妇乱子伦视频在线观看| 久久精品91蜜桃| 欧美av亚洲av综合av国产av| 丰满的人妻完整版| 露出奶头的视频| 丝袜人妻中文字幕| 丝袜人妻中文字幕| 精品无人区乱码1区二区| 一本精品99久久精品77| 欧美成人午夜精品| 美女大奶头视频| 久久中文字幕一级| 麻豆成人午夜福利视频| 日本三级黄在线观看| 黑人操中国人逼视频| 亚洲精品国产区一区二| 国产主播在线观看一区二区| 淫秽高清视频在线观看| 久久九九热精品免费| 中出人妻视频一区二区| 成人精品一区二区免费| 老汉色∧v一级毛片| 一进一出好大好爽视频| 色精品久久人妻99蜜桃| 亚洲成av人片免费观看| 亚洲男人的天堂狠狠| 久久久久久久久中文| 又紧又爽又黄一区二区| 国内揄拍国产精品人妻在线 | 中文字幕精品亚洲无线码一区 | 老司机在亚洲福利影院| 不卡一级毛片| 国产免费av片在线观看野外av| 欧美在线一区亚洲| 亚洲欧美激情综合另类| √禁漫天堂资源中文www| 亚洲国产日韩欧美精品在线观看 | 欧美不卡视频在线免费观看 | 777久久人妻少妇嫩草av网站| 午夜福利在线在线| АⅤ资源中文在线天堂| 亚洲国产精品成人综合色| 在线观看免费日韩欧美大片| 免费在线观看成人毛片| 久久人妻av系列| 免费无遮挡裸体视频| ponron亚洲| 男女床上黄色一级片免费看| 久久午夜亚洲精品久久| 啦啦啦观看免费观看视频高清| 精品国产亚洲在线| 久久久久久久久久黄片| 一区二区三区精品91| 欧洲精品卡2卡3卡4卡5卡区| 免费观看精品视频网站| 免费无遮挡裸体视频| 亚洲国产精品999在线| 久久香蕉国产精品| 国产亚洲精品久久久久5区| 两人在一起打扑克的视频| 国产日本99.免费观看| 最近最新中文字幕大全电影3 | 18禁观看日本| 黄色女人牲交| 91av网站免费观看| 十八禁人妻一区二区| 国产区一区二久久| 丁香六月欧美| netflix在线观看网站| 亚洲一区二区三区色噜噜| 国产精品综合久久久久久久免费| 男女下面进入的视频免费午夜 | 日韩三级视频一区二区三区| 90打野战视频偷拍视频| 欧美黑人巨大hd| 亚洲五月婷婷丁香| 国产av一区在线观看免费| 中文亚洲av片在线观看爽| 91成人精品电影| 亚洲成人精品中文字幕电影| 国产精品自产拍在线观看55亚洲| 免费在线观看黄色视频的| 人妻丰满熟妇av一区二区三区| 精品国产亚洲在线| 亚洲男人天堂网一区| www.www免费av| 天堂影院成人在线观看| 亚洲欧美精品综合一区二区三区| 国产高清激情床上av| 日本免费a在线| 很黄的视频免费| 色播在线永久视频| 欧美 亚洲 国产 日韩一| 伊人久久大香线蕉亚洲五| 午夜福利免费观看在线| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 国产一级毛片七仙女欲春2 | 久久久精品国产亚洲av高清涩受| 男男h啪啪无遮挡| 999精品在线视频| 国产黄a三级三级三级人| 一本一本综合久久| 中文字幕av电影在线播放| 日韩精品青青久久久久久| 俄罗斯特黄特色一大片| 国产精品亚洲av一区麻豆| 久久精品国产综合久久久| 亚洲中文字幕一区二区三区有码在线看 | 国产99白浆流出| 熟女电影av网| 久久精品国产清高在天天线| 久久精品国产亚洲av高清一级| 国产成人av教育| 天堂影院成人在线观看| 久久久久国产精品人妻aⅴ院| 91字幕亚洲| 国产91精品成人一区二区三区| 亚洲五月天丁香| 欧美成人免费av一区二区三区| 免费在线观看视频国产中文字幕亚洲| 嫩草影视91久久| 视频区欧美日本亚洲| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩高清在线视频| 少妇 在线观看| 精品乱码久久久久久99久播| 男人的好看免费观看在线视频 | 美女 人体艺术 gogo| 久久久久免费精品人妻一区二区 | 欧美黑人欧美精品刺激| 国产成人av教育| 美女国产高潮福利片在线看| 丰满的人妻完整版| 国产麻豆成人av免费视频| 搞女人的毛片| 国产av一区二区精品久久| 在线观看日韩欧美| 天天一区二区日本电影三级| 岛国视频午夜一区免费看| 老司机午夜十八禁免费视频| 国产免费男女视频| 国产又爽黄色视频| 十八禁人妻一区二区| 国内精品久久久久精免费| 亚洲国产毛片av蜜桃av| 久久久久九九精品影院| 伦理电影免费视频| 老汉色∧v一级毛片| 91av网站免费观看| 久久伊人香网站| 免费在线观看影片大全网站| 88av欧美| 在线观看日韩欧美| 在线十欧美十亚洲十日本专区| 日本 av在线| 中亚洲国语对白在线视频| 正在播放国产对白刺激| 黄网站色视频无遮挡免费观看| 久久久水蜜桃国产精品网| 久久这里只有精品19| 精品久久久久久久人妻蜜臀av| 国产在线观看jvid| cao死你这个sao货| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 免费在线观看完整版高清| 久久久久国内视频| 日韩大尺度精品在线看网址| 波多野结衣av一区二区av| 亚洲精华国产精华精| 日本一本二区三区精品| 老熟妇仑乱视频hdxx| 色综合站精品国产| 国产视频内射| 国产精品99久久99久久久不卡| 国产又色又爽无遮挡免费看| 1024视频免费在线观看| 国产99久久九九免费精品| www.自偷自拍.com| 亚洲av电影在线进入| 午夜福利成人在线免费观看| 成人特级黄色片久久久久久久| 国产精品二区激情视频| 精品熟女少妇八av免费久了| 一区二区三区国产精品乱码| 亚洲国产精品999在线| 国产精品久久久久久人妻精品电影| 少妇的丰满在线观看| 国产免费男女视频| 国产欧美日韩一区二区三| av福利片在线| 欧美乱码精品一区二区三区| 国内少妇人妻偷人精品xxx网站 | 亚洲第一欧美日韩一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲av五月六月丁香网| av欧美777| 久久精品91无色码中文字幕| 男女那种视频在线观看| 精品不卡国产一区二区三区| 男女做爰动态图高潮gif福利片| 免费看十八禁软件| 久久国产精品人妻蜜桃| 一区二区三区激情视频| АⅤ资源中文在线天堂| 久久久久久亚洲精品国产蜜桃av| 激情在线观看视频在线高清| 久久精品国产亚洲av香蕉五月| 欧美黑人巨大hd| 两个人看的免费小视频| 亚洲专区中文字幕在线| 变态另类丝袜制服| 国产日本99.免费观看| 日韩国内少妇激情av| 高清毛片免费观看视频网站| 久久草成人影院| 熟妇人妻久久中文字幕3abv| 亚洲成人精品中文字幕电影| 免费看美女性在线毛片视频| 欧美丝袜亚洲另类 | 欧美性猛交黑人性爽| av片东京热男人的天堂| 757午夜福利合集在线观看| 宅男免费午夜| 欧美乱妇无乱码| 精品久久久久久久人妻蜜臀av| 男女午夜视频在线观看| 免费搜索国产男女视频| 伦理电影免费视频| 好男人电影高清在线观看| 日韩欧美 国产精品| 在线国产一区二区在线| 18禁黄网站禁片免费观看直播| a在线观看视频网站| 久久久久免费精品人妻一区二区 | a级毛片a级免费在线| 国产单亲对白刺激| 国产精品99久久99久久久不卡| 一级a爱片免费观看的视频| 欧美激情高清一区二区三区| 最近最新中文字幕大全电影3 | 国产精品爽爽va在线观看网站 | 国产乱人伦免费视频| 亚洲久久久国产精品| 久久香蕉激情| 日本成人三级电影网站| av在线天堂中文字幕| 成人一区二区视频在线观看| 美女大奶头视频| 日本一区二区免费在线视频| 亚洲 欧美 日韩 在线 免费| 亚洲国产精品成人综合色| 看黄色毛片网站| 免费在线观看完整版高清| 色尼玛亚洲综合影院| 91成人精品电影| 久久久久九九精品影院| 日韩精品中文字幕看吧| 岛国视频午夜一区免费看| 欧美中文日本在线观看视频| 亚洲色图av天堂| АⅤ资源中文在线天堂| 十分钟在线观看高清视频www| 亚洲色图av天堂| 女生性感内裤真人,穿戴方法视频| 国产精品亚洲av一区麻豆| 变态另类成人亚洲欧美熟女| 变态另类丝袜制服| 成年女人毛片免费观看观看9| 精品少妇一区二区三区视频日本电影| 哪里可以看免费的av片| 日韩一卡2卡3卡4卡2021年| www日本在线高清视频| 欧美午夜高清在线| 午夜老司机福利片| 又黄又爽又免费观看的视频| 中文字幕人妻丝袜一区二区| aaaaa片日本免费| 又紧又爽又黄一区二区| 国产亚洲欧美在线一区二区| 精品无人区乱码1区二区| 亚洲精品在线观看二区| 母亲3免费完整高清在线观看| 国产成人啪精品午夜网站| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 国产人伦9x9x在线观看| 国产精品久久久久久亚洲av鲁大| 狠狠狠狠99中文字幕| 法律面前人人平等表现在哪些方面| 黄网站色视频无遮挡免费观看| 视频区欧美日本亚洲| 欧美日韩亚洲国产一区二区在线观看| 熟女电影av网| av在线天堂中文字幕| 久久久久精品国产欧美久久久| 97碰自拍视频| 日韩av在线大香蕉| 我的亚洲天堂| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看免费日韩欧美大片| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久黄片| 国产精品亚洲av一区麻豆| 色哟哟哟哟哟哟| 麻豆av在线久日| 搡老熟女国产l中国老女人| 亚洲,欧美精品.| 国产亚洲欧美98| 国产成人欧美| 日韩av在线大香蕉| 在线国产一区二区在线| 亚洲激情在线av| 久久久久国产一级毛片高清牌| 悠悠久久av| 欧美成人性av电影在线观看| 一级毛片女人18水好多| 国产av一区二区精品久久| 美女 人体艺术 gogo| 国产精品免费视频内射| 国产成人精品久久二区二区91| 人人妻人人看人人澡| 日韩 欧美 亚洲 中文字幕| 18禁黄网站禁片午夜丰满| 午夜精品久久久久久毛片777| 亚洲av中文字字幕乱码综合 | 母亲3免费完整高清在线观看| 9191精品国产免费久久| 欧美激情 高清一区二区三区| 首页视频小说图片口味搜索| 丝袜在线中文字幕| 亚洲国产精品成人综合色| 99久久99久久久精品蜜桃| 久久中文看片网| 免费看日本二区| 男人舔奶头视频| 精品欧美国产一区二区三| 亚洲欧美精品综合一区二区三区| 午夜福利免费观看在线| 夜夜躁狠狠躁天天躁| 少妇被粗大的猛进出69影院| 最新在线观看一区二区三区| 亚洲av五月六月丁香网| a级毛片a级免费在线| 精品国产乱子伦一区二区三区| 丁香欧美五月| 99热6这里只有精品| 亚洲一区二区三区色噜噜| 1024视频免费在线观看| 久久精品亚洲精品国产色婷小说| 在线天堂中文资源库| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 无人区码免费观看不卡| 熟女电影av网| 人成视频在线观看免费观看| netflix在线观看网站| 欧美黑人欧美精品刺激| 欧洲精品卡2卡3卡4卡5卡区| 久久人人精品亚洲av| 久久午夜亚洲精品久久| 国产精品亚洲一级av第二区| 国产成人精品久久二区二区91| 热re99久久国产66热| 黄片小视频在线播放| 亚洲精品国产一区二区精华液| 午夜影院日韩av| 国产麻豆成人av免费视频| 亚洲精品中文字幕在线视频| 久久精品91无色码中文字幕| 99国产精品一区二区三区| 日韩视频一区二区在线观看| 欧美一级a爱片免费观看看 | 久久精品人妻少妇| 久久性视频一级片| 少妇 在线观看| av视频在线观看入口| 精品高清国产在线一区| 国产精品九九99| 国产成人av激情在线播放| 久久天堂一区二区三区四区| 成年女人毛片免费观看观看9| 自线自在国产av| 后天国语完整版免费观看| 变态另类成人亚洲欧美熟女| 狠狠狠狠99中文字幕| 色综合婷婷激情| 亚洲久久久国产精品| 最近最新中文字幕大全电影3 | 在线免费观看的www视频| 午夜久久久在线观看| 搡老岳熟女国产| 99热6这里只有精品| www.自偷自拍.com| 最近在线观看免费完整版| 亚洲中文字幕日韩| 欧美zozozo另类| 熟妇人妻久久中文字幕3abv| 免费高清在线观看日韩| 久久国产亚洲av麻豆专区| 日本精品一区二区三区蜜桃| 高清在线国产一区| av欧美777| 一卡2卡三卡四卡精品乱码亚洲| av中文乱码字幕在线| 欧美人与性动交α欧美精品济南到| 一边摸一边做爽爽视频免费| 久久亚洲真实| 老熟妇乱子伦视频在线观看| 一级片免费观看大全| 精品福利观看| 色播亚洲综合网| 成年版毛片免费区| 免费搜索国产男女视频| 婷婷亚洲欧美| 亚洲五月天丁香|