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

    氣候變化背景下茶角胸葉甲潛在適生區(qū)預(yù)測(cè)

    2023-11-18 07:29:26楊藝帥楊學(xué)宇王玉生胡秋龍史子涵吉進(jìn)軍廖尹俊譚琳
    關(guān)鍵詞:環(huán)境變量適生區(qū)氣候變化

    楊藝帥,楊學(xué)宇,王玉生,胡秋龍,史子涵,吉進(jìn)軍,廖尹俊,譚琳*

    氣候變化背景下茶角胸葉甲潛在適生區(qū)預(yù)測(cè)

    楊藝帥1,楊學(xué)宇2,3,王玉生1,胡秋龍2,3,史子涵1,吉進(jìn)軍4,廖尹俊5,譚琳1*

    (1.湖南農(nóng)業(yè)大學(xué)植物保護(hù)學(xué)院,湖南 長(zhǎng)沙 410128;2.茶學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,湖南 長(zhǎng)沙 410128;3.湖南農(nóng)業(yè)大學(xué)園藝學(xué)院,湖南 長(zhǎng)沙 410128;4.湖南省茶業(yè)集團(tuán)股份有限公司,湖南 長(zhǎng)沙 410126;5.東北大學(xué)江河建筑學(xué)院,遼寧 沈陽(yáng) 110819)

    基于茶角胸葉甲(Lefèvre)截至2020年底的63個(gè)有效分布點(diǎn)和22個(gè)環(huán)境變量,分析篩選出影響茶角胸葉甲分布的8個(gè)重要環(huán)境變量,利用最大熵模型(MaxEnt)和ArcGIS軟件,對(duì)茶角胸葉甲2021—2100年氣候條件下的適生區(qū)進(jìn)行預(yù)測(cè)。結(jié)果表明:MaxEnt模型的受試者工作曲線面積(AUC)達(dá)0.992;在目前氣候條件下,茶角胸葉甲潛在的適生區(qū)總面積為6.69×105km2,主要位于中國(guó)南部和東南部地區(qū),除已報(bào)道發(fā)生過(guò)的地區(qū)外,還包括湖北、安徽、浙江、江蘇和臺(tái)灣;影響茶角胸葉甲潛在適生區(qū)分布的重要環(huán)境因子主要包括最冷季降水量、坡向和平均氣溫年較差,其中最冷季降水量是影響茶角胸葉甲潛在適生區(qū)分布范圍的最關(guān)鍵環(huán)境變量;在SSP126和SSP245情景下,茶角胸葉甲的總適生區(qū)面積在2061—2080年間達(dá)到峰值后下降;在SSP370情景下,茶角胸葉甲的總適生區(qū)面積在2041—2060年間達(dá)到峰值后下降;在SSP585情景下,茶角胸葉甲總適生區(qū)面積持續(xù)上升??傮w來(lái)看,茶角胸葉甲的適生區(qū)中心主要向西、向北和高緯度等內(nèi)陸地區(qū)轉(zhuǎn)移。

    茶角胸葉甲;最大熵模型;適生區(qū)預(yù)測(cè);氣候變化情景

    茶角胸葉甲(Lefèvre),又名黑足角胸葉甲[1],屬鞘翅目(Coleoptera)肖葉甲科(Eumolpidae),成蟲取食茶樹、油茶的嫩梢和芽葉,幼蟲在土下取食根系[2]。自20世紀(jì)80年代局部成災(zāi)以來(lái),該蟲不僅在南方茶區(qū)為害,且有迅速向北擴(kuò)散的趨勢(shì)[3–4]。以湖南為例,該蟲最初僅在湖南郴州茶區(qū)為害[3],但近年來(lái)在除湘西以外的全省境內(nèi)均有為害的報(bào)道[5]。另因其蟲口密度高[4],造成經(jīng)濟(jì)損失大[6],現(xiàn)已成為茶產(chǎn)業(yè)發(fā)展的隱患。茶角胸葉甲在南方茶區(qū)一年發(fā)生一代,以幼蟲在5~10 cm的表土層中越冬[7]。各地區(qū)的成蟲發(fā)生、為害盛期略有不同,如在廣東和福建以4月下旬至5月中旬為盛[1,8],而在湖南和江西則以5月中旬至6月中旬為盛[1,3]。茶角胸葉甲的發(fā)生與茶園降水量、溫度、相對(duì)濕度、坡向、海拔高度等密切相關(guān)。朱來(lái)佳[9]研究發(fā)現(xiàn),昭平縣茶角胸葉甲為害加重與降水量的增加密切相關(guān)。譚濟(jì)才等[3]研究發(fā)現(xiàn),溫度對(duì)茶角胸葉甲發(fā)生時(shí)間有明顯影響,相對(duì)濕度則影響其化蛹率和羽化率,東南向茶園的發(fā)生量比西北向茶園的大。謝振倫等[8]報(bào)道,山地及海拔較高的茶園茶角胸葉甲發(fā)生較重。由于茶角胸葉甲的3個(gè)蟲態(tài)生活在土壤中,隨著茶樹種植面積的不斷擴(kuò)大,跨地區(qū)調(diào)運(yùn)茶苗頻繁,茶苗帶蟲土調(diào)運(yùn)時(shí)均存在擴(kuò)散的風(fēng)險(xiǎn),加之潛在發(fā)生區(qū)域較少,給該蟲的監(jiān)測(cè)與防治帶來(lái)困難:因此,預(yù)測(cè)茶角胸葉甲的潛在分布區(qū),明確影響其分布的主要因素至關(guān)重要。

    唐繼洪等[10]運(yùn)用2041—2080年的氣候預(yù)測(cè)數(shù)據(jù),通過(guò)最大熵模型(MaxEnt)[11]對(duì)草地螟()的越冬區(qū)進(jìn)行了預(yù)測(cè),認(rèn)為草地螟越冬場(chǎng)所將會(huì)擴(kuò)大和北移。尹輝等[12]結(jié)合人類活動(dòng)干擾,運(yùn)用MaxEnt模型預(yù)測(cè)了2021—2100年駱駝刺(Shap.)的分布格局變化,認(rèn)為駱駝刺在西北干旱區(qū)生存并擴(kuò)張。姜明鑫等[13]利用MaxEnt模型預(yù)測(cè)茶小綠葉蟬(Matsuda)的潛在適生區(qū),結(jié)果表明,茶小綠葉蟬潛在適生區(qū)面積呈擴(kuò)大趨勢(shì)。

    以往有關(guān)茶角胸葉甲的研究多集中于其生物學(xué)特性與防治實(shí)踐[3,14],尚無(wú)其潛在適生區(qū)分布的研究。筆者收集茶角胸葉甲發(fā)生信息和氣候變化資料,采用ArcGIS軟件提取地形數(shù)據(jù)[15],導(dǎo)入MaxEnt模型,用以預(yù)測(cè)目前和2021—2100年氣候變化條件下茶角胸葉甲各級(jí)適生區(qū)面積和中心點(diǎn)的變化,以期為茶角胸葉甲的監(jiān)測(cè)、科學(xué)防治以及跨地域苗木調(diào)運(yùn)管理提供依據(jù)。

    1 數(shù)據(jù)來(lái)源與研究方法

    1.1 數(shù)據(jù)來(lái)源

    綜合統(tǒng)計(jì)已有文獻(xiàn)[8–9, 16],截至2020年底,茶角胸葉甲已在福建、廣東、廣西、貴州、湖南、江西等6省(自治區(qū))37個(gè)縣市(表1)發(fā)生、為害。結(jié)合實(shí)地調(diào)查情況,剔除重復(fù)報(bào)道和無(wú)效數(shù)據(jù)后,最終得到63處茶角胸葉甲的有效分布點(diǎn)信息。

    表1 茶角胸葉甲的分布

    環(huán)境數(shù)據(jù)均來(lái)源于世界氣候數(shù)據(jù)庫(kù)(https:// worldclim.org/),其中2021—2040年、2041—2060年、2061—2080年、2081—2100年氣候數(shù)據(jù)選取的是由國(guó)家氣候中心研究開發(fā)的BCC–CSM2–MR氣候系統(tǒng)模式[17]中的4種共享社會(huì)經(jīng)濟(jì)路徑(SSPs)情景,該情景是在典型濃度路徑的基礎(chǔ)上發(fā)展而來(lái)的,用以反映社會(huì)經(jīng)濟(jì)發(fā)展與氣候情景之間關(guān)聯(lián)的情景框架[18–19],包括SSP126、SSP245、SSP370、SSP585[15]。這4種路徑分別是假定2100年輻射強(qiáng)迫穩(wěn)定在2.6、4.5、7.0、8.5 W/m2的情景[20–21]。

    地形數(shù)據(jù)通過(guò)ArcGIS 10.4.1軟件從世界高程數(shù)據(jù)[15]中提取獲得,22個(gè)環(huán)境數(shù)據(jù)分辨率均為2.5。

    1.2 研究方法

    由于選用的氣候數(shù)據(jù)與地形數(shù)據(jù)的精度均為2.5,每個(gè)2.5×2.5的網(wǎng)格中只保留1個(gè)分布點(diǎn)[22]用于MaxEnt模型的構(gòu)建。

    為避免過(guò)度擬合,對(duì)獲取的環(huán)境變量數(shù)據(jù)進(jìn)行Pearson相關(guān)性分析,當(dāng)2項(xiàng)環(huán)境變量的相關(guān)系數(shù)≥0.8,即表現(xiàn)出較強(qiáng)的相關(guān)性時(shí),舍去其中對(duì)茶角胸葉甲分布貢獻(xiàn)率較小的環(huán)境變量,最終確定8個(gè)環(huán)境變量用于MaxEnt模型構(gòu)建(表2)。

    表2 構(gòu)建MaxEnt模型的8個(gè)環(huán)境變量

    參照馬世炎等[22]和劉學(xué)琴等[23]的方法,將茶角胸葉甲的分布數(shù)據(jù)、環(huán)境數(shù)據(jù)以及不同時(shí)間段的氣候數(shù)據(jù)導(dǎo)入MaxEnt3.4.1模型后,選用刀切法對(duì)權(quán)重進(jìn)行檢驗(yàn),設(shè)置受試者工作特征曲線(ROC),在“setting”中勾選“Random seed”,重復(fù)迭代類型,選取交叉驗(yàn)證的方法,重復(fù)10次,以此預(yù)測(cè)茶角胸葉甲的適生區(qū)分布。

    參照呂汝丹等[24]的方法,將茶甲胸葉甲適生區(qū)的預(yù)測(cè)結(jié)果以“Cloglog”格式和“asc”文件類型導(dǎo)出后,導(dǎo)入ArcGIS 10.4.1,利用重分類命令中的自然間斷點(diǎn)分級(jí)法,將茶角胸葉甲分布的適宜度分為非適生區(qū)[0~0.10),低度適生區(qū)[0.10~0.32)、中度適生區(qū)[0.32~0.62)、高度適生區(qū)[0.62~1.00] 4個(gè)適生等級(jí),得到茶角胸葉甲潛在適生區(qū)分布。

    參照姜明鑫等[13]的方法,采用ROC曲線對(duì)MaxEnt模型預(yù)測(cè)結(jié)果的精度進(jìn)行評(píng)估,用ROC曲線下面積值(AUC)表示。AUC值越大,表示選擇的環(huán)境變量與預(yù)測(cè)的分布模型之間的相關(guān)性越大,與隨機(jī)分布相距越遠(yuǎn),預(yù)測(cè)結(jié)果也就越準(zhǔn)確[20, 25]。0.8≤AUC<0.9時(shí),模型擬合良好;0.9≤AUC<1.0時(shí)模型擬合結(jié)果非常準(zhǔn)確[26]。

    2 結(jié)果與分析

    2.1 茶角胸葉甲的潛在適生區(qū)分布

    將63個(gè)茶角胸葉甲的有效分布點(diǎn)的地理坐標(biāo)導(dǎo)入MaxEnt模型,其ROC曲線(圖1)訓(xùn)練集的AUC值為0.986,測(cè)試集的AUC值為0.990。這說(shuō)明利用MaxEnt模型對(duì)茶角胸葉甲適生區(qū)進(jìn)行預(yù)測(cè)所得到的結(jié)果準(zhǔn)確。

    圖1 茶角胸葉甲MaxEnt模型的ROC曲線

    2.2 環(huán)境因子對(duì)茶角胸葉甲潛在適生區(qū)分布的影響

    利用刀切法對(duì)8個(gè)環(huán)境變量進(jìn)行分析,結(jié)果(圖2)表明:在僅考慮單一環(huán)境變量的情況下,對(duì)茶角胸葉甲分布增益值影響最大的3個(gè)環(huán)境變量大小依次為最冷季降水量、最干季降水量、平均氣溫年較差,說(shuō)明這3個(gè)變量具有對(duì)茶角胸葉甲適生區(qū)模型預(yù)測(cè)最有價(jià)值的信息,即對(duì)茶角胸葉甲的分布影響最大。海拔高度、最熱季平均溫度、最濕季平均溫度增益值均在0.8以上,對(duì)茶角胸葉甲的分布也有較大貢獻(xiàn)。當(dāng)缺少坡向變量時(shí),模型增益值減少最多,表明坡向?qū)τ诓杞切厝~甲的分布具有其他變量所沒(méi)有的信息,對(duì)模型預(yù)測(cè)結(jié)果的精確度具有無(wú)可替代的影響。

    環(huán)境變量與物種存在概率之間的響應(yīng)曲線能夠反映環(huán)境變量與物種之間的關(guān)系。當(dāng)茶角胸葉甲的存在概率大于0.2時(shí),此時(shí)對(duì)應(yīng)的環(huán)境變量值被認(rèn)為是茶角胸葉甲開始發(fā)生的閾值。由圖3可知,當(dāng)平均氣溫年較差為6.98 ℃時(shí),茶角胸葉甲開始發(fā)生;氣溫年較差為8.3 ℃時(shí),最適其生存;當(dāng)氣溫年較差達(dá)9.5 ℃后,不適合茶角胸葉甲的生存。當(dāng)最冷季降水量達(dá)到141.0 mm時(shí),茶角胸葉甲發(fā)生為害;當(dāng)最冷季降水量為196.2 mm時(shí),最適其生存;當(dāng)最冷季降水量超過(guò)349.5 mm時(shí),不適茶角胸葉甲生存。坡度的上升對(duì)茶角胸葉甲的存在概率有負(fù)面影響:茶角胸葉甲在各坡向均能存在,且西北坡(坡向賦值292.5~337.5)最適其生存,其次為北坡(0~22.5,337.5~360)和平地茶園(–1),其存在概率自北(0~22.5)向東北(22.5~67.5)到西北(292.5~337.5)逐步上升,自西北(292.5~337.5)向北(337.5~360)存在概率下降。

    圖2 單一環(huán)境變量下正規(guī)化訓(xùn)練增益值

    圖3 主要環(huán)境變量的響應(yīng)曲線

    2.3 2020年氣候條件下茶角胸葉甲適生區(qū)分布的預(yù)測(cè)

    利用MaxEnt模型模擬茶角胸葉甲在目前氣候條件下的潛在適生區(qū)分布。茶角胸葉甲的適生區(qū)總面積為6.69×105km2,其中高適生區(qū)面積為1.41×105km2,主要分布在江西、湖南東部和南部、廣西東北部、廣東北部、福建中部和西南部茶區(qū),湖北、安徽、浙江有零星分布;中適生區(qū)面積為1.93×105km2,主要分布在江西、浙江、湖南中部和西南部、廣西東北部、廣東北部、福建西部、湖北東部茶區(qū),安徽省內(nèi)有零星分布;低適生區(qū)面積為3.35×105km2。除上述茶區(qū)有分布外,貴州、江蘇和臺(tái)灣茶區(qū)也有分布。

    2.4 氣候變化對(duì)茶角胸葉甲潛在分布的影響

    利用MaxEnt模型對(duì)4個(gè)時(shí)間段(2021—2040年、2041—2060年、2061—2080年、2081—2100年)、4種氣候變化(SSP126、SSP245、SSP370、SSP585)情景下的茶角胸葉甲潛在適生區(qū)進(jìn)行預(yù)測(cè),獲得不同氣候變化條件下茶角胸葉甲的潛在適生區(qū)面積(圖4)。不同年份、不同路徑下的茶角胸葉甲潛在適生區(qū)相對(duì)于目前分布區(qū)發(fā)生了變化。在SSP126和SSP245情景下,2061—2080年茶角胸葉甲的總適生區(qū)面積在達(dá)到峰值后下降;在SSP370情景下,2041—2060年茶角胸葉甲的總適生區(qū)面積在達(dá)到峰值后下降;但在SSP585情景下,茶角胸葉甲總適生區(qū)面積將持續(xù)增加。

    圖4 不同氣候變化情景下茶角胸葉甲的適生區(qū)面積

    2021—2040年、SSP585情景下,茶角胸葉甲總適生區(qū)面積減少至2.18×105km2,是現(xiàn)有總適生區(qū)面積的32.6%,但在2061—2080年、SSP245情景下,總適生區(qū)面積達(dá)到9.059×106km2,是現(xiàn)有總適生區(qū)面積的13.6倍。

    2021—2040年、SSP585情景下,茶角胸葉甲高適生區(qū)面積減少至4.00×102km2,是現(xiàn)有總適生區(qū)面積的0.3%,但在2061—2080年、SSP245情景下,高適生區(qū)面積為2.974×106km2,是現(xiàn)有高適生區(qū)面積的21.1倍。

    2021—2040年、SSP585情景下,茶角胸葉甲中適生區(qū)面積減少至2.0×104km2,是現(xiàn)有中適生區(qū)面積的10.4%,但在2061—2080年、SSP245情景下,中適生區(qū)面積為2.956×106km2,是現(xiàn)有中適生區(qū)面積的15.3倍。

    2021—2040年、SSP585情景下,茶角胸葉甲低適生區(qū)面積減少至1.98×105km2,是現(xiàn)有低適生區(qū)面積的59.1%,但在2061—2080年、SSP245情景下,高適生區(qū)面積為3.129×106km2,是現(xiàn)有低適生區(qū)面積的9.3倍。

    對(duì)比目前氣候條件和不同氣候變化情景下茶角胸葉甲的各適生區(qū)地點(diǎn)(表3)發(fā)現(xiàn),在不同氣候變化情景下,茶角胸葉甲的高度適生區(qū)中心移動(dòng)的趨勢(shì)不一致,但大多向西部、北部地區(qū)移動(dòng),而茶角胸葉甲的中度適生區(qū)和低度適生區(qū)中心多有向西北方向、高緯度地區(qū)等內(nèi)陸移動(dòng)的趨勢(shì)。

    表3 不同氣候變化情景下茶角胸葉甲各適生區(qū)中心

    3 討論

    收集茶角胸葉甲發(fā)生信息,借助MaxEnt模型,對(duì)其目前氣候條件下適生區(qū)和4種氣候變化情景下的潛在適生區(qū)進(jìn)行了預(yù)測(cè)。結(jié)果表明,目前氣候條件下,茶角胸葉甲的潛在適生區(qū)總面積為6.69×105km2,主要分布在江西、湖南、湖北、廣西、廣東、福建、安徽、浙江、江蘇、貴州、臺(tái)灣,其中,湖北、安徽、浙江、江蘇和臺(tái)灣氣候條件適宜茶角胸葉甲生存,且可能定殖擴(kuò)散,因而在苗木調(diào)運(yùn)時(shí)需加強(qiáng)防控,避免引入造成危害。在SSP126、SSP245與SSP370情景下,茶角胸葉甲的總適生區(qū)面積達(dá)到峰值后下降;但在SSP585情景下,茶角胸葉甲總適生區(qū)面積持續(xù)上升。

    選取的22個(gè)環(huán)境變量中,最冷季降水量的貢獻(xiàn)率高達(dá)84.7%。茶角胸葉甲存在概率對(duì)最冷季降水量的響應(yīng)曲線先高后低,表明過(guò)濕或過(guò)干旱都會(huì)影響茶角胸葉甲的發(fā)生,與實(shí)際相印證[27–28]。此外,坡度也影響茶角胸葉甲的分布。平地與丘陵茶園一般比高山茶園發(fā)生多[27],這與茶角胸葉甲存在概率隨著坡度上升而減小一致。各坡向均可能存在茶角胸葉甲,且西北坡、北坡、平地最適其生存。而調(diào)查結(jié)果發(fā)現(xiàn)南坡茶園幼蟲發(fā)生量大,這可能是因?yàn)樵谀P湍M過(guò)程中,單從坡度這一要素來(lái)探討茶角胸葉甲存在概率,而實(shí)際調(diào)查的結(jié)論是從多變的復(fù)雜環(huán)境條件下獲得的。

    盡管選取了包括氣候、地形等共計(jì)22個(gè)環(huán)境變量進(jìn)行模擬,但考慮到茶角胸葉甲長(zhǎng)期居于地下,土壤的理化性質(zhì)、肥力水平、松緊程度甚至土壤的類型都可能對(duì)茶角胸葉甲的分布產(chǎn)生影響:因此,在未來(lái)研究茶角胸葉甲適生區(qū)時(shí),還需綜合考慮茶角胸葉甲的生物學(xué)信息、寄主分布、土壤理化性質(zhì)、茶園農(nóng)事管理、防控措施等更多的生物和非生物因素,使得預(yù)測(cè)結(jié)果與實(shí)際情況更加吻合。已收集了足夠多的茶角胸葉甲發(fā)生信息,但考慮到該蟲相關(guān)的公開報(bào)道次數(shù)及實(shí)地采樣的局限性,可能導(dǎo)致采集的樣點(diǎn)在空間分布上相對(duì)集中,可能會(huì)對(duì)預(yù)測(cè)結(jié)果產(chǎn)生一定影響,甚至導(dǎo)致極端值的出現(xiàn)[29]:因此,在尚未報(bào)道茶角胸葉甲為害的省份仍需加強(qiáng)監(jiān)測(cè)與預(yù)防,以獲得更好的防控效果。

    [1] 汪榮灶.黑足角胸葉甲發(fā)生與危害的初步考查[J].福建茶葉,1989,11(2):40–41.

    [2] 何學(xué)友,蔡守平,吳智才.油茶的一種新害蟲——黑足角胸葉甲[J].中國(guó)森林病蟲,2011,30(3):16–17.

    [3] 譚濟(jì)才,劉貴芳,王德興,等.茶角胸葉甲生物學(xué)特性及防治研究[J].湖南農(nóng)學(xué)院學(xué)報(bào),1986,12(4):51–60.

    [4] 汪榮灶,胡勇.江西黑足角胸葉甲成蟲的盛發(fā)期與防治[J].江西植保,2004,27(4):175.

    [5] 包強(qiáng),李冠華,肖蕾,等.4種植物源農(nóng)藥對(duì)茶角胸葉甲的室內(nèi)毒力測(cè)定[J].廣東茶業(yè),2020,171(3):12–15.

    [6] 林阿祥.茶角胸葉甲的發(fā)生與防治初報(bào)[J].茶葉科學(xué)簡(jiǎn)報(bào),1991,32(4):43–44.

    [7] 謝逸菲,何振,李密,等.茶角胸葉甲幼蟲在茶園土壤中的空間分布特征[J].湖南林業(yè)科技,2017,44(5):55–60.

    [8] 謝振倫,馬智華,朱俠慧.茶角胸葉甲的初步觀察[J].中國(guó)茶葉,1985,7(4):6–7.

    [9] 朱來(lái)佳.2011年昭平縣茶角胸葉甲重發(fā)生原因及對(duì)策[J].農(nóng)業(yè)研究與應(yīng)用,2013(2):64–67.

    [10] 唐繼洪,程云霞,羅禮智,等.基于Maxent模型的不同氣候變化情景下我國(guó)草地螟越冬區(qū)預(yù)測(cè)[J].生態(tài)學(xué)報(bào), 2017,37(14):4852–4863.

    [11] STEVEN J P,MIROSLAV D,Robert E S.Maxent software for modeling species niches and distributions (Version 3.4.1)[EB/OL]. [2022–08–19].http://biodivers ityinformatics. amnh.org/open_source/maxent/.

    [12] 尹輝,田聰,馬倩倩,等.氣候變化和人類活動(dòng)干擾下駱駝刺潛在分布格局變化特征[J].生態(tài)學(xué)報(bào),2022,42(18):7349–7361.

    [13] 姜明鑫,鐘文玉,胡海琴,等.氣候變化背景下茶小綠葉蟬在中國(guó)的潛在適生區(qū)預(yù)測(cè)[J].生態(tài)學(xué)雜志,2022,41(10):2008–2016.

    [14] 李先文,譚濟(jì)才,柏曉勇,等.幾種藥劑對(duì)茶角胸葉甲的室內(nèi)殺蟲活性測(cè)定及田間藥效試驗(yàn)[J].現(xiàn)代農(nóng)藥,2008,7(3):44–47.

    [15] 張華,趙浩翔,徐存剛.氣候變化背景下孑遺植物桫欏在中國(guó)的潛在地理分布[J].生態(tài)學(xué)雜志,2021,40(4):968–979.

    [16] 曹潘榮,劉克斌,劉春燕,等.茶角胸葉甲侵害對(duì)嶺頭單樅茶鮮葉芳香物質(zhì)的影響[J].應(yīng)用生態(tài)學(xué)報(bào),2006,17(11):2098–2101.

    [17] WU T,LU Y,F(xiàn)ANG Y,et al.The Beijing climate center climate system model (BCC-CSM): the main progress from CMIP5 to CMIP6[J].Geoscientific Model Development,2019,12(4):1573–1600.

    [18] 張杰,曹麗格,李修倉(cāng),等.IPCC AR5中社會(huì)經(jīng)濟(jì)新情景(SSPs)研究的最新進(jìn)展[J].氣候變化研究進(jìn)展,2013,9(3):225–228.

    [19] MOSS R H,EDMONDS J A,HIBBARD K A,et al.The next generation of scenarios for climate change research and assessment[J].Nature,2010,463:747–756.

    [20] 姜彤,呂嫣冉,黃金龍,等.CMIP6模式新情景(SSP–RCP)概述及其在淮河流域的應(yīng)用[J].氣象科技進(jìn)展,2020,10(5):102–109.

    [21] 張麗霞,陳曉龍,辛?xí)愿瑁瓹MIP6情景模式比較計(jì)劃(ScenarioMIP)概況與評(píng)述[J].氣候變化研究進(jìn)展,2019,15(5):519–525.

    [22] 馬世炎,于洪春,趙奎軍,等.基于MaxEnt模型的大豆蚜全球潛在地理分布分析[J].昆蟲學(xué)報(bào),2022,65(5):630–637.

    [23] 劉學(xué)琴,白明,賀達(dá)漢,等.基于MaxEnt模型的草原步甲屬潛在地理分布研究[J].生態(tài)學(xué)報(bào),2022,42(10): 4217–4224.

    [24] 呂汝丹,何健,劉慧杰,等.羽葉鐵線蓮的分布區(qū)與生態(tài)位模型分析[J].北京林業(yè)大學(xué)學(xué)報(bào),2019,41(2):70–79.

    [25] 張偉萍,胡云云,李智華,等.氣候變化情景下祁連圓柏在青海省的適宜分布區(qū)預(yù)測(cè)[J].應(yīng)用生態(tài)學(xué)報(bào),2021,32(7):2514–2524.

    [26] SWETS J A.Measuring the accuracy of diagnostic systems[J].Science,1988,240:1285–1293.

    [27] 譚濟(jì)才,劉貴芳,王德興,等.茶角胸葉甲的發(fā)生與防治[J].中國(guó)茶葉,1987,9(1):15–16.

    [28] 吳勇,林娟,黎健龍,等.英德市2019年茶角胸葉甲重發(fā)生的氣象成因淺析及其防治試驗(yàn)[J].廣東茶業(yè),2021(2):6–9.

    [29] 王艷君,高泰,石娟.基于MaxEnt模型對(duì)舞毒蛾全球適生區(qū)的預(yù)測(cè)及分析[J].北京林業(yè)大學(xué)學(xué)報(bào),2021,43(9):59–69.

    Prediction of the potential adaptive areas ofunder climate change scenarios

    YANG Yishuai1,YANG Xueyu2,3,WANG Yusheng1,HU Qiulong2,3,SHI Zihan1,JI Jinjun4,LIAO Yinjun5,TAN Lin1*

    (1.College of Plant Protection, Hunan Agricultural University, Changsha, Hunan 410128, China; 2.Key Laboratory of Tea Science of Ministry of Education, Changsha, Hunan 410128, China; 3.College of Horticulture, Hunan Agricultural University, Changsha, Hunan 410128, China; 4.Hunan Tea Group Corporation Limited, Changsha, Hunan 410126, China; 5.JangHo Architecture, Northeastern University, Shenyang, Liaoning 110819, China)

    By the end of 2020, based on 63 effective distribution records of the() in China and 22 environmental variables, eight environmental limiting factors affecting the distribution ofwere determined. The MaxEnt model and ArcGIS software were used to predict the potential adaptive areas of(from 2021 to 2100) under current and future climatic conditions in China. The results show that the receiver operating curve area (AUC) is 0.992 in MaxEnt model. At present climate conditions, the total area of adaptive area is 6.69×105km2, mainly located in the south and southeast of China. Besides the areas that had been infected by, its adaptive region includes Hubei, Anhui, Zhejiang, Jiangsu, and Taiwan provinces. The main environmental factors affecting the potential geographical distribution range ofwere precipitation of the coldest quarter, slope, and annual range of mean temperature. The precipitation of the coldest quarter was the most critical environmental variable affecting the potential adaptive areas range of. In scenarios SSP126 and SSP245, the total adaptive area ofpeaks between 2061 and 2080 and then declines. In scenario SSP370, the total adaptive area ofpeaks between 2041 and 2060 and then declines. However, in scenario SSP585, the total adaptive area ofincreases insistently. Overall, the center of adaptive areas predominantly shifts towards western, northern, and high-latitude inland areas.

    ; maximum entropy model(MaxEnt); adaptive area prediction; climate change scenarios

    Q969.48

    A

    1007–1032(2023)05–0581–07

    楊藝帥,楊學(xué)宇,王玉生,胡秋龍,史子涵,吉進(jìn)軍,廖尹俊,譚琳.氣候變化背景下茶角胸葉甲潛在適生區(qū)預(yù)測(cè)[J].湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2023,49(5):581–587.

    YANG Y S,YANG X Y,WANG Y S,HU Q L,SHI Z H,JI J J,LIAO Y J,TAN L.Prediction of the potential adaptive areas ofunder climate change scenarios[J].Journal of Hunan Agricultural University(Natural Sciences),2023,49(5):581–587.

    http://xb.hunau.edu.cn

    2023–03–31

    2023–08–15

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0202006–04)

    楊藝帥(1999—),男,湖南永州人,碩士研究生,主要從事農(nóng)業(yè)昆蟲與害蟲防治研究,519414001@qq.com;*通信作者,譚琳,博士,副教授,主要從事農(nóng)業(yè)昆蟲與害蟲防治研究,hqltanlin@163.com

    10.13331/j.cnki.jhau.2023.05.012

    責(zé)任編輯:羅慧敏

    英文編輯:羅維

    猜你喜歡
    環(huán)境變量適生區(qū)氣候變化
    基于最大熵模型的云南思茅松潛在分布區(qū)
    氣候變化下中國(guó)蒟蒻薯科箭根薯的地理分布格局預(yù)測(cè)
    未來(lái)氣候條件下當(dāng)歸適生區(qū)預(yù)測(cè)及時(shí)空變化分析
    《應(yīng)對(duì)氣候變化報(bào)告(2022)》發(fā)布
    從桌面右鍵菜單調(diào)用環(huán)境變量選項(xiàng)
    氣候變化下瀕危植物半日花在中國(guó)的潛在分布
    徹底弄懂Windows 10環(huán)境變量
    巴拉圭瓜多竹適生區(qū)分布研究
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    国产一级毛片在线| 亚洲av欧美aⅴ国产| 亚洲成色77777| 国产亚洲欧美在线一区二区| 永久免费av网站大全| 国产成人精品无人区| 汤姆久久久久久久影院中文字幕| 两人在一起打扑克的视频| 精品久久蜜臀av无| 免费一级毛片在线播放高清视频 | 老司机影院成人| 男男h啪啪无遮挡| 夜夜骑夜夜射夜夜干| 侵犯人妻中文字幕一二三四区| 丰满饥渴人妻一区二区三| 欧美精品一区二区免费开放| 久9热在线精品视频| 国产精品免费大片| www日本在线高清视频| 在线亚洲精品国产二区图片欧美| 在线精品无人区一区二区三| 久久免费观看电影| 2018国产大陆天天弄谢| 美国免费a级毛片| 蜜桃国产av成人99| 欧美av亚洲av综合av国产av| 国产在线视频一区二区| 超色免费av| e午夜精品久久久久久久| 91老司机精品| 99热全是精品| 欧美97在线视频| 国产精品.久久久| 亚洲精品日韩在线中文字幕| 午夜激情久久久久久久| 婷婷色综合www| 啦啦啦啦在线视频资源| 亚洲视频免费观看视频| 久久天躁狠狠躁夜夜2o2o | 欧美黑人精品巨大| 欧美精品人与动牲交sv欧美| 男女高潮啪啪啪动态图| 一级黄片播放器| 精品视频人人做人人爽| 国产女主播在线喷水免费视频网站| 亚洲国产av影院在线观看| 脱女人内裤的视频| 18禁观看日本| 午夜福利免费观看在线| 精品一区二区三区四区五区乱码 | 久久毛片免费看一区二区三区| 国产亚洲欧美在线一区二区| 天堂中文最新版在线下载| 精品少妇黑人巨大在线播放| 久久人妻熟女aⅴ| 久久人妻熟女aⅴ| 性色av乱码一区二区三区2| 久久天躁狠狠躁夜夜2o2o | 国产老妇伦熟女老妇高清| 天堂俺去俺来也www色官网| 曰老女人黄片| 麻豆乱淫一区二区| 亚洲三区欧美一区| 欧美成人午夜精品| 悠悠久久av| 丝瓜视频免费看黄片| 国产在线观看jvid| 久久亚洲国产成人精品v| 在线精品无人区一区二区三| 叶爱在线成人免费视频播放| 欧美日韩亚洲综合一区二区三区_| 国产精品 国内视频| 大香蕉久久成人网| 亚洲欧洲国产日韩| 国产日韩一区二区三区精品不卡| 99国产精品一区二区三区| 男人舔女人的私密视频| 久久久精品免费免费高清| 女性生殖器流出的白浆| 亚洲精品久久久久久婷婷小说| 两人在一起打扑克的视频| 欧美久久黑人一区二区| 国产片特级美女逼逼视频| 国产精品久久久久久精品电影小说| 成年人午夜在线观看视频| 18禁观看日本| 欧美黑人欧美精品刺激| 成人亚洲欧美一区二区av| 欧美日韩亚洲综合一区二区三区_| 欧美日韩福利视频一区二区| 国产高清videossex| 2021少妇久久久久久久久久久| 777久久人妻少妇嫩草av网站| 亚洲av成人不卡在线观看播放网 | 午夜免费成人在线视频| 亚洲,欧美精品.| 51午夜福利影视在线观看| 日韩大码丰满熟妇| 亚洲伊人久久精品综合| 亚洲国产精品一区三区| 高潮久久久久久久久久久不卡| 伦理电影免费视频| 免费一级毛片在线播放高清视频 | 高潮久久久久久久久久久不卡| 亚洲欧美一区二区三区久久| 97人妻天天添夜夜摸| 亚洲三区欧美一区| 午夜福利视频在线观看免费| 18禁国产床啪视频网站| 亚洲av电影在线进入| 国产精品 国内视频| 国产麻豆69| 大片电影免费在线观看免费| 青青草视频在线视频观看| 国产在视频线精品| 精品国产超薄肉色丝袜足j| 亚洲成人免费av在线播放| 欧美黑人精品巨大| 天天添夜夜摸| 麻豆av在线久日| 国产精品亚洲av一区麻豆| 高清不卡的av网站| 一边摸一边抽搐一进一出视频| 少妇粗大呻吟视频| 99久久人妻综合| 五月天丁香电影| av在线app专区| 在线观看人妻少妇| 国产精品秋霞免费鲁丝片| 国产熟女欧美一区二区| 免费不卡黄色视频| 亚洲精品国产色婷婷电影| 丁香六月欧美| 免费在线观看日本一区| 女人爽到高潮嗷嗷叫在线视频| 老司机在亚洲福利影院| 成年av动漫网址| 蜜桃在线观看..| 日本黄色日本黄色录像| 国产免费现黄频在线看| 交换朋友夫妻互换小说| 精品国产乱码久久久久久小说| √禁漫天堂资源中文www| 亚洲中文日韩欧美视频| 一二三四在线观看免费中文在| 亚洲精品一区蜜桃| 欧美人与善性xxx| 男女免费视频国产| 人人妻人人澡人人看| 免费av中文字幕在线| 亚洲av美国av| 成年人免费黄色播放视频| 国产亚洲精品久久久久5区| 精品熟女少妇八av免费久了| 日本猛色少妇xxxxx猛交久久| 日日爽夜夜爽网站| 欧美日韩精品网址| 好男人电影高清在线观看| 热99久久久久精品小说推荐| 丝袜美腿诱惑在线| 男人操女人黄网站| 夫妻午夜视频| 国语对白做爰xxxⅹ性视频网站| 一二三四在线观看免费中文在| 国产成人系列免费观看| 成年人免费黄色播放视频| 亚洲黑人精品在线| 交换朋友夫妻互换小说| 久久99一区二区三区| 色精品久久人妻99蜜桃| 亚洲国产av影院在线观看| 国产日韩一区二区三区精品不卡| 欧美 亚洲 国产 日韩一| 国产无遮挡羞羞视频在线观看| 乱人伦中国视频| 国产伦人伦偷精品视频| 丰满迷人的少妇在线观看| 亚洲国产成人一精品久久久| 国产精品久久久久久人妻精品电影 | 婷婷丁香在线五月| 每晚都被弄得嗷嗷叫到高潮| 免费观看人在逋| 国产精品 欧美亚洲| 欧美在线一区亚洲| 日本av手机在线免费观看| 黑人巨大精品欧美一区二区蜜桃| 日韩av不卡免费在线播放| 亚洲精品久久午夜乱码| 日韩电影二区| 亚洲欧美成人综合另类久久久| 中文字幕精品免费在线观看视频| 欧美黄色片欧美黄色片| 日韩中文字幕视频在线看片| 婷婷成人精品国产| 精品免费久久久久久久清纯 | 十八禁网站网址无遮挡| 熟女少妇亚洲综合色aaa.| 亚洲成人免费电影在线观看 | 美女国产高潮福利片在线看| 国产三级黄色录像| 国产片内射在线| 欧美 日韩 精品 国产| 精品国产一区二区三区久久久樱花| 国产激情久久老熟女| 热re99久久精品国产66热6| 国产成人精品在线电影| 久久国产精品人妻蜜桃| 国产精品麻豆人妻色哟哟久久| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 亚洲国产成人一精品久久久| 美女高潮到喷水免费观看| 十分钟在线观看高清视频www| 亚洲,一卡二卡三卡| 999久久久国产精品视频| 亚洲一码二码三码区别大吗| 亚洲成人手机| 人妻一区二区av| 国产一区二区在线观看av| 国产黄色视频一区二区在线观看| 成人免费观看视频高清| 9色porny在线观看| 涩涩av久久男人的天堂| 午夜av观看不卡| av电影中文网址| 熟女av电影| 夫妻午夜视频| 女人被躁到高潮嗷嗷叫费观| 亚洲黑人精品在线| 精品卡一卡二卡四卡免费| 国产主播在线观看一区二区 | 国产在线一区二区三区精| 久久久久久久久久久久大奶| 午夜福利视频精品| 人妻一区二区av| 一本—道久久a久久精品蜜桃钙片| 一边摸一边做爽爽视频免费| 久久人妻熟女aⅴ| 亚洲精品久久午夜乱码| 在线观看免费视频网站a站| 日韩伦理黄色片| 9热在线视频观看99| 九草在线视频观看| 国产精品一区二区免费欧美 | 51午夜福利影视在线观看| 人人妻人人添人人爽欧美一区卜| 日韩欧美一区视频在线观看| 亚洲av男天堂| 999久久久国产精品视频| 亚洲国产欧美一区二区综合| 一级毛片我不卡| 另类亚洲欧美激情| 久久性视频一级片| 在线观看免费高清a一片| cao死你这个sao货| 一级毛片黄色毛片免费观看视频| 一区二区av电影网| 久久ye,这里只有精品| 国产一级毛片在线| 中文字幕人妻丝袜一区二区| 大话2 男鬼变身卡| 久久这里只有精品19| 精品少妇内射三级| 精品人妻1区二区| 亚洲九九香蕉| 2018国产大陆天天弄谢| 欧美精品亚洲一区二区| 日韩电影二区| 一本久久精品| 少妇粗大呻吟视频| 色94色欧美一区二区| av片东京热男人的天堂| 超碰成人久久| 久久国产精品男人的天堂亚洲| 日韩av不卡免费在线播放| 久久影院123| 国产亚洲av高清不卡| 尾随美女入室| 国产成人一区二区三区免费视频网站 | 一级毛片电影观看| 国产成人欧美| 久久人妻熟女aⅴ| 国产精品一区二区免费欧美 | 涩涩av久久男人的天堂| 久久九九热精品免费| 女性被躁到高潮视频| videosex国产| 91老司机精品| 黄片小视频在线播放| 亚洲综合色网址| 在线观看www视频免费| 久久久国产精品麻豆| 少妇人妻久久综合中文| √禁漫天堂资源中文www| 黄色视频在线播放观看不卡| 欧美激情 高清一区二区三区| 午夜老司机福利片| 免费日韩欧美在线观看| 久久精品国产a三级三级三级| 国产免费现黄频在线看| 视频区图区小说| 99国产精品一区二区三区| 国产精品麻豆人妻色哟哟久久| 纵有疾风起免费观看全集完整版| 亚洲国产中文字幕在线视频| 亚洲三区欧美一区| 女人久久www免费人成看片| 女人精品久久久久毛片| a 毛片基地| 久久av网站| 后天国语完整版免费观看| 99九九在线精品视频| 欧美97在线视频| 在线观看免费日韩欧美大片| 9191精品国产免费久久| 国产又爽黄色视频| 亚洲五月婷婷丁香| 高清不卡的av网站| 青春草视频在线免费观看| 如日韩欧美国产精品一区二区三区| 99国产精品99久久久久| 美女中出高潮动态图| 亚洲专区国产一区二区| 人体艺术视频欧美日本| 欧美日韩国产mv在线观看视频| 日韩,欧美,国产一区二区三区| 在线观看国产h片| 亚洲欧美色中文字幕在线| 国产精品秋霞免费鲁丝片| 国产成人欧美| 久9热在线精品视频| 一边摸一边做爽爽视频免费| 视频在线观看一区二区三区| 五月天丁香电影| 国产熟女欧美一区二区| 午夜福利影视在线免费观看| 久久久国产精品麻豆| 看免费av毛片| 捣出白浆h1v1| 校园人妻丝袜中文字幕| 精品少妇黑人巨大在线播放| 黄色怎么调成土黄色| 国产不卡av网站在线观看| 午夜久久久在线观看| 国产亚洲av片在线观看秒播厂| 久久久精品国产亚洲av高清涩受| 老司机靠b影院| 大香蕉久久成人网| 啦啦啦啦在线视频资源| 国产精品一国产av| 亚洲国产欧美日韩在线播放| 制服诱惑二区| 每晚都被弄得嗷嗷叫到高潮| 欧美精品一区二区免费开放| 超碰97精品在线观看| 欧美亚洲 丝袜 人妻 在线| a级毛片黄视频| 成年美女黄网站色视频大全免费| 欧美性长视频在线观看| 51午夜福利影视在线观看| 最近手机中文字幕大全| 精品国产乱码久久久久久男人| 99久久99久久久精品蜜桃| 国产三级黄色录像| 精品少妇内射三级| 大陆偷拍与自拍| 亚洲国产欧美网| 色94色欧美一区二区| 中文字幕av电影在线播放| 欧美人与善性xxx| 欧美久久黑人一区二区| 99国产精品免费福利视频| 国产精品一区二区免费欧美 | 日韩一本色道免费dvd| 国产精品久久久久久精品电影小说| 大话2 男鬼变身卡| 午夜激情av网站| 老鸭窝网址在线观看| 欧美日韩精品网址| 视频在线观看一区二区三区| 一区二区三区激情视频| 一区在线观看完整版| 飞空精品影院首页| 亚洲av成人精品一二三区| 三上悠亚av全集在线观看| av福利片在线| 精品亚洲乱码少妇综合久久| 成人国产av品久久久| 无限看片的www在线观看| 夫妻性生交免费视频一级片| 亚洲精品国产av成人精品| 色网站视频免费| 午夜日韩欧美国产| 美女脱内裤让男人舔精品视频| 午夜91福利影院| 欧美日韩综合久久久久久| 色播在线永久视频| 天天添夜夜摸| 午夜福利影视在线免费观看| 国产福利在线免费观看视频| 欧美日韩av久久| 色精品久久人妻99蜜桃| 国产精品一二三区在线看| 91国产中文字幕| 国产一级毛片在线| 日本wwww免费看| www.自偷自拍.com| 国产精品亚洲av一区麻豆| 国产精品久久久av美女十八| 日韩av免费高清视频| 两个人免费观看高清视频| 久久ye,这里只有精品| 亚洲,欧美精品.| 色视频在线一区二区三区| 成年人免费黄色播放视频| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 王馨瑶露胸无遮挡在线观看| 性色av乱码一区二区三区2| 五月天丁香电影| 国产精品 国内视频| 一边摸一边做爽爽视频免费| 欧美人与性动交α欧美精品济南到| 欧美日韩亚洲高清精品| 99久久精品国产亚洲精品| 国产精品麻豆人妻色哟哟久久| 色综合欧美亚洲国产小说| 欧美97在线视频| 多毛熟女@视频| 欧美老熟妇乱子伦牲交| 日日爽夜夜爽网站| 香蕉丝袜av| a级毛片黄视频| 国产黄色免费在线视频| 晚上一个人看的免费电影| 国产成人一区二区三区免费视频网站 | 九草在线视频观看| 人妻人人澡人人爽人人| 男的添女的下面高潮视频| 久久人人爽人人片av| 国产一区二区三区av在线| 欧美日韩亚洲高清精品| 最新在线观看一区二区三区 | 国产麻豆69| tube8黄色片| 午夜久久久在线观看| 国产91精品成人一区二区三区 | 免费久久久久久久精品成人欧美视频| 欧美老熟妇乱子伦牲交| 国产免费视频播放在线视频| 国产精品.久久久| 欧美乱码精品一区二区三区| 国产精品久久久久久人妻精品电影 | 一级片'在线观看视频| 国产午夜精品一二区理论片| 50天的宝宝边吃奶边哭怎么回事| 午夜福利免费观看在线| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 中文字幕亚洲精品专区| 各种免费的搞黄视频| 日韩电影二区| 黄色a级毛片大全视频| av有码第一页| 热99国产精品久久久久久7| 亚洲激情五月婷婷啪啪| 久久久久久久大尺度免费视频| 99热网站在线观看| 久热这里只有精品99| 亚洲av电影在线观看一区二区三区| 国产黄色视频一区二区在线观看| 国产有黄有色有爽视频| 曰老女人黄片| 青草久久国产| 日本av手机在线免费观看| 亚洲av电影在线进入| 99久久综合免费| 天堂中文最新版在线下载| 狂野欧美激情性bbbbbb| 午夜福利乱码中文字幕| 50天的宝宝边吃奶边哭怎么回事| 大陆偷拍与自拍| 老司机在亚洲福利影院| 欧美人与性动交α欧美精品济南到| 十八禁高潮呻吟视频| 美女扒开内裤让男人捅视频| 男女床上黄色一级片免费看| 91九色精品人成在线观看| 欧美精品人与动牲交sv欧美| 老汉色∧v一级毛片| 久久这里只有精品19| 男人添女人高潮全过程视频| 亚洲情色 制服丝袜| 女性生殖器流出的白浆| 久久精品亚洲熟妇少妇任你| 嫁个100分男人电影在线观看 | 欧美黄色淫秽网站| 久久天堂一区二区三区四区| 99国产精品99久久久久| 天天操日日干夜夜撸| 久久精品亚洲熟妇少妇任你| 久久九九热精品免费| 精品亚洲乱码少妇综合久久| 欧美少妇被猛烈插入视频| 在线天堂中文资源库| 一本综合久久免费| 麻豆av在线久日| 久久久久久久国产电影| 国产亚洲欧美精品永久| 亚洲成av片中文字幕在线观看| 午夜福利视频精品| 伊人亚洲综合成人网| 国产一级毛片在线| 国产一区亚洲一区在线观看| 一本一本久久a久久精品综合妖精| 男男h啪啪无遮挡| 天天操日日干夜夜撸| 黄色 视频免费看| 亚洲伊人久久精品综合| 午夜福利乱码中文字幕| 日日摸夜夜添夜夜爱| 亚洲av日韩在线播放| 建设人人有责人人尽责人人享有的| 午夜老司机福利片| 成人午夜精彩视频在线观看| 91精品伊人久久大香线蕉| 欧美在线一区亚洲| 在线亚洲精品国产二区图片欧美| 中文精品一卡2卡3卡4更新| 欧美人与善性xxx| 亚洲欧美清纯卡通| 成人手机av| 免费黄频网站在线观看国产| 丰满迷人的少妇在线观看| 90打野战视频偷拍视频| 99国产精品99久久久久| 最近中文字幕2019免费版| 亚洲自偷自拍图片 自拍| 高清黄色对白视频在线免费看| 在现免费观看毛片| 亚洲一码二码三码区别大吗| 美女国产高潮福利片在线看| 国产免费现黄频在线看| 中文字幕av电影在线播放| 免费高清在线观看日韩| 国产精品熟女久久久久浪| 男的添女的下面高潮视频| 一级黄色大片毛片| 在线观看免费日韩欧美大片| 性色av一级| 欧美日韩视频精品一区| 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 亚洲精品国产区一区二| 波多野结衣一区麻豆| 在线亚洲精品国产二区图片欧美| 操出白浆在线播放| 乱人伦中国视频| 亚洲成国产人片在线观看| 午夜影院在线不卡| 午夜久久久在线观看| 日韩电影二区| 国产精品 欧美亚洲| 宅男免费午夜| 亚洲精品久久午夜乱码| 捣出白浆h1v1| 亚洲成人手机| 精品国产乱码久久久久久男人| 99国产精品一区二区三区| av网站免费在线观看视频| 麻豆乱淫一区二区| 青草久久国产| 天天躁夜夜躁狠狠躁躁| 国产97色在线日韩免费| 男女无遮挡免费网站观看| 一级毛片电影观看| 亚洲av美国av| 国产免费视频播放在线视频| 精品久久久久久久毛片微露脸 | 人妻人人澡人人爽人人| 丝袜人妻中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 久久青草综合色| 欧美精品人与动牲交sv欧美| 亚洲伊人色综图| 丰满迷人的少妇在线观看| 欧美xxⅹ黑人| 纯流量卡能插随身wifi吗| www.熟女人妻精品国产| 大陆偷拍与自拍| 亚洲av综合色区一区| 一级片'在线观看视频| 日韩熟女老妇一区二区性免费视频| 亚洲精品第二区| 一本—道久久a久久精品蜜桃钙片| 一二三四社区在线视频社区8| 丝袜脚勾引网站| 国产成人91sexporn| 国产成人一区二区在线| av不卡在线播放| 国产极品粉嫩免费观看在线| 亚洲视频免费观看视频| 免费在线观看视频国产中文字幕亚洲 | 久久这里只有精品19| 亚洲国产精品999| 国产精品熟女久久久久浪| 久久精品久久久久久噜噜老黄| 蜜桃在线观看..| 在现免费观看毛片| 少妇人妻 视频| 深夜精品福利| 国产片特级美女逼逼视频| 日韩制服骚丝袜av| 99re6热这里在线精品视频| 亚洲一区中文字幕在线|