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

    基于集合模型的草地貪夜蛾的潛在分布預(yù)測

    2021-03-16 13:22:34韓欣嬈朱耿平門永亮
    生物安全學(xué)報 2021年1期
    關(guān)鍵詞:種植區(qū)夜蛾草地

    韓欣嬈, 朱耿平, 門永亮, 楊 琢

    天津師范大學(xué)生命科學(xué)學(xué)院,天津市動植物抗性重點實驗室,天津 300380

    草地貪夜蛾Spodopterafrugiperda(J. E. Smith)自2018年底遷飛入侵我國以來,已擴散蔓延到18個省市區(qū),呈繼續(xù)北遷趨勢。為應(yīng)對其對玉米等糧食種植作物的影響,我國各地相繼啟動針對草地貪夜蛾暴發(fā)的應(yīng)對方案,但依然缺少對其分布的整體認識,特別是其適宜的生態(tài)空間和對玉米種植的影響,這種認識空缺會影響我國對草地貪夜蛾的整體防控。具體來說,哪些地方需要重點防控,哪些地方不需要防控,哪些地方需要針對保護玉米來防控,以及哪些地方需要針對非玉米作物來防控等問題亟待解決。本研究采用集合模型預(yù)測,將草地貪夜蛾適宜的生態(tài)空間和玉米種植范圍相整合,以期為當(dāng)前草地貪夜蛾的精準(zhǔn)防控提供依據(jù)。

    草地貪夜蛾,隸屬鱗翅目Lepidoptera夜蛾科Noctuidae,原產(chǎn)于美洲熱帶及亞熱帶地區(qū),廣泛分布在美洲大陸(Sparks,1979)。玉米為我國重要的糧食作物之一,是草地貪夜蛾的主要寄主和危害對象,聯(lián)合國糧食及農(nóng)業(yè)組織基于2018年12個非洲國家的糧食數(shù)據(jù),估算草地貪夜蛾對玉米造成的年損失多達1770萬t(秦譽嘉等,2020; 吳孔明等,2020)。此外,草地貪夜蛾食性廣,幼蟲可取食76科350多種植物,在沒有玉米的情況下可危害小麥、花生、大豆等作物。草地貪夜蛾生態(tài)位寬度較廣,11~30 ℃是其適宜的溫度范圍,在28 ℃條件下,30 d左右即可完成一個世代;具有很強的飛行能力,成蟲每晚能飛行100 km以上,目前尚缺乏有效的措施用以截斷其傳播(吳孔明等,2020)。

    生態(tài)位模型利用物種分布點所關(guān)聯(lián)的環(huán)境變量去預(yù)測物種的生態(tài)需求和潛在分布,被廣泛應(yīng)用于入侵生物學(xué)、保護生物學(xué)、全球變化生物學(xué)和疾病風(fēng)險傳播(朱耿平等,2013)。在生物入侵應(yīng)用中,生態(tài)位模型的預(yù)測出現(xiàn)了模型模擬結(jié)果和實際分布不符合、模型轉(zhuǎn)移能力低和生態(tài)位遷移等不恰當(dāng)結(jié)論,主要是由于模型構(gòu)建方法不當(dāng)(朱耿平等,2014)。當(dāng)前模型發(fā)展依然有很多問題亟待解決,例如將物種間相互作用整合至模型的方案、模擬預(yù)測不確定性,以及所模擬的生態(tài)位成分等(喬慧捷等,2013)。目前,已知有30多種生態(tài)位模型預(yù)測方法(Norbergetal.,2019),這些模型基于不同的前提假設(shè)和理論基礎(chǔ),采用不同數(shù)據(jù)方案和模型算法去推算生態(tài)需求和潛在分布,每種模型均有其獨特優(yōu)勢和不足,當(dāng)前沒有任何一個模型能夠較好地模擬所有物種的潛在分布(Norbergetal.,2019; Qiaoetal.,2015)。

    與其他生態(tài)模型類似,生態(tài)位模型預(yù)測存在不確定性,這種不確定性與物種分布點、環(huán)境變量、模型算法和參數(shù)等密切相關(guān),其中物種分布點和模型算法對模型不確定性的影響最大(Thibaudetal.,2014)。集合預(yù)測以單一模型預(yù)測為基礎(chǔ)構(gòu)建集合模型,能夠降低單一模型預(yù)測的不確定性,被廣泛應(yīng)用于氣候變化對物種分布影響的研究(Araújoetal.,2011),但在生物入侵種方面應(yīng)用相對較少。本研究在biomod 2平臺(Thuilleretal.,2009)構(gòu)建廣義相加模型(generalized additive models, GAM)、推進式回歸樹模型(push regression tree model, GBM)、廣義線性模型(general linear model, GLM)、隨機森林模型(random forest, RF)和最大熵模型(maximum entropy models, MaxEnt)等5種常用的生態(tài)位模型(Qiaoetal.,2015):其中GAM和GLM屬于統(tǒng)計模型,其模擬結(jié)果有時很難擬合物種分布點數(shù)據(jù),GBM、RF和MaxEnt屬于機械學(xué)習(xí)模型,其模擬結(jié)果對分布點數(shù)據(jù)擬合較高,但有時模型轉(zhuǎn)移能力較低。以此5種生態(tài)位模型為基礎(chǔ)構(gòu)建集合模型來預(yù)測草地貪夜蛾的潛在分布,以期通過集合模型預(yù)測降低單一模型預(yù)測的不確定性。

    Earlyetal.(2018)采用限制性環(huán)境因子和集合模型預(yù)測方案探討了草地貪夜蛾的潛在分布,然而其預(yù)測結(jié)果缺乏參考價值,特別是對中東和東亞地區(qū);秦譽嘉等(2019)和Balochetal.(2020)分別對草地貪夜蛾在中亞地區(qū)和我國的潛在分布進行了模擬,其單一模型預(yù)測(MaxEnt模型)存在不確定性;Wangetal.(2020)采用國際應(yīng)用生物科學(xué)中心(CABI, www.cabi.org/ISC/datasheet/29810)的分布數(shù)據(jù)構(gòu)建模型,CABI以國家或省州行政中心在宏觀尺度下繪制入侵物種的大致分布范圍,這些行政中心僅代表了物種在該地區(qū)被發(fā)現(xiàn),與物種具體分布點存在較大差別,并不能詳實地反映物種的真實分布。此外,上述模擬工作均未將玉米的空間分布進行整合對草地貪夜蛾進行風(fēng)險分析。本研究采用集合模型預(yù)測探討草地貪夜蛾的適應(yīng)生態(tài)空間,并將模擬結(jié)果和基于中分辨率成像光譜儀(moderate-resolution imaging spectroradiometer, MODIS)預(yù)測的玉米種植密度相整合,探討其對玉米等糧食作物種植的影響,以期更加精準(zhǔn)地對草地貪夜蛾進行風(fēng)險分析。

    1 材料與方法

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

    草地貪夜蛾世界范圍地理分布點數(shù)據(jù)來源于數(shù)據(jù)庫Global Biodiversity Information Facility (GBIF, www.gbif.org)、Biodiversity Information Serving Our Nation (BISON, https:∥bison.usgs.gov)、Integrated Digitised Biocollections (iDigBio, www.idigbio.org)和iNaturalist(www.inaturalist.org)。研究采用R語言數(shù)據(jù)包spocc(https:∥CRAN.R-project.org/package=spocc)在上述數(shù)據(jù)庫執(zhí)行查詢指令,并對錯誤點進行清除;通過查閱文獻(姜玉英等,2019a,2019b; 劉杰等,2019),國內(nèi)總計3433個分布記錄,這些物種分布點涵蓋了草地貪夜蛾當(dāng)前的分布范圍,可對其生態(tài)空間進行分析。物種分布點數(shù)據(jù)一般存在采樣偏差,在建立模型前需去除采樣偏差,研究在Wallace平臺使用spThin數(shù)據(jù)包(Kassetal.,2018)去除采樣偏差,使物種分布點間隔大于100 km,最終選取446個物種分布點用于模型構(gòu)建和評估。

    玉米種植密度數(shù)據(jù)參考了Monfredaetal. (2008)基于MODIS遙感數(shù)據(jù)對全球玉米的種植范圍和強度的模擬方法,該玉米種植強度基于National Agricultural Statistics(FAOSTAT)、Subnational Agricultural Statistics獲得地域土地數(shù)據(jù),計算出每個單元格的總耕地面積,然后將獲得的玉米作物面積除以總耕地面積來確定單位作物面積與總耕地面積之比,最后將作物數(shù)據(jù)集比與單位作物面積之比相乘來獲得作物種植面積比例,得到玉米種植強度(Monfredaetal.,2008; Ramankuttyetal.,2008)。

    環(huán)境變量的選擇主要考慮其對物種分布作用的限制性和環(huán)境變量間的相關(guān)性(Fanetal.,2019)。在MaxEnt模型中運用刀切法對環(huán)境變量進行重要性排序,然后根據(jù)環(huán)境變量之間的相關(guān)性,利用SDMToolbox工具包除去相關(guān)性較低或者較高(Peterson相關(guān)性值<0.8)的環(huán)境變量。鑒于本次模型轉(zhuǎn)移范圍相對于模型構(gòu)建區(qū)域較小,本文側(cè)重選擇了對草地貪夜蛾分布限制非常重要的環(huán)境因子(Fengetal.,2019),即年平均溫度(bio1)、晝夜溫差月均值(bio2)、最熱月份最高溫度(bio5)、最冷月份最低溫度(bio6)、年平均降水量(bio12)、最濕月份降水量(bio13)、最干月份降水量(bio14)和年平均輻射量(annual mean radiation),這些環(huán)境因子覆蓋了年平均和極端溫濕度,被證實能夠提高昆蟲模型預(yù)測的轉(zhuǎn)移能力(Zhuetal.,2017)。氣候因子來源于WorldClim (www.worldclim.org)的生物氣候?qū)W變量,時間跨度為1970—2000年的平均值,空間分辨率為5′,與分布點記錄相吻合。

    1.2 模型構(gòu)建和評估

    在物種分布點基礎(chǔ)上構(gòu)建最小凸多邊形(minimal convex polygon, MCP),然后將MCP做400 km緩沖,作為模型構(gòu)建區(qū)域。在該區(qū)域隨機采取10000個不存在的物種分布點,采用默認參數(shù),在biomod 2平臺構(gòu)建GAM、GBM、GLM、RF和MaxEnt生態(tài)位模型,然后將其投射至全球區(qū)域,具體模型算法參見畢迎鳳等(2013)。本文采用70%的物種分布點構(gòu)建模型,30%的點用于模型檢驗,模型檢驗采用partial ROC方案(朱耿平等,2017; Petersonetal.,2018)。與傳統(tǒng)ROC曲線和AUC值(area under the curve)相比,partial ROC方案能夠賦予物種分布點更多權(quán)重,同時將物種分布點可信度融入模型評估。partial ROC方案采用AUC 比率來評估模型,AUC比率>1,表明模型優(yōu)于隨機模型,AUC比率≤1表明模型差于隨機模型。最終模型構(gòu)建采用篩選后的所有物種分布點,集合模型預(yù)測結(jié)果基于5個單一模型預(yù)測的平均值。

    最后,研究基于集合模型預(yù)測的草地貪夜蛾的潛在分布和Monfredaetal. (2008)的玉米種植密度結(jié)果構(gòu)建二維地圖(bivariate map)。研究中所用底圖來源于自然資源部標(biāo)準(zhǔn)地圖服務(wù)系統(tǒng)(http:∥bzdt.ch.mnr.gov.cn/) 的世界地圖和中國地圖。本文使用二維地圖分別展示草地貪夜蛾的潛在分布和玉米種植密度在世界范圍和我國的分布情況,而且可將二者重疊部分和非重疊部分分別展示。

    2 結(jié)果與分析

    基于AUC比率測試單一模型對30%分布點的預(yù)測顯示,所有模型均優(yōu)于隨機模型(即AUC比率>1),其中GBM模型表現(xiàn)較好(AUC比率=1.38),其次為MaxEnt模型(AUC比率=1.31),GAM、GLM和RF模型表現(xiàn)相對較差(AUC比率分別為1.28、1.23、1.29)。

    5個單一模型的全球空間預(yù)測存在較大的差異,采用集合模型預(yù)測后顯示,草地貪夜蛾的適宜生態(tài)空間和潛在分布主要存在于美國中東部、中美洲、歐洲南部、東亞地區(qū),及南半球的中緯度地區(qū)(南緯25°~50°),包括巴西南部、阿根廷北部、非洲南部和澳洲東部(圖1A)。與草地貪夜蛾的潛在分布相比,玉米的種植范圍和密度相對較小,北半球玉米優(yōu)勢種植區(qū)主要在美國中北部、歐洲南部、中國西南、華北和東北地區(qū),南半球主要分布在南美洲和非洲南部(圖1B)。

    我國玉米分布廣泛,包含三大玉米密集種植區(qū)(Zhaoetal.,2018),即東北三省、華北平原地區(qū)和西南地區(qū)(圖2),具體玉米種植區(qū)包括東北華北春玉米區(qū)、黃淮海平原夏播玉米區(qū)、西北灌溉玉米區(qū)、西南山地玉米區(qū)、南方丘陵玉米區(qū)、青藏高原玉米區(qū)六大區(qū),這些區(qū)域均被草地貪夜蛾的潛在分布區(qū)所覆蓋(圖1)。在玉米密集區(qū)中,華北平原密集種植區(qū)的草地貪夜蛾適生性相對較高,其次為東北三省和西南密集種植區(qū)(圖1)。二維地圖顯示,草地貪夜蛾和玉米種植高度重疊區(qū)位于華北平原春玉米區(qū)和黃淮海平原夏播玉米區(qū),其次為東北平原春玉米區(qū)和西南山地玉米區(qū),這些地區(qū)需要重點防范,南方丘陵玉米區(qū)和西北灌溉玉米區(qū)受草地貪夜蛾的影響相對較小(圖2)。我國東南和華南地區(qū)是水稻的密集種植區(qū),這些非玉米種植區(qū)應(yīng)關(guān)注草地貪夜蛾對非玉米類農(nóng)作物的影響(圖2)。

    將預(yù)測模型與草地貪夜蛾實際發(fā)生狀況進行比較,模型預(yù)測結(jié)果很好地模擬了草地貪夜蛾的潛在分布,與實際發(fā)生狀況相吻合。針對中國來說,在西北灌溉區(qū)(如新疆、西藏等)實際少有發(fā)生,預(yù)測表現(xiàn)出重疊性降低,風(fēng)險較??;在華北平原種植區(qū)(如河北)、東北三省種植區(qū)表現(xiàn)重疊較高,實際發(fā)生量較大;部分地區(qū)雖氣候適合但不是玉米種植區(qū),不會造成草地貪夜蛾的廣泛發(fā)生。

    3 討論

    當(dāng)前有大于33種生態(tài)位模型,新的算法亦在不斷推出(Carlson,2020),這些模型在構(gòu)建過程中所基于的生態(tài)學(xué)理論和前提假設(shè)不相同,所得出的結(jié)果存在差異。GBM模型是一種在回歸樹上運用boosting的模型(Elithetal. ,2008);MaxEnt采用最大熵值算法來推算物種分布,是當(dāng)前最流行的生態(tài)位模型算法(Phillipsetal.,2006)。早期模型比較研究結(jié)果顯示MaxEnt、GBM和RF模型的預(yù)測優(yōu)于其他模型(Elithetal.,2006),與本研究結(jié)果類似。然而,最近模型系統(tǒng)比較研究顯示,沒有任何模型能適應(yīng)所有物種分布模擬(Norbergetal.,2019;Qiaoetal.,2015)。集合模型預(yù)測可能是當(dāng)前解決模型預(yù)測不確定性的重要途徑,然而,基于美國白蛾的預(yù)測顯示集合模型預(yù)測亦存在不確定性(Zhu & Peterson,2017)。

    圖1 世界范圍內(nèi)草地貪夜蛾的潛在分布和玉米種植密度[審圖號:GS(2021)923號]

    圖2 草地貪夜蛾的潛在分布和玉米種植密度在我國的二維分布圖[審圖號:GS(2021)923號]

    玉米是草地貪夜蛾重要的寄主植物,我國為玉米生產(chǎn)大國,入侵我國的草地貪夜蛾為“玉米型”,將玉米種植范圍和強度整合草地貪夜蛾的潛在分布進行分析十分必要。借鑒Earlyetal. (2018),本模擬采用了草地貪夜蛾的真實物種分布點和集合模型預(yù)方案,并將模擬結(jié)果和基于MODIS預(yù)測的玉米種植密度相整合,研究結(jié)果更加真實可靠。

    將草地貪夜蛾的潛在分布與其寄主玉米種植整合,可以加深對二者空間分布的整體認識,有利于對其進行風(fēng)險分析和綜合防治。從全球范圍來講,玉米的密集種植區(qū)完全被草地貪夜蛾的潛在分布所覆蓋,這些地區(qū)的玉米種植受到草地貪夜蛾的威脅,需重點防控,非玉米種植區(qū)應(yīng)監(jiān)測草地貪夜蛾對非玉米類農(nóng)作物的侵害。在我國,華北平原春玉米區(qū)和黃淮海平原夏播玉米區(qū)屬于草地貪夜蛾的高風(fēng)險區(qū),其次為東北三省和西南密集種植區(qū),南方丘陵玉米區(qū)和西北灌溉玉米區(qū)受草地貪夜蛾的影響相對較小,根據(jù)受草地貪夜蛾威脅程度不同而有區(qū)別地進行防控。我國東南和華南地區(qū)是草地貪夜蛾高適生區(qū),這些非玉米種植區(qū)應(yīng)監(jiān)控草地貪夜蛾對非玉米類農(nóng)作物的影響。

    將預(yù)測模型與草地貪夜蛾實際發(fā)生狀況進行比較,模型預(yù)測結(jié)果很好的模擬了草地貪夜蛾的潛在分布,與實際發(fā)生狀況相吻合。針對中國來說,在西北灌溉區(qū)(如新疆、西藏等)實際少有發(fā)生,預(yù)測表現(xiàn)出重疊性降低,風(fēng)險較??;在華北平原種植區(qū)(如河北)、東北三省種植區(qū)表現(xiàn)重疊較高,實際發(fā)生量較大;部分地區(qū)雖氣候適合但不是玉米種植區(qū),不會造成草地貪夜蛾的廣泛發(fā)生。

    猜你喜歡
    種植區(qū)夜蛾草地
    C市主要草莓種植區(qū)土壤重金屬鎘、鉛現(xiàn)狀調(diào)查
    不同種植區(qū)隴東苜蓿營養(yǎng)價值的比較研究
    中國飼料(2021年17期)2021-11-02 08:15:24
    侯馬市 采取果斷措施開展對草地貪夜蛾統(tǒng)防統(tǒng)治
    科學(xué)認知草地貪夜蛾 打贏防控攻堅戰(zhàn)
    草地上的事
    幼兒100(2020年31期)2020-11-18 03:42:00
    草地貪夜蛾的識別與防控
    草莓種植區(qū)土壤中典型有機氮化合物的分布及來源
    草地
    小太陽畫報(2018年6期)2018-05-14 17:19:28
    草地上
    河北昌黎縣葡萄種植區(qū)農(nóng)業(yè)地球化學(xué)特征
    亚洲国产最新在线播放| 成人漫画全彩无遮挡| 亚洲久久久国产精品| 少妇的逼水好多| av视频免费观看在线观看| 国产午夜精品一二区理论片| 国产精品国产av在线观看| 亚洲精品美女久久av网站| 中文字幕另类日韩欧美亚洲嫩草| 午夜免费观看性视频| 国产成人精品一,二区| 久久久久久久大尺度免费视频| 高清不卡的av网站| 精品午夜福利在线看| 国产成人91sexporn| 最黄视频免费看| 国产黄频视频在线观看| 丝袜美腿诱惑在线| 深夜精品福利| 亚洲精品美女久久久久99蜜臀 | 一本—道久久a久久精品蜜桃钙片| 国产精品秋霞免费鲁丝片| 久久久久精品久久久久真实原创| 最近最新中文字幕免费大全7| 亚洲av男天堂| 一边亲一边摸免费视频| 日韩三级伦理在线观看| 最新中文字幕久久久久| 最近最新中文字幕免费大全7| 黑人猛操日本美女一级片| 人人澡人人妻人| 久久精品国产a三级三级三级| 久久精品久久久久久噜噜老黄| 国产黄频视频在线观看| 丝袜人妻中文字幕| 电影成人av| 国产精品无大码| 亚洲精品自拍成人| 老司机亚洲免费影院| www.精华液| 欧美精品国产亚洲| 久久午夜综合久久蜜桃| 久久久久久久精品精品| 美女xxoo啪啪120秒动态图| 久久韩国三级中文字幕| 少妇人妻精品综合一区二区| 精品国产一区二区三区四区第35| 午夜福利网站1000一区二区三区| 少妇人妻久久综合中文| 2021少妇久久久久久久久久久| 国产野战对白在线观看| 免费观看av网站的网址| 一级片'在线观看视频| 婷婷色综合www| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲欧美一区二区av| 国产欧美日韩综合在线一区二区| 久久精品久久精品一区二区三区| 乱人伦中国视频| 欧美激情极品国产一区二区三区| 搡女人真爽免费视频火全软件| 观看美女的网站| 精品少妇一区二区三区视频日本电影 | 又粗又硬又长又爽又黄的视频| 一区二区三区乱码不卡18| 18禁动态无遮挡网站| 亚洲精品美女久久久久99蜜臀 | 精品一品国产午夜福利视频| 三级国产精品片| 在线亚洲精品国产二区图片欧美| 韩国高清视频一区二区三区| 亚洲国产成人一精品久久久| 亚洲美女搞黄在线观看| 韩国精品一区二区三区| 三上悠亚av全集在线观看| 老熟女久久久| 美女中出高潮动态图| 亚洲视频免费观看视频| 欧美日本中文国产一区发布| 亚洲av电影在线进入| 欧美xxⅹ黑人| 成人黄色视频免费在线看| 女人精品久久久久毛片| 超碰成人久久| 日本91视频免费播放| 午夜福利在线观看免费完整高清在| 精品久久久久久电影网| 久久国产亚洲av麻豆专区| 久久久久久人妻| 岛国毛片在线播放| 少妇人妻久久综合中文| 2021少妇久久久久久久久久久| 国产精品久久久av美女十八| 久久久亚洲精品成人影院| www.熟女人妻精品国产| 免费看av在线观看网站| 夫妻性生交免费视频一级片| 国产成人精品婷婷| 秋霞在线观看毛片| 国产成人精品福利久久| 亚洲激情五月婷婷啪啪| 一级毛片黄色毛片免费观看视频| 校园人妻丝袜中文字幕| 中文字幕av电影在线播放| 国产免费一区二区三区四区乱码| 99久久精品国产国产毛片| 国产精品99久久99久久久不卡 | 亚洲四区av| 久久久久久久亚洲中文字幕| 天天操日日干夜夜撸| 男女边吃奶边做爰视频| 国产一级毛片在线| 色网站视频免费| 爱豆传媒免费全集在线观看| 亚洲欧美一区二区三区国产| 国产精品二区激情视频| 免费在线观看黄色视频的| 亚洲 欧美一区二区三区| 欧美日韩亚洲高清精品| 成人免费观看视频高清| 国产成人免费观看mmmm| 国产精品99久久99久久久不卡 | 久久精品亚洲av国产电影网| 人人澡人人妻人| 99久国产av精品国产电影| 国产片特级美女逼逼视频| 男人爽女人下面视频在线观看| 丝袜脚勾引网站| 一区二区三区激情视频| 精品亚洲乱码少妇综合久久| 三上悠亚av全集在线观看| 天天躁狠狠躁夜夜躁狠狠躁| h视频一区二区三区| 一区二区日韩欧美中文字幕| www.av在线官网国产| 久久精品久久久久久久性| 久久国产亚洲av麻豆专区| 在线 av 中文字幕| 超碰97精品在线观看| 国产精品一国产av| 亚洲精品,欧美精品| freevideosex欧美| 国产麻豆69| 日本猛色少妇xxxxx猛交久久| 青青草视频在线视频观看| 一本久久精品| 亚洲国产最新在线播放| 大片免费播放器 马上看| 一级黄片播放器| 国产野战对白在线观看| 男女高潮啪啪啪动态图| 欧美精品一区二区免费开放| 国产极品粉嫩免费观看在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲美女视频黄频| 国产麻豆69| 性少妇av在线| 2018国产大陆天天弄谢| 天天躁日日躁夜夜躁夜夜| 婷婷色综合www| 超碰97精品在线观看| 涩涩av久久男人的天堂| 黄色一级大片看看| 在线看a的网站| 亚洲欧美成人综合另类久久久| 国产视频首页在线观看| 嫩草影院入口| av女优亚洲男人天堂| 亚洲精品在线美女| 亚洲国产日韩一区二区| 午夜久久久在线观看| 老司机影院成人| 亚洲精品国产av蜜桃| 波野结衣二区三区在线| 在线观看免费日韩欧美大片| 亚洲色图综合在线观看| 赤兔流量卡办理| 五月伊人婷婷丁香| 中文字幕精品免费在线观看视频| 国产毛片在线视频| 十分钟在线观看高清视频www| 亚洲国产精品成人久久小说| 久久精品人人爽人人爽视色| 欧美中文综合在线视频| 国产欧美日韩一区二区三区在线| 国产成人精品久久久久久| 国产精品一区二区在线不卡| 婷婷成人精品国产| 免费女性裸体啪啪无遮挡网站| www.熟女人妻精品国产| 极品少妇高潮喷水抽搐| 日产精品乱码卡一卡2卡三| 一二三四中文在线观看免费高清| 国产精品香港三级国产av潘金莲 | 日韩中文字幕视频在线看片| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的| 国产av一区二区精品久久| 久久99一区二区三区| 亚洲av男天堂| av不卡在线播放| 美女午夜性视频免费| 人人妻人人澡人人爽人人夜夜| 免费高清在线观看视频在线观看| 人人澡人人妻人| 欧美亚洲日本最大视频资源| 蜜桃在线观看..| 亚洲熟女精品中文字幕| 亚洲情色 制服丝袜| 国产精品嫩草影院av在线观看| 国产爽快片一区二区三区| 亚洲国产日韩一区二区| 26uuu在线亚洲综合色| 高清不卡的av网站| 91成人精品电影| 欧美人与善性xxx| 亚洲色图 男人天堂 中文字幕| 亚洲,一卡二卡三卡| 亚洲国产欧美日韩在线播放| 久久久久精品人妻al黑| 男女下面插进去视频免费观看| 夫妻午夜视频| 亚洲五月色婷婷综合| 国产免费又黄又爽又色| 国产免费现黄频在线看| 制服诱惑二区| 久久97久久精品| 老汉色∧v一级毛片| 如何舔出高潮| 亚洲精品久久成人aⅴ小说| 夫妻性生交免费视频一级片| 视频区图区小说| 啦啦啦中文免费视频观看日本| 777久久人妻少妇嫩草av网站| 亚洲少妇的诱惑av| 久久久欧美国产精品| 寂寞人妻少妇视频99o| 亚洲av.av天堂| 午夜福利影视在线免费观看| 欧美日韩亚洲高清精品| 午夜91福利影院| 人人澡人人妻人| 在线免费观看不下载黄p国产| 国产精品久久久av美女十八| 精品酒店卫生间| 一本—道久久a久久精品蜜桃钙片| 亚洲综合精品二区| 巨乳人妻的诱惑在线观看| 久久久国产精品麻豆| 黄片无遮挡物在线观看| 亚洲欧洲精品一区二区精品久久久 | 国产野战对白在线观看| 美女脱内裤让男人舔精品视频| 亚洲av.av天堂| 久久午夜福利片| 少妇的丰满在线观看| 嫩草影院入口| 欧美亚洲 丝袜 人妻 在线| 搡女人真爽免费视频火全软件| 国产成人免费无遮挡视频| 女人高潮潮喷娇喘18禁视频| 免费大片黄手机在线观看| 免费在线观看黄色视频的| 亚洲综合色网址| a 毛片基地| 少妇的逼水好多| 青草久久国产| 青春草亚洲视频在线观看| 一二三四中文在线观看免费高清| 亚洲四区av| 国产精品久久久久久久久免| 日韩视频在线欧美| 国产黄频视频在线观看| 国产成人av激情在线播放| 18禁动态无遮挡网站| 久久鲁丝午夜福利片| 婷婷色综合大香蕉| av女优亚洲男人天堂| 亚洲第一青青草原| 日韩av在线免费看完整版不卡| 久久人人爽av亚洲精品天堂| 亚洲欧美一区二区三区国产| 少妇的逼水好多| 婷婷色麻豆天堂久久| 免费久久久久久久精品成人欧美视频| 麻豆乱淫一区二区| 精品视频人人做人人爽| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片 在线播放| 国产成人a∨麻豆精品| 日韩制服丝袜自拍偷拍| 久久久久久久久久人人人人人人| 国产成人精品婷婷| 亚洲成人一二三区av| 99国产精品免费福利视频| 久久免费观看电影| 国产亚洲最大av| 日韩一区二区三区影片| 国产一区二区激情短视频 | 国产亚洲欧美精品永久| 人妻少妇偷人精品九色| 国产又色又爽无遮挡免| 国产精品秋霞免费鲁丝片| 在线 av 中文字幕| 日本欧美视频一区| 纵有疾风起免费观看全集完整版| 欧美日韩亚洲国产一区二区在线观看 | 国产精品亚洲av一区麻豆 | 啦啦啦在线免费观看视频4| videos熟女内射| av.在线天堂| 两个人看的免费小视频| 亚洲精品国产av蜜桃| 国产成人午夜福利电影在线观看| 欧美亚洲 丝袜 人妻 在线| 国产片内射在线| 国产精品嫩草影院av在线观看| 热99久久久久精品小说推荐| 日韩人妻精品一区2区三区| 成人漫画全彩无遮挡| 亚洲av综合色区一区| 国产成人午夜福利电影在线观看| 亚洲人成电影观看| 国产男女内射视频| 97人妻天天添夜夜摸| 少妇被粗大的猛进出69影院| 又大又黄又爽视频免费| 曰老女人黄片| 欧美日本中文国产一区发布| 高清av免费在线| 亚洲精品国产av蜜桃| 日韩一本色道免费dvd| 美女xxoo啪啪120秒动态图| 欧美 日韩 精品 国产| 看十八女毛片水多多多| 欧美精品国产亚洲| 五月天丁香电影| 国产日韩欧美在线精品| 老司机影院毛片| 男人操女人黄网站| 亚洲av电影在线观看一区二区三区| 自线自在国产av| 久久久亚洲精品成人影院| 中国国产av一级| 丰满饥渴人妻一区二区三| 日韩在线高清观看一区二区三区| 精品福利永久在线观看| 国产激情久久老熟女| 久久精品夜色国产| 国产精品二区激情视频| 免费播放大片免费观看视频在线观看| 亚洲成人手机| 韩国精品一区二区三区| 久久国产精品大桥未久av| 七月丁香在线播放| 嫩草影院入口| 久久精品国产亚洲av天美| 国产成人精品在线电影| 精品一区二区三区四区五区乱码 | 国产 精品1| 新久久久久国产一级毛片| 午夜福利乱码中文字幕| av在线观看视频网站免费| 欧美日韩视频高清一区二区三区二| 精品国产国语对白av| 麻豆精品久久久久久蜜桃| 在线观看免费视频网站a站| 国产在线免费精品| 婷婷色综合大香蕉| 国产精品免费大片| 婷婷成人精品国产| 欧美97在线视频| 欧美精品人与动牲交sv欧美| 最黄视频免费看| 亚洲精品国产色婷婷电影| 一本久久精品| 18禁观看日本| 成年动漫av网址| 久久人妻熟女aⅴ| 亚洲精品乱久久久久久| 丰满迷人的少妇在线观看| 午夜福利影视在线免费观看| 亚洲av综合色区一区| 亚洲精品久久成人aⅴ小说| 免费高清在线观看日韩| 亚洲美女黄色视频免费看| 又粗又硬又长又爽又黄的视频| 久久国产亚洲av麻豆专区| 热re99久久精品国产66热6| 男女无遮挡免费网站观看| 欧美日韩一区二区视频在线观看视频在线| 午夜av观看不卡| 午夜91福利影院| 国产精品蜜桃在线观看| 老司机影院毛片| 欧美人与性动交α欧美精品济南到 | 欧美日韩一区二区视频在线观看视频在线| 最新中文字幕久久久久| 制服人妻中文乱码| 97人妻天天添夜夜摸| 老女人水多毛片| 老鸭窝网址在线观看| 亚洲熟女精品中文字幕| 丝袜美足系列| 丝袜在线中文字幕| 99国产综合亚洲精品| 欧美日韩视频精品一区| 校园人妻丝袜中文字幕| 免费黄色在线免费观看| 午夜免费观看性视频| 亚洲欧美精品自产自拍| a级毛片在线看网站| 日韩一本色道免费dvd| 一级片免费观看大全| 1024视频免费在线观看| 汤姆久久久久久久影院中文字幕| 日韩av不卡免费在线播放| 成人毛片a级毛片在线播放| 在线看a的网站| 人成视频在线观看免费观看| 日韩一本色道免费dvd| 电影成人av| 久久99蜜桃精品久久| 欧美日韩视频精品一区| 搡女人真爽免费视频火全软件| 国产国语露脸激情在线看| av免费在线看不卡| 女人精品久久久久毛片| 亚洲精品久久午夜乱码| xxxhd国产人妻xxx| www日本在线高清视频| 国产精品嫩草影院av在线观看| 丰满饥渴人妻一区二区三| a级片在线免费高清观看视频| 天堂中文最新版在线下载| 日韩一本色道免费dvd| 久久精品亚洲av国产电影网| 有码 亚洲区| 亚洲色图综合在线观看| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看| 久久免费观看电影| 黑人猛操日本美女一级片| 亚洲国产色片| 午夜免费观看性视频| 人人妻人人澡人人看| 一区二区三区四区激情视频| 999精品在线视频| 色网站视频免费| 少妇人妻久久综合中文| 日韩欧美精品免费久久| av线在线观看网站| 伦理电影大哥的女人| 色婷婷av一区二区三区视频| 精品国产国语对白av| 精品人妻在线不人妻| 国产精品 国内视频| 久久久久久久久免费视频了| 在线观看www视频免费| 中文字幕人妻丝袜一区二区 | 久久精品亚洲av国产电影网| 丰满少妇做爰视频| 免费在线观看完整版高清| 亚洲国产av影院在线观看| 国产亚洲一区二区精品| 人人妻人人澡人人爽人人夜夜| 日本猛色少妇xxxxx猛交久久| 18禁动态无遮挡网站| 午夜福利网站1000一区二区三区| 91在线精品国自产拍蜜月| 久久久a久久爽久久v久久| 国产 一区精品| 99久国产av精品国产电影| 高清视频免费观看一区二区| 久久久久视频综合| 性色av一级| 麻豆av在线久日| 亚洲四区av| 99re6热这里在线精品视频| 亚洲av电影在线观看一区二区三区| 久久精品熟女亚洲av麻豆精品| 久久精品久久久久久噜噜老黄| 午夜激情av网站| 国产成人精品久久久久久| 精品久久久精品久久久| 亚洲国产日韩一区二区| 亚洲国产精品一区三区| 精品人妻在线不人妻| 一本—道久久a久久精品蜜桃钙片| 国产精品蜜桃在线观看| freevideosex欧美| 国产精品久久久久久精品电影小说| videosex国产| xxxhd国产人妻xxx| 日韩中文字幕视频在线看片| 亚洲国产欧美日韩在线播放| 18在线观看网站| 不卡av一区二区三区| 国产一区二区三区综合在线观看| 欧美少妇被猛烈插入视频| 亚洲伊人色综图| 日韩一本色道免费dvd| 丝袜美足系列| 国产精品国产av在线观看| 中国三级夫妇交换| 婷婷成人精品国产| 日产精品乱码卡一卡2卡三| 黄网站色视频无遮挡免费观看| 一级片免费观看大全| 人妻少妇偷人精品九色| 亚洲av日韩在线播放| av不卡在线播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品国产av成人精品| 色网站视频免费| 大片免费播放器 马上看| 女人久久www免费人成看片| av在线播放精品| 欧美日本中文国产一区发布| 日日爽夜夜爽网站| 捣出白浆h1v1| 观看美女的网站| 99re6热这里在线精品视频| 在线观看国产h片| 国产精品麻豆人妻色哟哟久久| 亚洲av成人精品一二三区| 热re99久久精品国产66热6| 电影成人av| 亚洲成人手机| 91精品三级在线观看| 国产片内射在线| 青春草视频在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品日韩在线中文字幕| 一级毛片我不卡| 日本黄色日本黄色录像| 亚洲内射少妇av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产av影院在线观看| 黄网站色视频无遮挡免费观看| 大话2 男鬼变身卡| 国产免费视频播放在线视频| 人妻 亚洲 视频| 有码 亚洲区| 欧美97在线视频| 菩萨蛮人人尽说江南好唐韦庄| 国产成人欧美| 久久免费观看电影| 春色校园在线视频观看| 99国产综合亚洲精品| 一二三四中文在线观看免费高清| 宅男免费午夜| 免费av中文字幕在线| 亚洲欧美色中文字幕在线| 菩萨蛮人人尽说江南好唐韦庄| 九九爱精品视频在线观看| 国产熟女午夜一区二区三区| 精品国产超薄肉色丝袜足j| 一级,二级,三级黄色视频| 老司机亚洲免费影院| 欧美激情高清一区二区三区 | 夫妻午夜视频| 欧美在线黄色| 日韩中文字幕视频在线看片| 91成人精品电影| 久久久久精品久久久久真实原创| 天堂俺去俺来也www色官网| 一级片免费观看大全| 在线 av 中文字幕| 狠狠婷婷综合久久久久久88av| 亚洲色图综合在线观看| 深夜精品福利| 色视频在线一区二区三区| 精品少妇黑人巨大在线播放| 夫妻午夜视频| 少妇的逼水好多| 啦啦啦在线观看免费高清www| 少妇精品久久久久久久| 亚洲美女视频黄频| 久久久久久久大尺度免费视频| 18在线观看网站| 欧美日韩精品网址| 欧美xxⅹ黑人| 免费在线观看视频国产中文字幕亚洲 | 五月天丁香电影| 中文字幕人妻丝袜制服| 亚洲精品国产av成人精品| 天堂中文最新版在线下载| 日韩欧美精品免费久久| 国产免费现黄频在线看| 国产1区2区3区精品| 捣出白浆h1v1| 在线观看美女被高潮喷水网站| 宅男免费午夜| 中文字幕最新亚洲高清| av视频免费观看在线观看| 国产男女内射视频| 亚洲人成网站在线观看播放| 国产亚洲午夜精品一区二区久久| 亚洲第一av免费看| 午夜激情久久久久久久| 精品国产一区二区三区四区第35| 亚洲av福利一区| 久久国产精品男人的天堂亚洲| 中文精品一卡2卡3卡4更新| 精品视频人人做人人爽| 又粗又硬又长又爽又黄的视频| 国产精品香港三级国产av潘金莲 | 韩国精品一区二区三区| 中文字幕人妻熟女乱码| 黄色 视频免费看| 亚洲色图 男人天堂 中文字幕| 国产精品久久久久久精品古装| 国产精品国产三级专区第一集| 少妇 在线观看|