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

    基于MaxEnt和DIVA-GIS的白蠟枯梢病菌在中國(guó)的潛在適生區(qū)分析

    2015-01-10 02:12:02羅志萍潘緒斌李來梅
    湖南農(nóng)業(yè)科學(xué) 2015年10期
    關(guān)鍵詞:白蠟適生區(qū)病菌

    羅志萍,潘緒斌,王 聰,李來梅

    (1.湖南出入境檢驗(yàn)檢疫局,湖南 長(zhǎng)沙 410004;2.中國(guó)檢驗(yàn)檢疫科學(xué)研究院,北京100176;3.湖南省微生物研究院,湖南 長(zhǎng)沙 410009)

    基于MaxEnt和DIVA-GIS的白蠟枯梢病菌在中國(guó)的潛在適生區(qū)分析

    羅志萍1,潘緒斌2,王 聰2,李來梅3

    (1.湖南出入境檢驗(yàn)檢疫局,湖南 長(zhǎng)沙 410004;2.中國(guó)檢驗(yàn)檢疫科學(xué)研究院,北京100176;3.湖南省微生物研究院,湖南 長(zhǎng)沙 410009)

    為了開展危險(xiǎn)性林木病害白蠟枯梢病的風(fēng)險(xiǎn)分析和檢疫工作,研究對(duì)白蠟枯梢病菌在中國(guó)的潛在適生區(qū)進(jìn)行了預(yù)測(cè)。利用MaxEnt軟件和DIVA-GIS數(shù)據(jù)庫(kù),選擇19個(gè)環(huán)境因子變量和2.5 min的空間分辨率,并對(duì)預(yù)測(cè)結(jié)果的可靠性進(jìn)行ROC曲線法分析評(píng)價(jià),同時(shí)采用Jackknife刀切法分析了各環(huán)境變量對(duì)預(yù)測(cè)結(jié)果的影響。結(jié)果表明:該病菌在我國(guó)局部地區(qū)有一定適生風(fēng)險(xiǎn);根據(jù)得到的潛在適生區(qū)預(yù)測(cè)圖,白蠟枯梢病菌在中國(guó)的高度適生區(qū)有臺(tái)灣中部,中度適生區(qū)有西藏西南邊陲林區(qū)、湖北神農(nóng)架和宜昌三峽林區(qū)、湖南武陵山和張家界林區(qū)、廣東、廣西、重慶、貴州、江西、福建的局部等地;各環(huán)境因子貢獻(xiàn)率分析表明最冷季降水量對(duì)白蠟枯梢病菌的適生區(qū)分布影響最大,研究預(yù)測(cè)結(jié)果AUC值為0.932,表明預(yù)測(cè)結(jié)果可信度高。

    白蠟樹枯梢??;白蠟鞘孢菌(Chalara fraxinea);適生性分析;MaxEnt

    白蠟樹枯梢病是由白蠟鞘孢菌Chalara fraxinea(有性階段Hymenoscyphus pseudoalbidus)引起的危險(xiǎn)性林木病害,主要靠繁殖材料和木材進(jìn)行遠(yuǎn)距離傳播。白蠟枯梢病在歐洲大陸擴(kuò)散蔓延,使歐洲各國(guó)遭受了巨大的經(jīng)濟(jì)損失。2007年該病被歐洲和地中海區(qū)域植物保護(hù)組織(EPPO)秘書處列入EPPO有害生物預(yù)警名單。2012年該病害在英國(guó)迅速擴(kuò)散,導(dǎo)致空前的林業(yè)疫情后,歐洲國(guó)家紛紛出臺(tái)禁令,禁止進(jìn)口和境內(nèi)運(yùn)輸白蠟樹種子、植株、樹木。據(jù)有關(guān)報(bào)道稱,只要有白蠟鞘孢菌寄主植物分布的地區(qū),該菌就能夠定殖。該菌潛在的寄主白蠟屬植物在中國(guó)分布約有20種以上,該菌一旦傳入我國(guó),定殖風(fēng)險(xiǎn)很大。

    為了順利開展白蠟枯梢病菌風(fēng)險(xiǎn)分析工作,研究利用MaxEnt軟件和DIVA-GIS數(shù)據(jù)庫(kù)對(duì)該有害生物在中國(guó)的潛在適生區(qū)進(jìn)行了科學(xué)預(yù)測(cè)。MaxEnt是利用最大熵模型(Maximum EntropyModel)通過已知物種分布地區(qū)和多種環(huán)境數(shù)據(jù)預(yù)測(cè)物種潛在分布范圍的一種軟件。該軟件是Phillips等以生態(tài)位理論為基礎(chǔ)提出來的,它可以用來分析驗(yàn)證物種已有的分布數(shù)據(jù),同時(shí)也可分析預(yù)測(cè)該物種可能分布范圍,近年來多用于入侵物種和珍稀瀕危野生物種的分布預(yù)測(cè)研究[1-6]。

    1 數(shù)據(jù)采集和整理

    1.1 數(shù)據(jù)來源

    1.1.1白蠟枯梢病菌分布發(fā)生數(shù)據(jù)通過查閱國(guó)內(nèi)外公開發(fā)表的文獻(xiàn)和國(guó)外官方網(wǎng)站信息[7-15],獲得白蠟枯梢病菌分布為害地區(qū)的大致信息。目前,在EPPO地區(qū)已經(jīng)鑒定發(fā)現(xiàn)C. fraxinea的國(guó)家有奧地利、比利時(shí)、捷克、芬蘭、法國(guó)、德國(guó)、匈牙利、意大利、立陶宛、荷蘭、挪威、波蘭、羅馬尼亞、俄國(guó)(加里寧格勒)、斯洛伐克、斯洛文尼亞、瑞典、英國(guó)。此外,觀察到該病害癥狀的國(guó)家有丹麥、愛沙尼亞、拉脫維亞、瑞士。亞州地區(qū),韓國(guó)曾有該病發(fā)生的報(bào)道。非洲、美洲和大洋洲均無發(fā)生記錄。將這些已有的分布信息盡可能轉(zhuǎn)換為經(jīng)緯度數(shù)據(jù)進(jìn)行保存。分別用該有害生物的Ash Dieback、Chalara fraxinea和Hymenoscyphus pseudoalbidus等名稱作為關(guān)鍵詞,在全球生物多樣性信息網(wǎng)絡(luò)(GBIF)和國(guó)際農(nóng)業(yè)與生物科學(xué)中心(CABI)中進(jìn)行搜索、對(duì)比和下載已有的分布記錄數(shù)據(jù),保存為Excel文件。經(jīng)過篩選甄別和轉(zhuǎn)換整理,最終獲得白蠟枯梢病有關(guān)的分布記錄(具體到經(jīng)緯度)1 812條。將這些數(shù)據(jù)按照MaxEnt軟件分析的要求轉(zhuǎn)換為csv格式文件保存。csv格式文件中每條記錄包含3個(gè)字段,分別為物種名、經(jīng)度和緯度,各字段間用逗號(hào)分隔。

    1.1.2環(huán)境和地圖數(shù)據(jù)采用WORLD-CLIMATE下載背景環(huán)境數(shù)據(jù),含19個(gè)環(huán)境因子變量,選擇2.5 min的空間分辨率。將氣候環(huán)境數(shù)據(jù)保存為獨(dú)立數(shù)據(jù)文件。在DIVA-GIS網(wǎng)站系統(tǒng)下載并處理最新的地理信息數(shù)據(jù)。

    1.2 數(shù)據(jù)處理

    MaxEnt軟件需要在Java環(huán)境中運(yùn)行。將整理好的環(huán)境數(shù)據(jù)和物種分布數(shù)據(jù)導(dǎo)入MaxEnt軟件,選擇適當(dāng)變量因子和參數(shù),運(yùn)行軟件,繪制出ASCII格式的白蠟枯梢病菌的世界預(yù)測(cè)分布圖,并對(duì)預(yù)測(cè)結(jié)果的可靠性進(jìn)行分析評(píng)價(jià),用Jackknife刀切法分析各個(gè)變量對(duì)預(yù)測(cè)結(jié)果的貢獻(xiàn)率,進(jìn)而確定影響物種分布的關(guān)鍵環(huán)境因子。將世界預(yù)測(cè)分布圖與DIVA-GIS中的中國(guó)地圖疊加切割處理,則得到白蠟枯梢病菌在中國(guó)的適生區(qū)預(yù)測(cè)圖。

    2 結(jié)果與分析

    2.1 預(yù)測(cè)結(jié)果

    在DIVA-GIS軟件中加入MaxEnt運(yùn)行得出的預(yù)測(cè)結(jié)果,將該柵格數(shù)據(jù)進(jìn)行分類處理。按照適生指數(shù)P從低到高,將適生區(qū)依次分為4個(gè)等級(jí),P<0.1為非適生區(qū),0.1≤P<0.3為低度適生區(qū),0.3≤P≤0.5為中度適生區(qū),P>0.5為高度適生區(qū)。最終生成白蠟枯梢病菌潛在的世界分布圖和中國(guó)潛在分布圖,分別見圖1和圖2。

    從圖1中可以看出,預(yù)測(cè)結(jié)果與目前歐洲大面積發(fā)生白蠟枯梢病菌的現(xiàn)狀相符,同時(shí)還預(yù)測(cè)了該菌在世界范圍內(nèi)潛在的適生區(qū)。

    由圖2可知,白蠟枯梢病菌在我國(guó)局部地區(qū)有一定適生風(fēng)險(xiǎn)。在我國(guó)發(fā)生潛在地區(qū)中,高度適生區(qū)有臺(tái)灣中部;中度適生區(qū)有位于喜馬拉雅山南坡的西藏西南邊陲林區(qū)(毗鄰尼泊爾)、湖北神農(nóng)架和宜昌三峽林區(qū)、湖南武陵山和張家界林區(qū)、廣東、廣西、重慶、貴州、江西、福建的局部地區(qū);低度適生區(qū)有西藏東部、云南、四川、甘肅、陜西、安徽、浙江、福建、湖南局部、湖北局部;非適生地區(qū)有新疆、內(nèi)蒙、青海、華北、東北地區(qū)。

    圖1 MaxEnt預(yù)測(cè)白蠟枯梢病菌世界潛在分布

    圖2 MaxEnt預(yù)測(cè)白蠟枯梢病菌中國(guó)潛在分布

    2.2 影響白蠟枯梢病菌適生的主要環(huán)境因子

    各環(huán)境變量對(duì)預(yù)測(cè)結(jié)果的Jackknife刀切法分析和貢獻(xiàn)率分別見圖3和表1。從表1和圖3中可以看出,對(duì)白蠟枯梢病菌預(yù)測(cè)結(jié)果影響較大的環(huán)境因子有bio 19、bio 15、bio 7和bio 8,分別代表最冷季降水量、降雨量變化方差、年溫變化范圍和最濕季平均氣溫。其中,最冷季降水量對(duì)預(yù)測(cè)結(jié)果的貢獻(xiàn)率最大,該氣候因子對(duì)白蠟枯梢病菌的適生區(qū)分布影響最大。

    2.3 預(yù)測(cè)結(jié)果的評(píng)價(jià)

    預(yù)測(cè)結(jié)果可靠性ROC曲線評(píng)價(jià)見圖4。在ROC曲線圖中,橫坐標(biāo)為假陽(yáng)性率,縱坐標(biāo)為靈敏度,曲線下面積為評(píng)價(jià)預(yù)測(cè)結(jié)果可靠性的AUC值,AUC取值范圍為0~1,值越大說明環(huán)境變量因子與該物種地理分布相關(guān)性越大,預(yù)測(cè)效果也越好。AUC值在0.5~0.7之間,認(rèn)為預(yù)測(cè)結(jié)果可靠性低;AUC值在0.7~0.9之間,認(rèn)為預(yù)測(cè)結(jié)果可靠性中等;AUC值在0.9以上則視為預(yù)測(cè)結(jié)果可靠性高。由圖4可知,該研究預(yù)測(cè)的白蠟枯梢病菌適生區(qū)AUC值為0.932,該預(yù)測(cè)結(jié)果高度可靠。

    表1 環(huán)境因子變量對(duì)白蠟枯梢病菌適生區(qū)預(yù)測(cè)結(jié)果的貢獻(xiàn)率

    圖3 白蠟枯梢病菌MaxEnt預(yù)測(cè)結(jié)果Jackknife刀切分析

    圖4 白蠟枯梢病菌MaxEnt預(yù)測(cè)結(jié)果的AUC精度評(píng)價(jià)

    3 討 論

    從預(yù)測(cè)結(jié)果來看,白蠟枯梢病菌主要在我國(guó)南方的原始森林地區(qū)有分布定植風(fēng)險(xiǎn),應(yīng)當(dāng)引起監(jiān)管和檢疫決策部門的警惕。

    理論上,用足夠全面的物種現(xiàn)有分布數(shù)據(jù)和詳盡的環(huán)境數(shù)據(jù)進(jìn)行該物種的潛在適生區(qū)預(yù)測(cè)是科學(xué)可靠的。目前可用的預(yù)測(cè)軟件很多,但MaxEnt軟件利用生態(tài)位最大熵原理設(shè)計(jì)運(yùn)算模式,綜合考慮了物種與復(fù)雜環(huán)境的關(guān)聯(lián)關(guān)系,得出的預(yù)測(cè)數(shù)據(jù)相較于其他分析軟件更科學(xué)更可靠,更適合植物病原菌的分布預(yù)測(cè)。但其前提條件是必須有充分詳盡的數(shù)據(jù)資料做支撐,需要盡可能多地搜集物種現(xiàn)有分布數(shù)據(jù)。前期數(shù)據(jù)搜集和整理工作的質(zhì)量直接影響到預(yù)測(cè)結(jié)果的可靠性。

    入侵生物適生區(qū)的預(yù)測(cè)工作是植物檢疫有害生物風(fēng)險(xiǎn)分析工作的基礎(chǔ),準(zhǔn)確可靠的預(yù)測(cè)結(jié)果可以給國(guó)際農(nóng)產(chǎn)品貿(mào)易檢疫準(zhǔn)入談判提供科學(xué)依據(jù)和技術(shù)支撐,也可以為農(nóng)業(yè)部門的有害生物預(yù)防和監(jiān)測(cè)工作提供指導(dǎo),我國(guó)應(yīng)該加強(qiáng)該項(xiàng)技術(shù)工作的開展。

    [1] 車 樂,曹 博,白成科,等. 基于MaxEnt 和ArcGIS對(duì)太白米的潛在分布預(yù)測(cè)及適宜性評(píng)價(jià)[J]. 生態(tài)學(xué)雜志,2014,33(6):1623-1628.

    [2] 朱耿平,劉國(guó)卿,卜文俊,等. 生態(tài)位模型的基本原理及其在生物多樣性保護(hù)中的應(yīng)用[J]. 生物多樣性,2013,21(1):90-98.

    [3] 郝朝運(yùn),譚樂和,范 睿,等. 利用最大熵模型預(yù)測(cè)藥用植物海南蒟的潛在地理布局[J].熱帶作物學(xué)報(bào),2011,32(8):1561-1566.

    [4] 王 穎,章桂明,楊偉東,等. 基于MaxEnt的大豆南北方莖潰瘍病菌在中國(guó)適生區(qū)的預(yù)測(cè)[J]. 植物檢疫,2009,(4):14-17.

    [5] 陳 晨,陳 娟,胡白石,等. 梨火疫病菌在中國(guó)的潛在分布及入侵風(fēng)險(xiǎn)分析[J]. 中國(guó)農(nóng)業(yè)科學(xué),2007,40(5):940-947.

    [6] 王運(yùn)生,謝丙炎,萬(wàn)方浩,等. Roc曲線分析在評(píng)價(jià)入侵物種分布模型中的應(yīng)用[J]. 生物多樣性,2007,15(4):365-372.

    [7] Kowalski T. Chalara fraxinea sp. nov. associated with dieback of ash (Fraxinus excelsior) in Poland[J]. Forest Pathology,2006,36:264-270.

    [8] European Plant Protection Organization (EPPO). Chalara fraxinea, Ash Dieback. (2007-3)[ EB/OL]. http://www.eppo.org/QUARANTINE/Alert_List/fungi/Chalara_fraxinea.htm,2012-11-06.

    [9] UK: Forestry Commission. Chalara dieback of ash (Chalara fraxinea) [EB/OL]. http://www.forestry.gov.uk/chalara,2012-12-25.

    [10] UK: Forestry Commission. The Plant Health (Forestry) (Amendment) Order 2012[EB/OL]. http://www.legislation.gov.uk/uksi/2012/2707/pdfs/ uksi_20122707_en.pdf/$FILE/uksi_20122707_en.pdf,2012-12-26.

    [11] Commonwealth Agricultural Bureau International (CABI). Hymenoscyphus albidus[EB/OL]. http://www.cabi.org/cpc/,2012-12-27.

    [12] UK: Forest Research. Rapid assessment of the need for a detailed Pest Risk Analysis for Chalara fraxinea[EB/OL]. http://www.fera.defra.gov. uk/plants/plantHealth/pestsDiseases/documents/chalaraFraxinea.pdf,2012-11-06.

    [13] Irish Government News Service(IGNS). Chalara disease found in young ash trees[EB/OL]. http://www.merrionstreet.ie/index.php/2012/10/ chalara-disease-found-in-young-ash-trees/?cat=12,2012-12-12.

    [14] Chalara ash dieback symptoms[EB/OL]. http://www.forestry.gov.uk/ forestry/INFD-92AHUK,2014-09-24.

    [15] Timmermann V,B?rja I,Hietala A M,et al. Ash dieback: pathogen spread and diurnal patterns of ascospore dispersal, with special emphasis on Norway[J]. EPPO/OEPP Bulletin,2011,41(1):14-20.

    (責(zé)任編輯:成 平)

    《湖南農(nóng)業(yè)科學(xué)》網(wǎng)上投稿系統(tǒng)開通公告

    為適應(yīng)期刊數(shù)字化出版和方便作者投稿,《湖南農(nóng)業(yè)科學(xué)》對(duì)原網(wǎng)站和投稿系統(tǒng)進(jìn)行了改版,并于2015年9月25日正式對(duì)外試運(yùn)行,網(wǎng)址:Http://www.hnnykx.com 。 網(wǎng)上投稿系統(tǒng)的啟用不僅有利于作者投稿、查詢、在線修改稿件信息,而且方便專家審稿,縮短審稿周期和數(shù)字化出版周期。在《湖南農(nóng)業(yè)科學(xué)》網(wǎng)上投稿系統(tǒng)試運(yùn)行期間,仍保留原有的投稿方式,投稿郵箱:hnnykx@vip.163.com。如有任何問題和建議,請(qǐng)及時(shí)向本刊編輯部垂詢和反饋。聯(lián)系電話:0731-84691322,0731-84693060。感謝大家對(duì)《湖南農(nóng)業(yè)科學(xué)》的支持與幫助!

    湖南農(nóng)業(yè)科學(xué)雜志社

    2015年9月

    Potential Distribution of Chalara fraxinea in China Predicted by MaxEnt and DIVA-GIS

    LUO Zhi-ping1,PAN Xu-bin2,WANG Cong2,LI Lai-mei3
    (1. Hunan Entry-Exit Inspection and Quarantine Bureau, Changsha 410004, PRC; 2. China Academy of Inspection and Quarantine, Beijing 100176, PRC; 3. Hunan Microbiology Institute , Changsha 410009, PRC)

    For the risk analysis of the disease ash dieback, this study predicted potential distribution regions in China of the fungusChalara fraxineaby using MaxEnt software and DIVA-GIS database. 19 environmental factors and 2.5 min spatial resolution were selected, and the ROC curve and Jackknife analysis method were used to evaluate the reliability of prediction results and the infuence of environmental variables on the prediction results, respectively. The results showed that C.fraxineacould distribute in some areas in China. According to the research fgure, the middle of Taiwan was the highly suitable area, and moderately suitable areas included those forest regions located in southwest Tibet, Hubei Shennongjia and Yichang Three Gorges, Hunan Wuling Mountain and Zhangjiajie, as well as some parts of Guangdong, Guangxi, Chongqing, Guizhou, Jiangxi, Fujian etc. The precipitation in the coldest quarter was the most important environmental factor which could infuence the potential distribution of C.fraxinea. The AUC value of the predication results was 0.932, which indicated that the results were highly reliable.

    ash dieback;Chalara fraxinea; potential distribution analysis; predication; MaxEnt

    S41-30

    :A

    :1006-060X(2015)10-0083-05

    10.16498/j.cnki.hnnykx.2015.10.025

    2015-09-12

    國(guó)家質(zhì)檢總局科技計(jì)劃項(xiàng)目(2014JK015)

    羅志萍(1972-),女,青海貴德縣人,高級(jí)農(nóng)藝師,主要從事進(jìn)出境植物檢疫和外來有害生物風(fēng)險(xiǎn)分析工作。

    猜你喜歡
    白蠟適生區(qū)病菌
    氣候變化下中國(guó)蒟蒻薯科箭根薯的地理分布格局預(yù)測(cè)
    未來氣候條件下當(dāng)歸適生區(qū)預(yù)測(cè)及時(shí)空變化分析
    氣候變化下瀕危植物半日花在中國(guó)的潛在分布
    白蠟種質(zhì)資源收集保存及應(yīng)用研究
    白蠟外齒莖蜂防治試驗(yàn)研究
    巴拉圭瓜多竹適生區(qū)分布研究
    小病菌影響鴉片戰(zhàn)爭(zhēng)
    特別健康(2018年3期)2018-07-04 00:40:24
    對(duì)節(jié)白蠟盆景欣賞(下)
    花卉(2016年12期)2016-12-10 08:06:48
    頭狀莖點(diǎn)霉病菌的新寄主高粱及病菌的檢疫鑒定(內(nèi)文第98~101頁(yè))圖版
    油茶炭疽病菌拮抗木霉菌的分離與篩選
    19禁男女啪啪无遮挡网站| 天堂网av新在线| 最好的美女福利视频网| 搡老妇女老女人老熟妇| 日本黄色片子视频| 黄色日韩在线| 狠狠狠狠99中文字幕| 欧美日韩瑟瑟在线播放| 我要搜黄色片| 成人永久免费在线观看视频| av在线天堂中文字幕| 亚洲中文日韩欧美视频| 99在线人妻在线中文字幕| 国产色婷婷99| 19禁男女啪啪无遮挡网站| 成人精品一区二区免费| or卡值多少钱| netflix在线观看网站| 国产精品精品国产色婷婷| 国产伦在线观看视频一区| 制服丝袜大香蕉在线| 久久久久国产精品人妻aⅴ院| 亚洲欧美日韩卡通动漫| 欧美不卡视频在线免费观看| 一个人观看的视频www高清免费观看| 免费无遮挡裸体视频| 18美女黄网站色大片免费观看| 亚洲中文字幕日韩| 国产成+人综合+亚洲专区| 一二三四社区在线视频社区8| 啦啦啦韩国在线观看视频| 综合色av麻豆| 欧美绝顶高潮抽搐喷水| 国产成年人精品一区二区| 有码 亚洲区| 日韩欧美精品v在线| 国产亚洲精品久久久久久毛片| 看片在线看免费视频| 色av中文字幕| 国产精品一区二区免费欧美| 国产精品女同一区二区软件 | 国产精品1区2区在线观看.| 男插女下体视频免费在线播放| 99久久久亚洲精品蜜臀av| 别揉我奶头~嗯~啊~动态视频| 欧美又色又爽又黄视频| 女人被狂操c到高潮| 男女做爰动态图高潮gif福利片| 亚洲精品日韩av片在线观看 | 午夜免费成人在线视频| av福利片在线观看| 一区二区三区国产精品乱码| 亚洲国产精品久久男人天堂| 欧美3d第一页| 国产免费av片在线观看野外av| 十八禁人妻一区二区| 亚洲中文字幕一区二区三区有码在线看| 好男人在线观看高清免费视频| 好男人电影高清在线观看| 校园春色视频在线观看| 亚洲av免费在线观看| 国产精品亚洲美女久久久| 亚洲精华国产精华精| 老熟妇仑乱视频hdxx| 婷婷丁香在线五月| 亚洲av免费在线观看| 国产亚洲精品一区二区www| 日韩欧美精品免费久久 | 欧美绝顶高潮抽搐喷水| 给我免费播放毛片高清在线观看| 97碰自拍视频| 国产精品永久免费网站| 久久久久久大精品| 麻豆成人av在线观看| 亚洲成人久久爱视频| 国产精品香港三级国产av潘金莲| 午夜精品一区二区三区免费看| 99国产综合亚洲精品| 欧美区成人在线视频| 两人在一起打扑克的视频| 日韩欧美一区二区三区在线观看| 99久久精品一区二区三区| 欧美zozozo另类| 12—13女人毛片做爰片一| 欧美极品一区二区三区四区| 亚洲精品美女久久久久99蜜臀| 久久精品亚洲精品国产色婷小说| 性色avwww在线观看| 2021天堂中文幕一二区在线观| av黄色大香蕉| 免费看美女性在线毛片视频| 麻豆成人av在线观看| 欧美丝袜亚洲另类 | 久久久色成人| 欧美激情在线99| 国产乱人伦免费视频| 国产成人av教育| 天天添夜夜摸| 麻豆成人午夜福利视频| 免费无遮挡裸体视频| 欧美黑人巨大hd| 精品人妻偷拍中文字幕| www国产在线视频色| 97超视频在线观看视频| 精品久久久久久久人妻蜜臀av| 亚洲av熟女| 欧美一区二区精品小视频在线| 熟女少妇亚洲综合色aaa.| 九色成人免费人妻av| 亚洲无线观看免费| 久久九九热精品免费| 十八禁人妻一区二区| 国产精品美女特级片免费视频播放器| 亚洲欧美日韩高清专用| 国产精品自产拍在线观看55亚洲| 一级a爱片免费观看的视频| 国产黄片美女视频| 国产熟女xx| 久久久久精品国产欧美久久久| 日本免费a在线| 日韩免费av在线播放| 精品人妻一区二区三区麻豆 | 国产成人欧美在线观看| 免费高清视频大片| 成年版毛片免费区| 国产日本99.免费观看| 亚洲av第一区精品v没综合| 精品久久久久久久末码| 免费人成视频x8x8入口观看| 久久精品综合一区二区三区| 国产伦精品一区二区三区视频9 | 精品99又大又爽又粗少妇毛片 | 国产日本99.免费观看| 国产精品三级大全| 精品一区二区三区视频在线观看免费| 欧美另类亚洲清纯唯美| 国产精品精品国产色婷婷| 12—13女人毛片做爰片一| 亚洲第一电影网av| 成人午夜高清在线视频| 日本五十路高清| 成人欧美大片| 青草久久国产| www.色视频.com| 国内精品一区二区在线观看| 国内精品久久久久久久电影| 亚洲av成人av| 女同久久另类99精品国产91| 老司机福利观看| 欧美日韩乱码在线| av中文乱码字幕在线| 国产精品98久久久久久宅男小说| 亚洲专区中文字幕在线| 国产蜜桃级精品一区二区三区| 九色成人免费人妻av| 我的老师免费观看完整版| 国产在视频线在精品| 亚洲av熟女| 国产真人三级小视频在线观看| 欧美乱码精品一区二区三区| 亚洲人与动物交配视频| 波多野结衣高清作品| 亚洲专区国产一区二区| 成人亚洲精品av一区二区| 男女之事视频高清在线观看| 国产三级在线视频| 色综合婷婷激情| 国产成人aa在线观看| 一二三四社区在线视频社区8| 日本黄色视频三级网站网址| 在线播放国产精品三级| 嫩草影院入口| 少妇丰满av| 欧美不卡视频在线免费观看| 国产精品久久久久久久久免 | 小蜜桃在线观看免费完整版高清| 给我免费播放毛片高清在线观看| 国产精品一及| 亚洲电影在线观看av| 最新美女视频免费是黄的| 老鸭窝网址在线观看| avwww免费| 一个人看视频在线观看www免费 | 亚洲中文字幕日韩| 中文在线观看免费www的网站| 亚洲五月天丁香| 他把我摸到了高潮在线观看| 国产精品1区2区在线观看.| 黄色丝袜av网址大全| 高清毛片免费观看视频网站| 亚洲av美国av| 亚洲中文字幕一区二区三区有码在线看| 最近最新中文字幕大全电影3| 男女视频在线观看网站免费| 国产欧美日韩一区二区三| 少妇高潮的动态图| 桃色一区二区三区在线观看| 我要搜黄色片| 中文字幕人成人乱码亚洲影| 网址你懂的国产日韩在线| 成熟少妇高潮喷水视频| 亚洲av熟女| www.色视频.com| 国产黄片美女视频| 亚洲avbb在线观看| 18美女黄网站色大片免费观看| 丰满人妻一区二区三区视频av | 美女大奶头视频| 两性午夜刺激爽爽歪歪视频在线观看| 91麻豆av在线| 一边摸一边抽搐一进一小说| 亚洲激情在线av| 日韩欧美 国产精品| 一级毛片高清免费大全| av片东京热男人的天堂| 一夜夜www| 欧美午夜高清在线| 久久精品91蜜桃| 美女被艹到高潮喷水动态| a级一级毛片免费在线观看| 午夜激情福利司机影院| 最后的刺客免费高清国语| 午夜福利欧美成人| 亚洲美女黄片视频| 男女下面进入的视频免费午夜| 悠悠久久av| 欧美激情久久久久久爽电影| 国产成人a区在线观看| 成年免费大片在线观看| 午夜免费激情av| 免费观看人在逋| 免费高清视频大片| 3wmmmm亚洲av在线观看| 日韩欧美 国产精品| 99久久九九国产精品国产免费| 午夜精品一区二区三区免费看| 又爽又黄无遮挡网站| av专区在线播放| 好男人电影高清在线观看| 69人妻影院| 真实男女啪啪啪动态图| 精品电影一区二区在线| 午夜亚洲福利在线播放| 久久伊人香网站| 麻豆久久精品国产亚洲av| 亚洲内射少妇av| 成年免费大片在线观看| 一本综合久久免费| 丁香欧美五月| 夜夜爽天天搞| 欧美一级a爱片免费观看看| 国产三级中文精品| 免费观看人在逋| 首页视频小说图片口味搜索| 精品福利观看| 免费av毛片视频| 亚洲国产色片| 一级毛片女人18水好多| 久久精品91蜜桃| 免费电影在线观看免费观看| 高清毛片免费观看视频网站| 国产在视频线在精品| 精品一区二区三区视频在线观看免费| 久久久成人免费电影| 天天躁日日操中文字幕| 亚洲国产色片| 人妻久久中文字幕网| 亚洲欧美日韩东京热| 亚洲中文字幕一区二区三区有码在线看| 最新中文字幕久久久久| 久久久久久久午夜电影| 精品欧美国产一区二区三| 国产欧美日韩精品一区二区| 身体一侧抽搐| 国内揄拍国产精品人妻在线| 中国美女看黄片| 国产精品永久免费网站| 成人国产综合亚洲| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 久久天躁狠狠躁夜夜2o2o| 女人十人毛片免费观看3o分钟| 特大巨黑吊av在线直播| netflix在线观看网站| 午夜两性在线视频| 精品人妻偷拍中文字幕| 19禁男女啪啪无遮挡网站| 欧美日韩中文字幕国产精品一区二区三区| 免费观看的影片在线观看| 91麻豆av在线| 欧美日韩国产亚洲二区| 亚洲自拍偷在线| 欧美zozozo另类| 最近最新中文字幕大全电影3| 动漫黄色视频在线观看| 99在线人妻在线中文字幕| 国内精品久久久久久久电影| 伊人久久精品亚洲午夜| 日韩欧美一区二区三区在线观看| 亚洲av免费高清在线观看| 在线国产一区二区在线| 一级a爱片免费观看的视频| a在线观看视频网站| 看免费av毛片| 最好的美女福利视频网| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 12—13女人毛片做爰片一| 久久久国产精品麻豆| 少妇的丰满在线观看| 三级男女做爰猛烈吃奶摸视频| 国产一区二区在线观看日韩 | 母亲3免费完整高清在线观看| 亚洲欧美精品综合久久99| 国产探花极品一区二区| 欧美bdsm另类| 中国美女看黄片| 最近最新中文字幕大全电影3| 国产亚洲精品久久久久久毛片| 91麻豆精品激情在线观看国产| www.www免费av| 免费人成在线观看视频色| 亚洲自拍偷在线| 亚洲国产精品成人综合色| 少妇人妻精品综合一区二区 | 日韩亚洲欧美综合| 亚洲人成网站在线播放欧美日韩| 69av精品久久久久久| 在线观看舔阴道视频| 国产一区二区亚洲精品在线观看| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 亚洲激情在线av| 亚洲不卡免费看| 色噜噜av男人的天堂激情| 成人鲁丝片一二三区免费| 99国产极品粉嫩在线观看| 国产精品美女特级片免费视频播放器| 中文字幕av成人在线电影| 1024手机看黄色片| 国产一区二区三区视频了| 国产成人影院久久av| av福利片在线观看| 18美女黄网站色大片免费观看| 欧美不卡视频在线免费观看| 嫩草影院入口| 欧美区成人在线视频| 精品久久久久久久末码| 在线观看免费视频日本深夜| 18禁裸乳无遮挡免费网站照片| 国产精品99久久99久久久不卡| 久久6这里有精品| 国产精品99久久久久久久久| 欧美日韩乱码在线| 欧美日韩瑟瑟在线播放| 两个人视频免费观看高清| 特级一级黄色大片| 天美传媒精品一区二区| 欧美激情久久久久久爽电影| 母亲3免费完整高清在线观看| 成人永久免费在线观看视频| 欧美日韩乱码在线| 日韩国内少妇激情av| 我要搜黄色片| 久久精品国产综合久久久| 有码 亚洲区| 99久久成人亚洲精品观看| 久久99热这里只有精品18| 99热这里只有精品一区| 久久精品国产自在天天线| 人妻久久中文字幕网| 最好的美女福利视频网| 欧美色视频一区免费| 亚洲av五月六月丁香网| 国产三级黄色录像| 亚洲久久久久久中文字幕| 伊人久久大香线蕉亚洲五| 久久亚洲真实| 亚洲中文日韩欧美视频| 日韩欧美三级三区| 欧美成人免费av一区二区三区| 国产免费一级a男人的天堂| 国产精品女同一区二区软件 | 亚洲精品在线观看二区| 欧美成人性av电影在线观看| 制服人妻中文乱码| 午夜影院日韩av| 国产精华一区二区三区| 亚洲中文日韩欧美视频| 少妇丰满av| 久久久国产成人免费| 婷婷亚洲欧美| 亚洲五月天丁香| 国产精品99久久99久久久不卡| 欧美成人a在线观看| 黄色成人免费大全| 久久久久九九精品影院| 色哟哟哟哟哟哟| 欧美日韩一级在线毛片| av天堂在线播放| 免费在线观看亚洲国产| 午夜免费成人在线视频| 欧美成人免费av一区二区三区| 最新中文字幕久久久久| 国产综合懂色| 国产一区二区在线av高清观看| 波多野结衣巨乳人妻| or卡值多少钱| 国产探花极品一区二区| 国产精品久久视频播放| 在线观看66精品国产| 午夜福利视频1000在线观看| 国产爱豆传媒在线观看| 国产亚洲欧美98| 日本与韩国留学比较| xxxwww97欧美| 日韩成人在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 88av欧美| 日日夜夜操网爽| 亚洲精品粉嫩美女一区| 亚洲午夜理论影院| 国产高清有码在线观看视频| 女同久久另类99精品国产91| 午夜精品久久久久久毛片777| 老熟妇仑乱视频hdxx| 中文字幕精品亚洲无线码一区| 亚洲精品在线观看二区| 亚洲av不卡在线观看| avwww免费| 中文在线观看免费www的网站| 国产精品一区二区三区四区久久| 国产激情偷乱视频一区二区| 听说在线观看完整版免费高清| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利在线观看吧| 国产极品精品免费视频能看的| 在线观看免费午夜福利视频| 久久国产乱子伦精品免费另类| 人妻久久中文字幕网| 他把我摸到了高潮在线观看| 99精品在免费线老司机午夜| 免费大片18禁| 亚洲专区国产一区二区| 亚洲国产色片| 日韩成人在线观看一区二区三区| 国产av一区在线观看免费| 午夜免费观看网址| 国产精品野战在线观看| 波多野结衣高清作品| 看免费av毛片| 国内毛片毛片毛片毛片毛片| 亚洲精品一卡2卡三卡4卡5卡| 美女cb高潮喷水在线观看| 日韩免费av在线播放| 亚洲一区高清亚洲精品| www.色视频.com| 国产色婷婷99| 亚洲性夜色夜夜综合| 色老头精品视频在线观看| a在线观看视频网站| 国内少妇人妻偷人精品xxx网站| 国产亚洲欧美98| 久久亚洲精品不卡| 看免费av毛片| 亚洲av一区综合| 白带黄色成豆腐渣| 久久久色成人| 无遮挡黄片免费观看| 色尼玛亚洲综合影院| 高清日韩中文字幕在线| 久久久久精品国产欧美久久久| 午夜a级毛片| 免费看日本二区| 人人妻人人澡欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线美女| 久久性视频一级片| 亚洲人成网站在线播放欧美日韩| 午夜福利在线观看吧| 欧美日韩精品网址| 少妇裸体淫交视频免费看高清| 中文字幕av在线有码专区| 国产精品,欧美在线| 欧美日韩亚洲国产一区二区在线观看| 免费看a级黄色片| 老司机午夜十八禁免费视频| 精品人妻偷拍中文字幕| 午夜免费激情av| netflix在线观看网站| 亚洲无线在线观看| 人妻丰满熟妇av一区二区三区| 国产主播在线观看一区二区| 老司机福利观看| 国产主播在线观看一区二区| 老司机福利观看| 日韩人妻高清精品专区| 欧美xxxx黑人xx丫x性爽| 日日摸夜夜添夜夜添小说| 国产精品香港三级国产av潘金莲| 亚洲成av人片免费观看| 特级一级黄色大片| 90打野战视频偷拍视频| 午夜视频国产福利| 有码 亚洲区| 99精品在免费线老司机午夜| 一二三四社区在线视频社区8| 亚洲欧美日韩高清在线视频| 一a级毛片在线观看| 精品无人区乱码1区二区| 精品不卡国产一区二区三区| 欧美日韩乱码在线| av专区在线播放| 国产欧美日韩一区二区三| bbb黄色大片| 国产极品精品免费视频能看的| 久久久久国内视频| 成人精品一区二区免费| 亚洲国产精品999在线| 亚洲精品久久国产高清桃花| 99热这里只有是精品50| 亚洲成人久久性| 听说在线观看完整版免费高清| 国产精品久久久久久久久免 | 久久6这里有精品| 黑人欧美特级aaaaaa片| 国产高潮美女av| 精品久久久久久久毛片微露脸| 亚洲成a人片在线一区二区| 欧美不卡视频在线免费观看| 婷婷六月久久综合丁香| 听说在线观看完整版免费高清| 国产亚洲欧美98| 国产老妇女一区| 亚洲内射少妇av| 欧美黄色淫秽网站| 国产乱人伦免费视频| 五月玫瑰六月丁香| 成年免费大片在线观看| 女同久久另类99精品国产91| 国产一区二区激情短视频| 免费看日本二区| 在线天堂最新版资源| 国产一区二区三区在线臀色熟女| 国产v大片淫在线免费观看| 亚洲人成网站在线播放欧美日韩| 亚洲国产欧美人成| 午夜福利成人在线免费观看| 国产99白浆流出| 日韩欧美免费精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩av在线大香蕉| 国产成人av激情在线播放| 日韩欧美精品v在线| 国产成人影院久久av| 9191精品国产免费久久| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 99久久九九国产精品国产免费| 有码 亚洲区| 精品久久久久久成人av| 人妻夜夜爽99麻豆av| 91在线精品国自产拍蜜月 | 亚洲国产精品sss在线观看| 亚洲五月婷婷丁香| 日本熟妇午夜| АⅤ资源中文在线天堂| 99久国产av精品| 国产精品综合久久久久久久免费| av视频在线观看入口| 国内精品久久久久精免费| 一a级毛片在线观看| 一进一出抽搐动态| 在线观看免费视频日本深夜| 天天添夜夜摸| 国产精品一区二区免费欧美| 成人三级黄色视频| 日韩有码中文字幕| 亚洲专区国产一区二区| 免费观看人在逋| 手机成人av网站| 免费人成视频x8x8入口观看| 美女被艹到高潮喷水动态| 一进一出好大好爽视频| 成人鲁丝片一二三区免费| 亚洲一区二区三区色噜噜| 国产高清有码在线观看视频| 免费一级毛片在线播放高清视频| 午夜a级毛片| 最新在线观看一区二区三区| 成年女人看的毛片在线观看| 久久6这里有精品| 免费av观看视频| 国产欧美日韩精品亚洲av| 99久久成人亚洲精品观看| 欧美av亚洲av综合av国产av| 51国产日韩欧美| 成年女人看的毛片在线观看| 99国产精品一区二区蜜桃av| 搡老熟女国产l中国老女人| 亚洲国产欧洲综合997久久,| 精品人妻一区二区三区麻豆 | 少妇的丰满在线观看| 一级黄色大片毛片| 日本精品一区二区三区蜜桃| 色综合欧美亚洲国产小说| 2021天堂中文幕一二区在线观| 亚洲欧美日韩高清专用| 老司机在亚洲福利影院| 宅男免费午夜| 最新在线观看一区二区三区| 99久久久亚洲精品蜜臀av| 久久久精品欧美日韩精品| 久久久久精品国产欧美久久久|