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

    外來入侵植物腫柄菊在云南的擴(kuò)散風(fēng)險(xiǎn)研究

    2021-08-16 00:50:10陳劍王四海朱楓袁春明
    廣西植物 2021年6期
    關(guān)鍵詞:擴(kuò)散環(huán)境變量

    陳劍 王四海 朱楓 袁春明

    摘 要: 預(yù)測外來入侵植物的潛在入侵范圍以及對不同環(huán)境的適應(yīng)性,能夠?yàn)橥鈦碇参锶肭痔峁╊A(yù)警和預(yù)防。該文通過實(shí)地調(diào)查腫柄菊(Tithonia diversifolia)在云南的入侵范圍,以及收集世界其他地區(qū)分布點(diǎn),結(jié)合12個(gè)關(guān)鍵環(huán)境變量,利用MaxEnt模型對腫柄菊適生區(qū)范圍和環(huán)境限制因子進(jìn)行分析。結(jié)果表明:在世界范圍適生區(qū)分布格局中,腫柄菊在云南無高適生區(qū),次適生區(qū)總面積占全省總面積的3.58%,中適生區(qū)占32.90%,低適生區(qū)占42.80%,非適生區(qū)占20.72%;有擴(kuò)散風(fēng)險(xiǎn)的區(qū)域主要為中適生區(qū),普遍分布于云南省26° N以南地區(qū);實(shí)地調(diào)查腫柄菊分布最高海拔2 200 m,在海拔1 000~1 500 m分布最為密集,而模型預(yù)測顯示海拔2 000 m以下都適宜腫柄菊生長,且500~2 000 m為最適宜海拔范圍;在12個(gè)環(huán)境變量中年降雨量、人類足跡、年均溫、海拔等環(huán)境變量對腫柄菊的分布起到了最重要的作用。綜上結(jié)果表明,目前在云南腫柄菊零星的分布點(diǎn)已到達(dá)其潛在分布區(qū)的邊緣,但已入侵區(qū)域還遠(yuǎn)小于潛在可入侵區(qū)域。為了預(yù)防腫柄菊在我國南方地區(qū)大面積擴(kuò)散帶來危害,有必要對其入侵和擴(kuò)散態(tài)勢進(jìn)行監(jiān)測。

    關(guān)鍵詞: 外來入侵植物, 腫柄菊, MaxEnt模型, 環(huán)境變量, 擴(kuò)散

    中圖分類號: Q948

    文獻(xiàn)標(biāo)識碼: A

    文章編號: 1000-3142(2021)05-0789-10

    Risk evaluation of Tithonia diversifolia dispersal in Yunnan Province,China

    CHEN Jian1,2, WANG Sihai1,2*, ZHU Feng3, YUAN Chunming 1,2

    ( 1. Yunnan Academy of Forestry and Grassland, Kunming 650201, China; 2. Yunnan Key Laboratory of Forest Plant Cultivation and Utilization /State Forestry Administration Key Laboratory of Yunnan Rare and Endangered Species Conservation and Propagation, Kunming 650201, China; 3. College of Geography and Ecotourism, Southwest Forestry University, Kunming 650224, China )

    Abstract:? Predicting the potential range of invasive alien plants and their adaptabilities to different environments could give early warning and prevention measures. The diverse bioclimatic types in Yunnan province would provide rich potential habitats for alien invasive plants. Tithonia diversifolia, an alien plant invading in tropical and sub-tropical regions of Yunnan, has quickly dispersed in recent years. In order to evaluate the risk of its invasion in different geographical environments, we investigated the invasion situation of T. diversifolia in Yunnan Province, and collected the distribution samples in the world. Using this species distribution coordinates and 12 main environmental variables, we analyzed the suitability of T. diversifolia to different environments and the dominating environmental factors by MaxEnt model. The results were as follows: Yunnan had not high-suitability area in the worldwide pattern of T. diversifolias distribution, meanwhile sub-, medium-, low-, un-suitability area accounted for 3.58%, 32.90%, 42.80% and 20.72% of the provincial area, respectively; Medium-suitability area was the main risk region, which widely interspersed in the south of 26° N in Yunnan; Field investigation showed T. diversifolia mainly distributed at the altitude of 1 000-1 500 m and its upper altitudinal limit was 2 200 m, while the model predicted that the altitude region below 2 000 m was suitable, especially between 500-2 000 m; The most important environmental factors affecting distribution of T. diversifolia were annual precipitation, human footprint, annual average temperature and altitude. In conclusion, although T. diversifolia has reached sporadically its suitable geographical limiting boundary at present, the invaded area is far less than the potential suitable area. It is necessary to keep monitoring T. diversifolias dispersal dynamic in the future, so that hazard caused by T. diversifolia could be controlled in time.

    Key words: alien invasive plant, Tithonia diversifolia, MaxEnt model, environmental variables, dispersal

    腫柄菊(Tithonia diversifolia)為菊科腫柄菊屬一年生草本植物,原產(chǎn)墨西哥及中美洲地區(qū),曾作為觀賞植物、綠肥和防治土壤侵蝕植物被廣泛引種到亞洲、非洲和澳洲的許多國家和地區(qū)。目前在東南亞、非洲、太平洋的一些熱帶和亞熱帶地區(qū)成為入侵雜草(Yang et al., 2012;Jama et al., 2000)。腫柄菊在我國的云南、廣西、廣東、福建、海南、臺灣等地已普遍分布(王四海等,2004; 羅瑛等,2009)。在云南,近年來隨著高速公路網(wǎng)的建設(shè),腫柄菊的擴(kuò)散速度明顯加快,有向更高緯度和海拔擴(kuò)散的趨勢(朱楓等,2018)。腫柄菊在熱帶和亞熱帶地區(qū)呈灌木或亞灌木狀生長,種子千粒重較大(4.5~6.5 g),種子的遠(yuǎn)距離傳播要依附于交通工具、人畜或其他媒介,這種具有特定路線和方式的傳播特點(diǎn)增加了其侵占新領(lǐng)地的時(shí)滯(王四海等,2008),從而可能導(dǎo)致腫柄菊的當(dāng)前實(shí)際分布和潛在適生區(qū)具有較大的空間差異。另外,云南由于受復(fù)雜地形和大氣環(huán)流的影響,各種氣候類型呈現(xiàn)出相互鑲嵌的特點(diǎn),影響植物分布的兩個(gè)主要因素(溫度和降雨)組合類型多樣(陳宗瑜,2001),在這種環(huán)境條件下腫柄菊的適應(yīng)性會(huì)隨著環(huán)境改變而不斷發(fā)生著變化(王四海等,2019)。當(dāng)前云南從北熱帶到北亞熱帶氣候類型均有腫柄菊分布,但是在不同的氣候區(qū)域群落面積分布極不均勻,在北亞熱帶僅為點(diǎn)狀的零星分布(朱楓等,2018)。擴(kuò)散時(shí)滯和環(huán)境適應(yīng)性都可能引起腫柄菊群落地理空間分布的不均勻,對腫柄菊在云南多變環(huán)境下的適應(yīng)性研究能確定其入侵潛能。

    腫柄菊在引種擴(kuò)散過程中已適應(yīng)多種生態(tài)氣候,在不同的地區(qū)限制其分布的環(huán)境因子會(huì)有差異。腫柄菊在100年前引入非洲的烏干達(dá),現(xiàn)在已擴(kuò)散到非洲的20多個(gè)國家,主要分布在撒哈拉以南降雨較多的區(qū)域,向北的繼續(xù)擴(kuò)散受到干旱環(huán)境的限制(Obiakara & Fourcade, 2018);在原產(chǎn)地墨西哥,腫柄菊在北部的分布會(huì)受到冬季霜凍影響(Lentz et al., 2008)。云南多樣化的地理環(huán)境可能會(huì)導(dǎo)致腫柄菊的適應(yīng)性在空間分布更為復(fù)雜。利用生態(tài)位模型或生境適應(yīng)性模型探討引種植物對新環(huán)境的適應(yīng)能力是近些年常使用的方法,通過模擬物種在環(huán)境中的生態(tài)位,然后把物種分布區(qū)投射到地理空間上,進(jìn)而實(shí)現(xiàn)預(yù)測物種的潛在分布范圍(Elith et al., 2010)。在眾多植物分布預(yù)測模型中,MaxEnt模型運(yùn)用最為普遍,常用來預(yù)測外來入侵植物潛在分布區(qū)(許仲林等,2015;劉曉彤等,2019)。本研究通過MaxEnt模型,結(jié)合氣候、土壤、海拔、人類影響指數(shù)(人類足跡)等多種環(huán)境因子對腫柄菊在云南的適應(yīng)性進(jìn)行分析,擬確定腫柄菊在不同區(qū)域的潛在適應(yīng)能力和將來的入侵范圍,以及影響腫柄菊分布的主要環(huán)境因子。

    1 材料與方法

    1.1 腫柄菊分布信息采集和分析

    根據(jù)腫柄菊易于沿公路擴(kuò)散和分布的特點(diǎn)(王四海等,2004),在2016年末和2017年初腫柄菊開花結(jié)實(shí)期,群落外貌特征具有較遠(yuǎn)距離能明顯識別的季節(jié),沿云南主要國道線和省道線沿路調(diào)查腫柄菊分布情況。依據(jù)群落為斑塊分布的特點(diǎn)(王四海等,2004),記錄斑塊(分布點(diǎn))的經(jīng)緯度、海拔信息,并采用目測法快速評估斑塊面積,調(diào)查范圍覆蓋腫柄菊在云南的分布區(qū)域。根據(jù)調(diào)查數(shù)據(jù),從分布最低點(diǎn)起把每50 m高差內(nèi)的腫柄菊群落面積進(jìn)行累加,所獲得的各海拔段群落面積用Oringin 9.1軟件擬合其分布曲線,獲得腫柄菊沿海拔分布規(guī)律。

    1.2 腫柄菊適生區(qū)分析

    (1)軟件及使用

    腫柄菊適生區(qū)分析利用MaxEnt V3.4.1和arcgis 10.2。MaxEnt用于物種適生區(qū)建模運(yùn)算,輸入數(shù)據(jù)為物種地理分布數(shù)據(jù)、海拔、氣候變量、土壤變量和人類足跡。模型重復(fù)運(yùn)行類別為交叉驗(yàn)證,樣本分布數(shù)據(jù)共分為10組,每一組數(shù)據(jù)均用于建模和交叉驗(yàn)證(Urbani et al., 2017),重復(fù)次數(shù)為5,運(yùn)算后輸出ASCII格式的物種適生區(qū)預(yù)測結(jié)果,用arcgis 10.2將其轉(zhuǎn)換成柵格圖層并按生境適宜度(0~1.00)進(jìn)行分級,獲得腫柄菊分布的5個(gè)適宜性等級,即高適生區(qū)(0.80~1.00)、次適生區(qū)(0.60~0.80)、中適生區(qū)(0.40~0.60)、低適生區(qū)(0.20~0.40)、非適生區(qū)(0~0.20),并按該分級標(biāo)準(zhǔn)對其進(jìn)行重分類后計(jì)算各級適生區(qū)面積。模型的預(yù)測準(zhǔn)確性以模型計(jì)算的ROC曲線的AUC值(曲線與橫坐標(biāo)所圍成的面積)進(jìn)行評價(jià),AUC 評價(jià)指標(biāo)(Swets, 1988):很差(≤0.060);較差(0.600~0.700);一般(0.700~0.800);良好(0.800~0.900);極好(0.900~1.000)。

    (2)環(huán)境數(shù)據(jù)來源與處理

    腫柄菊參與建模數(shù)據(jù)來自兩部分,第一部分為從全球生物多樣性信息網(wǎng)站(https://www.gbif.org/)、中國國家標(biāo)本資源平臺(www.nsii.org.cn)、中國數(shù)字植物標(biāo)本館(http://www.cvh.ac.cn/)等平臺上查詢世界范圍內(nèi)腫柄菊分布點(diǎn),并除去重復(fù)、錯(cuò)誤的分布點(diǎn);第二部分來自本研究實(shí)地采集的群落樣點(diǎn)。為防止模型過度擬合,以上兩部分分布點(diǎn)數(shù)據(jù)每5 km范圍內(nèi)只保留一個(gè)分布點(diǎn)(Veloz, 2009),最終共采用1 850個(gè)世界分布點(diǎn)(含云南省范圍內(nèi)的3個(gè)分布點(diǎn)),460個(gè)云南省實(shí)地采集分布點(diǎn),統(tǒng)一轉(zhuǎn)換為十進(jìn)制小數(shù)格式的.CSV文件用于建模。

    世界和中國云南行政邊界矢量圖來自國家地理信息系統(tǒng)網(wǎng)站(http://ngcc.sbsm.gov.cn/);19個(gè)典型氣候變量和海拔數(shù)據(jù)來自世界氣候數(shù)據(jù)庫(www.worldclim.org),氣候數(shù)據(jù)為1960—1990年共30 a的平均值,空間分辨率為30 arc sec。因19個(gè)氣候變量為非獨(dú)立變量,為提高模型模擬準(zhǔn)確度,用arcgis 10.0 多元分析工具各變量進(jìn)行相關(guān)性分析,對Pearson相關(guān)系數(shù)大于0.8的某一組數(shù)據(jù)集僅保留其中一個(gè)(Xu et al., 2019),最終使用的氣候變量為BIO1、BIO2、BIO10、BIO12、BIO15、BIO17、BIO18、BIO19等8個(gè)(表1);人類足跡數(shù)據(jù)來自哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)及應(yīng)用中心(https://sedac.ciesin.columbia.edu),空間分辨率為30 arc sec;土壤變量來自世界土壤數(shù)據(jù)庫(http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/),空間分辨率為30 arc sec,共34個(gè)土壤變量,參與MaxEnt建模后,僅列出參與建模重要值大于1的兩個(gè)變量S_GRAVEL和S_OC(表1)。以上所有環(huán)境數(shù)據(jù)均以云南省邊界進(jìn)行裁剪,并統(tǒng)一轉(zhuǎn)換為ASCII格式。最終參與建模的環(huán)境變量共12個(gè),各變量具體含義如表1。

    為了分析腫柄菊隨海拔升高的適應(yīng)性變化,把全省面積根據(jù)海拔段<500 m、500~1 000 m、1 000~2 000 m、2 000~2 500 m、2 500~3 000 m、>3 000 m進(jìn)行分類,根據(jù)MaxEnt適生區(qū)分析結(jié)果,計(jì)算在不同海拔段每類適生區(qū)面積占各自海拔段總面積的百分比。

    2 結(jié)果與分析

    2.1 腫柄菊在云南沿海拔分布現(xiàn)狀

    野外調(diào)查結(jié)果顯示, 當(dāng)前腫柄菊在云南分布最低海拔約在100 m,最高海拔為2 200 m,在1 000~1 500 m海拔段分布最為普遍(圖1)。共調(diào)查到的腫柄菊群落面積累計(jì)為1 965 548 m2,1 000~1 500 m范圍內(nèi)群落面積共計(jì)1 574 556 m2,占調(diào)查群落總面積的80.11%??梢娫? 000 m以下和1 500 m以上海拔段腫柄菊有著更多的潛在分布區(qū)。

    2.2 腫柄菊在云南適生區(qū)預(yù)測

    2.2.1 腫柄菊適生區(qū)分布格局 根據(jù)腫柄菊在全球適生狀況預(yù)測,在云南沒有高適生區(qū),次適生區(qū)主要零星分布在云南中部,至西南部和南部的區(qū)域,以及西北部的河谷地帶,面積約占全省總面積的3.58%;中適生區(qū)除滇西北和滇東北較高海拔區(qū)域外,大面積連續(xù)分布在云南的其余區(qū)域,占全省總面積的32.90%;低適生區(qū)間斷分布于全省各地,在云南中部和南部地區(qū)主要位于海拔相對較高區(qū)域,在滇西北和滇東北主要位于海拔相對較低區(qū)域,低適生區(qū)占全省總面積的42.80%;非適生區(qū)主要分布在滇西北和滇東北,占全省總面積的20.72%(表2,圖2)。腫柄菊在云南適生區(qū)域總體上隨海拔高差而變化,不同的適生區(qū)呈鑲嵌分布的特點(diǎn)。

    在500 m以下海拔段,低適生區(qū)占該海拔段的面積最大,為67.59%;其次為中適生區(qū),占29.43%;非適生區(qū)占2.94%。在500~1 000 m海拔段,次適生區(qū)、中適生區(qū)、低適生區(qū)和非適生區(qū)的面積分別占該海拔段總面積的4.28%、53.51%、38.09%和4.11%。在1 000~1 500 m海拔段,次適生區(qū)、中適生區(qū)、低適生區(qū)和非適生區(qū)的面積分別占該海拔段總面積的11.34%、49.71%、34.38%和4.56%(圖3)。在1 500~2 000 m海拔段,中適生區(qū)和低適生區(qū)在該海拔段所占面積變大,分別為49.71%和44.20%;次適生區(qū)和非適生區(qū)占比較小,分別為1.98%和6.04%。到了2 000~2 500 m海拔段,低適生區(qū)面積占比顯著變大,為73.50%。2 500~3000 m海拔段僅非適生區(qū)面積就高達(dá)74.83%,3 000 m以上海拔段非適生區(qū)面積已達(dá)100%(圖3)。在500~2 000 m海拔段次適生區(qū)和中適生區(qū)所占面積都在50%以上,該海拔段最適宜腫柄菊生長;2 000 m海拔段以上,從低適生區(qū)面積最大過渡到非適生區(qū)面積最大,大于2 500 m海拔高度腫柄菊就難以生長。500 m以下海拔段低適生區(qū)面積占比較大,這可能是因?yàn)樵谠颇?00 m以下地區(qū)主要是河谷深切區(qū)域,特殊地形限制了腫柄菊的生長。

    2.2.2 影響腫柄菊分布的主要環(huán)境因素 用于腫柄菊分布預(yù)測的12個(gè)環(huán)境變量中(表1),在腫柄菊適生區(qū)建模過程中,貢獻(xiàn)率大于1%的有6個(gè)環(huán)境因子,這6個(gè)環(huán)境因子的貢獻(xiàn)值占12個(gè)環(huán)境變量貢獻(xiàn)值總和的98.5%,貢獻(xiàn)值較多的環(huán)境因素為年降雨量、人類足跡、年均溫、海拔、最暖季度平均溫度、降雨季節(jié)性(表3)。從直接貢獻(xiàn)值看主要影響腫柄菊分布的環(huán)境因素主要是年降雨量、人類足跡、年均溫和海拔。置換重要值是將該特征變量值隨機(jī)置換為數(shù)模型運(yùn)算后獲得的值,可在建模貢獻(xiàn)率的基礎(chǔ)上對各變量貢獻(xiàn)進(jìn)行再次評估,該值越大,表明該變量影響越大。影響腫柄菊分布的置換重要值較高的為年均溫、人類足跡、最暖季度平均溫度、年降水量(表3)。

    綜合以上情況,年均溫、人類足跡、降水季節(jié)性、海拔、最暖季降水、最干季降水等6個(gè)環(huán)境變量對腫柄菊的分布起到了最重要的作用。圖4為這6個(gè)環(huán)境變量在單獨(dú)起作用時(shí)與腫柄菊出現(xiàn)概率的響應(yīng)曲線,當(dāng)各變量的值滿足年降水量>852 mm、人類足跡影響指數(shù)>10、年均溫17.2~23.7 ℃、海拔880~2 257 m、最暖季平均溫18.7~25.2 ℃、降雨季節(jié)性33.7~122.3時(shí),對腫柄菊的分布有最佳響應(yīng)(腫柄菊出現(xiàn)概率為0.5以上),也即為腫柄菊最適環(huán)境變量范圍。

    2.2.3 模型精度評價(jià) 圖5是檢測模型精度的ROC曲線(受試者工作特征曲線),AUC值為曲線與橫坐標(biāo)所圍的面積。AUC值越接近0.5(隨機(jī)模擬值),說明模型精度越低; AUC值越接近1, 則模型精度越高。本研究中重復(fù)測試的AUC平均值為0.883,平均標(biāo)準(zhǔn)差為0.002。說明建模預(yù)測效果良好。

    3 討論與結(jié)論

    3.1 腫柄菊的現(xiàn)有分布與潛在分布區(qū)

    選擇的12個(gè)環(huán)境變量中,年降雨量、人類足跡、年均溫、海拔、最暖季度平均溫度、降雨季節(jié)性等6個(gè)環(huán)境變量對腫柄菊的分布起到了最重要的作用。這種環(huán)境因素對腫柄菊分布的影響是基于對全球范圍的總體分析,在不同的區(qū)域影響腫柄菊分布的主導(dǎo)因素可能會(huì)有所不同。在云南北熱帶、南亞熱帶和中亞熱帶這些適合腫柄菊生長的氣候類型區(qū),傳播媒介(如公路網(wǎng))是影響腫柄菊分布的主導(dǎo)因素(朱楓等,2018);然而在北亞熱帶地區(qū),冬季的低溫顯著影響腫柄菊的結(jié)籽量,這一區(qū)域低的結(jié)籽量影響腫柄菊的進(jìn)一步擴(kuò)散(王四海等,2019)。在小尺度范圍的生境類型也影響著腫柄菊的分布(朱楓等,2018)。在未來腫柄菊擴(kuò)散過程中,由于局部影響分布主導(dǎo)因素的差異,即使適應(yīng)性相同的區(qū)域也可能存在擴(kuò)散蔓延速度不同。

    通過當(dāng)前腫柄菊在全球的分布點(diǎn)預(yù)測其潛在分布范圍,除了原產(chǎn)地中美洲和南美洲外,腫柄菊還可以大面積分布于非洲、亞洲,總體上高適生區(qū)分布范圍較少,預(yù)測分布區(qū)與已入侵區(qū)相對比,腫柄菊主要分布在次適生區(qū)和中適生區(qū)。云南當(dāng)前分布格局與世界的基本類似,已帶來嚴(yán)重生態(tài)危害的分布區(qū)域(王四海等,2004;朱楓等,2018;陳劍等,2020)也僅限于分布在次適生區(qū)和中適生區(qū)。由于腫柄菊在入侵的過程中氣候生態(tài)位較為保守(Obiakara & Fourcade, 2018), 在未來腫柄菊造成嚴(yán)重危害的區(qū)域應(yīng)主要是次適生區(qū)和中適生區(qū)。

    盡管當(dāng)前腫柄菊在云南分布的范圍很廣(跨越多于5個(gè)緯度,海拔變幅超過2 000 m),但只是在有些地區(qū)群落才高密度分布(朱楓等,2018),模型預(yù)測在500~2 000 m海拔段都較適宜腫柄菊生存,而當(dāng)前1 000~1 500 m海拔段腫柄菊分布最為密集。這表明現(xiàn)有的分布遠(yuǎn)沒有填滿潛在的分布區(qū)域,未來擴(kuò)散會(huì)在兩個(gè)層面進(jìn)行,一方面已經(jīng)有腫柄菊分布的地方,將會(huì)繼續(xù)擴(kuò)散和蔓延使群落分布更為密集;另一方面會(huì)向環(huán)境適宜區(qū)尚未有分布的區(qū)域擴(kuò)散,如沿怒江、瀾滄江、金沙江更為上游的河谷區(qū)域,以及云南東部的一些區(qū)域。

    3.2 腫柄菊在我國的分布潛力

    云南海拔高度從76.4 m到6 740 m,垂直高差達(dá)6 000 m以上。全省擁有從熱帶、亞熱帶、溫帶到寒帶等多種氣候類型,幾乎囊括我國所有的生態(tài)系統(tǒng)類型(楊宇明等,2008)。由于云南特殊的地勢地貌和多樣化的生物氣候類型,使得外來植物的各種氣候適應(yīng)性能力都能得到體現(xiàn)。本研究雖然以腫柄菊在云南的分布情況進(jìn)行分析,但其氣候的適應(yīng)性也可以作為在全國分布潛力的推斷依據(jù)。氣溫和降雨是影響植物分布的兩個(gè)主要因素,無論是模型預(yù)測還是當(dāng)前的實(shí)際分布情況,腫柄菊能分布最冷的區(qū)域?yàn)楸眮啛釒У貐^(qū),很顯然腫柄菊再向北更冷的區(qū)域擴(kuò)散時(shí)將會(huì)受到低溫的限制。腫柄菊在云南已有的分布區(qū)年降雨量最少區(qū)域僅為578 mm(王四海等,2019),我國600 mm等降雨量線大致沿青藏高原、黃土高原至華北地區(qū),這個(gè)范圍遠(yuǎn)超過腫柄菊所能向北分布的界限,從區(qū)域尺度看降雨不會(huì)是腫柄菊向北分布的限制因子。腫柄菊在云南能分布地區(qū)的最冷月份(一月)平均氣溫約8 ℃,如果以最冷月(一月)平均氣溫8 ℃推斷,腫柄菊在我國大致分布在廣西、廣東和福建及其以南地區(qū),這一分布范圍大致與世界較適宜分布區(qū)域沿我國南部地區(qū)分布范圍基本相當(dāng)。

    外來植物經(jīng)過對入侵地環(huán)境的長期適應(yīng)有時(shí)會(huì)發(fā)生生態(tài)位偏移,所以在某種程度上外來植物入侵范圍的精確預(yù)測具有困難(Guisan et al., 2014;李委濤等,2014)。原產(chǎn)于中美洲的紫莖澤蘭(Ageratina adenophora)、飛機(jī)草(Eupatorium odoratum)2種外來植物像腫柄菊一樣都是于20世紀(jì)30年代在云南建立了逃逸種群(薛紀(jì)如等,1979;王四海等,2008),當(dāng)前紫莖澤蘭已擴(kuò)散到西藏、四川、重慶和湖北等地(王翀等,2014),而飛機(jī)草僅限于在云南中南部、廣西西部和南部、廣東和福建南部等地區(qū)分布(楊波等,2009),兩種植物原產(chǎn)地和引入時(shí)間雖相似,但入侵潛力具有很大的不同。在21世紀(jì)初,飛機(jī)草在云南的入侵范圍與腫柄菊大致相同(王四海等,2004),但當(dāng)前腫柄菊入侵的海拔高度和向北的擴(kuò)散范圍都大于飛機(jī)草(朱楓等,2018),表現(xiàn)出較強(qiáng)的入侵態(tài)勢。紫莖澤蘭能耐受-5 ℃的低溫,飛機(jī)草在0 ℃時(shí)葉片和枝條可凍死(吳仁潤等,1984),腫柄菊根據(jù)目前分布情況推測在冬季可以忍耐短期0 ℃以下低溫,這三種植物在耐凍性方面的差異決定了它們潛在分布區(qū)的不同。另外,即使外來入侵植物沒有發(fā)生生態(tài)適應(yīng)性進(jìn)化,它們也可以占據(jù)與原產(chǎn)地不同的、非一致性的氣候生態(tài)位(Datta et al., 2019),這也可以解釋一些外來植物原產(chǎn)地相同,但入侵范圍可能不同的原因。綜合考慮到腫柄菊入侵氣候生態(tài)位保守性(Obiakara & Fourcade, 2018)、入侵時(shí)滯較長(王四海等,2008),以及全球氣溫變化因素,經(jīng)過較長時(shí)間的擴(kuò)散和適應(yīng),腫柄菊也有可能向北零星擴(kuò)散到一月平均溫度小于8 ℃的低適生區(qū)。雖已有標(biāo)本記載在廣西、廣東、福建、海南等地有腫柄菊分布,但對于其大面積密集分布并帶來危害性的報(bào)道并不多。為了預(yù)防腫柄菊在我國南方地區(qū)大面積擴(kuò)散帶來危害,有必要對其入侵和擴(kuò)散態(tài)勢進(jìn)行監(jiān)測(李委濤等,2014)。

    參考文獻(xiàn):

    CHEN J, WANG SH, YANG W, et al., 2020. Community dynamics of an alien invasive plant species Tithonia diversifolia? [J]. Chin J Ecol, 39(2):469-477.? [陳劍, 王四海, 楊衛(wèi), 等, 2020. 外來入侵植物腫柄菊群落動(dòng)態(tài)變化特征 [J]. 生態(tài)學(xué)雜志, 39(2): 469-477.]

    CHEN ZY, 2001. Climate of Yunnan? [M]. Beijing: China Meteorological Press.? [陳宗瑜, 2001. 云南氣候總論 [M]. 北京: 氣象出版社.]

    DATTA A, SCHWEIGER O, KHN I, 2019. Niche expansion of the invasive plant species Ageratina adenophora despite evolutionary constraints? [J]. J Biogeogr, 46: 1306-1315.

    ELITH J, KEARNEY M, PHILLIPS S, et al., 2010. The art of modeling range-shifting species? [J]. Mothods Ecol Evol? [J], 1: 330-342.

    GUISAN A, PETITPIERRE B, BROENNIMANN O, et al., 2014. Unifying niche shift studies: insights from biological invasions? [J]. Trend Ecol Evol, 29(5): 260-269.

    JAMA B, PALM CA, BURESH RJ, et al., 2000. Tithonia diversifolia as a green manure for soil fertility improvement in western Kenya: A review? [J]. Agrofor Syst, 49(2):201-221.

    LENTZ DL, BYE R, SNCHEZ-CORDERO V, 2008. Ecological niche modeling and distribution of wild sunflower in Mexico? [J]. Int J Plant Sci, 169(4): 541-549.

    LI WT, ZHENG YL, FENG YL, 2014. Differences in growth traits between Chromolaena odorata plants from invasive and native populations? [J]. Acta Ecol Sin, 34(23): 6890-6897.? [李委濤, 鄭玉龍, 馮玉龍, 2014.飛機(jī)草入侵種群與原產(chǎn)地種群生長性狀的差異 [J]. 生態(tài)學(xué)報(bào), 34(23): 6890-6897.]

    LIU XT, YUAN Q, NI J, 2019. Research advances in modeling plant species distribution in China? [J]. Chin J Plant Ecol, 43(4): 273-283.? [劉曉彤, 袁泉, 倪健, 2019. 中國植物分布模擬研究現(xiàn)狀 [J]. 植物生態(tài)學(xué)報(bào), 43(4): 273-283.]

    LUO Y, LIU Z, ZHAO JH, et al., 2009. Analysis of mineral elements and quality evaluation of Tithonia diversifolia(Hemsl.) A. Gray.? [J]. Trop Agric Eng, 33(3): 36-39.? [羅瑛, 劉壯, 趙君華, 等, 2009. 腫柄菊的礦質(zhì)營養(yǎng)元素分析及評價(jià) [J]. 熱帶農(nóng)業(yè)工程, 33(3):36-39.]

    OBIAKARA MC, FOURCADE Y, 2018. Climatic niche and potential distribution of Tithonia diversifolia (Hemsl.) A. Gray in Africa? [J]. PLoS ONE, 13(9): e0202421.

    SWETS JA, 1988. Measuring the accuracy of diagnostic systems? [J]. Science, 240(4857): 1285-1293.

    URBANI F, D′ALESSANDRO P, BIONDI M, 2017. Using Maximum Entropy Modeling (MaxEnt) to predict future trends in the distribution of high altitude endemic insects in response to climate change? [J]. B Insectol, 70(2): 189-200.

    VELOZ SD, 2009. Spatially autocorrelated sampling falsely inflates measures of accuracy for presence-only niche models [J]. J Biogeogr, 36(12): 2290-2299.

    WANG C, LIN HL, HE L, et al., 2014. Research on responses of Eupatorium adenophorums potential distribution to climate change? [J]. Acta Pratacult Sin, 23(4): 20-30.? [王翀, 林慧龍, 何蘭, 等, 2014. 紫莖澤蘭潛在分布對氣候變化響應(yīng)的研究 [J]. 草業(yè)科學(xué), 23(4): 20-30. ]

    WANG SH, CHEN J, LI NY, et al., 2019. Factors affecting seed production of Tithonia diversifolia, an invasive alien plant? [J]. Ecol Environ Sci, 28(7): 1369-1378.? [王四海, 陳劍, 李寧云, 等, 2019. 外來植物腫柄菊結(jié)籽量影響因素分析 [J]. 生態(tài)環(huán)境學(xué)報(bào), 28(7): 1369-1378.]

    WANG SH, SUN WB, CHENG X, 2004. Attributes of plant proliferation, geographic spread and the natural communities invaded by the naturalized alien plant species Tithonia diversifolia in Yunnan,China? [J]. Acta Ecol Sin, 24(3): 444-449.? [王四海, 孫衛(wèi)邦, 成曉, 2004. 逃逸外來植物腫柄菊在云南的生長繁殖特性、地理分布現(xiàn)狀及群落特征 [J]. 生態(tài)學(xué)報(bào), 24(3): 444-449.]

    WANG SH, SUN WB, CHENG X, et al., 2008. Reproductive characteristics of Tithonia diversifolia and its geographical spread in Yunnan Province of South-West China? [J]. Acta Ecol Sin, 28(3): 426-432.? [王四海, 孫衛(wèi)邦, 成曉, 等, 2008. 外來植物腫柄菊(Tithonia diversifolia)的繁殖特性及其地理擴(kuò)散 [J]. 生態(tài)學(xué)報(bào), 28(3): 426-432.]

    WU RR, ZHANG DY, LU XS, 1984. On the distribution, perniciousness and control of the two weeds,Chromolaena coelestina and C. odorata in Yunnan Province of China? [J]. Grassland Chin, 2: 17-22.? [吳仁潤, 張德銀, 盧欣石, 1984. 紫莖澤蘭和飛機(jī)草在云南的分布、危害與防治 [J]. 中國草原, 2: 17-22.]

    XU DP, ZHUO ZH, WANG RL, et al., 2019. Modeling the distribution of Zanthoxylum armatum in China with MaxEnt modeling? [J]. Global Ecol and Conserv, 19: e00691.

    XUE JR, DONG SC, YIN JH, et al., 1979. Adaptability and damage situation of Eupatorium adenophorum and control measures? [J]. Acta Bot Yunnan, 1(1): 106-114.? [薛紀(jì)如, 董世仁, 印嘉祜, 等, 1979. 紫莖澤蘭的適應(yīng)性能、危害情況及防除措施 [J]. 云南植物研究, 1(1): 106-114.]

    XU ZL, PENG HH, PENG SZ, 2015. The development and evaluation of species distribution models? [J]. Acta Ecol Sin, 35(2): 557-567.? [許仲林, 彭煥華, 彭守璋, 2015. 物種分布模型的發(fā)展及評價(jià)方法 [J]. 生態(tài)學(xué)報(bào), 35(2): 557-567.]

    YANG B, XUE YG, TANG XF, et al., 2009. Predicting potential geographic distribution of Eupatorium odoratum in China? [J]. Plant Protect, 35(4): 70-73.? [楊波, 薛躍規(guī), 唐小飛, 等, 2009. 外來入侵植物飛機(jī)草在中國的適生區(qū)預(yù)測 [J]. 植物保護(hù), 35(4): 70-73.]

    YANG J, TANG L, GUAN YL, et al., 2012. Genetic diversity of an alien invasive plant mexican sunflower (Tithonia diversifolia) in China? [J]. Weed Sci, 60(4): 552-557.

    YANG YM, WANG J, WANG JH, et al., 2008. Studies on the biodiversity and its conservation in Yunnan, China ?[M]. Beijing: Science Press.? [楊宇明, 王娟, 王建皓, 等, 2008. 云南生物多樣性及其保護(hù)研究 [M]. 北京: 科學(xué)出版社.]

    ZHU F, WANG SH, CHEN J, et al., 2018. Distribution characteristics and influencing factors of Tithonia diversifolia, an alien invasive plant species, in Yunnan Province, China? [J]. Chin J Ecol, 37(9): 27-34.? [朱楓, 王四海, 陳劍, 等, 2018. 外來入侵植物腫柄菊在云南的分布特征及其影響因素 [J]. 生態(tài)學(xué)雜志, 37(9): 27-34.]

    (責(zé)任編輯 何永艷)

    收稿日期:? 2020-03-06

    基金項(xiàng)目: 國家自然科學(xué)基金(31560092);云南省森林植物培育與開發(fā)利用重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目( ZZCX2016-10);云南省技術(shù)創(chuàng)新人才培養(yǎng)資助項(xiàng)目(2019HB066) [Supported by the National Natural Science Foundation of China (31560092); Open Fund of Yunnan Key Laboratory of Forest Plant Cultivation and Utilization (ZZCX2016-10); Talent Training Program on Technical Innovation of Yunnan Province (2019HB066)]。

    作者簡介: 陳劍(1978-),碩士,助理研究員,主要從事林業(yè)可持續(xù)發(fā)展研究,(E-mail)chenjianyaf@126.com。

    通信作者: 王四海,博士,副研究員,主要從事植物保護(hù)與利用研究,(E-mail)wangsh6688@163.com。

    猜你喜歡
    擴(kuò)散環(huán)境變量
    從桌面右鍵菜單調(diào)用環(huán)境變量選項(xiàng)
    徹底弄懂Windows 10環(huán)境變量
    混合無線多跳網(wǎng)絡(luò)安全性問題探討
    基于混沌理論的自適應(yīng)參數(shù)圖像加密算法
    科技資訊(2016年18期)2016-11-15 07:45:11
    雙船浮托安裝環(huán)境影響對比研究
    科技視界(2016年21期)2016-10-17 17:00:25
    基于三階段DEA—Malmquist模型的中國省域城鎮(zhèn)化效率測度及其收斂分析
    山東省產(chǎn)業(yè)集聚存在的問題與演進(jìn)路徑探討
    山東省港口發(fā)展效果評價(jià)
    一類具有時(shí)滯的非自治食物鏈捕食擴(kuò)散模型
    氣候變化約束下資源型企業(yè)最佳規(guī)模邊界的確定
    午夜日本视频在线| 成人午夜精彩视频在线观看| 久久草成人影院| 亚洲精品日本国产第一区| 毛片女人毛片| 成年人午夜在线观看视频 | 80岁老熟妇乱子伦牲交| 国产爱豆传媒在线观看| 亚洲内射少妇av| 高清欧美精品videossex| 成人亚洲精品av一区二区| 亚洲性久久影院| 91久久精品电影网| 又粗又硬又长又爽又黄的视频| 亚洲va在线va天堂va国产| 亚洲美女搞黄在线观看| 国产老妇女一区| 亚洲欧美成人综合另类久久久| 天美传媒精品一区二区| 久久久久久久大尺度免费视频| 亚洲,欧美,日韩| 日产精品乱码卡一卡2卡三| 国产精品av视频在线免费观看| 美女国产视频在线观看| a级毛片免费高清观看在线播放| 亚洲婷婷狠狠爱综合网| videos熟女内射| 国内精品一区二区在线观看| 插逼视频在线观看| 五月伊人婷婷丁香| 国产成人91sexporn| 免费看美女性在线毛片视频| 国产v大片淫在线免费观看| 色5月婷婷丁香| 日韩精品青青久久久久久| 2018国产大陆天天弄谢| 国产人妻一区二区三区在| 国产高清三级在线| 亚洲精品aⅴ在线观看| 国产精品一区二区三区四区久久| 午夜爱爱视频在线播放| 亚洲av中文字字幕乱码综合| 国产极品天堂在线| 日韩成人av中文字幕在线观看| 非洲黑人性xxxx精品又粗又长| 可以在线观看毛片的网站| 麻豆久久精品国产亚洲av| 国产成年人精品一区二区| 99久久精品国产国产毛片| 国产色爽女视频免费观看| 国产老妇女一区| 久久久久久久久大av| 国产精品久久久久久av不卡| 日日啪夜夜爽| 亚洲成人av在线免费| 国产成人免费观看mmmm| 国产亚洲精品久久久com| 国产 一区 欧美 日韩| 熟女电影av网| 丰满人妻一区二区三区视频av| 性色avwww在线观看| 尤物成人国产欧美一区二区三区| 永久网站在线| 久久久国产一区二区| 国产片特级美女逼逼视频| 夜夜爽夜夜爽视频| 2022亚洲国产成人精品| 少妇裸体淫交视频免费看高清| 国产69精品久久久久777片| 亚洲欧美成人精品一区二区| av免费观看日本| 国精品久久久久久国模美| 狠狠精品人妻久久久久久综合| www.av在线官网国产| 3wmmmm亚洲av在线观看| 成人无遮挡网站| 精品人妻一区二区三区麻豆| 成人毛片a级毛片在线播放| 床上黄色一级片| 国产亚洲精品久久久com| 国产v大片淫在线免费观看| 亚洲精品国产av成人精品| 七月丁香在线播放| 爱豆传媒免费全集在线观看| 久久人人爽人人爽人人片va| 国产高清国产精品国产三级 | 最新中文字幕久久久久| 欧美日韩国产mv在线观看视频 | 欧美 日韩 精品 国产| 欧美不卡视频在线免费观看| 亚洲伊人久久精品综合| 国产精品嫩草影院av在线观看| 午夜精品一区二区三区免费看| 欧美精品一区二区大全| 欧美一区二区亚洲| 中文字幕人妻熟人妻熟丝袜美| 免费观看av网站的网址| 国产视频首页在线观看| 久久久久久久国产电影| 高清午夜精品一区二区三区| 九色成人免费人妻av| 人妻少妇偷人精品九色| 熟妇人妻久久中文字幕3abv| 99久久中文字幕三级久久日本| 男插女下体视频免费在线播放| 高清在线视频一区二区三区| 18+在线观看网站| 肉色欧美久久久久久久蜜桃 | 亚洲四区av| 联通29元200g的流量卡| 精品人妻视频免费看| 女人十人毛片免费观看3o分钟| 国产久久久一区二区三区| 色网站视频免费| 亚洲av成人精品一二三区| 欧美一区二区亚洲| 少妇熟女欧美另类| 少妇的逼水好多| 久久99蜜桃精品久久| 欧美激情在线99| 高清在线视频一区二区三区| 色综合色国产| 亚洲国产精品专区欧美| 国产黄色视频一区二区在线观看| 国产成人freesex在线| 少妇被粗大猛烈的视频| 99久久精品国产国产毛片| 国产伦精品一区二区三区四那| 日产精品乱码卡一卡2卡三| 亚洲精品久久午夜乱码| 别揉我奶头 嗯啊视频| av.在线天堂| 精品久久久久久久久av| 亚洲精品中文字幕在线视频 | 国产综合精华液| 中文字幕制服av| 国产老妇女一区| 国产黄a三级三级三级人| 久久久久性生活片| 亚洲色图av天堂| 丰满人妻一区二区三区视频av| 久99久视频精品免费| 在线播放无遮挡| 亚洲乱码一区二区免费版| 在线观看免费高清a一片| 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| 18禁在线播放成人免费| av黄色大香蕉| 亚洲精品乱码久久久久久按摩| 日本免费a在线| 亚洲精品国产av蜜桃| 精品国产露脸久久av麻豆 | 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 亚洲欧美清纯卡通| 久久久久久久久久成人| 亚洲18禁久久av| 久久久久性生活片| 国产午夜福利久久久久久| 国产亚洲av嫩草精品影院| 国产成人精品一,二区| 亚洲av国产av综合av卡| 性色avwww在线观看| 国产在线男女| 26uuu在线亚洲综合色| 一本久久精品| 亚洲欧美成人精品一区二区| www.色视频.com| 简卡轻食公司| 日韩国内少妇激情av| 色吧在线观看| 国产中年淑女户外野战色| 天天一区二区日本电影三级| 久久久久久国产a免费观看| 午夜福利成人在线免费观看| 久久热精品热| 91久久精品国产一区二区成人| 乱码一卡2卡4卡精品| 亚洲精品亚洲一区二区| 亚洲国产最新在线播放| 97超碰精品成人国产| 久久精品夜夜夜夜夜久久蜜豆| 国内精品宾馆在线| 黄色配什么色好看| 91久久精品国产一区二区三区| 欧美高清性xxxxhd video| 午夜福利网站1000一区二区三区| 久久韩国三级中文字幕| 国产伦理片在线播放av一区| 我的女老师完整版在线观看| 成年免费大片在线观看| 午夜福利网站1000一区二区三区| 毛片一级片免费看久久久久| 久久亚洲国产成人精品v| www.色视频.com| 亚洲精品第二区| 国产精品一区二区在线观看99 | 国产日韩欧美在线精品| 久久久久精品久久久久真实原创| 亚洲精品第二区| 亚洲成色77777| 九色成人免费人妻av| 国产v大片淫在线免费观看| 欧美xxⅹ黑人| 搡老妇女老女人老熟妇| 国产在线男女| 99热网站在线观看| 亚洲精品中文字幕在线视频 | 18禁裸乳无遮挡免费网站照片| 午夜视频国产福利| 日韩不卡一区二区三区视频在线| 99热6这里只有精品| 久久久久久久久久人人人人人人| 亚洲av一区综合| 丝袜喷水一区| 午夜福利在线在线| 深爱激情五月婷婷| 亚洲国产成人一精品久久久| 人体艺术视频欧美日本| 九色成人免费人妻av| 麻豆av噜噜一区二区三区| 亚洲最大成人中文| 亚洲国产精品专区欧美| 2021少妇久久久久久久久久久| 成人综合一区亚洲| 日韩制服骚丝袜av| 欧美 日韩 精品 国产| 一区二区三区高清视频在线| 亚洲va在线va天堂va国产| 亚洲性久久影院| av又黄又爽大尺度在线免费看| 91精品国产九色| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜爱| 午夜免费观看性视频| 久久久成人免费电影| 成人性生交大片免费视频hd| 国产成人午夜福利电影在线观看| 欧美区成人在线视频| 日韩视频在线欧美| 国产女主播在线喷水免费视频网站 | 国产有黄有色有爽视频| 国产欧美另类精品又又久久亚洲欧美| 中文乱码字字幕精品一区二区三区 | 22中文网久久字幕| 嫩草影院新地址| 欧美成人午夜免费资源| 欧美日韩综合久久久久久| 精品午夜福利在线看| 人妻少妇偷人精品九色| 欧美潮喷喷水| 男的添女的下面高潮视频| 精品久久国产蜜桃| 国产成人aa在线观看| 少妇熟女欧美另类| 亚洲国产精品成人久久小说| 成年女人在线观看亚洲视频 | 欧美变态另类bdsm刘玥| 亚洲精品视频女| 久久久久久久午夜电影| 午夜爱爱视频在线播放| 亚洲三级黄色毛片| 精品熟女少妇av免费看| 美女黄网站色视频| 日韩欧美精品免费久久| 伊人久久国产一区二区| 九草在线视频观看| 2022亚洲国产成人精品| 在线观看美女被高潮喷水网站| 亚洲精品第二区| 熟妇人妻久久中文字幕3abv| 天天躁夜夜躁狠狠久久av| 噜噜噜噜噜久久久久久91| 啦啦啦啦在线视频资源| 日韩,欧美,国产一区二区三区| 3wmmmm亚洲av在线观看| 99久久精品国产国产毛片| 国产成人精品福利久久| 亚洲最大成人av| 国产乱人偷精品视频| 一区二区三区高清视频在线| 激情五月婷婷亚洲| 91精品国产九色| 国产伦在线观看视频一区| 日韩三级伦理在线观看| 色视频www国产| 国产毛片a区久久久久| 国产成人精品婷婷| 秋霞在线观看毛片| 国产av码专区亚洲av| 国产亚洲最大av| 国产91av在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 国产成人a区在线观看| 你懂的网址亚洲精品在线观看| 综合色av麻豆| 亚洲精品中文字幕在线视频 | 91精品伊人久久大香线蕉| 亚洲天堂国产精品一区在线| 成人毛片60女人毛片免费| 在线观看美女被高潮喷水网站| 成年av动漫网址| 日本黄大片高清| 噜噜噜噜噜久久久久久91| 免费看日本二区| 人妻制服诱惑在线中文字幕| 免费电影在线观看免费观看| 少妇熟女欧美另类| 岛国毛片在线播放| 精品午夜福利在线看| 亚洲国产精品成人久久小说| 99久久精品一区二区三区| 只有这里有精品99| 日日干狠狠操夜夜爽| 亚洲av在线观看美女高潮| 欧美最新免费一区二区三区| 亚洲av男天堂| 能在线免费观看的黄片| 免费av观看视频| 日韩欧美精品v在线| videossex国产| 丝袜喷水一区| 午夜福利在线观看免费完整高清在| 日韩三级伦理在线观看| 人妻少妇偷人精品九色| 欧美人与善性xxx| 日韩三级伦理在线观看| 亚洲一区高清亚洲精品| 成人亚洲精品一区在线观看 | 少妇猛男粗大的猛烈进出视频 | 国产乱人偷精品视频| av国产免费在线观看| 久久久久性生活片| 2021天堂中文幕一二区在线观| 最近的中文字幕免费完整| 男的添女的下面高潮视频| 日韩亚洲欧美综合| 大话2 男鬼变身卡| 老司机影院成人| 熟女电影av网| 一本久久精品| 日本免费在线观看一区| 亚洲最大成人手机在线| 永久免费av网站大全| 精品欧美国产一区二区三| 你懂的网址亚洲精品在线观看| 欧美日韩国产mv在线观看视频 | 免费黄网站久久成人精品| av.在线天堂| 一级二级三级毛片免费看| 如何舔出高潮| 亚洲国产高清在线一区二区三| 色综合亚洲欧美另类图片| 2021少妇久久久久久久久久久| 国产大屁股一区二区在线视频| 国产成人精品久久久久久| 床上黄色一级片| 建设人人有责人人尽责人人享有的 | 亚洲精品aⅴ在线观看| av一本久久久久| 女的被弄到高潮叫床怎么办| 人妻制服诱惑在线中文字幕| 亚洲一级一片aⅴ在线观看| 好男人在线观看高清免费视频| 亚洲在久久综合| 国产成人精品婷婷| 18禁动态无遮挡网站| 久久热精品热| 日韩强制内射视频| 十八禁国产超污无遮挡网站| 国产爱豆传媒在线观看| 99视频精品全部免费 在线| 国产亚洲最大av| 成人综合一区亚洲| 我的女老师完整版在线观看| 99热这里只有精品一区| 午夜福利在线在线| 大香蕉97超碰在线| 内射极品少妇av片p| 国产精品一区二区三区四区久久| 久热久热在线精品观看| 男女边摸边吃奶| 一级毛片我不卡| 国产成人freesex在线| 听说在线观看完整版免费高清| 少妇的逼水好多| 最近最新中文字幕免费大全7| 久久久久久久大尺度免费视频| 日韩精品青青久久久久久| 国内精品宾馆在线| 三级经典国产精品| 舔av片在线| 国产高清三级在线| 亚洲精品乱久久久久久| 中文字幕久久专区| or卡值多少钱| 熟妇人妻久久中文字幕3abv| 精品久久久久久久久av| 中文字幕亚洲精品专区| 国产在线一区二区三区精| 亚洲国产日韩欧美精品在线观看| 国产成年人精品一区二区| 又爽又黄a免费视频| 亚洲精品日韩av片在线观看| 国语对白做爰xxxⅹ性视频网站| 国产成人精品福利久久| 色吧在线观看| 成年人午夜在线观看视频 | 国产伦在线观看视频一区| 91aial.com中文字幕在线观看| 2021少妇久久久久久久久久久| 熟妇人妻久久中文字幕3abv| 色哟哟·www| 午夜激情福利司机影院| 欧美日韩一区二区视频在线观看视频在线 | 一个人看视频在线观看www免费| 欧美另类一区| 老师上课跳d突然被开到最大视频| 国产不卡一卡二| 秋霞在线观看毛片| 寂寞人妻少妇视频99o| 国产免费视频播放在线视频 | 青青草视频在线视频观看| 大又大粗又爽又黄少妇毛片口| 精品国产三级普通话版| 干丝袜人妻中文字幕| 国产黄色免费在线视频| 在线 av 中文字幕| 亚洲色图av天堂| 校园人妻丝袜中文字幕| 精品久久久久久久久久久久久| 国产黄频视频在线观看| 午夜爱爱视频在线播放| 成年女人看的毛片在线观看| 美女内射精品一级片tv| 一个人看视频在线观看www免费| 欧美性猛交╳xxx乱大交人| 可以在线观看毛片的网站| av又黄又爽大尺度在线免费看| 精品久久久久久久久久久久久| 精华霜和精华液先用哪个| 亚洲av电影不卡..在线观看| 亚州av有码| 国产成人精品久久久久久| 成人av在线播放网站| av在线观看视频网站免费| 99热全是精品| 日韩精品青青久久久久久| av女优亚洲男人天堂| 国产黄a三级三级三级人| 欧美日本视频| 久久这里只有精品中国| 欧美日韩亚洲高清精品| 男女啪啪激烈高潮av片| 久久久久精品性色| 网址你懂的国产日韩在线| 啦啦啦中文免费视频观看日本| 久久亚洲国产成人精品v| 一级毛片电影观看| 大又大粗又爽又黄少妇毛片口| 国产黄频视频在线观看| 国产国拍精品亚洲av在线观看| 综合色丁香网| 一个人看视频在线观看www免费| av在线播放精品| 亚洲在线自拍视频| 国产成年人精品一区二区| 免费观看精品视频网站| 亚洲精品乱久久久久久| 国产探花在线观看一区二区| 精品久久国产蜜桃| 国产激情偷乱视频一区二区| 男女边吃奶边做爰视频| 中文字幕人妻熟人妻熟丝袜美| 一级毛片 在线播放| 最近最新中文字幕大全电影3| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩av不卡免费在线播放| 精品人妻熟女av久视频| 国产高潮美女av| 亚洲国产色片| 肉色欧美久久久久久久蜜桃 | 国产精品一二三区在线看| 亚洲精品亚洲一区二区| 51国产日韩欧美| 天天躁夜夜躁狠狠久久av| 能在线免费观看的黄片| 国内少妇人妻偷人精品xxx网站| 亚洲国产色片| 中文天堂在线官网| 国产人妻一区二区三区在| 纵有疾风起免费观看全集完整版 | 国产精品人妻久久久久久| 激情五月婷婷亚洲| 国产午夜精品一二区理论片| 男女边吃奶边做爰视频| 赤兔流量卡办理| 国产一区二区亚洲精品在线观看| 又爽又黄a免费视频| 国产 亚洲一区二区三区 | 久久精品国产亚洲av天美| 婷婷色综合大香蕉| 18禁动态无遮挡网站| 男人舔女人下体高潮全视频| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 精品久久久久久成人av| 日韩人妻高清精品专区| 男女下面进入的视频免费午夜| 建设人人有责人人尽责人人享有的 | 亚洲国产精品sss在线观看| 天堂俺去俺来也www色官网 | 亚洲av一区综合| 国产色婷婷99| 成人国产麻豆网| 国产精品伦人一区二区| 亚洲性久久影院| 中文欧美无线码| 在线免费观看不下载黄p国产| 日韩欧美 国产精品| 国产人妻一区二区三区在| 日韩欧美 国产精品| 国产精品久久久久久精品电影小说 | 欧美精品国产亚洲| 可以在线观看毛片的网站| a级毛片免费高清观看在线播放| 国产午夜精品一二区理论片| 精品一区在线观看国产| 亚洲国产日韩欧美精品在线观看| 在线观看美女被高潮喷水网站| 久久久久久久久久成人| 亚洲精品久久午夜乱码| 精品酒店卫生间| av国产久精品久网站免费入址| 熟女电影av网| 国产精品嫩草影院av在线观看| 一区二区三区高清视频在线| 简卡轻食公司| 神马国产精品三级电影在线观看| 国产69精品久久久久777片| 美女主播在线视频| 狠狠精品人妻久久久久久综合| 亚洲怡红院男人天堂| av在线老鸭窝| 九色成人免费人妻av| 美女被艹到高潮喷水动态| 狂野欧美激情性xxxx在线观看| 免费av不卡在线播放| 99热全是精品| 国产精品国产三级国产av玫瑰| 亚洲,欧美,日韩| 国内精品宾馆在线| eeuss影院久久| 国产激情偷乱视频一区二区| 美女内射精品一级片tv| 国产精品国产三级专区第一集| 久久久久久久久久久丰满| 国产高清国产精品国产三级 | 色综合色国产| 97精品久久久久久久久久精品| 久久精品国产自在天天线| 天堂av国产一区二区熟女人妻| 美女内射精品一级片tv| 只有这里有精品99| 成人美女网站在线观看视频| 最近中文字幕高清免费大全6| 国产女主播在线喷水免费视频网站 | ponron亚洲| av线在线观看网站| 噜噜噜噜噜久久久久久91| 你懂的网址亚洲精品在线观看| 国产精品伦人一区二区| 一级a做视频免费观看| .国产精品久久| 色视频www国产| 国产成人精品婷婷| 一个人看的www免费观看视频| 人妻系列 视频| 国产一区亚洲一区在线观看| 国产毛片a区久久久久| 免费黄色在线免费观看| 最近最新中文字幕免费大全7| 99re6热这里在线精品视频| 亚洲国产色片| 亚洲欧洲国产日韩| 天堂网av新在线| 噜噜噜噜噜久久久久久91| 国产大屁股一区二区在线视频| 老司机影院毛片| 春色校园在线视频观看| 国产一级毛片七仙女欲春2| 简卡轻食公司| 男人舔奶头视频| 特大巨黑吊av在线直播| 我的女老师完整版在线观看| 国产精品人妻久久久久久| 天天躁日日操中文字幕| 国产一级毛片七仙女欲春2| 久99久视频精品免费| 中文字幕av成人在线电影| 麻豆精品久久久久久蜜桃| 亚洲精品日韩av片在线观看| 伊人久久国产一区二区| 天堂√8在线中文| 日韩成人伦理影院| 久久韩国三级中文字幕| 国产午夜福利久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 国产 一区 欧美 日韩| 成人特级av手机在线观看| 九草在线视频观看| 禁无遮挡网站| 午夜福利在线在线| 免费黄网站久久成人精品|