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

    基于Maxent的椰子致死性黃化植原體在中國的適生性分析

    2014-04-29 15:47:39曹學(xué)仁等
    熱帶作物學(xué)報 2014年11期
    關(guān)鍵詞:分布預(yù)測

    曹學(xué)仁等

    摘 要 利用Maxent生態(tài)位模型和ArcGIS軟件,并結(jié)合椰子樹(以下簡稱椰樹)在中國的種植分布,對椰子致死性黃化植原體在中國及中國椰樹種植區(qū)的適生性進行了預(yù)測。結(jié)果表明:中國主要的椰樹種植區(qū)均為椰子致死性黃化病的潛在適生區(qū),其中華南沿海地區(qū)及海南大部分地區(qū)屬于中風險區(qū)和高風險區(qū)。經(jīng)ROC曲線分析法驗證,Maxent模型的AUC值為0.992,表明預(yù)測結(jié)果具有較高的可信度。

    關(guān)鍵詞 椰子致死性黃化植原體;Maxent;分布;預(yù)測

    中圖分類號 S432.1 文獻標識碼 A

    Abstract The potential distribution regions of coconut lethal yellowing phytoplasma in China and in coconut growth regions of China were predicted by using Maxent model and ArcGIS,combined with the distribution of coconut in China. The results showed that most coconut growth regions in China were the potential distribution regions for the pathogen. Coastal areas of south China and most regions in Hainan were at moderate to high risk of the pathogen establishment. The checked result of receiver operating characteristic(ROC)curve analysis showed that the area under curve(AUC)value of Maxent model reached 0.992,indicating that the predicted results was fairly satisfactory.

    Key words Coconut lethal yellowing phytoplasma;Maxent;Distribution;Prediction

    doi 10.3969/j.issn.1000-2561.2014.11.027

    由椰子致死性黃化植原體(Coconut lethal yellowing phytoplasma)引起的椰子致死性黃化?。↙ethal Yellowing,LY)是一種毀滅性病害。病害傳播快、致病性高,從開始出現(xiàn)癥狀至全株死亡,一般只需要3~6個月,各齡椰樹均可感病[1]。除侵染椰樹外,椰子致死性黃化植原體還可造成油棕、棗椰、馬尼拉棕等35種棕櫚植物的死亡[2]。在牙買加,1961年有600萬株椰樹感病,至1981年約80%的椰樹死亡;多哥在1964年約50%的椰樹(近6萬株)被毀,美國佛羅里達州的椰樹也由于該病而造成嚴重的損失[3-4]。

    椰子致死性黃化病于19世紀末在牙買加首次被發(fā)現(xiàn)[1],目前主要分布于拉丁美洲的加勒比海地區(qū)和非洲等地[5],中國還未見有關(guān)發(fā)現(xiàn)此病害的報道[4],但在海關(guān)檢疫過程中發(fā)現(xiàn)過此病原。1992年農(nóng)業(yè)部發(fā)布的《中華人民共和國進境植物檢疫危險性病、蟲、雜草名錄》中,椰子致死性黃化植原體被列為一類有害生物,2007年發(fā)布的《中華人民共和國進境植物檢疫性有害生物名錄》中,該有害生物也被列為檢疫對象。朱輝等[4]遵循聯(lián)合國糧農(nóng)組織(FAO)和國家植物保護公約(IPCC)的有害生物分析準則,對椰子致死性黃化植原體傳入中國的風險進行了研究,結(jié)果發(fā)現(xiàn)椰子致死性黃化植原體的風險性R值為2.27,屬于高度危險性有害生物。因此,開展椰子致死性黃化病在中國的適生性研究是十分有必要的。

    目前,研究人員一般采用生態(tài)位模型進行適生區(qū)預(yù)測,該模型的基本原理是根據(jù)每種生物特殊的生存環(huán)境,即生態(tài)位要求,從目標物種的已知分布區(qū)出發(fā),利用數(shù)學(xué)模型歸納或模擬其生態(tài)位需求,然后將其投射到目標地區(qū)即可得到目標物種的適生區(qū)分布[6]。Maxent是近幾年開發(fā)出來的物種地理分布預(yù)測軟件,它根據(jù)物種現(xiàn)實分布點和現(xiàn)實分布地區(qū)的環(huán)境變量運算得出預(yù)測模型,再利用此模型模擬目標物種在目標地區(qū)的可能分布情況[7]。國內(nèi)相關(guān)研究者已經(jīng)利用Maxent軟件對甜菜孢囊線蟲[8]、大豆南北方莖潰瘍病菌[9]、玉米霜霉病[10]、麥瘟病[11]、橡膠樹棒孢霉落葉病[12]等植物病害在中國的適生區(qū)進行了預(yù)測。本研究利用Maxent軟件并結(jié)合椰樹的分布情況對椰子致死性黃化植原體在中國的潛在適生區(qū)進行預(yù)測,可為進一步開展該病害在中國的分布區(qū)域測報及制定有效的檢疫措施和防治決策提供依據(jù)。

    1 材料與方法

    1.1 材料

    1.1.1 軟件來源 所用Maxent軟件在MAXENT主頁(http://www.cs.princeton.edu/~schapire/maxent/)下載,版本為3.3.3k;GIS軟件采用ESRI公司開發(fā)的ArcMap 9.3。

    1.1.2 環(huán)境數(shù)據(jù) 從WorldClim(http://www.worldclim.org/)免費下載19個生物氣候變量,數(shù)據(jù)分辨率為5 arc-minutes。

    1.1.3 地圖數(shù)據(jù) 從國家基礎(chǔ)地理信息系統(tǒng)(http://nfgis.nsdi.gov.cn/)下載1 ∶ 400萬的中國行政區(qū)劃圖,將其作為中國分析的底圖。

    1.1.4 椰子致死性黃化病分布數(shù)據(jù)的收集與處理 椰子致死性黃化病分布數(shù)據(jù)的獲取主要根據(jù)國內(nèi)外公開發(fā)表的相關(guān)論文,目前該病害主要分布在拉美地區(qū)的開曼群島、巴哈馬群島、伯利茲、危地馬拉、古巴、多米尼加、海地、洪都拉斯、墨西哥、美國佛羅里達州南部地區(qū)及非洲的加納、喀麥隆、尼日利亞、坦桑尼亞、肯尼亞、莫桑比克等地[4-5,13-16]。根據(jù)Maxent軟件要求,將椰子致死性黃化病實際分布點按物種名、分布點經(jīng)度和分布點緯度的順序儲存成后綴名為csv格式的文件,其中東經(jīng)和北緯為正,西經(jīng)和南緯為負。

    1.1.5 椰樹種植分布 中國椰樹種植分布在海南,廣東、廣西、云南、臺灣等地也有零星種植[17-18]。

    1.2 方法

    1.2.1 模型分析方法 將分布數(shù)據(jù)和環(huán)境數(shù)據(jù)導(dǎo)入Maxent,隨機選取25%的分布點作為測試集(test data),剩余的作為訓(xùn)練集(training data),其他參數(shù)均為軟件默認值,Maxent軟件分析結(jié)果的輸出格式為ASCII柵格圖層,給出的是待測物種在全球的適生性概率,取值范圍在0~1之間。將預(yù)測結(jié)果導(dǎo)入GIS軟件ArcMap 9.3后與1 ∶ 400萬的中國行政區(qū)劃圖相疊加,利用空間分析工具Spatial Analyst將中國的預(yù)測結(jié)果從全球預(yù)測結(jié)果中切割出來,得到椰子致死性黃化植原體在中國的預(yù)測分布圖,同時與中國椰樹種植分布圖疊加可得到該病菌在中國椰樹種植區(qū)的預(yù)測分布圖。

    參考已有文獻[10-11],將風險等級劃分為高風險區(qū)(>0.25)、中風險區(qū)(0.05~0.25)、低風險區(qū)(0.0001~0.05)和基本不發(fā)生區(qū)(0~0.000 1),使預(yù)測結(jié)果分級顯示。

    1.2.2 模型預(yù)測效果檢驗 對風險分析的結(jié)果進行精度檢驗時,采用ROC(receiver operating characteristic)曲線分析法[19]。ROC曲線分析法最開始是應(yīng)用于雷達信號分析,后來在醫(yī)學(xué)顯像領(lǐng)域應(yīng)用非常廣泛,常用于診斷試驗評價。ROC曲線是以真陽性率為縱坐標、假陽性率為橫坐標所形成的曲線,其中真陽性率是指實際有分布且被預(yù)測為有分布的概率,假陽性率是指實際沒有分布但被預(yù)測為有分布的概率。ROC曲線下的面積即AUC值(area under curve)可以用來反映模型預(yù)測結(jié)果的可靠性,AUC值越大表示與隨機分布相距越遠,環(huán)境變量與預(yù)測的物種地理分布模型之間的相關(guān)性越大,即模型預(yù)測效果越好[19-20]。

    Maxent軟件在預(yù)測結(jié)果中直接繪制出了ROC曲線并給出了模型的AUC值,為模型預(yù)測效果的判斷提供了方便。

    2 結(jié)果與分析

    2.1 Maxent模型預(yù)測結(jié)果

    利用GIS將Maxent運行的結(jié)果與中國底圖疊加,得到椰子致死性黃化植原體在中國的適生性分布結(jié)果(圖1)。結(jié)果表明,中國華南沿海地區(qū)及海南大部分地區(qū)屬于中風險區(qū)和高風險區(qū),長江中下游地區(qū)、華南北部、貴州和重慶全部、四川東部、云南大部、河南中部和南部、山東南部等屬于低風險區(qū),其他地區(qū)均為基本不發(fā)生區(qū)。

    2.2 椰樹種植區(qū)適生分布范圍

    由于椰子致死性黃化植原體屬于專性寄生病原物,只能在有寄主分布的地區(qū)才會引起病害,椰樹是該病原物的主要寄主,因此結(jié)合中國椰樹的種植分布范圍,可以得到該病原物在中國椰樹種植區(qū)的潛在分布(圖2)。從圖中可以看出,中國主要的椰樹種植區(qū)都是椰子致死性黃化病的潛在適生區(qū),其中海南的北部和東北部及廣東的雷州半島是高風險區(qū),海南中北部、西部和南部、東南部、廣西的南部和廣東除雷州半島外的沿海地區(qū)均為中風險區(qū),海南中南部、云南南部、臺灣大部分地區(qū)均為低風險區(qū)。

    2.3 生物氣候變量對適生分布范圍的影響

    關(guān)于各生物氣候變量因子對椰子致死性黃化植原體分布影響的大小,本研究采用Maxent軟件中的刀切法,結(jié)果見圖3。所選19個生物氣候變量中,以溫度變化反差(bio4)、年溫變化范圍(bio7)、最冷月份最低溫度(bio6)和最冷季度平均溫度(bio11)對病害分布范圍預(yù)測結(jié)果的影響較大,說明溫度特別是最低溫度對病原的分布有重要的影響。

    2.4 精度檢驗

    采用ROC曲線分析法對應(yīng)用Maxent軟件預(yù)測的椰子致死性黃化植原體適生性分布結(jié)果進行精度檢驗(圖4)。由圖中可知,Maxent模型的AUC值為0.992(非常接近1),比隨機分布模型的AUC(0.500)大,表明預(yù)測結(jié)果具有較高的精度。

    3 討論與結(jié)論

    本研究利用Maxent生態(tài)位模型和ArcGIS軟件,對椰子致死性黃化植原體在中國的潛在分布進行了研究。結(jié)果表明,中國云貴、華南、長江中下游、東南沿海和華北南部地區(qū)均是該病原的潛在適生區(qū),特別是廣西大部分地區(qū)、海南大部分地區(qū)、廣東西部、南部和東南沿海及福建東南沿海和臺灣部分沿海地區(qū)屬于高風險區(qū)和中風險區(qū)(圖1)。包括油棕、棗椰、馬尼拉棕等在內(nèi)的棕櫚科植物在中國主要種植于云南、廣西、廣東、海南和臺灣等地,而椰子致死性黃化植原體在中國的中風險區(qū)和高風險區(qū)均位于這些地區(qū),因此一旦該病原傳入中國,將給中國棕櫚科植物帶來重大損失。植原體病菌是一種專性寄生生物,其適生區(qū)域不但由環(huán)境條件決定, 寄主的分布也是決定其定殖擴散的重要因素,預(yù)測結(jié)果和國內(nèi)寄主范圍的分布情況比較吻合也驗證了預(yù)測結(jié)果的可靠性。研究結(jié)果經(jīng)ROC曲線分析法驗證可知,Maxent模型的AUC值為0.992(非常接近1),表明預(yù)測結(jié)果具有較高的精度。

    椰樹是中國重要的熱帶經(jīng)濟作物,對椰子致死性黃化植原體在中國椰樹種植區(qū)的潛在分布的研究結(jié)果表明,中國主要的椰樹種植區(qū)都是椰子致死性黃化病的高風險區(qū)和中風險區(qū)(圖2),因此對于這種尚缺乏有效的防治技術(shù)的毀滅性病害,加強對該病原的檢疫及監(jiān)測對確保中國椰子產(chǎn)業(yè)的健康發(fā)展十分重要。

    與已報道的椰子致死性黃化植原體傳入中國的風險分析研究[4]相比,本研究利用Maxent軟件對椰子致死性黃化植原體在中國的潛在適生區(qū)進行研究,可以了解椰子致死性黃化植原體對中國不同地區(qū)的潛在威脅,為決策者制定相應(yīng)的管理措施提供技術(shù)支撐。

    目前利用模型預(yù)測的植物病原生物包括真菌[9-12]、細菌[21]、 線蟲[8, 22]和病毒[23]等,關(guān)于植原體適生區(qū)預(yù)測的研究報道很少,這可能是由于植原體的生物學(xué)參數(shù)難以準確獲得,而Climex、Dymex等模型需要根據(jù)病原的生物學(xué)特性調(diào)整參數(shù)來預(yù)測其潛在的地理分布。但包括Maxent、Bioclim、GARP等在內(nèi)的模型則不需要任何病原的生物學(xué)資料,只需知道物種已有的分布,然后根據(jù)環(huán)境變量預(yù)測該物種在其他地區(qū)的分布情況。此外,研究發(fā)現(xiàn)Maxent模型的結(jié)果要優(yōu)于GARP、Bioclim等模型,特別是在物種分布數(shù)據(jù)不全的情況下,Maxent模型仍然能得到較為滿意的結(jié)果[19]。對相似穿孔線蟲和飛機草在中國的適生區(qū)預(yù)測研究發(fā)現(xiàn),Maxent模型預(yù)測結(jié)果較GARP好[22,24]。因此本研究選用Maxent模型對椰子致死性黃化植原體在中國的潛在分布情況進行研究,這也為其他植原體的適生性研究提供借鑒。

    本研究所依據(jù)的環(huán)境變量僅包括WORLDCLIM中的19個生物氣候變量,沒有考慮地形、土壤類型等對病害分布的影響;此外,全球氣候變化應(yīng)該也會對椰子致死性黃化植原體的分布產(chǎn)生影響。雖然本研究是對椰子致死性黃化植原體在中國及中國椰樹種植區(qū)的潛在分布的初步研究結(jié)果,但從一定程度上表明該病害對中國椰子產(chǎn)業(yè)的潛在威脅巨大。

    參考文獻

    [1] Eden-Green S J. History, distribution and present status of lethal yellowing-like disease of palms[C]//In Proceedings of an International Workshop on Lethal Yellowing-like Diseases of Coconut. Elmina, Ghana, 1995, 11: 9-25.

    [2] Harrison N A, Cordova I, Richardson P, et a1. Detection and diagnosis of lethal yellowing[M]//Oropoza C, Verdeil J L, Ashburlier G R, et a1. Current Advances in Coconut Biotechnology. Dordrecht: Kluwer Academic Publishers, 1999: 183-196.

    [3] 余鳳玉, 覃偉權(quán), 朱 輝, 等. 椰子致死性黃化病研究進展[J]. 中國熱帶農(nóng)業(yè), 2008(1): 42-44.

    [4]朱 輝, 余鳳玉, 吳多揚, 等. 椰子致死性黃化植原體傳入中國的風險性分析[J]. 江西農(nóng)業(yè)學(xué)報, 2010, 22(11): 84-87.

    [5] Bonnot F, de Franqueville H, Lourenco, E. Spatial and spatiotemporal pattern analysis of coconut lethal yellowing in Mozambique[J]. Phytopathology, 2010, 100(4): 300-312.

    [6] 孫文濤, 劉雅婷. 生物入侵風險分析的研究進展[J]. 中國農(nóng)學(xué)通報, 2010, 26(7): 233-236.

    [7] Phillips S J, Anderson R P, Schapire R E. Maximum entropy modeling of species geographic distributions[J]. Ecological Modelling, 2006, 190: 231-259.

    [8]李建中, 彭德良, 劉淑艷. 潛在外來入侵甜菜孢囊線蟲在中國的適生性風險分析[J]. 植物保護, 2008, 34(5): 90-94.

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

    [10]趙文娟, 陳 林, 丁克堅, 等. 利用MAXENT預(yù)測玉米霜霉病在中國的適生區(qū)[J]. 植物保護, 2009, 35(2): 32-38.

    [11]曹學(xué)仁, 陳 林, 周益林, 等. 基于MaxEnt的麥瘟病在全球及中國的潛在分布區(qū)預(yù)測[J]. 植物保護, 2011, 37(3): 80-83.

    [12]曲偉偉, 李志紅, 黃貴修, 等. 利用MAXENT預(yù)測橡膠樹棒孢霉落葉病在中國的適生區(qū)[J]. 植物保護, 2011, 37(4): 52-57.

    [13] Ashbumer G R, CórdovaⅡ, Oropeza C M, et a1. First report of coconut lethal yellowing disease in Honduras[J]. Plant Disease, 1996, 80(8): 960.

    [14] Harrison N A, Bourne C M, Cox R L, et al. DNA probes for detection of mycoplasmalike organisms associated with lethal yellowing disease of palms in Florida[J]. Phytopathology, 1992, 82(2): 216-224.

    [15] Cardena R, Villanueva M A, Santamaria J M, et al. Presence in Yucatan of mycoplasmalike organisms in Cocos nucifera palms showing lethal yellowing disease symptoms[J]. Canadian Journal of Plant Pathology, 1991, 13(2): 135-138.

    [16] Mpunami A, Tymon A, Jones P, et al. Genetic diversity in the coconut lethal yellowing disease phytoplasmas of East Africa[J]. Plant Pathology, 1999, 48(1): 109-114.

    [17] 馮美利, 李 杰, 王必尊. 我國椰子綜合研究進展概述[J]. 中國熱帶農(nóng)業(yè), 2007(5): 30-31.

    [18] 陳豪軍, 陳永森, 林聰田, 等. 基于OpenModeller-GARP的椰子在中國的適生區(qū)預(yù)測[J]. 南方農(nóng)業(yè)學(xué)報, 2012, 43(10): 1 525-1 529.

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

    [20] Hanley J A, Mcneil B J. The meaning and use of the area under a Receiver Operating Characteristic(ROC)curve[J]. Radiology, 1982, 143(1): 29-36.

    [21]徐 進, 陳 林, 許景生, 等. 香蕉細菌性枯萎病菌在中國的潛在適生區(qū)域[J]. 植物保護學(xué)報, 2008, 35(3): 233-238.

    [22]王運生, 謝丙炎, 萬方浩, 等. 相似穿孔線蟲在中國的適生區(qū)預(yù)測[J]. 中國農(nóng)業(yè)科學(xué), 2007, 40(11): 2 502-2 506.

    [23]鄭 耘, 楊偉東, 章桂明, 等. 李痘病毒適生區(qū)預(yù)測及模型選擇[J]. 植物檢疫, 2009, 23(3): 1-4.

    [24] 楊 波, 薛躍規(guī), 唐小飛, 等. 外來入侵植物飛機草在中國的適生區(qū)預(yù)測[J]. 植物保護, 2009, 35(4): 70-73.

    猜你喜歡
    分布預(yù)測
    無可預(yù)測
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測卷(A卷)
    選修2-2期中考試預(yù)測卷(B卷)
    選修2—2期中考試預(yù)測卷(A卷)
    不可預(yù)測
    不必預(yù)測未來,只需把握現(xiàn)在
    大葉千斤拔活性成分分布及積累動態(tài)
    28例醫(yī)療糾紛起訴案件特點分析
    “一帶一路”沿線直接投資分布與挑戰(zhàn)應(yīng)對
    剩余油分布狀況研究綜述
    大型黄色视频在线免费观看| 天天一区二区日本电影三级 | 亚洲精品一区av在线观看| 老鸭窝网址在线观看| 波多野结衣av一区二区av| 少妇熟女aⅴ在线视频| 国产高清videossex| 欧美日韩乱码在线| 国产精品国产高清国产av| 国产精品 欧美亚洲| 久久中文字幕一级| 亚洲第一av免费看| 巨乳人妻的诱惑在线观看| 精品一区二区三区视频在线观看免费| av福利片在线| 又大又爽又粗| 免费人成视频x8x8入口观看| 亚洲激情在线av| 1024香蕉在线观看| 久久久久久国产a免费观看| or卡值多少钱| 亚洲欧美一区二区三区黑人| 久久久久久大精品| а√天堂www在线а√下载| av网站免费在线观看视频| 一本久久中文字幕| 别揉我奶头~嗯~啊~动态视频| 国内精品久久久久久久电影| 日本撒尿小便嘘嘘汇集6| 人妻久久中文字幕网| 亚洲国产精品久久男人天堂| 久久久久久久久免费视频了| 热re99久久国产66热| 国产成+人综合+亚洲专区| 欧美不卡视频在线免费观看 | 日韩欧美在线二视频| 国产精品一区二区免费欧美| 成人亚洲精品av一区二区| 国产极品粉嫩免费观看在线| 国产精品亚洲av一区麻豆| 熟女少妇亚洲综合色aaa.| 国产aⅴ精品一区二区三区波| 两个人看的免费小视频| 久久国产亚洲av麻豆专区| 精品久久久精品久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 97人妻天天添夜夜摸| 在线观看www视频免费| 日本免费一区二区三区高清不卡 | 少妇的丰满在线观看| 法律面前人人平等表现在哪些方面| 淫妇啪啪啪对白视频| 香蕉久久夜色| 变态另类成人亚洲欧美熟女 | 国产一区二区激情短视频| 亚洲av第一区精品v没综合| 十八禁人妻一区二区| 国产极品粉嫩免费观看在线| 亚洲avbb在线观看| 国内精品久久久久精免费| 91九色精品人成在线观看| 国产精品99久久99久久久不卡| 亚洲一区二区三区色噜噜| 免费看a级黄色片| 美女免费视频网站| 一二三四社区在线视频社区8| 黄色视频,在线免费观看| 人成视频在线观看免费观看| 欧美日韩一级在线毛片| 91麻豆av在线| 国产成人av激情在线播放| 亚洲自偷自拍图片 自拍| 久久香蕉精品热| 精品人妻1区二区| 午夜福利一区二区在线看| 久久久久久国产a免费观看| 免费搜索国产男女视频| 中文字幕人妻熟女乱码| 大码成人一级视频| 黑人欧美特级aaaaaa片| 无限看片的www在线观看| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久精品电影 | 亚洲色图 男人天堂 中文字幕| 亚洲国产中文字幕在线视频| 国产免费男女视频| 亚洲专区中文字幕在线| 国产精品一区二区三区四区久久 | 99国产综合亚洲精品| tocl精华| 成人18禁在线播放| av电影中文网址| 日本免费一区二区三区高清不卡 | 精品国内亚洲2022精品成人| 国产精品1区2区在线观看.| 欧美日本亚洲视频在线播放| 欧美黑人欧美精品刺激| 最近最新中文字幕大全免费视频| 三级毛片av免费| 亚洲欧美日韩高清在线视频| 欧美日韩精品网址| 国产精品久久久av美女十八| 久久久久久久午夜电影| 给我免费播放毛片高清在线观看| 成人18禁在线播放| 欧美不卡视频在线免费观看 | 久久精品亚洲精品国产色婷小说| 男人舔女人的私密视频| 国产91精品成人一区二区三区| 两个人视频免费观看高清| 免费在线观看亚洲国产| 久久青草综合色| 精品少妇一区二区三区视频日本电影| 久久人人精品亚洲av| 99在线人妻在线中文字幕| 给我免费播放毛片高清在线观看| 黄色丝袜av网址大全| 成年版毛片免费区| 丁香欧美五月| 嫩草影院精品99| 欧美日韩福利视频一区二区| 精品一区二区三区av网在线观看| 亚洲人成网站在线播放欧美日韩| av欧美777| 性色av乱码一区二区三区2| 国产成人av教育| 九色国产91popny在线| 国产精品国产高清国产av| 久久人人精品亚洲av| 99在线人妻在线中文字幕| 高清毛片免费观看视频网站| 禁无遮挡网站| 国产精华一区二区三区| 久久精品国产亚洲av高清一级| 搡老妇女老女人老熟妇| 久久久精品国产亚洲av高清涩受| 一级毛片精品| 在线观看免费视频日本深夜| 国产真人三级小视频在线观看| 亚洲精品美女久久av网站| av在线天堂中文字幕| 久久久久久久午夜电影| 啪啪无遮挡十八禁网站| bbb黄色大片| 国内精品久久久久久久电影| 亚洲激情在线av| 97人妻精品一区二区三区麻豆 | 亚洲国产日韩欧美精品在线观看 | 视频在线观看一区二区三区| 国产成人精品在线电影| 亚洲性夜色夜夜综合| 国产亚洲精品久久久久5区| 久久人妻熟女aⅴ| 后天国语完整版免费观看| 少妇的丰满在线观看| 啦啦啦观看免费观看视频高清 | 极品教师在线免费播放| 日韩大码丰满熟妇| 很黄的视频免费| 丰满人妻熟妇乱又伦精品不卡| 国产免费男女视频| 麻豆久久精品国产亚洲av| 久久久久亚洲av毛片大全| 免费在线观看完整版高清| а√天堂www在线а√下载| 非洲黑人性xxxx精品又粗又长| 最近最新中文字幕大全电影3 | av免费在线观看网站| 一二三四社区在线视频社区8| 欧美色欧美亚洲另类二区 | 如日韩欧美国产精品一区二区三区| 亚洲中文日韩欧美视频| 国产一区在线观看成人免费| a在线观看视频网站| 成人欧美大片| 国产不卡一卡二| 免费女性裸体啪啪无遮挡网站| 久久亚洲真实| aaaaa片日本免费| 亚洲欧美激情综合另类| 久久人妻av系列| 国产精品1区2区在线观看.| 亚洲片人在线观看| 亚洲国产精品成人综合色| 黄色丝袜av网址大全| 亚洲三区欧美一区| 国产高清视频在线播放一区| 宅男免费午夜| 午夜亚洲福利在线播放| 午夜精品在线福利| 亚洲精品国产色婷婷电影| 欧美乱色亚洲激情| 一级a爱片免费观看的视频| 免费高清在线观看日韩| 天天添夜夜摸| 亚洲精华国产精华精| 亚洲第一欧美日韩一区二区三区| 91精品国产国语对白视频| 亚洲欧美激情在线| 一区二区日韩欧美中文字幕| 999久久久精品免费观看国产| 搡老熟女国产l中国老女人| 国产视频一区二区在线看| 午夜福利欧美成人| 日韩三级视频一区二区三区| 好男人在线观看高清免费视频 | 国产精品自产拍在线观看55亚洲| 久99久视频精品免费| 国产成人啪精品午夜网站| 好看av亚洲va欧美ⅴa在| ponron亚洲| 在线观看免费视频网站a站| 亚洲精品中文字幕一二三四区| 精品熟女少妇八av免费久了| 9色porny在线观看| 国产一级毛片七仙女欲春2 | 黄色视频不卡| 国产欧美日韩一区二区精品| 国内精品久久久久精免费| 亚洲专区中文字幕在线| 天堂动漫精品| 国产精品 欧美亚洲| 亚洲成人久久性| 国产免费av片在线观看野外av| av视频在线观看入口| 亚洲色图 男人天堂 中文字幕| 欧美色视频一区免费| 欧美成人免费av一区二区三区| 亚洲欧美精品综合久久99| 国产成人免费无遮挡视频| 欧美日韩瑟瑟在线播放| 午夜两性在线视频| 精品国产一区二区三区四区第35| 欧美激情 高清一区二区三区| 成在线人永久免费视频| 午夜精品在线福利| 制服人妻中文乱码| 亚洲色图av天堂| 久久久国产欧美日韩av| 精品久久久久久久毛片微露脸| 国产av又大| 国产成人av激情在线播放| 老司机在亚洲福利影院| 国产精品香港三级国产av潘金莲| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 精品国产乱子伦一区二区三区| 国产精品亚洲一级av第二区| 国产一区二区三区在线臀色熟女| 啦啦啦韩国在线观看视频| 男女下面进入的视频免费午夜 | 亚洲免费av在线视频| 亚洲精品在线美女| 亚洲精华国产精华精| 91九色精品人成在线观看| 国产精品一区二区三区四区久久 | 欧美成人免费av一区二区三区| 精品电影一区二区在线| 亚洲成av人片免费观看| 曰老女人黄片| 国产高清有码在线观看视频 | 性少妇av在线| 成人国产一区最新在线观看| 在线观看日韩欧美| 色老头精品视频在线观看| 视频区欧美日本亚洲| 国产精品一区二区免费欧美| 国产一卡二卡三卡精品| 久久伊人香网站| 无遮挡黄片免费观看| 一区在线观看完整版| 国产精品1区2区在线观看.| 制服丝袜大香蕉在线| 日日夜夜操网爽| 国产免费av片在线观看野外av| 国产区一区二久久| 婷婷精品国产亚洲av在线| 亚洲色图 男人天堂 中文字幕| 国产精品一区二区精品视频观看| 狠狠狠狠99中文字幕| 看片在线看免费视频| 亚洲av美国av| 日韩欧美国产在线观看| 女人高潮潮喷娇喘18禁视频| 无人区码免费观看不卡| 精品电影一区二区在线| 精品熟女少妇八av免费久了| 999精品在线视频| 久久亚洲真实| 成年版毛片免费区| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 亚洲av第一区精品v没综合| 亚洲av成人一区二区三| 露出奶头的视频| 欧美最黄视频在线播放免费| 亚洲精品美女久久久久99蜜臀| 午夜久久久在线观看| 欧美av亚洲av综合av国产av| 夜夜夜夜夜久久久久| 亚洲专区中文字幕在线| 国产极品粉嫩免费观看在线| 18禁观看日本| 在线观看日韩欧美| 老司机靠b影院| 在线观看午夜福利视频| 国产精品免费视频内射| 最近最新免费中文字幕在线| 亚洲色图av天堂| 一夜夜www| 国产日韩一区二区三区精品不卡| 久久久久亚洲av毛片大全| tocl精华| 亚洲国产高清在线一区二区三 | 日韩大码丰满熟妇| 9191精品国产免费久久| 亚洲精品美女久久久久99蜜臀| 美女扒开内裤让男人捅视频| 午夜福利高清视频| 国产精品久久视频播放| 亚洲专区国产一区二区| 亚洲精品久久国产高清桃花| 在线观看66精品国产| 亚洲七黄色美女视频| 麻豆一二三区av精品| 深夜精品福利| 久久久久久久久免费视频了| 国产精品二区激情视频| 欧美久久黑人一区二区| 美女高潮喷水抽搐中文字幕| 精品第一国产精品| 欧美激情 高清一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 国产又色又爽无遮挡免费看| 人妻久久中文字幕网| 亚洲国产毛片av蜜桃av| 91麻豆av在线| 一区福利在线观看| 成人18禁在线播放| 午夜免费观看网址| 亚洲男人的天堂狠狠| 国产精品亚洲一级av第二区| 欧美老熟妇乱子伦牲交| 淫秽高清视频在线观看| 好男人在线观看高清免费视频 | 搞女人的毛片| 日韩三级视频一区二区三区| 久久久久国产一级毛片高清牌| 国产视频一区二区在线看| 国产一区二区三区在线臀色熟女| 欧洲精品卡2卡3卡4卡5卡区| 999久久久国产精品视频| 日韩av在线大香蕉| 欧美一级a爱片免费观看看 | 国产成人免费无遮挡视频| 99re在线观看精品视频| 久久久国产成人免费| 色婷婷久久久亚洲欧美| 国产欧美日韩精品亚洲av| 色av中文字幕| 波多野结衣高清无吗| 一级,二级,三级黄色视频| 色综合婷婷激情| 国产精品一区二区三区四区久久 | 亚洲午夜理论影院| 免费看美女性在线毛片视频| 香蕉久久夜色| 日日干狠狠操夜夜爽| 中文字幕精品免费在线观看视频| 精品免费久久久久久久清纯| 亚洲av五月六月丁香网| 亚洲五月天丁香| 亚洲精品一卡2卡三卡4卡5卡| 一级片免费观看大全| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 欧美精品亚洲一区二区| 一边摸一边抽搐一进一小说| 精品国产亚洲在线| 男人操女人黄网站| 97碰自拍视频| 国产亚洲欧美98| 久久久国产欧美日韩av| 一卡2卡三卡四卡精品乱码亚洲| 国产在线观看jvid| 日韩欧美在线二视频| av视频免费观看在线观看| 久久伊人香网站| 日韩欧美一区二区三区在线观看| 国产精品一区二区免费欧美| 亚洲精品粉嫩美女一区| 看黄色毛片网站| 大陆偷拍与自拍| 一本大道久久a久久精品| 精品国产乱码久久久久久男人| 熟妇人妻久久中文字幕3abv| 一个人观看的视频www高清免费观看 | 亚洲精品一区av在线观看| 国产精品秋霞免费鲁丝片| 日本免费一区二区三区高清不卡 | 精品一品国产午夜福利视频| 午夜久久久久精精品| 亚洲激情在线av| 久久久久久久久免费视频了| 琪琪午夜伦伦电影理论片6080| 亚洲国产欧美网| 国产精品精品国产色婷婷| 成人永久免费在线观看视频| 国产精品久久久久久人妻精品电影| 色综合婷婷激情| 亚洲九九香蕉| 精品福利观看| 成人手机av| 每晚都被弄得嗷嗷叫到高潮| 黄色片一级片一级黄色片| 午夜免费成人在线视频| 成熟少妇高潮喷水视频| 在线永久观看黄色视频| 日韩大尺度精品在线看网址 | 美女国产高潮福利片在线看| 99riav亚洲国产免费| 色哟哟哟哟哟哟| 在线观看日韩欧美| 色在线成人网| 日韩欧美在线二视频| 亚洲成人国产一区在线观看| 一级毛片精品| 亚洲成人精品中文字幕电影| 午夜成年电影在线免费观看| 亚洲,欧美精品.| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品在线电影| 757午夜福利合集在线观看| 99riav亚洲国产免费| 午夜福利免费观看在线| 亚洲av电影不卡..在线观看| 最好的美女福利视频网| 久99久视频精品免费| 国产午夜福利久久久久久| 这个男人来自地球电影免费观看| 国产精品乱码一区二三区的特点 | 高清在线国产一区| 男人舔女人下体高潮全视频| 国产区一区二久久| 嫩草影视91久久| 91麻豆精品激情在线观看国产| 精品欧美一区二区三区在线| 老汉色∧v一级毛片| 国产精品久久久久久人妻精品电影| 一区二区三区国产精品乱码| 好男人电影高清在线观看| 最近最新中文字幕大全免费视频| 99国产精品一区二区三区| 国产高清有码在线观看视频 | 欧美在线黄色| 国产麻豆69| 国产免费男女视频| 亚洲片人在线观看| 亚洲人成伊人成综合网2020| 亚洲精品在线美女| 丝袜人妻中文字幕| 成年版毛片免费区| 日本a在线网址| 午夜福利影视在线免费观看| 国产一区二区激情短视频| 高潮久久久久久久久久久不卡| 在线永久观看黄色视频| 亚洲男人的天堂狠狠| 一级a爱片免费观看的视频| 亚洲免费av在线视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜两性在线视频| 这个男人来自地球电影免费观看| 国产精品98久久久久久宅男小说| 免费人成视频x8x8入口观看| 少妇熟女aⅴ在线视频| 日本免费a在线| 好看av亚洲va欧美ⅴa在| 黑丝袜美女国产一区| 国产成人免费无遮挡视频| 女生性感内裤真人,穿戴方法视频| 99精品在免费线老司机午夜| 日韩三级视频一区二区三区| 怎么达到女性高潮| 国产一区二区三区综合在线观看| 精品国产乱子伦一区二区三区| 国产欧美日韩精品亚洲av| 丁香六月欧美| 淫秽高清视频在线观看| 午夜日韩欧美国产| 欧美日韩中文字幕国产精品一区二区三区 | 久久久久精品国产欧美久久久| 欧美黑人精品巨大| 欧美成人免费av一区二区三区| 国产精品久久久av美女十八| 欧美黄色淫秽网站| 亚洲美女黄片视频| 女人被狂操c到高潮| 在线视频色国产色| 久久精品国产99精品国产亚洲性色 | 亚洲第一av免费看| 老司机午夜福利在线观看视频| 欧美av亚洲av综合av国产av| 久久久国产成人精品二区| 天天躁狠狠躁夜夜躁狠狠躁| 很黄的视频免费| 久久人人精品亚洲av| 国产成人欧美在线观看| 村上凉子中文字幕在线| 亚洲人成伊人成综合网2020| 男女午夜视频在线观看| 久久草成人影院| 亚洲男人天堂网一区| 国产亚洲精品一区二区www| www.www免费av| 亚洲伊人色综图| 丰满人妻熟妇乱又伦精品不卡| 琪琪午夜伦伦电影理论片6080| 亚洲午夜精品一区,二区,三区| 国产熟女午夜一区二区三区| 日本欧美视频一区| 99精品在免费线老司机午夜| 欧美在线一区亚洲| 悠悠久久av| 可以免费在线观看a视频的电影网站| 免费在线观看完整版高清| 香蕉久久夜色| 午夜福利视频1000在线观看 | 好男人在线观看高清免费视频 | 精品久久久久久,| 老司机福利观看| 一a级毛片在线观看| 精品久久久久久成人av| 久久亚洲精品不卡| 黄色丝袜av网址大全| 在线永久观看黄色视频| 淫妇啪啪啪对白视频| 久久久久久大精品| 男男h啪啪无遮挡| 久久久精品欧美日韩精品| 精品国产国语对白av| 一夜夜www| 精品福利观看| 午夜福利视频1000在线观看 | 国产成人av教育| 国产男靠女视频免费网站| 日韩一卡2卡3卡4卡2021年| 亚洲自拍偷在线| www.精华液| 午夜福利,免费看| 91老司机精品| 国产精品一区二区三区四区久久 | 欧美丝袜亚洲另类 | 久久国产精品男人的天堂亚洲| 亚洲avbb在线观看| 免费av毛片视频| 精品国产乱子伦一区二区三区| 欧美激情极品国产一区二区三区| 久久久久久久久中文| 成在线人永久免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 波多野结衣巨乳人妻| 欧美日韩亚洲综合一区二区三区_| 色婷婷久久久亚洲欧美| 又大又爽又粗| 99久久久亚洲精品蜜臀av| 国产成人av教育| 免费少妇av软件| 国产人伦9x9x在线观看| 亚洲av日韩精品久久久久久密| 亚洲专区国产一区二区| 欧美丝袜亚洲另类 | 男人的好看免费观看在线视频 | 午夜精品国产一区二区电影| www.精华液| 嫩草影院精品99| 亚洲欧美激情在线| 老司机午夜福利在线观看视频| 黄色视频不卡| 一级a爱视频在线免费观看| 亚洲人成伊人成综合网2020| 日韩欧美三级三区| www日本在线高清视频| 操出白浆在线播放| 国产一卡二卡三卡精品| 不卡av一区二区三区| 一级作爱视频免费观看| 国产免费av片在线观看野外av| 亚洲一区二区三区不卡视频| 亚洲一卡2卡3卡4卡5卡精品中文| 久久婷婷人人爽人人干人人爱 | 99国产精品一区二区蜜桃av| www.自偷自拍.com| 亚洲欧洲精品一区二区精品久久久| 久久久精品欧美日韩精品| 精品一区二区三区视频在线观看免费| 色av中文字幕| 久久久精品欧美日韩精品| 国产一级毛片七仙女欲春2 | 在线观看免费视频日本深夜| 中文字幕色久视频| 女性生殖器流出的白浆| 中文亚洲av片在线观看爽| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品国产精品久久久不卡| 又紧又爽又黄一区二区| 成年女人毛片免费观看观看9| 亚洲中文av在线| 亚洲 欧美 日韩 在线 免费| 亚洲精品国产一区二区精华液| 大陆偷拍与自拍| 少妇粗大呻吟视频| 久久午夜综合久久蜜桃| 亚洲欧美一区二区三区黑人|