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

    基于最大熵模型預(yù)測太白貝母的潛在分布

    2021-10-13 05:04:46鄒萌蔣瑞平李佳倫段淇洋廖海周嘉裕
    中國中醫(yī)藥信息雜志 2021年9期
    關(guān)鍵詞:環(huán)境模型

    鄒萌,蔣瑞平,李佳倫,段淇洋,廖海,周嘉裕

    西南交通大學(xué)生命科學(xué)與工程學(xué)院,四川 成都 610031

    太白貝母Fritillaria taipaiensisP.Y.Li 為百合科貝母屬植物,是川貝母藥材的基原植物之一,以鱗莖入藥[1]。太白貝母含有生物堿成分,具有潤肺止咳、化痰平喘等功效,臨床廣泛應(yīng)用于肺熱燥咳、肺陰 虛等。由于人為大量采挖和生態(tài)環(huán)境的破壞,太白貝母的野生資源受到嚴(yán)重破壞,已瀕臨枯竭[2]。鑒于太白貝母野生資源所面臨的嚴(yán)峻形勢,《陜西省中藥材保護和發(fā)展實施方案(2016-2020)》(陜政辦發(fā)〔2016〕9 號)已將太白貝母列入中藥材重點保護品種目錄。為緩解太白貝母需求與野生資源匱乏的矛盾,對太白貝母開展引種與人工栽培是一種可行的替代方法。

    有關(guān)太白貝母的栽培報道,始載于清光緒《大寧縣志》“貝母,銀廠坪所產(chǎn)為佳”?!按髮幙h”指今重慶市巫溪縣地區(qū),說明從清代起,當(dāng)?shù)匾殉蔀樘棕惸傅脑耘鄥^(qū)。其后,陜西、甘肅、重慶與四川等地也陸續(xù)有太白貝母的野生資源及栽培區(qū)的報道[3],但盲目引種及擴大栽培區(qū)會嚴(yán)重影響太白貝母生產(chǎn)的合理布局,甚至破壞太白貝母的道地性,導(dǎo)致藥材品質(zhì)下降。因此,開展太白貝母生態(tài)環(huán)境適應(yīng)性研究具有重大現(xiàn)實意義。

    本課題組通過2019 年6 月-2020 年9 月野外實地勘查、查詢相關(guān)文獻與網(wǎng)絡(luò)數(shù)據(jù)庫信息,整理出太白貝母在中國的地理分布區(qū)資料,利用最大熵(MaxEnt)模型和地理信息系統(tǒng)ArcGIS 軟件分析影響太白貝母地理分布的主要環(huán)境因子及其影響程度,并以此為依據(jù)分類預(yù)測太白貝母不同程度適生區(qū)劃圖,分析其潛在地理分布格局,以期為太白貝母資源的可持續(xù)利用提供參考。

    1 數(shù)據(jù)來源與方法

    1.1 太白貝母分布數(shù)據(jù)

    太白貝母的分布數(shù)據(jù)主要來源于3 個部分:①實驗室人員野外實地考察,以及中國數(shù)字植物標(biāo)本館(http://www.cvh.org.cv/)標(biāo)本植物的分布地經(jīng)緯度信息;②通過查閱文獻中的分布區(qū)數(shù)據(jù)整理得到經(jīng)緯度信息,其中對僅有物種名及物種詳細(xì)分布地名稱的,通過地球在線(https://www.earthol.com/)查找具體分布點并記錄;③通過全球生物多樣性信息平臺(https://www.gbif.org/)整理太白貝母在中國范圍內(nèi)分布點,并找到相應(yīng)的經(jīng)緯度信息。最終整理得到中國境內(nèi)115 份包含20 個地點的分布點數(shù)據(jù),將各分布點數(shù)據(jù)統(tǒng)一轉(zhuǎn)換為十進制并以.csv 格式保存。詳見表1。

    表1 太白貝母分布數(shù)據(jù)

    1.2 環(huán)境變量

    從世界氣象數(shù)據(jù)(https://www.worldclim.org/)下載32 個環(huán)境變量因子,涵蓋19 個基礎(chǔ)生物氣候變量(bio1~19),以及1、5、6、7 月的最低溫度(tmin1、tmin5、tmin6、tmin7)、最高溫度(tmax1、tmax5、tmax6、tmax7)、平均溫度(tavg1、tavg5、tavg6、tavg7)和海拔(alt)。采用ChinaMap 中4p 格式的1∶400 萬中華人民共和國地圖和中華人民共和國行政區(qū)劃圖(https://www.csdn.net/)作為底圖得到太白貝母在中國境內(nèi)的分布點。詳見圖1。

    圖1 太白貝母分布點

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

    1.3.1 環(huán)境因子選取

    由于環(huán)境因子間有一定的相關(guān)性,所以本研究對32 個環(huán)境因子進行相關(guān)性檢驗后用于MaxEnt 模型。利用ArcGIS12.0 軟件提取各采樣點的點插值,以SPSS20.0 軟件Spearman 相關(guān)系數(shù)查驗相關(guān)性,選取相關(guān)性<0.8 的環(huán)境因子,若2 個因子相關(guān)性≥0.8,只選擇任意1 個因子運行模型。

    1.3.2 最大熵模型預(yù)測

    運行MaxEnt3.4.1 軟件,在樣本欄(samples)添加上述.csv 文件,在環(huán)境圖層欄(environmental layers)添加所選取的低相關(guān)性環(huán)境因子,采用刀切法(jackknife)預(yù)測各環(huán)境因子的貢獻率,制作響應(yīng)曲線(response curve)并生成初步適生區(qū)預(yù)測圖像。預(yù)測模型的準(zhǔn)確性是利用受試者工作特征(ROC)曲線下方面積(AUC)值判定,其值越接近1,說明與隨機分布相距越遠(yuǎn),且環(huán)境因子與預(yù)測的物種地理分布之間的相關(guān)性越大,即MaxEnt 模型預(yù)測的結(jié)果越準(zhǔn)確。

    本研究采用Roberto、Natural breaks 等研究方法,選取25%分布點作為測試集與75%分布點作為訓(xùn)練集,模型運算次數(shù)設(shè)定為500 次,重復(fù)運算設(shè)定為10 次,其余參數(shù)均為默認(rèn)模型的自動設(shè)置。

    1.4 適應(yīng)性分布區(qū)劃分

    將MaxEnt 模型運行得到的結(jié)果導(dǎo)入ArcGIS12.0軟件,利用重分類進行疊加分析和地圖制作,采用人工分級方法劃分出太白貝母適應(yīng)性分布等級,得到主要生態(tài)因子影響下的適應(yīng)性區(qū)劃圖。根據(jù)MaxEnt 模型運行結(jié)果中各環(huán)境因子的響應(yīng)曲線得出主要因子適宜值范圍及太白貝母適生區(qū)的生態(tài)特征。

    2 結(jié)果與分析

    2.1 預(yù)測的有效性

    近年來,ROC 曲線分析法在預(yù)測物種潛在適生區(qū)模型評價中應(yīng)用廣泛[4]。該方法是對曲線下方面積進行計算,即得到AUC 值,由此判斷模型模擬的精準(zhǔn)度。一般認(rèn)為,AUC<0.6 時判斷價值較低,AUC>0.8 時判斷價值較高,AUC=1 時為理想分布[5]。本研究在對各環(huán)境因子重復(fù)運行10 次后得到訓(xùn)練數(shù)據(jù)集AUC=0.958,說明MaxEnt 模型對太白貝母潛在分布區(qū)的預(yù)測結(jié)果較好,可信度非常高。見圖2。

    圖2 太白貝母分布MaxEnt 模型ROC 曲線及AUC 值

    2.2 環(huán)境因子貢獻率

    jackknife 是在預(yù)測物種潛在分布的適生區(qū)和分析環(huán)境因子對所預(yù)測結(jié)果的影響程度大小,以及通過運算確定對物種分布貢獻率最大的環(huán)境因子過程中最常用的方法。本研究根據(jù)jackknife 分析得到的預(yù)測結(jié)果顯示,各環(huán)境因子在對太白貝母適應(yīng)度的影響所占的權(quán)重不同。見圖3。

    圖3 環(huán)境因子在MaxEnt 模型中預(yù)測太白貝母分布刀切圖

    通過模型篩選得到對太白貝母潛在分布貢獻率>0 的環(huán)境因子有6 個,即晝夜溫差均值(41.21%)、海拔(26.42%)、年降水量(18.32%)、溫度季節(jié)性變化的標(biāo)準(zhǔn)差(13.43%)、最濕月降水量(0.63%)、6月最低氣溫(0.12%)。其中晝夜溫差均值和海拔的貢獻率最高,溫度、海拔與降水是影響太白貝母分布區(qū)的重要環(huán)境因素,見表2。

    表2 對太白貝母潛在分布貢獻率>0 的6 個環(huán)境因子

    2.3 主要環(huán)境因子

    通過運行MaxEnt 模型可以得到所選取各環(huán)境因子的單變量響應(yīng)曲線,從而分析其最適區(qū)間值。為明確主要環(huán)境因子的參數(shù)值與太白貝母存在概率之間的關(guān)系,本研究在MaxEnt 模型運行結(jié)果基礎(chǔ)上繪制僅包含單一環(huán)境因子的響應(yīng)曲線,其中包含bio2、alt、bio12 與bio4。結(jié)果表明,太白貝母的潛在分布概率隨環(huán)境的改變波動值較大。選擇存在概率>0.2 的區(qū)間作為其適生區(qū)范圍。

    晝夜溫差均值和海拔影響最大,其中當(dāng)晝夜溫差均值約為8.1 ℃時太白貝母適生概率最大,達到66.28%,晝夜溫差2.7~11.9 ℃為其適生區(qū)范圍。當(dāng)海拔為241 m 時達到太白貝母適生范圍;241 m<海拔<1 318 m時呈增長趨勢;1 318 m時到達最適概率,約為65.7%;1 318 m<海拔<3 828 m 時適生概率逐漸降低。因此,在晝夜溫差約為8.1 ℃及海拔約為1 318 m 時,最適宜太白貝母生長。詳見圖4、圖5。

    圖4 太白貝母晝夜溫差均值響應(yīng)曲線

    圖5 太白貝母海拔響應(yīng)曲線

    2.4 太白貝母的適應(yīng)性劃分

    將MaxEnt 模型的預(yù)測結(jié)果導(dǎo)入ArcGIS 軟件中運行,對總體適生區(qū)進行重分類,確定把適生等級分為4 類:0<適應(yīng)指數(shù)≤0.06 為非適生區(qū);0.06<適應(yīng)指數(shù)≤0.23 為低適生區(qū),主要分布于四川、甘肅南部、云南南部、浙江南部,在青海東南部、新疆北部和西藏西部有零星分布;0.23<適應(yīng)指數(shù)≤0.51 為適生區(qū),主要分布于甘肅南部、貴州西北部、四川中部及東北部、重慶南部、陜西南部及湖北西部;0.51<適生指數(shù)<0.86 為高適生區(qū),主要分布于甘肅南部、貴州西北部、四川中部及東北部、重慶南部、陜西南部及湖北西部等地區(qū)。見圖6。潛在分布區(qū)在各?。▍^(qū))的分布面積統(tǒng)計情況見表3。

    圖6 太白貝母潛在分布核心區(qū)域

    表3 太白貝母潛在分布區(qū)域面積

    3 討論

    利用生態(tài)位模型的運行對物種的潛在分布區(qū)進行預(yù)測,進而對物種建立保護區(qū)及規(guī)劃人工栽培地域是目前廣泛應(yīng)用的方法。MaxEnt 模型運行簡單,預(yù)測結(jié)果的準(zhǔn)確性高且可信度較強,已被廣泛用于物種的潛在分布區(qū)預(yù)測及適應(yīng)性評價[6]。

    本研究結(jié)合MaxEnt 模型和ArcGIS,采用的分布點數(shù)據(jù)涵蓋了太白貝母的現(xiàn)有分布區(qū),結(jié)果發(fā)現(xiàn)太白貝母主要集中在陜西、甘肅、重慶及四川等地,與陳士林等[7]關(guān)于太白貝母最大生境區(qū)域的研究結(jié)果一致。

    基于收集分布點的32 個環(huán)境因子數(shù)據(jù),本研究對太白貝母的潛在分布區(qū)進行預(yù)測及適應(yīng)性等級劃分。結(jié)果MaxEnt 模型ROC 曲線訓(xùn)練集AUC=0.958,表示該模型預(yù)測準(zhǔn)確度達到較高水平,可用于太白貝母潛在分布區(qū)的預(yù)測。較之現(xiàn)有分布區(qū)[8],潛在分布區(qū)增加了貴州、湖南及湖北的部分地區(qū),其中貴州尤其值得關(guān)注,不僅與重慶接壤,且其總面積的76%可作為太白貝母的潛在分布區(qū),表明貴州極有可能成為未來太白貝母產(chǎn)業(yè)的重點發(fā)展區(qū)域,需提前開展合理布局與規(guī)劃。潛在分布區(qū)中,高適生區(qū)是太白貝母生長與繁殖的最適區(qū)域,主要分布于甘肅南部、貴州西北部、四川中部和東北部、重慶南部、陜西南部及湖北西部等地區(qū)。這些區(qū)域集中了太白貝母的一些傳統(tǒng)產(chǎn)區(qū),如陜西太白縣晝夜溫差大,氣候較寒冷,一直是太白貝母種植最適宜的地區(qū)[9],印證了本研究預(yù)測結(jié)果的準(zhǔn)確性。太白貝母的潛在核心分布區(qū)與川貝母Fritillaria cirrhosaD.Don、暗紫貝母Fritillaria unibracteataP.K.Hsiao &K.C.Hsia 和甘肅貝母Fritillaria przewalskiiMaxim.ex Batal 等品種的分布區(qū)有一定范圍的重疊,其中在甘肅南部與甘肅貝母Fritillaria przewalskiiMaxim.ex Batal 有重疊,在四川省與暗紫貝母Fritillaria unibracteataP.K.Hsiao &K.C.Hsia 和川貝母Fritillaria cirrhosaD.Don 重疊,在湖北西部和重慶市與湖北貝母Fritillaria hupehensisHsiao et K.C.Hsia 有重疊。這些重疊區(qū)的存在,一方面體現(xiàn)出不同品種貝母所需的生態(tài)環(huán)境因子高度相似,另一方面也表明這些區(qū)域擁有較好的生態(tài)優(yōu)越性,盡管出現(xiàn)了不同貝母的空間生態(tài)位重疊,但彼此間卻不產(chǎn)生生存競爭。

    在影響太白貝母分布的主導(dǎo)環(huán)境因子中,氣候因子對太白貝母的影響大于土壤與地形因子。晝夜溫差均值是影響太白貝母分布的最重要環(huán)境因子,當(dāng)晝夜溫差均值為8.1 ℃時,最有利于太白貝母生長過程中營養(yǎng)物質(zhì)的累積,并促進太白貝母鱗莖的生長。有報道,等溫性是影響川貝母Fritillaria cirrhosaD.Don分布的最重要環(huán)境因子[10]。太白貝母與川貝母Fritillaria cirrhosaD.Don 具有不同的最重要環(huán)境因子,推測其原因可能在于太白貝母屬于川貝母類群中較為獨立的品種。太白貝母的獨立性也體現(xiàn)在DNA上,其ITS2 等條形碼序列與其他品種川貝母有較大差別[11],在進化樹上處于不同分支,表明太白貝母較早從川貝母類群中分化出來,導(dǎo)致其DNA 與所處環(huán)境與其他川貝母品種產(chǎn)生較大差別。由此,我們推測植物的生長環(huán)境適應(yīng)性可能作為植物系統(tǒng)進化研究的一個參數(shù),但有待于進一步驗證。

    本研究顯示,海拔是影響太白貝母分布的第2 個重要因素,最適海拔為1 318 m。研究顯示,不同海拔高度影響藥材品質(zhì),海拔梯度造成的水分、溫度及土壤肥力等環(huán)境差異性會對植物生長發(fā)育產(chǎn)生極大的影響,進而導(dǎo)致植物生理和生態(tài)適應(yīng)性的變化[12],如生長在海拔1 500 m 的重慶產(chǎn)桔梗藥材質(zhì)量最優(yōu)[13]。由于暗紫貝母等其余川貝母品種的最低海拔在3 000 m 左右[14],因此同一地區(qū)不同海拔可種植不同類型的川貝母品種,如低處種植太白貝母,高處種植川貝母,如此則能有效利用生境資源,實現(xiàn)空間利用最大化。

    目前川貝母的野生資源不容樂觀,而栽培技術(shù)已較為成熟,所以對太白貝母的適生區(qū)進行預(yù)測并以此指導(dǎo)進行合理種植與栽培十分必要。相較其他品種,太白貝母有望成為低海拔地區(qū)川貝母栽培的優(yōu)良候選品種。在合理規(guī)劃前提下建立太白貝母中藥材生產(chǎn)質(zhì)量管理規(guī)范(GAP)生產(chǎn)基地,對太白貝母資源的可持續(xù)利用有重要意義。

    猜你喜歡
    環(huán)境模型
    一半模型
    長期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學(xué)習(xí)的虛擬仿真環(huán)境
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    孕期遠(yuǎn)離容易致畸的環(huán)境
    不能改變環(huán)境,那就改變心境
    環(huán)境
    孕期遠(yuǎn)離容易致畸的環(huán)境
    3D打印中的模型分割與打包
    亚洲欧美激情在线| 亚洲国产精品999在线| 婷婷精品国产亚洲av在线| 亚洲伊人色综图| 波多野结衣av一区二区av| 国产又爽黄色视频| 亚洲第一青青草原| 美女高潮喷水抽搐中文字幕| 在线天堂中文资源库| 免费高清在线观看日韩| 国产高清有码在线观看视频 | 91在线观看av| 亚洲色图av天堂| 国产午夜福利久久久久久| 97碰自拍视频| 免费看a级黄色片| 亚洲国产看品久久| 欧美老熟妇乱子伦牲交| 久久久久久久午夜电影| 满18在线观看网站| 给我免费播放毛片高清在线观看| 欧美色视频一区免费| 亚洲免费av在线视频| 国产精品 欧美亚洲| 国产精品电影一区二区三区| 欧美绝顶高潮抽搐喷水| 在线观看免费视频日本深夜| 大型黄色视频在线免费观看| 精品一品国产午夜福利视频| 亚洲,欧美精品.| 亚洲av电影在线进入| 国产高清有码在线观看视频 | 黑人巨大精品欧美一区二区mp4| 国产精品综合久久久久久久免费 | 亚洲av第一区精品v没综合| 免费看十八禁软件| 亚洲国产欧美一区二区综合| 在线观看午夜福利视频| 午夜老司机福利片| 美女扒开内裤让男人捅视频| 50天的宝宝边吃奶边哭怎么回事| 丝袜美足系列| 美女国产高潮福利片在线看| 免费在线观看影片大全网站| 免费少妇av软件| 亚洲男人天堂网一区| 国产av精品麻豆| 十八禁人妻一区二区| 99久久精品国产亚洲精品| 成人18禁在线播放| 中文字幕高清在线视频| 又黄又爽又免费观看的视频| 不卡一级毛片| 美女大奶头视频| 国产欧美日韩一区二区三| 欧美日韩中文字幕国产精品一区二区三区 | 中文字幕另类日韩欧美亚洲嫩草| 久久精品成人免费网站| 亚洲欧美激情综合另类| 国产片内射在线| 午夜日韩欧美国产| 精品少妇一区二区三区视频日本电影| 午夜精品国产一区二区电影| 国产又色又爽无遮挡免费看| 免费高清视频大片| 日本a在线网址| 如日韩欧美国产精品一区二区三区| 国产亚洲欧美98| 精品第一国产精品| 老司机深夜福利视频在线观看| 久久精品91蜜桃| 午夜影院日韩av| 亚洲国产毛片av蜜桃av| 亚洲激情在线av| 99国产综合亚洲精品| www.www免费av| 女警被强在线播放| 日本vs欧美在线观看视频| 久久久久国产一级毛片高清牌| xxx96com| 一本大道久久a久久精品| 亚洲av五月六月丁香网| 老熟妇仑乱视频hdxx| 国产成人系列免费观看| 国产成年人精品一区二区| 99国产综合亚洲精品| 成人手机av| 在线av久久热| 黑人欧美特级aaaaaa片| 女人爽到高潮嗷嗷叫在线视频| 夜夜看夜夜爽夜夜摸| 亚洲自拍偷在线| 国产伦人伦偷精品视频| 国产人伦9x9x在线观看| 大型黄色视频在线免费观看| 亚洲九九香蕉| 成人18禁在线播放| 91字幕亚洲| 日韩视频一区二区在线观看| 欧美中文综合在线视频| 每晚都被弄得嗷嗷叫到高潮| 免费观看人在逋| 中出人妻视频一区二区| 成人特级黄色片久久久久久久| 免费高清视频大片| 精品欧美一区二区三区在线| 日本黄色视频三级网站网址| 国产精品精品国产色婷婷| 欧美中文日本在线观看视频| 可以在线观看的亚洲视频| 丝袜美足系列| 一区二区日韩欧美中文字幕| 99香蕉大伊视频| 色综合婷婷激情| 久久国产乱子伦精品免费另类| 国产高清有码在线观看视频 | 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产在线观看| 国产精品久久久久久亚洲av鲁大| 免费搜索国产男女视频| 国产激情欧美一区二区| 国产真人三级小视频在线观看| 日本 欧美在线| 国产在线精品亚洲第一网站| 国产男靠女视频免费网站| 日本a在线网址| 成人18禁在线播放| 亚洲一区二区三区色噜噜| 成人免费观看视频高清| www.精华液| 亚洲人成网站在线播放欧美日韩| 啦啦啦 在线观看视频| 亚洲人成伊人成综合网2020| 亚洲专区国产一区二区| 亚洲无线在线观看| 欧洲精品卡2卡3卡4卡5卡区| 每晚都被弄得嗷嗷叫到高潮| 国产av精品麻豆| 亚洲avbb在线观看| 亚洲三区欧美一区| 国内毛片毛片毛片毛片毛片| 9色porny在线观看| 国产精品久久久av美女十八| tocl精华| 少妇熟女aⅴ在线视频| 午夜福利一区二区在线看| 制服丝袜大香蕉在线| 一区在线观看完整版| 日韩欧美在线二视频| 精品久久久精品久久久| 亚洲精品中文字幕在线视频| 国产不卡一卡二| 露出奶头的视频| a级毛片在线看网站| 非洲黑人性xxxx精品又粗又长| 久久国产精品人妻蜜桃| 老熟妇乱子伦视频在线观看| 精品久久久久久,| 国产一区二区三区在线臀色熟女| 午夜精品久久久久久毛片777| 99精品在免费线老司机午夜| avwww免费| 真人一进一出gif抽搐免费| 亚洲国产欧美网| 美女免费视频网站| 69精品国产乱码久久久| 欧美国产精品va在线观看不卡| 欧美午夜高清在线| 两人在一起打扑克的视频| av免费在线观看网站| 91精品三级在线观看| av片东京热男人的天堂| 色哟哟哟哟哟哟| 久久人妻福利社区极品人妻图片| 怎么达到女性高潮| 国产主播在线观看一区二区| 国产欧美日韩一区二区精品| 欧美不卡视频在线免费观看 | 色综合亚洲欧美另类图片| 久久香蕉激情| av免费在线观看网站| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| 国产精品免费视频内射| 叶爱在线成人免费视频播放| 国产精品一区二区在线不卡| 97人妻精品一区二区三区麻豆 | 桃色一区二区三区在线观看| 少妇 在线观看| 亚洲熟妇熟女久久| 久久久国产成人精品二区| 久久国产精品影院| 搡老妇女老女人老熟妇| 麻豆一二三区av精品| 男人舔女人下体高潮全视频| x7x7x7水蜜桃| 久久婷婷成人综合色麻豆| 少妇裸体淫交视频免费看高清 | 真人一进一出gif抽搐免费| 国产伦一二天堂av在线观看| 日韩精品免费视频一区二区三区| 久久久久久久精品吃奶| 久久久久国内视频| 性欧美人与动物交配| 香蕉久久夜色| 国产野战对白在线观看| 91精品三级在线观看| 91国产中文字幕| 午夜久久久久精精品| 国语自产精品视频在线第100页| 中文字幕人成人乱码亚洲影| 18禁美女被吸乳视频| 欧美中文日本在线观看视频| 天天一区二区日本电影三级 | 又大又爽又粗| 性欧美人与动物交配| 看黄色毛片网站| 亚洲 国产 在线| 国产欧美日韩精品亚洲av| a在线观看视频网站| 国产av一区在线观看免费| 一二三四在线观看免费中文在| 日本欧美视频一区| 欧美激情高清一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产精品电影一区二区三区| 午夜精品国产一区二区电影| or卡值多少钱| 国产单亲对白刺激| 熟妇人妻久久中文字幕3abv| 午夜免费成人在线视频| 久久久久久国产a免费观看| 99在线视频只有这里精品首页| 午夜久久久久精精品| 黑人巨大精品欧美一区二区mp4| 欧美激情 高清一区二区三区| 久久精品成人免费网站| 高清在线国产一区| 一二三四在线观看免费中文在| 亚洲视频免费观看视频| 欧美成人一区二区免费高清观看 | 97人妻天天添夜夜摸| 亚洲国产精品久久男人天堂| www日本在线高清视频| 国产av一区二区精品久久| 久久这里只有精品19| 变态另类成人亚洲欧美熟女 | 亚洲专区国产一区二区| 满18在线观看网站| 亚洲专区字幕在线| 午夜福利免费观看在线| 一级毛片女人18水好多| 在线免费观看的www视频| 首页视频小说图片口味搜索| 精品电影一区二区在线| 国产亚洲av嫩草精品影院| 国产伦一二天堂av在线观看| 欧美日本中文国产一区发布| 国内久久婷婷六月综合欲色啪| 18禁观看日本| 99riav亚洲国产免费| 亚洲午夜理论影院| 在线免费观看的www视频| 国产单亲对白刺激| 日韩精品青青久久久久久| 久久精品成人免费网站| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲天堂国产精品一区在线| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 9191精品国产免费久久| 在线天堂中文资源库| √禁漫天堂资源中文www| 18禁观看日本| 亚洲色图 男人天堂 中文字幕| 99国产精品一区二区三区| 国产欧美日韩精品亚洲av| 亚洲男人的天堂狠狠| 此物有八面人人有两片| 免费搜索国产男女视频| www.熟女人妻精品国产| 久久国产精品影院| 欧美在线黄色| 大型av网站在线播放| 国产av一区在线观看免费| 不卡av一区二区三区| 久久久国产欧美日韩av| 国产亚洲精品久久久久5区| 性欧美人与动物交配| 少妇熟女aⅴ在线视频| 日韩欧美一区视频在线观看| 一级片免费观看大全| 黄色丝袜av网址大全| 黄色成人免费大全| 国产私拍福利视频在线观看| 久久久久九九精品影院| 久久人人爽av亚洲精品天堂| 亚洲 欧美 日韩 在线 免费| 欧美激情 高清一区二区三区| 两性夫妻黄色片| 中文字幕久久专区| 此物有八面人人有两片| 免费看美女性在线毛片视频| 亚洲欧美精品综合久久99| 在线十欧美十亚洲十日本专区| 久久国产乱子伦精品免费另类| 免费女性裸体啪啪无遮挡网站| 午夜福利成人在线免费观看| 午夜激情av网站| 91麻豆精品激情在线观看国产| 91国产中文字幕| 成人欧美大片| 村上凉子中文字幕在线| 日韩欧美免费精品| 极品人妻少妇av视频| 一区二区三区国产精品乱码| 国产高清视频在线播放一区| or卡值多少钱| 可以免费在线观看a视频的电影网站| 在线观看免费视频日本深夜| 一区福利在线观看| 一个人观看的视频www高清免费观看 | 涩涩av久久男人的天堂| 香蕉丝袜av| 亚洲avbb在线观看| 欧美亚洲日本最大视频资源| 国产三级黄色录像| 欧美一级毛片孕妇| 亚洲电影在线观看av| 色综合站精品国产| 久久人人精品亚洲av| 大陆偷拍与自拍| 中出人妻视频一区二区| 女同久久另类99精品国产91| av有码第一页| 午夜福利成人在线免费观看| 老熟妇仑乱视频hdxx| 黄频高清免费视频| 国内精品久久久久精免费| 日韩成人在线观看一区二区三区| 亚洲精品久久国产高清桃花| 亚洲国产欧美网| 日本精品一区二区三区蜜桃| 桃色一区二区三区在线观看| 久久精品91蜜桃| 免费女性裸体啪啪无遮挡网站| 我的亚洲天堂| 免费av毛片视频| 后天国语完整版免费观看| 在线观看66精品国产| 国产av精品麻豆| 亚洲欧美日韩高清在线视频| 91成人精品电影| 69精品国产乱码久久久| 久久香蕉精品热| 久久久久久久精品吃奶| 男男h啪啪无遮挡| 精品乱码久久久久久99久播| 午夜精品久久久久久毛片777| 亚洲一区二区三区不卡视频| 男女下面进入的视频免费午夜 | 一级片免费观看大全| 亚洲男人的天堂狠狠| 韩国精品一区二区三区| 欧美成人性av电影在线观看| 国产成人免费无遮挡视频| 国产在线观看jvid| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 香蕉丝袜av| 丝袜美足系列| 亚洲精品美女久久av网站| 国产精品一区二区在线不卡| 国产成人系列免费观看| 精品无人区乱码1区二区| 亚洲熟女毛片儿| 亚洲精品中文字幕一二三四区| 男男h啪啪无遮挡| 精品乱码久久久久久99久播| www.精华液| 乱人伦中国视频| 真人一进一出gif抽搐免费| 美女 人体艺术 gogo| 咕卡用的链子| 国产精品日韩av在线免费观看 | 国产亚洲精品第一综合不卡| 99riav亚洲国产免费| 不卡av一区二区三区| 91成年电影在线观看| 精品久久久久久久久久免费视频| 亚洲自拍偷在线| 人人妻人人澡人人看| 老汉色∧v一级毛片| 18美女黄网站色大片免费观看| 搞女人的毛片| 90打野战视频偷拍视频| 人人澡人人妻人| 男女下面进入的视频免费午夜 | 国产精品久久视频播放| 咕卡用的链子| 在线天堂中文资源库| 欧美国产日韩亚洲一区| 亚洲成人精品中文字幕电影| 久久久国产成人免费| 少妇 在线观看| 男女之事视频高清在线观看| 欧美最黄视频在线播放免费| 在线永久观看黄色视频| 亚洲一区高清亚洲精品| 午夜视频精品福利| 免费看十八禁软件| 两性午夜刺激爽爽歪歪视频在线观看 | www.精华液| 日韩欧美国产一区二区入口| 日韩大码丰满熟妇| www日本在线高清视频| 免费高清视频大片| 夜夜夜夜夜久久久久| 午夜激情av网站| 久久精品人人爽人人爽视色| 淫妇啪啪啪对白视频| 国产高清激情床上av| 亚洲一区高清亚洲精品| 乱人伦中国视频| 免费在线观看黄色视频的| 亚洲人成网站在线播放欧美日韩| 日本 欧美在线| 国语自产精品视频在线第100页| 久久久久久大精品| 成人欧美大片| 国产精品二区激情视频| 丰满人妻熟妇乱又伦精品不卡| 国产又爽黄色视频| 欧美黄色淫秽网站| 曰老女人黄片| 亚洲精品美女久久av网站| 高清在线国产一区| 99re在线观看精品视频| aaaaa片日本免费| 每晚都被弄得嗷嗷叫到高潮| a级毛片在线看网站| 亚洲va日本ⅴa欧美va伊人久久| 亚洲熟妇中文字幕五十中出| 制服人妻中文乱码| 国产精品久久久av美女十八| 国产国语露脸激情在线看| 久久婷婷成人综合色麻豆| 久久精品成人免费网站| 搞女人的毛片| 国产精品自产拍在线观看55亚洲| 麻豆一二三区av精品| 国产一区二区三区综合在线观看| 亚洲avbb在线观看| 欧美精品啪啪一区二区三区| 欧美精品啪啪一区二区三区| 琪琪午夜伦伦电影理论片6080| av天堂久久9| 婷婷六月久久综合丁香| 夜夜夜夜夜久久久久| 此物有八面人人有两片| 热99re8久久精品国产| 悠悠久久av| 给我免费播放毛片高清在线观看| 亚洲三区欧美一区| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 亚洲性夜色夜夜综合| 中文字幕久久专区| 色婷婷久久久亚洲欧美| 国产精品av久久久久免费| 免费一级毛片在线播放高清视频 | 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 国产精品98久久久久久宅男小说| 性色av乱码一区二区三区2| 精品高清国产在线一区| 国产精品,欧美在线| www.精华液| 国产又色又爽无遮挡免费看| 欧美丝袜亚洲另类 | 国产成人免费无遮挡视频| 亚洲av美国av| 午夜成年电影在线免费观看| 巨乳人妻的诱惑在线观看| 久久久国产成人精品二区| 久久人人97超碰香蕉20202| 免费看十八禁软件| 一边摸一边抽搐一进一出视频| 一夜夜www| 99国产精品一区二区三区| 久久中文字幕一级| 免费一级毛片在线播放高清视频 | 人妻丰满熟妇av一区二区三区| 长腿黑丝高跟| 在线国产一区二区在线| 精品国产乱码久久久久久男人| 一本大道久久a久久精品| 亚洲九九香蕉| 免费在线观看视频国产中文字幕亚洲| 免费女性裸体啪啪无遮挡网站| 中文字幕精品免费在线观看视频| 中文字幕另类日韩欧美亚洲嫩草| 看片在线看免费视频| 久久久久久久午夜电影| 韩国av一区二区三区四区| 91成年电影在线观看| 精品国产国语对白av| 国产一区二区在线av高清观看| 黄色毛片三级朝国网站| 淫秽高清视频在线观看| 99久久99久久久精品蜜桃| 极品教师在线免费播放| 黑人操中国人逼视频| 国产精品av久久久久免费| av电影中文网址| 性色av乱码一区二区三区2| 久久精品亚洲熟妇少妇任你| 18禁美女被吸乳视频| 女性被躁到高潮视频| 国产精品久久久久久人妻精品电影| 岛国在线观看网站| 人成视频在线观看免费观看| 国语自产精品视频在线第100页| 久久天躁狠狠躁夜夜2o2o| 欧美成人午夜精品| 国产成人影院久久av| 精品一区二区三区视频在线观看免费| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 不卡一级毛片| 男人操女人黄网站| 午夜视频精品福利| 午夜老司机福利片| 涩涩av久久男人的天堂| 中文字幕色久视频| 久久狼人影院| 日韩欧美在线二视频| 午夜免费观看网址| 午夜免费鲁丝| АⅤ资源中文在线天堂| 日日夜夜操网爽| 午夜福利免费观看在线| 亚洲一区二区三区不卡视频| 狠狠狠狠99中文字幕| 亚洲午夜精品一区,二区,三区| a级毛片在线看网站| 久久香蕉精品热| 男人操女人黄网站| 999久久久国产精品视频| x7x7x7水蜜桃| 在线视频色国产色| 亚洲av电影不卡..在线观看| av视频免费观看在线观看| 久久热在线av| 女生性感内裤真人,穿戴方法视频| 热re99久久国产66热| 日韩欧美免费精品| 欧美激情久久久久久爽电影 | 757午夜福利合集在线观看| 夜夜爽天天搞| 国产亚洲精品久久久久久毛片| tocl精华| 乱人伦中国视频| 狂野欧美激情性xxxx| 黄色女人牲交| 一级毛片女人18水好多| 天天添夜夜摸| 激情在线观看视频在线高清| 这个男人来自地球电影免费观看| 不卡一级毛片| 精品一区二区三区四区五区乱码| 午夜福利18| 成年女人毛片免费观看观看9| 中文字幕人成人乱码亚洲影| 久久久久精品国产欧美久久久| 十八禁网站免费在线| 热re99久久国产66热| 91av网站免费观看| 在线十欧美十亚洲十日本专区| 国产精品久久视频播放| 亚洲 欧美 日韩 在线 免费| 国产国语露脸激情在线看| 亚洲精品一区av在线观看| 91麻豆av在线| 免费av毛片视频| 日日爽夜夜爽网站| 黄色视频,在线免费观看| 老熟妇乱子伦视频在线观看| www.www免费av| 好男人电影高清在线观看| 欧美丝袜亚洲另类 | 国产成人系列免费观看| 黄色丝袜av网址大全| 18禁国产床啪视频网站| 看免费av毛片| 成人亚洲精品一区在线观看| 岛国在线观看网站| 在线观看午夜福利视频| 精品电影一区二区在线| 操美女的视频在线观看| 亚洲成人精品中文字幕电影| 日本在线视频免费播放| 99热只有精品国产| 91精品国产国语对白视频| 久久久国产精品麻豆| 欧美性长视频在线观看| 中国美女看黄片| 精品电影一区二区在线| 欧美中文日本在线观看视频| 日本 av在线| 亚洲精品国产色婷婷电影| 国产91精品成人一区二区三区| 怎么达到女性高潮| 久久精品人人爽人人爽视色|