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

    基于MaxEnt與GIS的黃精生態(tài)適宜分布區(qū)研究

    2025-06-23 00:00:00蔣絲絲葉菊何旭光汪成馬鈺潔唐媛媛馬敏
    草地學(xué)報(bào) 2025年5期
    關(guān)鍵詞:模型

    中圖分類號(hào):Q948.13 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1007-0435(2025)05-1622-09

    Abstract:To evaluate the distribution of the traditional Chinese medicine Polygonati Rhizoma in China and the main environmental factors affcting its distribution by means of ecological niche modellng,a total of 679 geographic locations of Polygonati Rhizoma distribution were collcted,and combined with 19 climatic and 4 topographic factors to simulate and analyze current and future habitat suitability by using MaxEnt algorithm and GIS technology,respectively. Currently,the total suitable area for Polygonati Rhizoma is 2142500km2 , mainly including Zhejiang,Jiangxi,eastern Yunnan,southern Sichuan,Chongqing and Hunan. Under future climate scenarios,the centre of mass migration is located in southern Shaanxi,northern Henan and Shanxi provinces,and the areas of both its high suitable and more suitable zones show an expanding trend,while the area of low suitable zones decreases and then increases.The key environmental factors affcting the distribution of Polygonati Rhizoma were elevation( 200-3000m ,mean temperature of driest quarter( -5°C~ 5°C ),precipitation of the driest month ( 19-215mm ) and precipitation of the warmest quarter (480-1500 mm).The MaxEnt model used in this study predicted accurate and reliable results ( .AUC=0.824 ,which provides a certain reference basis for the study of wild resources conservation and introduction of cultivation of Polygonati Rhizoma.

    Key words:Polygonati Rhizoma ;MaxEnt;GIS;Ecological suitability;Climatic factor

    中藥黃精為百合科黃精屬植物滇黃精(Polygonatum kingianum Coll. et Hemsl.)、黃精(P.sibiricum Red.)或多花黃精(P.cyrtonemaHua)的干燥根莖,味甘、性平,歸脾、肺、腎經(jīng),具有補(bǔ)氣養(yǎng)陰、潤肺、健脾、益腎等作用?,F(xiàn)代藥學(xué)研究證明,黃精富含多糖、生物堿、黃酮、三萜、甾體皂苷、木脂素及揮發(fā)油等活性物質(zhì),具有降血糖、提高免疫力、抗腫瘤、抗炎、抗衰老和抗氧化等藥理活性[1-2]。黃精是我國大宗藥材,也是民間常用食材,具有非常廣泛的應(yīng)用價(jià)值和開發(fā)前景。近年來,隨著黃精產(chǎn)業(yè)的不斷興起,野生資源已不能滿足開發(fā)利用現(xiàn)狀[3],四川、貴州、云南、安徽和湖北等省已開展了人工種植工程[4-5]。然而,黃精在大規(guī)模栽培過程中,存在著品種及其適生區(qū)篩選不當(dāng)?shù)膯栴},導(dǎo)致品種混亂、產(chǎn)量不高和品質(zhì)不佳等現(xiàn)象[6]。目前,研究物種生態(tài)適應(yīng)性的方法較多,如最大熵模型(Maximum Entropy,MaxEnt)、規(guī)則集遺傳算法(GeneticAlgorithm forRule-set Production,GARP)和生物氣候(BiocliMaticPrediction,BIO-CLIM)模型等[7-8]。其中,最大熵模型可通過較小的樣本量得到較高的預(yù)測精度,還能處理連續(xù)變量和分類變量,常用于動(dòng)植物適生區(qū)篩選及分布范圍變遷方面的研究,具有速度快、效率高、受樣本數(shù)量影響較小的優(yōu)點(diǎn),并且可以對環(huán)境變量進(jìn)行刀切法來得到各個(gè)環(huán)境變量的重要值,是目前表現(xiàn)最好、應(yīng)用最廣的物種分布模型[9-11]。通過MaxEnt模型與環(huán)境變量結(jié)合分析,發(fā)現(xiàn)意大利、法國和中國的秦嶺淮河一帶是黑腹尼虎天牛的最適分布區(qū)[12];常悅等[13]預(yù)測了加拿大一枝黃花在中國的適生區(qū),為其入侵動(dòng)態(tài)提供了早期預(yù)警機(jī)制;李政升等[14]利用MaxEnt模型預(yù)測了中亞紫菀木在當(dāng)前和未來的潛在適生區(qū),并得出影響其分布的主導(dǎo)氣候因子;宋晴等[15]也利用此模型研究興安盟草原灌叢植被潛在適生區(qū)分布情況,結(jié)果表明當(dāng)溫室氣體排放量持續(xù)增加時(shí),可能會(huì)引起興安盟地區(qū)水分匱缺,導(dǎo)致草原灌叢高適生區(qū)面積大幅降低,部分高適生區(qū)向中、低適生區(qū)轉(zhuǎn)化。將模型MaxEnt與地理信息系統(tǒng)(GIS)結(jié)合運(yùn)用,相比其他模型運(yùn)算速度快、穩(wěn)定性好、預(yù)測精度高[16]。因此,本研究利用MaxEnt模型、結(jié)合GIS技術(shù),研究黃精資源分布、適生區(qū)域和適宜生境條件的基礎(chǔ)上,探究當(dāng)前和未來氣候情境下,黃精適宜性分布位點(diǎn)與環(huán)境變量之間的關(guān)系,以及潛在適生區(qū)的演變情況,為其生態(tài)適宜性、道地性及規(guī)范化人工栽培技術(shù)提供科學(xué)依據(jù)與理論支撐。

    1 材料與方法

    1. 1 數(shù)據(jù)來源

    1.1.1樣點(diǎn)信息黃精的分布數(shù)據(jù)來源于網(wǎng)絡(luò)數(shù)據(jù)庫和文獻(xiàn)數(shù)據(jù)。通過文獻(xiàn)調(diào)查并收集《中國植物志》、中國數(shù)字植物標(biāo)本館(CVH,https://www.cvh.ac.cn/)和中國國家標(biāo)本平臺(tái)(NSII,https://www.nsii.org.cn)獲取黃精的地理分布數(shù)據(jù)。以“黃精\"為檢索詞在中國知網(wǎng)(CNKI,https://www.cnki.net/)檢索在2000一2023年公開發(fā)表的關(guān)于黃精的文獻(xiàn),記錄經(jīng)緯度。對有詳細(xì)位置但無經(jīng)緯度的標(biāo)本使用百度地圖軟件獲取其經(jīng)緯度信息。

    1.1.2 環(huán)境信息通過全球氣候數(shù)據(jù)庫(https://www.worldclim.org/)獲取環(huán)境數(shù)據(jù),通過世界土壤數(shù)據(jù)庫(http://www.fao.org/soils-portal/)獲取土壤數(shù)據(jù),并將上述數(shù)據(jù)信息(19項(xiàng)氣候因子和4個(gè)地形因子)轉(zhuǎn)換存儲(chǔ)為.asc文件格式。

    未來氣候數(shù)據(jù)模式的選擇來自“第六次國際耦合模式比較計(jì)劃(CMIP6)中國(北京)氣候中心中等分辨率氣候系統(tǒng)模式(BCC ~ CSM2\~MR)”,體現(xiàn)碳排放影響差異的SSP1-2.6(可持續(xù)發(fā)展路徑)、SSP2-4.5(中間發(fā)展路徑)SSP3-7.0(區(qū)域競爭發(fā)展路徑)和SSP5-8.5(化石燃料常規(guī)發(fā)展路徑)4種共享社會(huì)經(jīng)濟(jì)路徑[7]。

    1.1.3地理信息地圖數(shù)據(jù)從自然資源部標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站(http://bzdt.ch.mnr.gov.cn)下載的審圖號(hào)為GS(2024)0650號(hào)的標(biāo)準(zhǔn)地圖作為分析底圖。1.1.4數(shù)據(jù)分析軟件本研究采用生態(tài)分析軟件MaxEnt模型(3.4.4版本)、地圖處理GIS技術(shù)(10.7版本)和數(shù)據(jù)分析軟件SPSS(28版本)。

    1. 2 研究方法

    1.2.1MaxEnt模型參數(shù)設(shè)置使用MaxEnt3.4.4軟件導(dǎo)入黃精地理位置坐標(biāo),隨機(jī)以 25% 的分布點(diǎn)用于測試集(Testdata)進(jìn)行模型驗(yàn)證,剩余75% 的分布樣本作為訓(xùn)練集(Trainingdata)進(jìn)行模型構(gòu)建,設(shè)置最大迭代次數(shù)為 1×106 ,啟用響應(yīng)曲線,ROC曲線和刀切法(Jackknife),選擇刀切法來檢測變量,其他參數(shù)保持默認(rèn)設(shè)定[18]。利用響應(yīng)曲線來評價(jià)生態(tài)因子適宜度范圍,通過ROC曲線及其面積(AUC)評估模型準(zhǔn)確性,并運(yùn)用刀切法檢驗(yàn)生態(tài)因子的權(quán)重,以Loigstic輸出形式展現(xiàn)[19]。

    1.2.2MaxEnt模型精確度檢驗(yàn)受試者工作特征曲線(ROC曲線)中,將假陽性率作為橫坐標(biāo),將真陽性率作為縱坐標(biāo),依據(jù)受試者在特定刺激條件下,采用不同判斷標(biāo)準(zhǔn)得出的不同結(jié)果作出曲線[20]。應(yīng)用ROC曲線分析法對MaxEnt軟件預(yù)測的適宜區(qū)進(jìn)行精度檢驗(yàn),得到AUC值,按照上述方法,建構(gòu)基于最大熵理論的模型并經(jīng)過10次迭代實(shí)驗(yàn)。

    1.2.3適生區(qū)等級(jí)劃分進(jìn)行10次MaxEnt模型預(yù)測,從中選取AUC結(jié)果最大的一組,計(jì)算唯一值,將得到的asc.格式文件導(dǎo)入GIS軟件中,運(yùn)用“柵格分析工具-重分類\"工具,通過自然間斷點(diǎn)法劃分黃精生境適宜指數(shù),并加載到中國行政區(qū)域地圖中,將潛在生長區(qū)域分為四個(gè)不同的生境適應(yīng)性級(jí)別,得到預(yù)測生境適宜區(qū)域圖。

    2 結(jié)果與分析

    2.1 分布樣點(diǎn)信息

    為了減少空間上的過度擬合,將收集到的黃精分布點(diǎn)數(shù)據(jù),利用GIS剔除技術(shù),去除重復(fù)信息,以最小可查行政單位為準(zhǔn),篩選后共有679份黃精樣品分布信息可用于分析,并將其導(dǎo)入中國地圖,分布點(diǎn)見圖1。

    圖1黃精分布點(diǎn)的地理位置Fig.l Geographic location of the distribution points of the PolygonatiRhizoma

    海南省、臺(tái)灣省、廣東省南部、福建省南部及新疆自治區(qū)、西藏自治區(qū)和青海省交界的廣大地區(qū)無黃精的分布點(diǎn);黑龍江省、內(nèi)蒙古自治區(qū)、甘肅省、青海省東部、西藏自治區(qū)東南部和廣西省分布點(diǎn)較少;西南地區(qū)、湖南省、江西省、浙江省和安徽省分布點(diǎn)多且密集。由此可知,黃精在我國的分布極為廣泛,北至黑龍江、南至廣東、東至山東、西至新疆的31個(gè)省、市、自治區(qū)均有其分布。原因可能是中藥黃精的人藥品種較多,且各品種的生態(tài)適應(yīng)性均較強(qiáng)。其中,滇黃精主要分布在南部地區(qū),包括湖南、廣西、重慶、四川、貴州和云南,適宜海拔 700~ 3600m[21] 。多花黃精主要分布于四川、貴州、湖南、湖北、河南(南部和西部)江西、安徽、江蘇(南部)、浙江、福建、廣東(中部和北部)和廣西(北部),適宜海拔 500~2100m[22] 。黃精主要分布于北京、天津、河北、山西、內(nèi)蒙古、遼寧、黑龍江、江蘇、浙江、安徽、江西、山東、河南、湖南、湖北、廣西、重慶、四川、貴州、云南、陜西、甘肅、青海和寧夏,適宜海拔800~2800m[23] 。從各品種資源的分布情況來看,我國黃精資源可分為南北兩大類。其中,北黃精以大興安嶺南部、東北平原、內(nèi)蒙古高原和賀蘭山為分布中心,其原植物為黃精;南黃精以云貴高原和江南丘陵地帶為分布中心,其原植物為滇黃精和多花黃精。進(jìn)一步從分布范圍和資源量來看,南黃精的分布范圍更廣泛,具有更高的常見度和群集度。

    2.2主導(dǎo)生態(tài)因子的選擇

    MaxEnt模型通過物種已知地理分布點(diǎn)和相關(guān)環(huán)境變量預(yù)測物種的潛在分布,鑒于模型預(yù)測的變量之間存在較高的多重共線性,可能會(huì)導(dǎo)致構(gòu)建的模型出現(xiàn)數(shù)據(jù)的過擬合現(xiàn)象[24]。根據(jù)MaxEnt預(yù)實(shí)驗(yàn)得到貢獻(xiàn)率以及通過GIS技術(shù)的“數(shù)據(jù)管理工具-柵格-柵格處理-重采樣\"工具提取的數(shù)據(jù),進(jìn)行相關(guān)性分析。貢獻(xiàn)率是指每個(gè)主成分解釋的方差比例,貢獻(xiàn)率越高,說明該主成分在描述數(shù)據(jù)方面越重要。利用SPSS軟件計(jì)算19個(gè)氣候因子和4個(gè)地形因子的相關(guān)系數(shù),如果2個(gè)變量的相關(guān)系數(shù)大于0.8,則確認(rèn)變量之間存在相關(guān)性。進(jìn)一步結(jié)合變量對黃精分布的貢獻(xiàn)率,保留貢獻(xiàn)率較高的一個(gè)變量[25],結(jié)果如圖2所示。

    紅色代表正相關(guān),藍(lán)色代表負(fù)相關(guān),且相關(guān)性系數(shù)越大,顏色越深。根據(jù)相關(guān)性分析的結(jié)果進(jìn)行篩選,去除相關(guān)性大于0.8的生態(tài)因子,得到9個(gè)環(huán)境變量,包括海拔、坡度、等溫性、方面正弦、方面余弦、月平均晝夜溫差、最干季度平均氣溫、最干月降水量和最暖季度降水量。海拔、坡度、方面正弦和方面余弦是地理方向四個(gè)代表因子。其中,將方面正弦類比為“非水平度”的度量,即地形表面偏離水平面的程度,當(dāng)?shù)匦伪砻鏋樗矫鏁r(shí),方面余弦為1,方面正弦為0;當(dāng)?shù)匦伪砻鏋榇怪泵鏁r(shí),方面余弦為0,方面正弦為1。

    2.3模型預(yù)測精度

    為保證預(yù)測結(jié)果與實(shí)際數(shù)據(jù)的擬合度和精確性,進(jìn)行模型預(yù)測精度的檢驗(yàn),結(jié)果見圖3。AUC是ROC曲線以下的面積,用來衡量預(yù)測的精確性,其取值范圍為 ,數(shù)值越大表明預(yù)測的可靠性越高。AUC值的評價(jià)標(biāo)準(zhǔn): 0.5~0.6 表示預(yù)測失敗, 0.6~0.7 較差, 0.7~0.8 一般, 0.8~0.9 較好,0.9~1.0 非常好[26]。由圖3可見,黑色線代表AUC值為0.5,等于0.5時(shí),真實(shí)性最低,無應(yīng)用的價(jià)值;紅色線代表選擇訓(xùn)練數(shù)據(jù)AUC值為0.824,藍(lán)色線代表測試數(shù)據(jù)AUC值為0.846,表明MaxEnt模型對黃精在我國潛在分布區(qū)的預(yù)測效果較好。

    2.4生態(tài)因子適宜性分析

    運(yùn)用最大熵模型估算的過程中,通常用“交叉刀切法驗(yàn)證技術(shù)\"分析環(huán)境變量影響潛在分布的重要性,某一環(huán)境變量的貢獻(xiàn)率越高,其影響程度就越大。結(jié)果見表1。最暖季度降水量、最干月降水量、最干季度平均氣溫和海拔的貢獻(xiàn)率分別是47.7%,24.7%,14.5% 和 6.9% ,說明這4個(gè)因子是影響黃精分布的主導(dǎo)氣候因子;進(jìn)一步歸類分析得出,制約其潛在分布格局的首要因子是降水,其次是溫度,再次是海拔,這與文國衛(wèi)等2研究結(jié)果一致。一般情況下,降水量和氣溫隨海拔高度的不斷增加而遞減,即海拔通過影響降水量和氣溫而影響植物的分布,因此,影響黃精分布的關(guān)鍵因子可歸納為降水量和氣溫。

    圖3黃精潛在分布預(yù)測ROC曲線Fig.3PredictiveROC curves for the potential distribution of Polygonati Rhizoma

    貢獻(xiàn)率由高到低依次是:最暖季度降水量 gt; 最干月降水量 gt; 最干季度平均氣溫 gt; 海拔 gt; 坡度 gt; 月平均晝夜溫差gt;等溫性 gt; 方面余弦 gt; 方面正弦,且上述9個(gè)因子對黃精生境適宜性的累計(jì)貢獻(xiàn)率高達(dá) 100% ,因而選擇這些因子進(jìn)行刀切法分析,結(jié)果見圖4和圖5。

    圖4刀切法評價(jià)環(huán)境變量重要性的結(jié)果
    圖5黃精主導(dǎo)環(huán)境變量的響應(yīng)曲線Fig.5Response curves of environmental variables dominated by Polygonati Rhizoma

    圖4為訓(xùn)練增益的刀切圖評價(jià)結(jié)果,可見,最暖季度降水量(bio18)和最干月降水量(bio14)對模型準(zhǔn)確性的提升作用較為顯著,表明二者相較于其它變量,在獨(dú)立使用時(shí)能提供更多有效信息,對黃精適生環(huán)境的判定影響最大。另外,當(dāng)去掉方面正弦或方面余弦時(shí),對模型的預(yù)測能力影響不大,說明它們包含的有效信息相對較少。進(jìn)一步繪制模型輸出邏輯值與環(huán)境因子之間的響應(yīng)曲線28,根據(jù)貢獻(xiàn)率得到,最暖季度降水量、最干月降水量、最干季度平均氣溫和海拔4個(gè)環(huán)境因子的貢獻(xiàn)率較高,所以選擇這4個(gè)環(huán)境因子繪制曲線,結(jié)果均為單峰值,見圖5。

    根據(jù)黃精對環(huán)境變量的響應(yīng)情況,設(shè)置存在概率閾值為0.7,當(dāng)存在概率大于0.7時(shí),對應(yīng)的環(huán)境變量區(qū)間適合黃精的生長。可見,最暖季度降水量范圍為: 480~1500mm ,且降水量為 1000mm 時(shí),峰值最高。最干月降水量范圍為: 19~215mm ,且降水量在 60mm 時(shí),達(dá)到峰值后一直保持不變。最干季度平均氣溫范圍為: -5°C~5°C ,且平均氣溫為0°C 時(shí),峰值最高。海拔范圍為: 200~3000m ,且海拔為 1000m 時(shí),峰值最高。上述參數(shù)為黃精的規(guī)模化栽培區(qū)域篩選及栽培技術(shù)提供了可靠的數(shù)據(jù)支撐。

    2.5當(dāng)前情景下的地理分布

    本文采用可視性圖分析了中藥黃精的生境適宜區(qū),共分為4個(gè)區(qū)域,即非適生區(qū) (0~0.5) 、低適生區(qū) (0.5~0.7) 、中適生區(qū) (0.7~0.9) 、高適生區(qū)(大于0.9),結(jié)果見圖6和圖7。

    表1MaxEnt模型運(yùn)算的環(huán)境變量貢獻(xiàn)率( AUC=0.824 0Table1 Contribution of environmental variables to theMaxEntmodellingexercise(AUC =0,824 )
    圖6模型預(yù)測當(dāng)前氣候條件下黃精適宜分布區(qū)域 Fig. 6Model predictions of suitable distribution areas for PolygonatiRhizoma undercurrentclimatic conditions
    圖7當(dāng)前條件下黃精適生區(qū)面積占比Fig.7Percentage of the area suitable for Polygonati Rhizomaundercurrent conditions

    黃精適宜生境總面積為 ,約占我國國土面積的 22.3% 。其中,高適生區(qū)集中在浙江、湖南、江西等地區(qū),面積為 161500km2 ;中適生區(qū)集中在云南東部、四川南部、貴州東部等地區(qū),面積為 719300km2 ;低適生區(qū)集中在云南西部、貴州西部、四川東北部等地區(qū),面積為1261700km2 ;高、中和低適生區(qū)分別占總適生區(qū)的 7.5%.33.6% 和 58.9% 。非適生區(qū)集中在西北和東北地區(qū),如新疆、西藏、青海、內(nèi)蒙古、寧夏、黑龍江和吉林,這些地區(qū)雖然有黃精分布,但由于寒冷和干旱等氣候因素,生存環(huán)境并不理想[29],分布量極少且零散,可視為黃精生長發(fā)育的不適宜區(qū),面積為 7457500km2 。

    2.6 氣候變化情景下的地理分布

    未來氣候變化情景包括四種共享社會(huì)經(jīng)濟(jì)路徑情景:SSP1-2.6,SSP2-4.5,SSP3-7.0和SSP5-8.5[30] 。SSP1-2.6情景表示低脆弱性、低減緩壓力和低輻射強(qiáng)迫的綜合影響;SSP2-4.5情景表示中等社會(huì)脆弱性與中等輻射強(qiáng)迫的組合;SSP3-7.0情景表示高社會(huì)脆弱性與相對高的人為輻射強(qiáng)迫的組合;SSP5-8.5情景表示高強(qiáng)迫情景[31]。本文在上述情景模式下,預(yù)測分析了黃精未來16個(gè)時(shí)期的潛在分布,結(jié)果見表2和圖8。

    從表2可知,未來氣候變化情景下,高適生區(qū)面積呈現(xiàn)擴(kuò)大的趨勢,在SSP1-2.62021—2040情景下,高適生區(qū)面積最大,高達(dá) 377900km2 ,相較于現(xiàn)在,增加了 216400km2 。中適生區(qū)面積整體呈現(xiàn)擴(kuò)增趨勢,在SSP5-8.02081一2100情景下減少了17100km2 ,在 SSP1-2.62081-2100 情景下增大至 。低適生區(qū)面積出現(xiàn)先減少再增加的趨勢,在SSP1-2.62021—2040情景下,減少面積最多,減少值為 517700km2 ;在 SSP3–7.02081-2100 情景下,增加面積最多,增加值為 402900km2 。非適生區(qū)面積總體上呈現(xiàn)減少趨勢,在SSP3-7.02081一2100情景下,減少面積最多,減少值為479000km2 ;在SSP1-2.62021-2040情景下,非適生區(qū)面積高達(dá) 7586600km2 ??傮w來看,不同情景模式下,未來四個(gè)時(shí)期相較于當(dāng)前,黃精的總適生區(qū)面積呈現(xiàn)增加趨勢。

    從圖8可知,在SSP3-7.02081—2100情境下,低適生區(qū)面積變化最大,占比最高,占總面積的17.34% ;在SSP1-2.62021—2040情景下,中適生區(qū)面積、高適生區(qū)面積變化最大,占比最高,分別占總面積的 9.29% 和 3.94% 。

    表2不同時(shí)期黃精適生區(qū)面積Table2 The size of the suitablearea for PolygonatiRhizoma atdifferent times

    2.7 適宜生境質(zhì)心的轉(zhuǎn)移

    本文將黃精潛在適生分布區(qū)以幾何中心點(diǎn)來定義,以此方式來模擬在不同時(shí)期氣候情景下質(zhì)心遷移變化[32-33],結(jié)果見圖9和表3。

    在SSP1-2.6情景下,2041-2060時(shí)期相較于現(xiàn)代質(zhì)心位置向東遷移了 1205.09km;2021-2040 時(shí)期質(zhì)心位置向東北遷移了 1112.21km;2061-2080 時(shí)期和2081—2100時(shí)期質(zhì)心位置向西北遷移了79.17和 364.99km 。在SSP2-4.5情景下,2021—2040時(shí)期質(zhì)心位置向東北遷移 143.48km ;2041—2060時(shí)期和2061—2080時(shí)期質(zhì)心位置向西北遷移了349.24和 514.50km ;2081—2100時(shí)期質(zhì)心位置向東北遷移了 525.81km 。在SSP3-7.0情景下,

    2021—2040時(shí)期 時(shí)期質(zhì)心位置向西北遷移了 146.34km,391.57km;2061-2080E 寸期、2081—2100時(shí)期質(zhì)心位置向東北遷移了 689.93km 、845.17km 。在SSP5-8.5情景下, 2021-2040 時(shí)期質(zhì)心位置向西北遷移了 137.26km;2041-2060 時(shí)期 時(shí)期、 2081-2100 時(shí)期質(zhì)心位置向東北遷移了418.63,724.00和 1175.96km 。

    總體來看,從當(dāng)代到未來16個(gè)時(shí)期,均有不同的遷移位置。跟現(xiàn)代質(zhì)心位置相比,先是向東部遷移,再是向西北遷移,最后向東北遷移,表明未來的總體質(zhì)心遷移中心偏東北方向。其中,SSP1-2.6情境下,2041一2060時(shí)期質(zhì)心位置遷移的距離最多,2061—2080時(shí)期質(zhì)心位置遷移距離最小。

    表3不同時(shí)期氣候情景下黃精適生區(qū)質(zhì)心遷移距離Table3Distanceofcentre-of-massmigrationinasuitableareaforPolygonatiRhizomaunderdiferentclimaticscenarioindif ferent time periods
    圖8不同時(shí)期黃精適生區(qū)面積百分比Fig.8Percentage of suitablearea forPolygonati Rhizoma atdifferentperiods
    圖9不同時(shí)期氣候情景下黃精適生區(qū)質(zhì)心遷移位置Fig.9Position of the centre of mass migration in the Polygo nati Rhizoma habitat under different climate scenarios for different time periods

    3討論

    模型的精度取決于樣本的覆蓋度、區(qū)域和樣本量的多少,AUC值是衡量模型精度的指標(biāo)[34]。本研究運(yùn)用MaxEnt模型和GIS技術(shù)分析黃精地理分布現(xiàn)狀與生物環(huán)境變量的關(guān)系,篩選了主要環(huán)境變量,確定了未來氣候情境下對其地理變化的影響,充分保障預(yù)測的可靠性。收集的679個(gè)分布數(shù)據(jù)基本覆蓋黃精的分布范圍,得到AUC值為0.824,并將預(yù)測的高適生區(qū)與實(shí)際收集的分布位點(diǎn)進(jìn)行比對,發(fā)現(xiàn)分布范圍基本一致。利用MaxEnt模型預(yù)測黃精潛在分布的性能出色,預(yù)測結(jié)果精度高,有

    效避免了過擬合現(xiàn)象。

    現(xiàn)階段,黃精的高適生區(qū)主要包括湖南省、浙江省和江西省,而東北和西北地區(qū)不適宜黃精的生長。章鵬飛等[35采用“物種分布模型\"預(yù)測了多花黃精(P.cyrtonemaHua)的潛在分布區(qū),其高適生區(qū)也包括湖南和江西,而中、低適生區(qū)分布范圍存在較大差異,這是由于多花黃精是中藥黃精的基源植物之一,二者適生區(qū)范圍不能相互混淆。在雞頭黃精(P.sibiricumRed.)的潛在適生區(qū)研究中,朱寶潔等[36篩選的適生區(qū)與本研究篩選的地域存在較多的重疊性。在滇黃精(P.kingianumColl.etHemsl.)潛在分布和適宜性研究中,姚馨等[3得到其高度適生區(qū)主要位于云南高原中部、南部、東部和烏蒙山以東的貴州高原局部,中適生區(qū)主要分布于云南、四川和貴州省境內(nèi),少數(shù)位于廣西省境內(nèi),低適生區(qū)位于四川和貴州兩地。上述研究結(jié)果在中、高適生區(qū)分布及氣候類型上均與本研究結(jié)果基本一致。與現(xiàn)階段情況相比,未來黃精的質(zhì)心遷移中心位于陜西南部、河南北部和山西省。這可能是由于溫室效應(yīng)的不斷加劇,上述地區(qū)的未來降水量和溫度有所增加,更適于黃精的生長發(fā)育,從而成為未來環(huán)境模式下的高適生區(qū)。但是,該質(zhì)心遷移中心的變遷規(guī)律、影響因素及機(jī)制尚需進(jìn)一步深入研究。

    4結(jié)論

    本研究基于679個(gè)黃精分布點(diǎn)數(shù)據(jù)和篩選得到的9個(gè)環(huán)境因子數(shù)據(jù),利用優(yōu)化后的MaxEnt模型預(yù)測未來4種氣候情景不同時(shí)期黃精在國內(nèi)的潛在分布及其變化,以及影響其分布的主要環(huán)境因子,為黃精的生態(tài)適宜區(qū)域和分布提供了重要線索。降水量、氣溫和海拔是影響其分布的主要環(huán)境因子,當(dāng)前高適生區(qū)位于浙江、湖南、江西,未來氣候變化下質(zhì)心中心依次先向東、再向西北、后向東北遷移。

    參考文獻(xiàn)

    [1]袁怡菁,王秋紅.黃精化學(xué)成分、藥理作用研究進(jìn)展及質(zhì)量標(biāo)志物預(yù)測分析[J].中醫(yī)藥信息,2024,41(2):72-80,86

    [2]王慧,袁德培,曾楚華,等.黃精的藥理作用及臨床應(yīng)用研究進(jìn)展[J].湖北民族學(xué)院學(xué)報(bào)(醫(yī)學(xué)版),2017,34(2):58-60,64

    [3]蘇文田,劉躍鈞,蔣燕鋒,等.黃精產(chǎn)業(yè)發(fā)展現(xiàn)狀與可持續(xù)發(fā)展的建議[J].中國中藥雜志,2018,43(13):2831-2835

    [4]何俊玲,康琪,彭騰,等.專利視角下中藥黃精全產(chǎn)業(yè)鏈的發(fā)展對策研究[J].中藥與臨床,2023,14(6):78-82,96

    [5]張娜娜,潘慧云,侯芯芯,等.黃精研究現(xiàn)狀與可持續(xù)發(fā)展探析[J].安徽農(nóng)學(xué)通報(bào),2022,28(10):24-26

    [6]龍艷,麻生富,張國俊,等.黃精病蟲害防治研究進(jìn)展[J].耕作與栽培,2023,43(5):63-66

    [7]ZHONG Q,ZHANG JE,DITOMMASO A,et al. PredictinginvasionsofWedeliatrilobata(L.)Hitchc.withMaxentandGARPmodels[J].Journal of PlantResearch,2015,128(5):763-75

    [8]SEMWAL P D,PANDEY A,GORE P G,et al. Habitat prediction mapping using BioClim model for prioritizing germplasmcollection and conservation of an aquatic cash crop‘makhana’(Euryale ferox Salisb.)in India[J]. Genetic Resources andCrop Evolution,2021,68(8):3445-3456

    [9]GEBREWAHIDY,ABREHES,MERESAE,etal.Currentand future predicting potential areas of Oxytenanthera Abys-sinica(A. Richard) using MaxEnt model under climate changein Northern Ethiopia[J].Ecological Processes,202O,9(1):1-15

    [10]YANG J,HUANG Y,JIANG X,et al.Potential geographicaldistribution of the edangred.plant Isoetes under human activi-ties using MaxEnt and GARP[J]. Global Ecology and Conser-vation,2022,38:eo2186

    [11]郎顯鵬,樊如月,李青豐.中國北方地區(qū)蒙古韭潛在適生區(qū)分析[J].草地學(xué)報(bào),2023,31(11):3525-3534

    [12]王坤,石娟,梁特.基于MaxEnt和ArcGIS的黑腹尼虎天牛潛在地理分布預(yù)測[J].植物保護(hù)學(xué)報(bào),2023,50(6):1508-1517

    [13]常悅,李玉靜,李周旭,等.基于Maxent和ArcGIS加拿大一枝黃花適生區(qū)預(yù)測[J].農(nóng)學(xué)學(xué)報(bào),2024,14(3):40-47

    [14]李政升,馬玉壽,李有鑫,等,氣候變化下青藏高原中亞紫菀木潛在適生區(qū)的時(shí)空動(dòng)態(tài)變化[J].生態(tài)學(xué)雜志,2024,43(6):1566-1575

    [15]宋晴,付鴻莉,王鐵梅,等.興安盟草原灌叢植被潛在適生區(qū)分布模擬分析[J].草地學(xué)報(bào),2024,32(2):579-587

    [16]尹芳,朱家政,孟文睿,等.基于MaxEnt模型的陜西省菊芋種植潛在適宜區(qū)分析[J].中國農(nóng)業(yè)氣象,2024,45(3):268-280

    [17]TONGWEN W,YIXIONG L,YONGJIE F,et al. The beijingclimate center climate system model(BCC-CSM):the mainpIUgICSS HUI CIVIII J tU CIVII U LJ」. GCUSUICIuIIU IVIUUCIDevelopment,2019,12(4):1573-1600

    [18]唐楊欣,皮杰,劉新華,等.基于最大熵模型預(yù)測氣候變化下河蜆在中國的潛在分布[J].生態(tài)學(xué)報(bào),2023,43(10):4250-4259

    [19]王曉帆,段雨萱,金露露,等.基于優(yōu)化的最大熵模型預(yù)測中國高山櫟組植物的歷史、現(xiàn)狀與未來分布變化[J].生態(tài)學(xué)報(bào),2023,43(16):6590-6604

    [20]郭旭東,席少陽,楊莉萍,等.基于優(yōu)化MaxEnt模型的我國淫羊藿潛在適生區(qū)預(yù)測研究[J].中國中醫(yī)藥信息雜志,2024,31(3):1-7

    [21]田霞,吳楠,續(xù)若冰,等.滇黃精化學(xué)成分與提取、應(yīng)用研究進(jìn)展[J].中國食品工業(yè),2022(18):99-101

    [22]章登停,楊健,程銘恩,等.基于高光譜數(shù)據(jù)的多花黃精產(chǎn)地識(shí)別研究[J].中國中藥雜志,2023,48(16):4347-4361

    [23]羅霜,王劍波,陳柯宇,等.不同產(chǎn)地藥食同源黃精的成分及體外抗氧化研究[J].中國調(diào)味品,2023,48(6):59-65

    [24]AOUINTIH,MOUTAHIR H,TOUHAMI I,et al. Observedand predicted geographic distribution of Acer monspesulanumL.using the MaxEnt model in the context of climate change[J].Forests,2022,13(12):2049

    [25]王麗麗,楊采青,王瑛,等.全球入侵物種馬鈴薯塊莖蛾生態(tài)位轉(zhuǎn)移及適生區(qū)擴(kuò)展[J].應(yīng)用生態(tài)學(xué)報(bào),2024,35(3):797-805

    [26]HANLEYJA,MCNEIL B J. The meaning and use of the areaundera receiver operating characteristic(ROC)curve[J].Radiology,1982,143(1):29-36

    [27]文國衛(wèi),黃秋良,呂增偉,等.氣候變化情境下木荷潛在地理分布及生態(tài)適宜性分析[J].生態(tài)學(xué)報(bào),2023,43(16):6617-6626

    [28]董蘇君,馬松梅,張丹,等.新疆三種獨(dú)尾草的適宜分布及其對未來氣候變化的響應(yīng)[J].生態(tài)學(xué)報(bào),2022,42(21):8809-8817

    [29]胡穎,程詩明,柳新紅,等.藥用黃精遺傳資源評價(jià)、選育和栽培管理技術(shù)研究進(jìn)展[J].浙江林業(yè)科技,2023,43(6):117-121

    [30]謝孟,張學(xué)星,羅燕,等.基于MaxEnt模型的云南干熱河谷適生樹種選擇[J].生態(tài)學(xué)報(bào),2024,44(9):3689-3707

    [31]張彥靜,斯琴,胡潔,等.氣候變化情景下裸冠菊在中國的潛在適生區(qū)分布預(yù)測[J].生態(tài)學(xué)報(bào),2023,43(21):8852-8864

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

    [33]王彥珺,蘇鑫宇,朱成豪,等.基于MaxEnt模型和ArcGIS的金槐生態(tài)適宜性區(qū)劃研究[J].中國中醫(yī)藥信息雜志,2024,31(6):1-6

    [34]賀一鳴,王馳,王海濤,等.氣候變化對蒙古獲潛在適生區(qū)的影響[J].草地學(xué)報(bào),2023,31(2):540-550

    [35]章鵬飛,張虹,張小波,等.多花黃精生態(tài)適宜性區(qū)劃研究[J].中國中藥雜志,2020,45(13):3073-3078

    [36]朱寶潔,張瑾,龍飛,等.氣候背景下雞頭黃精適宜區(qū)預(yù)測[J].中藥與臨床,2023,14(3):6-9

    [37]姚馨,張金渝,萬清清,等.滇黃精的潛在分布與氣候適宜性分析[J].熱帶亞熱帶植物學(xué)報(bào),2018,26(5):439-448

    (責(zé)任編輯彭露茜)

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    欧美精品亚洲一区二区| 香蕉国产在线看| 久久久久久久午夜电影| 老司机靠b影院| 999久久久国产精品视频| 欧美日本亚洲视频在线播放| 美女大奶头视频| 亚洲片人在线观看| 亚洲熟妇中文字幕五十中出| 久久香蕉国产精品| 亚洲性夜色夜夜综合| 色播亚洲综合网| 国产熟女午夜一区二区三区| 国产亚洲精品第一综合不卡| 国产精品99久久99久久久不卡| 超碰成人久久| 亚洲精品在线美女| 757午夜福利合集在线观看| 久久久久国产一级毛片高清牌| 久久久久久久午夜电影| 色精品久久人妻99蜜桃| 很黄的视频免费| 国产97色在线日韩免费| 麻豆成人午夜福利视频| 亚洲精华国产精华精| 午夜精品在线福利| 特大巨黑吊av在线直播| 国产成人av激情在线播放| 久久九九热精品免费| 国产精品久久久av美女十八| 男人舔女人的私密视频| 嫩草影视91久久| 此物有八面人人有两片| 18禁观看日本| 亚洲一区二区三区色噜噜| 小说图片视频综合网站| 91九色精品人成在线观看| 亚洲精品av麻豆狂野| 国产精品野战在线观看| 国产av一区在线观看免费| 国产精品一区二区免费欧美| 亚洲国产中文字幕在线视频| 最近最新中文字幕大全电影3| 成人亚洲精品av一区二区| 欧美成人免费av一区二区三区| 日本撒尿小便嘘嘘汇集6| 又爽又黄无遮挡网站| 老司机福利观看| 国产精品亚洲av一区麻豆| 久久久久国产一级毛片高清牌| 成熟少妇高潮喷水视频| 精品欧美一区二区三区在线| 熟女少妇亚洲综合色aaa.| 欧美日韩亚洲国产一区二区在线观看| 久久久国产成人精品二区| 女人被狂操c到高潮| 少妇的丰满在线观看| 久久草成人影院| 国产精品久久电影中文字幕| 在线播放国产精品三级| 久久热在线av| 欧美最黄视频在线播放免费| 精品第一国产精品| 99热这里只有是精品50| 免费观看精品视频网站| 亚洲精品av麻豆狂野| 色老头精品视频在线观看| 757午夜福利合集在线观看| 激情在线观看视频在线高清| 久久久久国产精品人妻aⅴ院| 一二三四在线观看免费中文在| 无限看片的www在线观看| 成年版毛片免费区| 曰老女人黄片| 成人国产综合亚洲| 岛国在线观看网站| 国产亚洲av高清不卡| 国产精品影院久久| 美女午夜性视频免费| 久久久久久免费高清国产稀缺| 91字幕亚洲| 午夜福利高清视频| 午夜两性在线视频| 久久精品91蜜桃| 国产高清视频在线播放一区| 麻豆久久精品国产亚洲av| 身体一侧抽搐| 97碰自拍视频| 国产三级在线视频| 精品国产超薄肉色丝袜足j| 9191精品国产免费久久| 蜜桃久久精品国产亚洲av| 日本一二三区视频观看| 久久婷婷人人爽人人干人人爱| 日韩成人在线观看一区二区三区| 少妇被粗大的猛进出69影院| 国产黄a三级三级三级人| 久久香蕉国产精品| 亚洲精品国产精品久久久不卡| 精品免费久久久久久久清纯| 国产av又大| 中出人妻视频一区二区| 天堂av国产一区二区熟女人妻 | 欧美乱码精品一区二区三区| 精品少妇一区二区三区视频日本电影| 午夜福利高清视频| 国产精品久久久久久精品电影| 国产一区二区激情短视频| 亚洲国产精品成人综合色| 老司机午夜十八禁免费视频| xxx96com| 婷婷亚洲欧美| 人人妻,人人澡人人爽秒播| 国产一级毛片七仙女欲春2| √禁漫天堂资源中文www| 国产探花在线观看一区二区| 欧美性长视频在线观看| 亚洲乱码一区二区免费版| 99热只有精品国产| 亚洲av熟女| 日本一二三区视频观看| 久久久国产成人免费| 无遮挡黄片免费观看| 又紧又爽又黄一区二区| 亚洲一区中文字幕在线| 亚洲国产高清在线一区二区三| 亚洲精品中文字幕在线视频| 亚洲乱码一区二区免费版| 一进一出好大好爽视频| 欧美性猛交╳xxx乱大交人| 99riav亚洲国产免费| www国产在线视频色| 日韩大尺度精品在线看网址| 国产免费男女视频| 每晚都被弄得嗷嗷叫到高潮| 一级a爱片免费观看的视频| 99在线人妻在线中文字幕| 免费人成视频x8x8入口观看| 国产区一区二久久| 一级毛片高清免费大全| 国产高清视频在线播放一区| 日本 av在线| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址| 色噜噜av男人的天堂激情| 免费看十八禁软件| 亚洲国产日韩欧美精品在线观看 | 老司机午夜十八禁免费视频| 在线免费观看的www视频| 18禁黄网站禁片免费观看直播| 免费观看人在逋| 免费看日本二区| 久久精品亚洲精品国产色婷小说| 国产亚洲精品第一综合不卡| 黄频高清免费视频| 成人三级做爰电影| 成人高潮视频无遮挡免费网站| 99久久久亚洲精品蜜臀av| 中文字幕人成人乱码亚洲影| 日韩三级视频一区二区三区| 男女视频在线观看网站免费 | av福利片在线观看| 一区福利在线观看| 99精品欧美一区二区三区四区| 国产黄色小视频在线观看| 嫁个100分男人电影在线观看| 精品久久蜜臀av无| 亚洲精品av麻豆狂野| 全区人妻精品视频| 欧美在线黄色| 色综合站精品国产| 香蕉国产在线看| 日本在线视频免费播放| 欧美绝顶高潮抽搐喷水| ponron亚洲| 国产精品一区二区三区四区久久| 亚洲精品久久国产高清桃花| 一卡2卡三卡四卡精品乱码亚洲| 久久精品亚洲精品国产色婷小说| 一a级毛片在线观看| 午夜福利视频1000在线观看| 国产真实乱freesex| 精品日产1卡2卡| 99热这里只有是精品50| 99久久99久久久精品蜜桃| 此物有八面人人有两片| 国产精品爽爽va在线观看网站| 日韩有码中文字幕| 两个人看的免费小视频| 日韩国内少妇激情av| 成人国产一区最新在线观看| 国产黄a三级三级三级人| 亚洲欧美日韩高清专用| 成年免费大片在线观看| 老司机在亚洲福利影院| 在线永久观看黄色视频| 少妇裸体淫交视频免费看高清 | 不卡av一区二区三区| 久久久久亚洲av毛片大全| 久久精品综合一区二区三区| 久久午夜亚洲精品久久| 亚洲 欧美一区二区三区| 亚洲一区二区三区色噜噜| 午夜a级毛片| 欧美成人性av电影在线观看| av福利片在线| 超碰成人久久| 久久久国产欧美日韩av| 免费观看人在逋| 99久久国产精品久久久| 久久 成人 亚洲| 一区二区三区高清视频在线| 一进一出好大好爽视频| 两个人视频免费观看高清| 免费看日本二区| 一本综合久久免费| 999久久久国产精品视频| 18禁国产床啪视频网站| 国产成人aa在线观看| 成人一区二区视频在线观看| 精品国产超薄肉色丝袜足j| 99精品欧美一区二区三区四区| 最近视频中文字幕2019在线8| 日韩有码中文字幕| 中文字幕久久专区| 亚洲中文字幕日韩| 天堂动漫精品| 在线看三级毛片| 免费一级毛片在线播放高清视频| 一进一出抽搐gif免费好疼| 99精品欧美一区二区三区四区| 国产麻豆成人av免费视频| 香蕉久久夜色| 国产三级中文精品| 欧美性猛交黑人性爽| 欧美成狂野欧美在线观看| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品一区av在线观看| av视频在线观看入口| 久久久久国产一级毛片高清牌| 在线十欧美十亚洲十日本专区| av中文乱码字幕在线| 亚洲性夜色夜夜综合| 国产精品亚洲av一区麻豆| 欧美激情久久久久久爽电影| 国产精品久久久久久精品电影| 成人手机av| 窝窝影院91人妻| 12—13女人毛片做爰片一| 中出人妻视频一区二区| 亚洲午夜精品一区,二区,三区| 一级a爱片免费观看的视频| 久久久久久久精品吃奶| 午夜成年电影在线免费观看| 国产精品免费视频内射| 国产精品,欧美在线| 在线观看午夜福利视频| 淫妇啪啪啪对白视频| 一级黄色大片毛片| 91九色精品人成在线观看| 天天躁夜夜躁狠狠躁躁| 亚洲性夜色夜夜综合| 一本精品99久久精品77| www日本黄色视频网| 精品一区二区三区视频在线观看免费| 国产av在哪里看| a级毛片a级免费在线| 免费在线观看亚洲国产| 丁香六月欧美| 波多野结衣高清无吗| 精品国产乱码久久久久久男人| 亚洲第一欧美日韩一区二区三区| 又黄又爽又免费观看的视频| 婷婷六月久久综合丁香| 一本久久中文字幕| 国产激情欧美一区二区| 男人舔奶头视频| 久久人妻福利社区极品人妻图片| 国产单亲对白刺激| 亚洲欧美精品综合一区二区三区| 久久精品国产亚洲av香蕉五月| 丝袜人妻中文字幕| 久久国产精品人妻蜜桃| 欧美+亚洲+日韩+国产| 在线观看午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 在线观看美女被高潮喷水网站 | 超碰成人久久| 90打野战视频偷拍视频| 国产成人av激情在线播放| 日韩欧美精品v在线| 亚洲成av人片免费观看| 欧美日韩瑟瑟在线播放| 久久久久久国产a免费观看| 久久伊人香网站| 人人妻,人人澡人人爽秒播| 久久九九热精品免费| 亚洲国产精品合色在线| av有码第一页| 制服人妻中文乱码| 午夜免费观看网址| 日本在线视频免费播放| 精品久久久久久,| 久久这里只有精品19| 亚洲中文字幕一区二区三区有码在线看 | 午夜a级毛片| 欧美色欧美亚洲另类二区| 国内少妇人妻偷人精品xxx网站 | 国内精品久久久久精免费| 精品国产美女av久久久久小说| 成熟少妇高潮喷水视频| 欧美av亚洲av综合av国产av| 欧美色欧美亚洲另类二区| 日韩免费av在线播放| 在线观看免费日韩欧美大片| 国产精品一区二区精品视频观看| 两个人视频免费观看高清| 亚洲一区二区三区不卡视频| or卡值多少钱| 天堂av国产一区二区熟女人妻 | 亚洲欧洲精品一区二区精品久久久| 成人三级黄色视频| 亚洲第一欧美日韩一区二区三区| 美女免费视频网站| 精品午夜福利视频在线观看一区| 男女午夜视频在线观看| 欧美3d第一页| 无遮挡黄片免费观看| 欧美色欧美亚洲另类二区| 91字幕亚洲| 午夜精品久久久久久毛片777| 小说图片视频综合网站| 高清毛片免费观看视频网站| 女同久久另类99精品国产91| 俺也久久电影网| 变态另类成人亚洲欧美熟女| 夜夜爽天天搞| 免费无遮挡裸体视频| 黄频高清免费视频| 欧美日韩瑟瑟在线播放| 男女视频在线观看网站免费 | 桃红色精品国产亚洲av| 黑人巨大精品欧美一区二区mp4| 99久久无色码亚洲精品果冻| 亚洲欧美精品综合一区二区三区| 69av精品久久久久久| 日本熟妇午夜| 精品免费久久久久久久清纯| 婷婷亚洲欧美| 老鸭窝网址在线观看| 国产精品乱码一区二三区的特点| 日本a在线网址| 91九色精品人成在线观看| or卡值多少钱| 亚洲中文日韩欧美视频| 长腿黑丝高跟| 国产精品久久电影中文字幕| 久久国产精品影院| 两个人看的免费小视频| ponron亚洲| 一本久久中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲男人天堂网一区| 最新在线观看一区二区三区| 亚洲人成网站在线播放欧美日韩| 午夜福利免费观看在线| 99久久国产精品久久久| 中文字幕av在线有码专区| 男女午夜视频在线观看| 久久天躁狠狠躁夜夜2o2o| 中文字幕久久专区| 亚洲美女黄片视频| 男女视频在线观看网站免费 | 无遮挡黄片免费观看| 男人的好看免费观看在线视频 | 91成年电影在线观看| www.精华液| 91字幕亚洲| 18美女黄网站色大片免费观看| 黄色毛片三级朝国网站| 99re在线观看精品视频| 国产成人aa在线观看| av福利片在线观看| 香蕉丝袜av| 日韩大尺度精品在线看网址| 亚洲男人天堂网一区| 淫妇啪啪啪对白视频| 一二三四社区在线视频社区8| 看免费av毛片| a级毛片a级免费在线| 91在线观看av| 狠狠狠狠99中文字幕| 国产精品乱码一区二三区的特点| 日本 欧美在线| 黑人巨大精品欧美一区二区mp4| 一本久久中文字幕| 一进一出抽搐gif免费好疼| 女警被强在线播放| 亚洲乱码一区二区免费版| 久久久国产成人精品二区| 久久国产乱子伦精品免费另类| 精品久久蜜臀av无| 少妇被粗大的猛进出69影院| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区在线av高清观看| 精品高清国产在线一区| 国产黄色小视频在线观看| www日本在线高清视频| 亚洲av熟女| 黄色成人免费大全| 午夜福利视频1000在线观看| 亚洲专区中文字幕在线| 亚洲成a人片在线一区二区| 日本一二三区视频观看| 韩国av一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 一本综合久久免费| 一级片免费观看大全| 一边摸一边做爽爽视频免费| 久久久久国内视频| 日韩国内少妇激情av| 神马国产精品三级电影在线观看 | 九色成人免费人妻av| 色尼玛亚洲综合影院| 搡老岳熟女国产| 国产亚洲av嫩草精品影院| 午夜日韩欧美国产| 搡老岳熟女国产| 禁无遮挡网站| 色综合婷婷激情| 国产精品一区二区三区四区久久| 精品久久蜜臀av无| 亚洲精品在线美女| 亚洲精品美女久久久久99蜜臀| 久99久视频精品免费| 人人妻人人看人人澡| 五月伊人婷婷丁香| 色精品久久人妻99蜜桃| 成人欧美大片| 亚洲精品在线观看二区| 久久精品国产亚洲av高清一级| 最近最新中文字幕大全电影3| 人妻夜夜爽99麻豆av| 国产一区二区三区在线臀色熟女| www国产在线视频色| 一区二区三区高清视频在线| 18禁美女被吸乳视频| 怎么达到女性高潮| 91字幕亚洲| 黄色视频,在线免费观看| 欧美另类亚洲清纯唯美| 日韩成人在线观看一区二区三区| 国产成人影院久久av| 久久久久久亚洲精品国产蜜桃av| 国产三级中文精品| 亚洲成人久久爱视频| 一边摸一边做爽爽视频免费| 久久亚洲真实| 无限看片的www在线观看| 久久久国产精品麻豆| 日日爽夜夜爽网站| 国产欧美日韩精品亚洲av| 两个人看的免费小视频| 亚洲精品色激情综合| 久久热在线av| 日韩欧美三级三区| 国产又黄又爽又无遮挡在线| 国产欧美日韩一区二区三| 国产人伦9x9x在线观看| 国产精品一及| 久久热在线av| 国产精品一区二区三区四区久久| 黄色a级毛片大全视频| 日本精品一区二区三区蜜桃| www.精华液| 精华霜和精华液先用哪个| 97超级碰碰碰精品色视频在线观看| 国产av在哪里看| 亚洲人成伊人成综合网2020| 天天躁夜夜躁狠狠躁躁| 18禁裸乳无遮挡免费网站照片| 久久 成人 亚洲| 日韩欧美一区二区三区在线观看| 久久久久久亚洲精品国产蜜桃av| 97碰自拍视频| 久久久久国内视频| 国产精品免费视频内射| 亚洲aⅴ乱码一区二区在线播放 | 深夜精品福利| 国产亚洲欧美98| 在线a可以看的网站| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 女生性感内裤真人,穿戴方法视频| 亚洲第一欧美日韩一区二区三区| 亚洲中文字幕日韩| 一个人免费在线观看电影 | 日韩大尺度精品在线看网址| 五月伊人婷婷丁香| 亚洲国产中文字幕在线视频| 人人妻人人澡欧美一区二区| 女警被强在线播放| 午夜福利高清视频| 在线观看美女被高潮喷水网站 | 亚洲人成网站高清观看| 宅男免费午夜| 亚洲在线自拍视频| 在线免费观看的www视频| 国产av一区在线观看免费| 天天躁夜夜躁狠狠躁躁| 好男人电影高清在线观看| 一本一本综合久久| 色哟哟哟哟哟哟| 狂野欧美白嫩少妇大欣赏| 亚洲精品一区av在线观看| 欧美一级毛片孕妇| 波多野结衣高清无吗| 国产精华一区二区三区| 亚洲熟女毛片儿| 亚洲性夜色夜夜综合| 欧美日韩国产亚洲二区| 国产精品av久久久久免费| 特级一级黄色大片| 国产黄a三级三级三级人| 亚洲精品国产精品久久久不卡| 日本撒尿小便嘘嘘汇集6| 色哟哟哟哟哟哟| 99riav亚洲国产免费| 亚洲国产精品合色在线| 一个人免费在线观看的高清视频| 欧美乱码精品一区二区三区| 国产一区二区在线av高清观看| 一本精品99久久精品77| 两个人免费观看高清视频| www.精华液| 亚洲美女视频黄频| 国产单亲对白刺激| 老司机福利观看| 亚洲国产精品sss在线观看| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av在线| 欧美色欧美亚洲另类二区| 最好的美女福利视频网| 亚洲第一欧美日韩一区二区三区| 欧美av亚洲av综合av国产av| 久久草成人影院| 国产主播在线观看一区二区| 蜜桃久久精品国产亚洲av| 久久婷婷人人爽人人干人人爱| 欧美国产日韩亚洲一区| 香蕉av资源在线| 久久久久久久久免费视频了| 老汉色av国产亚洲站长工具| 99精品在免费线老司机午夜| 精品久久久久久久毛片微露脸| 一本大道久久a久久精品| 久久中文字幕一级| 日韩国内少妇激情av| 国产伦人伦偷精品视频| 国产伦在线观看视频一区| 国产亚洲精品久久久久5区| 曰老女人黄片| 嫁个100分男人电影在线观看| 人妻夜夜爽99麻豆av| 欧美色欧美亚洲另类二区| 美女免费视频网站| 中文在线观看免费www的网站 | 啦啦啦韩国在线观看视频| 两个人视频免费观看高清| 亚洲专区中文字幕在线| 伦理电影免费视频| 免费观看人在逋| 少妇粗大呻吟视频| 国产精品永久免费网站| 中文字幕熟女人妻在线| 国产精品亚洲一级av第二区| www国产在线视频色| 亚洲精品粉嫩美女一区| 一级黄色大片毛片| 国产区一区二久久| 国产高清视频在线观看网站| 亚洲av五月六月丁香网| 狂野欧美白嫩少妇大欣赏| 国产91精品成人一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看 | 一本大道久久a久久精品| 最近在线观看免费完整版| 婷婷六月久久综合丁香| 亚洲一区二区三区不卡视频| 午夜精品久久久久久毛片777| 国产av一区二区精品久久| 国产午夜精品论理片| 两性夫妻黄色片| 久久亚洲精品不卡| 亚洲av成人一区二区三| 一本一本综合久久| 两个人免费观看高清视频| 久久精品亚洲精品国产色婷小说| 麻豆成人午夜福利视频| 他把我摸到了高潮在线观看| 不卡av一区二区三区| 国产成人av教育| 岛国视频午夜一区免费看| 成年版毛片免费区| 久久久久久亚洲精品国产蜜桃av| 亚洲18禁久久av| 欧美不卡视频在线免费观看 | 国产高清激情床上av| 欧美日韩国产亚洲二区| 亚洲人成电影免费在线| 国产精品一及| 亚洲午夜精品一区,二区,三区| 99久久精品国产亚洲精品| 变态另类丝袜制服|