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

    美國(guó)雜草風(fēng)險(xiǎn)評(píng)估方法對(duì)我國(guó)入侵植物的可應(yīng)用性*

    2020-06-02 09:07:22何善勇溫俊寶印麗萍
    林業(yè)科學(xué) 2020年4期
    關(guān)鍵詞:生活型植物評(píng)估

    何善勇 徐 飛 張 寧 溫俊寶 印麗萍

    (1.北京林業(yè)大學(xué)林木有害生物防治北京市重點(diǎn)實(shí)驗(yàn)室 北京 100083;2.上海海關(guān) 上海 200135)

    生物入侵在全世界范圍內(nèi)造成了嚴(yán)重的社會(huì)、經(jīng)濟(jì)和環(huán)境影響,已成為生態(tài)學(xué)研究的焦點(diǎn)之一(Pheloungetal., 1999)。相比其他類別的入侵生物,入侵植物數(shù)量更多,全球入侵物種數(shù)據(jù)庫(kù)(global invasive species database,GISD)收錄的891種外來(lái)入侵生物中,458種(51.4%)為外來(lái)入侵植物(Invasive Species Specialist Group, 2017);全球100種最危險(xiǎn)入侵生物中,入侵植物數(shù)量最多,占比高達(dá)37%(Luqueetal., 2014)。當(dāng)外來(lái)入侵植物在引入地區(qū)生存定殖后,由于缺乏天敵、競(jìng)爭(zhēng)者等外界控制因素,加之自身在繁殖、適應(yīng)能力、資源利用等多方面的優(yōu)勢(shì)(孫士國(guó)等, 2018),極易擴(kuò)散蔓延,形成大面積單優(yōu)群落,占據(jù)主導(dǎo)地位,危害本土動(dòng)植物,給經(jīng)濟(jì)、生態(tài)系統(tǒng)帶來(lái)災(zāi)難性后果(張震等, 2018),甚至對(duì)引入地區(qū)的社會(huì)、文化和人類健康造成嚴(yán)重危害(楊期和等, 2002)。

    植物貿(mào)易引入是植物入侵的主要途徑之一(趙宇翔等, 2015)。自2001年加入世界貿(mào)易組織以來(lái),我國(guó)植物的主動(dòng)引入十分頻繁,以林木引種為例,2005—2014年10年間,我國(guó)總計(jì)從國(guó)外主動(dòng)引入植物1 500余種(趙宇翔, 2017);與之相對(duì)應(yīng),我國(guó)外來(lái)入侵植物的種類開始急劇增加,數(shù)量由90種(2002年)劇增至291種(2017年)(李振宇等, 2002;賀紅, 2017),每年因外來(lái)植物入侵造成的經(jīng)濟(jì)損失高達(dá)1 000億元人民幣(賀紅, 2017)。雖然我國(guó)從20世紀(jì)80年代就已實(shí)行對(duì)外開放,但植物貿(mào)易真正飛速發(fā)展始于21世紀(jì)初葉(Normile, 2004),因此今后很長(zhǎng)時(shí)期我國(guó)外來(lái)植物的入侵形勢(shì)將十分嚴(yán)峻。

    風(fēng)險(xiǎn)評(píng)估是預(yù)防外來(lái)入侵植物最經(jīng)濟(jì)有效的手段(Lodgeetal., 2006)。我國(guó)已研究并建立了多種入侵植物多指標(biāo)綜合風(fēng)險(xiǎn)評(píng)估體系,如廈門、溫州、江浙滬地區(qū)、華東地區(qū)外來(lái)植物入侵風(fēng)險(xiǎn)評(píng)估體系(歐健等, 2006;馮幼義等, 2010;李明麗, 2010;顧慧, 2014),以及外來(lái)入侵植物雜草化風(fēng)險(xiǎn)“五階評(píng)估法”(強(qiáng)勝, 2009)、外來(lái)植物入侵風(fēng)險(xiǎn)評(píng)價(jià)指標(biāo)體系(謝國(guó)文等, 2010)、外來(lái)樹種生物入侵風(fēng)險(xiǎn)評(píng)價(jià)體系(鄭勇奇等, 2014)等,但從植物內(nèi)在屬性(生物生態(tài)學(xué)特性)建立風(fēng)險(xiǎn)評(píng)估方法用以評(píng)判植物自身是否為潛在入侵植物的研究和應(yīng)用則相對(duì)較少。鄭勇奇等(2014)建立的外來(lái)樹種生物入侵風(fēng)險(xiǎn)評(píng)價(jià)體系中,包含對(duì)樹種自身入侵性的評(píng)估方法,但其指標(biāo)主要與喬木類植物有關(guān)。Chen等(2015)基于6種植物生物生態(tài)學(xué)特性構(gòu)建多種機(jī)器學(xué)習(xí)算法入侵植物評(píng)估模型,總體評(píng)估準(zhǔn)確率為66%~82%,但其使用時(shí)需要采用相關(guān)的軟件或相關(guān)計(jì)算機(jī)語(yǔ)言,不便于廣泛應(yīng)用。國(guó)際植物保護(hù)公約(international plant protection convention,IPPC)要求對(duì)植物自身是否為有害植物(包括“雜草”和“入侵植物”)展開風(fēng)險(xiǎn)評(píng)估(ISPM 2, 2007;ISPM 11, 2013),2017年國(guó)家林業(yè)局發(fā)布的《境外林木引種檢疫審批風(fēng)險(xiǎn)評(píng)估管理規(guī)范》(中國(guó)林業(yè)網(wǎng), 2018)明確要求對(duì)外來(lái)植物自身的入侵性進(jìn)行評(píng)估,但沒(méi)有給出具體的參考方法,因此,在我國(guó)林木引種檢疫管理工作中,迫切需要評(píng)判效果較好、操作方便的入侵植物風(fēng)險(xiǎn)評(píng)估方法,用于評(píng)判擬引入植物是否為入侵植物,以降低外來(lái)入侵植物隨植物貿(mào)易引入的風(fēng)險(xiǎn)。美國(guó)雜草風(fēng)險(xiǎn)評(píng)估方法(United States weed risk assessment system, USWRA)用于評(píng)估植物自身是否為潛在入侵植物(對(duì)此,Koop等(2012)指出該方法標(biāo)題中的“Weed”和評(píng)估結(jié)果中的“Invader”等同,含義均指能在自然或人工環(huán)境中繁殖擴(kuò)散并造成經(jīng)濟(jì)或環(huán)境損害的非本土植物。本研究按中文翻譯習(xí)慣,將“Weed”翻譯為“雜草”,但其類別包含所有植物種類,不僅限于中文字面所指的草本植物),是在全世界使用最為廣泛的澳大利亞雜草風(fēng)險(xiǎn)評(píng)估方法(Australian weed risk assessment system, AWRA)基礎(chǔ)上建立而成的,并且提高了對(duì)非入侵植物的評(píng)估準(zhǔn)確率(Koopetal., 2012)。為此,本研究擬采用我國(guó)已知身份的外來(lái)植物對(duì)USWRA進(jìn)行測(cè)試,分析該方法對(duì)我國(guó)外來(lái)植物是否具有較高的評(píng)估效果,評(píng)判該方法是否適用于我國(guó)林木引種檢疫管理工作,為建立我國(guó)入侵植物風(fēng)險(xiǎn)評(píng)估方法奠定研究基礎(chǔ)。

    1 材料與方法

    USWRA將受試植物評(píng)估為嚴(yán)重入侵植物、一般入侵植物或非入侵植物(Koopetal., 2012)。本研究采用USWRA對(duì)我國(guó)已知身份的外來(lái)入侵植物(包括嚴(yán)重入侵和一般入侵)和外來(lái)非入侵植物進(jìn)行逐一評(píng)估,整理并分析USWRA對(duì)各類植物的評(píng)估效果(拒絕率、接受率和再評(píng)估率),采用ROC曲線對(duì)USWRA的評(píng)估效果進(jìn)行檢驗(yàn)。

    1.1 受試植物的選擇 外來(lái)入侵植物:從馬金雙(2013)主編的《中國(guó)入侵植物名錄》和黃冠勝(2014)主編的《中國(guó)外來(lái)生物入侵與檢疫防范》2部著作中,選出55種外來(lái)入侵植物(前書48種,后書7種)進(jìn)行測(cè)試,其中嚴(yán)重入侵植物19種、一般入侵植物36種。本次測(cè)試中,將《中國(guó)入侵植物名錄》一書中的全國(guó)性入侵植物劃分為嚴(yán)重入侵植物,局部入侵植物劃分為一般入侵植物,來(lái)自《中國(guó)外來(lái)生物入侵與檢疫防范》中的7種入侵植物按《中國(guó)入侵植物名錄》的定義歸為一般入侵植物。

    受試入侵植物選擇時(shí)主要從傳入時(shí)間和分類地位2方面設(shè)置篩選條件。傳入時(shí)間方面,外來(lái)植物要成為入侵植物往往需要幾十年時(shí)間(Williamson, 1996),但具體年份未知,本研究根據(jù)Nishida等(2009)的測(cè)試,確定傳入時(shí)間至少為40年。由于我國(guó)近代之前(即1840年以前)植物貿(mào)易很少(趙宇翔, 2017),因此受試外來(lái)入侵植物的傳入時(shí)間主要為近代以后(1840年以后)且至少已傳入40年。分類地位方面,根據(jù)傳入時(shí)間要求,從2部著作中整理出48科189種外來(lái)入侵植物,從每科中隨機(jī)抽取1種用于測(cè)試;為保證受試植物生活型更加豐富,增加7種木本和藤本植物,分屬豆科(Fabacea)、菊科(Asteraceae)和西番蓮科(Passiflorweae),最終選出55種受試外來(lái)入侵植物。

    外來(lái)非入侵植物:從何家慶(2012)編著的《中國(guó)外來(lái)植物》中選出76種外來(lái)非入侵植物進(jìn)行測(cè)試。非入侵植物主要按照未在我國(guó)造成經(jīng)濟(jì)或環(huán)境危害這一標(biāo)準(zhǔn)進(jìn)行篩選(Daehleretal., 2004)。篩選時(shí),按照上述外來(lái)入侵植物的傳入時(shí)間要求,從《中國(guó)外來(lái)植物》一書中整理出75科187種外來(lái)非入侵植物,從每科中隨機(jī)抽取1種用于測(cè)試,天南星科(Araceae)的另一種藤本植物被額外選擇,最終選出76種受試外來(lái)非入侵植物。

    1.2 信息來(lái)源 植物特性信息主要來(lái)源于相關(guān)的期刊、著作和網(wǎng)絡(luò)數(shù)據(jù)庫(kù)資源,其中植物信息網(wǎng)站主要包括作物保護(hù)大全(crop protection compendium,CPC)、入侵物種大全(invasive species compendium,ISC)、中國(guó)國(guó)家有害生物檢疫信息系統(tǒng)、全球入侵物種數(shù)據(jù)庫(kù)(global invasive species database,GISD)、夏威夷-太平洋雜草風(fēng)險(xiǎn)評(píng)估數(shù)據(jù)庫(kù)(Hawaii-Pacific weed risk assessment)、Plant For A Future(PFAF)、種子信息數(shù)據(jù)庫(kù)(seed information database,SID)、植物染色體數(shù)目索引數(shù)據(jù)庫(kù)(index to plant chromosome numbers,IPCN)、熱帶有用植物數(shù)據(jù)庫(kù)(useful tropical plants database)、美國(guó)農(nóng)業(yè)部植物信息數(shù)據(jù)庫(kù)(USDA、NRCS、The PLANTS Database)、美國(guó)食品藥品管理局有毒植物數(shù)據(jù)庫(kù)(FDA poisonous plant database)、加拿大有毒植物信息系統(tǒng)(Canadian poisonous plants information system)、全球雜草大全(global compendium of weeds)、中國(guó)在線植物志等。

    1.3 美國(guó)雜草風(fēng)險(xiǎn)評(píng)估方法簡(jiǎn)介 采用經(jīng)典的分類方法——邏輯斯蒂回歸建立評(píng)判公式,其實(shí)質(zhì)為線性回歸模型,只是通過(guò)Sigmoid函數(shù)將該線性模型轉(zhuǎn)化為解決分類問(wèn)題的模型,其結(jié)果取值范圍為[0,1],即評(píng)判對(duì)象成為某一類別的概率。USWRA包括定殖/擴(kuò)散潛力(establishment/spread potential, ES)和后果影響潛力(impact potential, IMP)2個(gè)指標(biāo)體系。ES描述受試植物在風(fēng)險(xiǎn)地區(qū)逃逸、繁殖和擴(kuò)散的潛力,IMP描述受試植物對(duì)自然環(huán)境、人工生產(chǎn)系統(tǒng)(如農(nóng)田、果園等)、人為景觀(如花園、城市草坪等)造成直接或間接影響的潛力,二者總體均為判別式打分指標(biāo),主要由“是否問(wèn)題”構(gòu)成,其他還包括一些“多選擇問(wèn)題”。ES考慮了植物生物地理學(xué)、生物生態(tài)學(xué)、環(huán)境適應(yīng)能力等多方面屬性,具體包括23個(gè)評(píng)判因素;IMP考慮了植物對(duì)風(fēng)險(xiǎn)地區(qū)生物、環(huán)境、經(jīng)濟(jì)產(chǎn)品、人類、社會(huì)等多方面綜合影響,具體包括18個(gè)評(píng)判因素。評(píng)估時(shí)對(duì)各項(xiàng)因素進(jìn)行逐項(xiàng)評(píng)判打分求和,即分別求得ES和IMP的評(píng)估分?jǐn)?shù)。最終USWRA采用邏輯斯蒂回歸方法建立了關(guān)于ES和IMP的3個(gè)分類公式,分別用于計(jì)算受試植物成為非入侵植物(non-invader)、一般入侵植物(minor invader)和嚴(yán)重入侵植物(major invader)的概率(即Pnon-invader、Pminor invader和Pmajor invader):

    (1)

    (2)

    (3)

    當(dāng)Pnon-invader≥44.9%時(shí),受試植物被評(píng)估為非入侵植物;當(dāng)Pmajor invader≥38.8%時(shí),受試植物被評(píng)估為嚴(yán)重入侵植物;其他情況則被評(píng)定為一般入侵植物,需再評(píng)估(Koopetal., 2012;APHIS-PPQ, 2015)。

    此類測(cè)試研究主要采用拒絕率、接受率和再評(píng)估率對(duì)評(píng)估結(jié)果進(jìn)行統(tǒng)計(jì),即對(duì)嚴(yán)重入侵植物予以拒絕,非入侵植物予以接受,一般入侵植物則需再評(píng)估(Daehleretal.,2000; Daehleretal., 2004; Gordonetal., 2008)。

    1.4 USWRA評(píng)估效果的檢驗(yàn) 利用受試者工作特征曲線(receiver operating characteristic curve,ROC曲線)分別檢驗(yàn)USWRA對(duì)嚴(yán)重入侵植物和非入侵植物的評(píng)估效果,以及將一般入侵植物和嚴(yán)重入侵植物歸為入侵植物后USWRA對(duì)入侵植物和非入侵植物的評(píng)估效果。ROC曲線檢驗(yàn)在統(tǒng)計(jì)分析軟件SPSS22中操作完成。

    ROC曲線是在以假陽(yáng)性率為橫坐標(biāo)、真陽(yáng)性率為縱坐標(biāo)形成的坐標(biāo)系中根據(jù)評(píng)估結(jié)果繪制而成的一條曲線,線下的面積稱為AUC(area under curve)值,以AUC值大小反映診斷試驗(yàn)的價(jià)值。一般認(rèn)為AUC值為0.5~0.7時(shí)診斷價(jià)值較低,為0.7~0.9時(shí)診斷價(jià)值中等,大于0.9時(shí)診斷價(jià)值高(王運(yùn)生等, 2007)。

    1.5 不同生活型植物成為入侵植物的差異分析 1)不同生活型植物成為入侵性植物比例的差異比較。結(jié)合《中國(guó)入侵植物名錄》《中國(guó)外來(lái)生物入侵與檢疫規(guī)范》《中國(guó)外來(lái)植物》以及《中國(guó)植物志》(網(wǎng)絡(luò)版)的記錄,整理匯總受試外來(lái)植物的生活型。采用Iman等(1980)在Friedman非參數(shù)檢驗(yàn)方法基礎(chǔ)上(Friedman, 1937;1940)提出的新檢驗(yàn)方法(以下仍統(tǒng)稱該方法為“Friedman檢驗(yàn)”)檢驗(yàn)本研究選擇的不同生活型植物中成為入侵植物的比例是否存在顯著差異:

    (4)

    (5)

    式中:N為入侵層級(jí),本研究分為嚴(yán)重入侵和一般入侵2個(gè)層級(jí);k為生活型數(shù)量。

    當(dāng)FF>Fα[(k-1),(k-1)(N-1)](α最低為0.1水平)時(shí),認(rèn)為存在顯著差異。

    當(dāng)檢驗(yàn)結(jié)果顯示具有顯著差異(至少α=0.1水平上),采用Nemenyi(1963)方法進(jìn)行事后多重比較,但用于比較的關(guān)鍵差異值(critical difference,CD),按照Bonferroni-Dunn檢驗(yàn)提供的關(guān)鍵值(critical values)(表1)進(jìn)行計(jì)算(Dunn, 1961; Dem?ar, 2006),當(dāng)不同生活型植物成為2類入侵植物的平均秩值之差大于關(guān)鍵差異值CD時(shí),則二者生活型植物間存在顯著差異:

    (6)

    2)不同生活型植物綜合風(fēng)險(xiǎn)值大小的差異比較。USWRA中的線性回歸方程為“0.235 6×ES+0.601 9×IMP”,其計(jì)算結(jié)果被稱為綜合風(fēng)險(xiǎn)值(composite risk scores,以下簡(jiǎn)稱“風(fēng)險(xiǎn)值”),用以描述受試植物成為入侵植物的風(fēng)險(xiǎn)大小,其值越大,則表明植物越可能成為入侵植物(Koopetal.,2012)。據(jù)此公式計(jì)算本研究中各受試植物的風(fēng)險(xiǎn)值,利用單因素方差分析和LSD多重比較方法(操作在SPSS 22軟件中進(jìn)行)比較不同生活型植物之間風(fēng)險(xiǎn)值的差異,以初步了解我國(guó)外來(lái)植物的風(fēng)險(xiǎn)特點(diǎn)。

    表1 Bonferroni-Dunn檢驗(yàn)關(guān)鍵值

    2 結(jié)果與分析

    2.1 USWRA評(píng)估效果 131種受試外來(lái)植物的評(píng)估結(jié)果詳見附件1。結(jié)果顯示,所有受試植物定殖/擴(kuò)散潛力(ES)得分范圍為-9~17,其中嚴(yán)重入侵植物為4~17,平均為10.68;一般入侵植物為0~16,平均為6.67;非入侵植物為-9~10,平均為-1.47。所有受試植物后果影響潛力(IMP)得分范圍為1~2.5,其中嚴(yán)重入侵植物為1.1~2.5,平均為1.84;一般入侵植物為1.0~2.4,平均為1.66;非入侵植物為1~1.8,平均為1.06。

    單因素方差分析結(jié)果顯示,在ES得分上,非入侵植物、一般入侵植物和嚴(yán)重入侵植物呈現(xiàn)出兩兩極顯著差異;IMP得分上,非入侵植物分別與一般入侵植物和嚴(yán)重入侵植物呈現(xiàn)出極顯著差異,但一般入侵植物和嚴(yán)重入侵植物之間未呈現(xiàn)出顯著差異(表2)。

    對(duì)比USWRA評(píng)估受試植物得到的預(yù)測(cè)類別和受試植物在我國(guó)的已知身份,統(tǒng)計(jì)整理出USWRA對(duì)入侵植物和非入侵植物相應(yīng)的拒絕率、接受率和再評(píng)估率(表3)。結(jié)果顯示,USWRA對(duì)非入侵植物的識(shí)別效果最高,能準(zhǔn)確識(shí)別出84.2%的非入侵植物,同時(shí)確保沒(méi)有一種非入侵植物被拒絕;對(duì)嚴(yán)重入侵植物的識(shí)別效果次之,能準(zhǔn)確識(shí)別出68.4%的嚴(yán)重入侵植物,同時(shí)沒(méi)有一種嚴(yán)重入侵植物被接受;對(duì)一般入侵植物的識(shí)別效果較復(fù)雜,包含拒絕、接受和再評(píng)估3種結(jié)果,再評(píng)估率最高,達(dá)66.7%,拒絕率次之,為27.8%,同時(shí)有5.6%的一般入侵植物被接受。將一般入侵植物和嚴(yán)重入侵植物同時(shí)歸入入侵植物類,USWRA對(duì)入侵植物的評(píng)估效果顯示,雖然對(duì)入侵植物的拒絕率較低,僅為41.8%,但其接受的入侵植物并不多,僅為3.6%,只是無(wú)法評(píng)估的數(shù)量較大,占54.5%。

    表2 美國(guó)雜草風(fēng)險(xiǎn)評(píng)估方法對(duì)我國(guó)外來(lái)植物的評(píng)估結(jié)果分?jǐn)?shù)比較①

    ①上標(biāo)大寫字母A、B、C表示在0.01水平上差異顯著。The superscript capital letters A, B and C indicate significant differences at 0.01 level.

    表3 美國(guó)雜草風(fēng)險(xiǎn)評(píng)估方法對(duì)我國(guó)外來(lái)植物的評(píng)估效果

    2.2 ROC檢驗(yàn)結(jié)果 ROC檢驗(yàn)結(jié)果顯示,由USWRA對(duì)嚴(yán)重入侵植物和非入侵植物評(píng)估效果形成的曲線面積AUC值為0.990±0.007(SE),極顯著大于0.5(P<0.000 1);由USWRA對(duì)入侵植物和非入侵植物評(píng)估效果形成的曲線面積AUC值為0.963±0.014(SE),極顯著大于0.5(P<0.000 1)。具體如圖1所示。

    圖1 美國(guó)雜草風(fēng)險(xiǎn)評(píng)估方法對(duì)我國(guó)外來(lái)植物評(píng)估結(jié)果的ROC檢驗(yàn)

    2.3 不同生活型植物成為入侵植物的差異比較 1)不同生活型植物成為入侵植物的比例差異。各受試植物的生活型詳見附件1,生活型數(shù)量和比例見表4。本研究受試植物包括草本(66種)、喬木(29種)、灌木(19種)、藤本(12種)、竹類(3種)、肉質(zhì)(1種)、蕨類(1種)7種生活型,其中,草本類植物數(shù)量最多且占比最大,肉質(zhì)、蕨類和竹類植物數(shù)量最少且占比最低,喬木、灌木、藤本植物的數(shù)量和占比處于中間位置。

    表4 不同生活型植物中成為入侵植物的比例及秩值

    表5 不同生活型植物成為入侵植物比例間的多重比較結(jié)果

    2)不同生活型植物風(fēng)險(xiǎn)值大小的差異比較。草本、喬木、灌木和藤本4種生活型植物的風(fēng)險(xiǎn)值之間存在顯著差異(P=0.043),平均風(fēng)險(xiǎn)值從高到低依次為草本>藤本>灌木>喬木。草本和藤本植物的風(fēng)險(xiǎn)值之間無(wú)顯著差異,但二者風(fēng)險(xiǎn)值均顯著高于喬木(α=0.1),灌木和其他生活型植物的風(fēng)險(xiǎn)值之間未呈現(xiàn)出顯著差異。

    表6 不同生活型植物的風(fēng)險(xiǎn)值結(jié)果比較①

    ①a、b:0.1水平上差異顯著Significant difference at 0.1 level.

    3 討論與建議

    3.1 對(duì)嚴(yán)重入侵植物和非入侵植物評(píng)估效果的分析 本次測(cè)試中,USWRA正確拒絕了68.4%的嚴(yán)重入侵植物,沒(méi)有錯(cuò)誤接受一種嚴(yán)重入侵植物,其余31.6%需再評(píng)估。嚴(yán)重入侵植物的正確拒絕率略低但接近70%,且不會(huì)導(dǎo)致嚴(yán)重入侵植物被接受,因此綜合認(rèn)為USWRA對(duì)嚴(yán)重入侵植物具有有效的評(píng)估效果。USWRA對(duì)嚴(yán)重入侵植物拒絕率略低,主要是因?yàn)闇y(cè)試將外來(lái)植物視為尚未進(jìn)入我國(guó)境內(nèi),未評(píng)價(jià)在其他地區(qū)的入侵歷史指標(biāo)。本研究中受試外來(lái)植物的傳入時(shí)間基本在19世紀(jì)40年代至20世紀(jì)50年代之間,很難查詢到這一傳入時(shí)間或之前植物在國(guó)外其他地區(qū)的入侵歷史信息,測(cè)試時(shí)未評(píng)價(jià)USWRA中的ES-1、ES-3、IMP-N6、IMP-A4和IMP-P6指標(biāo)(Koopetal., 2012);而植物在其他地區(qū)的入侵歷史,甚至可以單獨(dú)用于評(píng)估(Mack, 1996),所以這一信息在評(píng)估方法中占有一定比重。USWRA中ES指標(biāo)分?jǐn)?shù)為-25~32分,其中其他地區(qū)的入侵歷史分?jǐn)?shù)為0~6分,IMP指標(biāo)分?jǐn)?shù)為1~5.1分,其中其他地區(qū)的入侵歷史分?jǐn)?shù)為0~1.6分。ES得分對(duì)USWRA的評(píng)估結(jié)果影響較大(Koopetal., 2012),本次測(cè)試中嚴(yán)重入侵植物ES平均得分為10.68,IMP平均得分為1.84,據(jù)此計(jì)算出植物成為嚴(yán)重入侵植物的平均概率為37.5%(低于閾值38.8%),由此推算ES僅需增加0.42分,就能保證成為嚴(yán)重入侵植物的概率超過(guò)閾值38.8%。而在USWRA的ES指標(biāo)中,如果外來(lái)植物在其他地區(qū)有發(fā)生歷史,則ES將賦值為5分(見ES-1的賦值方式)。Dawson等(2009)研究發(fā)現(xiàn),78%能自主擴(kuò)散的外來(lái)植物,在其他地區(qū)具入侵歷史,若據(jù)此比例推算(本研究受試的嚴(yán)重入侵植物不僅能自主擴(kuò)散,而且實(shí)際造成了嚴(yán)重危害),當(dāng)評(píng)估植物在其他地區(qū)的入侵歷史指標(biāo)后,USWRA對(duì)嚴(yán)重入侵植物的正確拒絕率將至少上升至78%,同時(shí)大幅度降低再評(píng)估率,并保持零錯(cuò)誤接受率,總體評(píng)估效果得到較高提升。因此在實(shí)際應(yīng)用時(shí),當(dāng)獲取到植物在其他地區(qū)的入侵歷史信息,其評(píng)估效果將更加高效。

    USWRA正確接受了84.2%的非入侵植物,沒(méi)有錯(cuò)誤拒絕一種非入侵植物,其余15.8%需再評(píng)估,據(jù)此認(rèn)為USWRA對(duì)非入侵植物的評(píng)估效果較好。對(duì)于非入侵植物,本研究認(rèn)為未評(píng)判植物在其他地區(qū)的入侵歷史,幾乎不會(huì)對(duì)評(píng)估效果產(chǎn)生影響。本次選擇的76種非入侵植物中,60種屬于觀賞類植物,其他或?yàn)槭卟嘶驗(yàn)楣麡浠驗(yàn)樗幱玫?,人工馴化程度均較高,能形成自我繁殖種群的比例較少,這大大降低了其曾在其他地區(qū)成為入侵植物的可能性。Dawson等(2009)研究發(fā)現(xiàn),在東非雨林地帶能“生存”(surviving)下來(lái)的外來(lái)植物在其他地區(qū)具有入侵歷史者占18%,由此推斷那些不能形成單獨(dú)種群的外來(lái)植物,具有入侵歷史的比例會(huì)更低。

    3.2 對(duì)一般入侵植物評(píng)估效果的分析 USWRA對(duì)一般入侵植物的評(píng)估結(jié)果較低,拒絕率僅為27.8%,而未識(shí)別率達(dá)到66.7%,錯(cuò)誤率為5.6%,表現(xiàn)出未識(shí)別率較高、拒絕率較低的情況。

    首先,假設(shè)一般入侵植物在實(shí)際中應(yīng)為嚴(yán)重入侵植物而予以拒絕。當(dāng)將一般入侵植物視為嚴(yán)重入侵植物,并與嚴(yán)重入侵植物合并為入侵植物后,結(jié)果顯示USWRA對(duì)入侵植物的拒絕率較低,為42%。造成拒絕率低的原因,可能有以下2種:一是前述由于受到植物其他地區(qū)入侵歷史缺失的影響;二可能是因?yàn)樵撘话闳肭种参飳?shí)際中應(yīng)為非入侵植物,而原劃分標(biāo)準(zhǔn)將其錯(cuò)誤地歸為一般入侵植物(馬金雙, 2013)。在《中國(guó)入侵植物名錄》中,外來(lái)入侵植物第4個(gè)等級(jí)為一般入侵類,其定義的主要依據(jù)是“危害不大或不明顯”,本研究認(rèn)為這一評(píng)判規(guī)則容易引起主觀偏差,從而導(dǎo)致非入侵植物被歸入到一般入侵植物類中,或者將一般入侵植物被歸入到非入侵植物類中。一般入侵植物和非入侵植物之間不是離散型,而是一個(gè)連續(xù)狀態(tài),二者之間具有共同的特征(Gordonetal., 2008),由此導(dǎo)致評(píng)估方法對(duì)一般入侵植物無(wú)法判別。從另一方面看,本次測(cè)試中無(wú)法識(shí)別的一般入侵植物,提示《中國(guó)入侵植物名錄》中對(duì)于一般入侵植物(IV)的歸類或許需要重新調(diào)整,而閆小玲等(2014)也認(rèn)為該書中入侵植物的等級(jí)應(yīng)根據(jù)實(shí)際情況做出動(dòng)態(tài)調(diào)整。

    其次,一般入侵植物與非入侵植物以及嚴(yán)重入侵植物之間的重疊性,導(dǎo)致無(wú)法真正定義一般入侵植物(Gordonetal., 2008),即一般入侵植物應(yīng)為嚴(yán)重入侵植物而予以拒絕,還是應(yīng)為非入侵植物而被接受,尚無(wú)法給出準(zhǔn)確的評(píng)判。因此,雖然本次評(píng)估結(jié)果顯示一般入侵植物拒絕率低,但這并不代表USWRA的準(zhǔn)確率低,而是由于一般入侵植物的特殊性造成的。如果以對(duì)一般入侵植物的評(píng)估效率認(rèn)為USWRA效率偏低,會(huì)導(dǎo)致結(jié)論會(huì)有所偏差。因此本研究認(rèn)為,應(yīng)更加注重考察USWRA對(duì)嚴(yán)重入侵植物和非入侵植物的評(píng)估效果。

    3.3 不同生活型植物成為入侵植物的比例 本研究所選擇的4種不同生活型植物,非參數(shù)檢驗(yàn)結(jié)果認(rèn)為各自成為入侵植物的比例存在差異,只是進(jìn)一步多重比較結(jié)果無(wú)法找出具有顯著差異者。分析原因認(rèn)為主要與樣本選擇方式有關(guān),本研究主要致力于檢驗(yàn)USWRA對(duì)我國(guó)外來(lái)植物的評(píng)估效果,未充分考慮生活型的樣本數(shù)量,灌木(19種)和藤本(12種)植物數(shù)量相對(duì)較少;但另一方面,USWRA對(duì)不同生活型植物的風(fēng)險(xiǎn)值結(jié)果表明,喬木植物在我國(guó)成為入侵植物的風(fēng)險(xiǎn)值顯著低于草本植物。由此初步表明,我國(guó)面臨的植物入侵風(fēng)險(xiǎn)主要可能來(lái)自于草本植物,與閆小玲等(2014)研究結(jié)果相似,即我國(guó)外來(lái)入侵植物中超過(guò)80%為草本植物,而喬木植物僅占4%。因此,USWRA的綜合風(fēng)險(xiǎn)值也能較好指示植物成為入侵植物的風(fēng)險(xiǎn)大小,從另一側(cè)面表明該方法具有較高的評(píng)判效果。

    3.4 USWRA在我國(guó)的應(yīng)用與發(fā)展建議 USWRA的指標(biāo)主要為植物特性(Koopetal., 2012),其實(shí)質(zhì)屬于風(fēng)險(xiǎn)評(píng)估中的歸類方法,用以識(shí)別出潛在的入侵植物和非入侵植物(ISPM 2, 2007)。因此在實(shí)際應(yīng)用時(shí),不能僅根據(jù)USWRA評(píng)估結(jié)果即對(duì)受試植物做出允許入境的管理決定,還應(yīng)重點(diǎn)評(píng)估植物攜帶其他有害生物的風(fēng)險(xiǎn)。

    我國(guó)林木引種需求較高(趙宇翔, 2017),但主動(dòng)引入往往造成植物入侵。目前我國(guó)6種重要的林業(yè)外來(lái)入侵植物中有4種屬于主動(dòng)引種傳入(趙宇翔等, 2015),因此有效評(píng)估擬引入林木自身是否可能為入侵植物,將降低入侵植物的傳入風(fēng)險(xiǎn),提高風(fēng)險(xiǎn)評(píng)估工作的效率。本次測(cè)試結(jié)果認(rèn)為,USWRA可以作為一種參考方法應(yīng)用于我國(guó)林木引種的風(fēng)險(xiǎn)評(píng)估工作。為更好地應(yīng)用USWRA以及促進(jìn)我國(guó)入侵植物風(fēng)險(xiǎn)評(píng)估方法的研究,提出以下建議:

    1)應(yīng)用USWRA時(shí),可對(duì)USWRA中的ES-8“是否為禾本科草本植物”和ES-17“自然傳播方式的數(shù)量”進(jìn)行修改(Koopetal., 2012)。建議在ES-8中加入菊科草本植物,因?yàn)榫湛祁惾肭种参镌谖覈?guó)發(fā)生數(shù)量最多且危害嚴(yán)重(閆小玲等, 2014)。對(duì)于ES-17指標(biāo),其賦值方式是按自然傳播方式的總數(shù)量賦分(即0種=-4;1種=-2;2種=0;3種=2;4或5種=4),這就要求評(píng)估者必須全部掌握受試植物5種自然傳播方式的信息才能做出準(zhǔn)確評(píng)判,但有些傳播方式難于查找甚至無(wú)相關(guān)信息,如動(dòng)物體內(nèi)傳播等。因此建議,可在保持ES-17指標(biāo)分值不變(4分)的情況下,對(duì)5種自然傳播方式實(shí)施逐項(xiàng)賦值(平均或權(quán)重),相加后得到ES-17的評(píng)估分?jǐn)?shù)。

    2)加強(qiáng)開展多類群物種的特性比較研究。從受試植物來(lái)看,我國(guó)外來(lái)草本植物占據(jù)大多數(shù)[與Weber等(2008)的研究結(jié)果相似],且草本和藤本植物的風(fēng)險(xiǎn)值較高,喬木植物的風(fēng)險(xiǎn)值相對(duì)較低,初步認(rèn)為草本和藤本植物對(duì)我國(guó)的入侵風(fēng)險(xiǎn)最高,而喬木植物入侵風(fēng)險(xiǎn)較低。這與美國(guó)的入侵植物生活型特點(diǎn)不盡相同,因?yàn)閺腢SWRA的指標(biāo)可知(ES-5、7、9),藤本、水生和木本植物對(duì)美國(guó)具有較高的入侵風(fēng)險(xiǎn)(Koopetal., 2012)。因此建議,加強(qiáng)開展多類群物種間特性的比較研究,如外來(lái)入侵植物與外來(lái)非入侵植物、外來(lái)入侵植物與本土非入侵植物等(Py?eketal.,2007;2008),系統(tǒng)篩選出與我國(guó)入侵植物高度相關(guān)的特性,建立適合我國(guó)特點(diǎn)的入侵植物風(fēng)險(xiǎn)評(píng)估歸類方法。

    3)對(duì)于USWRA評(píng)定為再評(píng)估的植物,可利用相關(guān)的再評(píng)估規(guī)則進(jìn)行最終判斷,如Koop等(2012)為USWRA補(bǔ)充建立的再評(píng)估規(guī)則,或Daehler等(2004)建立的判別樹再評(píng)估規(guī)則。

    猜你喜歡
    生活型植物評(píng)估
    明堂山主要藥用植物資源調(diào)查分析
    論數(shù)學(xué)教學(xué)中有效運(yùn)用數(shù)學(xué)模型
    成才之路(2017年21期)2017-07-27 20:52:40
    哦,不怕,不怕
    寬闊水國(guó)家級(jí)自然保護(hù)區(qū)珙桐群落特征及演替趨勢(shì)分析
    將植物穿身上
    植物罷工啦?
    關(guān)于初中數(shù)學(xué)作業(yè)布置的幾點(diǎn)思考
    評(píng)估依據(jù)
    植物也瘋狂
    立法后評(píng)估:且行且盡善
    浙江人大(2014年5期)2014-03-20 16:20:25
    av视频在线观看入口| 久久国内精品自在自线图片| 一二三四中文在线观看免费高清| 国产精品国产三级专区第一集| 中文资源天堂在线| 好男人在线观看高清免费视频| 高清av免费在线| 国产高清有码在线观看视频| 国产精品国产三级国产专区5o | 国产精品福利在线免费观看| 久久久久久大精品| 亚洲人与动物交配视频| 免费人成在线观看视频色| 精品酒店卫生间| 久久久亚洲精品成人影院| www.色视频.com| 91在线精品国自产拍蜜月| 水蜜桃什么品种好| 久久这里有精品视频免费| 精品不卡国产一区二区三区| 久久这里只有精品中国| av卡一久久| 久久精品国产99精品国产亚洲性色| 久久这里有精品视频免费| 最近中文字幕2019免费版| 国产一区二区在线观看日韩| 两个人视频免费观看高清| 秋霞在线观看毛片| 91狼人影院| 久久热精品热| 一级毛片电影观看 | 国产毛片a区久久久久| 久久热精品热| 国产高清国产精品国产三级 | 久久99精品国语久久久| 成人漫画全彩无遮挡| 简卡轻食公司| 免费一级毛片在线播放高清视频| 国产三级中文精品| 最近视频中文字幕2019在线8| 欧美色视频一区免费| 国产大屁股一区二区在线视频| 岛国毛片在线播放| 久久久久精品久久久久真实原创| 丝袜美腿在线中文| 色尼玛亚洲综合影院| 亚洲欧美日韩东京热| 老司机影院毛片| 天堂影院成人在线观看| 成人美女网站在线观看视频| 国产视频首页在线观看| 国产精品久久久久久久电影| 蜜桃久久精品国产亚洲av| 插阴视频在线观看视频| 舔av片在线| 国产伦精品一区二区三区视频9| 在线观看av片永久免费下载| 高清午夜精品一区二区三区| 啦啦啦观看免费观看视频高清| 欧美极品一区二区三区四区| 国产精品伦人一区二区| 亚洲精品456在线播放app| 国产精品人妻久久久久久| 边亲边吃奶的免费视频| 男女国产视频网站| 亚洲综合色惰| 老司机影院成人| 一级爰片在线观看| 天天躁夜夜躁狠狠久久av| 久久精品影院6| 亚洲最大成人手机在线| av专区在线播放| 九草在线视频观看| 国产大屁股一区二区在线视频| 夜夜爽夜夜爽视频| 国产精品一及| 少妇熟女aⅴ在线视频| 国产精华一区二区三区| 最近最新中文字幕免费大全7| 我的老师免费观看完整版| 日韩一区二区视频免费看| 日韩欧美国产在线观看| 美女高潮的动态| 国产真实乱freesex| 亚洲成人精品中文字幕电影| 国产av码专区亚洲av| av黄色大香蕉| 特大巨黑吊av在线直播| 可以在线观看毛片的网站| 中文亚洲av片在线观看爽| 国产一区亚洲一区在线观看| 青春草国产在线视频| 尾随美女入室| 国产伦精品一区二区三区视频9| 日韩欧美精品免费久久| 国产激情偷乱视频一区二区| 亚洲精品成人久久久久久| 女人十人毛片免费观看3o分钟| 亚洲av不卡在线观看| 岛国毛片在线播放| 国产男人的电影天堂91| 青春草国产在线视频| 国产成人freesex在线| 色噜噜av男人的天堂激情| 亚洲欧美日韩卡通动漫| 黄色欧美视频在线观看| 国产三级中文精品| 性插视频无遮挡在线免费观看| 精品久久久久久久久久久久久| 高清午夜精品一区二区三区| 有码 亚洲区| 午夜激情欧美在线| 三级经典国产精品| 变态另类丝袜制服| 中文字幕制服av| 亚洲性久久影院| 日韩高清综合在线| 精品一区二区三区人妻视频| 国产免费福利视频在线观看| 亚洲精品日韩在线中文字幕| 日韩一本色道免费dvd| 99视频精品全部免费 在线| 亚洲成av人片在线播放无| 永久网站在线| 久久99蜜桃精品久久| 精品一区二区免费观看| 久久久久久久午夜电影| 日本av手机在线免费观看| 日韩亚洲欧美综合| 免费在线观看成人毛片| 国产69精品久久久久777片| 日本一本二区三区精品| 91精品一卡2卡3卡4卡| 亚洲无线观看免费| 久久精品91蜜桃| 色哟哟·www| 国产视频内射| 国产精品久久久久久久电影| 久久久久久大精品| 亚洲天堂国产精品一区在线| 只有这里有精品99| 99久久无色码亚洲精品果冻| 热99re8久久精品国产| videossex国产| 欧美日韩综合久久久久久| 亚洲国产高清在线一区二区三| 国产精品一区二区性色av| 中文天堂在线官网| 久久久久精品久久久久真实原创| av又黄又爽大尺度在线免费看 | 国产一级毛片在线| 亚洲中文字幕一区二区三区有码在线看| 色哟哟·www| 婷婷色av中文字幕| 国产精品一区二区三区四区免费观看| 汤姆久久久久久久影院中文字幕 | 天堂av国产一区二区熟女人妻| 国产精品综合久久久久久久免费| 国产成人a∨麻豆精品| 天堂网av新在线| 婷婷色麻豆天堂久久 | 插逼视频在线观看| 少妇裸体淫交视频免费看高清| 免费黄网站久久成人精品| 人体艺术视频欧美日本| 青青草视频在线视频观看| 亚洲精品一区蜜桃| 日韩一本色道免费dvd| 两个人视频免费观看高清| 男女那种视频在线观看| 精华霜和精华液先用哪个| 禁无遮挡网站| 国产成人freesex在线| 国产精品国产三级国产av玫瑰| 乱系列少妇在线播放| 亚洲在久久综合| 美女黄网站色视频| 一级黄片播放器| 欧美不卡视频在线免费观看| 亚洲美女视频黄频| 久久精品夜夜夜夜夜久久蜜豆| 精品国产一区二区三区久久久樱花 | 亚洲人与动物交配视频| 啦啦啦啦在线视频资源| 日韩一区二区视频免费看| 人体艺术视频欧美日本| 国产成人免费观看mmmm| 精品无人区乱码1区二区| 看片在线看免费视频| 国产精品电影一区二区三区| 国产成人a∨麻豆精品| 日韩亚洲欧美综合| av卡一久久| 成人二区视频| 91狼人影院| 亚洲性久久影院| 国产精品1区2区在线观看.| 哪个播放器可以免费观看大片| 亚洲av一区综合| 99热6这里只有精品| 91精品一卡2卡3卡4卡| 久久99蜜桃精品久久| 国产精品爽爽va在线观看网站| 99久久精品国产国产毛片| 蜜桃久久精品国产亚洲av| 精品国内亚洲2022精品成人| 亚洲成人av在线免费| 欧美日韩一区二区视频在线观看视频在线 | 搡女人真爽免费视频火全软件| 中国国产av一级| 色尼玛亚洲综合影院| 中文字幕av成人在线电影| 一二三四中文在线观看免费高清| 日日摸夜夜添夜夜添av毛片| 精品午夜福利在线看| 国产精品久久久久久久久免| 亚洲精品乱码久久久久久按摩| 久久韩国三级中文字幕| 国产伦在线观看视频一区| 国产成人91sexporn| 精品无人区乱码1区二区| 夫妻性生交免费视频一级片| 69人妻影院| 国产精品电影一区二区三区| 国产成人免费观看mmmm| 精品久久久久久成人av| 国产午夜精品久久久久久一区二区三区| 97热精品久久久久久| 91精品一卡2卡3卡4卡| 最近最新中文字幕免费大全7| 色视频www国产| 永久网站在线| 日韩高清综合在线| 女人久久www免费人成看片 | 国产精品.久久久| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| 亚洲欧美精品专区久久| 简卡轻食公司| 亚洲第一区二区三区不卡| 午夜亚洲福利在线播放| 18+在线观看网站| av免费在线看不卡| 午夜福利成人在线免费观看| 亚洲乱码一区二区免费版| 男人和女人高潮做爰伦理| 伊人久久精品亚洲午夜| 青春草亚洲视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 欧美精品一区二区大全| 亚洲国产欧美在线一区| 欧美日韩综合久久久久久| 少妇被粗大猛烈的视频| 亚洲18禁久久av| 一级二级三级毛片免费看| 国产老妇伦熟女老妇高清| 国产色爽女视频免费观看| 久久国产乱子免费精品| 麻豆精品久久久久久蜜桃| 永久网站在线| 18禁裸乳无遮挡免费网站照片| 亚洲无线观看免费| 国产视频首页在线观看| a级一级毛片免费在线观看| www.色视频.com| 成人二区视频| 久久99蜜桃精品久久| av女优亚洲男人天堂| 久久久久久久久久久免费av| 欧美性猛交黑人性爽| 免费黄网站久久成人精品| 久久久久久九九精品二区国产| 超碰av人人做人人爽久久| 亚洲欧美日韩高清专用| 国产精品一区二区三区四区久久| 日本免费a在线| 精品久久久久久成人av| 边亲边吃奶的免费视频| АⅤ资源中文在线天堂| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 免费播放大片免费观看视频在线观看 | 精品久久久久久电影网 | 亚洲成色77777| 日韩制服骚丝袜av| 亚洲四区av| 老女人水多毛片| 91久久精品电影网| 九九热线精品视视频播放| 最近的中文字幕免费完整| 久久这里有精品视频免费| 久久欧美精品欧美久久欧美| 在线观看美女被高潮喷水网站| 日日干狠狠操夜夜爽| 国产91av在线免费观看| 日日摸夜夜添夜夜添av毛片| 男人舔奶头视频| 一区二区三区四区激情视频| 一级av片app| 亚洲国产高清在线一区二区三| 免费搜索国产男女视频| 亚洲国产精品合色在线| 狠狠狠狠99中文字幕| 久久亚洲国产成人精品v| 亚洲国产欧洲综合997久久,| 国产乱人偷精品视频| 日本色播在线视频| 一级av片app| 日产精品乱码卡一卡2卡三| 人妻少妇偷人精品九色| 成人毛片60女人毛片免费| 最近手机中文字幕大全| 亚洲国产高清在线一区二区三| 乱人视频在线观看| 国产亚洲精品av在线| 国产私拍福利视频在线观看| 欧美成人午夜免费资源| 草草在线视频免费看| 成人二区视频| 97人妻精品一区二区三区麻豆| 欧美xxxx性猛交bbbb| 精品一区二区三区人妻视频| 亚洲精品国产av成人精品| 国产精品一区二区在线观看99 | 18禁在线无遮挡免费观看视频| 国产免费男女视频| 免费观看性生交大片5| 欧美日韩国产亚洲二区| 亚洲av电影在线观看一区二区三区 | 国产片特级美女逼逼视频| 美女脱内裤让男人舔精品视频| 丝袜美腿在线中文| 不卡视频在线观看欧美| av视频在线观看入口| 国产在线男女| 久久久精品94久久精品| 天堂av国产一区二区熟女人妻| 亚洲av成人av| av专区在线播放| 变态另类丝袜制服| 国产在线男女| 中国国产av一级| 人人妻人人看人人澡| 男女边吃奶边做爰视频| 小蜜桃在线观看免费完整版高清| 色哟哟·www| 一级毛片aaaaaa免费看小| 国产成年人精品一区二区| 久久久a久久爽久久v久久| 亚洲综合色惰| 国产一区有黄有色的免费视频 | 激情 狠狠 欧美| 赤兔流量卡办理| 国产一区二区亚洲精品在线观看| 免费观看精品视频网站| 国产一区二区亚洲精品在线观看| 成人一区二区视频在线观看| 五月玫瑰六月丁香| 秋霞伦理黄片| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 国产成人精品婷婷| a级一级毛片免费在线观看| 有码 亚洲区| 日韩一区二区视频免费看| 精品久久久久久久久av| 18+在线观看网站| 国产精品爽爽va在线观看网站| 中文字幕av在线有码专区| 久久精品综合一区二区三区| 欧美高清性xxxxhd video| 免费av不卡在线播放| 春色校园在线视频观看| 日韩一区二区视频免费看| 亚洲精品aⅴ在线观看| av视频在线观看入口| 内地一区二区视频在线| 熟女电影av网| 99九九线精品视频在线观看视频| 麻豆av噜噜一区二区三区| 熟妇人妻久久中文字幕3abv| 男人舔女人下体高潮全视频| 国语对白做爰xxxⅹ性视频网站| 最近最新中文字幕大全电影3| 成人国产麻豆网| 久久草成人影院| 日日啪夜夜撸| 性插视频无遮挡在线免费观看| h日本视频在线播放| 一区二区三区免费毛片| 亚洲欧美精品自产自拍| 亚洲国产色片| 日产精品乱码卡一卡2卡三| 日本黄色片子视频| 色噜噜av男人的天堂激情| 高清毛片免费看| 99热网站在线观看| 欧美一级a爱片免费观看看| 高清av免费在线| 成人午夜精彩视频在线观看| 嫩草影院精品99| av在线亚洲专区| 亚洲av一区综合| 九草在线视频观看| 一级毛片我不卡| 美女脱内裤让男人舔精品视频| 一个人免费在线观看电影| 午夜视频国产福利| 久久国产乱子免费精品| 看十八女毛片水多多多| 噜噜噜噜噜久久久久久91| 国产久久久一区二区三区| 亚洲精品乱码久久久久久按摩| 久久精品久久精品一区二区三区| 美女被艹到高潮喷水动态| 日韩制服骚丝袜av| 免费在线观看成人毛片| 欧美日韩一区二区视频在线观看视频在线 | 国产精品综合久久久久久久免费| 丝袜美腿在线中文| 秋霞伦理黄片| 丰满人妻一区二区三区视频av| 国产片特级美女逼逼视频| 岛国毛片在线播放| 国产视频首页在线观看| 老女人水多毛片| 深夜a级毛片| 精品欧美国产一区二区三| av黄色大香蕉| 久久99蜜桃精品久久| 亚洲欧美日韩高清专用| 99久久九九国产精品国产免费| 国产一区二区亚洲精品在线观看| 亚洲国产精品久久男人天堂| 亚洲自拍偷在线| 最近2019中文字幕mv第一页| 国产av码专区亚洲av| 99国产精品一区二区蜜桃av| 欧美+日韩+精品| 久久人人爽人人爽人人片va| 国产成人精品一,二区| 日本猛色少妇xxxxx猛交久久| 日韩在线高清观看一区二区三区| 校园人妻丝袜中文字幕| 美女脱内裤让男人舔精品视频| 日韩亚洲欧美综合| 欧美日本亚洲视频在线播放| 午夜福利成人在线免费观看| 国产精品,欧美在线| 国产亚洲av嫩草精品影院| 身体一侧抽搐| 中文字幕精品亚洲无线码一区| 亚洲欧美日韩无卡精品| 亚洲精品一区蜜桃| 麻豆国产97在线/欧美| 国产伦一二天堂av在线观看| 69av精品久久久久久| 插逼视频在线观看| 欧美xxxx性猛交bbbb| 少妇人妻一区二区三区视频| 精品久久久久久久久久久久久| 长腿黑丝高跟| 搞女人的毛片| 亚洲精品国产av成人精品| 午夜福利高清视频| 如何舔出高潮| 黑人高潮一二区| 日韩欧美在线乱码| 一个人观看的视频www高清免费观看| 日本免费在线观看一区| 成人性生交大片免费视频hd| 国产午夜精品久久久久久一区二区三区| 亚洲不卡免费看| 色5月婷婷丁香| 国产免费视频播放在线视频 | 国产伦精品一区二区三区视频9| 天堂网av新在线| 国产精品精品国产色婷婷| 成年版毛片免费区| 在线免费十八禁| 亚洲成人精品中文字幕电影| 中文在线观看免费www的网站| 有码 亚洲区| 日本熟妇午夜| 99久久精品国产国产毛片| 神马国产精品三级电影在线观看| 九九爱精品视频在线观看| 国产男人的电影天堂91| 特级一级黄色大片| 午夜免费激情av| 久久久午夜欧美精品| 老司机福利观看| 国产精品人妻久久久影院| 91精品一卡2卡3卡4卡| 日本爱情动作片www.在线观看| 欧美日韩综合久久久久久| 国产黄色视频一区二区在线观看 | 精品久久久噜噜| 青春草亚洲视频在线观看| 免费av不卡在线播放| 能在线免费观看的黄片| 久久精品久久久久久噜噜老黄 | 超碰av人人做人人爽久久| 小说图片视频综合网站| 日韩 亚洲 欧美在线| 免费观看a级毛片全部| 亚洲精品乱码久久久v下载方式| 亚洲av成人av| 人妻制服诱惑在线中文字幕| 久久久久网色| 亚洲国产精品sss在线观看| 观看免费一级毛片| 免费大片18禁| 日本-黄色视频高清免费观看| 久久99热6这里只有精品| 99热网站在线观看| 亚洲成av人片在线播放无| 免费大片18禁| 久久久久国产网址| 丝袜美腿在线中文| 女人久久www免费人成看片 | 大香蕉久久网| 国产精品久久视频播放| 国语自产精品视频在线第100页| 久久韩国三级中文字幕| 网址你懂的国产日韩在线| 搡老妇女老女人老熟妇| 91狼人影院| 国产亚洲91精品色在线| 色尼玛亚洲综合影院| 美女国产视频在线观看| 欧美成人午夜免费资源| 26uuu在线亚洲综合色| 亚洲av不卡在线观看| 男女下面进入的视频免费午夜| 精品久久久久久成人av| 日本wwww免费看| 国产又黄又爽又无遮挡在线| 亚州av有码| 亚洲在线自拍视频| 亚洲精品色激情综合| 十八禁国产超污无遮挡网站| 亚洲av免费在线观看| 欧美97在线视频| 色综合站精品国产| 一个人看的www免费观看视频| 午夜爱爱视频在线播放| 亚洲美女视频黄频| kizo精华| 一本久久精品| 精品人妻视频免费看| 国产日韩欧美在线精品| 亚洲欧美日韩卡通动漫| 日韩人妻高清精品专区| 亚洲国产精品国产精品| 欧美性猛交╳xxx乱大交人| 夜夜看夜夜爽夜夜摸| 成人二区视频| 少妇的逼好多水| 女人久久www免费人成看片 | 深夜a级毛片| 欧美97在线视频| 亚洲国产精品成人综合色| 高清午夜精品一区二区三区| 久久久国产成人免费| 日韩一本色道免费dvd| 99热这里只有是精品在线观看| 国产精品福利在线免费观看| 亚洲av不卡在线观看| 精品久久久久久电影网 | 国产一区亚洲一区在线观看| av专区在线播放| 18禁在线播放成人免费| 免费看光身美女| 99九九线精品视频在线观看视频| 久久亚洲国产成人精品v| 视频中文字幕在线观看| 国产精品乱码一区二三区的特点| 丝袜美腿在线中文| 插逼视频在线观看| 午夜爱爱视频在线播放| 久久99精品国语久久久| 欧美丝袜亚洲另类| 哪个播放器可以免费观看大片| 亚洲国产精品久久男人天堂| 久久久成人免费电影| 亚洲av成人精品一区久久| 国产精品日韩av在线免费观看| 国产亚洲精品av在线| 欧美zozozo另类| 精品熟女少妇av免费看| 午夜福利成人在线免费观看| 亚洲精品国产成人久久av| 最近中文字幕2019免费版| 国产精品伦人一区二区| 国产一区二区亚洲精品在线观看| 欧美区成人在线视频| 男人狂女人下面高潮的视频| 亚洲国产欧美人成| 中文亚洲av片在线观看爽| 极品教师在线视频| 久久久精品94久久精品| 日本免费在线观看一区| 欧美另类亚洲清纯唯美| 舔av片在线| 两个人的视频大全免费| 又爽又黄a免费视频| 最近中文字幕高清免费大全6| 又粗又爽又猛毛片免费看| 国产一区有黄有色的免费视频 | 国内精品宾馆在线| 亚洲欧美成人综合另类久久久 | 一边摸一边抽搐一进一小说| 日韩一本色道免费dvd| 中文字幕制服av|