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

    北美刺龍葵在中國的適生區(qū)預(yù)測

    2016-04-08 16:09:14冼曉青萬方浩
    生物安全學(xué)報(bào) 2016年2期
    關(guān)鍵詞:適生區(qū)

    王 瑞, 冼曉青, 萬方浩,2*

    1中國農(nóng)業(yè)科學(xué)院植物保護(hù)研究所,植物病蟲害生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,北京 100193;2青島農(nóng)業(yè)大學(xué)農(nóng)學(xué)與植物保護(hù)學(xué)院,山東 青島 266109

    ?

    北美刺龍葵在中國的適生區(qū)預(yù)測

    王瑞1, 冼曉青1, 萬方浩1,2*

    1中國農(nóng)業(yè)科學(xué)院植物保護(hù)研究所,植物病蟲害生物學(xué)國家重點(diǎn)實(shí)驗(yàn)室,北京 100193;2青島農(nóng)業(yè)大學(xué)農(nóng)學(xué)與植物保護(hù)學(xué)院,山東 青島 266109

    摘要:【背景】北美刺龍葵是一種全球廣泛分布的惡性雜草,已被列入我國進(jìn)境檢疫性有害生物名單。近年來北美刺龍葵不斷隨進(jìn)口貨物傳入我國,明確其傳入途徑和適生區(qū)對(duì)控制其入侵具有重要意義?!痉椒ā坎捎肎IS、空間統(tǒng)計(jì)學(xué)、Maxent生態(tài)位模型等方法分析了北美刺龍葵的傳入途徑與潛在分布區(qū),并通過ROC分析法對(duì)模型進(jìn)行檢驗(yàn)?!窘Y(jié)果】跨區(qū)域農(nóng)產(chǎn)品貿(mào)易是北美刺龍葵全球擴(kuò)散的驅(qū)動(dòng)力與傳入我國的主要途徑。生態(tài)模型預(yù)測結(jié)果表明,北美刺龍葵在我國具有廣闊的適生區(qū),除黑龍江、吉林、內(nèi)蒙古、青海、甘肅、西藏、四川西北部以外的區(qū)域都是其在我國的適生區(qū),其中高風(fēng)險(xiǎn)區(qū)主要集中在東部和南部沿海、西南邊境和新疆的部分地區(qū)。AUC值為0.789,表明本研究建立的Maxent模型的預(yù)測能力較強(qiáng),能夠很好地?cái)M合物種已知分布的環(huán)境生態(tài)位。【結(jié)論與意義】北美刺龍葵在我國的傳入風(fēng)險(xiǎn)極高?;诒泵来听埧谖覈闹饕獋魅胪緩脚c潛在擴(kuò)散媒介的時(shí)空分布,劃定了重點(diǎn)監(jiān)測的區(qū)域,建議對(duì)適生區(qū)內(nèi)極易傳入的高風(fēng)險(xiǎn)區(qū)如港口、機(jī)場、物流中轉(zhuǎn)站、加工廠等開展早期監(jiān)測預(yù)警,以預(yù)防其再次入侵與進(jìn)一步擴(kuò)散蔓延。

    關(guān)鍵詞:北美刺龍葵; 適生區(qū); 入侵植物; 生態(tài)位模型

    隨著全球化進(jìn)程的發(fā)展,國際間頻繁的貿(mào)易、旅游等活動(dòng)極大地促進(jìn)了外來有害生物入侵(Andersonetal.,2015)。生物入侵已成為當(dāng)前世界各國關(guān)注的重要環(huán)境和生物安全問題。我國是世界上受到外來有害生物入侵危害最嚴(yán)重的國家之一,已入侵的物種不斷地?cái)U(kuò)散蔓延,新的危險(xiǎn)性有害生物不斷地傳入(萬方浩等,2009)。由于外來生物能夠自我繁殖,一旦成功入侵就很難根除,因此,風(fēng)險(xiǎn)評(píng)估和早期預(yù)警是預(yù)防外來有害生物傳入和擴(kuò)散的最有效途徑(萬方浩等,2010; Nataleetal.,2012)。

    北美刺龍葵SolanumcarolinenseL.屬茄科Solanaceae茄屬Solanum,是多年生草本植物(Anderson,1999),原產(chǎn)美國東南部伊利諾伊州、馬薩諸塞州、佛羅里達(dá)州和德克薩斯州等地區(qū)。當(dāng)前,北美刺龍葵已擴(kuò)散至美國全境以及大洋洲、歐洲、中南美洲、亞洲、非洲的36個(gè)國家和地區(qū)(Bassett & Munro,1986; Nicholsetal.,1992; Wise,2007; Wise & Sacchi,1996 )(圖1)。北美刺龍葵入侵定植后極易通過分泌化感物質(zhì)等途徑形成單優(yōu)勢群落,對(duì)當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)和生物多樣性保護(hù)構(gòu)成極大威脅,如入侵農(nóng)田后可造成農(nóng)作物減產(chǎn)35%~60%(Beeleretal.,2004; Frank,1990; Halpernetal.,2014)。北美刺龍葵還是農(nóng)作物病、蟲害的寄主,對(duì)牲畜和人類有毒(Bassett & Munro,1986)。北美刺龍葵具有極強(qiáng)的繁殖能力,每株每年可產(chǎn)生1500~7200粒種子,種子可隨風(fēng)力、水流、動(dòng)物、交通工具等途徑進(jìn)行自然和人為因素主導(dǎo)的擴(kuò)散蔓延(潘揚(yáng)等,2014; Cipollini & Levey,1997; Miyazaki,2008; Miyazaki & Ito,2004; Urakawa & Koide,2004a、2004b)。種子擴(kuò)散到新地區(qū)后,其休眠特性(休眠期可達(dá)10年)能極大地提高定植能力(Onenetal.,2006)??傊泵来听埧且环N正在全球擴(kuò)散蔓延的危險(xiǎn)性有毒植物,已被許多國家列為重點(diǎn)防控的檢疫性有害生物。

    近年來,北美刺龍葵已入侵我國并呈擴(kuò)散蔓延趨勢。北美刺龍葵最早于2006年在浙江被發(fā)現(xiàn),此后相繼在臺(tái)州市椒江區(qū)上大陳鎮(zhèn)和溫州市海島等地區(qū)發(fā)現(xiàn)(陳模舜,2008; 李根有等,2006)。同時(shí),北美刺龍葵的種子還在不斷地隨貿(mào)易等途徑傳入我國。因此,急需制定早期監(jiān)測預(yù)警措施,阻止其再次入侵,抑制已入侵種群的進(jìn)一步擴(kuò)散蔓延。

    1材料與方法

    1.1地理分布數(shù)據(jù)的搜集與整理

    北美刺龍葵在全球的分布數(shù)據(jù)主要通過GBIF(http:∥www.gbif.org/)數(shù)據(jù)庫和已發(fā)表的文獻(xiàn)獲取(Bassett & Munro,1986; Follak & Strauss,2010; Imaizumietal.,2006)。對(duì)下載自GBIF的數(shù)據(jù)進(jìn)行整理核對(duì),刪除信息不完整、重復(fù)記錄等數(shù)據(jù);對(duì)信息記錄完整但缺少地理坐標(biāo)的數(shù)據(jù)通過Google earth等手動(dòng)查詢地理坐標(biāo),同時(shí)對(duì)已經(jīng)包含地理坐標(biāo)的數(shù)據(jù)進(jìn)行核對(duì)與校正;最后把分布點(diǎn)地理坐標(biāo)存儲(chǔ)為適生區(qū)預(yù)測和時(shí)空擴(kuò)散動(dòng)態(tài)統(tǒng)計(jì)分析需要的格式。

    1.2全球擴(kuò)散的時(shí)空動(dòng)態(tài)和地理分布格局

    根據(jù)各個(gè)國家的發(fā)生時(shí)間統(tǒng)計(jì)北美刺龍葵分布點(diǎn)發(fā)生的時(shí)空動(dòng)態(tài),利用分布點(diǎn)、地統(tǒng)計(jì)方法分析其在不同空間尺度上(全球和廣泛分布國家日本)的分布模式,以判斷空間地理分布格局,最終解析其傳入和擴(kuò)散機(jī)制。

    1.3傳入我國的主要途徑分析

    通過查閱近5年來國家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局發(fā)布的進(jìn)境植物檢疫截獲有害生物報(bào)表(http:∥search.aqsiq.gov.cn/was/search)統(tǒng)計(jì)北美刺龍葵傳入我國的動(dòng)態(tài)和途徑,對(duì)極易攜帶北美刺龍葵種子的進(jìn)口貿(mào)易接受港口開展排查,確定其在我國的分布現(xiàn)狀。

    1.4氣候生態(tài)位模擬及潛在適生區(qū)的預(yù)測

    探明外來物種傳入我國的適生區(qū)對(duì)制定早期預(yù)防預(yù)警措施具有重要意義(萬方浩等,2010; Nataleetal.,2012; Wang & Wang,2006)。目前,常用的生態(tài)位模型有BIOCLIM、BIOMAPPER、DIVA、DOMAIN/、廣義相加模型(generalized additive model, GAM)、GARP、廣義線性模型(generalized linear model, GLM)、最大熵模型(Maxent)。其中,Maxent生態(tài)位模型在外來植物適生區(qū)預(yù)測中應(yīng)用最為廣泛且準(zhǔn)確性高(Elith,2014; Petersonetal.,2011)。該模型的基本原理是根據(jù)樣本信息對(duì)未知分布做出推斷,即利用物種的已知分布數(shù)據(jù)和分布地的環(huán)境數(shù)據(jù)判斷物種的生態(tài)需求,并預(yù)測物種的潛在分布,最終的模型可再被投影到地理空間,形成預(yù)測適合物種分布地區(qū)的數(shù)字地圖(Petersonetal.,2011; Phillipsetal.,2006)。本研究基于此模型,采用Maxent 3.3.3版(http:∥www.cs.princeton.edu/~schapire/maxent/)預(yù)測北美刺龍葵在中國的適生區(qū)。環(huán)境數(shù)據(jù)集來自WORLDCLIM數(shù)據(jù)集,包括海拔、氣溫、降雨等20個(gè)環(huán)境因子(http:∥www.worldclim.org),空間分辨率為5 min。

    選取合適的分布數(shù)據(jù)構(gòu)建預(yù)測模型,理論上原產(chǎn)地的分布數(shù)據(jù)是最適合構(gòu)建模型預(yù)測物種的實(shí)際生態(tài)位,盡管近年來一些研究發(fā)現(xiàn),部分外來物種在入侵地會(huì)侵占新的生態(tài)位,但生態(tài)位漂移并不是普遍現(xiàn)象(王瑞和萬方浩,2010; Petersonetal.,2011; Petitpierreetal.,2012)。因此,本研究采用原產(chǎn)地分布數(shù)據(jù)構(gòu)建生態(tài)位模型預(yù)測北美刺龍葵在全球和中國的適生區(qū)預(yù)測。

    選取與北美刺龍葵分布相關(guān)的環(huán)境數(shù)據(jù)模擬其適生的環(huán)境生態(tài)位。刀切法(jackknife)分析已被廣泛用來判斷每一個(gè)因子對(duì)預(yù)測結(jié)果的影響。其原理是在模型運(yùn)算過程中通過系統(tǒng)省略一個(gè)環(huán)境變量,使用其余的變量進(jìn)行建模預(yù)測,然后通過分析變量的有無對(duì)運(yùn)算結(jié)果的影響以判斷影響物種分布的關(guān)鍵因素(Petersonetal.,2003)。本研究采用Maxent軟件自帶的jackknife程序選擇最優(yōu)環(huán)境變量。具體運(yùn)算如下,在參數(shù)設(shè)置中勾選刀切法對(duì)話框,通過設(shè)置數(shù)據(jù)對(duì)話框,把用于模型運(yùn)算的分布數(shù)據(jù)平分為2個(gè)部分:50%為構(gòu)建模型的訓(xùn)練數(shù)據(jù)(training data),余下的50%用于模型的檢驗(yàn)(外部檢驗(yàn)數(shù)據(jù),extrinsic test data)。

    采用受試者工作特征曲線(receiver operating characteristic curve, ROC)評(píng)估上述篩選的環(huán)境參數(shù)和分布數(shù)據(jù)構(gòu)建的模型模擬北美刺龍葵環(huán)境生態(tài)位的準(zhǔn)確性與可行性。ROC曲線下的面積(area under curve, AUC)的大小作為模型預(yù)測準(zhǔn)確度的衡量指標(biāo)。AUC值位于0.5~0.6為失敗,0.6~0.7較差,0.7~0.8一般,0.8~0.9良好,0.9~1.0優(yōu)秀。AUC值越接近1,說明與隨機(jī)分布相距越遠(yuǎn),環(huán)境變量與預(yù)測物種地理分布之間的相關(guān)性越大,模型預(yù)測結(jié)果越準(zhǔn)確(王運(yùn)生等,2007)。

    把上述構(gòu)建的最優(yōu)模型投影到中國預(yù)測其適生區(qū)并生成適生性風(fēng)險(xiǎn)等級(jí)圖。在最終的預(yù)測分析中,已知的分布數(shù)據(jù)全部用于模型分析,其他參數(shù)為軟件默認(rèn)值,分析結(jié)果以Cumulative格式、ASCII類型文件輸出。將預(yù)測生成的柵格文件在ArcGIS 9.3中進(jìn)行重分類。按照模型運(yùn)算數(shù)據(jù)中已知分布點(diǎn)的最小值(minimum training presence)劃分適生區(qū)和非適生區(qū);利用幾何分級(jí)法(geometrical interval)將適生區(qū)的預(yù)測值標(biāo)記為由橙黃色到深紅色的漸變色,以直觀顯示其適生程度的變化趨勢,把上述文件輸出為圖片,做出適生區(qū)預(yù)測圖。

    2結(jié)果與分析

    2.1北美刺龍葵在全球的擴(kuò)散動(dòng)態(tài)與分布現(xiàn)狀

    北美刺龍葵最早于1901年擴(kuò)散到加拿大安大略省。此后,相繼在新西蘭(1934年)、日本(1950年)、韓國(1969年)、巴西(1974年)、英國(1975年)、克羅地亞(1978年)、芬蘭(1981年)、荷蘭(1982年)、挪威(1987)、澳大利亞(1992)、西班牙(1997年)、奧地利(2005年)、中國(2006年)、意大利(2008年)、德國(2010年)等國家發(fā)現(xiàn)該植物的分布記錄。此外,印度、尼泊爾、孟加拉也有分布報(bào)道,但入侵時(shí)間和分布地點(diǎn)不詳(圖1)。北美刺龍葵在全球的分布點(diǎn)呈指數(shù)增長趨勢(圖2)。

    北美刺龍葵在歐亞已入侵國家的最早分布點(diǎn)均位于港口等人類活動(dòng)頻繁的地區(qū),如日本70%的貿(mào)易港口都有其分布記錄,橫濱港、神戶港和岡山港3個(gè)主要貿(mào)易港口是主要聚集區(qū)(圖3)。

    2.2北美刺龍葵在我國的分布、傳入頻率與途徑

    目前,我國的北美刺龍葵主要分布在浙江省臺(tái)州市和溫州市境內(nèi)。根據(jù)進(jìn)境植物檢疫截獲有害生物月度通報(bào)數(shù)據(jù)(2013年11月—2015年8月)顯示,近2年的截獲次數(shù)呈指數(shù)增加,截至2015年8月累計(jì)截獲115次(圖4),其中99.6%源于進(jìn)口糧谷中攜帶的該植物種子。河北秦皇島、廣西防城港、山東濱州、日照嵐山、浙江寧波口岸均從美國進(jìn)口的大豆、高粱、小麥中截獲了該植物的種子。

    2.3北美刺龍葵在我國的適生區(qū)

    北美刺龍葵生態(tài)位模型環(huán)境參數(shù)刀切法分析表明,年平均溫度(bio1)、晝夜溫差與年溫差比值(bio3)、最濕季平均溫度(bio8)、最冷季平均溫度(bio11)、最濕月份降雨量(bio13)、最干月份降雨量 (bio14)、降雨量的季節(jié)變化(bio15)對(duì)模型的貢獻(xiàn)率最高,是影響北美刺龍葵分布的關(guān)鍵因素,因而被用來建立最終模型預(yù)測北美刺龍葵的適生區(qū)。利用上述篩選的環(huán)境變量和原產(chǎn)地分布數(shù)據(jù)構(gòu)建模型,將模型模擬的環(huán)境生態(tài)位投影到全球,AUC值為0.789,遠(yuǎn)高于隨機(jī)預(yù)測值0.5,表明本研究建立的Maxent模型的預(yù)測能力較強(qiáng),能夠很好地?cái)M合物種已知分布的環(huán)境生態(tài)位。

    把該模型投影到全球與中國,并按照適生性等級(jí)的劃分方法對(duì)預(yù)測值進(jìn)行分級(jí)以生成適生性等級(jí)圖,把區(qū)劃后的適生性等級(jí)圖與行政區(qū)劃疊加,生成北美刺龍葵在全球和中國的適生性風(fēng)險(xiǎn)預(yù)測圖(圖5、6)。全球預(yù)測結(jié)果顯示,中南美洲、歐洲中南部、東南亞、非洲等地都是其適生區(qū),且其在全球的已知分布點(diǎn)均被預(yù)測為適生區(qū)(圖5)。這也間接顯示北美刺龍葵原產(chǎn)地?cái)?shù)據(jù)模擬的氣候生態(tài)位能夠用來預(yù)測其適生區(qū)。

    北美刺龍葵在中國的適生區(qū)預(yù)測圖顯示,除黑龍江、吉林、內(nèi)蒙古、青海、甘肅、西藏、四川西北部以外的區(qū)域都是其在我國的適生區(qū),其中高風(fēng)險(xiǎn)區(qū)主要集中在東部和南部沿海、西南邊境和新疆的部分地區(qū)(圖6)。

    3討論

    3.1北美刺龍葵在全球擴(kuò)散的驅(qū)動(dòng)力與傳入我國的主要途徑推斷

    北美刺龍葵在全球的擴(kuò)散動(dòng)態(tài)表明其在原產(chǎn)地以外地區(qū)的最早發(fā)生點(diǎn)主要位于貿(mào)易等人類活動(dòng)頻繁的港口、交通樞紐等區(qū)域。北美刺龍葵是一種有毒植物,目前并沒有發(fā)現(xiàn)其在全球范圍內(nèi)的引種記錄?;诒泵来听埧谌虻姆植几窬旨皶r(shí)空動(dòng)態(tài),有理由推斷跨區(qū)域的貿(mào)易活動(dòng)中的無意攜帶是其在全球擴(kuò)散蔓延的主要途徑與驅(qū)動(dòng)力。此外,北美刺龍葵在日本的地理分布格局表明,多次傳入是其成功入侵與廣泛分布危害的主要根源。Imaizumietal.(2006)的研究表明,北美刺龍葵在日本的東京地區(qū)主要沿高野河擴(kuò)散蔓延。此外,已經(jīng)定植的種群也在隨人類活動(dòng)沿道路進(jìn)行長距離的跳躍式擴(kuò)散。北美刺龍葵具有極強(qiáng)的繁殖能力,且具有種子和營養(yǎng)體雙重繁殖方式等特性,均有助于種群的定植、增長和擴(kuò)散蔓延。因此,國際貿(mào)易活動(dòng)中攜帶的種子是導(dǎo)致其全球擴(kuò)張的主要驅(qū)動(dòng)力,定植后跨區(qū)域的人類活動(dòng)和隨不同媒介擴(kuò)散的潛力會(huì)促進(jìn)其進(jìn)一步的擴(kuò)散蔓延。

    我國已于2007年把北美刺龍葵列為進(jìn)境檢疫性有害生物和重點(diǎn)監(jiān)測的對(duì)象(吳海榮等,2008)。盡管如此,北美刺龍葵還是成功入侵我國,目前主要分布在我國東部沿海地區(qū),應(yīng)該是隨著人類活動(dòng)無意攜帶傳入我國。隨著我國一帶一路戰(zhàn)略的發(fā)展和自由貿(mào)易區(qū)的建設(shè),對(duì)外交流的區(qū)域、貿(mào)易量和頻次勢必會(huì)快速增長,這也為北美刺龍葵等有害生物進(jìn)入我國更多區(qū)域提供了更多機(jī)會(huì),如近年來北美刺龍葵的截獲次數(shù)就呈指數(shù)增長趨勢。目前,我國對(duì)外開放的一類口岸有253個(gè),其中水運(yùn)口岸133個(gè)、鐵路口岸17個(gè)、公路口岸47個(gè)、航空口岸56個(gè),廣泛分布于我國的東部和南部沿海、西部和北部邊境線和長江流域??诎兜倪@種離散式的地理分布格局也可為北美刺龍葵等外來有害生物通過貿(mào)易途徑同時(shí)、多次傳入我國不同地區(qū)提供了機(jī)會(huì),這些通過貿(mào)易攜帶的有害生物還會(huì)隨著本地運(yùn)輸而擴(kuò)散到更多地區(qū)。同時(shí),已入侵的種群也會(huì)不斷地繁殖、增長和擴(kuò)散蔓延。總之,北美刺龍葵在我國傳入和擴(kuò)散的風(fēng)險(xiǎn)極高,急需發(fā)展早期監(jiān)測預(yù)警措施預(yù)防其再次傳入,抑制其進(jìn)一步擴(kuò)散。

    3.2加強(qiáng)北美刺龍葵的早期監(jiān)測預(yù)警

    探明外來物種的傳入和擴(kuò)散途徑,根據(jù)傳入媒介的空間分布,結(jié)合影響種群定植和繁殖生長的環(huán)境條件,劃定物種重點(diǎn)監(jiān)測的區(qū)域,通過系統(tǒng)監(jiān)測達(dá)到早發(fā)現(xiàn)、早根除的目標(biāo)是控制有害生物入侵和擴(kuò)散蔓延的最有效途徑(Andersonetal.,2015; Nataleetal.,2012; Wangetal.,2011)。

    基于北美刺龍葵在我國的傳入與擴(kuò)散途徑,結(jié)合其在我國的預(yù)測適生區(qū),本研究繪制了北美刺龍葵在中國的早期監(jiān)測預(yù)警圖,以便開展系統(tǒng)監(jiān)測與制定針對(duì)性防控措施(圖6)。首先,抑制已入侵種群的進(jìn)一步擴(kuò)散蔓延。由于北美刺龍葵具有極強(qiáng)的繁殖能力,其種子可隨水流和人類活動(dòng)擴(kuò)散,因此應(yīng)該加緊對(duì)已入侵種群10 km范圍內(nèi)的河流、道路的監(jiān)測,同時(shí)采取人工拔除、化學(xué)防治等手段控制已發(fā)生地區(qū)種群的增長,根除新發(fā)現(xiàn)的植株。其次,抑制繁殖體的再次入侵?;诒泵来听埧饕S進(jìn)口糧谷攜帶傳入,根據(jù)進(jìn)口糧谷在我國的主要運(yùn)轉(zhuǎn)途徑,確定重點(diǎn)監(jiān)測港口、中轉(zhuǎn)站和加工/存儲(chǔ)地,結(jié)合北美刺龍葵在我國的環(huán)境適應(yīng)性決定重點(diǎn)監(jiān)測區(qū)域。建議對(duì)這些地區(qū)開展系統(tǒng)監(jiān)測,以預(yù)防其進(jìn)一步入侵。

    參考文獻(xiàn)

    陳模舜, 2008. 臺(tái)州市外來植物區(qū)系組成與分布特征的研究. 臺(tái)州學(xué)院學(xué)報(bào), 30(6): 43-48.

    李根有, 金水虎, 哀建國, 2006. 浙江省有害植物種類、特點(diǎn)及防治. 浙江林學(xué)院學(xué)報(bào), 23(6): 614-624.

    潘揚(yáng), 羅芳, 魯長虎, 2014. 脊椎動(dòng)物傳播植物肉質(zhì)果中的次生物質(zhì)及其生態(tài)作用. 生態(tài)學(xué)報(bào), 34(10): 2490-2497.

    萬方浩, 郭建英, 張峰, 2009. 中國生物入侵研究. 北京: 科學(xué)出版社.

    萬方浩, 彭德良, 王瑞, 2010. 生物入侵: 預(yù)警篇. 北京: 科學(xué)出版社.

    王瑞, 萬方浩, 2010. 外來入侵植物意大利蒼耳在我國適生區(qū)預(yù)測. 草業(yè)學(xué)報(bào), 19(6): 222-230.

    王運(yùn)生, 謝丙炎, 萬方浩, 肖啟明, 戴良英, 2007. ROC曲線分析在評(píng)價(jià)入侵物種分布模型中的應(yīng)用. 生物多樣性, 15(4): 365-372.

    吳海榮, 鐘國強(qiáng), 胡學(xué)難, 王定國, 2008. 淺析我國新頒布進(jìn)境檢疫雜草名錄的特點(diǎn). 植物檢疫, 22(4): 231-233.

    Anderson L G, Rocliffe S, Haddaway N R and Dunn A M, 2015. The role of tourism and recreation in the spread of non-native species: a systematic review and meta-analysis.PLoSONE, 10(10): e0140833.

    Anderson W P, 1999.PerennialWeeds:CharacteristicsandIdentificationofSelectedHerbaceousSpecies. Ames, IA, USA: Iowa State University Press.

    Bassett I J and Munro D B, 1986. The biology of Canadian weeds. 78.SolanumcarolinenseL. andSolanumrostratumDunal.CanadianJournalofPlantScience, 66(4): 977-991.

    Beeler J E, Rhodes G N, Bates G E, Main C L and Mueller T C, 2004. Horsenettle (Solanumcarolinense) control in tall fescue (Festucaarundinacea) and clover (Trifoliumsp.) pastures with mixtures of 2,4-D and picloram.WeedTechnology, 18(4): 1091-1095.

    Cipollini M L and Levey D J, 1997. Antifungal activity ofSolanumfruit glycoalkaloids: implications for frugivory and seed dispersal.Ecology, 78(3): 799-809.

    Elith J, 2014. Predicting distributions of invasive species (Chapter 6)∥Walshe T R, Robinson A, Nunn M and Burgman M A.InvasiveSpecies:RiskAssessmentandManagement. Cambridge: Cambridge University Press: 1-28.

    Follak S and Strauss G, 2010. Potential distribution and management of the invasive weedSolanumcarolinensein Central Europe.WeedResearch, 50(6): 544-552.

    Frank J R, 1990. Influence of horsenettle (Solanumcarolinense) on snapbean (Phaseolusvulgaris).WeedScience, 38(3): 220-223.

    Halpern S L, Bednar D, Chisholm A and Underwood N, 2014. Plant-mediated effects of host plant density on a specialist herbivore ofSolanumcarolinense.EcologicalEntomology, 39(2): 217-225.

    Imaizumi T, Kurokawa S, Ito M, Auld B and Wang G X, 2006. Population structure ofSolanumcarolinensealong the Takano River in Kyoto, Japan as determined by amplified fragment length polymorphism analysis.WeedResearch, 46(3): 219-225.

    Miyazaki K, 2008. Root system architecture and its relationship to the vegetative reproduction function in horsenettle (Solanumcarolinense).WeedBiologyandManagement, 8(2): 97-103.

    Miyazaki K and Ito M, 2004. Root system structure and shoot arrangement of 1-year-oldSolanumcarolinenseL.WeedBiologyandManagement, 4(2): 122-125.

    Natale E, Zalba M S, Reinoso H and Damilano G, 2012. Assessing invasion process through pathway and vector analysis: case of saltcedar (Tamarixspp.).ManagementofBiologicalInvasions, 3(1): 37-44.

    Nichols R L, Cardina J, Lynch R L, Minton N A and Wells H D, 1992. Insects, nematodes, and pathogens associated with horsenettle (Solanumcarolinense) in bermudagrass (Cynodondactylon) pastures.WeedScience, 40(2): 320-325.

    Onen H, Misako I and Imaizumi T, 2006. Horsenettle (SolanumcarolinenseL.) plants emerged at different times after corn (ZeamaysL.) planting.WeedBiologyandManagement, 6(1): 55-58.

    Peterson A T, Papes M and Kluza D A, 2003. Predicting the potential invasive distribution of four alien plant species in north America.WeedScience, 51(6): 863-868.

    Peterson A T, Soberon J, Pearson R G, Anderson R P, Martinez-Meyer E, Nakamura M and Araujo M B, 2011.EcologicalNichesandGeographicDistributions. Princeton, NJ: Princeton University Press.

    Petitpierre B, Kueffer C, Broennimann O, Randin C, Daehler C and Guisan A, 2012. Climatic niche shifts are rare among terrestrial plant invaders.Science, 335: 1344-1348.

    Phillips S J, Anderson R P and Schapire R E, 2006. Maximum entropy modeling of species geographic distributions.EcologicalModelling, 190(3-4): 231-259.

    Urakawa S and Koide I, 2004a. Horsenettle (SolanumcarolinenseL.) diffusion through rotary tilling in invaded corn fields.GrasslandScience, 50(2): 194-200.

    Urakawa S and Koide I, 2004b. Growth characteristics and establishment of seeds of horsenettle (SolanumcarolinenseL.) in the cornfield.GrasslandScience, 50(1): 64-70.

    Wang R, Wang J F, Qiu Z J, Meng B, Wan F H and Wang Y Z, 2011. Multiple mechanisms underlie rapid expansion of an invasive alien plant.NewPhytologist, 191(3): 828-839.

    Wang R and Wang Y Z, 2006. Invasion dynamics and potential spread of the invasive alien plant speciesAgeratinaadenophora(Asteraceae) in China.DiversityandDistributions, 12(4): 397-408.

    Wise M J, 2007. The herbivores ofSolanumcarolinense(Horsenettle) in northern Virginia: natural history and damage assessment.SoutheasternNaturalist, 6(3): 505-522.

    Wise M J and Sacchi C F, 1996. Impact of two specialist insect herbivores on reproduction of horse nettle,Solanumcarolinense.Oecologia, 108(2): 328-337.

    (責(zé)任編輯:郭瑩)

    Predicting the potential invasive distribution ofSolanumcarolinensein China

    Rui WANG1, Xiao-qing XIAN1, Fang-hao WAN1,2*

    1StateKeyLaboratoryforBiologyandPlantDiseasesandInsectPests,InstituteofPlantProtection,ChineseAcademyofAgriculturalSciences,Beijing100193,China;2SchoolofAgronomyandPlantProtection,QingdaoAgriculturalUniversity,Qingdao,Shandong266109,China

    Abstract:【Background】Solanumcarolinenseis a notorious weed in America where it is native and other countries in Europe and Asia where it is introduced. In recent decade,S.carloninensehas been declared a harmful quarantine weed in China since it invaded southern China as its seeds were frequently intercepted from imported commodities. It is necessary to identify areas at risk of an invasion and develop an efficient management strategy. 【Method】 GIS and spatial statistics were used to reconstruct worldwide invasion dynamics and identify its introduction vectors. MaxEnt niche modeling was applied to predict the potential distribution ofS.carolinensein China and the model performance was evaluated by Receiver Operator Characteristic (ROC) curve. 【Result】 Its historical introduction from native America to other continents showed that the trade of commodities containing seeds promoted its invasion and spread across the world. The area under the ROC curve (AUC) was 0.789, which indicated that climatic niche models on the basis of native data had a highly predictive ability. And thus the modeled climatic niches were projected on to the landscape of the China. With the exception of Heilongjiang, Jilin, Inner Mongolia, Qinghai, Gansu, Tibet, northwestern parts of Sichuan, most areas in China were potential areas forS.carolinenseinvasion. 【Conclusion and significance】 High introduction risk with increasing trade coupled with large potentially available distribution area may lead to increased invasion in the future. An early monitoring and eradication program should be adopted in the highly vulnerable areas, including coastal harbors, airports, transfer stations and processing factories of the imported cereal grains or fodder crops, to prevent further invasion and spread.

    Key words:Solanumcarolinense; potential invasive distribution; invasive plant; ecological niche model

    收稿日期(Received):2015-05-09接受日期(Accepted): 2015-06-21

    基金項(xiàng)目:國家自然科學(xué)基金(31471827、31000883); 環(huán)保公益性行業(yè)專項(xiàng)(201409061); 國家科技支撐計(jì)劃課題(2015BAD08B03); 山東省“泰山學(xué)者”建設(shè)工程專項(xiàng)

    作者簡介:王瑞, 男, 副研究員。 研究方向: 入侵物種的時(shí)空擴(kuò)散機(jī)制與早期監(jiān)測預(yù)警。 E-mail: wangrcaas@163.com *通訊作者(Author of correspondence), E-mail: wanfanghao@caas.cn

    DOI:10. 3969/j.issn.2095-1787.2016.02.005

    猜你喜歡
    適生區(qū)
    氣候變化背景下藥用植物牡丹的潛在適生區(qū)分析
    氣候變化對(duì)落葉松人工林在中國適生區(qū)分布的影響
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預(yù)測
    未來氣候條件下當(dāng)歸適生區(qū)預(yù)測及時(shí)空變化分析
    當(dāng)前及未來氣候下四川省草地貪夜蛾適生區(qū)分布
    基于物種分布模型的畫稿溪國家級(jí)自然保護(hù)區(qū)桫欏保護(hù)現(xiàn)狀及影響因素*
    基于MaxEnt和ArcGIS的烏檀生長適應(yīng)性區(qū)劃研究
    基于組合模型預(yù)測蛇足石杉在中國的適生區(qū)
    不同氣候條件下桫欏在中國的潛在適生區(qū)分布
    氣候變化下瀕危植物半日花在中國的潛在分布
    啦啦啦韩国在线观看视频| 欧美激情国产日韩精品一区| 九九久久精品国产亚洲av麻豆| 亚洲午夜理论影院| 非洲黑人性xxxx精品又粗又长| 一个人观看的视频www高清免费观看| 99久久精品热视频| 亚洲最大成人中文| 日本 欧美在线| 级片在线观看| 国产乱人视频| 免费搜索国产男女视频| 国产亚洲精品av在线| 变态另类丝袜制服| 国产极品精品免费视频能看的| 极品教师在线视频| 国产真实伦视频高清在线观看 | 99热精品在线国产| 久久久色成人| 综合色av麻豆| 成人毛片a级毛片在线播放| 久久久国产成人精品二区| 热99在线观看视频| 亚洲性久久影院| 最近最新免费中文字幕在线| 免费电影在线观看免费观看| 亚洲国产高清在线一区二区三| 日日摸夜夜添夜夜添小说| 性欧美人与动物交配| 男女之事视频高清在线观看| 成人av一区二区三区在线看| 在线国产一区二区在线| 午夜福利在线观看吧| 国产 一区精品| 三级国产精品欧美在线观看| 搞女人的毛片| 国产亚洲欧美98| 午夜免费激情av| 天堂动漫精品| 身体一侧抽搐| 国产在视频线在精品| 午夜精品久久久久久毛片777| 国产精品三级大全| 久久精品久久久久久噜噜老黄 | 日韩强制内射视频| 此物有八面人人有两片| 国产探花极品一区二区| 九九爱精品视频在线观看| 精品欧美国产一区二区三| 国产精品自产拍在线观看55亚洲| 亚洲精品日韩av片在线观看| 精品一区二区三区av网在线观看| 99国产精品一区二区蜜桃av| 人人妻人人澡欧美一区二区| 欧美日韩综合久久久久久 | 国产成人影院久久av| 18禁在线播放成人免费| 亚洲欧美日韩卡通动漫| 99精品久久久久人妻精品| 91精品国产九色| 级片在线观看| 日韩在线高清观看一区二区三区 | 欧洲精品卡2卡3卡4卡5卡区| 毛片一级片免费看久久久久 | 亚洲av第一区精品v没综合| 色噜噜av男人的天堂激情| 亚洲成人久久性| 国产白丝娇喘喷水9色精品| 午夜精品一区二区三区免费看| 亚洲人成网站在线播| 性欧美人与动物交配| 欧美成人a在线观看| 美女黄网站色视频| 身体一侧抽搐| 国产精品乱码一区二三区的特点| 可以在线观看的亚洲视频| 日韩强制内射视频| 国产高清视频在线观看网站| 三级男女做爰猛烈吃奶摸视频| 国产av一区在线观看免费| 亚洲一级一片aⅴ在线观看| 亚洲国产精品合色在线| 高清在线国产一区| 欧美极品一区二区三区四区| 亚洲avbb在线观看| 国产精品永久免费网站| bbb黄色大片| 国产麻豆成人av免费视频| 午夜爱爱视频在线播放| 女的被弄到高潮叫床怎么办 | av天堂在线播放| av在线蜜桃| 午夜福利高清视频| 亚洲精品在线观看二区| av天堂在线播放| 国产精品无大码| 国产精品福利在线免费观看| 村上凉子中文字幕在线| 欧美日韩精品成人综合77777| 久久精品夜夜夜夜夜久久蜜豆| 国产精品av视频在线免费观看| 国产一区二区在线观看日韩| 九九久久精品国产亚洲av麻豆| 亚洲自偷自拍三级| 日韩欧美三级三区| АⅤ资源中文在线天堂| 久久精品国产清高在天天线| 在线免费观看的www视频| 精品午夜福利在线看| 女生性感内裤真人,穿戴方法视频| 制服丝袜大香蕉在线| 色播亚洲综合网| 久久久久国内视频| 精品人妻熟女av久视频| 国产精品久久电影中文字幕| 国产av一区在线观看免费| 免费一级毛片在线播放高清视频| 色综合色国产| av.在线天堂| 色哟哟哟哟哟哟| 97碰自拍视频| 无遮挡黄片免费观看| 国产精品久久久久久精品电影| 亚洲狠狠婷婷综合久久图片| 色播亚洲综合网| 久久亚洲真实| 亚洲最大成人手机在线| 日本欧美国产在线视频| 亚洲欧美激情综合另类| 美女 人体艺术 gogo| 久久久久性生活片| 韩国av在线不卡| 日本黄色片子视频| 久久草成人影院| 俺也久久电影网| 日本黄色片子视频| 蜜桃亚洲精品一区二区三区| 精品国内亚洲2022精品成人| 成人精品一区二区免费| 国产成人一区二区在线| 国产高清不卡午夜福利| 久久人人精品亚洲av| 久久精品国产亚洲av香蕉五月| 久久99热这里只有精品18| 无遮挡黄片免费观看| 久久精品夜夜夜夜夜久久蜜豆| 日本与韩国留学比较| 色在线成人网| 国产私拍福利视频在线观看| АⅤ资源中文在线天堂| 变态另类丝袜制服| 香蕉av资源在线| 神马国产精品三级电影在线观看| 亚洲,欧美,日韩| 国产探花在线观看一区二区| 日韩欧美在线二视频| 亚洲精品456在线播放app | 夜夜爽天天搞| 我的老师免费观看完整版| 99久久久亚洲精品蜜臀av| 成人av一区二区三区在线看| 午夜福利18| 亚洲国产精品sss在线观看| 搡老岳熟女国产| 在线观看一区二区三区| 亚洲一级一片aⅴ在线观看| 欧美一区二区国产精品久久精品| 天堂√8在线中文| 少妇裸体淫交视频免费看高清| 成人特级av手机在线观看| 久久久久久伊人网av| 国产精品av视频在线免费观看| 日韩中字成人| 男人舔女人下体高潮全视频| 精品人妻视频免费看| 成人国产一区最新在线观看| 国产爱豆传媒在线观看| 免费看a级黄色片| 国产爱豆传媒在线观看| 91久久精品电影网| 99久久久亚洲精品蜜臀av| 日本-黄色视频高清免费观看| 国产精品98久久久久久宅男小说| 网址你懂的国产日韩在线| 特级一级黄色大片| 看黄色毛片网站| 亚洲在线观看片| 国产在线精品亚洲第一网站| 中国美女看黄片| 亚洲最大成人手机在线| 在现免费观看毛片| 最近最新免费中文字幕在线| 中亚洲国语对白在线视频| 丰满乱子伦码专区| 琪琪午夜伦伦电影理论片6080| 国产av不卡久久| 丰满的人妻完整版| 亚洲18禁久久av| 亚洲无线在线观看| 欧美性感艳星| 中国美女看黄片| 99久久精品国产国产毛片| 成人无遮挡网站| 好男人在线观看高清免费视频| 国产视频内射| 蜜桃久久精品国产亚洲av| 人妻丰满熟妇av一区二区三区| 欧美人与善性xxx| 成年女人看的毛片在线观看| 亚洲成人精品中文字幕电影| 一区福利在线观看| 九九热线精品视视频播放| 少妇猛男粗大的猛烈进出视频 | 亚洲av.av天堂| 五月玫瑰六月丁香| 有码 亚洲区| 日日干狠狠操夜夜爽| 国产精华一区二区三区| 天堂av国产一区二区熟女人妻| 国产精品一区二区性色av| 久久99热这里只有精品18| 亚洲国产欧美人成| 国产精品爽爽va在线观看网站| 99九九线精品视频在线观看视频| 超碰av人人做人人爽久久| 亚洲av熟女| 可以在线观看的亚洲视频| 日本黄色片子视频| 亚洲性久久影院| 日韩欧美在线二视频| 国产精品永久免费网站| 久久久精品大字幕| 国产在线男女| 色哟哟哟哟哟哟| 91久久精品电影网| 3wmmmm亚洲av在线观看| 色综合亚洲欧美另类图片| 综合色av麻豆| 97热精品久久久久久| 成人av在线播放网站| 老师上课跳d突然被开到最大视频| 又紧又爽又黄一区二区| 色在线成人网| 国产精品国产三级国产av玫瑰| 国产精品人妻久久久影院| 国产精品无大码| 69人妻影院| 夜夜夜夜夜久久久久| 亚洲精华国产精华液的使用体验 | 久久天躁狠狠躁夜夜2o2o| 春色校园在线视频观看| 一级黄片播放器| 女的被弄到高潮叫床怎么办 | 美女大奶头视频| 欧美黑人欧美精品刺激| 成人av在线播放网站| 校园春色视频在线观看| 内射极品少妇av片p| 国产女主播在线喷水免费视频网站 | 亚洲av中文字字幕乱码综合| 内地一区二区视频在线| 窝窝影院91人妻| 日本成人三级电影网站| 国产一级毛片七仙女欲春2| 悠悠久久av| 波多野结衣高清无吗| 久久精品人妻少妇| 无遮挡黄片免费观看| 99久久精品热视频| 88av欧美| 在线观看美女被高潮喷水网站| 国产亚洲精品久久久com| 成年免费大片在线观看| 精品久久久久久成人av| 国产aⅴ精品一区二区三区波| 国产精品野战在线观看| 亚洲欧美日韩东京热| 一个人看视频在线观看www免费| 国产午夜精品论理片| 99视频精品全部免费 在线| 男女啪啪激烈高潮av片| 精品人妻偷拍中文字幕| 久久久成人免费电影| 免费看美女性在线毛片视频| 国产成人影院久久av| 男人和女人高潮做爰伦理| 一级av片app| 麻豆成人午夜福利视频| 国产v大片淫在线免费观看| 校园春色视频在线观看| 色综合婷婷激情| 黄色日韩在线| 一区二区三区高清视频在线| 午夜福利18| 亚洲精品粉嫩美女一区| 日本一二三区视频观看| 免费人成视频x8x8入口观看| 女人被狂操c到高潮| 真实男女啪啪啪动态图| 精品久久久久久久久久免费视频| 韩国av一区二区三区四区| 一个人看的www免费观看视频| 22中文网久久字幕| 91久久精品电影网| 国产一区二区三区av在线 | 久久6这里有精品| 日本a在线网址| 成人无遮挡网站| 男人和女人高潮做爰伦理| 成人二区视频| 精品久久久久久成人av| 国产精品电影一区二区三区| 在线观看66精品国产| 免费无遮挡裸体视频| 亚洲精品一卡2卡三卡4卡5卡| 国产熟女欧美一区二区| 国产主播在线观看一区二区| 国产亚洲av嫩草精品影院| 99热这里只有是精品在线观看| 99久久成人亚洲精品观看| 国产高清三级在线| 久久久久国产精品人妻aⅴ院| 人妻夜夜爽99麻豆av| av天堂在线播放| 波多野结衣高清作品| 丰满乱子伦码专区| 久久国内精品自在自线图片| 校园春色视频在线观看| 成人永久免费在线观看视频| 午夜福利在线在线| 国产精品人妻久久久影院| 午夜福利视频1000在线观看| 日本黄色片子视频| 久久久久久久久久久丰满 | 欧美日韩黄片免| 日韩欧美 国产精品| 大型黄色视频在线免费观看| 亚洲一级一片aⅴ在线观看| 久久久久久大精品| 成人av在线播放网站| 久久人妻av系列| 国产精品不卡视频一区二区| 亚洲中文字幕日韩| 精品久久久久久久久亚洲 | 热99re8久久精品国产| 亚洲av不卡在线观看| 日韩精品中文字幕看吧| 日韩精品青青久久久久久| 深夜a级毛片| 精品久久国产蜜桃| 午夜老司机福利剧场| 麻豆国产av国片精品| 九九爱精品视频在线观看| 国产伦精品一区二区三区四那| 全区人妻精品视频| 一区福利在线观看| 最后的刺客免费高清国语| 午夜激情福利司机影院| 国产av麻豆久久久久久久| 熟妇人妻久久中文字幕3abv| 国产精品一区二区三区四区免费观看 | 嫩草影院精品99| 欧美一区二区亚洲| 亚洲欧美日韩卡通动漫| 女人被狂操c到高潮| 国产亚洲精品av在线| 免费观看在线日韩| 精品久久久久久,| 国产男人的电影天堂91| 欧美日韩乱码在线| 亚洲欧美日韩卡通动漫| 美女cb高潮喷水在线观看| 热99re8久久精品国产| 天堂网av新在线| 在线观看免费视频日本深夜| 男女下面进入的视频免费午夜| 我要看日韩黄色一级片| 亚洲av免费高清在线观看| 中文在线观看免费www的网站| 在线a可以看的网站| 日韩一本色道免费dvd| 狂野欧美激情性xxxx在线观看| 99在线视频只有这里精品首页| 亚洲成a人片在线一区二区| 麻豆一二三区av精品| 国产黄片美女视频| 少妇丰满av| 一个人观看的视频www高清免费观看| 性色avwww在线观看| 人妻久久中文字幕网| 亚洲,欧美,日韩| 国产主播在线观看一区二区| 欧美成人a在线观看| 色尼玛亚洲综合影院| av天堂中文字幕网| 国产欧美日韩精品亚洲av| 成人性生交大片免费视频hd| 联通29元200g的流量卡| 国产伦在线观看视频一区| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产自在天天线| 成人无遮挡网站| 中文字幕精品亚洲无线码一区| 国产成人a区在线观看| 亚洲人成网站在线播放欧美日韩| 国产成年人精品一区二区| 1024手机看黄色片| 日本一本二区三区精品| 午夜福利在线观看免费完整高清在 | 久久精品人妻少妇| 日本精品一区二区三区蜜桃| videossex国产| 国产aⅴ精品一区二区三区波| 男人舔奶头视频| 久久久久久久精品吃奶| 欧美一级a爱片免费观看看| 国产精品无大码| 亚洲五月天丁香| 小蜜桃在线观看免费完整版高清| 午夜影院日韩av| 免费搜索国产男女视频| 日韩欧美在线二视频| 国产极品精品免费视频能看的| 国产欧美日韩精品亚洲av| 两人在一起打扑克的视频| 欧美不卡视频在线免费观看| 88av欧美| 无人区码免费观看不卡| 亚洲,欧美,日韩| 国产高清视频在线观看网站| 成人无遮挡网站| 伦理电影大哥的女人| 精品免费久久久久久久清纯| 日日摸夜夜添夜夜添小说| 草草在线视频免费看| 人妻少妇偷人精品九色| 在线观看66精品国产| 亚洲精品粉嫩美女一区| 老司机福利观看| 在线观看免费视频日本深夜| 国产男人的电影天堂91| 啦啦啦啦在线视频资源| 熟女电影av网| 亚洲av日韩精品久久久久久密| 中文字幕免费在线视频6| 特级一级黄色大片| 国产免费av片在线观看野外av| 国产亚洲精品久久久com| 亚洲va日本ⅴa欧美va伊人久久| 国产精品一及| 国产视频内射| 欧美最新免费一区二区三区| 女生性感内裤真人,穿戴方法视频| 成人二区视频| 亚洲电影在线观看av| 国产精品女同一区二区软件 | 波多野结衣巨乳人妻| avwww免费| 一级黄片播放器| 欧美日韩瑟瑟在线播放| 国产男人的电影天堂91| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品综合一区在线观看| 久久婷婷人人爽人人干人人爱| 国产亚洲精品久久久com| 欧美人与善性xxx| 午夜福利在线在线| 精品人妻偷拍中文字幕| 亚洲七黄色美女视频| 窝窝影院91人妻| 色噜噜av男人的天堂激情| 桃色一区二区三区在线观看| 12—13女人毛片做爰片一| 欧美成人a在线观看| 亚洲av中文字字幕乱码综合| 男人舔女人下体高潮全视频| 能在线免费观看的黄片| 国产精品98久久久久久宅男小说| 日韩强制内射视频| 免费观看精品视频网站| 3wmmmm亚洲av在线观看| 久久久色成人| 熟女电影av网| 亚洲五月天丁香| 欧美+日韩+精品| 国产精品人妻久久久影院| 女人被狂操c到高潮| 身体一侧抽搐| 观看美女的网站| 午夜免费男女啪啪视频观看 | 一区二区三区高清视频在线| 亚洲av免费高清在线观看| 999久久久精品免费观看国产| 国产亚洲av嫩草精品影院| 成人国产综合亚洲| 校园人妻丝袜中文字幕| 色5月婷婷丁香| 长腿黑丝高跟| 国产av一区在线观看免费| 村上凉子中文字幕在线| 在现免费观看毛片| 婷婷精品国产亚洲av| 狠狠狠狠99中文字幕| 久久久久久久久中文| 88av欧美| 午夜久久久久精精品| 草草在线视频免费看| 成熟少妇高潮喷水视频| 国产精品福利在线免费观看| 精品欧美国产一区二区三| 国产白丝娇喘喷水9色精品| 观看美女的网站| 国产午夜精品论理片| 久久久久久国产a免费观看| 九色国产91popny在线| 国产私拍福利视频在线观看| 中文字幕免费在线视频6| 少妇人妻一区二区三区视频| 国产成人av教育| 国产免费一级a男人的天堂| 在线观看午夜福利视频| 伊人久久精品亚洲午夜| 91久久精品国产一区二区三区| 联通29元200g的流量卡| 深夜精品福利| 亚洲久久久久久中文字幕| 美女免费视频网站| 麻豆精品久久久久久蜜桃| 人妻夜夜爽99麻豆av| 免费搜索国产男女视频| 亚洲最大成人中文| 制服丝袜大香蕉在线| 在线国产一区二区在线| 国产精品三级大全| 欧美最黄视频在线播放免费| 久久久久久久久久黄片| 九色国产91popny在线| 美女黄网站色视频| 国产色爽女视频免费观看| 丰满的人妻完整版| 久久精品国产亚洲av天美| 成人三级黄色视频| 国产精品美女特级片免费视频播放器| 亚洲人成网站在线播放欧美日韩| 夜夜爽天天搞| 亚洲欧美日韩东京热| 99九九线精品视频在线观看视频| 亚洲国产日韩欧美精品在线观看| 91久久精品国产一区二区三区| 真人做人爱边吃奶动态| 亚洲熟妇熟女久久| 国产色婷婷99| 长腿黑丝高跟| 午夜精品一区二区三区免费看| 免费不卡的大黄色大毛片视频在线观看 | 成人欧美大片| 波多野结衣高清作品| 国产高清三级在线| 久久久久久久久大av| 一卡2卡三卡四卡精品乱码亚洲| 精华霜和精华液先用哪个| 欧美成人免费av一区二区三区| 国产一区二区在线观看日韩| 91精品国产九色| 午夜福利欧美成人| 在线播放国产精品三级| 国产精品,欧美在线| 国产高清三级在线| 免费高清视频大片| 国产单亲对白刺激| eeuss影院久久| 国产单亲对白刺激| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 精华霜和精华液先用哪个| 一卡2卡三卡四卡精品乱码亚洲| 一区福利在线观看| 精品人妻一区二区三区麻豆 | 国产精品1区2区在线观看.| 国产精品久久久久久av不卡| 在线a可以看的网站| a在线观看视频网站| 国产极品精品免费视频能看的| 日韩亚洲欧美综合| 国产乱人视频| 午夜激情欧美在线| 自拍偷自拍亚洲精品老妇| 免费看av在线观看网站| av国产免费在线观看| 久久久久精品国产欧美久久久| 国产成年人精品一区二区| 免费观看精品视频网站| 中文字幕熟女人妻在线| 熟妇人妻久久中文字幕3abv| 69av精品久久久久久| 色综合色国产| 男人的好看免费观看在线视频| 两个人的视频大全免费| 18禁黄网站禁片午夜丰满| 最好的美女福利视频网| 又爽又黄无遮挡网站| 男人舔女人下体高潮全视频| 九九爱精品视频在线观看| 色吧在线观看| 成人欧美大片| 免费大片18禁| 网址你懂的国产日韩在线| 性色avwww在线观看| 午夜福利在线观看免费完整高清在 | 成人无遮挡网站| 日本-黄色视频高清免费观看| 久久久国产成人免费| 欧美+亚洲+日韩+国产| 性欧美人与动物交配| 日韩精品有码人妻一区|