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

    落葉松枯梢病在中國的適生性

    2024-05-06 03:20:12張秀蕓伍文慧梁英梅
    生態(tài)學(xué)報(bào) 2024年7期
    關(guān)鍵詞:環(huán)境變量適生區(qū)落葉松

    張秀蕓,伍文慧,梁英梅

    1北京林業(yè)大學(xué)省部共建森林培育與保護(hù)教育部重點(diǎn)實(shí)驗(yàn)室,北京 100083 2 北京林業(yè)大學(xué)博物館,北京 100083

    落葉松喜光、耐寒、適應(yīng)性強(qiáng)且生長迅速,自然分布于中國東北、華北、內(nèi)蒙古及西南林區(qū),是大面積人工造林的主要樹種。由落葉松新殼梭孢菌(Neofusicoccumlaricinum(Sawada) Hattori &Nakash)引致的落葉松枯梢病[1]主要為害落葉松幼林,引起新梢枯死,連年受害造成樹冠禿頂,對中國落葉松人工林的營造構(gòu)成嚴(yán)重的威脅。落葉松枯梢病于1938年在日本北海道發(fā)現(xiàn),此后在朝鮮、菲律賓、英國、加拿大等國有不同程度的分布和危害[2]。在我國,該病1973年在吉林省首次被發(fā)現(xiàn),其后相繼在東北地區(qū)蔓延并造成巨大經(jīng)濟(jì)損失,現(xiàn)已在我國多個(gè)省市區(qū)有所報(bào)道。該病是一種危險(xiǎn)性傳染病,一旦入侵落葉松林將給樹木生長帶來毀滅性的損失[3]。為此,該病被列入農(nóng)村農(nóng)業(yè)部等六部門發(fā)布的“重點(diǎn)管理外來入侵物種名錄”中的五種植物病害之一。落葉松枯梢病的病原菌在林間主要通過氣流傳播,長距離跨區(qū)域傳播主要靠苗木調(diào)運(yùn),而病害的發(fā)生與林地氣候環(huán)境有密切關(guān)系。因此,基于中國落葉松林的分布特點(diǎn),結(jié)合氣候環(huán)境變化,分析該病害的擴(kuò)散風(fēng)險(xiǎn)及其分布格局,對控制該病害的擴(kuò)散蔓延具有指導(dǎo)意義。

    全球氣溫在過去的100年里上升了大約0.6℃[4],根據(jù)第六次國際耦合模式比較計(jì)劃 (CMIP6)到21世紀(jì)末全球氣候變化還會不斷加劇。氣候變化主要在大尺度上決定物種分布并使得物種分布范圍發(fā)生變化,如隨著氣候變暖許多陸生生物的分布范圍向高海拔和高緯度轉(zhuǎn)移[5—6]。氣候變化也可以通過影響植物病原及其寄主植物從而改變病害的空間分布、發(fā)病率和嚴(yán)重程度。在過去30年里,有超過100篇的文獻(xiàn)及綜述關(guān)注到氣候變化對植物病原和植物疾病風(fēng)險(xiǎn)的影響[7]。氣候變化對不同植物病害的適生范圍變化的影響存在差異,如引起植物病害適生范圍縮小,增加或轉(zhuǎn)移等[8]。而落葉松枯梢病如何響應(yīng)未來氣候變化尚不明確,因此,本研究通過構(gòu)建生態(tài)位模型對其進(jìn)行合理預(yù)測。

    生態(tài)位模型(ENMs)被廣泛應(yīng)用于基于物種分布和環(huán)境變量的特定機(jī)器算法來評估物種的生態(tài)位[9],從而確定物種潛在分布區(qū)。目前常用的生態(tài)位模型包括GARP、BIOCLIM、DOMAIN及MaxEnt[10]。其中最大熵模型(MaxEnt)[11]是基于最大熵理論的統(tǒng)計(jì)模型[12],具有所需樣本量要求不嚴(yán)苛、預(yù)測結(jié)果具有較高精度和可信度的優(yōu)點(diǎn),被廣泛應(yīng)用于入侵生物風(fēng)險(xiǎn)預(yù)測評估[13]、珍稀動植物生境適宜性評價(jià)[14—15]以及植物病蟲害風(fēng)險(xiǎn)預(yù)測等方面。如通過預(yù)測柑橘黃龍病[16](CandidatusLiberibacter asiaticus andCandidatusLiberibacter africanus)、柑橘輪斑病[17](Pseudofabraeacitricarpa)及蘋果腐爛病[18](Valsamali)的潛在分布區(qū),為柑橘和蘋果的生產(chǎn)提供重要保障。此外,通過MaxEnt分析松材線蟲病[19](Bursaphelenchusxylophilus)和松針紅斑病[20](Dothistromapini)的發(fā)生風(fēng)險(xiǎn),為及時(shí)預(yù)防和管理松樹病害提供了指導(dǎo)。落葉松枯梢病作為國家重點(diǎn)管理的兩種林木病害之一,對其在中國的流行及適生性分析還很缺乏,因此對落葉松枯梢病在中國的潛在適生區(qū)進(jìn)行預(yù)測和分析顯得格外重要和迫切。在本研究中,基于落葉松枯梢病分布數(shù)據(jù)和環(huán)境數(shù)據(jù),利用ENMTools、R、MaxEnt、ArcGIS、SDMTools等技術(shù),確定影響落葉松枯梢病分布的的重要環(huán)境因子,明確當(dāng)前氣候下落葉松枯梢病在中國的潛在風(fēng)險(xiǎn)區(qū)、氣候變化對落葉松枯梢病分布區(qū)的影響及其擴(kuò)散風(fēng)險(xiǎn),以期為落葉松枯梢病未分布區(qū)的檢疫和早期預(yù)警提供理論依據(jù),以及時(shí)阻止病害擴(kuò)散蔓延。

    1 材料和方法

    1.1 分布點(diǎn)數(shù)據(jù)

    落葉松枯梢病的地理分布數(shù)據(jù)來自課題組于2021—2022年間在黑龍江塔河縣、內(nèi)蒙古阿里河林場、河北省塞罕壩、遼寧省海城市、遼寧省撫順縣、吉林省敦化市、吉林省延吉市、寧夏六盤山、陜西寧陜縣等地對落葉松枯梢病發(fā)生情況開展實(shí)地調(diào)查的基礎(chǔ)上,結(jié)合相關(guān)文獻(xiàn)和資料中記錄的分布信息以及相關(guān)林業(yè)部門關(guān)于落葉松枯梢病的分布數(shù)據(jù)獲得。為避免由于病害分布點(diǎn)聚集而導(dǎo)致的模型過度擬合,利用ENMTools[21]依據(jù)氣候圖層分辨率篩選分布點(diǎn)數(shù)據(jù)[14],得到用于模型模擬的落葉松枯梢病地理分布點(diǎn)63個(gè)。最后將數(shù)據(jù)在Excel中按照物種名稱、經(jīng)度、緯度整理,以csv格式保存?zhèn)溆谩F渲?東經(jīng)和北緯為正,西經(jīng)和南緯為負(fù)[22]。

    1.2 環(huán)境變量數(shù)據(jù)

    與落葉松枯梢病發(fā)生流行具有密切聯(lián)系的環(huán)境因素主要是溫度和降水,而19個(gè)生物氣候變量在描述年度趨勢和季節(jié)性的溫度和降水時(shí)具有重要生物學(xué)意義。因此環(huán)境變量選擇來自WorldClim數(shù)據(jù)庫(https://www.worldclim.org/)中的19個(gè)氣候因子。未來數(shù)據(jù)選擇第六次國際耦合模式比較計(jì)劃(CMIP6)中BCC-CSM2-MR模型的預(yù)測數(shù)據(jù)[14],選用新排放情景共享社會經(jīng)濟(jì)路徑(SSPs)中的3種情景,即綠色發(fā)展路徑(ssp126)、中等發(fā)展路徑(ssp245)、高速發(fā)展路徑(ssp585)情景進(jìn)行分析模擬[21],探究在未來不同發(fā)展路徑下落葉松枯梢病在中國的潛在地理分布。數(shù)據(jù)空間分辨率為5arc-minutes。

    將數(shù)據(jù)利用ArcGIS進(jìn)行掩膜提取、格式轉(zhuǎn)化(ASCII)和采樣[23]等處理,用于模型模擬和變量間相關(guān)性大小的計(jì)算。為避免環(huán)境數(shù)據(jù)間過度擬合,對19個(gè)環(huán)境變量進(jìn)行篩選。首先利用63個(gè)落葉松枯梢病分布點(diǎn)和19個(gè)環(huán)境變量在MaxEnt模型中運(yùn)行一次,利用刀切法得到19個(gè)環(huán)境變量的貢獻(xiàn)率,然后用R計(jì)算各環(huán)境變量間的相關(guān)性大小。當(dāng)變量間的相關(guān)系數(shù)絕對值大于0.85時(shí),去除貢獻(xiàn)率小的變量[24—25],最終篩選出10個(gè)環(huán)境變量用于模型模擬。

    1.3 最大熵模型構(gòu)建及精度檢驗(yàn)

    將落葉松枯梢病分布點(diǎn)數(shù)據(jù)(csv)和環(huán)境變量(asc)加載到MaxEnt中分析落葉松枯梢病在中國的潛在適生區(qū)。模型參數(shù)設(shè)置如下:隨機(jī)選取75%和25%的物種分布數(shù)據(jù)作為訓(xùn)練集(Training data)和測試集(Testing data),輸出格式(Output format)選擇Logistic,重復(fù)訓(xùn)練次數(shù)(Replicates)10次,重復(fù)迭代方式(Replicated run type)選擇Subsample。用刀切法計(jì)算環(huán)境變量對預(yù)測模型的貢獻(xiàn)率,運(yùn)行創(chuàng)建響應(yīng)曲線(Create response curves)。此外,MaxEnt模型的準(zhǔn)確性和預(yù)測精度受調(diào)控倍頻和特征組合兩個(gè)參數(shù)的影響較大[26],本研究通過R包kuenm優(yōu)化參數(shù)避免模型過度擬合[27]。在此測試了MaxEnt模型中40個(gè)調(diào)控倍頻[0—4]和29種特征組合[L,Q,P,T,H,LQ,LP,LT,LH,QP,QT,QH,PT,PH,TH,LQP,LQT,LQH,LPT,LPH,QPT,QPH、QTH,PTH,LQPT,LQPH,LQTH,LPTH,LQPTH]等1160種參數(shù)組合,最后選擇delta AICc(The minimum information criterion AICc value. delta. AICc)最小值對應(yīng)的參數(shù)組合建立MaxEnt模型[28—29]。

    利用受試者工作特征曲線下面積(AUC)評價(jià)預(yù)測結(jié)果的精確度,當(dāng)AUC值為0.5—0.6時(shí),模型預(yù)測視為失敗;AUC值為0.6—0.7時(shí),模型預(yù)測視為較差;AUC值為0.7—0.8時(shí),模型預(yù)測視為一般;AUC值為0.8—0.9時(shí),模型預(yù)測視為好;AUC值大于0.9時(shí),模型預(yù)測視為非常好[30]。

    1.4 ArcGIS處理

    將MaxEnt模型運(yùn)行的結(jié)果導(dǎo)入ArcGIS中,進(jìn)行格式轉(zhuǎn)換和重分類[24]。采用等間距分割法[31]將落葉松枯梢病的適生等級劃分為四類,分別為非適生區(qū)(0—0.2)、低適生區(qū)(0.2—0.4)、中適生區(qū)(0.4—0.6)、高適生區(qū)(0.6—1)。繼而對各氣候模式下的結(jié)果圖層利用適生區(qū)柵格屬性表Attribute Table中的圖形比例換算得到各適生等級的面積占比。使用ArcGIS插件SDMTools中的質(zhì)心變化(Centroid Change)功能來模擬未來適生區(qū)與當(dāng)前適生區(qū)相比的遷移方向與路徑[32]。

    2 結(jié)果與分析

    2.1 模型優(yōu)化參數(shù)與精度評估

    基于篩選出的63個(gè)落葉松枯梢病分布點(diǎn)和10個(gè)環(huán)境變量,運(yùn)用R語言kuenm包篩選MaxEnt模型最優(yōu)參數(shù),結(jié)果如下(表1)。落葉松枯梢病在中國的適生區(qū)預(yù)測模型的最優(yōu)參數(shù)組合為: 組合特征=LQ、調(diào)控倍數(shù)=0.8,在此參數(shù)組合下最小信息準(zhǔn)則AICc最小為0,低于默認(rèn)值。

    表1 MaxEnt模型kuenm優(yōu)化評價(jià)指標(biāo)Table 1 Evaluation metrics of MaxEnt model generated by kuenm

    利用AUC值評估模型精度。對不同年份和不同氣候情景下落葉松枯梢病的潛在分布區(qū)進(jìn)行預(yù)測。結(jié)果表明訓(xùn)練集AUC值和測試集AUC值均大于0.9,顯著高于隨機(jī)模型AUC值(0.5)(表2)。說明該模型預(yù)測結(jié)果有較高的準(zhǔn)確度和可信度,即落葉松枯梢病分布模擬結(jié)果與可能的實(shí)際分布范圍存在高度一致性。

    表2 MaxEnt模型預(yù)測精度AUC值Table 2 AUC of prediction accuracy of MaxEnt model

    2.2 影響落葉松枯梢病分布的主要環(huán)境變量

    依據(jù)模型運(yùn)行結(jié)果(表3),對落葉松枯梢病分布貢獻(xiàn)率較大的環(huán)境變量為年平均氣溫(bio1)、最暖季度降水量(bio18)、等溫性(bio3)、溫度季節(jié)性變化標(biāo)準(zhǔn)差(bio4)、最干季度平均溫度(bio9),其累計(jì)貢獻(xiàn)率為94.4%;根據(jù)各環(huán)境因子的AUC值(圖1)表明最暖季度平均溫度(bio10)、年降水量(bio12)、年平均氣溫(bio1)、最濕季度降水量(bio16)、等溫性(bio3)對落葉松枯梢病在中國的潛在適生分布預(yù)測影響較大。綜合以上分析表明年平均氣溫(bio1)、最暖季度降水量(bio18)、年降水量(bio12)、最暖季度平均溫度(bio10)在預(yù)測落葉松枯梢病中國潛在地理分布中起著主導(dǎo)作用。

    圖1 MaxEnt模型預(yù)測環(huán)境變量因子的AUCFig.1 AUC value of environmental variable factors predicted by MaxEnt model AUC:受試者工作特征曲線下面積 Area under the receiver operating characteristic curve ;bio1:年平均氣溫Annual Mean Temperature;bio10:最暖季度平均溫度Mean Temperature of Warmest Quarter;bio12:年降水量Annual Precipitation;bio16:最濕季度降水量Precipitation of the wettest quarter;bio18:最暖季度降水量Precipitation of Warmest Quarter;bio19:最冷季度降水量Precipitation of Coldest Quarter;bio3:等溫性Isothermality;bio4:溫度季節(jié)性變化標(biāo)準(zhǔn)差Temperature Seasonality; bio8:最濕季度平均溫度Mean Temperature of Wettest Quarter;bio9:最干季度平均溫度Mean Temperature of Driest Quarter

    表3 主要的的環(huán)境變量貢獻(xiàn)率Table 3 Percent contribution of major environmental variables

    2.3 現(xiàn)代氣候情境下落葉松枯梢病在中國的潛在適生分布區(qū)

    利用ArcGIS將MaxEnt模型預(yù)測出的結(jié)果進(jìn)行可視化處理并統(tǒng)計(jì)各適生等級面積占比(表4;圖2)。結(jié)果表明在當(dāng)前氣候情境下落葉松枯梢病中國的潛在適生區(qū)面積占中國國土總面積的18.02%,主要分布在中國的遼寧、吉林、黑龍江的東部、內(nèi)蒙古的東北部、陜西西部、甘肅南部。此外,在河北、北京、山西、湖北、甘肅、青海、陜西、寧夏、四川、貴州等地也有零星分布。其中高適生區(qū)占3.1%,集中在吉林東部、遼寧東部、黑龍江東南部及中部局部地區(qū);中適生區(qū)占5.5%,主要位于遼寧、吉林的中部、黑龍江中東部、內(nèi)蒙古東部的局部地區(qū);低適生區(qū)占9.4%,主要位于遼寧西部、吉林西部及黑龍江西部、內(nèi)蒙古東北部、山東東部、河北北部、甘肅東南部。

    圖2 現(xiàn)代情境下落葉松枯梢病在中國的潛在適生分布區(qū)Fig.2 Potential suitable area of shoot blight of larch in China under current climate文中地圖來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(https://www.resdc.cn/)

    表4 未來氣候情景下落葉松枯梢病在中國的潛在適生區(qū)面積所占比例/%Table 4 Potential suitable area of shoot blight of larch in China under future scenarios

    2.4 未來氣候條件下落葉松枯梢病的潛在分布區(qū)變化

    未來氣候情境下落葉松枯梢病在中國的潛在適生區(qū)面積和當(dāng)前氣候相比呈減少趨勢(表4),其質(zhì)心有由東北向華北、西南蔓延的趨勢。

    在綠色發(fā)展模式下(圖3),到2021—2040年、2041—2060年、2061—2080年、2081—2100年落葉松枯梢病在中國的潛在適生區(qū)面積占比依次減少到15.56%、15.47%、16.59%、17.70%。相較現(xiàn)代氣候情景,未來落葉松枯梢病在中國吉林、遼寧東部的適生性下降,部分高適生區(qū)轉(zhuǎn)為中適生區(qū)。在中等程度發(fā)展模式下(圖4),到2021—2040年、2041—2060年、2061—2080年、2081—2100年落葉松枯梢病在中國的潛在適生區(qū)面積占比減少到16.01%、14.71%、14.66%、13.80%。其中高適生區(qū)仍集中在中國東北地區(qū),但其適生性整體下降。在吉林、遼寧西部、黑龍江西北部及內(nèi)蒙古東部部分地區(qū)的適生性下降,部分地區(qū)由低適生區(qū)轉(zhuǎn)為非適生區(qū)。在高化石燃料消耗的快速發(fā)展模式下(圖5),到2021—2040年、2041—2060年、2061—2080年、2081—2100年落葉松枯梢病在中國的潛在適生區(qū)面積占比下降到18.10%、15.12%、11.66%、10.87%。該模式下,落葉松枯梢病在中國東北地區(qū)的適生性較中等程度發(fā)展模式進(jìn)一步下降,在中國西南、西北局部地區(qū)適生性增加,部分非適生區(qū)轉(zhuǎn)為適生區(qū)。

    圖3 未來氣候情景(ssp126)下落葉松枯梢病在中國的潛在適生分布區(qū)Fig.3 Potential suitable area of shoot blight of larch in China under future scenario (ssp126)

    圖4 未來氣候情景(ssp245)下落葉松枯梢病在中國的潛在適生分布區(qū)Fig.4 Potential suitable area of shoot blight of larch in China under future scenario (ssp245)

    圖5 未來氣候情景(ssp585)下落葉松枯梢病在中國的潛在適生分布區(qū)Fig.5 Potential suitable area of shoot blight of larch in China under future scenario (ssp585)

    綜合幾種不同的氣候情景來看,到21世紀(jì)末落葉松枯梢病在中國的潛在分布區(qū)面積均會減少。但值得注意的是三種氣候變化情境下其面積變化趨勢不同,在綠色發(fā)展模式下,落葉松枯梢病的在中國的潛在適生區(qū)相較當(dāng)前先大幅度減少,后逐漸增加。此外,在該模式下落葉松枯梢病在中國的潛在適生區(qū)面積變幅最大,說明在全球氣候變化的背景下,落葉松枯梢病對以綠色發(fā)展為主的濃度發(fā)展路徑更為敏感。在中等程度發(fā)展模式和高化石燃料消耗的快速發(fā)展模式下,相較當(dāng)前氣候落葉松枯梢病在中國的潛在適生區(qū)面積一直呈減少態(tài)勢。根據(jù)落葉松枯梢病在中國潛在適生區(qū)的質(zhì)心變化發(fā)現(xiàn)(圖6),在較長時(shí)間段內(nèi)落葉松枯梢病的適生范圍將由中國東北向華北、西南地區(qū)擴(kuò)散蔓延。

    圖6 未來氣候模式下落葉松枯梢病的適生區(qū)質(zhì)心變化趨勢Fig.6 Centroid migration and change of shoot blight of larch suitable area under future scenarios

    3 討論

    3.1 影響落葉松枯梢病分布的環(huán)境變量

    在篩選出的對落葉松枯梢病分布影響較大的10個(gè)環(huán)境變量中,其中與溫度相關(guān)的環(huán)境變量的貢獻(xiàn)率達(dá)到67.6%,與降水相關(guān)的環(huán)境變量的貢獻(xiàn)率為58.3%,最近的一項(xiàng)研究也表明溫度相較降水對真菌病害的分布有更重要的影響[33]。本研究表明影響落葉松枯梢病在中國分布的主導(dǎo)環(huán)境變量包括年平均氣溫、最暖季度降水量、年降水量、最暖季度平均溫度。由此可看出該病對最暖季度的溫濕條件響應(yīng)較為顯著,而此時(shí)的落葉松新梢也由于木質(zhì)化程度低更容易受病原菌侵染而導(dǎo)致病害發(fā)生,這與實(shí)際情況中落葉松枯梢病在中國的中高風(fēng)險(xiǎn)區(qū)的7、8月份為發(fā)病高峰相吻合。此外,本研究通過環(huán)境因子響應(yīng)曲線發(fā)現(xiàn)適宜落葉松枯梢病發(fā)生的最適最暖季度平均溫度范圍為10.1—24.0℃,這接近先前研究適宜落葉松枯梢病菌孢子的最適萌發(fā)溫度范圍15—27℃[34],說明環(huán)境變量的最佳取值范圍與此前的研究結(jié)果基本一致。除溫度外,降水在落葉松枯梢病傳播擴(kuò)散中發(fā)揮著重要的作用,落葉松枯梢病菌的子囊座和分生孢子器在降雨后吸水釋放孢子,后經(jīng)雨水淋洗、反濺而擴(kuò)大病害侵染范圍。因此建議在最暖季度尤其是大量降雨前噴灑藥劑對該病害進(jìn)行防治。

    3.2 氣候變化下落葉松枯梢病在中國的潛在適生區(qū)變化

    非生物、生物因素和物種遷移能力在長期進(jìn)化過程中影響了物種分布,不同的歷史時(shí)期物種分布范圍不同,氣候變化也會導(dǎo)致物種適宜區(qū)的位置發(fā)生變化。對于植物病害而言,氣候變化會使植物病害向高緯度轉(zhuǎn)移[35],但在本研究中落葉松枯梢病的適生范圍并沒有表現(xiàn)出向高緯度轉(zhuǎn)移的趨勢。這可能是由于主要環(huán)境變量的值沒有超過影響落葉松枯梢病適宜發(fā)生的閾值。

    研究采用新的排放情景-共享社會經(jīng)濟(jì)路徑(SSPs)中的3種情景預(yù)測未來落葉松枯梢病在中國的潛在分布及轉(zhuǎn)移趨勢,發(fā)現(xiàn)落葉松枯梢病在未來三種氣候模式下的響應(yīng)結(jié)果略有差異。就其適生區(qū)面積變化而言,在綠色發(fā)展模式下,落葉松枯梢病的潛在適生面積先大幅度減少后逐漸增加,到本世紀(jì)末其適生面積占中國國土總面積的17.70%;而在中度發(fā)展模式和高化石燃料消耗的快速發(fā)展模式下,其適生區(qū)面積均持續(xù)減少,到本世紀(jì)末其適生區(qū)面積減少到中國國土總面積的13.80%和10.78%。此外,相較綠色發(fā)展模式,在中度發(fā)展模式和高化石燃料消耗的快速發(fā)展模式下,落葉松枯梢病的高適生區(qū)呈明顯的破碎化,表明病害有集中爆發(fā)的趨勢。不同氣候模式下適生區(qū)面積變化的差異可能是由于綠色發(fā)展模式下二氧化碳的排放量較少,使溫度上升幅度較低造成的。比較未來三種氣候模式下落葉松枯梢病的質(zhì)心轉(zhuǎn)移軌跡發(fā)現(xiàn),質(zhì)心均位于內(nèi)蒙古自治區(qū)內(nèi),但綠色發(fā)展模式下其轉(zhuǎn)移趨勢偏向華北地區(qū),而中度發(fā)展模式和高化石燃料消耗的快速發(fā)展模式下其軌跡更明顯向西南地區(qū)轉(zhuǎn)移。整體上質(zhì)心轉(zhuǎn)移軌跡是由東北地區(qū)轉(zhuǎn)向華北、西南地區(qū),這可能與中國不同地區(qū)溫室氣體的排放程度不同使年均溫上升存在差異有關(guān),其中在北方地區(qū)特別是東北地區(qū),年均溫上升幅度高于南方地區(qū)[36]。因此在向華北及西南地區(qū)調(diào)運(yùn)落葉松苗木時(shí),有關(guān)檢疫部門應(yīng)加強(qiáng)地區(qū)檢疫預(yù)警工作。

    3.3 研究的局限性

    在滿足病原和寄主同時(shí)存在的情況下,適宜的溫、濕度是病害流行的必要條件。在前人研究的基礎(chǔ)上,本研究基于落葉松枯梢病分布數(shù)據(jù)和氣候數(shù)據(jù),利用MaxEnt模型對落葉松枯梢病在中國的潛在適生區(qū)進(jìn)行了預(yù)測。而在實(shí)際環(huán)境中,落葉松枯梢病的分布和適宜程度也會受到其他復(fù)雜因素的影響,如風(fēng)力、人為干擾、海拔、樹齡、寄主品種及物種間相互作用[37]等。落葉松枯梢病菌的孢子通過風(fēng)力傳播且落葉松苗木在運(yùn)輸過程中不確定的人為干擾會影響病害擴(kuò)散傳播。由于這些因素對病害的影響機(jī)制不明確或難以量化而影響模型模擬效果,因此該預(yù)測結(jié)果可作為判斷落葉松枯梢病與環(huán)境變量之間關(guān)系的參考,但不能完全概括二者之間的關(guān)系。此外,研究結(jié)果旨在預(yù)測落葉松枯梢病潛在適生分布區(qū)并非其實(shí)際分布,故即使生態(tài)位模型預(yù)測該物種分布可能性很高,也不代表其在此區(qū)域有實(shí)際分布[38]。

    寄主的存在是研究病害適生區(qū)的前提,該研究對落葉松枯梢病及其寄主在中國的分布情況進(jìn)行調(diào)查,結(jié)果顯示落葉松枯梢病的分布與預(yù)測結(jié)果較吻合,且落葉松枯梢病的寄主分布涵蓋了預(yù)測結(jié)果中落葉松枯梢病的中高適生分布區(qū)(圖7),說明了MaxEnt模型在林木病害預(yù)測方面的可行性。落葉松枯梢病主要為害的落葉松屬植物種類在中國分布較為廣泛,有利于病害在中國的傳播擴(kuò)散。而落葉松枯梢病長距離傳播主要通過調(diào)運(yùn)感病落葉松苗木,因此本研究為落葉松枯梢病的監(jiān)測預(yù)警提供了參考依據(jù)。此外,落葉松枯梢病發(fā)展蔓延迅速,其林內(nèi)擴(kuò)散途徑主要包括孢子借風(fēng)力雨水?dāng)U大侵染范圍而加重病情、菌絲在病株體內(nèi)擴(kuò)散蔓延引起部分新稍和老枝發(fā)病,加之林木病害相較作物病害其寄主具多年生、受人為干擾小的特性,因此該病一旦發(fā)生便很難控制,這使得預(yù)測結(jié)果更具有現(xiàn)實(shí)意義。

    圖7 落葉松枯梢病寄主植物在中國的分布Fig.7 Distribution of the main hosts of shoot blight of larch in China

    4 結(jié)論

    本研究利用MaxEnt生態(tài)位模型預(yù)測了落葉松枯梢病在中國當(dāng)前和未來氣候下的潛在適生分布區(qū)并明確了影響落葉松枯梢病發(fā)生的主要環(huán)境變量。根據(jù)分析結(jié)果發(fā)現(xiàn)年平均氣溫、最暖季度降水量、年降水量、最暖季度平均溫度與落葉松枯梢病的發(fā)生流行有密切的聯(lián)系。當(dāng)前氣候情境下落葉松枯梢病的適生區(qū)主要位于中國東北地區(qū)、西南及西北局部地區(qū)。未來不同氣候情境下落葉松枯梢病在中國的潛在適生面積出現(xiàn)不同程度的縮小,其適生范圍有由東北向華北、西南地區(qū)擴(kuò)散的趨勢。落葉松枯梢病危害嚴(yán)重,其寄主落葉松在我國又分布廣泛,因此本研究通過對落葉松枯梢病在中國的潛在適生區(qū)進(jìn)行預(yù)測和分析,可為林業(yè)相關(guān)檢疫部門對落葉松枯梢病的防控與檢疫工作提供參考依據(jù)。

    猜你喜歡
    環(huán)境變量適生區(qū)落葉松
    基于最大熵模型的云南思茅松潛在分布區(qū)
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預(yù)測
    未來氣候條件下當(dāng)歸適生區(qū)預(yù)測及時(shí)空變化分析
    山西落葉松雜交良種逾10萬畝
    從桌面右鍵菜單調(diào)用環(huán)境變量選項(xiàng)
    氣候變化下瀕危植物半日花在中國的潛在分布
    徹底弄懂Windows 10環(huán)境變量
    長白落葉松離體再生體系的建立
    巴拉圭瓜多竹適生區(qū)分布研究
    神奇的落葉松提取物
    国产精品一区二区在线不卡| 日韩av在线免费看完整版不卡| 中文乱码字字幕精品一区二区三区| 嫩草影院入口| 欧美日韩视频高清一区二区三区二| 精品久久久久久久久av| 午夜免费鲁丝| 观看美女的网站| 午夜91福利影院| 老司机影院毛片| 成年女人在线观看亚洲视频| 日日啪夜夜爽| av在线app专区| 丝袜美足系列| 日本av手机在线免费观看| 午夜福利网站1000一区二区三区| 最近中文字幕2019免费版| 国产在线视频一区二区| 免费观看av网站的网址| 母亲3免费完整高清在线观看 | 久久av网站| 亚洲国产精品一区三区| 国产精品欧美亚洲77777| 天堂8中文在线网| 成人无遮挡网站| a级毛片黄视频| 如何舔出高潮| 国产色婷婷99| 国产色婷婷99| 岛国毛片在线播放| 国产成人一区二区在线| 国产成人午夜福利电影在线观看| 国产一级毛片在线| 九九爱精品视频在线观看| 丝袜在线中文字幕| videos熟女内射| 秋霞伦理黄片| 伦精品一区二区三区| videosex国产| 在线天堂最新版资源| 成人手机av| 亚洲伊人久久精品综合| 大香蕉久久成人网| 街头女战士在线观看网站| 亚洲av中文av极速乱| 欧美精品一区二区大全| 街头女战士在线观看网站| 成人国产av品久久久| 简卡轻食公司| 2018国产大陆天天弄谢| 国产精品一国产av| 日韩人妻高清精品专区| 高清黄色对白视频在线免费看| 少妇的逼水好多| 中文字幕免费在线视频6| 狂野欧美激情性bbbbbb| 国内精品宾馆在线| 亚洲成人手机| 欧美日韩视频高清一区二区三区二| 日本爱情动作片www.在线观看| 亚洲精品色激情综合| 久久精品国产自在天天线| 国产精品无大码| 男女啪啪激烈高潮av片| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久成人| 伦精品一区二区三区| 91精品三级在线观看| av在线app专区| 亚洲丝袜综合中文字幕| videos熟女内射| 天堂俺去俺来也www色官网| 亚洲四区av| 最近的中文字幕免费完整| 少妇熟女欧美另类| 欧美日韩成人在线一区二区| 纵有疾风起免费观看全集完整版| 在线观看人妻少妇| 美女中出高潮动态图| 三级国产精品片| 综合色丁香网| 天天躁夜夜躁狠狠久久av| av国产精品久久久久影院| 校园人妻丝袜中文字幕| av国产久精品久网站免费入址| 国产成人精品在线电影| 国产精品嫩草影院av在线观看| 精品人妻熟女毛片av久久网站| 在线 av 中文字幕| 色吧在线观看| 亚洲国产av影院在线观看| 免费观看av网站的网址| 国产熟女欧美一区二区| 另类亚洲欧美激情| 免费黄网站久久成人精品| 少妇人妻精品综合一区二区| 久久精品久久精品一区二区三区| 成人国语在线视频| 岛国毛片在线播放| 亚洲久久久国产精品| 免费人成在线观看视频色| 欧美丝袜亚洲另类| 亚洲精品第二区| 精品一区二区三区视频在线| 天堂8中文在线网| 精品国产乱码久久久久久小说| 永久免费av网站大全| 亚洲精品自拍成人| 在线天堂最新版资源| 91精品三级在线观看| 亚洲美女搞黄在线观看| 老司机亚洲免费影院| 99久久综合免费| 中文字幕精品免费在线观看视频 | 99热这里只有是精品在线观看| 老司机亚洲免费影院| a级毛片黄视频| 国产精品99久久99久久久不卡 | 久久久久久久国产电影| 婷婷色av中文字幕| 欧美成人午夜免费资源| 乱人伦中国视频| 久久国产精品男人的天堂亚洲 | av视频免费观看在线观看| 高清在线视频一区二区三区| 午夜福利视频在线观看免费| 欧美精品高潮呻吟av久久| 久久狼人影院| 国产成人午夜福利电影在线观看| 在线观看人妻少妇| 又黄又爽又刺激的免费视频.| av在线老鸭窝| 国产精品麻豆人妻色哟哟久久| 亚洲五月色婷婷综合| 亚洲综合精品二区| 国产片特级美女逼逼视频| 97精品久久久久久久久久精品| 亚洲国产色片| 熟妇人妻不卡中文字幕| 两个人的视频大全免费| 欧美三级亚洲精品| 久久人妻熟女aⅴ| 黑丝袜美女国产一区| 色94色欧美一区二区| 女的被弄到高潮叫床怎么办| 久久精品久久久久久久性| 国产精品久久久久成人av| 日韩欧美精品免费久久| 搡老乐熟女国产| 有码 亚洲区| 国产国语露脸激情在线看| 男女免费视频国产| 老熟女久久久| 在线观看美女被高潮喷水网站| 我的女老师完整版在线观看| 国产精品国产三级专区第一集| 满18在线观看网站| 久久久久久久久久成人| 亚洲av在线观看美女高潮| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产一区二区三区久久久樱花| 欧美成人精品欧美一级黄| 毛片一级片免费看久久久久| 91在线精品国自产拍蜜月| 久久人人爽人人爽人人片va| 97超视频在线观看视频| 极品人妻少妇av视频| 国产精品秋霞免费鲁丝片| 男人操女人黄网站| 午夜日本视频在线| 狠狠婷婷综合久久久久久88av| 在线观看一区二区三区激情| 国产精品久久久久久精品古装| 最近的中文字幕免费完整| 午夜免费鲁丝| 欧美bdsm另类| 熟女电影av网| 女人久久www免费人成看片| 日韩熟女老妇一区二区性免费视频| 亚洲性久久影院| 香蕉精品网在线| 精品99又大又爽又粗少妇毛片| 亚洲av不卡在线观看| 99九九在线精品视频| 成人无遮挡网站| 王馨瑶露胸无遮挡在线观看| 一区二区三区精品91| 亚洲无线观看免费| 成年女人在线观看亚洲视频| 精品人妻熟女av久视频| 欧美 亚洲 国产 日韩一| 色婷婷久久久亚洲欧美| 国产精品嫩草影院av在线观看| 黄色视频在线播放观看不卡| 一级毛片我不卡| 一本大道久久a久久精品| 大片电影免费在线观看免费| 国产一区二区在线观看av| 亚洲一区二区三区欧美精品| 99精国产麻豆久久婷婷| 久久久精品免费免费高清| 精品亚洲成国产av| 久久韩国三级中文字幕| 91午夜精品亚洲一区二区三区| 国产有黄有色有爽视频| videosex国产| 看免费成人av毛片| 久久久国产精品麻豆| 大香蕉97超碰在线| 亚洲欧洲精品一区二区精品久久久 | 一本大道久久a久久精品| 国国产精品蜜臀av免费| 91久久精品国产一区二区三区| 我要看黄色一级片免费的| 午夜激情av网站| 免费日韩欧美在线观看| 校园人妻丝袜中文字幕| 欧美性感艳星| 日本av免费视频播放| 国产av码专区亚洲av| 视频在线观看一区二区三区| 成人免费观看视频高清| 久久精品国产亚洲av涩爱| 欧美另类一区| 能在线免费看毛片的网站| 国产在线一区二区三区精| 夫妻性生交免费视频一级片| 日本色播在线视频| 99热国产这里只有精品6| 91久久精品国产一区二区三区| 亚洲精品国产色婷婷电影| 日日爽夜夜爽网站| 一区二区三区乱码不卡18| 欧美日韩国产mv在线观看视频| 国产综合精华液| 精品一区在线观看国产| 高清不卡的av网站| 嫩草影院入口| 99热6这里只有精品| 欧美 日韩 精品 国产| 丰满饥渴人妻一区二区三| 久久精品国产a三级三级三级| 一个人看视频在线观看www免费| 久久久久久久久久人人人人人人| 午夜福利网站1000一区二区三区| 91aial.com中文字幕在线观看| 日韩人妻高清精品专区| av福利片在线| 日韩av免费高清视频| 久久久久久久精品精品| 在线观看美女被高潮喷水网站| 黑人高潮一二区| 国产成人91sexporn| 国产一区二区在线观看av| a级毛片在线看网站| 最新中文字幕久久久久| 日本-黄色视频高清免费观看| .国产精品久久| 国产精品一区二区三区四区免费观看| 99视频精品全部免费 在线| 国产免费现黄频在线看| 亚洲精品aⅴ在线观看| 亚洲精品,欧美精品| 18禁裸乳无遮挡动漫免费视频| 久久久精品94久久精品| 精品少妇久久久久久888优播| 男男h啪啪无遮挡| 中文乱码字字幕精品一区二区三区| 亚洲五月色婷婷综合| 欧美成人午夜免费资源| 一区二区日韩欧美中文字幕 | 黑人猛操日本美女一级片| 亚洲国产欧美在线一区| 校园人妻丝袜中文字幕| 亚洲精品国产av成人精品| 另类精品久久| 久久久a久久爽久久v久久| 少妇高潮的动态图| av专区在线播放| 亚洲性久久影院| 男的添女的下面高潮视频| 汤姆久久久久久久影院中文字幕| 极品人妻少妇av视频| 夫妻午夜视频| 99久久中文字幕三级久久日本| 国产男女内射视频| 黑人欧美特级aaaaaa片| 毛片一级片免费看久久久久| 精品一区二区三区视频在线| 九九爱精品视频在线观看| 中文字幕免费在线视频6| 校园人妻丝袜中文字幕| 黑人欧美特级aaaaaa片| 国产成人免费观看mmmm| av免费在线看不卡| 最近的中文字幕免费完整| 高清视频免费观看一区二区| 亚洲在久久综合| av视频免费观看在线观看| 午夜免费男女啪啪视频观看| 久久久久久久久久人人人人人人| 欧美bdsm另类| 亚洲精品久久午夜乱码| 黄色怎么调成土黄色| 97在线视频观看| 好男人视频免费观看在线| 97精品久久久久久久久久精品| av在线app专区| 不卡视频在线观看欧美| 精品一区二区三卡| 乱人伦中国视频| 久久精品国产亚洲av天美| 日本欧美视频一区| 亚洲av二区三区四区| 性高湖久久久久久久久免费观看| 国产男女超爽视频在线观看| 午夜日本视频在线| 99热国产这里只有精品6| 一本大道久久a久久精品| 少妇的逼水好多| 少妇被粗大猛烈的视频| 大香蕉97超碰在线| 曰老女人黄片| 国产 精品1| 在线精品无人区一区二区三| 国产欧美日韩综合在线一区二区| 五月天丁香电影| 亚洲性久久影院| 如日韩欧美国产精品一区二区三区 | 欧美性感艳星| 精品国产国语对白av| 91成人精品电影| 亚洲一区二区三区欧美精品| 啦啦啦视频在线资源免费观看| 亚洲图色成人| av网站免费在线观看视频| 国产成人av激情在线播放 | 满18在线观看网站| 久久精品久久精品一区二区三区| 国产日韩欧美在线精品| 国产精品无大码| 久久国内精品自在自线图片| 18+在线观看网站| 七月丁香在线播放| 国产亚洲一区二区精品| 国产高清不卡午夜福利| 黑人猛操日本美女一级片| 午夜精品国产一区二区电影| 国产永久视频网站| 毛片一级片免费看久久久久| 日韩欧美一区视频在线观看| 少妇丰满av| 亚洲成人手机| 97在线人人人人妻| 在线观看一区二区三区激情| 国产精品一国产av| 一二三四中文在线观看免费高清| a级毛色黄片| 婷婷色av中文字幕| 人妻制服诱惑在线中文字幕| 亚洲精华国产精华液的使用体验| 国产淫语在线视频| 日本免费在线观看一区| 热re99久久精品国产66热6| av女优亚洲男人天堂| 少妇 在线观看| 最近2019中文字幕mv第一页| 日韩一区二区三区影片| 亚洲欧美成人精品一区二区| 国产成人一区二区在线| 永久免费av网站大全| xxx大片免费视频| 美女cb高潮喷水在线观看| 在线观看一区二区三区激情| 久久狼人影院| 日本免费在线观看一区| 草草在线视频免费看| 啦啦啦中文免费视频观看日本| 在线免费观看不下载黄p国产| 亚洲av不卡在线观看| kizo精华| 视频在线观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 成人手机av| 久久久久久伊人网av| 国产视频首页在线观看| 99九九在线精品视频| 亚洲婷婷狠狠爱综合网| 亚洲成人一二三区av| 91精品国产九色| 边亲边吃奶的免费视频| 91在线精品国自产拍蜜月| 18+在线观看网站| 女人精品久久久久毛片| 国产精品一二三区在线看| 天堂中文最新版在线下载| 黑人欧美特级aaaaaa片| 麻豆乱淫一区二区| 免费少妇av软件| 亚洲久久久国产精品| 观看av在线不卡| av天堂久久9| 乱人伦中国视频| 高清在线视频一区二区三区| 人人妻人人爽人人添夜夜欢视频| 精品酒店卫生间| 亚洲色图 男人天堂 中文字幕 | 人妻系列 视频| 丝袜脚勾引网站| av线在线观看网站| h视频一区二区三区| 精品视频人人做人人爽| 久久av网站| 国产有黄有色有爽视频| 大话2 男鬼变身卡| 亚洲精品,欧美精品| 国产精品三级大全| 精品少妇内射三级| 国产在线一区二区三区精| 久久久久国产网址| 日本wwww免费看| 中国三级夫妇交换| 一级片'在线观看视频| 亚洲五月色婷婷综合| 亚洲久久久国产精品| 我的老师免费观看完整版| 最新的欧美精品一区二区| 久久精品国产自在天天线| 中文字幕精品免费在线观看视频 | 少妇被粗大猛烈的视频| 亚洲av.av天堂| 国产精品免费大片| 国产 精品1| 亚洲丝袜综合中文字幕| 大香蕉久久网| 在线看a的网站| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产av成人精品| 亚洲精品乱码久久久v下载方式| 蜜桃在线观看..| 日韩在线高清观看一区二区三区| 亚洲成人av在线免费| 久久99蜜桃精品久久| 国产免费一区二区三区四区乱码| 国产伦理片在线播放av一区| 18在线观看网站| 777米奇影视久久| 国产免费视频播放在线视频| 在线观看www视频免费| 亚洲中文av在线| 欧美bdsm另类| 91在线精品国自产拍蜜月| 91精品国产九色| 国产免费现黄频在线看| 777米奇影视久久| av卡一久久| 欧美亚洲日本最大视频资源| 国产又色又爽无遮挡免| 五月伊人婷婷丁香| 亚洲精品av麻豆狂野| 中文乱码字字幕精品一区二区三区| 九九爱精品视频在线观看| 国产精品.久久久| 国产精品 国内视频| 国产亚洲欧美精品永久| 女性被躁到高潮视频| 国产免费又黄又爽又色| 这个男人来自地球电影免费观看 | 激情五月婷婷亚洲| 久久午夜综合久久蜜桃| 午夜福利网站1000一区二区三区| 一级毛片我不卡| 秋霞伦理黄片| 高清黄色对白视频在线免费看| 一区在线观看完整版| 91久久精品国产一区二区成人| 超碰97精品在线观看| 国产毛片在线视频| 国产一区二区在线观看日韩| 久久精品久久精品一区二区三区| 久久精品国产亚洲av天美| av.在线天堂| 久久ye,这里只有精品| 妹子高潮喷水视频| 亚洲经典国产精华液单| 国产一区二区在线观看av| 久久青草综合色| 熟妇人妻不卡中文字幕| 十八禁网站网址无遮挡| 只有这里有精品99| 天美传媒精品一区二区| 中文字幕人妻熟人妻熟丝袜美| 天天躁夜夜躁狠狠久久av| 成人18禁高潮啪啪吃奶动态图 | 久久久久网色| 精品一品国产午夜福利视频| 国产女主播在线喷水免费视频网站| 亚洲欧洲精品一区二区精品久久久 | videos熟女内射| 国产欧美亚洲国产| 久久久久久久久大av| 高清不卡的av网站| 免费黄网站久久成人精品| 欧美日韩综合久久久久久| av网站免费在线观看视频| 久久精品国产自在天天线| 国产精品久久久久久av不卡| 亚洲美女搞黄在线观看| 交换朋友夫妻互换小说| 国产成人av激情在线播放 | xxx大片免费视频| 国产无遮挡羞羞视频在线观看| 国产有黄有色有爽视频| 亚洲色图综合在线观看| 久久鲁丝午夜福利片| 一区二区三区四区激情视频| 日日摸夜夜添夜夜添av毛片| 91aial.com中文字幕在线观看| 久久久久网色| 国产成人精品久久久久久| 一本大道久久a久久精品| 亚洲av二区三区四区| 日日摸夜夜添夜夜爱| 看十八女毛片水多多多| 在线观看美女被高潮喷水网站| 免费人成在线观看视频色| 夫妻午夜视频| 特大巨黑吊av在线直播| 午夜福利网站1000一区二区三区| 成人国语在线视频| 最近的中文字幕免费完整| 久久这里有精品视频免费| 2022亚洲国产成人精品| 在线观看免费高清a一片| 日韩一本色道免费dvd| 亚洲内射少妇av| 建设人人有责人人尽责人人享有的| 国产日韩欧美在线精品| 肉色欧美久久久久久久蜜桃| 久久久久视频综合| 国产精品嫩草影院av在线观看| 欧美激情 高清一区二区三区| 国产精品99久久久久久久久| 国产极品粉嫩免费观看在线 | 中文字幕av电影在线播放| 中文字幕精品免费在线观看视频 | 成人午夜精彩视频在线观看| 最新中文字幕久久久久| 亚洲,一卡二卡三卡| 久久精品国产a三级三级三级| 亚洲av中文av极速乱| 亚洲欧美日韩卡通动漫| 99热这里只有精品一区| av视频免费观看在线观看| 18禁动态无遮挡网站| 亚洲精品乱码久久久v下载方式| 成年人午夜在线观看视频| 赤兔流量卡办理| 久久久久久伊人网av| 欧美日韩av久久| 亚洲第一av免费看| 91午夜精品亚洲一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 51国产日韩欧美| 午夜福利,免费看| 成人国产av品久久久| 国产69精品久久久久777片| 国产免费视频播放在线视频| 亚洲高清免费不卡视频| 久久这里有精品视频免费| 免费播放大片免费观看视频在线观看| 最后的刺客免费高清国语| 国产淫语在线视频| 国产永久视频网站| 日本午夜av视频| av在线播放精品| 性色av一级| 国产熟女欧美一区二区| 狂野欧美激情性xxxx在线观看| 狂野欧美激情性bbbbbb| 久久亚洲国产成人精品v| 国产成人一区二区在线| 国产熟女欧美一区二区| 这个男人来自地球电影免费观看 | 国产精品人妻久久久影院| 国产精品一二三区在线看| 99国产综合亚洲精品| 色吧在线观看| 丰满饥渴人妻一区二区三| 精品一区二区三卡| 免费高清在线观看视频在线观看| 91久久精品国产一区二区三区| 久久韩国三级中文字幕| 在线观看一区二区三区激情| 观看av在线不卡| 插阴视频在线观看视频| 国产伦理片在线播放av一区| 日本午夜av视频| 国产精品嫩草影院av在线观看| 91精品一卡2卡3卡4卡| 人人澡人人妻人| 免费大片黄手机在线观看| 91久久精品国产一区二区成人| 久久精品人人爽人人爽视色| 一本久久精品| 91久久精品电影网| 曰老女人黄片| 国产极品天堂在线| 国产精品蜜桃在线观看| 国产国语露脸激情在线看| 日本91视频免费播放| 免费久久久久久久精品成人欧美视频 | 久久久久久久久大av| 日本91视频免费播放| 精品久久国产蜜桃|