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

    氣候變化對燕山地區(qū)北蒼術(shù)適生區(qū)分布的影響

    2023-11-03 06:26:38王冶薛忠財楊越常佳寧趙春穎李云峰曾凡然
    中國中醫(yī)藥信息雜志 2023年11期
    關(guān)鍵詞:適生區(qū)燕山蒼術(shù)

    王冶 ,薛忠財 ,楊越 ,常佳寧 ,趙春穎 ,李云峰 ,曾凡然

    1.河北民族師范學院資源與環(huán)境科學學院,河北 承德 067000;2.河北省山區(qū)地質(zhì)環(huán)境重點實驗室,河北 承德 067000;3.冀北山區(qū)土壤與特色植物營養(yǎng)創(chuàng)新研究中心,河北 承德 067000;4.承德醫(yī)學院河北省中藥研究與開發(fā)重點實驗室,河北 承德 067000;5.承德筆架山生態(tài)農(nóng)業(yè)科技發(fā)展有限公司,河北 承德 067000

    物種分布區(qū)特征是物種重要的生態(tài)學和進化學特征[1],對制定物種保育策略具有重要參考價值[2],諸多因子均可能影響物種分布區(qū)特征[3-5],其中氣候變化已成為影響生態(tài)系統(tǒng)和生物類群的最主要因素之一[6]。氣候變化導致我國一些野生藥用植物資源面臨巨大威脅[7-8],預測未來氣候變化對物種地理分布的影響,有利于指導瀕危物種保育策略的制定。目前用于預測物種地理分布的常見模型已較為成熟[9],其中最大熵(MaxEnt)模型基于最大熵原理,對物種的潛在分布進行預測,因其預測效果好、預測覆蓋面積大等特點[10],被廣泛用于分析未來氣候?qū)ξ锓N分布的影響、野生動植物保護、物種適宜性等[11]。北蒼術(shù)Atractylodes chinensis(DC.) Koidz.是菊科蒼術(shù)屬多年生草本植物[12],其根莖是中藥蒼術(shù)的主要來源[13]。燕山地區(qū)所產(chǎn)北蒼術(shù)“朱砂點密,香氣濃郁”,因此作為該藥材的道地產(chǎn)區(qū)[14]。燕山地區(qū)地形地貌、生態(tài)系統(tǒng)相對復雜,氣候變化對北蒼術(shù)適生區(qū)的潛在影響尚不明確。因此,研究氣候變化影響下的北蒼術(shù)適生區(qū)分布規(guī)律十分必要。

    本研究基于實地調(diào)查和查閱文獻獲取燕山地區(qū)北蒼術(shù)地理分布信息,應(yīng)用MaxEnt模型和地理信息系統(tǒng)軟件ArcGIS,分析當前燕山地區(qū)北蒼術(shù)適生區(qū)分布,預測不同氣候情景下未來不同時期潛在適生區(qū)變化,為北蒼術(shù)引種馴化、資源保護和產(chǎn)業(yè)發(fā)展提供理論依據(jù)。

    1 資料與方法

    1.1 研究區(qū)概況與樣點分布

    本研究區(qū)域為廣義燕山地區(qū),位于壩上高原以南、河北平原以北、白河谷地以東、遼河平原以西的山地(93°41'~43°19'N、116°60'~112°13'E),調(diào)查區(qū)域總面積1.3×105km2,包含北京、張家口、唐山、承德、秦皇島、赤峰、葫蘆島、朝陽、錦州等市的全部或部分區(qū)域。燕山地區(qū)處于溫帶大陸性季風氣候區(qū),年均溫8.7~12.3 ℃,年降水量487.6 mm左右,作為京津冀地區(qū)的天然屏障和水源涵養(yǎng)地,海拔梯度較大,涵蓋森林、灌叢、草地、農(nóng)田、濕地等多樣的生境類型。

    野生北蒼術(shù)地理分布信息收集:①本團隊于2020年8月和2021年7月赴燕山地區(qū)進行實地考察,基于主要環(huán)境因子在燕山地區(qū)的空間分布特征進行布點,收集北蒼術(shù)地理分布信息共36個;②查閱中國數(shù)字植物標本館(https://www.cvh.ac.cn/)、中國國家標本資源平臺(http://www.nsii.org.cn/)及相關(guān)文獻[15-17],經(jīng)過數(shù)據(jù)篩選與核查,收集記錄信息河北12 個、內(nèi)蒙古5個、遼寧4個,利用ArcGIS10.0軟件,基于燕山地區(qū)范圍進行分布點信息校對,除去無效記錄,剩余13個。最終確定燕山地區(qū)野生北蒼術(shù)地理分布信息共49 個(見圖1)。本文涉及的地圖來自自然資源部標準地圖服務(wù)系統(tǒng)(http://bzdt.ch.mnr.gov.cn/)。

    圖1 燕山地區(qū)北蒼術(shù)樣點分布

    1.2 環(huán)境因子數(shù)據(jù)

    選擇生物氣候變量、地形地貌變量及土壤屬性變量等30個環(huán)境因子數(shù)據(jù)(見表1),19個溫度、降水相關(guān)的生物氣候變量數(shù)據(jù)[18]來自世界氣候數(shù)據(jù)庫(WorldClim,Version2.0,http://www.worldclim.org),包括當前時期(1970-2000 年)及未來3 個時期(2021-2040年、2041-2060年、2061-2080年),并基于聯(lián)合國政府間氣候變化專門委員會第六次評估報告,選擇全球氣候模型預測的溫室氣體低排放情景(SSP1-2.6)和高排放情景(SSP5-8.5)氣候數(shù)據(jù)。地形地貌變量中數(shù)字高程模型(DEM)數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺(http://www.gscloud.cn),坡度和坡向數(shù)據(jù)基于DEM,采用ArcGIS10.0軟件表面分析工具計算獲取。土壤屬性變量(0~30 cm)數(shù)據(jù)來自世界土壤屬性數(shù)據(jù)庫(https://geodata.pku.edu.cn/)。以上數(shù)據(jù)空間分辨率均統(tǒng)一為2.5 km。

    表1 30個環(huán)境因子

    1.3 模型構(gòu)建與精度評價

    采用MaxEnt3.4.1 軟件將北蒼術(shù)樣點分布信息(CSV格式)與不同時期的環(huán)境因子(ASCⅡ格式)數(shù)據(jù)輸入軟件中進行建模運算,刀切法計算各環(huán)境因子的貢獻率,隨機選擇75%數(shù)據(jù)作為訓練集,25%數(shù)據(jù)作為測試集,設(shè)置10次重復,同時設(shè)定最大迭代次數(shù)為10 000,其他參數(shù)為默認值[19]。利用ArcGIS10.0軟件進行作圖。采用受試者工作特征(ROC)曲線下面積(AUC)對預測結(jié)果進行評價[20]:AUC≤0.6表明預測結(jié)果失敗,0.6<AUC≤0.7表明預測結(jié)果較差,0.7<AUC≤0.8表明預測結(jié)果一般,0.8<AUC≤0.9表明預測結(jié)果良好,0.9<AUC≤1.0表明預測結(jié)果優(yōu)秀。AUC取值范圍為0~1,越接近1表明模型預測越準確[21]。

    1.4 適生區(qū)劃分

    基于模型結(jié)果,采用ArcGIS10.0軟件重分類工具對適生區(qū)進行等級劃分。根據(jù)概率P將適生區(qū)劃分為4個等級:高適生區(qū)(P≥0.8)、中適生區(qū)(0.6≤P<0.8)、低適生區(qū)(0.4≤P<0.6)和非適生區(qū)(P<0.4)[22]。

    1.5 時空變化分析

    采用重心遷移模型分析不同時期環(huán)境因子變化影響下,燕山地區(qū)北蒼術(shù)適生區(qū)重心遷移軌跡,從而反映其時空格局變化趨勢。重心遷移模型廣泛應(yīng)用在物種遷移、景觀格局演變等空間過程與格局的研究中[7,23-24],是反映對象空間演變過程的有效方法,其計算公式為:

    式中,zi為第i個平面空間單元的AUC值,(xi,yi)為第i個平面空間單元的投影坐標值,點(xˉ,yˉ)為北蒼術(shù)適生區(qū)地理分布重心的投影坐標值。重心移動的距離計算公式為:

    式中,dm表示地理分布重心移動的距離;yt+m和yt分別表示在時期t+m和t的地理分布重心投影y坐標;xt+m和xt分別表示在時期t+m和t的地理分布重心投影x坐標。

    2 結(jié)果與分析

    2.1 模型評價

    當前時期北蒼術(shù)MaxEnt模型AUC平均值為0.913;未來3個時期(2021-2040年、2041-2060年、2061-2080年)不同溫室氣體排放情景下,MaxEnt模型AUC平均值均大于0.9,表明模型預測結(jié)果優(yōu)秀,可用于北蒼術(shù)適生區(qū)分布的預測研究。見圖2。

    圖2 不同時期北蒼術(shù)生態(tài)分布預測模型AUC比較

    2.2 環(huán)境因子篩選

    影響北蒼術(shù)生長的主要環(huán)境因子為溫度季節(jié)性變化(46.6%),其次為高程、降水量季節(jié)性變化、最濕季度降水量、年均降水量,累積貢獻率達73.4%。主要環(huán)境因子適宜值和貢獻率見表2。

    表2 影響北蒼術(shù)分布主要環(huán)境因子適宜值與貢獻率

    2.3 適生區(qū)及其動態(tài)變化

    2.3.1 當前時期適生區(qū)分布

    當前時期燕山地區(qū)北蒼術(shù)適生區(qū)分布見圖3??傔m生區(qū)面積約7.54×104km2,占燕山地區(qū)總面積的59.10%,其中高適生區(qū)主要集中在承德西南部、張家口東部、秦皇島西北部和唐山北部,總面積約1.06×104km2,占燕山地區(qū)總面積的8.29%;中適生區(qū)主要集中在承德東部、秦皇島東部和赤峰西南部,總面積約2.12×104km2,占燕山地區(qū)總面積的16.61%;低適生區(qū)主要集中在朝陽西部、葫蘆島西部和北京北部,總面積約4.36×104km2,占燕山地區(qū)總面積的34.19%。綜上表明,燕山地區(qū)北蒼術(shù)適生區(qū)主要以承德、秦皇島和張家口為核心,高、中、低適生區(qū)相連并呈發(fā)散狀向東北方向分布。本研究北蒼術(shù)地理分布信息雖無位于燕山地區(qū)西北部的分布點,但根據(jù)模型預測結(jié)果可以看出,該區(qū)域也適合北蒼術(shù)生長。

    圖3 當前時期燕山地區(qū)北蒼術(shù)適生區(qū)分布

    2.3.2 未來氣候變化對適生區(qū)分布的影響

    未來不同時期2種氣候情景下燕山地區(qū)北蒼術(shù)適生區(qū)分布見圖4,適生區(qū)面積變化見表3。結(jié)果表明,在SSP1-2.6情景下,隨著時間變化,適生區(qū)面積增加的區(qū)域主要集中在朝陽中、西部和葫蘆島南部,而承德北部與赤峰南部適生區(qū)面積減少。2021-2040年適生區(qū)面積減少2 358.40 km2,占4.16%;2041-2060年適生區(qū)面積減少3 849.21 km2,占6.79%;2061-2080年適生區(qū)面積減少5 599.63 km2,占9.88%。在SSP5-8.5情景下,承德中北部與赤峰南部適生區(qū)面積隨時間增加,而朝陽、葫蘆島、秦皇島、唐山地區(qū)適生區(qū)面積均減少。2021-2040年適生區(qū)面積增加1 056.04 km2,占1.86%;2041-2060年適生區(qū)面積增加2 887.32 km2,占5.09%;2061-2080年適生區(qū)面積增加5 610.13 km2,占9.90%。綜上表明,SSP1-2.6情景導致的氣候變化使燕山地區(qū)北蒼術(shù)適生區(qū)面積呈退化趨勢,主要體現(xiàn)為高、中適生區(qū)向低、非適生區(qū)轉(zhuǎn)化,而SSP5-8.5情景導致的氣候變化使燕山地區(qū)北蒼術(shù)適生區(qū)呈擴張趨勢,主要體現(xiàn)為低、非適生區(qū)向高、中適生區(qū)轉(zhuǎn)化。

    表3 未來不同氣候情景下燕山地區(qū)北蒼術(shù)適生區(qū)面積轉(zhuǎn)移矩陣

    圖4 未來不同氣候情景下燕山地區(qū)北蒼術(shù)適生區(qū)變化

    2.4 氣候變化對適生區(qū)空間遷移的影響

    2.4.1 適生區(qū)水平分布變化

    未來2種氣候情景下,燕山地區(qū)北蒼術(shù)重心遷移路線見圖5。結(jié)果表明,當前適生區(qū)重心位于承德中部,SSP1-2.6情景下,適生區(qū)重心總體向燕山地區(qū)中、東部遷移,該區(qū)域為內(nèi)蒙古高原與遼河平原的過渡帶,最終遷至承德與朝陽交界處,各時期遷移距離分別為21.02、26.34、36.16 km,總體遷移距離為83.52 km。SSP5-8.5情景下,適生區(qū)重心總體向燕山地區(qū)北部遷移,該區(qū)域為山區(qū)與內(nèi)蒙古高原,最終遷至赤峰南部地區(qū),各時期遷移距離分別為26.27、27.27、46.94 km,總體遷移距離為100.44 km。

    圖5 不同氣候情景下燕山地區(qū)北蒼術(shù)適生區(qū)重心遷移路線

    2.4.2 適生區(qū)垂向分布變化

    提取2種氣候情景下,不同時期的燕山地區(qū)北蒼術(shù)適生區(qū)海拔數(shù)據(jù),進行高度為300 m的頻率分析。結(jié)果表明,各時期適生區(qū)面積垂向分布規(guī)律一致,隨海拔升高,適生區(qū)面積呈先增加后減少的特征,主要分布海拔區(qū)間為300~900 m。見圖6。

    圖6 不同氣候情景下燕山地區(qū)北蒼術(shù)適生區(qū)垂向分布變化

    對燕山地區(qū)北蒼術(shù)適生區(qū)垂向分布變化統(tǒng)計后發(fā)現(xiàn):到2061-2080年,SSP1-2.6情景下,適生區(qū)面積在海拔區(qū)間0~300 m 累積增加30.12%,在海拔區(qū)間300~600、600~900 m,適生區(qū)面積分別累積減少19.06%、7.56%,其他海拔區(qū)間變化不顯著;SSP5-8.5情景下,在海拔區(qū)間0~300 m,適生區(qū)面積累積減少24.96%,在海拔區(qū)間300~600、600~900 m,適生區(qū)面積累積增加10.14%、14.57%,其他海拔區(qū)間變化不顯著??傮w表明,在SSP1-2.6情景下,燕山地區(qū)北蒼術(shù)適生區(qū)垂向分布呈現(xiàn)由高海拔向低海拔區(qū)域遷移的趨勢,平均海拔下降88 m;在SSP5-8.5情景下,呈現(xiàn)由低海拔向高海拔區(qū)域遷移的趨勢,平均海拔上升119 m。

    3 討論

    根據(jù)環(huán)境因子貢獻率結(jié)果,對燕山地區(qū)北蒼術(shù)影響最大的氣候變量為溫度季節(jié)性變化,地形地貌變量主要為高程,二者累積貢獻率達58.6%。溫度季節(jié)性變化在1 100~1 350為最適宜條件,即全年4個季度溫度標準偏差在11~13.5 ℃最有利于北蒼術(shù)生長,與標本資料、實地調(diào)查獲取及相關(guān)文獻[25-26]基本一致;北蒼術(shù)最適海拔為300~900 m,此結(jié)果與中國數(shù)字植物標本館統(tǒng)計數(shù)據(jù)及《中國植物志》[12]記載較為一致。燕山地區(qū)屬溫帶大陸性季風氣候區(qū),具有鮮明的季度氣溫差異,且由于坡向、海拔等地理因素[27],氣溫呈現(xiàn)顯著的東西差異,西部地區(qū)溫度季節(jié)性變化為1 050~1 457,東部地區(qū)溫度季節(jié)性變化為890~980,表明西部地區(qū)具有更大的季節(jié)性氣溫變化。且燕山地區(qū)的東、西部地形地貌差異顯著,西部地區(qū)以高原山地為主,海拔在200~2 100 m,東部以丘陵平原為主,海拔30~180 m,總體呈現(xiàn)西高東低的地形特征。北蒼術(shù)適生區(qū)集中在燕山西部地區(qū),主要環(huán)境因子特征與燕山氣溫變化及地形特征高度吻合,表明模型計算的主要環(huán)境因子結(jié)果合理。

    當前條件下,燕山地區(qū)北蒼術(shù)適生區(qū)主要集中在承德、秦皇島和張家口東部地區(qū),與文獻記載[15,28]一致。SSP1-2.6情景下,適生區(qū)分布有所退化,地理分布重心向東遷移,總體東向遷移至海拔較低的內(nèi)蒙古高原與遼河平原過渡帶;SSP5-8.5情景下,適生區(qū)分布有所擴張,地理分布重心向東北遷移,總體向更高海拔的內(nèi)蒙古高原方向遷移??赡苁菧厥覛怏w低排放情景下,低海拔地區(qū)溫度下降,達到了適生條件,導致北蒼術(shù)向海拔較低的平原方向遷移;溫室氣體高排放情景下,氣溫升高使北蒼術(shù)原適生區(qū)物種豐富度增加,出現(xiàn)生態(tài)位重疊,進而產(chǎn)生種間競爭,且高海拔山區(qū)溫度有所升高,導致北蒼術(shù)向海拔更高的高原方向遷移。趙文龍等[7]研究表明,氣候變暖使紅花綠絨蒿低海拔適生區(qū)出現(xiàn)退化,高海拔適生區(qū)面積擴大;郭彥龍等[29]研究發(fā)現(xiàn),氣候變暖使桃兒七適生區(qū)的平均海拔逐漸升高。

    本研究顯示,溫度季節(jié)性變化與高程是影響燕山地區(qū)北蒼術(shù)生長的主要環(huán)境因子,當前總適生區(qū)面積約7.54×104km2,適生區(qū)較為廣泛。未來隨著氣候變化,北蒼術(shù)分布范圍將出現(xiàn)向不同海拔地區(qū)遷移的現(xiàn)象,因此在北蒼術(shù)人工馴化和引種區(qū)域的海拔選擇上,應(yīng)注意氣候變化,防止適生區(qū)退化。本研究結(jié)果有助于了解氣候變化對燕山地區(qū)北蒼術(shù)適生區(qū)分布的影響,為北蒼術(shù)資源可持續(xù)利用及產(chǎn)業(yè)發(fā)展提供參考。

    猜你喜歡
    適生區(qū)燕山蒼術(shù)
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預測
    未來氣候條件下當歸適生區(qū)預測及時空變化分析
    蒼術(shù)行情疲軟運行
    燕山水庫
    氣候變化下瀕危植物半日花在中國的潛在分布
    巴拉圭瓜多竹適生區(qū)分布研究
    竹子學報(2018年2期)2018-10-25 02:02:12
    GC法同時測定不同種類蒼術(shù)中6種成分
    中成藥(2018年5期)2018-06-06 03:11:59
    關(guān)蒼術(shù)根莖中蒼術(shù)酮含量及變化規(guī)律研究
    燕山水庫
    家在燕山
    詩選刊(2015年6期)2015-10-26 09:47:34
    在线免费观看的www视频| 一区福利在线观看| 国产欧美日韩一区二区三| 国产免费一级a男人的天堂| 久久精品国产亚洲av香蕉五月| 久久99热这里只有精品18| 久久久久久国产a免费观看| 女生性感内裤真人,穿戴方法视频| 欧美日韩乱码在线| 麻豆国产av国片精品| 国产久久久一区二区三区| 亚洲国产欧美网| 每晚都被弄得嗷嗷叫到高潮| 黄片小视频在线播放| 18禁在线播放成人免费| 操出白浆在线播放| 99热精品在线国产| av在线蜜桃| 一进一出好大好爽视频| 一二三四社区在线视频社区8| 嫩草影院入口| 亚洲不卡免费看| 内射极品少妇av片p| 国产精品影院久久| 国产一区二区三区在线臀色熟女| 亚洲精品久久国产高清桃花| 欧美一级毛片孕妇| 国产午夜精品论理片| 亚洲国产高清在线一区二区三| 久久这里只有精品中国| 久久精品91无色码中文字幕| aaaaa片日本免费| 成年女人毛片免费观看观看9| 啦啦啦韩国在线观看视频| 国产av一区在线观看免费| 欧美高清成人免费视频www| 日本免费一区二区三区高清不卡| 麻豆国产97在线/欧美| 欧美色视频一区免费| 尤物成人国产欧美一区二区三区| 嫁个100分男人电影在线观看| 欧美另类亚洲清纯唯美| 少妇高潮的动态图| 国产精品久久久人人做人人爽| 亚洲七黄色美女视频| 最新在线观看一区二区三区| 91在线观看av| 999久久久精品免费观看国产| 欧美在线黄色| 丝袜美腿在线中文| 欧美极品一区二区三区四区| 欧美激情久久久久久爽电影| 久久久成人免费电影| 啪啪无遮挡十八禁网站| 国产又黄又爽又无遮挡在线| av在线蜜桃| 哪里可以看免费的av片| 非洲黑人性xxxx精品又粗又长| 97超级碰碰碰精品色视频在线观看| 成人欧美大片| 久久香蕉精品热| 国产一区在线观看成人免费| 一个人观看的视频www高清免费观看| 在线观看免费午夜福利视频| 香蕉久久夜色| 国内少妇人妻偷人精品xxx网站| 长腿黑丝高跟| 国产美女午夜福利| 国产免费男女视频| eeuss影院久久| 成人午夜高清在线视频| 亚洲天堂国产精品一区在线| 不卡一级毛片| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 老熟妇乱子伦视频在线观看| 欧美午夜高清在线| 亚洲人成电影免费在线| 一级毛片高清免费大全| 国产精品久久久久久亚洲av鲁大| 国产精品女同一区二区软件 | tocl精华| 国产亚洲精品av在线| 欧美日韩黄片免| 老汉色av国产亚洲站长工具| 99久久精品国产亚洲精品| 91九色精品人成在线观看| 日韩成人在线观看一区二区三区| 久久亚洲真实| 亚洲不卡免费看| 欧美区成人在线视频| 亚洲,欧美精品.| 18禁裸乳无遮挡免费网站照片| 国产精品 欧美亚洲| 免费无遮挡裸体视频| 91av网一区二区| 午夜亚洲福利在线播放| 美女 人体艺术 gogo| 精品一区二区三区av网在线观看| 18+在线观看网站| 日本五十路高清| 18禁在线播放成人免费| 亚洲精品乱码久久久v下载方式 | 久久久色成人| 看黄色毛片网站| 99久久九九国产精品国产免费| 日本 av在线| 超碰av人人做人人爽久久 | 亚洲精品亚洲一区二区| 一个人免费在线观看的高清视频| 亚洲无线在线观看| 一级作爱视频免费观看| 亚洲精品456在线播放app | 午夜视频国产福利| 成人三级黄色视频| 18禁国产床啪视频网站| tocl精华| 久久精品国产综合久久久| 伊人久久大香线蕉亚洲五| 久久久久久国产a免费观看| 久久亚洲真实| 欧美色欧美亚洲另类二区| 国产欧美日韩精品一区二区| 亚洲乱码一区二区免费版| 制服丝袜大香蕉在线| 亚洲精品日韩av片在线观看 | 午夜精品一区二区三区免费看| 成人特级黄色片久久久久久久| 国产99白浆流出| 美女黄网站色视频| e午夜精品久久久久久久| 国产免费男女视频| 欧美在线黄色| 国产淫片久久久久久久久 | 少妇高潮的动态图| АⅤ资源中文在线天堂| 免费搜索国产男女视频| 女警被强在线播放| 国产午夜精品论理片| 亚洲精品成人久久久久久| 女人十人毛片免费观看3o分钟| 成人18禁在线播放| 久久久久国产精品人妻aⅴ院| 免费看日本二区| 757午夜福利合集在线观看| 国产精品1区2区在线观看.| 欧美最新免费一区二区三区 | 精品国产三级普通话版| 亚洲成av人片免费观看| 国产日本99.免费观看| 一级毛片女人18水好多| 最好的美女福利视频网| 免费大片18禁| 男女那种视频在线观看| 久久人人精品亚洲av| 国产高清视频在线观看网站| 欧美黄色淫秽网站| 九九在线视频观看精品| 九九热线精品视视频播放| 国产单亲对白刺激| 久久6这里有精品| www.熟女人妻精品国产| 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 日日夜夜操网爽| 国产男靠女视频免费网站| 淫妇啪啪啪对白视频| 亚洲精品在线美女| 69人妻影院| 久久人妻av系列| 两个人的视频大全免费| 亚洲aⅴ乱码一区二区在线播放| 男女那种视频在线观看| 丁香欧美五月| 亚洲成av人片免费观看| 又黄又粗又硬又大视频| 男女床上黄色一级片免费看| 国产美女午夜福利| 午夜免费男女啪啪视频观看 | 蜜桃亚洲精品一区二区三区| 一区二区三区高清视频在线| 精品国产美女av久久久久小说| 久久午夜亚洲精品久久| 一个人看视频在线观看www免费 | 在线国产一区二区在线| 一a级毛片在线观看| 69人妻影院| 久久伊人香网站| 夜夜看夜夜爽夜夜摸| 麻豆国产97在线/欧美| 精品午夜福利视频在线观看一区| 久久精品国产亚洲av涩爱 | 亚洲国产精品成人综合色| 亚洲欧美精品综合久久99| 又紧又爽又黄一区二区| 国产黄色小视频在线观看| 国产色爽女视频免费观看| 日韩精品青青久久久久久| a级毛片a级免费在线| 色av中文字幕| 91久久精品电影网| 日日干狠狠操夜夜爽| 精品久久久久久成人av| 91字幕亚洲| 日韩av在线大香蕉| av专区在线播放| 一a级毛片在线观看| 毛片女人毛片| 欧美激情久久久久久爽电影| 亚洲狠狠婷婷综合久久图片| 黄色丝袜av网址大全| 久久久久久人人人人人| 最近最新免费中文字幕在线| 久久精品国产亚洲av涩爱 | 给我免费播放毛片高清在线观看| 欧美日韩国产亚洲二区| 色吧在线观看| 久久精品国产99精品国产亚洲性色| 全区人妻精品视频| 亚洲国产精品成人综合色| 国产色爽女视频免费观看| 亚洲自拍偷在线| 欧美黑人巨大hd| 国产av麻豆久久久久久久| 一个人免费在线观看的高清视频| 性色avwww在线观看| 级片在线观看| 一本一本综合久久| 美女免费视频网站| 欧美激情在线99| 日本一本二区三区精品| 精品电影一区二区在线| 偷拍熟女少妇极品色| 色av中文字幕| 在线a可以看的网站| 黄片小视频在线播放| 真人一进一出gif抽搐免费| 在线播放国产精品三级| 精品不卡国产一区二区三区| 婷婷六月久久综合丁香| 最后的刺客免费高清国语| 国产激情欧美一区二区| 欧美一区二区精品小视频在线| 亚洲狠狠婷婷综合久久图片| 国产日本99.免费观看| 欧美日本视频| 国产爱豆传媒在线观看| 国产一区二区激情短视频| 国产真实乱freesex| 国产精品电影一区二区三区| 欧美不卡视频在线免费观看| 国产不卡一卡二| 99国产综合亚洲精品| 国产麻豆成人av免费视频| 久久久精品大字幕| 精品人妻偷拍中文字幕| 久久婷婷人人爽人人干人人爱| 观看美女的网站| 国内精品久久久久久久电影| 看免费av毛片| 精品一区二区三区视频在线观看免费| avwww免费| 成人午夜高清在线视频| 蜜桃久久精品国产亚洲av| 日本黄色视频三级网站网址| 性色av乱码一区二区三区2| 欧美一级a爱片免费观看看| 日本 av在线| x7x7x7水蜜桃| 黄色成人免费大全| 国产伦人伦偷精品视频| 18+在线观看网站| 亚洲精品色激情综合| 色av中文字幕| 午夜影院日韩av| 欧美成狂野欧美在线观看| 网址你懂的国产日韩在线| av在线天堂中文字幕| 国内久久婷婷六月综合欲色啪| 伊人久久大香线蕉亚洲五| 午夜福利视频1000在线观看| 国产精品1区2区在线观看.| 国产一区二区亚洲精品在线观看| 久久久精品大字幕| 国产高清激情床上av| av女优亚洲男人天堂| 亚洲av二区三区四区| 亚洲av一区综合| 免费一级毛片在线播放高清视频| 变态另类丝袜制服| 亚洲av电影不卡..在线观看| 午夜福利18| 成年人黄色毛片网站| av福利片在线观看| 欧美成人a在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 乱人视频在线观看| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 婷婷丁香在线五月| 91av网一区二区| 一级黄片播放器| 色哟哟哟哟哟哟| 日日干狠狠操夜夜爽| 一级毛片高清免费大全| 老汉色av国产亚洲站长工具| 少妇的逼好多水| 18禁黄网站禁片午夜丰满| 国产极品精品免费视频能看的| 天天躁日日操中文字幕| 免费无遮挡裸体视频| 亚洲无线在线观看| 久久人妻av系列| eeuss影院久久| e午夜精品久久久久久久| 大型黄色视频在线免费观看| 好男人在线观看高清免费视频| 在线a可以看的网站| 91久久精品国产一区二区成人 | 日本与韩国留学比较| 亚洲欧美日韩高清专用| 身体一侧抽搐| 国产精品久久久久久人妻精品电影| 午夜视频国产福利| 欧美一区二区国产精品久久精品| 成人精品一区二区免费| 3wmmmm亚洲av在线观看| 女人被狂操c到高潮| 欧美一区二区精品小视频在线| 日韩欧美精品免费久久 | 久久99热这里只有精品18| 人妻久久中文字幕网| 欧美日本亚洲视频在线播放| 在线免费观看不下载黄p国产 | 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人综合色| 国产熟女xx| 女生性感内裤真人,穿戴方法视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品免费一区二区三区在线| 欧美极品一区二区三区四区| 热99在线观看视频| 日本在线视频免费播放| 久久国产乱子伦精品免费另类| 女警被强在线播放| 日韩精品青青久久久久久| 一本一本综合久久| 日韩欧美 国产精品| 麻豆国产av国片精品| 国产99白浆流出| 亚洲精华国产精华精| 老司机在亚洲福利影院| 国模一区二区三区四区视频| 亚洲成人免费电影在线观看| 一区二区三区免费毛片| 中文亚洲av片在线观看爽| 色播亚洲综合网| 国产探花极品一区二区| 69av精品久久久久久| 国产免费男女视频| 我要搜黄色片| 特级一级黄色大片| 99国产精品一区二区蜜桃av| 欧美+日韩+精品| 国产精品女同一区二区软件 | 丁香六月欧美| 国产成人福利小说| 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 内射极品少妇av片p| 观看美女的网站| 久久精品国产亚洲av涩爱 | www国产在线视频色| 波多野结衣高清无吗| 午夜精品久久久久久毛片777| 亚洲成人久久性| 嫁个100分男人电影在线观看| 熟女电影av网| 99在线视频只有这里精品首页| 亚洲中文字幕日韩| 最近在线观看免费完整版| 国内精品美女久久久久久| 1024手机看黄色片| 成人国产一区最新在线观看| 桃色一区二区三区在线观看| 久久人人精品亚洲av| 国产日本99.免费观看| 男女午夜视频在线观看| www国产在线视频色| 午夜福利18| 桃色一区二区三区在线观看| 久久人人精品亚洲av| 一级黄片播放器| 2021天堂中文幕一二区在线观| 国产精品亚洲av一区麻豆| 亚洲精品日韩av片在线观看 | 国产男靠女视频免费网站| 少妇人妻精品综合一区二区 | 国产91精品成人一区二区三区| 亚洲男人的天堂狠狠| eeuss影院久久| 蜜桃久久精品国产亚洲av| 黄色女人牲交| 国产av一区在线观看免费| 婷婷六月久久综合丁香| av黄色大香蕉| 天堂√8在线中文| 国产在线精品亚洲第一网站| 亚洲精品色激情综合| 看免费av毛片| 最近视频中文字幕2019在线8| 欧美乱色亚洲激情| 狠狠狠狠99中文字幕| 激情在线观看视频在线高清| 国产精品久久视频播放| 国产老妇女一区| 亚洲av免费高清在线观看| 三级国产精品欧美在线观看| 午夜激情福利司机影院| 一级a爱片免费观看的视频| 亚洲精品456在线播放app | 日韩成人在线观看一区二区三区| 亚洲第一欧美日韩一区二区三区| 午夜福利欧美成人| 亚洲不卡免费看| 国产精品免费一区二区三区在线| 夜夜躁狠狠躁天天躁| 免费高清视频大片| 亚洲av中文字字幕乱码综合| 欧美一区二区亚洲| h日本视频在线播放| 久久久久九九精品影院| 精品一区二区三区人妻视频| 99热6这里只有精品| 少妇人妻一区二区三区视频| www.熟女人妻精品国产| АⅤ资源中文在线天堂| 中文字幕高清在线视频| 制服丝袜大香蕉在线| 国产精品乱码一区二三区的特点| 成人无遮挡网站| 国产午夜福利久久久久久| 最近在线观看免费完整版| 男人舔奶头视频| www国产在线视频色| 3wmmmm亚洲av在线观看| 九色国产91popny在线| 天堂动漫精品| 国产激情偷乱视频一区二区| 性色av乱码一区二区三区2| h日本视频在线播放| 免费观看的影片在线观看| 久久久久久人人人人人| 少妇的丰满在线观看| 99久久精品热视频| 国产激情欧美一区二区| 天天添夜夜摸| 欧美不卡视频在线免费观看| 男人和女人高潮做爰伦理| 亚洲av成人不卡在线观看播放网| 国产精品亚洲av一区麻豆| 草草在线视频免费看| 天堂网av新在线| 国产一区二区三区在线臀色熟女| 69人妻影院| 国产精品99久久久久久久久| 色精品久久人妻99蜜桃| 老熟妇仑乱视频hdxx| e午夜精品久久久久久久| av欧美777| 国产成人a区在线观看| 亚洲av熟女| 成人精品一区二区免费| 久久亚洲精品不卡| 亚洲av电影不卡..在线观看| 2021天堂中文幕一二区在线观| 欧美色欧美亚洲另类二区| 精品人妻1区二区| 国产日本99.免费观看| 婷婷精品国产亚洲av在线| 成年人黄色毛片网站| 久久精品91蜜桃| 观看美女的网站| 在线国产一区二区在线| 九色成人免费人妻av| 欧美最黄视频在线播放免费| 亚洲,欧美精品.| 一级毛片高清免费大全| 日本a在线网址| 久久香蕉精品热| 国产成人影院久久av| 嫩草影视91久久| 国产成人av教育| 亚洲专区中文字幕在线| 一级作爱视频免费观看| 夜夜看夜夜爽夜夜摸| 啪啪无遮挡十八禁网站| 少妇熟女aⅴ在线视频| 国产一区二区三区在线臀色熟女| 国产精品一区二区免费欧美| 国产91精品成人一区二区三区| 级片在线观看| 国产一区二区激情短视频| 国产伦在线观看视频一区| 久久久久国内视频| 男人舔奶头视频| 国产精品久久电影中文字幕| 男女床上黄色一级片免费看| 国产免费男女视频| 午夜老司机福利剧场| 好男人电影高清在线观看| 欧美另类亚洲清纯唯美| 麻豆国产97在线/欧美| www日本在线高清视频| 桃红色精品国产亚洲av| 国产野战对白在线观看| 国产aⅴ精品一区二区三区波| 国产精品三级大全| 成人国产综合亚洲| 亚洲av免费高清在线观看| 国产午夜精品久久久久久一区二区三区 | 国产伦人伦偷精品视频| 国产一级毛片七仙女欲春2| 国产精品亚洲av一区麻豆| 一进一出好大好爽视频| 可以在线观看毛片的网站| 亚洲av一区综合| 精品国产超薄肉色丝袜足j| 美女高潮的动态| 亚洲,欧美精品.| 国产伦在线观看视频一区| 国模一区二区三区四区视频| 亚洲精品在线观看二区| 熟女少妇亚洲综合色aaa.| 欧美黑人欧美精品刺激| 国产伦一二天堂av在线观看| 最近最新免费中文字幕在线| 97碰自拍视频| 法律面前人人平等表现在哪些方面| 亚洲国产欧洲综合997久久,| 国产激情偷乱视频一区二区| 51午夜福利影视在线观看| 精品国产超薄肉色丝袜足j| 国产单亲对白刺激| 欧美在线黄色| av专区在线播放| 天天躁日日操中文字幕| 久久久久亚洲av毛片大全| 午夜福利成人在线免费观看| 中文字幕人妻丝袜一区二区| 国产免费男女视频| 一边摸一边抽搐一进一小说| 97碰自拍视频| 国产精品久久久久久久久免 | 国产精品自产拍在线观看55亚洲| 国产精品1区2区在线观看.| 操出白浆在线播放| avwww免费| 少妇的逼水好多| 夜夜爽天天搞| 大型黄色视频在线免费观看| 欧美黄色淫秽网站| 亚洲激情在线av| 怎么达到女性高潮| 亚洲av一区综合| av在线蜜桃| 亚洲久久久久久中文字幕| 淫秽高清视频在线观看| 成人国产综合亚洲| 国产欧美日韩一区二区精品| 精品国产三级普通话版| 免费av观看视频| 黄片大片在线免费观看| 少妇人妻一区二区三区视频| 高潮久久久久久久久久久不卡| 18禁美女被吸乳视频| 欧美最新免费一区二区三区 | 亚洲国产精品国产精品| 国产精品精品国产色婷婷| 午夜福利在线观看免费完整高清在| 嘟嘟电影网在线观看| 国产成人福利小说| 国产精品99久久久久久久久| 2021少妇久久久久久久久久久| eeuss影院久久| 亚洲在线自拍视频| 国精品久久久久久国模美| 亚洲精品乱码久久久久久按摩| av黄色大香蕉| 丰满乱子伦码专区| 免费看光身美女| 99久国产av精品国产电影| 看非洲黑人一级黄片| 精品久久久久久久人妻蜜臀av| av卡一久久| 男女啪啪激烈高潮av片| 亚洲久久久久久中文字幕| 国产黄a三级三级三级人| 亚洲成色77777| 天堂俺去俺来也www色官网 | 中文在线观看免费www的网站| 麻豆国产97在线/欧美| 免费大片黄手机在线观看| 99热6这里只有精品| 国产精品人妻久久久影院| 搞女人的毛片| 最近2019中文字幕mv第一页| 日韩欧美三级三区| 伦理电影大哥的女人| 国产淫片久久久久久久久| 免费播放大片免费观看视频在线观看| 国产精品不卡视频一区二区| 午夜福利视频1000在线观看| 国产不卡一卡二| 久久久久免费精品人妻一区二区|