• <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片在线观看爽 | 国产精品免费视频内射| 成熟少妇高潮喷水视频| 性少妇av在线| 99精品久久久久人妻精品| 久久天躁狠狠躁夜夜2o2o| 黄色视频,在线免费观看| 9热在线视频观看99| 又黄又粗又硬又大视频| av一本久久久久| 多毛熟女@视频| 国产激情久久老熟女| 99久久精品国产亚洲精品| 91国产中文字幕| 久久精品亚洲精品国产色婷小说| 色尼玛亚洲综合影院| 亚洲精品国产色婷婷电影| 校园春色视频在线观看| 麻豆国产av国片精品| 久久午夜亚洲精品久久| 制服诱惑二区| 如日韩欧美国产精品一区二区三区| 中亚洲国语对白在线视频| 在线观看午夜福利视频| 欧美人与性动交α欧美精品济南到| 日本wwww免费看| 成年女人毛片免费观看观看9 | 咕卡用的链子| 国产精品99久久99久久久不卡| 久久人人97超碰香蕉20202| 在线免费观看的www视频| 精品国产乱码久久久久久男人| 国产国语露脸激情在线看| 精品久久久久久久久久免费视频 | 97人妻天天添夜夜摸| 热re99久久精品国产66热6| 日韩欧美国产一区二区入口| 男女午夜视频在线观看| 成人国产一区最新在线观看| 午夜福利视频在线观看免费| 99re在线观看精品视频| 18禁国产床啪视频网站| 十八禁网站免费在线| 男女床上黄色一级片免费看| 久久久久久久午夜电影 | 亚洲一码二码三码区别大吗| 1024视频免费在线观看| 人人妻人人澡人人看| 午夜久久久在线观看| 女人爽到高潮嗷嗷叫在线视频| 日韩免费av在线播放| 人妻一区二区av| 丝袜人妻中文字幕| 我的亚洲天堂| 亚洲精品中文字幕在线视频| 99热只有精品国产| av福利片在线| 午夜日韩欧美国产| 91成人精品电影| 人人妻人人澡人人爽人人夜夜| 在线观看日韩欧美| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 99国产精品一区二区蜜桃av | 午夜免费观看网址| 香蕉国产在线看| xxx96com| 亚洲精品中文字幕在线视频| 色94色欧美一区二区| 欧美日韩一级在线毛片| 看免费av毛片| 精品人妻1区二区| 露出奶头的视频| 99久久国产精品久久久| 亚洲中文日韩欧美视频| 国产欧美日韩一区二区三区在线| 俄罗斯特黄特色一大片| 99久久精品国产亚洲精品| 成人精品一区二区免费| 九色亚洲精品在线播放| 亚洲色图av天堂| 亚洲三区欧美一区| 久久久精品区二区三区| 国产在视频线精品| av电影中文网址| 亚洲在线自拍视频| 成人国语在线视频| 国产精品美女特级片免费视频播放器 | 在线观看日韩欧美| 性少妇av在线| 成人影院久久| 一本综合久久免费| 精品一区二区三区视频在线观看免费 | 亚洲成国产人片在线观看| 久久人妻av系列| 搡老熟女国产l中国老女人| 在线观看66精品国产| 精品电影一区二区在线| 一a级毛片在线观看| 亚洲熟妇熟女久久| 高清黄色对白视频在线免费看| 男人操女人黄网站| 亚洲精品乱久久久久久| 欧美久久黑人一区二区| 国产精品香港三级国产av潘金莲| 欧美日韩黄片免| 极品少妇高潮喷水抽搐| 99re在线观看精品视频| 亚洲成av片中文字幕在线观看| 91麻豆精品激情在线观看国产 | 亚洲中文av在线| 高清毛片免费观看视频网站 | 最新的欧美精品一区二区| 免费在线观看黄色视频的| 国产亚洲欧美98| 日日夜夜操网爽| 麻豆乱淫一区二区| 亚洲成人免费电影在线观看| 十八禁高潮呻吟视频| 久久精品熟女亚洲av麻豆精品| 中文字幕人妻丝袜一区二区| 亚洲欧美一区二区三区黑人| 99热国产这里只有精品6| 国产精品 欧美亚洲| 真人做人爱边吃奶动态| 国产精品一区二区精品视频观看| 久久人人97超碰香蕉20202| 国产91精品成人一区二区三区| 露出奶头的视频| 久久精品亚洲av国产电影网| 91麻豆av在线| 真人做人爱边吃奶动态| 国产精品99久久99久久久不卡| 午夜精品国产一区二区电影| 国产淫语在线视频| 叶爱在线成人免费视频播放| 欧美成狂野欧美在线观看| 高清av免费在线| 亚洲在线自拍视频| 美女国产高潮福利片在线看| 久久久久久久国产电影| 久热这里只有精品99| 女性生殖器流出的白浆| 国产91精品成人一区二区三区| 黑丝袜美女国产一区| 午夜福利在线观看吧| 国产欧美日韩一区二区精品| 国产精华一区二区三区| 91国产中文字幕| 波多野结衣一区麻豆| 久久性视频一级片| 久久精品成人免费网站| 久久人妻福利社区极品人妻图片| 天堂√8在线中文| 人妻 亚洲 视频| 韩国精品一区二区三区| 精品国产超薄肉色丝袜足j| 亚洲三区欧美一区| 999久久久精品免费观看国产| 热re99久久精品国产66热6| 国产在线观看jvid| 亚洲av电影在线进入| 色精品久久人妻99蜜桃| 精品一区二区三区四区五区乱码| 午夜日韩欧美国产| 91字幕亚洲| 久久精品国产清高在天天线| 亚洲av熟女| 久久久久久亚洲精品国产蜜桃av| 久久国产亚洲av麻豆专区| 国产无遮挡羞羞视频在线观看| 国产在线观看jvid| 久久99一区二区三区| 国产精品亚洲av一区麻豆| 一区二区三区国产精品乱码| 超碰97精品在线观看| 人妻 亚洲 视频| 亚洲国产欧美一区二区综合| 日韩一卡2卡3卡4卡2021年| 国产日韩一区二区三区精品不卡| 国产精品98久久久久久宅男小说| 亚洲精品在线观看二区| 岛国毛片在线播放| 后天国语完整版免费观看| 国产人伦9x9x在线观看| 成年版毛片免费区| 精品国产一区二区三区四区第35| а√天堂www在线а√下载 | 美女视频免费永久观看网站| 精品国产一区二区三区四区第35| 亚洲视频免费观看视频| 亚洲av成人av| 国产精品一区二区在线观看99| 五月开心婷婷网| 中亚洲国语对白在线视频| 老鸭窝网址在线观看| 午夜福利免费观看在线| 日韩熟女老妇一区二区性免费视频| 国产成人精品无人区| 国产精品一区二区在线观看99| 母亲3免费完整高清在线观看| 久久精品国产99精品国产亚洲性色 | 欧美一级毛片孕妇| 在线观看免费高清a一片| 久热爱精品视频在线9| 午夜视频精品福利| 国产一区二区三区视频了| 少妇粗大呻吟视频| 国产99白浆流出| 国产精品1区2区在线观看. | 精品亚洲成a人片在线观看| av有码第一页| 欧美激情高清一区二区三区| 满18在线观看网站| 很黄的视频免费| 男女高潮啪啪啪动态图| 午夜91福利影院| 伊人久久大香线蕉亚洲五| 久久精品人人爽人人爽视色| 视频区图区小说| 国产精品电影一区二区三区 | cao死你这个sao货| 亚洲精品乱久久久久久| 搡老熟女国产l中国老女人| 嫩草影视91久久| 岛国毛片在线播放| 亚洲熟女精品中文字幕| 国内久久婷婷六月综合欲色啪| 97人妻天天添夜夜摸| 久久亚洲精品不卡| 国产精品久久久久久人妻精品电影| 免费不卡黄色视频| 99精品在免费线老司机午夜| 老汉色∧v一级毛片| 亚洲精品美女久久av网站| 久久香蕉激情| av国产精品久久久久影院| 亚洲精品在线美女| 午夜福利,免费看| 18在线观看网站| 国产精品一区二区在线不卡| 97人妻天天添夜夜摸| 91麻豆精品激情在线观看国产 | 国产精品二区激情视频| 亚洲精品美女久久久久99蜜臀| 99re6热这里在线精品视频| 亚洲人成77777在线视频| 满18在线观看网站| 国产精华一区二区三区| 久久久精品国产亚洲av高清涩受| 欧美精品av麻豆av| 亚洲欧美一区二区三区久久| 老司机亚洲免费影院| 大型黄色视频在线免费观看| 国产精品影院久久| 午夜福利视频在线观看免费| 日韩 欧美 亚洲 中文字幕| 国产精品久久久人人做人人爽| 午夜福利,免费看| 一级,二级,三级黄色视频| 亚洲综合色网址| 欧美不卡视频在线免费观看 | 99久久精品国产亚洲精品| 国产在线精品亚洲第一网站| 国产成人免费观看mmmm| av不卡在线播放| 亚洲欧美精品综合一区二区三区| 人妻 亚洲 视频| 一区二区三区精品91| 人人澡人人妻人| 黑人巨大精品欧美一区二区mp4| 精品国产美女av久久久久小说| 成人手机av| 欧美国产精品一级二级三级| 无限看片的www在线观看| 亚洲伊人色综图| 狠狠狠狠99中文字幕| 很黄的视频免费| 91精品国产国语对白视频| 国产不卡一卡二| 欧美色视频一区免费| av国产精品久久久久影院| 亚洲专区中文字幕在线| 天堂中文最新版在线下载| 久久ye,这里只有精品| 国产深夜福利视频在线观看| 在线看a的网站| 老司机福利观看| 午夜两性在线视频| 一区福利在线观看| 成在线人永久免费视频| 欧美成人免费av一区二区三区 | 99久久人妻综合| 日本黄色日本黄色录像| 91麻豆av在线| 国产视频一区二区在线看| 午夜福利,免费看| 免费在线观看黄色视频的| 激情在线观看视频在线高清 | 老司机午夜福利在线观看视频| 精品一区二区三区视频在线观看免费 | 成人国语在线视频| 久久香蕉精品热| 美女高潮到喷水免费观看| 老司机午夜福利在线观看视频| 999精品在线视频| 男女之事视频高清在线观看| www.精华液| 婷婷成人精品国产| 亚洲九九香蕉| 日本wwww免费看| 午夜成年电影在线免费观看| 男人舔女人的私密视频| 18禁美女被吸乳视频| 中文字幕最新亚洲高清| av中文乱码字幕在线| 国产精品 国内视频| 久久国产精品人妻蜜桃| av中文乱码字幕在线| 亚洲五月婷婷丁香| 亚洲国产欧美网| 午夜影院日韩av| 精品福利永久在线观看| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看黄色视频的| 成人国语在线视频| 好看av亚洲va欧美ⅴa在| 亚洲精品在线观看二区| 99精品欧美一区二区三区四区| 母亲3免费完整高清在线观看| 两人在一起打扑克的视频| 亚洲在线自拍视频| 国产精品香港三级国产av潘金莲| 国产三级黄色录像| 大片电影免费在线观看免费| 动漫黄色视频在线观看| 亚洲成av片中文字幕在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 色综合婷婷激情| 日韩欧美一区二区三区在线观看 | 天堂动漫精品| 国产成人精品久久二区二区91| 国产精品欧美亚洲77777| 韩国av一区二区三区四区| 啦啦啦免费观看视频1| 中文字幕人妻丝袜一区二区| 交换朋友夫妻互换小说| 女人精品久久久久毛片| 日日爽夜夜爽网站| 日本黄色视频三级网站网址 | 中文字幕人妻熟女乱码| 国产精品久久电影中文字幕 | 日韩成人在线观看一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 欧美精品人与动牲交sv欧美| 亚洲情色 制服丝袜| 天天躁夜夜躁狠狠躁躁| 人成视频在线观看免费观看| 窝窝影院91人妻| 搡老岳熟女国产| 欧美乱妇无乱码| 成人永久免费在线观看视频| 亚洲片人在线观看| 欧美精品一区二区免费开放| 91老司机精品| 国产真人三级小视频在线观看| 香蕉国产在线看| 一本大道久久a久久精品| 久久精品aⅴ一区二区三区四区| 国产欧美日韩综合在线一区二区| a级毛片黄视频| 亚洲熟妇中文字幕五十中出 | 久久天躁狠狠躁夜夜2o2o| 黄色片一级片一级黄色片| 日日摸夜夜添夜夜添小说| 午夜精品国产一区二区电影| 少妇裸体淫交视频免费看高清 | 视频区欧美日本亚洲| 一区二区三区精品91| 夜夜躁狠狠躁天天躁| 涩涩av久久男人的天堂| 久久人妻av系列| 黄色毛片三级朝国网站| 91成人精品电影| 男女之事视频高清在线观看| 多毛熟女@视频| 国产精品一区二区精品视频观看| 国产免费av片在线观看野外av| 国产成人精品久久二区二区91| 91大片在线观看| 国产99白浆流出| a级片在线免费高清观看视频| 在线观看免费日韩欧美大片| 欧美大码av| 亚洲五月色婷婷综合| 国产精品国产av在线观看| а√天堂www在线а√下载 | 国产精品98久久久久久宅男小说| 黄片小视频在线播放| 香蕉国产在线看| 色播在线永久视频| ponron亚洲| 日韩视频一区二区在线观看| 一级片'在线观看视频| 精品国产超薄肉色丝袜足j| 黑丝袜美女国产一区| 欧美中文综合在线视频| 久久人人97超碰香蕉20202| 亚洲五月天丁香| 一级毛片高清免费大全| 成年动漫av网址| 免费观看a级毛片全部| 午夜激情av网站| 亚洲性夜色夜夜综合| 精品第一国产精品| 亚洲国产看品久久| 在线观看免费午夜福利视频| 好看av亚洲va欧美ⅴa在| 午夜免费观看网址| 亚洲欧美激情综合另类| 女人被躁到高潮嗷嗷叫费观| 1024视频免费在线观看| 国产精品乱码一区二三区的特点 | x7x7x7水蜜桃| 高清欧美精品videossex| 久久久久精品人妻al黑| 成人三级做爰电影| 热re99久久精品国产66热6| 亚洲七黄色美女视频| 亚洲一区二区三区不卡视频| 久久精品人人爽人人爽视色| 成人18禁高潮啪啪吃奶动态图| 无遮挡黄片免费观看| 国产欧美日韩一区二区精品| 男女午夜视频在线观看| 国产成人精品在线电影| 久久ye,这里只有精品| 极品少妇高潮喷水抽搐| 亚洲色图 男人天堂 中文字幕| 亚洲专区字幕在线| 精品国产美女av久久久久小说| 三级毛片av免费| 97人妻天天添夜夜摸| 成人免费观看视频高清| 久久人妻av系列| 亚洲片人在线观看| 女人被狂操c到高潮| 一区二区三区激情视频| 久久中文字幕一级| 法律面前人人平等表现在哪些方面| 中文欧美无线码| 日韩 欧美 亚洲 中文字幕| a在线观看视频网站| 亚洲av片天天在线观看| 亚洲精品久久成人aⅴ小说| 制服诱惑二区| 亚洲五月色婷婷综合| 亚洲欧美一区二区三区黑人| 热re99久久国产66热| 久久国产精品男人的天堂亚洲| 国产成人精品久久二区二区91| 脱女人内裤的视频| 91麻豆精品激情在线观看国产 | 老鸭窝网址在线观看| 久久精品亚洲精品国产色婷小说| 久久香蕉国产精品| 精品国产乱子伦一区二区三区| 久久精品国产清高在天天线| 丝瓜视频免费看黄片| 日韩大码丰满熟妇| 又黄又爽又免费观看的视频| 国产精品一区二区在线观看99| xxxhd国产人妻xxx| 久久久精品免费免费高清| 可以免费在线观看a视频的电影网站| 国产成人欧美| 麻豆av在线久日| 女性被躁到高潮视频| 亚洲av日韩精品久久久久久密| 久久ye,这里只有精品| 国产在线观看jvid| 女警被强在线播放| 午夜两性在线视频| 性色av乱码一区二区三区2| 久久精品国产综合久久久| 欧美日韩中文字幕国产精品一区二区三区 | 欧美成人免费av一区二区三区 | 国产有黄有色有爽视频| 久久人妻av系列| 日韩精品免费视频一区二区三区| 亚洲色图 男人天堂 中文字幕| 日本黄色视频三级网站网址 | 高清欧美精品videossex| 人人妻人人澡人人爽人人夜夜| 欧美精品av麻豆av| 免费看a级黄色片| 两性夫妻黄色片| 国产精品免费视频内射| 中文字幕另类日韩欧美亚洲嫩草| 夜夜躁狠狠躁天天躁| 久久精品国产亚洲av香蕉五月 | 麻豆乱淫一区二区| 老鸭窝网址在线观看| 亚洲va日本ⅴa欧美va伊人久久| 免费观看人在逋| 女人爽到高潮嗷嗷叫在线视频| 人人妻人人澡人人爽人人夜夜| 看黄色毛片网站| 国产欧美日韩一区二区三| 成人三级做爰电影| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 免费日韩欧美在线观看| 国产精品九九99| 国产日韩一区二区三区精品不卡| 免费黄频网站在线观看国产| 国产成人av教育| 亚洲国产欧美日韩在线播放| 亚洲人成电影观看| 亚洲专区国产一区二区| 亚洲国产欧美一区二区综合| 国产精品一区二区精品视频观看| 51午夜福利影视在线观看| 国产免费男女视频| 大陆偷拍与自拍| 看黄色毛片网站| videosex国产| 日日摸夜夜添夜夜添小说| 日韩欧美一区二区三区在线观看 | 看片在线看免费视频| 亚洲精品美女久久av网站| 国产aⅴ精品一区二区三区波| 久久精品国产清高在天天线| 18禁黄网站禁片午夜丰满| 欧美精品啪啪一区二区三区| 国产精品国产av在线观看| 久久午夜亚洲精品久久| 热re99久久国产66热| 久久精品熟女亚洲av麻豆精品| 成人手机av| 美国免费a级毛片| 黄片小视频在线播放| 午夜精品久久久久久毛片777| 欧美日韩国产mv在线观看视频| 成年动漫av网址| 欧美日韩视频精品一区| 色播在线永久视频| 亚洲成人免费av在线播放| cao死你这个sao货| 丰满人妻熟妇乱又伦精品不卡| 妹子高潮喷水视频| 国产精品一区二区在线观看99| 久久午夜综合久久蜜桃| 久久久久久久国产电影| 香蕉丝袜av| 久久久久久久久免费视频了| av天堂在线播放| 欧美乱妇无乱码| 激情在线观看视频在线高清 | 中文字幕人妻熟女乱码| 一二三四在线观看免费中文在| 黑丝袜美女国产一区| 国产精品一区二区免费欧美| 在线永久观看黄色视频| 在线av久久热| 女人高潮潮喷娇喘18禁视频| 十分钟在线观看高清视频www| 国产野战对白在线观看| 欧美黄色淫秽网站| 一区福利在线观看| 欧美国产精品一级二级三级| 久热这里只有精品99| 高清在线国产一区| 在线永久观看黄色视频| 国产精品一区二区在线不卡| 人妻一区二区av| 国产99久久九九免费精品| 国产欧美亚洲国产| 在线观看免费高清a一片| 亚洲性夜色夜夜综合| av不卡在线播放| 国产精品 国内视频| 在线观看www视频免费| 久99久视频精品免费| 新久久久久国产一级毛片| 亚洲精品美女久久久久99蜜臀| 国产精品一区二区在线不卡| 国产午夜精品久久久久久| 午夜91福利影院| 人成视频在线观看免费观看| 国产亚洲精品久久久久久毛片 | 亚洲 欧美一区二区三区| 免费高清在线观看日韩| 一本综合久久免费| 一级片免费观看大全| 中文亚洲av片在线观看爽 | 亚洲中文字幕日韩| 国产成人影院久久av| 国产精品98久久久久久宅男小说| 久久精品成人免费网站| 在线观看www视频免费| 久久久精品国产亚洲av高清涩受| 国产真人三级小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 看免费av毛片| 国产精品二区激情视频| 亚洲av成人一区二区三| 狠狠婷婷综合久久久久久88av| 老司机在亚洲福利影院| 男人的好看免费观看在线视频 | 亚洲 欧美一区二区三区| 女人久久www免费人成看片| 国产免费男女视频| 久久久久国内视频| 精品电影一区二区在线|