• <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)保法官
    氣候變化:法官的作用
    成人漫画全彩无遮挡| 只有这里有精品99| 尤物成人国产欧美一区二区三区| 在线观看国产h片| 欧美三级亚洲精品| 精品一区二区三卡| 99热全是精品| 最近中文字幕高清免费大全6| 久久99热这里只有精品18| 男的添女的下面高潮视频| 国产伦在线观看视频一区| 中文字幕制服av| 综合色丁香网| h日本视频在线播放| 成人综合一区亚洲| 国产 精品1| 日韩三级伦理在线观看| 日韩一本色道免费dvd| 日韩视频在线欧美| 日本黄色片子视频| 能在线免费看毛片的网站| 日韩av免费高清视频| 国产精品福利在线免费观看| 免费看不卡的av| 欧美成人午夜免费资源| 国产精品熟女久久久久浪| 国产高清有码在线观看视频| 免费av毛片视频| 午夜福利高清视频| 在线精品无人区一区二区三 | 国产免费一区二区三区四区乱码| 一级片'在线观看视频| 成人亚洲欧美一区二区av| 国产在线男女| 一本色道久久久久久精品综合| 亚洲国产精品999| 日日啪夜夜爽| 日韩一区二区三区影片| 午夜亚洲福利在线播放| 99久久精品热视频| 亚洲真实伦在线观看| 欧美区成人在线视频| 亚洲天堂国产精品一区在线| 一本久久精品| 欧美日韩视频高清一区二区三区二| 亚洲综合色惰| 观看免费一级毛片| 亚洲国产精品国产精品| 亚洲av免费高清在线观看| 免费av毛片视频| 国内少妇人妻偷人精品xxx网站| 亚洲,欧美,日韩| 男女那种视频在线观看| 久久久久精品性色| 少妇猛男粗大的猛烈进出视频 | 欧美xxxx性猛交bbbb| 国产综合精华液| 亚洲精品乱久久久久久| av在线观看视频网站免费| 国产精品一区二区在线观看99| 一级毛片aaaaaa免费看小| 性色av一级| 日韩欧美精品免费久久| 国国产精品蜜臀av免费| 久久久久精品久久久久真实原创| 在线a可以看的网站| 中国美白少妇内射xxxbb| 又爽又黄a免费视频| 亚洲成人久久爱视频| 久久韩国三级中文字幕| 欧美bdsm另类| 大话2 男鬼变身卡| 97超碰精品成人国产| 亚洲国产成人一精品久久久| 国语对白做爰xxxⅹ性视频网站| 成人特级av手机在线观看| 国产精品秋霞免费鲁丝片| 22中文网久久字幕| 人人妻人人看人人澡| 国产91av在线免费观看| 熟女人妻精品中文字幕| 一级毛片我不卡| 中文精品一卡2卡3卡4更新| 亚洲av国产av综合av卡| av福利片在线观看| 日本熟妇午夜| 亚洲丝袜综合中文字幕| 18禁动态无遮挡网站| 亚洲精品成人av观看孕妇| 好男人视频免费观看在线| 午夜激情久久久久久久| 看非洲黑人一级黄片| 国产女主播在线喷水免费视频网站| 久久精品国产亚洲av涩爱| 久久精品国产a三级三级三级| 色视频在线一区二区三区| 亚洲天堂av无毛| av专区在线播放| 香蕉精品网在线| 女人十人毛片免费观看3o分钟| 欧美精品国产亚洲| 亚洲天堂av无毛| 午夜福利在线观看免费完整高清在| 亚洲成色77777| 十八禁网站网址无遮挡 | 亚洲内射少妇av| 亚洲av免费在线观看| 国产国拍精品亚洲av在线观看| 成人国产麻豆网| 免费看a级黄色片| 日韩中字成人| 亚洲欧美精品专区久久| 亚洲欧美精品专区久久| 永久网站在线| 熟女电影av网| 高清午夜精品一区二区三区| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 自拍欧美九色日韩亚洲蝌蚪91 | 久久精品夜色国产| 亚洲精品日本国产第一区| 日本一二三区视频观看| 日本猛色少妇xxxxx猛交久久| 国产精品福利在线免费观看| 2022亚洲国产成人精品| 日韩精品有码人妻一区| 亚州av有码| 国产综合精华液| 少妇猛男粗大的猛烈进出视频 | 夜夜看夜夜爽夜夜摸| 男女边摸边吃奶| 中国三级夫妇交换| 18禁在线无遮挡免费观看视频| 精品少妇久久久久久888优播| 久久国内精品自在自线图片| 亚洲婷婷狠狠爱综合网| 亚洲国产精品成人综合色| 能在线免费看毛片的网站| 国语对白做爰xxxⅹ性视频网站| 一区二区三区免费毛片| 18禁在线无遮挡免费观看视频| 午夜福利视频精品| 亚洲精品久久午夜乱码| 亚洲国产精品成人综合色| 少妇猛男粗大的猛烈进出视频 | 男女无遮挡免费网站观看| 特大巨黑吊av在线直播| 亚洲精品成人av观看孕妇| 男插女下体视频免费在线播放| 欧美日韩国产mv在线观看视频 | 春色校园在线视频观看| 国产成人a∨麻豆精品| 黄色怎么调成土黄色| 舔av片在线| 91久久精品国产一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 2022亚洲国产成人精品| 日韩国内少妇激情av| 午夜激情福利司机影院| 国产精品一区www在线观看| 国产av国产精品国产| 午夜免费鲁丝| 国产男人的电影天堂91| 亚洲国产欧美人成| 国产欧美另类精品又又久久亚洲欧美| 日韩一本色道免费dvd| 国产有黄有色有爽视频| 99久久精品一区二区三区| 国产亚洲av片在线观看秒播厂| 观看免费一级毛片| 99热这里只有是精品在线观看| 亚洲欧美中文字幕日韩二区| 国产人妻一区二区三区在| 丝袜脚勾引网站| 插阴视频在线观看视频| 午夜免费鲁丝| 偷拍熟女少妇极品色| 欧美少妇被猛烈插入视频| 亚洲天堂av无毛| 国产精品偷伦视频观看了| 国内精品美女久久久久久| 永久免费av网站大全| 亚洲激情五月婷婷啪啪| freevideosex欧美| 日韩精品有码人妻一区| 久久国产乱子免费精品| 少妇人妻一区二区三区视频| 欧美日韩视频精品一区| 精品一区二区免费观看| 少妇人妻 视频| 在线观看三级黄色| 丰满乱子伦码专区| 一级二级三级毛片免费看| 国精品久久久久久国模美| 一级毛片我不卡| 在线 av 中文字幕| 免费在线观看成人毛片| 国产精品伦人一区二区| 免费黄网站久久成人精品| 日日啪夜夜爽| videos熟女内射| 免费少妇av软件| 日韩一区二区视频免费看| 日韩电影二区| 亚洲在线观看片| 久热这里只有精品99| 国产在线男女| 伦精品一区二区三区| 国产成人freesex在线| 少妇猛男粗大的猛烈进出视频 | 人妻系列 视频| 国产精品一二三区在线看| 99热这里只有精品一区| 女的被弄到高潮叫床怎么办| 韩国高清视频一区二区三区| 欧美 日韩 精品 国产| 久热久热在线精品观看| 日本三级黄在线观看| 精品亚洲乱码少妇综合久久| 久久久久久久亚洲中文字幕| 国产伦精品一区二区三区四那| 欧美老熟妇乱子伦牲交| 午夜视频国产福利| 欧美老熟妇乱子伦牲交| 黄色配什么色好看| 国产高清不卡午夜福利| 男人舔奶头视频| 久久久成人免费电影| 精品一区二区免费观看| 亚洲怡红院男人天堂| 亚州av有码| 精品午夜福利在线看| 一级二级三级毛片免费看| 欧美激情在线99| 大香蕉97超碰在线| 狂野欧美激情性bbbbbb| 亚洲人成网站在线播| 女的被弄到高潮叫床怎么办| 女人被狂操c到高潮| 亚洲美女搞黄在线观看| 国产 精品1| 中文字幕制服av| 观看免费一级毛片| 看黄色毛片网站| 在线播放无遮挡| 亚洲最大成人中文| av福利片在线观看| 一级毛片黄色毛片免费观看视频| 国产精品一区二区性色av| 国产免费一级a男人的天堂| 欧美精品人与动牲交sv欧美| 日本wwww免费看| 舔av片在线| 国产免费视频播放在线视频| 久久久久精品性色| 久久人人爽人人片av| 极品教师在线视频| 老女人水多毛片| 18+在线观看网站| 在线精品无人区一区二区三 | 亚洲精品一区蜜桃| 亚洲三级黄色毛片| 国产精品99久久久久久久久| 亚洲最大成人中文| 日本黄大片高清| 欧美bdsm另类| 九色成人免费人妻av| 啦啦啦中文免费视频观看日本| 日韩亚洲欧美综合| 女的被弄到高潮叫床怎么办| 免费人成在线观看视频色| 在线精品无人区一区二区三 | 国产毛片a区久久久久| 18禁裸乳无遮挡动漫免费视频 | 国产精品久久久久久精品电影小说 | 亚洲第一区二区三区不卡| 日韩人妻高清精品专区| 美女主播在线视频| 国产精品一及| 欧美一区二区亚洲| 丰满少妇做爰视频| 好男人在线观看高清免费视频| 国产乱来视频区| 国产综合懂色| 国产一区二区三区综合在线观看 | 亚洲高清免费不卡视频| 国产精品熟女久久久久浪| 韩国高清视频一区二区三区| 国产色婷婷99| 精品酒店卫生间| 永久免费av网站大全| 免费电影在线观看免费观看| 啦啦啦中文免费视频观看日本| av女优亚洲男人天堂| 晚上一个人看的免费电影| 一级毛片aaaaaa免费看小| 久久女婷五月综合色啪小说 | xxx大片免费视频| 亚洲色图综合在线观看| 麻豆成人午夜福利视频| 成人欧美大片| 国产成人精品久久久久久| 视频区图区小说| 激情 狠狠 欧美| 丝袜脚勾引网站| 中文在线观看免费www的网站| 五月开心婷婷网| 久久久久久久久大av| 激情 狠狠 欧美| 午夜福利视频1000在线观看| 99热这里只有精品一区| 国产伦在线观看视频一区| 亚洲一区二区三区欧美精品 | 日本三级黄在线观看| 午夜免费男女啪啪视频观看| 我的老师免费观看完整版| 美女国产视频在线观看| 国精品久久久久久国模美| 极品教师在线视频| 国产免费又黄又爽又色| 欧美最新免费一区二区三区| 深夜a级毛片| 观看美女的网站| 午夜激情福利司机影院| 国产探花在线观看一区二区| 欧美三级亚洲精品| 九色成人免费人妻av| 久久99热6这里只有精品| 日韩视频在线欧美| 纵有疾风起免费观看全集完整版| 男女边吃奶边做爰视频| 午夜激情福利司机影院| av在线亚洲专区| 最新中文字幕久久久久| 欧美成人a在线观看| 99热这里只有精品一区| 国产精品久久久久久久电影| 男女国产视频网站| 亚洲精品国产av蜜桃| 欧美三级亚洲精品| 国国产精品蜜臀av免费| 秋霞在线观看毛片| 免费看光身美女| 日本午夜av视频| 亚洲欧美日韩卡通动漫| 成人鲁丝片一二三区免费| 乱码一卡2卡4卡精品| 伦理电影大哥的女人| 久久精品国产自在天天线| 日本一二三区视频观看| 免费黄频网站在线观看国产| 91精品伊人久久大香线蕉| 国产在视频线精品| 天堂网av新在线| 制服丝袜香蕉在线| 欧美精品国产亚洲| 欧美xxxx性猛交bbbb| 久久久久久久久久成人| 亚洲成人精品中文字幕电影| 亚洲精品影视一区二区三区av| 国产在线一区二区三区精| 亚洲精品aⅴ在线观看| 校园人妻丝袜中文字幕| 国产乱人偷精品视频| 大话2 男鬼变身卡| 国产亚洲91精品色在线| 久久精品夜色国产| 亚洲av中文字字幕乱码综合| 亚洲精品亚洲一区二区| 青春草亚洲视频在线观看| 国产精品久久久久久精品古装| 内射极品少妇av片p| 97人妻精品一区二区三区麻豆| 一区二区三区四区激情视频| 亚洲久久久久久中文字幕| av天堂中文字幕网| 日韩大片免费观看网站| 男人和女人高潮做爰伦理| 午夜免费观看性视频| 亚洲自拍偷在线| 色播亚洲综合网| av在线老鸭窝| 久久久午夜欧美精品| 欧美国产精品一级二级三级 | 日本午夜av视频| 久久韩国三级中文字幕| 欧美97在线视频| 亚洲自拍偷在线| 丰满少妇做爰视频| 日韩电影二区| 国产真实伦视频高清在线观看| 女人十人毛片免费观看3o分钟| 国产亚洲av片在线观看秒播厂| 九九在线视频观看精品| 国产亚洲91精品色在线| 国产中年淑女户外野战色| 久久精品国产亚洲av涩爱| 26uuu在线亚洲综合色| 欧美xxxx黑人xx丫x性爽| 精品视频人人做人人爽| 亚洲精品久久久久久婷婷小说| 精品少妇久久久久久888优播| 别揉我奶头 嗯啊视频| 国产精品一区二区性色av| 嫩草影院精品99| 免费看日本二区| 国产精品一区二区三区四区免费观看| www.色视频.com| 国产精品嫩草影院av在线观看| 丝袜喷水一区| 午夜日本视频在线| 天美传媒精品一区二区| 中文字幕制服av| 一级毛片黄色毛片免费观看视频| 少妇的逼好多水| 麻豆精品久久久久久蜜桃| 国产精品一区二区在线观看99| 成人国产麻豆网| 亚洲欧美中文字幕日韩二区| 亚洲欧美日韩东京热| 日韩大片免费观看网站| 51国产日韩欧美| 99re6热这里在线精品视频| av在线观看视频网站免费| 毛片一级片免费看久久久久| 黄色欧美视频在线观看| av在线观看视频网站免费| 六月丁香七月| 久久精品熟女亚洲av麻豆精品| 日产精品乱码卡一卡2卡三| 久久久久久久精品精品| av在线app专区| 成年女人看的毛片在线观看| 男女国产视频网站| 精品99又大又爽又粗少妇毛片| 麻豆久久精品国产亚洲av| 久久影院123| 高清在线视频一区二区三区| 男的添女的下面高潮视频| 日本爱情动作片www.在线观看| 精品人妻一区二区三区麻豆| 99视频精品全部免费 在线| 亚洲av福利一区| 国产69精品久久久久777片| av福利片在线观看| 国产精品爽爽va在线观看网站| 高清毛片免费看| 精品人妻偷拍中文字幕| 六月丁香七月| 欧美高清成人免费视频www| av免费观看日本| 丝袜脚勾引网站| 国产成人精品福利久久| 中文精品一卡2卡3卡4更新| 美女高潮的动态| 在线a可以看的网站| 国产亚洲最大av| 中文字幕制服av| 成人亚洲欧美一区二区av| 国产成人福利小说| 香蕉精品网在线| 日本-黄色视频高清免费观看| 亚洲国产欧美人成| 国产精品国产三级国产专区5o| 亚洲最大成人中文| 午夜激情久久久久久久| 国产 一区 欧美 日韩| 国产高潮美女av| 国产免费视频播放在线视频| 一级爰片在线观看| 亚洲欧美精品专区久久| 最近手机中文字幕大全| 男女无遮挡免费网站观看| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 伊人久久精品亚洲午夜| 97精品久久久久久久久久精品| 少妇裸体淫交视频免费看高清| 精品亚洲乱码少妇综合久久| 国产精品一区二区三区四区免费观看| 大片电影免费在线观看免费| 亚洲图色成人| 日韩一本色道免费dvd| 国产精品久久久久久av不卡| 午夜精品一区二区三区免费看| 婷婷色综合大香蕉| 99久国产av精品国产电影| 国产欧美亚洲国产| 夫妻午夜视频| 精品亚洲乱码少妇综合久久| 国产一区二区三区av在线| 2021少妇久久久久久久久久久| 黄色配什么色好看| 丰满乱子伦码专区| 丝袜喷水一区| eeuss影院久久| 麻豆精品久久久久久蜜桃| 肉色欧美久久久久久久蜜桃 | 韩国av在线不卡| 高清午夜精品一区二区三区| 色视频在线一区二区三区| 免费人成在线观看视频色| 丰满人妻一区二区三区视频av| 成人亚洲精品av一区二区| 亚洲av国产av综合av卡| 国产人妻一区二区三区在| 99精国产麻豆久久婷婷| 免费观看性生交大片5| 人人妻人人爽人人添夜夜欢视频 | 99热全是精品| 22中文网久久字幕| 免费观看的影片在线观看| 女人久久www免费人成看片| 18禁在线播放成人免费| 久久久精品免费免费高清| 国产亚洲91精品色在线| 伊人久久国产一区二区| 免费播放大片免费观看视频在线观看| 韩国高清视频一区二区三区| 国产女主播在线喷水免费视频网站| 亚洲四区av| 亚洲国产高清在线一区二区三| 国产精品久久久久久av不卡| 伦理电影大哥的女人| 一级二级三级毛片免费看| 亚洲欧美成人精品一区二区| 亚洲国产欧美人成| 一级毛片aaaaaa免费看小| 亚洲国产欧美在线一区| 国产色婷婷99| 欧美精品国产亚洲| 99热这里只有是精品50| 亚洲av免费在线观看| 国产亚洲av片在线观看秒播厂| 菩萨蛮人人尽说江南好唐韦庄| 国产精品伦人一区二区| 国产午夜精品久久久久久一区二区三区| 在线 av 中文字幕| 人体艺术视频欧美日本| 久久6这里有精品| 97超碰精品成人国产| 国产av不卡久久| 乱码一卡2卡4卡精品| 久久久久精品久久久久真实原创| 麻豆乱淫一区二区| 肉色欧美久久久久久久蜜桃 | 日韩在线高清观看一区二区三区| 最近最新中文字幕免费大全7| 欧美 日韩 精品 国产| 国产亚洲5aaaaa淫片| 国产精品久久久久久精品电影小说 | 国精品久久久久久国模美| 少妇的逼水好多| av又黄又爽大尺度在线免费看| 内地一区二区视频在线| 日产精品乱码卡一卡2卡三| 日韩亚洲欧美综合| 色网站视频免费| 久久99热这里只频精品6学生| 国产在线男女| 中文精品一卡2卡3卡4更新| 99久久人妻综合| 五月伊人婷婷丁香| 国产乱来视频区| 91久久精品电影网| 国产亚洲av嫩草精品影院| 亚洲经典国产精华液单| 啦啦啦中文免费视频观看日本| 亚洲欧美成人精品一区二区| 日韩av免费高清视频| av黄色大香蕉| 亚洲欧美成人综合另类久久久| 色视频www国产| 国产精品国产三级专区第一集| 日韩 亚洲 欧美在线| 三级国产精品欧美在线观看| 能在线免费看毛片的网站| 五月开心婷婷网| 国精品久久久久久国模美| 国产伦在线观看视频一区| 国产高清国产精品国产三级 | 观看美女的网站| 91精品一卡2卡3卡4卡| 亚洲最大成人av| 小蜜桃在线观看免费完整版高清| 一级a做视频免费观看| 国产探花在线观看一区二区| 亚洲自偷自拍三级| 亚洲av成人精品一区久久| 高清日韩中文字幕在线| 乱码一卡2卡4卡精品| 国产成人免费观看mmmm| av线在线观看网站| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 国产精品无大码| 韩国高清视频一区二区三区| 神马国产精品三级电影在线观看| 午夜爱爱视频在线播放| 欧美精品一区二区大全| 国产 精品1| 男女边摸边吃奶| 一级黄片播放器| 一区二区三区精品91| 三级国产精品片| 日韩av在线免费看完整版不卡| 午夜激情福利司机影院| 婷婷色麻豆天堂久久| 麻豆精品久久久久久蜜桃| 日韩欧美精品v在线| 亚洲国产精品国产精品| 高清日韩中文字幕在线| 欧美97在线视频| 免费大片18禁| 美女脱内裤让男人舔精品视频|