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

    基于Maxent和ARCGIS的醉魚草屬植物適生區(qū)預(yù)測

    2021-06-10 05:44:12歐陽晟秦思思陳王賦顏玉娟陳月華
    生態(tài)科學(xué) 2021年3期
    關(guān)鍵詞:適生區(qū)氣候因子季節(jié)

    歐陽晟, 秦思思, 陳王賦, 顏玉娟, 陳月華

    基于Maxent和ARCGIS的醉魚草屬植物適生區(qū)預(yù)測

    歐陽晟, 秦思思, 陳王賦, 顏玉娟*, 陳月華

    中南林業(yè)科技大學(xué)風(fēng)景園林學(xué)院, 長沙 410004

    為分析園林觀賞植物醉魚草屬在中國的潛在適生區(qū), 和環(huán)境因子對其在地理分布上的影響, 基于世界氣候數(shù)據(jù)網(wǎng)提供的19個氣候因子數(shù)據(jù), 以及中國地區(qū)內(nèi)128份可查證的分布點, 運用最大熵模型及相關(guān)生物地理軟件, 模擬醉魚草屬植物在中國的潛在分布區(qū), 并通過刀切法得出影響其分布的主導(dǎo)氣候因子。結(jié)果表明: 醉魚草屬植物潛在適生區(qū)為西南云貴地區(qū)和東南江浙一帶; 暖季降雨量、冷季均溫、溫度季節(jié)變異系數(shù)、降雨季節(jié)變異系數(shù)是影響醉魚草屬植物分布的主導(dǎo)因子, 其適宜生長在暖季降水量為200—2500 mm, 冷季均溫為-20—20 ℃的地區(qū); 醉魚草適生區(qū)與年降水量具有一定的相關(guān)性, 其分布邊界與400 mm年等降水量線大致重合。

    適生區(qū)預(yù)測; 最大熵模型; ARCGIS; 醉魚草屬

    0 前言

    醉魚草屬隸屬于馬錢科(參照恩格勒分類系統(tǒng)), 廣泛分布于江浙、兩湖兩廣及云貴等中國南方各省區(qū), 常生于山地、鄉(xiāng)村路緣、河畔林間, 多為灌木, 適應(yīng)性強。本研究涉及的種有醉魚草()、大葉醉魚草()、巴東醉魚草()、喉藥醉魚草()、短序醉魚草()、大花醉魚草()、云南醉魚草()、紫花醉魚草()、滇川醉魚草()等, 花大色艷, 芬芳四溢, 具有較高的觀賞價值, 在園林綠化中可用來點綴草地、護坡綠化、點綴山石、美化院落, 還可用來布置成花海景觀等, 是優(yōu)良的園林綠化植物, 醉魚草屬植物的花、葉及根都可作為藥材使用。國內(nèi)學(xué)者對醉魚草屬植物的研究主要集中在其化學(xué)成分測定及藥用價值方面, 也有對其園林應(yīng)用方面的研究, 對醉魚草屬植物適生區(qū)的預(yù)測在醫(yī)學(xué)和園林應(yīng)用領(lǐng)域都具有重要意義[1–5]。

    Maxent模型又叫最大熵模型, 是近年來新開發(fā)的一款軟件, 常用于預(yù)測生物的潛在適生區(qū)[6–9]。ArcGIS(地理信息系統(tǒng))作為近年來主流空間信息分析軟件, 具有強大的空間分析能力, 在眾多行業(yè)得到廣泛運用。根據(jù)Jaynes提出的最大熵理論, 在已知條件下, 事物在熵最大時最接近它的真實狀態(tài)[10–12]。

    1 材料與方法

    1.1 數(shù)據(jù)來源

    醉魚草屬植物的分布數(shù)據(jù)來源于http: //www. cvh.org.cn(中國數(shù)字標本館)、http://www.gbif.org/(全球生物多樣性信息平臺)、http://www.papc.cn/ (中國自然保護區(qū)資源平臺)、《中國植物志》及各地方植物志, 另外, 查詢醉魚草野外調(diào)查數(shù)據(jù)和文獻資料。初步獲得醉魚草屬植物標本信息270個, 直接使用其中具有準確經(jīng)緯度信息的標本, 對于已知詳細分布點的標本, 通過Google Earth等地圖軟件查詢其經(jīng)緯度信息, 最終整理出128份醉魚草屬植物地理分布信息, 主要發(fā)布在中國南方各省市, 具體分布如下(圖1)。

    海拔圖層獲取于DIVA-GIS軟件網(wǎng)站數(shù)據(jù)下載中心(http: //www.diva-gis.org), 植物覆蓋圖層數(shù)據(jù)來自于網(wǎng)站(http: //www.iscgm.org.), 公共氣候數(shù)據(jù)獲取于WorldClim網(wǎng)站(http: //www.Worldclim.org), 采用2.5 min數(shù)據(jù)庫, 共包括19個氣候因子[13–15]。分析用的底圖源于國家基礎(chǔ)地理信息系統(tǒng)官網(wǎng)(http: //nfgis.Nsdi.gov.cn /nfgis/chinese /c_xz.htm), 比例尺為1: 400萬。

    圖1 醉魚草屬植物標本分布圖

    Figure 1 Distribution map ofsamples in China

    1.2 相關(guān)軟件

    本研究所使用的軟件主要有Excel表格、DIVA-GIS 7.5、Maxent(最大熵模型)、ARCGIS 10.2(地理信息系統(tǒng))。Excel表格主要用于進行坐標數(shù)據(jù)的整理, DIVA-GIS用于繪制物種分布點地圖及氣候數(shù)據(jù)的格式轉(zhuǎn)換, Maxent用于分析預(yù)測醉魚草屬植物適生區(qū)間, ARCGIS用于繪制預(yù)測分布圖及空間分析。

    1.3 數(shù)據(jù)處理

    1.3.1 繪制物種分布點地圖

    利用DIVA-GIS 7.5軟件處理所獲得的醉魚草屬植物經(jīng)緯度數(shù)據(jù), 繪制而得到醉魚草屬植物在中國境內(nèi)的地理分布圖。再以海拔和植物覆蓋圖層為底圖, 在ARCGIS 10.2軟件中繪制出醉魚草屬植物的海拔分布圖及分布地的植物覆蓋圖[16–17]。

    1.3.2 分析預(yù)測醉魚草屬植物適生區(qū)間

    根據(jù)Maxent模型需要, 將所得醉魚草屬植物的經(jīng)緯度坐標數(shù)據(jù)通過Excel表格轉(zhuǎn)換“植物名+經(jīng)度+緯度”的csv格式, 導(dǎo)入Maxent模型中的Sample菜單下; 將WorldClim網(wǎng)站下載的19個氣候因子數(shù)據(jù), 包括年平均溫(bio1)、晝夜溫月均值(bio2)、晝夜溫差與年溫差比值(bio3)、溫度變化方差(bio4)、最熱月份最高溫(bio5)、最冷月份最低溫(bio6)、年溫變化范圍(bio7)、最濕季節(jié)平均溫度(bio8)、最干季節(jié)平均溫度(bio9)、最暖季節(jié)平均溫度(bio10)、最冷季節(jié)平均溫度(bio11)、年平均濕度(bio12)、最濕月濕度(bio13)、最干月濕度(bio14)、濕度變化方差(bio15)、最濕季節(jié)濕度(bio16)、最干季節(jié)濕度(bio17)、最暖季節(jié)平均濕度(bio18)、最冷季節(jié)平均濕度(bio19), 通過DIVA-GIS 7.5軟件的Data菜單下Export Grdfile-Multiple files, 將19個氣候因子的grd文件全部轉(zhuǎn)換為ASCII 格式, 導(dǎo)入到Environmental layers菜單下。在隨機狀態(tài)下, 選擇75%的數(shù)據(jù)作為訓(xùn)練集組建模型, 剩下25%的數(shù)據(jù)作為測試集用于驗證, 勾選Create response curves、Make pictures of prediction、Do Jackknife to measure variable importance復(fù)選框, 在setting菜單下characteristic curve設(shè)為10次重復(fù), 生成10個預(yù)測隨機模型, 其他參數(shù)均為軟件默認值。

    1.3.3 預(yù)測分布圖及空間分析

    模擬結(jié)果以ASCII格式輸出, 再轉(zhuǎn)入ARCGIS 10.2軟件中加載Maxent運行結(jié)果, 利用arctoobox工具箱, 將ASCII文件轉(zhuǎn)為光柵文件類型, 選擇輸出文件類型為float, 繪制出預(yù)測結(jié)果的分布圖[18–19]。

    2 結(jié)果與分析

    2.1 模型精度分析

    將實際沒有該物種分布, 而被預(yù)測為陽性的概率(假陽性率)為橫坐標; 實際有該物種分布, 且預(yù)測為陽性的概率(真陽性率)為縱坐標, 繪制成受試者工作特征曲線(Reciever Operating Characteristic Curve, ROC曲線), 選取受試者工作特征曲線對模型的預(yù)測精度進行評估, 其具體指標是采用ROC曲線與橫坐標所圍成的面積(AUC值)來判斷模型預(yù)測的優(yōu)良程度(表1), AUC值越接近1, 表明與隨機分布相距越遠, 環(huán)境變量與預(yù)測結(jié)果越準確[20]。

    利用Maxent模型預(yù)測醉魚草的AUC值為0.990(圖2), 表明最大熵模型對醉魚草屬植物在中國適生區(qū)的預(yù)測結(jié)果極佳。

    表1 AUC值與預(yù)測結(jié)果間的關(guān)系

    圖2 利用Maxent模型預(yù)測醉魚草的ROC曲線

    Figure 2 Maxent Model predicting the ROC curve of the

    2.2 主導(dǎo)因子分析

    刀切法常用于分析各個環(huán)境變量對預(yù)測結(jié)果的影響程度, 運用刀切法分析影響醉魚草屬植物分布的主要環(huán)境變量時, 其原理是在運算過程中, 對環(huán)境變量進行依次省略, 再用剩余的變量建立模型, 對省略變量與遺漏誤差之間的相關(guān)性進行分析, 如果遺漏誤差因某個環(huán)境變量的缺失而顯著升高, 則表明該環(huán)境參數(shù)對模型的預(yù)測結(jié)果影響突出[21]。

    將氣候因子導(dǎo)入Maxent模型后, 各氣候因子對醉魚草屬植物分布的相對貢獻率見表2。

    從表2可知, 相對貢獻率大于10%的氣候因子有暖季降雨量(Bio18)、溫度季節(jié)變異系數(shù)(Bio4)、冷季均溫(Bio11)、降雨季節(jié)變異系數(shù)(Bio15), 屬主導(dǎo)因子, 對醉魚草屬植物分布概率的影響見圖3。其中Bio18對醉魚草屬植物適生區(qū)影響最大, 貢獻率為53.2%。其次, Bio4和Bio11以15.7%和12.4%的貢獻率居第二位與第三位, 說明溫度變化條件也對醉魚草屬植物的分布產(chǎn)生了一定的影響。

    從圖3(a)可以看出: 當暖季降雨量大于200 mm時, Buddleja分布的概率隨降雨量的增加而顯著提升; 當降雨量600 mm時, 其分布概率仍隨著降雨量的增加上升速度減慢; 當降雨量達到2500 mm時, 其分布概率達到平衡, 所以, 暖季降雨量大于2500 mm時將不利于醉魚草屬植物的生長。圖3(b)顯示, 溫度季節(jié)變異系數(shù)大于0后, 醉魚草屬植物分布概率開始迅速上升, 在4000—13000范圍內(nèi)變化緩慢, 大于15000后趨近平衡。圖3(c)顯示: 冷季均溫在-20—20 ℃范圍內(nèi)適合醉魚草屬植物的生長, 在均溫為0℃時達到最大分布概率, 低于-25 ℃和高于25 ℃均不利于其生長。圖3(d)顯示, 降雨季節(jié)變化系數(shù)(Bio15)在85時最較適合醉魚草屬植物的生長。

    2.3 地理分布及預(yù)測結(jié)果分析

    2.3.1 地理分布

    從海拔分布圖(圖4)和植被覆蓋圖(圖5)可以看出南方地區(qū)的各種地形地段, 在熱帶氣候地區(qū)少有分布, 高海拔寒冷地區(qū)難以生存, 標本數(shù)據(jù)主要分布于云貴高原、四川盆地、東南丘陵及長江中下游地區(qū), 海拔跨度從200 m到2700 m, 其中250—500 m(東南丘陵)和1000—2000 m(云貴高原)分布相對集中。生長地植被類型主要為亞熱帶常綠闊葉林及常綠針葉林, 多于殼斗科、樟科、木蘭科、松柏科、金縷梅等喬木植物混生, 常見的混生喬木有栲樹、香樟、木荷、馬尾松、杉木等; 常與草本和灌木伴生, 其中以禾本科及菊科植物居多。

    表2 Maxent模型中19氣候因子相對貢獻率

    圖3 主導(dǎo)因子對醉魚草屬植物分布概率的影響曲線. 暖季降雨量(a) 、溫度季節(jié)變異系數(shù)(b)、冷季均溫(c)、降雨季節(jié)變異系數(shù)(d)

    Figure 3 The curve of probability of the dominant factor released. Response curve of precipitation of warmest quarter(a), response curve of temperature seasonality(b), response curve of mean temperature of coldest quarter(c), response curve of precipitation seasonality(d)

    圖4 醉魚草屬植物分布點的海拔圖

    Figure 4 Altitudes in the geographical distribution map of

    圖5 醉魚草屬分布點的植物覆蓋圖

    Figure 5 landcover type in the geographical of distribution map of

    醉魚草屬植物集中分布的地區(qū)多處于亞熱帶季風(fēng)氣候區(qū), 臨近太平洋, 夏季容易受到海上濕熱氣流的影響而產(chǎn)生充沛的降水, 滿足常綠闊葉植物的水分和溫度的需求。在生長條件滿足的情況下, 醉魚草屬植物可生長于多種植被環(huán)境中: 既能于灌木草本間雜混生, 也能處于高大喬木的遮蔽之下, 或水邊路旁。與菊科植物、蕨類植物、禾本科植物等都能形成良好的植物群落。

    2.3.2 預(yù)測結(jié)果分析

    將Maxent預(yù)測結(jié)果導(dǎo)入ARCGIS后得到醉魚草屬植物適生區(qū)間的預(yù)測結(jié)果(圖6), 形成以云貴高原及四川盆地為中心和以長江中下游地區(qū)為中心的“一大一小”的分布格局, 其中主要分布的省份為云南、貴州、四川、重慶(一大)和安徽、浙江、福建(一小); 此外, 甘肅、湖南、湖北、江西、江蘇、廣西等地區(qū)也較為適合其生長繁殖。貴州省為預(yù)測分布面積最大的省區(qū), 大部分北方地區(qū)及熱帶地區(qū)不適合醉魚草屬植物的生長。

    預(yù)測結(jié)果表明, 醉魚草屬植物在我國的適生范圍主要集中在西南、華東地區(qū)。其中高度適生區(qū)包括貴州絕大部分, 四川、云南、重慶大部分以及廣西北部, 福建西北及東南部, 江蘇南部。適生區(qū)包括廣西、廣東、湖南、湖北、江西、浙江、安徽、甘肅大部分地區(qū), 西藏和寧夏也有零星分布。與當前分布點不同的是, 甘肅、寧夏、湖北、安徽將成為醉魚草屬植物新的適生區(qū)間。新的適生區(qū)間集中在30°N附近, 較當前分布地區(qū)有所北移, 區(qū)間跨度更為廣闊。

    預(yù)測中高度適生區(qū)間及適生區(qū)間主要集中在20—40°N, 使得醉魚草可作為這些地區(qū)的鄉(xiāng)土草本而進行開發(fā), 通過花壇、花境等應(yīng)用形式進行園林栽培, 而醉魚草屬植物在這些地區(qū)極強的適應(yīng)性也可減少其養(yǎng)護管理費用, 降低成本。

    2.3.3 適生區(qū)主導(dǎo)因子與地理特征預(yù)測

    利用刀切法對19個氣候因子分析后, 得到4個關(guān)鍵影響因子(主導(dǎo)因子): 暖季降雨量、溫度季節(jié)變異系數(shù)、冷季均溫、降雨季節(jié)變異系數(shù)。當暖季降雨量小于200 mm時, 醉魚草植物的分布概率極低, 大于200 mm后, 開始增長(圖3(a))。結(jié)合我國年等降水量線和醉魚草植物標本分布圖(圖7, 年等降水量線圖來自中國國家地理網(wǎng)), 發(fā)現(xiàn)現(xiàn)有醉魚草屬植物在中國分布地區(qū)的邊界與400 mm年等降水量大致重合; 醉魚草分布較密集的地區(qū)與800 mm年等降水量線相吻合。從等降水量線的意義看, 它既是我國地理上的分界線, 也是氣候區(qū)的分界線; 400 mm年等降水量線是我國半濕潤與半干旱地區(qū)的分界線, 也是東部季風(fēng)區(qū)與西北干旱半干旱區(qū)的分界線; 800 mm等降水量線作為傳統(tǒng)意義上的南北分界線, 是溫帶季風(fēng)氣候與亞熱帶季風(fēng)氣候的分界線。由圖3(c)得知,當冷季均溫在0℃的時候, 分布概率達到最大, 這與秦嶺—淮河的中國南北分界線大致重合。根據(jù)以上氣候特征可知: 降水量和溫度對醉魚草的生長影響較其他氣候因子大, 其難以在干旱沙漠地區(qū)大面積生長繁殖, 比較適合生長于南方濕潤溫暖的氣候。在培育過程中要注意水分和溫度的控制, 夏季給予其充足的水分, 北方地區(qū)需要對其加以防護措施, 以有利其度過冬季的嚴寒, 得到更好的栽培效果[22-27]。

    圖6 利用Maxent模型預(yù)測的醉魚草屬植物在中國的潛在適生區(qū)

    Figure 6 A prediction of the potential distribution range ofusing Maxent model

    圖7 醉魚草屬植物標本分布與中國年等降水量線

    Figure 7specimen and annual precipitation line of china

    3 討論

    根據(jù)醉魚草屬植物標本分布圖(圖1)及中國植物志等網(wǎng)站、文獻及相關(guān)信息的查詢, 醉魚草現(xiàn)階段主要分布于中國南方各省, 為典型的亞熱帶季風(fēng)氣候, 此氣候類型容易在南北緯25—35 ℃亞熱帶大陸東部, 夏季炎熱多雨, 冬季溫暖少雨, 冬季的平均溫度可在0 ℃以上, 這與預(yù)測結(jié)果中的主導(dǎo)因子相吻合: 夏季的高溫多雨使得暖季降雨量充足; 四季明顯的氣候使得氣溫的季節(jié)變化顯著; 冬季寒冷, 均溫能維持在0 ℃以上; 夏季多雨, 冬季少雨, 降雨量的季節(jié)變化十分明顯。這些都表明, 醉魚草屬植物的生長和環(huán)境氣候條件息息相關(guān), 其更適合生存于亞熱帶季風(fēng)氣候區(qū)。在預(yù)測結(jié)果中, 醉魚草植物的分布范圍有往北移的趨勢, 甘肅、江蘇、河南、河北等一些非熱帶季風(fēng)氣候區(qū)也成為醉魚草生存概率極高的區(qū)域, 這與全球氣候變暖相吻合。

    當前醉魚草屬植物主要運用于南方各城市中, 因氣候條件限制, 很少被配置于北方的城市園林。預(yù)測結(jié)果顯示, 醉魚草屬植物的適生區(qū)間將逐步向北擴展, 中原一帶將成為其新的適生地段。醉魚草生長較快, 易為栽植, 成活率高, 生命力強, 植株不高, 是良好的花壇植物, 還可作切花使用, 具有很好的發(fā)展應(yīng)用前景, 可先作為南方的鄉(xiāng)土物種進行開發(fā)栽培, 在南方地區(qū)進行推廣培育, 等到氣候條件適宜, 再引種北方進行繁育。

    [1] 中國科學(xué)院《中國植物志》編輯委員會.中國植物志[M].北京: 科學(xué)出版社, 2004.

    [2] 顏玉娟, 趙虎. 園林地被植物的選擇與應(yīng)用[J]. 湖北林業(yè)科技, 2005, (133): 54–57.

    [3] 謝紅梅. 湖南醉魚草屬植物資源評價及園林應(yīng)用研究[D].長沙: 中南林業(yè)科技大學(xué), 2017.

    [4] Zhang Wei, Li Zheng, Xu Fengqing, et al. Mimengosides J and K: two new neuroprotective triterpenoids from the fruits of[J]. Journal of Asian Natural Products Research, 2018(11): 1–9.

    [5] 張玉梅. 醉魚草扦插繁殖技術(shù)及其生理生化研究[D]. 長沙: 中南林業(yè)科技大學(xué), 2018.

    [6] PHILLIPS S J, ANDERSON R P, SCHAPRISE R E. Maximum entropy modeling of species geographic distributions[J]. Ecological Modelling, 2006, 190(3/4): 231–259.

    [7] KUMAR S, STOHLGREN T J. Maxent modeling for predicting suitable habitat for threatened and endangered treein New Caledonia[J]. Journal of Ecology and the Natural Environment, 2009, 1(4): 94–98.

    [8] LOZIE R J D, ANIELLO P, HICKERSON M J. Predicting the distribution of Sasquatch in western north America: Anything goes with ecological niche modeling[J]. Journal of Biogeography, 2009, 36(9): 1623–1627.

    [9] 曹向鋒, 錢國良, 胡白石, 等. 采用生態(tài)位模型預(yù)測黃頂菊在中國的潛在適生區(qū)[J]. 應(yīng)用生態(tài)學(xué)報, 2010, 21(12): 3063–3069.

    [10] 龔維, 夏青, 陳紅鋒, 等. 珍稀樹種伯樂樹的潛在適生區(qū)預(yù)測[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報, 2015, 36(4): 98–104.

    [11] 邵揚, 葉丹, 歐陽臻, 等. 薄荷的生境適宜性區(qū)劃及品質(zhì)區(qū)劃研究[J].中國中藥雜志, 2016, 41(17): 3169–3175.

    [12] 謝春平. 南方紅豆杉分布生態(tài)適應(yīng)性分析[J]. 熱帶地理, 2014, 34(3): 359–365.

    [13] 田聰, 穆尼熱, 朱忠艷, 等. 基于DIVA-GIS快速獲取氣候信息的方法[J]. 農(nóng)業(yè)學(xué)報, 2015, 5(5): 109–113.

    [14] 杜中修, 鄭勇奇, 魯法典, 等. 火炬樹在我國的適生區(qū)預(yù)測[J]. 廣東農(nóng)業(yè)科學(xué), 2011, (2): 50–52.

    [15] 孔令斌, 林偉, 李志紅, 等. 基于CLIMEX和DIVA-GIS的瓜實蠅潛在地理分布預(yù)測[J]. 植物保護學(xué)報, 2008, 35(2): 148–154.

    [16] 謝春平. 基于DIVA-GIS生物地理分布圖的繪制[J]. 湖北農(nóng)業(yè)科學(xué), 2011, 50(11): 2345–2348.

    [17] 毛俐慧, 李垚, 劉暢, 等. 基于Maxent模型預(yù)測細葉小羽蘚在中國的潛在分布區(qū)[J]. 生態(tài)學(xué)雜志, 2017, 36(1): 54–60.

    [18] 熊中人, 張曉晨, 鄒旭, 等. 中國天山花楸適生區(qū)預(yù)測及其對氣候變化的響應(yīng)[J]. 生態(tài)科學(xué), 2019, 38(5): 44–51.

    [19] 閆宇航, 岑云峰, 張鵬巖, 等. 基于MaxEnt模型的中國馬尾松分布格局及未來變化[J]. 生態(tài)學(xué)雜志, 2019, 38(9): 2896–2901.

    [20] 柳曉燕, 李俊生, 趙彩云, 等.基于MAXENT, 模型和ArcGIS預(yù)測豚草在中國的適生區(qū)[J].植物保護學(xué)報, 2016, 43(6): 1041–1048

    [21] 劉建宏, 熊小真, 潘涌智, 等.滇柏大痣小蜂Roques & Pan的潛在分布區(qū)預(yù)測[J].中國農(nóng)學(xué)通報, 2011, 27(10): 39–43

    [22] MENON S, CHOUDHURY B I, KHAN M L, et al. Ecological niche modeling and local knowledge predict new populations of, a critically endan-gered tree species[J]. Endangered Species Research, 2010, 11(2): 175–181.

    [23] BOURG N A, MCSHEA W J, GILL D E. Putting a CART before the search: Successful habitat prediction for a rare forest herb[J]. Ecology, 2005, 86(10): 2793–2804.

    [24] VANAGAS G. Receiver operating characteristic curves and comparison of cardiac surgery risk stratification systems[J]. Interactive Cardiovascular and Thoracic Surgery, 2004, 3(2): 319–322.

    [25] ROBERTO M, ZAMORA R. Predictive modeling of microhabitats for endemic birds in south Chilean temperate forestsusing Maximum entropy[J]. Ecological Informatics, 2011, 6: 364–370.

    [26] 秦譽嘉, 藍帥, 盧國彩, 等. 氣候變化條件下櫻桃繞實蠅在中國的潛在地理分布預(yù)測[J]. 植物保護學(xué)報2019, 46(1): 63–70.

    [27] 張興旺, 李垚, 謝艷萍, 等. 氣候變化對黃山花楸潛在地理分布的影響[J]. 植物資源與環(huán)境學(xué)報, 2018, 27(4): 31–34.

    Prediction of potential distribution ofbased on the Maxent and ARCGIS

    OUYANG Sheng, QIN Sisi, CHEN Wangfu, YAN Yujuan*, CHEN Yuehua**

    College of Landscape Architecture, Central South University of Forestry and Technology, Changsha 410004, China

    For the analysis of the potential distribution area of the ornamental plantin China, and the impact of environmental factors on the geographical distribution, based on 19climatic factors data provided by the World Climate Data Network and 128 recorded distribution points in China, we used the maximum entropy (Maxent) model and related biogeographic software to simulate the potential distribution area ofin China, and used the Jackknife method to obtain the dominant climatic factors affecting their distribution. The results showed thatthe potential suitable areas ofwere Yunnan Guizhou region in Southwest China area and Jiangsu and Zhejiang in Southeast China. Warm season rainfall, uniform temperature in cold season, seasonal variation coefficient of temperature, and seasonal variation coefficient of rainfall were affecting the dominant factor for the distribution of. Warm season rainfall 200-2500 mm and uniform temperature in cold season -20 to 20℃ were a suitable growing area. There was a certain correlation between the suitable area of grass and annual precipitation. The distribution boundary roughly coincided with the 400 mm annual isoprecipitation line.

    potential distribution; maximum entropy model; DIVA-GIS;

    歐陽晟, 秦思思, 陳王賦, 顏玉娟, 陳月華. 基于Maxent和ARCGIS的醉魚草屬植物適生區(qū)預(yù)測[J]. 生態(tài)科學(xué), 2021, 40(3): 165–173.

    OUYANG Sheng, Qin Sisi, CHEN Wangfu, YAN Yujuan, CHEN Yuehua. Prediction of potential distribution ofbased on the Maxent and ARCGIS[J]. Ecological Science, 2021, 40(3): 165–173.

    10.14108/j.cnki.1008-8873.2021.03.020

    S157.2

    A

    1008-8873(2021)03-165-09

    2019-11-13;

    2020-01-23

    國家林業(yè)局“十三五”重點學(xué)科項目([2015]44); 湖南省教育廳“十二五”重點學(xué)科項目(2011-76); 湖南高校一流本科專業(yè)(湘教通〔2019〕138號)(園林專業(yè))

    歐陽晟(1995—), 男, 湖南衡陽人, 碩士研究生, 研究方向為風(fēng)景園林規(guī)劃與設(shè)計, E-mail: 959041792@qq.com

    顏玉娟, 女, 四川達州人, 副教授、博士, 主要研究領(lǐng)域為植物景觀規(guī)劃與設(shè)計、園林植物資源與應(yīng)用、城市綠地系統(tǒng)規(guī)劃、植物景觀維護、森林公園植物景觀規(guī)劃, E-mail:1964610285@qq.com

    猜你喜歡
    適生區(qū)氣候因子季節(jié)
    青海省果洛州高原鼠兔宜生區(qū)分布及其對氣候變化的響應(yīng)
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預(yù)測
    未來氣候條件下當歸適生區(qū)預(yù)測及時空變化分析
    氣候變化下瀕危植物半日花在中國的潛在分布
    我喜歡的季節(jié)7
    季節(jié)蠕變
    英語文摘(2019年5期)2019-07-13 05:50:06
    季節(jié)的變換
    幼兒畫刊(2018年10期)2018-10-27 05:44:36
    巴拉圭瓜多竹適生區(qū)分布研究
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    花的季節(jié)
    日韩电影二区| 9热在线视频观看99| 久久久久久久久免费视频了| 在线 av 中文字幕| 久久国产精品男人的天堂亚洲| av有码第一页| 久久这里只有精品19| 一级片'在线观看视频| 少妇猛男粗大的猛烈进出视频| 亚洲成人一二三区av| 久久毛片免费看一区二区三区| 99国产综合亚洲精品| 午夜老司机福利片| 免费在线观看完整版高清| 黄频高清免费视频| 美女扒开内裤让男人捅视频| 超碰成人久久| 美女中出高潮动态图| 国产熟女欧美一区二区| 大码成人一级视频| 欧美 亚洲 国产 日韩一| 欧美成人精品欧美一级黄| 午夜福利,免费看| 国产成人精品久久久久久| av片东京热男人的天堂| 不卡av一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 19禁男女啪啪无遮挡网站| 香蕉国产在线看| 两个人看的免费小视频| 晚上一个人看的免费电影| 国产淫语在线视频| 午夜日本视频在线| 人人澡人人妻人| 亚洲熟女毛片儿| 国产成人午夜福利电影在线观看| 1024视频免费在线观看| 天堂8中文在线网| 精品少妇一区二区三区视频日本电影 | 久久国产亚洲av麻豆专区| 性色av一级| 国产精品国产av在线观看| 母亲3免费完整高清在线观看| xxx大片免费视频| 超碰成人久久| 日韩 欧美 亚洲 中文字幕| 91精品三级在线观看| 91精品三级在线观看| 美女中出高潮动态图| 飞空精品影院首页| 人人澡人人妻人| 日韩一区二区三区影片| 尾随美女入室| 波多野结衣一区麻豆| xxxhd国产人妻xxx| 老熟女久久久| 满18在线观看网站| 大片电影免费在线观看免费| 韩国av在线不卡| 亚洲精品一二三| 一区二区三区精品91| 日本猛色少妇xxxxx猛交久久| 在线 av 中文字幕| 宅男免费午夜| 中文欧美无线码| 黄片播放在线免费| 熟妇人妻不卡中文字幕| 中文字幕人妻丝袜制服| 一区二区av电影网| 国产精品女同一区二区软件| 好男人视频免费观看在线| 欧美国产精品va在线观看不卡| 亚洲三区欧美一区| 国产精品免费视频内射| 国产片特级美女逼逼视频| 国产又爽黄色视频| 欧美日韩成人在线一区二区| 日日爽夜夜爽网站| 丁香六月天网| av网站免费在线观看视频| 国产精品三级大全| 国产精品久久久久久精品古装| 欧美黑人欧美精品刺激| 日韩制服骚丝袜av| 老司机亚洲免费影院| 国产又爽黄色视频| 亚洲七黄色美女视频| 又大又黄又爽视频免费| 日韩伦理黄色片| 日韩中文字幕视频在线看片| 日本爱情动作片www.在线观看| 色婷婷久久久亚洲欧美| 9热在线视频观看99| 亚洲成人一二三区av| 久久久精品免费免费高清| 国产精品久久久av美女十八| 老汉色av国产亚洲站长工具| 亚洲国产av新网站| 老司机深夜福利视频在线观看 | 色婷婷久久久亚洲欧美| 欧美老熟妇乱子伦牲交| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲最大av| 日本av手机在线免费观看| av电影中文网址| 综合色丁香网| 卡戴珊不雅视频在线播放| 热99国产精品久久久久久7| 国产片特级美女逼逼视频| 又黄又粗又硬又大视频| 一区福利在线观看| 成人亚洲精品一区在线观看| 日韩制服骚丝袜av| 少妇 在线观看| 赤兔流量卡办理| 国产精品国产三级专区第一集| 国产高清国产精品国产三级| 精品国产超薄肉色丝袜足j| 欧美成人精品欧美一级黄| 最近手机中文字幕大全| 多毛熟女@视频| 亚洲精品第二区| 午夜福利,免费看| 日韩伦理黄色片| 咕卡用的链子| 两个人免费观看高清视频| 这个男人来自地球电影免费观看 | 亚洲成av片中文字幕在线观看| 欧美日韩一级在线毛片| 久久久久网色| 91精品三级在线观看| 一区二区av电影网| 免费黄色在线免费观看| 夫妻性生交免费视频一级片| 中文字幕最新亚洲高清| 日日啪夜夜爽| 在线看a的网站| 亚洲美女视频黄频| 亚洲国产最新在线播放| 亚洲精华国产精华液的使用体验| 日韩制服骚丝袜av| 日韩一区二区三区影片| 欧美激情极品国产一区二区三区| 可以免费在线观看a视频的电影网站 | 亚洲成色77777| 18禁裸乳无遮挡动漫免费视频| 男人爽女人下面视频在线观看| 青春草视频在线免费观看| 男女无遮挡免费网站观看| 中国国产av一级| 日本av手机在线免费观看| 久久人人爽av亚洲精品天堂| 啦啦啦视频在线资源免费观看| 一二三四在线观看免费中文在| 欧美老熟妇乱子伦牲交| 国产成人欧美在线观看 | 人人妻,人人澡人人爽秒播 | 亚洲欧美成人综合另类久久久| 日韩免费高清中文字幕av| 日韩电影二区| 久久这里只有精品19| 五月天丁香电影| 两个人免费观看高清视频| 亚洲天堂av无毛| 精品久久久精品久久久| 一本一本久久a久久精品综合妖精| 国产极品粉嫩免费观看在线| 不卡av一区二区三区| 18禁国产床啪视频网站| 日本vs欧美在线观看视频| 黄色视频在线播放观看不卡| 精品人妻熟女毛片av久久网站| 日韩制服骚丝袜av| 精品国产超薄肉色丝袜足j| 国产男女超爽视频在线观看| √禁漫天堂资源中文www| 亚洲av综合色区一区| 国产av精品麻豆| 国产精品欧美亚洲77777| 纯流量卡能插随身wifi吗| 黄频高清免费视频| 久久 成人 亚洲| 免费少妇av软件| 超色免费av| 99九九在线精品视频| 成人午夜精彩视频在线观看| 日本一区二区免费在线视频| 老司机亚洲免费影院| 99re6热这里在线精品视频| e午夜精品久久久久久久| 777久久人妻少妇嫩草av网站| 老汉色av国产亚洲站长工具| 黑人猛操日本美女一级片| 久久久久久人人人人人| 久久免费观看电影| 欧美人与性动交α欧美软件| 成年人免费黄色播放视频| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜精品| 黄频高清免费视频| 国产成人a∨麻豆精品| 精品视频人人做人人爽| 精品一区二区三区四区五区乱码 | 在线亚洲精品国产二区图片欧美| 老鸭窝网址在线观看| 啦啦啦中文免费视频观看日本| 国产又爽黄色视频| 国产麻豆69| 一二三四在线观看免费中文在| 亚洲中文av在线| 美女扒开内裤让男人捅视频| 日日啪夜夜爽| 一区在线观看完整版| 欧美国产精品va在线观看不卡| 国产成人一区二区在线| 又粗又硬又长又爽又黄的视频| 久久精品久久久久久久性| 欧美日韩一区二区视频在线观看视频在线| 久久青草综合色| 日韩电影二区| 最近最新中文字幕免费大全7| 麻豆精品久久久久久蜜桃| 一本一本久久a久久精品综合妖精| 如日韩欧美国产精品一区二区三区| 高清av免费在线| 久久精品国产综合久久久| 超碰97精品在线观看| 欧美黑人精品巨大| 免费高清在线观看日韩| 69精品国产乱码久久久| 久久久久久久久免费视频了| 丝袜脚勾引网站| 尾随美女入室| 欧美国产精品一级二级三级| 欧美精品高潮呻吟av久久| 两性夫妻黄色片| 我要看黄色一级片免费的| 99国产精品免费福利视频| 中文精品一卡2卡3卡4更新| 蜜桃国产av成人99| 男女边摸边吃奶| 丝袜美足系列| 91精品三级在线观看| 2018国产大陆天天弄谢| 成年av动漫网址| 满18在线观看网站| 亚洲精品国产av成人精品| 久久av网站| 狂野欧美激情性xxxx| 女的被弄到高潮叫床怎么办| 电影成人av| 热re99久久国产66热| 国产淫语在线视频| 日韩中文字幕视频在线看片| 97在线人人人人妻| 成人国产麻豆网| 建设人人有责人人尽责人人享有的| 99精国产麻豆久久婷婷| 国产成人免费无遮挡视频| 亚洲精品久久成人aⅴ小说| 少妇人妻 视频| 黄网站色视频无遮挡免费观看| 黄片小视频在线播放| 欧美日本中文国产一区发布| 国产亚洲av高清不卡| 免费黄网站久久成人精品| 99国产精品免费福利视频| 最黄视频免费看| 亚洲国产最新在线播放| 色吧在线观看| 天堂8中文在线网| 久久久精品区二区三区| 欧美人与善性xxx| 亚洲美女黄色视频免费看| 日日摸夜夜添夜夜爱| 亚洲图色成人| 免费高清在线观看日韩| 亚洲精品日韩在线中文字幕| 男女下面插进去视频免费观看| 精品国产超薄肉色丝袜足j| 精品亚洲成国产av| 一级,二级,三级黄色视频| 久久久国产一区二区| 亚洲欧美精品综合一区二区三区| h视频一区二区三区| 色网站视频免费| 国产精品亚洲av一区麻豆 | 青春草视频在线免费观看| 1024视频免费在线观看| 免费久久久久久久精品成人欧美视频| 欧美日韩av久久| 赤兔流量卡办理| 如日韩欧美国产精品一区二区三区| 亚洲国产精品一区三区| 国产精品国产三级国产专区5o| 美女中出高潮动态图| 电影成人av| av一本久久久久| 午夜激情av网站| 亚洲第一青青草原| 亚洲av国产av综合av卡| 色视频在线一区二区三区| 日本欧美视频一区| 国产麻豆69| 欧美 日韩 精品 国产| 亚洲精品自拍成人| 下体分泌物呈黄色| 美女大奶头黄色视频| 夜夜骑夜夜射夜夜干| 午夜福利免费观看在线| 国产成人免费无遮挡视频| 蜜桃在线观看..| 亚洲国产中文字幕在线视频| 在线观看一区二区三区激情| 日韩中文字幕视频在线看片| 在线观看一区二区三区激情| 大码成人一级视频| 亚洲成人免费av在线播放| 日韩伦理黄色片| 一区二区三区激情视频| 国产极品天堂在线| 亚洲欧美日韩另类电影网站| 亚洲激情五月婷婷啪啪| 黄色视频在线播放观看不卡| 人成视频在线观看免费观看| 国产片特级美女逼逼视频| 精品国产乱码久久久久久男人| 国产成人免费无遮挡视频| 十八禁网站网址无遮挡| 制服人妻中文乱码| 国产不卡av网站在线观看| 久久久国产一区二区| 国产野战对白在线观看| 国产国语露脸激情在线看| 妹子高潮喷水视频| 亚洲精品aⅴ在线观看| 亚洲精品国产av成人精品| 国产男女超爽视频在线观看| 一级毛片我不卡| 日韩熟女老妇一区二区性免费视频| 天美传媒精品一区二区| 国产一区二区三区综合在线观看| 人人妻人人添人人爽欧美一区卜| 亚洲精品美女久久久久99蜜臀 | 色婷婷久久久亚洲欧美| 亚洲第一青青草原| 午夜老司机福利片| a级毛片黄视频| 午夜福利,免费看| 男人操女人黄网站| 日日摸夜夜添夜夜爱| 欧美精品高潮呻吟av久久| 亚洲美女黄色视频免费看| 国产av国产精品国产| 男女午夜视频在线观看| av网站免费在线观看视频| 免费观看性生交大片5| 精品国产一区二区三区四区第35| 青草久久国产| 性色av一级| 一区二区三区精品91| 亚洲三区欧美一区| 亚洲色图 男人天堂 中文字幕| 麻豆乱淫一区二区| 亚洲精品,欧美精品| 少妇人妻精品综合一区二区| a级毛片在线看网站| 成年人午夜在线观看视频| 成年美女黄网站色视频大全免费| 国产精品嫩草影院av在线观看| 免费观看人在逋| 久久精品国产综合久久久| 老司机影院毛片| 亚洲婷婷狠狠爱综合网| 你懂的网址亚洲精品在线观看| 国产在视频线精品| 成人毛片60女人毛片免费| 2021少妇久久久久久久久久久| 叶爱在线成人免费视频播放| 国产精品嫩草影院av在线观看| 久久99一区二区三区| a级毛片黄视频| 欧美精品一区二区大全| 99热全是精品| 午夜av观看不卡| 国产精品久久久久久精品古装| 国产精品偷伦视频观看了| 天堂俺去俺来也www色官网| 视频在线观看一区二区三区| 性色av一级| 成人黄色视频免费在线看| 999精品在线视频| 美女主播在线视频| 亚洲av成人不卡在线观看播放网 | 国产精品无大码| 丰满饥渴人妻一区二区三| 大香蕉久久成人网| 国产精品久久久久久久久免| 一区二区av电影网| 亚洲国产毛片av蜜桃av| 亚洲av男天堂| 无遮挡黄片免费观看| 精品久久蜜臀av无| h视频一区二区三区| 久久久久精品国产欧美久久久 | 日韩一区二区三区影片| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻,人人澡人人爽秒播 | 青春草亚洲视频在线观看| 免费在线观看完整版高清| 黄色怎么调成土黄色| 亚洲国产日韩一区二区| 欧美日韩国产mv在线观看视频| 肉色欧美久久久久久久蜜桃| 亚洲一区中文字幕在线| 男女边摸边吃奶| 丰满饥渴人妻一区二区三| 国产熟女午夜一区二区三区| www.自偷自拍.com| 久久久久视频综合| 秋霞伦理黄片| 色视频在线一区二区三区| 欧美激情极品国产一区二区三区| 一级毛片我不卡| 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久 | 搡老岳熟女国产| 久久国产精品大桥未久av| 熟女少妇亚洲综合色aaa.| 国产爽快片一区二区三区| 久久久亚洲精品成人影院| 亚洲婷婷狠狠爱综合网| 人人妻人人爽人人添夜夜欢视频| 99热国产这里只有精品6| 国产精品一区二区精品视频观看| 在线天堂最新版资源| 久久青草综合色| 另类亚洲欧美激情| 色吧在线观看| 黑人巨大精品欧美一区二区蜜桃| 午夜av观看不卡| 黄网站色视频无遮挡免费观看| 成人亚洲精品一区在线观看| √禁漫天堂资源中文www| 韩国精品一区二区三区| 黄色 视频免费看| 捣出白浆h1v1| 大片电影免费在线观看免费| 亚洲精品国产av蜜桃| 亚洲精品视频女| 一区二区三区乱码不卡18| 亚洲国产精品999| 亚洲av福利一区| 考比视频在线观看| 我的亚洲天堂| 99精国产麻豆久久婷婷| 老司机在亚洲福利影院| 国产黄频视频在线观看| 亚洲美女搞黄在线观看| 大码成人一级视频| 亚洲,欧美精品.| 高清不卡的av网站| 亚洲欧美精品自产自拍| 精品少妇一区二区三区视频日本电影 | 亚洲av福利一区| 久久久精品免费免费高清| 观看美女的网站| 五月开心婷婷网| 亚洲精品国产av成人精品| 老鸭窝网址在线观看| 久久久久久免费高清国产稀缺| 成人国产麻豆网| 免费少妇av软件| 久久午夜综合久久蜜桃| 人人妻人人澡人人爽人人夜夜| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久久久大奶| 飞空精品影院首页| videosex国产| av电影中文网址| 亚洲精品aⅴ在线观看| 韩国高清视频一区二区三区| 人人妻人人澡人人看| 久久国产精品男人的天堂亚洲| 久久久久人妻精品一区果冻| 久久久久久久久久久免费av| 免费日韩欧美在线观看| 亚洲精品,欧美精品| 黑人巨大精品欧美一区二区蜜桃| 一级片'在线观看视频| 高清视频免费观看一区二区| 午夜91福利影院| 一级片免费观看大全| 老司机靠b影院| 免费不卡黄色视频| 午夜福利一区二区在线看| 日本爱情动作片www.在线观看| 国产精品av久久久久免费| 亚洲欧美一区二区三区国产| 日本色播在线视频| 国产亚洲精品第一综合不卡| 色综合欧美亚洲国产小说| 亚洲男人天堂网一区| 亚洲国产最新在线播放| 欧美精品一区二区大全| 亚洲国产欧美日韩在线播放| 狂野欧美激情性xxxx| 久久亚洲国产成人精品v| 色播在线永久视频| 男女免费视频国产| 韩国高清视频一区二区三区| 亚洲av男天堂| 国产一卡二卡三卡精品 | 亚洲精品在线美女| 又大又爽又粗| 麻豆乱淫一区二区| 日韩av不卡免费在线播放| 国产精品99久久99久久久不卡 | 免费在线观看视频国产中文字幕亚洲 | 色婷婷av一区二区三区视频| 男人舔女人的私密视频| 1024视频免费在线观看| 中国三级夫妇交换| 男女午夜视频在线观看| e午夜精品久久久久久久| 久久久国产欧美日韩av| 亚洲av男天堂| 午夜福利视频在线观看免费| 国产成人欧美在线观看 | 国精品久久久久久国模美| 午夜福利,免费看| 成人18禁高潮啪啪吃奶动态图| 一区福利在线观看| 美女扒开内裤让男人捅视频| 日韩精品免费视频一区二区三区| av网站在线播放免费| 考比视频在线观看| √禁漫天堂资源中文www| 久久婷婷青草| 男女国产视频网站| 精品人妻在线不人妻| 久久ye,这里只有精品| 国产不卡av网站在线观看| 黄片播放在线免费| 久久鲁丝午夜福利片| 男男h啪啪无遮挡| 丝袜美腿诱惑在线| 在线 av 中文字幕| 色精品久久人妻99蜜桃| 黄网站色视频无遮挡免费观看| 桃花免费在线播放| a级片在线免费高清观看视频| 久久精品国产综合久久久| 十八禁人妻一区二区| 精品亚洲乱码少妇综合久久| 在线观看国产h片| 欧美日韩综合久久久久久| 麻豆乱淫一区二区| 色吧在线观看| 在线精品无人区一区二区三| 欧美国产精品va在线观看不卡| 婷婷色综合大香蕉| 一区福利在线观看| 午夜日韩欧美国产| 大香蕉久久成人网| 免费少妇av软件| 中文字幕av电影在线播放| a级毛片在线看网站| 国产乱人偷精品视频| 成人影院久久| 国产av码专区亚洲av| 国产精品久久久久久精品古装| 午夜激情久久久久久久| 女人久久www免费人成看片| 狂野欧美激情性bbbbbb| 亚洲av男天堂| 十八禁网站网址无遮挡| 午夜激情av网站| 国产成人91sexporn| a级片在线免费高清观看视频| 免费在线观看视频国产中文字幕亚洲 | 赤兔流量卡办理| 女人爽到高潮嗷嗷叫在线视频| 如日韩欧美国产精品一区二区三区| 亚洲国产看品久久| 午夜福利视频精品| av在线观看视频网站免费| 亚洲成人免费av在线播放| 久久精品久久精品一区二区三区| 天天操日日干夜夜撸| 女人被躁到高潮嗷嗷叫费观| 成人国产av品久久久| www.av在线官网国产| www日本在线高清视频| 少妇被粗大猛烈的视频| 国产精品国产av在线观看| 一本大道久久a久久精品| 中文字幕人妻熟女乱码| 国产精品国产av在线观看| 高清欧美精品videossex| 免费在线观看黄色视频的| 亚洲成人国产一区在线观看 | 免费少妇av软件| 男男h啪啪无遮挡| 久久精品亚洲av国产电影网| 美女中出高潮动态图| 欧美日韩国产mv在线观看视频| 亚洲色图 男人天堂 中文字幕| 精品人妻在线不人妻| 婷婷色综合www| 又大又黄又爽视频免费| 欧美av亚洲av综合av国产av | av在线老鸭窝| 成人亚洲欧美一区二区av| 国产精品久久久久成人av|