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

    氣候變化條件下苦參在我國潛在分布區(qū)的預(yù)測分析

    2022-04-08 01:08:39張濤胡菀賈天嬌趙三增孔丹宇劉毅
    廣西植物 2022年3期
    關(guān)鍵詞:苦參氣候變化

    張濤 胡菀 賈天嬌 趙三增 孔丹宇 劉毅

    摘 要:為了解氣候變化情景下苦參在中國的潛在分布區(qū)變化,探討生物氣候因子與苦參適宜分布格局的關(guān)系。該文通過收集苦參的地理分布點并結(jié)合19項生態(tài)因子,運用最大熵模型(MaxEnt)和地理信息系統(tǒng)(ArcGIS)對末次盛冰期、當(dāng)前氣候、未來氣候三種氣候情景下苦參在我國適生區(qū)的分布格局進行模擬,并分析影響苦參生長的主導(dǎo)生態(tài)因子。結(jié)果表明:(1)當(dāng)前氣候條件下,最暖季度平均溫度(bio10)、最濕季度降水量(bio16)是影響苦參分布的主導(dǎo)氣候因子。(2)苦參的適宜生境占我國國土總面積的35.90%,高適生區(qū)主要分布在我國800 mm等降水量線與400 mm等降水量線之間的地區(qū)。(3)在末次盛冰期至未來氣候情境下,苦參在我國的適宜生境面積逐漸減少,呈現(xiàn)遞減的趨勢,且整體有向高緯度地區(qū)移動的趨勢。(4)當(dāng)前氣候條件下苦參的適宜生境面積較末次盛冰期減少0.3%,未來氣候條件下苦參的適宜生境面積較當(dāng)前氣候減少0.75%。全球氣候變暖對苦參潛在分布區(qū)具有一定負(fù)面影響,造成苦參的適宜生境縮減、實際生態(tài)位變窄,一定程度上不利于苦參的生長。該研究結(jié)果為苦參的野生保育及人工栽培的合理布局提供了理論依據(jù)。

    關(guān)鍵詞: 苦參, MaxEnt, 氣候變化, 生態(tài)因子, 潛在分布區(qū)

    中圖分類號:Q948.2

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

    文章編號:1000-3142(2022)03-0349-14

    Prediction of potential distribution of Sophora flavescens

    in China under climate change

    ZHANG Tao1, HU Wan1, JIA Tianjiao1, ZHAO Sanzeng2, KONG Danyu1, LIU Yi1*

    ( 1. Lushan Botanical Garden, Jiangxi Province and Chinese Academy of Sciences, Lushan 332900, Jiangxi, China;

    2. College of Tropical Crops, Hainan University, Haikou 570228, China )

    Abstract:In order to understand the changes of potential distribution of Sophora flavescens in China under the climate change and to explore the relationship between bioclimatic factors and the suitable distribution pattern of S. flavescens, we investigated the distribution pattern of S. flavescens in China under the three climate scenarios (last glacial maximum, current climate, and future climate), and analyzed the dominant ecological factors affecting the growth of S. flavescens with Maximum Entropy Model (MaxEnt) and Geographic Information System (ArcGIS) by collecting the geographical distribution points of S. flavescens combining with 19 ecological factors. The results were as follows: (1) Under the current climate conditions, the average temperature in the warmest quarter (bio10) and the precipitation in the wettest quarter (bio16) were the dominant climate factors affecting the distribution of S. flavescens. (2) The suitable habitat of S. flavescens accounted for 35.90% of the total land area of China. The optimum areas of S. flavescens were mainly distributed between the 800 mm isopyet and the 400 mm isopyet in China. (3) From the last glacial maximum to the future climate, the suitable habitat area of S. flavescens in China gradually decreased, and showed a trend of moving to high latitude. (4) The suitable habitat area of S. flavescens under the current climate was 0.3% less than that in the last glacial maximum, and the suitable habitat area of S. flavescens under the future climate was 0.75% less than that in the current climate. In conclusion, global warming plays a negative role in the potential distribution area of S. flavescens, and leads to the reduction of suitable habitat and the narrowing of actual niche of S. flavescens, which is not conducive to the growth of S. flavescens to a certain extent. The results of this study provide the theoretical basis for the rational distribution of wild conservation and artificial cultivation of S. flavescens.

    Key words: Sophora flavescens, MaxEnt, climate change, ecological factor, potential distribution area

    苦參(Sophora flavescens)為豆科(Legumi-nosae)槐屬(Sophora)多年生落葉亞灌木,稀呈灌木狀或草本,別名地槐、白莖地骨、山槐、野槐等(中國科學(xué)院中國植物志編輯委員會,1994)。《中國藥典》記載臨床上常以苦參的干燥根入藥,具有清熱燥濕、殺蟲、利尿等功效(國家藥典委員會,2020)?,F(xiàn)代藥理學(xué)研究表明,苦參具有抗腫瘤、抗心律失常、調(diào)節(jié)免疫、抗菌、抗病毒等多種藥理活性,特別是對腫瘤細(xì)胞有顯著的抑制作用(鐘贛生,2016)。除用于中藥方劑外,苦參還是制藥工業(yè)中的常用制劑、日用品、生物農(nóng)藥和獸藥的主要原料。但是,隨著我國中藥產(chǎn)業(yè)的迅速發(fā)展,對中藥材的需求量日益增加,野生資源已不能滿足中藥材產(chǎn)業(yè)的發(fā)展需求??鄥⒆鳛槌S弥兴帲猛緩V泛, 用量逐年遞增,導(dǎo)致野生苦參資源嚴(yán)重匱乏(紀(jì)瑛等,2011)。因此,厘清苦參的全國分布格局,分析其生態(tài)適宜性意義重大。傳統(tǒng)的資源調(diào)查方法雖是開展野外調(diào)研,但野外調(diào)查獲得的數(shù)據(jù)不足以充分反映物種的整體分布格局。而藥用植物資源的分布與地理環(huán)境有著直接關(guān)系,受到諸多環(huán)境因子的制約,以及人類活動的影響,是在特定空間下的產(chǎn)物,具有空間信息特征(郭蘭萍等,2005)。物種能夠生存和繁衍后代的所有條件的總和稱為生態(tài)位,是定義物種的最小分布單元(喬慧捷等,2013)。生態(tài)位模型可利用物種已知的分布數(shù)據(jù)和相關(guān)環(huán)境變量,根據(jù)統(tǒng)計學(xué)的運算方法推算被研究物種的生態(tài)需求,將此運算結(jié)果投射至不同的空間和時間中用以預(yù)測被研究物種的潛在分布(朱耿平等,2013)。

    近年來,生態(tài)位模型越來越多的應(yīng)用于入侵生物學(xué)、保護生物學(xué)、氣候變化對物種分布的影響以及傳染病空間傳播等研究領(lǐng)域中(朱耿平等,2013)?;诓煌睦碚摶A(chǔ)、基礎(chǔ)數(shù)據(jù)和分析方式,已開發(fā)出許多成熟的生態(tài)位模型算法,主要的模型有生物氣候分析系統(tǒng)(BIOCLIM)、生態(tài)因子分析模型(ENFA)、基于規(guī)則集的遺傳算法模型(GARP)、最大熵模型(MaxEnt)等,其中MaxEnt模型是目前應(yīng)用最為廣泛的生態(tài)位模型。通過對比16種生態(tài)位模型的預(yù)測能力,結(jié)果發(fā)現(xiàn)MaxEnt 模型的預(yù)測能力最高(Elith et al., 2006)。MaxEnt模型通過物種的分布數(shù)據(jù)和環(huán)境數(shù)據(jù),找出物種分布的最大熵,從而對物種的分布進行預(yù)測,相較于其他常見的生態(tài)位模型,最大熵模型精度優(yōu)于其他類型的模型,并且在物種樣本量很小甚至樣本量數(shù)據(jù)殘缺的情境下,仍具有良好的預(yù)測效果和可信度(王運生等,2007;邢丁亮等,2011;車樂等,2014)。該模型的穩(wěn)定性好且預(yù)測的結(jié)果與物種的實際分布基本吻合(劉艷梅等,2018)。MaxEnt 模型作為評價物種棲息地生境質(zhì)量中最常用的模型,被廣泛運用于農(nóng)作物適宜區(qū)預(yù)測、動植物潛在生境評價、外來入侵物種風(fēng)險評估和藥用植物潛在生境分布等眾多研究領(lǐng)域。從魚類(Kim et al.,2020)到兩棲類(黃勇杰等,2017),從昆蟲(Kamel et al., 2021; Saha et al., 2021; Ji et al., 2021)到鳥類(Jha KK & Jha R, 2021),從植物類(Liu et al., 2021; Boral et al., 2021; Zeng et al., 2021)到大型獸類(Rozhnov et al., 2021; 溫平等, 2021)。從全球的大尺度范圍到物種廊道的小尺度范圍都取得了較好的預(yù)測結(jié)果。在藥用植物研究領(lǐng)域已被成功運用于如何首烏(Fallopia multiflora)(潘石玉等,2016)、甘遂(Euphorbia kansui) (盧有媛等,2018)、秦艽(Gentiana macrophylla)(檀逸虹等,2020)、京大戟(Euphorbia pekinensis)(裴蘇婷等,2021)等藥用植物的潛在分布預(yù)測分析,且預(yù)測效果良好。因此,本研究選用MaxEnt 模型來進行氣候變化條件下苦參在我國潛在分布區(qū)的模擬。

    當(dāng)前苦參的人工栽培已成為保護苦參野生資源,滿足中藥產(chǎn)業(yè)對苦參大量需求的最重要途徑之一,而對苦參進行在全國范圍內(nèi)生態(tài)適宜性區(qū)劃的研究尚未見有報道。因此,通過對苦參進行生態(tài)位分析,厘清苦參在全國的分布格局,掌握其在全國的生態(tài)適宜性分布區(qū)域情況,可為苦參種質(zhì)資源調(diào)查、野生種源保護以及分子譜系地理學(xué)等研究提供基礎(chǔ)資料。同時,可以在全國范圍內(nèi)劃分合理的道地產(chǎn)區(qū),指導(dǎo)苦參在全國的產(chǎn)業(yè)發(fā)展規(guī)劃,建立符合GAP標(biāo)準(zhǔn)的苦參原藥材種植基地,實現(xiàn)質(zhì)量穩(wěn)定的道地藥材規(guī)范化、規(guī)?;a(chǎn)。本研究以中國為研究區(qū)域,以苦參為研究對象,選取19個生物氣候因子作為環(huán)境變量,依托中國數(shù)字標(biāo)本館(https://www.cvh.ac.cn/) 收集的苦參分布信息,采用MaxEnt模型和ArcGIS空間分析技術(shù),通過模型模擬預(yù)測該物種在三個不同歷史時期的潛在分布區(qū)和適宜生境,擬探討以下問題:(1)預(yù)測分析苦參的潛在分布區(qū)域的時空分布格局;(2)分析苦參潛在適宜分布格局與環(huán)境因子之間的聯(lián)系;(3)提出苦參在我國范圍內(nèi)的保護和開發(fā)利用的科學(xué)依據(jù)。

    1 材料與方法

    1.1 樣本信息

    查閱中國數(shù)字標(biāo)本館(https://www.cvh.ac.cn/) 已收集苦參標(biāo)本的分布信息,刪除重復(fù)和記錄不詳細(xì)樣本分布信息,對有經(jīng)緯度記載的數(shù)據(jù)直接使用,僅有位置記錄的分布點以最小可查行政單位為準(zhǔn),通過百度地圖坐標(biāo)拾取工具獲得經(jīng)緯度數(shù)據(jù),在ArcGIS中對獲得的數(shù)據(jù)進行篩選,一個柵格中只保留一個分布信息,確保分布點間的距離均大于1 km,避免過度擬合造成的誤差,最終獲得725份苦參樣本分布信息,詳見圖1。

    1.2 環(huán)境數(shù)據(jù)獲取

    本研究使用的三個歷史時期(末次盛冰期、當(dāng)前時期、未來時期)氣候數(shù)據(jù)均下載自世界氣候數(shù)據(jù)庫(WorldClim,http: //worldclim.org) ,包括19 個生物氣候變量(bio1~bio19)(表1)。末次盛冰期選取符合中國氣候類型的CCSM4, 選擇最小分辨率2.5′(5 km×5 km);當(dāng)前氣候年份選擇1970—2000(version 2.1);未來年份選擇2070,未來氣候選取IPCC第5份報告中描述的未來四種氣候變暖情景中代表高排放情景的RCP8.5(高)。當(dāng)前和未來的氣候圖層分辨率均為30″(1 km × 1 km)。中國底圖下載自國家基礎(chǔ)地理信息系統(tǒng)網(wǎng)站(http://www.ngcc.cn/ngcc/)。

    1.3 MaxEnt 模型設(shè)定及評價

    MaxEnt建模參數(shù)設(shè)定:選用MaxEnt軟件(V3.3.3K),隨機測試比例20%,最大迭代次數(shù)1 × 104,設(shè)置響應(yīng)曲線、ROC曲線和刀切法。以響應(yīng)曲線評價生態(tài)因子適宜度范圍,以ROC曲線及曲線下面積(即AUC值)評價模型精度,以刀切法檢驗生態(tài)因子的權(quán)重。

    1.4 環(huán)境變量篩選

    選取模型輸出結(jié)果中貢獻率大、權(quán)重高的環(huán)境變量為主導(dǎo)環(huán)境變量。

    1.5 生態(tài)適宜區(qū)劃分采用GIS

    首先將建模得到的輸出結(jié)果中的ASCII文件通過ArcGIS中的格式轉(zhuǎn)化工具轉(zhuǎn)化成柵格文件,然后通過重分類工具將圖層進行重分類,從而進行適生等級的劃分。輸出結(jié)果顯示的是不同地區(qū)的生境適宜性,用P表示,以概率的形式體現(xiàn),取值0~1。參考2007年聯(lián)合國政府間氣候變化專門委員會(IPCC)第四次評估報告對于不確定性的等級劃分,本文對研究對象在全國的適生等級進行如下劃分:P<0.33為非適生;0.33≤P<0.66為低適生;P≥0.66為高度適生。

    2 結(jié)果與分析

    2.1 MaxEnt 模型預(yù)測精度

    MaxEnt模型預(yù)測精度以ROC曲線下面積(即AUC值)為衡量指標(biāo),AUC值越大表明模型預(yù)測精度越高。評價標(biāo)準(zhǔn): 0.5~0.6即模型預(yù)測失敗;0.6~0.7較差;0.7~0.8一般;0.8~0.9好;0.9~1非常好?;诳鄥?25個樣本分布點,利用MaxEnt模型在過去、當(dāng)前和未來三種氣候模式下對苦參適生區(qū)進行預(yù)測,結(jié)果顯示過去、現(xiàn)在、未來三種情境的氣候條件下,模型運算的測試數(shù)據(jù)集和訓(xùn)練數(shù)據(jù)集的AUC平均值分別為0.858(圖2:A)、0.857(圖2:B)、0.863(圖2:C),表明由MaxEnt模型預(yù)測的苦參生態(tài)適宜性具有較好的準(zhǔn)確度和可信度,可用于苦參潛在適生區(qū)分布的研究。

    2.2 苦參潛在分布區(qū)主導(dǎo)氣候因子

    基于刀切法(Jackknife)來檢驗?zāi)P瓦\算結(jié)果,反映出不同環(huán)境變量對于分布增益的貢獻大小。當(dāng)前氣候條件下19個環(huán)境變量的貢獻率排列前11(貢獻率>1%)依次為bio16(最濕季度降水量,mm)、bio10(最暖季度平均溫度,℃)、bio6(最冷月份最低溫度,℃)、bio15(降水量季變異系數(shù))、bio9(最干季度平均溫度,℃)、bio11(最冷季度平均溫度,℃)、bio4(溫度季節(jié)性變化標(biāo)準(zhǔn)差)、bio8(最濕季度平均溫度,℃)、bio5(最暖月份最高溫度,℃)、bio14(最干月份降水量,mm)、bio13(最濕月份降水量,mm),累計貢獻率達96.8%(表2)。綜合19個環(huán)境變量的相關(guān)性分析結(jié)果發(fā)現(xiàn),bio16(最濕季度降水量,mm)和bio10(最暖季度平均溫度,℃)是影響當(dāng)代苦參適宜生境分布的主導(dǎo)氣候因子(累計貢獻率達72.4%),其次是bio6(最冷月份最低溫度,℃)、bio15(降水量季變異系數(shù))、bio9(最干季度平均溫度,℃)、bio11(最冷季度平均溫度,℃)、bio4(溫度季節(jié)性變化標(biāo)準(zhǔn)差)、bio8(最濕季度平均溫度,℃)、bio5(最暖月份最高溫度,℃)、bio14(最干月份降水量,mm)、bio13(最濕月份降水量,mm),但貢獻率較前兩個?。ū?)。為進一步明確苦參在當(dāng)前氣候條件下主導(dǎo)氣候因子與存在概率間的關(guān)系,分別將前兩個氣候因子導(dǎo)入MaxEnt模型中單獨建模,繪制單變量響應(yīng)曲線(圖3),隨后計算兩個主導(dǎo)氣候因子的閾值。存在概率大于50%時,最濕季度降水量(bio16)最適宜取值為400~800 mm,最暖季度平均溫度(bio10)值為20~26 ℃。

    2.3 氣候變化條件下苦參在中國的適生度分布

    在氣候變化背景下,利用MaxEnt模型預(yù)測了我國苦參在末次盛冰期、當(dāng)前氣候、未來氣候(2070)三個適宜性生境分布區(qū)域。將建模得到的輸出結(jié)果中的ASCII文件通過ArcGIS中的格式轉(zhuǎn)化工具轉(zhuǎn)化成柵格文件,繪制不同歷史時期的苦參在全國分布的適生度圖(圖4, 圖5, 圖6),圖中的色帶表示適生度,取值由0至1,適生度取值越大表明該區(qū)域苦參適宜分布的概率越大,該區(qū)域即為苦參的潛在分布區(qū)域。

    2.4 氣候變化對我國苦參適宜性生境分布的影響

    在氣候變化背景下,利用MaxEnt模型預(yù)測了我國苦參在末次盛冰期、當(dāng)前氣候、未來氣候(2070)三個適宜性生境分布區(qū)域。將三個歷史時期的適生度(圖4, 圖5, 圖6)通過重分類工具,對圖層進行如下劃分:P<0.33為非適生區(qū),0.33≤P<0.66為低適生區(qū),P≥0.66為高度適生區(qū),得到三個歷史時期的氣候情境下苦參在全國的適宜性生境分布,如(圖7, 圖8, 圖9)所示,綠色為非適生區(qū)、 紫色為低適生區(qū)、 紅色為高適生區(qū), 適生區(qū)

    包括低適生區(qū)與高適生區(qū)。

    根據(jù)模型結(jié)果,從末次盛冰期開始到未來氣候(2070,RCP8.5)三種氣候情境下,苦參在我國的適宜生境面積逐漸減少,呈現(xiàn)遞減的趨勢,且整體有向高緯度地區(qū)移動的趨勢,但變化范圍較小;高適宜生境面積先減少再增加,變化范圍同樣較小。具體情況是(1)末次盛冰期苦參的適宜生境面積占研究區(qū)域總面積的36.20%,高適宜生境面積占研究區(qū)域總面積的4.22%。(2)當(dāng)前氣候條件下苦參的適宜生境面積占研究區(qū)域總面積的35.90%,較末次盛冰期減少0.3%;當(dāng)前氣候條件下苦參的高適宜生境面積占研究區(qū)域總面積的2.43%,較末次盛冰期減少1.78%。(3)未來氣候條件下苦參的適宜生境面積占我國總面積的35.15%,較當(dāng)前氣候減少0.75%;未來氣候條件下苦參的高適宜生境面積占研究區(qū)域總面積的2.87%,較當(dāng)前氣候增加0.44%(圖7,圖8,圖9,表3)。

    2.5 苦參在全國適宜分布區(qū)的面積統(tǒng)計

    在ArcGIS中將三個不同歷史時期氣候情境下的適生度重分類,進行幾何計算,統(tǒng)計各圖層的面積。末次盛冰期、當(dāng)前氣候及未來氣候下的適生總面積分別為348.621 5×105、345.280 1×105、337.867 6 × 105? km2 (表3)。

    3 討論與結(jié)論

    3.1 氣候變化條件下苦參在我國潛在分布區(qū)模擬

    目前對模型精度評價應(yīng)用最廣泛的方法是ROC曲線法(AUC法),AUC值被公認(rèn)為生態(tài)位模型評價領(lǐng)域的最佳評價指標(biāo)。AUC的取值范圍在0.5~1, 值越接近于1,模型精度越高。本研究基于苦參的725個地理分布點和19個氣候因子,采用MaxEnt 模型和地理信息系統(tǒng)(ArcGIS)對苦參在末次盛冰期、當(dāng)前氣候、未來氣候(2070,RCP8.5)三種氣候情景下在我國的適生區(qū)分布格局進行預(yù)測。受試者工作特征曲線顯示,三種情境的氣候條件下,MaxEnt 模型訓(xùn)練集和測試集的AUC平均值均大于0.85,顯著高于隨機預(yù)測的AUC值(0.5)。這說明此次模型預(yù)測的地理分布與苦參實際分布擬合度較高,擬合效果優(yōu)秀,預(yù)測結(jié)果可信度高。

    模型結(jié)果顯示,從末次盛冰期開始到未來氣候(2070,RCP8.5)三種氣候情境下,苦參的適宜生境區(qū)域主要集中分布在我國大興安嶺-張家口-蘭州-拉薩-喜馬拉雅山東南端一線的東南部地區(qū)。末次盛冰期至當(dāng)前氣候情境下,苦參在呼倫貝爾高原和三江平原的適宜生境分布范圍變小,在東北平原和橫斷山脈有向北移動的趨勢。當(dāng)前氣候情境至未來氣候情境下, 苦參在呼倫貝爾高原、東北平原和長白山區(qū)域適宜生境面積向北擴張增加。

    當(dāng)前氣候情境下,高適宜生境區(qū)域主要集中分布在遼東半島、 河北東北部與遼寧西南部的燕山區(qū)域、河北與山西交接的太行山區(qū)域、陜西南部、甘肅西南部、湖北西部、四川和重慶北部等地區(qū);山東中部和膠東半島、江蘇南部、安徽西南部大別山區(qū)以及江西中部等地區(qū)有零星分布。末次盛冰期至當(dāng)前氣候情境下,陜北高原、山西中部的呂梁山區(qū)和東南部的太行山、中條山和膠東半島、遼東半島等地區(qū)的高適宜生境區(qū)域減少。當(dāng)前氣候情境至未來氣候下,燕山和遼東半島區(qū)域的高適宜生境分布區(qū)向北擴展,大巴山區(qū)和鄂西地區(qū)的高適宜生境分布區(qū)域減少,山東東南部、陜西中部、貴州中部、重慶南部、江西北部、浙江和福建沿海地區(qū)以及臺灣北部等地區(qū)將發(fā)展成為新的高適宜生境分布區(qū)域。

    模型的模擬結(jié)果可以指導(dǎo)以后的實際調(diào)研方向和范圍,進一步探索和豐富苦參的現(xiàn)實分布區(qū)域。

    3.2 氣候變化對苦參潛在分布區(qū)域的時空分布格局的影響

    氣候變化是人類當(dāng)前面臨的最為嚴(yán)峻的環(huán)境問題,政府間氣候變化專門委員會第五次評估報告(intergovernmental panel on climate change,IPCC)指出,相比于1850—1900年,21世紀(jì)末全球氣溫增幅可能會超過1.5 ℃,在RCP6.0和RCP8.5情景下,溫度升高可能會超過2 ℃(Alexander, 2013)。氣候變化對植物的生長發(fā)育、地理分布及種群數(shù)量大小等都將會產(chǎn)生極大的影響( Alan et al., 2006)。越來越多的研究發(fā)現(xiàn),氣候變暖可能會導(dǎo)致植物適生區(qū)范圍減?。▌⒀蟮?,2009),使得植物向高海拔、高緯度地區(qū)遷移( Wardle et al., 1992;吳軍等,2011) 。然而,氣候變化對不同植物的影響不一致。因此,研究植物對氣候變化的響應(yīng)對保護生物多樣性具有十分重要的意義。

    本研究在氣候變化背景下,利用MaxEnt模型預(yù)測了苦參在末次盛冰期、當(dāng)前氣候、未來氣候(2070,RCP8.5)三個適宜性生境分布區(qū)域。根據(jù)模型結(jié)果,從末次盛冰期開始到未來氣候(2070,RCP8.5)三種氣候情境下,苦參在我國的適宜生境面積逐漸減少,呈現(xiàn)遞減的趨勢,且整體有向高緯度地區(qū)移動的趨勢。當(dāng)前氣候條件下苦參的適宜生境面積較末次盛冰期減少0.3%,未來氣候條件下苦參的適宜生境面積較當(dāng)前氣候減少0.75%。這可能由于全球氣候變暖,海平面上升,低緯度地區(qū)變熱,高緯度地區(qū)逐漸變暖,使得原有的較低緯度適生區(qū)變成了低適生區(qū)或者非適生區(qū),原有的高緯度非適生區(qū)域變成了適生區(qū)域。全球氣候變暖對苦參潛在分布區(qū)具有強烈的影響,暖期讓苦參實際生態(tài)位變窄,不利于苦參生存。

    氣候變暖情景下,苦參高適生區(qū)面積將減少,說明氣候變暖對苦參的擴張有消極影響,使得非適生區(qū)域增加。在氣候變暖情境下,可能導(dǎo)致北半球中高緯度地區(qū)降水強度增加,同時全國各地區(qū)極端高溫指數(shù)顯著增加,進而使得苦參的非適生區(qū)域增加,高適生區(qū)域減少且向高緯度地區(qū)移動(趙宗慈等,2008)。這符合苦參喜涼、不耐澇的生物學(xué)特性。

    預(yù)計到21世紀(jì)末氣候變暖情景下,中國中西部地區(qū)夏季降水量將顯著減少,這導(dǎo)致苦參原本不適宜生長的區(qū)域更加不適宜其擴張(梁玉蓮和延曉冬,2016)。這與本研究模型運算得到三種氣候情境下苦參非適生區(qū)域面積逐漸增加的結(jié)果一致,也與Wardle等(1992)得到的氣候變暖會使得植物向高海拔、高緯度地區(qū)遷移的研究結(jié)論一致。

    3.3 生物氣候因子與苦參適宜分布格局的關(guān)系

    MaxEnt模型的刀切法檢驗和主導(dǎo)因子響應(yīng)曲線分析表明,苦參的潛在適生區(qū)分布影響較大的氣候因子是最暖季度平均溫度(bio10)、最濕季度降水量(bio16),其最適宜區(qū)間分別為20~26 ℃、400~800 mm??梢姡瑴囟群徒邓渴怯绊懣鄥⒎植嫉闹鲗?dǎo)因子。溫度和降水作為重要的氣候因子顯著影響著植物的生理生態(tài)特征、分布范圍、多樣性和豐富度等。

    苦參莖葉和根系均在7、8月份進入快速生長的階段,在8月下旬種子乳熟,隨后進入種子脫水階段,9月上旬種子完全成熟,莢果開裂,種子掉落,完成生殖生長。這一階段正好與我國大部分地區(qū)最暖季度重合,而20~26 ℃的最暖季度平均溫度有利于苦參種子的乳熟和脫水,從而影響其在地理尺度上的分布格局。因此,在苦參的生產(chǎn)實踐中,在選取作為選育優(yōu)質(zhì)種源的種植基地時要注意當(dāng)?shù)販囟?,?yīng)優(yōu)先選擇最暖季度平均溫度的適宜取值在20~26 ℃之間的地域。

    相對于干旱脅迫,苦參植株受澇害影響更嚴(yán)重,高溫高濕條件下更容易引起苦參根腐病等病害的發(fā)生(張文龍等,2018)。適宜程度的降水可以提供充足的水分供應(yīng),促進植物生長。但水分過多卻嚴(yán)重阻礙植物生長產(chǎn)生,土壤中過量水分導(dǎo)致土壤透氣性變差,影響根系正常呼吸作用,不利于植物生長發(fā)育。同時過多水分引起的高濕度有利于病原體的快速繁殖,形成嚴(yán)重病害,威脅苦參的存活,從而影響其在地理尺度上的分布格局 (Xin et al., 2016)。

    我國疆域遼闊,經(jīng)緯跨度大,地形落差顯著,一些特定的等降水量線反而成為中國自然區(qū)劃的重要參考指標(biāo)。例如:800 mm等降水量線是中國南方濕潤區(qū)和北方半濕潤區(qū)的重要分界線;400 mm等降水量線將中國劃分為東部濕潤半濕潤區(qū)、西部干旱半干旱區(qū)兩部分;而200 mm等降水量線以東為半干旱區(qū),以西為干旱區(qū)(車彥軍等,2020)。我國降水的時空分布極不均勻,從時間分布上看,夏季降水較多,我國大部分地區(qū)夏季降水占全年降水量在40%以上,春、秋季次之,冬季最少;從空間上看東南沿海地區(qū)降水豐沛,西北內(nèi)陸地區(qū)降水稀少(肖潺等,2015)。

    由模型輸出結(jié)果可知,三個歷史時期氣候情境下苦參在我國高適生區(qū)主要分布在我國800 mm等降水量線與400 mm等降水量線之間。未來氣候條件下苦參的適宜生境面積較當(dāng)前氣候減少0.75%,這與未來排放情景下,我國400 mm等降水量線將發(fā)生西移和北抬的研究結(jié)果相一致(高艷紅等,2020)。因此,在苦參的生產(chǎn)實踐中,在選取高品質(zhì)藥材的種植基地時要注意當(dāng)?shù)亟邓?,尤其要重點關(guān)注最濕季度降水量的適宜取值在400~800 mm的地域。

    但是,物種存在概率對環(huán)境變量的反饋曲線反映的是單一環(huán)境變量的作用,苦參生命活動受到各種環(huán)境變量(包括氣候因子、天敵種類、植被條件等)的綜合影響,因此,此結(jié)果可作為判斷苦參與環(huán)境變量之間關(guān)系的參考。

    3.4 苦參在我國范圍內(nèi)的保護和開發(fā)利用的建議

    研究植物分布格局對氣候變化的響應(yīng), 了解植物的氣候需求及其與物種地理分布間的關(guān)系, 對于揭示物種的形成、遷移和擴散歷史, 提出合理的物種多樣性保護策略及合理利用、生態(tài)恢復(fù)與建設(shè)具有重要意義(王綺等,2020)??茖W(xué)預(yù)測藥用植物適宜分布區(qū),并規(guī)劃建立野生種質(zhì)資源庫、自然保護區(qū),已成為保護藥用植物資源的一種有效保護途徑(Yan et al., 2020)。

    本研究通過模型預(yù)測發(fā)現(xiàn),在當(dāng)前氣候條件下苦參在全國高適宜生境區(qū)域主要集中分布在遼東半島、河北東北部與遼寧西南部的燕山區(qū)域、河北與山西交接的太行山區(qū)域、陜西南部、甘肅西南部、湖北西部、四川和重慶北部等地區(qū);山東中部和膠東半島、江蘇南部、安徽西南部大別山區(qū)以及江西中部等地區(qū)有零星分布。這些地區(qū)為苦參的生態(tài)適宜分布區(qū),也是野生苦參資源的重點分布區(qū)域,是未來苦參重點保育區(qū)。

    氣候變暖情景下,當(dāng)前氣候條件下苦參在我國的適宜生境面積較末次盛冰期減少0.3%,未來氣候條件下苦參的適宜生境面積較當(dāng)前氣候減少0.75%。未來氣候條件下,受氣候變化影響,四川、重慶、貴州、湖北和江西將成為苦參適宜分布的高風(fēng)險地區(qū),這些地區(qū)的野生苦參資源在未來要加強保護或遷地保育等。

    氣候變化情景下,苦參在我國相對穩(wěn)定的高適生區(qū)主要分布在河北西部和東北部山區(qū)、內(nèi)蒙古西南部、遼寧南部等地區(qū),在這些地區(qū)合理優(yōu)化布局等措施,進一步加強保護力度,建立自然保護區(qū),禁止濫采濫挖,同時,建立廣泛的苦參種植和生產(chǎn)加工等基地。

    本研究首次在全國范圍內(nèi)對苦參進行了潛在分布區(qū)和適宜生境進行了預(yù)測和分析,得出全球氣候變暖對苦參潛在分布區(qū)具有一定負(fù)面影響的結(jié)論,即全球氣候變暖讓苦參適生生境縮減,實際生態(tài)位變窄,一定程度上不利于苦參的生長?;诒狙芯拷Y(jié)果,在苦參的生產(chǎn)實踐中,應(yīng)綜合考慮溫度和降水這兩個生態(tài)因子,即在選取作為選育優(yōu)質(zhì)種源的種植基地時要注意當(dāng)?shù)販囟?,?yīng)優(yōu)先選擇最暖季度平均溫度的適宜取值在20~26 ℃之間的地域;在選取高品質(zhì)藥材的種植基地時要注意當(dāng)?shù)亟邓浚绕湟攸c關(guān)注最濕季度降水量的適宜取值在400~800 mm的地域。

    參考文獻:

    ALAN PJ, BUSTAMANTE MR, COLOMA LA, et al., 2006. Widespread amphibian extinctions from epidemic disease driven by global warming [J]. Nature, 439: 161-167.

    ALEXANDER L, ALLEN S, LEEBINDOFF N, 2013. Climate Change 2013: The Physical Science Basis [M]//STOCKER TF, QIN D,? PLATINER GK, et al. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom: Cambridge University Press: 5-6.

    BORAL D, MOKTAN S, 2021. Predictive distributionmodeling of Swertia bimaculata in Darjeeling-Sikkim Eastern Himalaya using MaxEnt: current and future scenarios[J]. Ecol Proc, 10(1): 1-16.

    Editorial Committee of Flora of China, Chinese Academy of Sciences, 1994. Flora Republicae Popularis Sinicae (Vol. 40) [M]. Beijing: Science Press: 81-82. [中國科學(xué)院中國植物志編輯委員會, 1994. 中國植物志(第四十卷)[M]. 北京: 科學(xué)出版社: 81-82.]

    ELITH J, GRAHAM CH, ANDERSON RP, et al. , 2006. Novel methods improve prediction of species’ distributions from occurrence data[J]. Ecography, 29(2): 119-151.

    GAO YH, XU JW, ZHANG M, et al., 2020. Advances of the 400 mm isohyet and wetness and dryness changes in the mainland of China[J]. Adv Earth Sci, 35(11): 1101-1112. [高艷紅, 許建偉, 張萌, 等, 2020. 中國400 mm 等降水量線和干濕變化研究進展[J]. 地球科學(xué)進展. 35(11): 1101-1112.]

    GUO LP, HUANG LQ, JIANG YX, 2005. Application of “3S” technologies in continuous usage of traditional Chinese medicine resources [J]. Chin J Chin Mat Med, 30(18): 1397-1400. [郭蘭萍, 黃璐琦, 蔣有緒, 2005. “3S”技術(shù)在中藥資源可持續(xù)利用中的應(yīng)用[J]. 中國中藥雜志, 30(18): 1397-1400.]

    HUANG YJ, LU JB, WANG FT, 2017. Predicting the potential geographical distribution of Hainan Odorous Frog (Odorrana hainanensis) in Hainan Province by MaxEnt[J]. Chin J Zool, 52(1): 30-41. [黃勇杰, 盧佳斌, 王鋒堂, 等, 2017. 基于MaxEnt模型預(yù)測海南島海南臭蛙的潛在地理分布[J]. 動物學(xué)雜志, 52(1): 30-41.]

    JHA KK, JHA R, 2021. Study of vulture habitat suitability and impact of climate change incentral India Using MaxEnt[J]. J Resour Ecol, 12(1), 30-42.

    JI W, GAO G, WEI JF, 2021. Potential global distribution of Daktulosphaira vitifoliae under climate change based on MaxEnt [J]. Insect, 12(4): 347.

    JI Y, DU YB, WANG B, et al., 2011. Study on growth dynamics of Sophora flavescens from different germplasm sources [J]. Pratac Sci, 28 (5): 802-806 [紀(jì)瑛, 杜彥斌, 王斌, 等, 2011. 不同來源地苦參生長動態(tài)的研究[J]. 草業(yè)科學(xué), 28(5): 802-806.]

    JU L, CAO B, BAI CK, et al., 2014. Predictive distribution and habitat suitability assessment of Notholirion bulbuliferum based on MaxEnt and ArcGIS [J]. Chin J Ecol, 33(6): 1623-1628. [車樂, 曹博, 白成科, 等, 2014. 基于MaxEnt和ArcGIS對太白米的潛在分布預(yù)測及適宜性評價[J]. 生態(tài)學(xué)雜志, 33(6): 1623-1628.]

    JU YJ, GUAN XC, WANG SJ, et al., 2020. Spatial analysis of annual precipitation lines of 800 mm in the eastern monsoon of China[J]. Plateau Meteorl, 39(5): 997-1006. [車彥軍, 管小春, 王世金, 等, 2020. 中國東部季風(fēng)區(qū)800 mm等降水量線的空間變化分析[J]. 高原氣象, 39(5): 997-1006.]

    KAMEL M, BREAM AS, MOURSY MM, et al., 2021. Predicting the geographic distribution habitats of Schizomyia buboniae (Diptera: Cecidomyiidae) and its host plant Deverra tortuosa (Apiaceae) in Egypt by using MaxEnt modeling[J]. J Basic Appl Zool, 82(1): 1-13.

    KIM Z, SHIM T, KOO YM, et al., 2020. Predicting the impact of climate change on freshwater fish distribution by incorporating water flow rate and quality variables[J]. Sustainability, 12(23): 10001.

    LIANG YL, YAN XD, 2016. Prediction of climate change over China and uncertainty analysis during the 21st Century under RCPs [J]. J Tropical Meteorol, 32(2): 183-192. [梁玉蓮, 延曉冬, 2016. RCPs情景下中國21世紀(jì)氣候變化預(yù)估及不確定性分析[J]. 熱帶氣象學(xué)報, 32(2): 183-192.]

    LIU L, GUAN LL, ZHAO HX, et al., 2021.Modeling habitat suitability of Houttuynia cordata Thunb (Ceercao) using MaxEnt under climate change in China[J]. Ecol Inform, 2021, 63: 101324.

    LIU YM, ZHOU SD, XIE DF, et al., 2018. Potential distribution of? Fritillaria unibracteata predicted by the MaxEnt model[J]. Guihaia, 38(3): 352-360. [劉艷梅, 周頌東, 謝登峰, 等, 2018. 基于最大熵模型(MaxEnt)預(yù)測暗紫貝母的潛在分布[J]. 廣西植物, 38(3): 352-360.]

    LIU Y,ZHANG J, YANG WQ, 2009. Responses of alpine biodiversity to climate change[J]. Biodivers Sci, 17(1): 88-96.[劉洋, 張健, 楊萬勤, 2009. 高山生物多樣性對氣候變化響應(yīng)的研究進展[J]. 生物多樣性, 17(1): 88-96.]

    LU WH, 1959. Study on materia medica of? Sophora flavescens [J]. J Chengdu Coll Traditional Chin Med, (3): 49-53. [陸聞鴻, 1959. 苦參的本草學(xué)研究[J]. 成都中醫(yī)學(xué)院學(xué)報, (3): 49-53.]

    LU YY, GUO S, YAN H, et al., 2018. Ecological suitability study on Euphorbia kansui[J]. Mod Chin Med, 20(12): 1471-1475. [盧有媛, 郭盛, 嚴(yán)輝, 等, 2018. 甘遂生態(tài)適宜性區(qū)劃研究[J]. 中國現(xiàn)代中藥, 20(12): 1471-1475.]

    National Pharmacopoeia Commission, 2020. Pharmacopoeia of the People’s Republic of China (Volume I, 2020 Edition) [M]. Beijing: China Pharmaceutical Science and Technology Press: 211-212. [國家藥典委員會, 2020. 中華人民共和國藥典(2020年版一部)[M]. 北京: 中國醫(yī)藥科技出版社: 211-212.]

    PAN SY, ZHU ZH, YAO TH, et al., 2016. Predicting suitable habitats of medicinal plant Fallopia multiflora in China under climate change scenarios[J]. J NW A & F Univ (Nat Sci Ed), 44(1): 192-198. [潘石玉, 朱志紅, 姚天華, 等, 2016. 氣候變化背景下藥用植物何首烏在中國適生區(qū)分布預(yù)測[J]. 西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版), 44(1): 192-198.]

    PEI ST, ZHANG H, ZHANG PF, et al., 2021. Study on ecological suitability of Euphorbiae Pekinensis Rupr. [J]. Chin J Inf Tradit Chin Med, 28(5): 1-4. [裴蘇婷, 張虹, 章鵬飛, 等, 2021. 京大戟生態(tài)適宜性區(qū)劃研究[J]. 中國中醫(yī)藥信息雜志, 28(5): 1-4.]

    QIAO HJ, HU JH, HUANG JH, 2013. Theoretical basis, future directions, and challenges for ecological nichemodels [J]. Sci Sin (Vitae), 43(11): 915-927. [喬慧捷, 胡軍華, 黃繼紅, 2013. 生態(tài)位模型的理論基礎(chǔ)、發(fā)展方向與挑戰(zhàn)[J]. 中國科學(xué): 生命科學(xué), 43(11): 915-927.]

    ROZHNOV VV, PSHEGUSOV RH, HERNANDEZ-BIANLO JA, et al., 2021. MaxEnt modeling for predicting suitable habitats in the North Caucasus (Russian Part) for Persian Leopard (P. p. ciscaucasica) based on GPS data from collared and released animals[J]. Izv, Atmos Oceanic Phys, 56(9): 1090-1106.

    SAHA A, RAHMAN S, ALAM S, 2021. Modeling current and future potential distributions of desert locust Schistocerca gregaria (Forskl) under climate change scenarios using Maxent. [J]. J Asia-Pac Biodivers (prepublish), 14: 399-409.

    TAN YH, ZHANG XJ, YUAN SS, et al., 2020. Prediction of the ecological suitability of Gentiana macrophylla Pall. under scenarios of global climate change [J]. Chin J Ecol, 39(11): 3766-3773. [檀逸虹, 張喜娟, 原樹生, 等, 2020. 全球氣候變化背景下秦艽生態(tài)適宜性預(yù)測[J]. 生態(tài)學(xué)雜志, 39(11): 3766-3773.]

    WEN P, PENG K, CHEN X, et al., 2021. Habitat suitability evaluation of four ungulate species in BaishuiheNational Nature Reserve based on maxent model[J]. J Sichuan For Sci Technol, 42(1): 70-75. [溫平, 彭科, 陳旭, 等, 2021. 基于MaxEnt模型的白水河國家級保護區(qū)有蹄類保護動物生境適宜性評價[J]. 四川林業(yè)科技, 42(1): 70-75.]

    WANG Q, FAN BG, ZHAO GH, 2020. Prediction of potential distribution area of Corylus mandshurica in China under climate change[J]. Chin J Ecol, 39(11): 3774-3784. [王綺, 樊保國, 趙光華, 2020. 氣候變化下毛榛在中國的潛在適生區(qū)預(yù)測[J]. 生態(tài)學(xué)雜志, 39(11): 3774-3784.]

    WANG YS, XIE BY, WAN FH, et al., 2007. Application of ROC curve analysis in evaluating the performance of alien species’ potential distribution models[J]. Biodivers Sci, 15(4): 365-372. [王運生, 謝丙炎, 萬方浩, 等, 2007. ROC曲線分析在評價入侵物種分布模型中的應(yīng)用[J]. 生物多樣性, 15(4): 365-372.]

    WARDLE P, COLEMAN MC, 1992. Evidence for rising upper limits of four native New Zealand forest trees [J]. New Zealand J Bot, 30: 303-314.

    WU J, XU HG, CHEN L, 2011. A review of impacts of climate change on species[J]. J Ecol Rural Environ, 27(4): 1-6. [吳軍, 徐海根, 陳煉, 2011. 氣候變化對物種影響研究綜述[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 27(4): 1-6.]

    XIAO C, YU RC, YUAN WH, et al., 2015. Spatial and temporal differences characteristics of the rainy season of mainland China [J]. Acta Meteorol Sin, 73 (1): 84-92. [肖潺, 宇如聰, 原韋華, 等, 2015. 中國大陸雨季時空差異特征分析[J]. 氣象學(xué)報, 73(1): 84-92.]

    XIN XF, NOMUR K, AUNG K, et al., 2016. Bacteria establish an aqueous living space in plants crucial for virulence[J]. Nature, 539, 524-529.

    XING DL, HAO ZQ, 2011. The principle of maximum entropy and its applications in ecology[J]. Biodivers Sci, 19(3): 295-302. [邢丁亮, 郝占慶, 2011. 最大熵原理及其在生態(tài)學(xué)研究中的應(yīng)用[J]. 生物多樣性, 19(3): 295-302.]

    YAN HY, FENG L, ZHAO YF, et al., 2020. Prediction of the spatial distribution of Alternanthera philoxeroides in China based on ArcGIS and MaxEnt[J]. Glob Ecol Conservation, 21: 1-8.

    ZENG JR, LI CM, LIU JZ, et al., 2021. Ecological assessment of current and future Pogostemon cablin Benth. potential planting regions in China based on maxent and ArcGIS models [J]. J Appl Res Med Aromat Plants (prepublish), 24(9): 100308.

    ZHANG WL, DU FQ, HUANG XL, et al., 2018. The outbreak regularity and control methods of the main diseases and insect pests for artificially planted Sophora flavescens [J]. Lishizhen Med Mat Med Res, 29 (9): 2241-2243 [張文龍, 杜富強, 黃旭龍, 等, 2018. 貴州特色藥材苦參人工栽培常見病蟲害發(fā)生規(guī)律及防治措施[J]. 時珍國醫(yī)國藥, 29(9): 2241-2243.]

    ZHAO ZC, LUO Y, JIANG Y, et al., 2008. Assessment and prediction of precipitation and droughts/floods changes over the world and in China[J]. Sci Technol Rev, (6): 28-33. [趙宗慈, 羅勇, 江瀅, 等, 2008. 全球和中國降水、旱澇變化的檢測評估[J]. 科技導(dǎo)報, (6): 28-33.]

    ZHONG GS, 2016. Traditional Chinese medicine [M]. 4th ed. Beijing: China Traditional Chinese Medicine Press: 98 . [鐘贛生, 2016. 中藥學(xué) [M]. 4版. 北京: 中國中醫(yī)藥出版社: 98.]

    ZHU GP, LIU GQ, BU WJ, et al., 2013. Ecological niche modeling and its applications in biodiversity conservation [J]. Biodivers Sci, 21(1): 90-98. [朱耿平, 劉國卿, 卜文俊, 等, 2013. 生態(tài)位模型的基本原理及其在生物多樣性保護中的應(yīng)用[J]. 生物多樣性, 21(1): 90-98.]

    (責(zé)任編輯 蔣巧媛)

    收稿日期:2021-09-09

    基金項目:國家自然科學(xué)基金(41961009)[Supported by National Natural Science Foundation of China(41961009)]。

    第一作者: 張濤(1988-),碩士,助理研究員,主要從事藥用植物種質(zhì)資源保護與開發(fā)利用研究,(E-mail)zhangtao361@126.com。

    通信作者:劉毅,博士,副研究員,主要從事藥用與功能植物學(xué)研究,(E-mail)yiliu609@outlook.com。

    猜你喜歡
    苦參氣候變化
    北極熊如何適應(yīng)氣候變化
    英語世界(2023年10期)2023-11-17 09:19:12
    Efficacy of Kushen decoction (苦參湯) on high-fat-diet-induced hyperlipidemia in rats
    以苦參為主治療心律失常的療效觀察
    《應(yīng)對氣候變化報告(2022)》發(fā)布
    苦參百部液治蟲咬皮炎
    保護生物多樣性 積極應(yīng)對氣候變化
    大劑苦參治不寐
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    苦參凝膠用于宮頸 LEEP 術(shù)后180例臨床觀察
    黄片播放在线免费| 精品人妻1区二区| 日韩人妻精品一区2区三区| 乱人伦中国视频| 国产成人免费无遮挡视频| 亚洲熟女精品中文字幕| 日韩 欧美 亚洲 中文字幕| 久热爱精品视频在线9| 操美女的视频在线观看| 51午夜福利影视在线观看| 女人高潮潮喷娇喘18禁视频| 国产成人av激情在线播放| 成人黄色视频免费在线看| 国产又色又爽无遮挡免| 老司机亚洲免费影院| 久久精品久久久久久噜噜老黄| 我的亚洲天堂| 成人手机av| 人人妻人人爽人人添夜夜欢视频| 日日夜夜操网爽| 国产97色在线日韩免费| 国产亚洲欧美在线一区二区| 超碰97精品在线观看| 人妻人人澡人人爽人人| 人人妻人人澡人人看| 女人久久www免费人成看片| 国产深夜福利视频在线观看| 久久午夜综合久久蜜桃| 麻豆国产av国片精品| 黄频高清免费视频| 国产精品国产av在线观看| 男人爽女人下面视频在线观看| 国产亚洲欧美精品永久| 国产熟女欧美一区二区| 日韩精品免费视频一区二区三区| 建设人人有责人人尽责人人享有的| 成年女人毛片免费观看观看9 | 精品第一国产精品| 中文字幕人妻丝袜制服| 精品一区二区三卡| 亚洲国产毛片av蜜桃av| 丝瓜视频免费看黄片| 伊人久久大香线蕉亚洲五| 国产野战对白在线观看| 欧美 日韩 精品 国产| 精品少妇一区二区三区视频日本电影| 中文字幕高清在线视频| 啦啦啦在线观看免费高清www| 亚洲欧洲国产日韩| 国产成人影院久久av| 亚洲欧美精品自产自拍| 国产成人精品久久二区二区91| 国产成人精品在线电影| 大香蕉久久网| 看十八女毛片水多多多| 亚洲人成77777在线视频| 成在线人永久免费视频| 在线观看免费午夜福利视频| 免费在线观看完整版高清| 国产男女超爽视频在线观看| 国产福利在线免费观看视频| 免费在线观看黄色视频的| 国产视频首页在线观看| 亚洲自偷自拍图片 自拍| 亚洲国产欧美日韩在线播放| 久久精品熟女亚洲av麻豆精品| 久久精品人人爽人人爽视色| 欧美在线一区亚洲| 高清黄色对白视频在线免费看| 免费观看av网站的网址| 亚洲av日韩在线播放| 国产精品 国内视频| 18在线观看网站| 免费观看a级毛片全部| 亚洲欧美成人综合另类久久久| 操出白浆在线播放| 亚洲国产av新网站| 国产一区亚洲一区在线观看| 亚洲国产中文字幕在线视频| av一本久久久久| 欧美日韩视频高清一区二区三区二| 精品第一国产精品| 久久青草综合色| 色播在线永久视频| 午夜91福利影院| 国产精品熟女久久久久浪| 高清欧美精品videossex| 亚洲欧美清纯卡通| 老司机亚洲免费影院| 一级毛片女人18水好多 | 久久久久国产精品人妻一区二区| 韩国高清视频一区二区三区| 蜜桃国产av成人99| 亚洲专区国产一区二区| 国产高清国产精品国产三级| 性少妇av在线| 在线观看免费日韩欧美大片| 亚洲色图 男人天堂 中文字幕| 国产精品秋霞免费鲁丝片| 自拍欧美九色日韩亚洲蝌蚪91| 777米奇影视久久| √禁漫天堂资源中文www| 美女扒开内裤让男人捅视频| 中文字幕最新亚洲高清| 黄色毛片三级朝国网站| 热99国产精品久久久久久7| 日韩中文字幕欧美一区二区 | 国产亚洲精品久久久久5区| 免费女性裸体啪啪无遮挡网站| 色精品久久人妻99蜜桃| 日本午夜av视频| 国产精品秋霞免费鲁丝片| e午夜精品久久久久久久| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 亚洲中文日韩欧美视频| 不卡av一区二区三区| 一级黄片播放器| 飞空精品影院首页| 晚上一个人看的免费电影| 无遮挡黄片免费观看| 青春草视频在线免费观看| 国产成人系列免费观看| 美女午夜性视频免费| 日韩,欧美,国产一区二区三区| 视频区欧美日本亚洲| www.999成人在线观看| 亚洲av美国av| www.999成人在线观看| 岛国毛片在线播放| 亚洲欧美精品自产自拍| 少妇人妻 视频| 午夜老司机福利片| 七月丁香在线播放| 中文精品一卡2卡3卡4更新| 91麻豆av在线| av国产精品久久久久影院| 亚洲av片天天在线观看| 亚洲欧美一区二区三区国产| 婷婷成人精品国产| 无遮挡黄片免费观看| 蜜桃在线观看..| 精品久久久久久久毛片微露脸 | 一区二区三区激情视频| 99久久人妻综合| 久久久欧美国产精品| 黄频高清免费视频| 欧美人与性动交α欧美精品济南到| 波多野结衣一区麻豆| 久久天躁狠狠躁夜夜2o2o | 黑人欧美特级aaaaaa片| 菩萨蛮人人尽说江南好唐韦庄| 久久这里只有精品19| 老司机午夜十八禁免费视频| 99国产精品免费福利视频| 老汉色∧v一级毛片| 男人舔女人的私密视频| 国产真人三级小视频在线观看| 1024香蕉在线观看| 美女脱内裤让男人舔精品视频| 最新的欧美精品一区二区| 欧美日韩一级在线毛片| 一本—道久久a久久精品蜜桃钙片| 男女高潮啪啪啪动态图| 日日摸夜夜添夜夜爱| 婷婷色综合www| 下体分泌物呈黄色| 三上悠亚av全集在线观看| av在线app专区| 男女无遮挡免费网站观看| videosex国产| 亚洲激情五月婷婷啪啪| 一本色道久久久久久精品综合| 亚洲欧洲精品一区二区精品久久久| 欧美 亚洲 国产 日韩一| av线在线观看网站| 欧美黑人欧美精品刺激| 永久免费av网站大全| 精品国产一区二区久久| 午夜免费鲁丝| 中文字幕精品免费在线观看视频| 国产精品一区二区精品视频观看| 少妇粗大呻吟视频| 熟女av电影| 免费女性裸体啪啪无遮挡网站| 日本色播在线视频| www.熟女人妻精品国产| 一区二区三区乱码不卡18| 欧美成人午夜精品| 国产免费又黄又爽又色| 女人久久www免费人成看片| 青草久久国产| 欧美精品av麻豆av| 免费在线观看视频国产中文字幕亚洲 | 日韩熟女老妇一区二区性免费视频| 国产伦人伦偷精品视频| 成人黄色视频免费在线看| 少妇的丰满在线观看| 国产av国产精品国产| 国产老妇伦熟女老妇高清| 欧美+亚洲+日韩+国产| 久久毛片免费看一区二区三区| 成年人免费黄色播放视频| 天天添夜夜摸| 精品国产乱码久久久久久小说| 国产男女内射视频| 婷婷色综合www| 亚洲av在线观看美女高潮| 免费在线观看完整版高清| 亚洲欧美精品自产自拍| 国产成人av教育| 精品少妇久久久久久888优播| 亚洲国产精品成人久久小说| 国产野战对白在线观看| 一二三四社区在线视频社区8| 亚洲一区二区三区欧美精品| 国产97色在线日韩免费| 大片免费播放器 马上看| 亚洲av综合色区一区| 热99久久久久精品小说推荐| 国产精品 国内视频| 大片电影免费在线观看免费| 色播在线永久视频| 自线自在国产av| 国产一区二区激情短视频 | 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品一区二区在线不卡| 欧美人与性动交α欧美软件| 国产日韩欧美亚洲二区| 欧美激情极品国产一区二区三区| 亚洲人成电影观看| www.熟女人妻精品国产| 久久精品亚洲熟妇少妇任你| 亚洲三区欧美一区| 精品一区二区三卡| 成年美女黄网站色视频大全免费| 国产精品 国内视频| 国产成人av教育| 免费在线观看日本一区| 九色亚洲精品在线播放| 久久亚洲精品不卡| 免费观看av网站的网址| 亚洲欧洲国产日韩| 脱女人内裤的视频| 97人妻天天添夜夜摸| 在线观看人妻少妇| 欧美激情极品国产一区二区三区| 亚洲五月婷婷丁香| 又粗又硬又长又爽又黄的视频| 一级毛片 在线播放| 欧美黑人精品巨大| 国产熟女欧美一区二区| 可以免费在线观看a视频的电影网站| 日本一区二区免费在线视频| 免费观看人在逋| 亚洲一区二区三区欧美精品| 国产精品人妻久久久影院| 99国产精品99久久久久| 久久精品久久久久久久性| 美女午夜性视频免费| 啦啦啦 在线观看视频| 国产黄色免费在线视频| 色综合欧美亚洲国产小说| 男男h啪啪无遮挡| 亚洲熟女毛片儿| 亚洲熟女精品中文字幕| 狂野欧美激情性xxxx| 99国产精品99久久久久| 国产成人欧美| 精品少妇黑人巨大在线播放| 亚洲av欧美aⅴ国产| 欧美日韩一级在线毛片| 成人影院久久| 亚洲欧美中文字幕日韩二区| 国产爽快片一区二区三区| 深夜精品福利| cao死你这个sao货| 桃花免费在线播放| 精品一区二区三区av网在线观看 | 日韩欧美一区视频在线观看| 国产精品一区二区在线观看99| 婷婷色综合大香蕉| 夫妻午夜视频| 在现免费观看毛片| 国产日韩欧美亚洲二区| 精品国产一区二区三区久久久樱花| 观看av在线不卡| 国产成人av激情在线播放| 久热爱精品视频在线9| 在线观看免费视频网站a站| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 久久国产亚洲av麻豆专区| 啦啦啦啦在线视频资源| 男女床上黄色一级片免费看| 后天国语完整版免费观看| 成人18禁高潮啪啪吃奶动态图| 亚洲欧美清纯卡通| 黑人巨大精品欧美一区二区蜜桃| 高清视频免费观看一区二区| 18在线观看网站| 在线观看国产h片| 999精品在线视频| 国产欧美亚洲国产| 一区二区三区精品91| www.精华液| 在现免费观看毛片| av欧美777| 老鸭窝网址在线观看| 国产亚洲午夜精品一区二区久久| 亚洲国产最新在线播放| 免费黄频网站在线观看国产| 国产精品av久久久久免费| 超色免费av| xxxhd国产人妻xxx| 日本五十路高清| 国产成人精品久久二区二区91| svipshipincom国产片| 少妇 在线观看| 国产主播在线观看一区二区 | 国产成人a∨麻豆精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲成色77777| 亚洲一区二区三区欧美精品| 精品熟女少妇八av免费久了| 在线观看免费高清a一片| 不卡av一区二区三区| 国产在线一区二区三区精| 五月开心婷婷网| 最新在线观看一区二区三区 | 日本五十路高清| 国产亚洲欧美在线一区二区| 少妇人妻久久综合中文| netflix在线观看网站| 韩国精品一区二区三区| 亚洲国产精品一区二区三区在线| 国产激情久久老熟女| 天堂8中文在线网| 999久久久国产精品视频| 国产免费福利视频在线观看| 亚洲精品乱久久久久久| 午夜日韩欧美国产| 97在线人人人人妻| 国产一区二区激情短视频 | 极品人妻少妇av视频| 曰老女人黄片| 在线精品无人区一区二区三| 精品一区二区三区四区五区乱码 | 99香蕉大伊视频| 久久国产精品人妻蜜桃| 777久久人妻少妇嫩草av网站| 男人添女人高潮全过程视频| 久久精品国产亚洲av高清一级| av又黄又爽大尺度在线免费看| 亚洲午夜精品一区,二区,三区| 午夜免费观看性视频| 一区二区三区精品91| 精品高清国产在线一区| 丁香六月天网| 又大又爽又粗| 亚洲精品久久午夜乱码| 精品亚洲成a人片在线观看| 国产av精品麻豆| 大片电影免费在线观看免费| 久久久久久免费高清国产稀缺| 捣出白浆h1v1| 少妇人妻久久综合中文| 下体分泌物呈黄色| 啦啦啦 在线观看视频| 少妇人妻久久综合中文| 国产片特级美女逼逼视频| 欧美黄色片欧美黄色片| 国产熟女午夜一区二区三区| 99国产精品99久久久久| 国产真人三级小视频在线观看| 男女边摸边吃奶| 色视频在线一区二区三区| 欧美成人午夜精品| 看十八女毛片水多多多| 啦啦啦在线观看免费高清www| 亚洲国产精品成人久久小说| 赤兔流量卡办理| 亚洲欧美成人综合另类久久久| 欧美日韩av久久| 日韩,欧美,国产一区二区三区| 亚洲国产av新网站| 亚洲精品日本国产第一区| 美女午夜性视频免费| 久久精品亚洲av国产电影网| 国产99久久九九免费精品| 波多野结衣av一区二区av| 精品福利永久在线观看| 五月天丁香电影| 亚洲成人国产一区在线观看 | 国产精品久久久久成人av| 久热这里只有精品99| 免费观看a级毛片全部| 99国产精品一区二区三区| 国产午夜精品一二区理论片| 99热国产这里只有精品6| 老司机深夜福利视频在线观看 | 人妻一区二区av| 久久久精品国产亚洲av高清涩受| 女性被躁到高潮视频| 欧美黄色片欧美黄色片| 国产成人欧美在线观看 | av线在线观看网站| 纵有疾风起免费观看全集完整版| 国产成人啪精品午夜网站| 999精品在线视频| 亚洲欧美一区二区三区黑人| 午夜免费鲁丝| 巨乳人妻的诱惑在线观看| 99香蕉大伊视频| 亚洲午夜精品一区,二区,三区| 国产在线免费精品| 又粗又硬又长又爽又黄的视频| 韩国高清视频一区二区三区| 性高湖久久久久久久久免费观看| 亚洲精品国产一区二区精华液| 建设人人有责人人尽责人人享有的| 国产av精品麻豆| 91老司机精品| 亚洲av成人不卡在线观看播放网 | 成年美女黄网站色视频大全免费| 久久午夜综合久久蜜桃| 97在线人人人人妻| 18在线观看网站| 国产精品久久久久成人av| 一二三四社区在线视频社区8| 高清av免费在线| 成年av动漫网址| 汤姆久久久久久久影院中文字幕| 亚洲国产最新在线播放| 欧美变态另类bdsm刘玥| 欧美日韩福利视频一区二区| 曰老女人黄片| 90打野战视频偷拍视频| 久久ye,这里只有精品| 精品人妻一区二区三区麻豆| 又大又爽又粗| 一级毛片女人18水好多 | 少妇人妻 视频| 秋霞在线观看毛片| 一级毛片我不卡| 性少妇av在线| 老汉色∧v一级毛片| 国产成人精品久久久久久| 好男人视频免费观看在线| 一本一本久久a久久精品综合妖精| 国产一区亚洲一区在线观看| 欧美日韩亚洲高清精品| 电影成人av| 午夜91福利影院| 国产精品二区激情视频| 自拍欧美九色日韩亚洲蝌蚪91| 91字幕亚洲| 国产色视频综合| 午夜激情av网站| 色网站视频免费| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 亚洲免费av在线视频| 晚上一个人看的免费电影| 久久99精品国语久久久| 侵犯人妻中文字幕一二三四区| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品美女久久久久99蜜臀 | 国产精品香港三级国产av潘金莲 | 免费黄频网站在线观看国产| 亚洲av成人精品一二三区| 狠狠精品人妻久久久久久综合| 亚洲精品自拍成人| 丰满少妇做爰视频| 亚洲精品美女久久久久99蜜臀 | 免费少妇av软件| 最黄视频免费看| 在线观看www视频免费| 日韩一本色道免费dvd| 不卡av一区二区三区| 性高湖久久久久久久久免费观看| 久久久久久久久免费视频了| 涩涩av久久男人的天堂| 在线观看免费午夜福利视频| 在线观看www视频免费| 欧美精品一区二区大全| 中文字幕精品免费在线观看视频| 久久久久精品人妻al黑| 国产精品免费大片| 国产在线一区二区三区精| 中文字幕人妻熟女乱码| 在线观看www视频免费| 中国美女看黄片| 在现免费观看毛片| 亚洲成色77777| 人人妻人人添人人爽欧美一区卜| 国产欧美日韩一区二区三 | cao死你这个sao货| 欧美日韩国产mv在线观看视频| 一二三四社区在线视频社区8| 国产精品一区二区免费欧美 | 亚洲久久久国产精品| 青春草亚洲视频在线观看| 国产淫语在线视频| 各种免费的搞黄视频| av天堂在线播放| 久9热在线精品视频| 国产视频一区二区在线看| 国产女主播在线喷水免费视频网站| netflix在线观看网站| 岛国毛片在线播放| 看免费av毛片| 十八禁网站网址无遮挡| 欧美精品av麻豆av| 在线观看免费视频网站a站| 亚洲成国产人片在线观看| 丝瓜视频免费看黄片| 亚洲人成电影免费在线| 亚洲精品中文字幕在线视频| 久久九九热精品免费| 别揉我奶头~嗯~啊~动态视频 | 在现免费观看毛片| 精品久久久久久久毛片微露脸 | 免费不卡黄色视频| 丁香六月天网| 国产成人欧美在线观看 | 国产av一区二区精品久久| 欧美黄色淫秽网站| 国产精品久久久久久精品古装| 国产精品偷伦视频观看了| 亚洲精品第二区| 亚洲一区二区三区欧美精品| 老司机影院成人| 丝袜美腿诱惑在线| 真人做人爱边吃奶动态| 亚洲精品美女久久av网站| 制服人妻中文乱码| 亚洲视频免费观看视频| 男女国产视频网站| av视频免费观看在线观看| 国产视频首页在线观看| 亚洲男人天堂网一区| 久久久久久久久免费视频了| 国产一区二区 视频在线| 男女无遮挡免费网站观看| 熟女少妇亚洲综合色aaa.| 欧美97在线视频| 十分钟在线观看高清视频www| 国产野战对白在线观看| 老司机午夜十八禁免费视频| 久久人人97超碰香蕉20202| bbb黄色大片| 丝瓜视频免费看黄片| 国产精品秋霞免费鲁丝片| 国产高清视频在线播放一区 | 欧美日韩一级在线毛片| 91字幕亚洲| av片东京热男人的天堂| 国产在线观看jvid| 国产97色在线日韩免费| 少妇精品久久久久久久| 成年女人毛片免费观看观看9 | 久久人人爽人人片av| 欧美精品人与动牲交sv欧美| 欧美日韩av久久| 麻豆av在线久日| 亚洲欧美日韩另类电影网站| 十八禁人妻一区二区| 国产视频一区二区在线看| 国产精品偷伦视频观看了| 性色av一级| 日本av免费视频播放| 欧美成人精品欧美一级黄| 亚洲人成电影免费在线| 天天躁夜夜躁狠狠躁躁| 久久久精品国产亚洲av高清涩受| 国产亚洲av片在线观看秒播厂| 大片电影免费在线观看免费| 亚洲精品久久成人aⅴ小说| 国产欧美日韩综合在线一区二区| 国产视频首页在线观看| 美女午夜性视频免费| 欧美日韩亚洲国产一区二区在线观看 | xxx大片免费视频| 国产日韩欧美在线精品| 免费久久久久久久精品成人欧美视频| 久热这里只有精品99| 久久影院123| www.av在线官网国产| 99久久综合免费| 大话2 男鬼变身卡| 久久久久国产精品人妻一区二区| 青草久久国产| 亚洲专区中文字幕在线| 亚洲免费av在线视频| 国产免费又黄又爽又色| 国产成人啪精品午夜网站| 一本一本久久a久久精品综合妖精| av在线播放精品| 最近中文字幕2019免费版| av不卡在线播放| 人成视频在线观看免费观看| 黄网站色视频无遮挡免费观看| 观看av在线不卡| 久久精品国产a三级三级三级| 在线av久久热| 亚洲精品久久成人aⅴ小说| 欧美精品av麻豆av| 在线天堂中文资源库| 高清欧美精品videossex| 超色免费av| av天堂在线播放| 欧美少妇被猛烈插入视频| 一本一本久久a久久精品综合妖精|