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

    華中櫻桃適生區(qū)模擬和生態(tài)特征分析

    2022-11-29 13:52:02董京京王賢榮伊賢貴
    關鍵詞:適生區(qū)華中櫻桃

    董京京,陳 潔,楊 宏,李 蒙,王賢榮,伊賢貴

    (南京林業(yè)大學,南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心,南京林業(yè)大學生物與環(huán)境學院,南京林業(yè)大學櫻花研究中心,江蘇 南京 210037)

    氣候變化影響全球環(huán)境變化,而物種對環(huán)境需要作出相應的反應,因此氣候變化可以通過改變生態(tài)系統(tǒng)影響物種的地理分布[1]。氣候變化可以促進或抑制植物生長發(fā)育,進而影響其適生區(qū)。使用科學的方法研究氣候和植物的相互影響,探討在氣候變化下植物的分布規(guī)律,有利于更加深入了解氣候變化情況下植物生態(tài)系統(tǒng)的適應性,對植物的引種開發(fā)、資源布局和保護也有著重要的現(xiàn)實意義[2-3]。近年來,基于生態(tài)位的理念,物種分布模型隨著全球信息系統(tǒng)技術的快速發(fā)展在生物地理學和基礎生態(tài)學研究方面有重要作用,為在全球氣候變化下物種潛在分布區(qū)的預測提供了一種新途徑[4]。常用的物種分布模型有DOMAIN、MaxEnt、Bioclim、ENFA、MARS、ANN等[5-6],其中2004年由Phillips等[7]提出MaxEnt模型具有預測精度高、操作方便、使用性強等特點,被廣泛應用于研究物種潛在適生區(qū)的預測。

    華中櫻桃[Cerasusconradinae(Koehne) Yü et Li],隸屬于薔薇科(Rosaceae)櫻屬(CerasusMill.),為中國特有櫻屬植物資源,分布于中高海拔500~2 100 m區(qū)域的溝邊林中[8-9]。華中櫻桃樹形高大美觀,花期長,可觀花觀葉觀果,或作為造林生態(tài)樹種[10],還可作為優(yōu)秀的櫻花苗木砧木材料,具有極高的觀賞和利用價值[11]。但目前華中櫻桃處于半野生狀態(tài),尚待研究和發(fā)掘利用。在此之前,華中櫻桃的研究主要集中在親緣關系、資源調查、群落結構等方面:付濤等[12]運用 SSR分子標記技術對華中櫻桃的親緣關系進行探討;伊賢貴[13]對武夷山櫻屬資源進行調查,揭示了華中櫻桃在廣大的中部及東部都可直接引種栽植;黎錄松[14]對武夷山自然保護區(qū)華中櫻桃的群落結構進行研究,間接揭示了其在群落中為優(yōu)勢種但有漸漸衰退的趨勢。然而目前基于生態(tài)地理學方面華中櫻桃的環(huán)境適生區(qū)預測研究鮮見報道。本研究基于華中櫻桃在中國的標本記錄,結合野外分布居群的地理信息定位和氣候數(shù)據(jù),運用MaxEnt(Version 3.4.1)模型和Arc-GIS(Version 10.6)軟件重建華中櫻桃的地理分布格局,模擬華中櫻桃在中國當代、2050s、2070s潛在的適生分布區(qū)域,并進一步研究在全球氣候變化下氣候因子與其分布格局的相互關系,以期為華中櫻桃的資源利用和種質保護提供理論依據(jù)[15]。

    1 材料與方法

    1.1 華中櫻桃地理分布數(shù)據(jù)收集

    華中櫻桃地理分布信息的收集主要通過查詢物種標本數(shù)據(jù)庫,包括檢索國家標本資源共享平臺內共建的中國數(shù)字植物標本館(http://www.cvh.org.cn/)、全球生物多樣性信息網(http://www.gbif.org/)、有關文獻數(shù)據(jù)庫(如www.cnki.net)等,搜索到24所標本館的554份標本,用Excel記錄、整理。為獲取精確的標本信息,刪除了地理信息模糊以及人為種植或重復的標本信息,并準確逐一對照下載的標本信息數(shù)據(jù),確保標本地理位置準確無誤[16]。運用百度地圖對有明確地點但無經緯度的地名進行精準定位,為防止樣點過于集中而產生太多誤差,在每個2.5′×2.5′網格中只取1個離中心點距離最短的分布點,結合華中櫻桃目前的資源分布進行實地調查野外分布居群的地理信息定位數(shù)據(jù),最終搜集到有效的201個華中櫻桃地理分布點數(shù)據(jù)[17],并轉化成MaxEnt模型所需要的.csv格式保存。

    1.2 環(huán)境因子來源及主導環(huán)境因子分析

    環(huán)境因子數(shù)據(jù)主要下載于世界氣候數(shù)據(jù)庫WorldClim Version 2.1,January 2020(http://www.worldclim.org/),當代(1970—2000年)、2050s(2041—2060年)、2070s(2061—2080年)等時期的19個氣候因子(.tiff)文件,空間分辨率2.5′,未來氣候數(shù)據(jù)(2050s、2070s)采用通用大氣環(huán)流模型CCSM4,模擬發(fā)布了4種未來溫室氣體排放濃度,RCP指典型濃度路徑(representative concentration pathway),RCP 2.6、RCP 4.5、RCP 6.0、RCP 8.5分別代表溫室氣體排放最低、中低、中高、最高等情境[18]。下載后需將已獲得的全球氣候數(shù)據(jù)用Arc-GIS軟件把.tiff格式的氣候數(shù)據(jù)文件轉換成MaxEnt模型所能識別的.asc格式文件并用掩膜剪裁工具剪裁為中國范圍。

    根據(jù)田聰?shù)萚19]基于DIVA-GIS(Version 7.5)軟件獲取氣候信息的方法提取華中櫻桃實際有效分布點的氣候數(shù)據(jù),用SPSS 26.0軟件對基于華中櫻桃抽取的環(huán)境變量進行主成分分析,篩選出決定當下影響華中櫻桃地理分布的限制氣候因子[20-21],與各個華中櫻桃地理分布點的海拔、經緯度進行Pearson相關性分析。

    1.3 適生區(qū)預測模型構建

    將19個氣候因子(表1,按貢獻率列舉)的數(shù)據(jù).asc格式以及華中櫻桃標本信息數(shù)據(jù).csv格式(數(shù)據(jù)依次物種名、經度、緯度)導入MaxEnt模型,設置測試集25%的華中櫻桃分布點來建立模型,訓練集75%的華中櫻桃分布點來測試模型,勾選刀切法、繪制預測分布圖以及環(huán)境響應曲線,其他選項均保持默認,以重復運算10次來保證模型檢測的準確性。

    表1 MaxEnt模型預測各個環(huán)境變量的貢獻率

    以中國省級行政區(qū)劃矢量地圖為底圖,導入華中櫻桃地理分布數(shù)據(jù),繪制實際地理分布圖。中國省級行政區(qū)劃矢量地圖(地圖比例尺為1∶4 000 000)下載于國家基礎地理信息中心(http://www.ngcc.cn/ngcc/)。根據(jù)MaxEnt模型模擬結果,將數(shù)據(jù)導入Arc-GIS軟件中,并轉化成柵格數(shù)據(jù),再利用Reclassify工具,完成華中櫻桃當前、2050s以及2070s氣候情境下適生區(qū)分布的模擬預測。利用Arc-GIS軟件對華中櫻桃潛在適生區(qū)預測結果進行劃分,Arc-GIS軟件的等級劃分方法中選擇平均間距法(equal interval)將適宜度的邏輯值(p)從0~1劃分為5個等級:不適宜區(qū)(0,0.2]、低適宜區(qū)(0.2,0.4]、中適宜區(qū)(0.4,0.6]、高適宜區(qū)(0.6,0.8]、極高適生區(qū)(0.8,1.0][22-23]。

    1.4 預測模型精確度檢測方法

    利用受試者特征工作曲線(receiver operating characteristic curve,ROC)預測效果,該曲線以橫坐標為假陽性率(即1-特異率),縱坐標為真陽性率(即靈敏度,為1-遺漏率)繪制而成,與橫坐標圍成的面積值AUC(crea under receiver operating characteristic curve)來評價預測的模型精度結果,為1-遺漏率,其結果范圍為[0,1][24],參考王運生等[25]的方法設置模型檢測標準(以模型精度區(qū)分):失敗(0.5,0.6]、較差(0.6,0.7]、一般(0.7,0.8]、較好(0.8,0.9]、極好(0.9,1.0]。

    2 結果與分析

    2.1 華中櫻桃分布區(qū)預測模型的精確度檢測結果

    通過MaxEnt模型檢測的10組數(shù)據(jù),其中運算結果最大的一組訓練集AUC為0.945,測試集的為0.949,根據(jù)受試者工作曲線(ROC)的模型評價標準可知,MaxEnt模型基于氣候數(shù)據(jù)對華中櫻桃潛在分布區(qū)的檢測達到極好的標準(圖1)。

    圖1 MaxEnt 模型中ROC曲線預測結果

    2.2 華中櫻桃地理分布和當代潛在適生區(qū)預測結果

    通過檢索標本以及查閱文獻,表明華中櫻桃水平地理分布范圍為98°21′0″~119°33′36″E、22°43′48″~34°17′24″N,主要分布于西南、華中、華東地區(qū)的四川、貴州、湖南、重慶、江西、云南等地,東至浙江天目山一帶、西至云南怒江傈僳族自治州、南至云南普洱市西盟佤族自治縣、北至陜西雞峰山。

    根據(jù)華中櫻桃的201個地理分布點數(shù)據(jù)以及19個環(huán)境因子數(shù)據(jù),使用MaxEnt模型和Arc-GIS軟件預測其當代潛在適生區(qū)(圖2)。由圖2可知,華中櫻桃在當代環(huán)境因子影響下潛在分布區(qū)范圍為93°~122°E、22°~35°N,總適宜區(qū)面積預測共170.00萬km2。其中極高適宜區(qū)面積1.46萬km2,主要集中在重慶武隆縣、忠縣,面積占適宜區(qū)比重最小。高適宜區(qū)面積共28.43萬km2,分布于兩個核心區(qū)域(0.6

    底圖審圖號:GS(2019)1822。下同。

    2.3 華中櫻桃適生區(qū)主導環(huán)境因子分析

    由MaxEnt最大熵模型輸出(表1)可知,影響華中櫻桃潛在分布區(qū)的環(huán)境貢獻率依次為bio 12、bio 7、bio 4、bio 17,累計貢獻率達79.2%;置換重要值排前4的依次是bio 12、bio 4、bio 9、bio 19。

    利用刀切法對環(huán)境主導因子分析可知(圖3)正規(guī)化訓練增益值前4依次是bio 7、bio 12、bio 6、bio 11,測試增益值前4依次是bio 7,bio 16、bio 12、bio 18,AUC值前4依次是bio 16、bio 7、bio 12、bio 17,綜上所述可知影響華中櫻桃適生區(qū)分布的環(huán)境因子有bio 4、bio 6、bio 7、bio 9、bio 11、bio 12、bio 16、bio 17、bio 18、bio 19。

    藍色表示單一變量,綠色表示除該變量外的其他變量組合,紅色表示所有變量。Blue represents a single variable, green represents combinations of variables other than that variable, and red represents all variables.

    通過DIVA-GIS軟件獲取10個環(huán)境變量實際地點的氣候數(shù)據(jù),對其進行主成分分析,結果表明:前3個主成分方差貢獻率分別是45.22%、30.68%、16.38%,累計貢獻率達92.27%,集中了10個氣候變量的絕大多數(shù)信息(表2)。由前3個主成分相對于10個氣候變量的得分系數(shù)(表3)可知,在第1主成分載荷較大的是最濕季降水量(bio 16)、年降水量(bio 12),可代表降水因子對地理分布的主要影響,在第2主成分載荷較大是溫度季節(jié)變化方差(bio 4)與溫度年較差(bio 7),可代表溫度因子對地理分布的重要影響(表3)。

    表2 主成分前3個方差解釋

    表3 前3個主成分相對于10個氣候變量的得分系數(shù)

    2.4 華中櫻桃分布范圍主導氣候因子與海拔及經緯度的相關性

    由MaxEnt模型繪制華中櫻桃存在概率與bio 16、bio 12、bio 4、bio 7的響應曲線結果見圖4。結果表明,將華中櫻桃存在概率大于0.5定為環(huán)境變量的適宜范圍,則最濕季降水量(bio 16)范圍為500~800 mm,年均降水量(bio 12)的范圍為1 050~1 700 mm,溫度季節(jié)變化方差(bio 4)的范圍為650~840(標準差×100),溫度年較差(bio 7)的范圍為26.0~30.5 ℃。

    圖4 華中櫻桃存在概率與主導環(huán)境因子響應曲線

    對主導限制氣候因子與海拔、經緯度進行Pearson相關性分析,通過單尾檢測相關系數(shù)(R)絕對值大于0.50結果表明:bio 7與bio 4及緯度呈極顯著正相關(P< 0.01),相關系數(shù)分別為0.806、0.601;bio 16與bio 12呈極顯著正相關,與緯度呈極顯著負相關(P< 0.01),相關性系數(shù)分別為0.820、-0.516;bio 4與經度呈極顯著正相關,與海拔呈極顯著負相關(P< 0.01),相關系數(shù)分別為0.685、-0.601;bio 12與經度呈極顯著正相關(P< 0.01),相關系數(shù)為0.756(表4)。

    2.5 未來氣候變化下華中櫻桃潛在適生區(qū)預測

    選取2050s、2070s數(shù)據(jù)采用通用大氣環(huán)流模型CCSM4對華中櫻桃潛在適生區(qū)進行模擬,得到不同氣候情境下潛在適生區(qū)的分布面積(表5)。

    表5 不同氣候情境下華中櫻桃的適生面積

    由表5可見,在2050s RCP 2.6、RCP 4.5溫室氣體排放濃度較低情境下,華中櫻桃總適宜區(qū)面積呈現(xiàn)先減小后逐漸增大狀態(tài),且在RCP 4.5濃度時達到總適宜區(qū)面積最大狀態(tài),在RCP 6.0、RCP 8.5溫室氣體排放濃度較高時總適宜區(qū)面積逐漸減小,整體而言,2050s RCP 8.5氣候情境下的總適宜區(qū)面積比當代總適宜區(qū)面積減小1.63萬km2,高適宜區(qū)和極高適宜區(qū)整體面積增加,低適宜區(qū)和中適宜區(qū)整體面積減小。從表5、圖2、圖5可知,極高適宜區(qū)在2050s逐漸向東北方向擴張,在RCP 2.6情況下極高適宜區(qū)核心區(qū)域主要集中在湖北恩施、鶴峰、咸豐、巴東附近等地;在RCP 4.5氣候情境下,極高適宜區(qū)面積比當代面積增加最多,可達3.22萬km2;在RCP 6.0氣候情境下,四川雅安附近區(qū)域出現(xiàn)一定范圍的極高適生區(qū),且有逐漸向高緯度地區(qū)遷移的趨勢;在RCP 8.5氣候情境下,極高適宜區(qū)核心區(qū)域除湖北外其他部分更破碎化,重慶、貴州、湖南、江西、福建、浙江、安徽均有分布。高適宜區(qū)面積在不同氣候情境下穩(wěn)定增加,RCP 4.5時面積整體最大,相比于當代,面積增加10.54萬km2,分布趨勢逐漸向低緯度擴展,覆蓋重慶和貴州絕大多數(shù)區(qū)域。中適宜區(qū)和低適宜區(qū)與當代分布地區(qū)沒有明顯變化,基本趨于穩(wěn)定狀態(tài)。

    在2070s RCP 4.5氣候情境下華中櫻桃總適宜區(qū)相比于當代面積減小,其他氣候情境下總適宜區(qū)面積皆穩(wěn)定增加,在RCP 6.0氣候情境下總適宜區(qū)面積達到最大,比當代面積增加10.03萬km2。極高適宜區(qū)和高適宜區(qū)隨著氣候的變暖面積皆穩(wěn)定增長。中適宜區(qū)在RCP 2.6氣候情境下較當代面積增大,其他氣候情境下面積均小于當代。低適宜區(qū)在RCP 8.5時面積比當代小,其他氣候情境下面積均大于當代。由表5、圖2、圖5可知,在2070s氣候條件下極高適宜區(qū)在RCP 6.0氣候情境時面積最大,相較當代,增加近4.3倍且依然有向東北方向擴張的趨勢,高適宜區(qū)在RCP 2.6氣候情境下相比于當代逐漸向重慶周圍聚攏且在浙江和福建交界處出現(xiàn)一定面積,在RCP 8.5氣候情境下高適宜區(qū)繼續(xù)向四周偏向低緯度擴張,面積最大比當代增加12.46萬km2。

    圖5 基于MaxEnt模型在2050s和2070s RCP 2.6、RCP 4.5、RCP 6.0、RCP 8.5氣候情境下華中櫻桃潛在適宜區(qū)

    3 討 論

    華中櫻桃為中國特有櫻屬植物,分布范圍跨度廣,溫帶和亞熱帶均有分布,適應氣候能力強[26]。筆者通過MaxEnt模型結合Arc-GIS地理信息系統(tǒng)對華中櫻桃的生態(tài)位進行模擬,基于地理信息技術DIVA-GIS提取出潛在分布區(qū)所在生境的10個環(huán)境因子,利用PCA主成分分析和相關性分析,去除變量因子之間的共線性影響,最終明確了最濕季降水量(bio 16)、年降水量(bio 12)是限制華中櫻桃當代潛在地理分布的首要因子,溫度季節(jié)變化方差(bio 4)和溫度年較差(bio 7)是重要的氣候因子,與前人研究櫻屬植物的生態(tài)地理學結論具有相似性。如王華辰等[27]通過研究中國特有雪落櫻桃(Cerasusxueluoensis)[28]潛在適生區(qū)和氣候特征的關系時對其生態(tài)適應性特征進行定性定量分析表明降水是主導氣候因子;朱弘等[20]通過研究浙閩櫻桃(Cerasusschneideriana)地理分布的氣候限制因子時表明年降水量(bio 12)、最濕季降雨量(bio 16)、最暖季降雨量(bio 18)為代表的水分因子對其分布有主要影響。華中櫻桃的主導氣候限制因子水分相比溫度影響更大,推測與華中櫻桃分布區(qū)氣候帶受北亞熱帶和中亞熱帶濕潤季風區(qū)的影響有關。

    朱淑霞等[29]通過預測當代氣候變化下高盆櫻桃(Cerasuscerasoides)和鐘花櫻桃(Cerasuscampanulata)的地理分布,結果表明分布中心為云貴高原、武夷山脈、南嶺山脈,在未來氣候(CCM3)變暖情況下,適生區(qū)面積將減小。張興旺等[30]通過研究氣候變化對黃山花楸(Sorbusamabilis)潛在分布區(qū)的影響時發(fā)現(xiàn)在未來不同氣候情境下各個適宜區(qū)面積均較當代面積減小且有向高海拔區(qū)域遷移趨勢。本研究結果表明,華中櫻桃在當代潛在適生區(qū)分布于重慶、貴州、湖南、四川、武漢、浙江、福建、江西、安徽、河南、云南、廣西、廣州、西藏以及中國臺灣,在江蘇、山東和上海少量地區(qū)也有分布趨勢,核心分布區(qū)主要集中于西南和華中地區(qū)。未來分布區(qū)面積隨著氣候逐漸變暖,在2050s總適宜區(qū)相較于當代面積呈減少狀態(tài),但在2070s總適宜面積增加,且極高適宜區(qū)有向高緯度和東北方向擴散的趨勢,高適宜區(qū)有向低緯度方向擴散的趨勢,但核心適宜區(qū)依然保持在西南和華中地區(qū)。江志紅等[31]對21世紀氣候變化(IPCC-AR4)進行評估表示溫度增幅在21世紀中葉(2021—2050年)不同氣候情境下差異不大,到21世紀末期(2071—2100年)預估氣候變暖差異大,中國區(qū)域年平均溫度升高2.5 ℃,降水量增大12%;姜大膀等[32]對東亞氣候情境預測研究時揭示溫度增幅在東北、西部、華中地區(qū)表現(xiàn)出明顯變化,同時導致大陸降水全區(qū)域性增大。綜上預測,在2070s中國區(qū)域尤其以華中、西南地區(qū)降雨量增加,溫度增幅變化明顯,華中櫻桃在2070s總適宜區(qū)面積增大的預測結果與上述氣候綜合因子的變化趨勢相關。本研究基于收集華中櫻桃地理分布點數(shù)據(jù)對其潛在適生區(qū)和生態(tài)特征進行分析,為華中櫻桃的種質保護、親緣地理與資源利用的后續(xù)研究奠定了重要基礎。

    猜你喜歡
    適生區(qū)華中櫻桃
    氣候變化下中國蒟蒻薯科箭根薯的地理分布格局預測
    未來氣候條件下當歸適生區(qū)預測及時空變化分析
    華中要塞:義陽三關
    櫻桃肉,讓年味飄香
    金橋(2022年2期)2022-03-02 05:43:04
    華中建筑2021年總目錄
    華中建筑(2021年12期)2022-01-17 02:08:42
    新四軍華中抗戰(zhàn)
    氣候變化下瀕危植物半日花在中國的潛在分布
    明年或激增40%?華中3萬多噸加州鱸市場誰能笑到最后?
    當代水產(2019年2期)2019-05-16 02:46:14
    巴拉圭瓜多竹適生區(qū)分布研究
    竹子學報(2018年2期)2018-10-25 02:02:12
    一樹櫻桃?guī)в昙t
    人大建設(2018年5期)2018-08-16 07:09:08
    亚洲丝袜综合中文字幕| 成人午夜精彩视频在线观看| 久久久久久人人人人人| 国产精品熟女久久久久浪| 成人影院久久| 男女国产视频网站| 亚洲综合色网址| 亚洲欧美成人综合另类久久久| 26uuu在线亚洲综合色| av免费观看日本| 亚洲天堂av无毛| 这个男人来自地球电影免费观看 | 黄网站色视频无遮挡免费观看| 亚洲成人一二三区av| 久久精品国产a三级三级三级| 久久久久国产网址| 最新中文字幕久久久久| 看免费av毛片| 亚洲美女搞黄在线观看| 免费av不卡在线播放| 观看av在线不卡| 男女高潮啪啪啪动态图| 九色亚洲精品在线播放| 日韩电影二区| 日韩一本色道免费dvd| 国产男人的电影天堂91| 色婷婷久久久亚洲欧美| 国产日韩欧美亚洲二区| 亚洲第一区二区三区不卡| 日本vs欧美在线观看视频| 黄色毛片三级朝国网站| 久久久久久人人人人人| 大香蕉久久网| 纵有疾风起免费观看全集完整版| 国产男女超爽视频在线观看| av天堂久久9| 看十八女毛片水多多多| 午夜免费观看性视频| 免费黄频网站在线观看国产| 自线自在国产av| 亚洲伊人色综图| 亚洲美女搞黄在线观看| 久久精品国产a三级三级三级| 男女无遮挡免费网站观看| 久热这里只有精品99| 色5月婷婷丁香| 久久久久精品人妻al黑| 在线精品无人区一区二区三| 国产熟女欧美一区二区| 国产精品 国内视频| 久久久国产一区二区| 国内精品宾馆在线| 国产精品国产av在线观看| 一级片免费观看大全| 国产亚洲欧美精品永久| 一二三四在线观看免费中文在 | 99热国产这里只有精品6| 久久99热这里只频精品6学生| 国产69精品久久久久777片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 18在线观看网站| 九色成人免费人妻av| 国产精品嫩草影院av在线观看| 又黄又爽又刺激的免费视频.| 国产av码专区亚洲av| 美女国产视频在线观看| 热99国产精品久久久久久7| 丝袜喷水一区| 美女大奶头黄色视频| 国产男女内射视频| 麻豆精品久久久久久蜜桃| kizo精华| 看非洲黑人一级黄片| 久久97久久精品| 亚洲国产欧美日韩在线播放| 免费高清在线观看日韩| 国产成人精品婷婷| 久久人人爽人人片av| 亚洲国产精品999| 大香蕉久久成人网| 国产av国产精品国产| 女人被躁到高潮嗷嗷叫费观| 九九在线视频观看精品| 国产激情久久老熟女| 久久久久精品性色| 蜜桃在线观看..| 丝袜人妻中文字幕| 2018国产大陆天天弄谢| 两个人看的免费小视频| 夫妻午夜视频| 777米奇影视久久| 亚洲av福利一区| 精品一区二区免费观看| 在线观看免费日韩欧美大片| 国产女主播在线喷水免费视频网站| 亚洲一区二区三区欧美精品| 男女下面插进去视频免费观看 | 热re99久久国产66热| 亚洲国产精品一区三区| av不卡在线播放| 成人二区视频| 久久精品aⅴ一区二区三区四区 | 久久久久精品性色| 国产成人欧美| 免费看不卡的av| 大香蕉久久网| 视频在线观看一区二区三区| 日韩精品有码人妻一区| a级毛片在线看网站| 久久精品国产亚洲av涩爱| 久久久国产欧美日韩av| 亚洲国产欧美日韩在线播放| 青春草国产在线视频| 国产 精品1| 亚洲国产精品成人久久小说| 免费看av在线观看网站| 在线观看人妻少妇| 在线观看人妻少妇| 深夜精品福利| 国产av精品麻豆| 亚洲av国产av综合av卡| 精品国产一区二区久久| 国产精品三级大全| 在线观看免费高清a一片| 国产精品99久久99久久久不卡 | 亚洲精品日本国产第一区| 国产又色又爽无遮挡免| 亚洲国产av影院在线观看| 人妻人人澡人人爽人人| 国产伦理片在线播放av一区| 最新中文字幕久久久久| 国产亚洲精品久久久com| 考比视频在线观看| 精品一区二区三区四区五区乱码 | 欧美最新免费一区二区三区| 777米奇影视久久| 丝袜美足系列| 成人综合一区亚洲| 人妻 亚洲 视频| 男女午夜视频在线观看 | 精品99又大又爽又粗少妇毛片| 成人二区视频| 亚洲精品456在线播放app| 国产精品蜜桃在线观看| 亚洲精品456在线播放app| 18禁观看日本| 精品99又大又爽又粗少妇毛片| 高清毛片免费看| 丝瓜视频免费看黄片| 精品99又大又爽又粗少妇毛片| 97在线人人人人妻| av黄色大香蕉| 国产xxxxx性猛交| 成年女人在线观看亚洲视频| 免费人成在线观看视频色| 少妇被粗大的猛进出69影院 | 欧美精品人与动牲交sv欧美| 纵有疾风起免费观看全集完整版| 91午夜精品亚洲一区二区三区| 男女啪啪激烈高潮av片| 1024视频免费在线观看| 久久99精品国语久久久| 看非洲黑人一级黄片| 亚洲国产精品专区欧美| 亚洲内射少妇av| 成人漫画全彩无遮挡| 亚洲 欧美一区二区三区| 丰满少妇做爰视频| 午夜91福利影院| a级毛片黄视频| 三级国产精品片| 欧美日韩亚洲高清精品| 精品人妻一区二区三区麻豆| 午夜福利乱码中文字幕| 久久精品久久精品一区二区三区| 午夜福利视频精品| 美国免费a级毛片| 曰老女人黄片| 国产乱来视频区| 婷婷色av中文字幕| av在线app专区| 男人添女人高潮全过程视频| 午夜免费男女啪啪视频观看| 伦精品一区二区三区| 性高湖久久久久久久久免费观看| 成年美女黄网站色视频大全免费| av网站免费在线观看视频| 日产精品乱码卡一卡2卡三| 观看av在线不卡| 成人国语在线视频| 日韩av不卡免费在线播放| 精品人妻偷拍中文字幕| 日本av免费视频播放| 日韩精品有码人妻一区| 亚洲伊人久久精品综合| 免费女性裸体啪啪无遮挡网站| 波野结衣二区三区在线| 亚洲人成网站在线观看播放| 91精品三级在线观看| 夫妻性生交免费视频一级片| 日本免费在线观看一区| 国产福利在线免费观看视频| 国产精品国产三级专区第一集| 日韩伦理黄色片| 国产精品国产三级国产av玫瑰| 少妇的逼好多水| 亚洲精品久久午夜乱码| 最后的刺客免费高清国语| 久久久久网色| 一边摸一边做爽爽视频免费| 欧美 日韩 精品 国产| 欧美精品高潮呻吟av久久| 国产又色又爽无遮挡免| 国产一区二区三区综合在线观看 | 色视频在线一区二区三区| 国产精品久久久久久av不卡| 男人爽女人下面视频在线观看| 国产亚洲一区二区精品| 国产在线一区二区三区精| 秋霞伦理黄片| 大香蕉久久网| 久久这里只有精品19| 毛片一级片免费看久久久久| 国产欧美日韩一区二区三区在线| 亚洲精品一区蜜桃| 亚洲国产看品久久| 欧美人与善性xxx| 久久国产精品大桥未久av| 亚洲av国产av综合av卡| 亚洲欧洲国产日韩| 中文精品一卡2卡3卡4更新| 中国国产av一级| 丝袜美足系列| 国产精品一区www在线观看| 香蕉精品网在线| 国产一区有黄有色的免费视频| 日韩成人伦理影院| 精品国产露脸久久av麻豆| 久久久亚洲精品成人影院| 99香蕉大伊视频| 黑人高潮一二区| 人人澡人人妻人| 一本大道久久a久久精品| 黄网站色视频无遮挡免费观看| av黄色大香蕉| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品一区二区三区视频在线| 国产精品久久久久久久久免| 欧美精品高潮呻吟av久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜精品国产一区二区电影| 亚洲国产毛片av蜜桃av| 亚洲精品美女久久久久99蜜臀 | 免费高清在线观看视频在线观看| 最近最新中文字幕免费大全7| 国产精品一区二区在线不卡| 精品久久久精品久久久| 99热这里只有是精品在线观看| 欧美精品高潮呻吟av久久| 亚洲在久久综合| 午夜福利在线观看免费完整高清在| 亚洲天堂av无毛| 夫妻午夜视频| 日韩免费高清中文字幕av| 久久久久久人人人人人| 亚洲国产看品久久| 久久久久国产网址| 国产精品 国内视频| 国产成人免费观看mmmm| 久久综合国产亚洲精品| 女性被躁到高潮视频| 老熟女久久久| 久久久久久久久久成人| 高清黄色对白视频在线免费看| 精品国产一区二区久久| 成人18禁高潮啪啪吃奶动态图| 国产精品嫩草影院av在线观看| 99久久人妻综合| 国产毛片在线视频| 美女大奶头黄色视频| 桃花免费在线播放| 亚洲精品av麻豆狂野| 最近的中文字幕免费完整| 国产高清三级在线| 国产白丝娇喘喷水9色精品| 久久99精品国语久久久| 亚洲综合色网址| 国产精品久久久久久av不卡| 国产无遮挡羞羞视频在线观看| 久久国产精品男人的天堂亚洲 | 制服人妻中文乱码| 一本大道久久a久久精品| 久久久久久人人人人人| 色5月婷婷丁香| 赤兔流量卡办理| 纵有疾风起免费观看全集完整版| 欧美精品人与动牲交sv欧美| 亚洲精品乱码久久久久久按摩| 少妇被粗大猛烈的视频| 2018国产大陆天天弄谢| 女人精品久久久久毛片| 精品第一国产精品| 久久久精品94久久精品| 9191精品国产免费久久| 飞空精品影院首页| 久久亚洲国产成人精品v| 精品人妻偷拍中文字幕| 18禁裸乳无遮挡动漫免费视频| 好男人视频免费观看在线| 国产精品一二三区在线看| 欧美+日韩+精品| 国产永久视频网站| 日韩三级伦理在线观看| 18禁在线无遮挡免费观看视频| 另类精品久久| 久久久精品免费免费高清| 精品熟女少妇av免费看| 精品一区二区三区四区五区乱码 | 亚洲国产av新网站| 少妇猛男粗大的猛烈进出视频| 国产亚洲一区二区精品| 日韩在线高清观看一区二区三区| 欧美精品一区二区大全| 精品久久蜜臀av无| 久久国产亚洲av麻豆专区| 欧美3d第一页| 熟女av电影| 久久99精品国语久久久| 国产男女内射视频| 日韩免费高清中文字幕av| 国产欧美日韩一区二区三区在线| 不卡视频在线观看欧美| 在线观看免费日韩欧美大片| 精品一区在线观看国产| 一二三四在线观看免费中文在 | 国产精品 国内视频| 精品久久国产蜜桃| 国产精品一区www在线观看| 国产淫语在线视频| 久久久久久久久久久久大奶| 有码 亚洲区| 日韩 亚洲 欧美在线| 在线观看免费视频网站a站| 寂寞人妻少妇视频99o| 又粗又硬又长又爽又黄的视频| 欧美国产精品va在线观看不卡| 日本91视频免费播放| 免费观看在线日韩| 一级毛片 在线播放| 日韩欧美精品免费久久| 国产亚洲av片在线观看秒播厂| 夫妻性生交免费视频一级片| 精品一区二区免费观看| 精品国产国语对白av| 国产日韩欧美亚洲二区| 岛国毛片在线播放| 少妇被粗大猛烈的视频| 高清不卡的av网站| 欧美成人精品欧美一级黄| 久久久久国产网址| 最黄视频免费看| 久久青草综合色| 亚洲成人一二三区av| 国产精品不卡视频一区二区| 91午夜精品亚洲一区二区三区| 欧美少妇被猛烈插入视频| a 毛片基地| 国产午夜精品一二区理论片| 人成视频在线观看免费观看| 又粗又硬又长又爽又黄的视频| 亚洲av免费高清在线观看| 午夜福利视频在线观看免费| 午夜福利在线观看免费完整高清在| 精品午夜福利在线看| 国内精品宾馆在线| 久久久久久久亚洲中文字幕| 久久人妻熟女aⅴ| 欧美精品国产亚洲| 人妻少妇偷人精品九色| 少妇的逼水好多| 欧美精品一区二区免费开放| 人人妻人人添人人爽欧美一区卜| 性色av一级| 黄色一级大片看看| 亚洲欧美成人综合另类久久久| videossex国产| 久久久a久久爽久久v久久| 日本黄大片高清| 大片电影免费在线观看免费| 国产精品久久久久久精品电影小说| 两个人免费观看高清视频| 久久久亚洲精品成人影院| 国产日韩一区二区三区精品不卡| 成人二区视频| 国产男女超爽视频在线观看| 亚洲一级一片aⅴ在线观看| 久久久久网色| 国产不卡av网站在线观看| 国产 一区精品| 美女福利国产在线| 久久久a久久爽久久v久久| 永久免费av网站大全| 激情五月婷婷亚洲| 999精品在线视频| 亚洲人成77777在线视频| 又粗又硬又长又爽又黄的视频| 51国产日韩欧美| 91精品三级在线观看| 亚洲欧美成人综合另类久久久| 人妻人人澡人人爽人人| 男人添女人高潮全过程视频| 这个男人来自地球电影免费观看 | 国产欧美亚洲国产| 亚洲av男天堂| 成年人午夜在线观看视频| 宅男免费午夜| 最新中文字幕久久久久| 一本大道久久a久久精品| 国产av一区二区精品久久| 午夜av观看不卡| av卡一久久| 亚洲三级黄色毛片| 国产欧美日韩一区二区三区在线| 亚洲国产日韩一区二区| 菩萨蛮人人尽说江南好唐韦庄| 中文乱码字字幕精品一区二区三区| 久久亚洲国产成人精品v| 一区二区av电影网| 大香蕉97超碰在线| kizo精华| a级毛片在线看网站| 黑人高潮一二区| 亚洲人与动物交配视频| 丁香六月天网| 日本黄大片高清| 国产成人av激情在线播放| 精品亚洲成国产av| 天天躁夜夜躁狠狠久久av| 免费观看无遮挡的男女| 国产色婷婷99| 日本色播在线视频| 在现免费观看毛片| 国产一区二区三区综合在线观看 | 亚洲国产精品国产精品| 美女大奶头黄色视频| 一区二区三区四区激情视频| 国产一区有黄有色的免费视频| 国产黄频视频在线观看| 亚洲高清免费不卡视频| 国产av一区二区精品久久| av国产精品久久久久影院| 成年av动漫网址| 美女大奶头黄色视频| 国产1区2区3区精品| 成人综合一区亚洲| 国产精品一国产av| 国内精品宾馆在线| 久久久久久久久久成人| 日韩中文字幕视频在线看片| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 美女xxoo啪啪120秒动态图| 丝袜喷水一区| 制服诱惑二区| 国内精品宾馆在线| 十八禁网站网址无遮挡| 亚洲av福利一区| 欧美成人精品欧美一级黄| 纵有疾风起免费观看全集完整版| 观看美女的网站| 天天影视国产精品| 日本午夜av视频| a级毛片黄视频| 久久鲁丝午夜福利片| 2018国产大陆天天弄谢| 青春草国产在线视频| 插逼视频在线观看| 丰满乱子伦码专区| 久久久久久久国产电影| 一级片免费观看大全| 美女内射精品一级片tv| 国产又色又爽无遮挡免| 最近最新中文字幕免费大全7| 老熟女久久久| 日韩伦理黄色片| 亚洲av综合色区一区| 日本-黄色视频高清免费观看| 少妇猛男粗大的猛烈进出视频| 女性生殖器流出的白浆| 卡戴珊不雅视频在线播放| 午夜福利视频在线观看免费| 五月伊人婷婷丁香| 国产亚洲一区二区精品| 亚洲av国产av综合av卡| 日韩 亚洲 欧美在线| 最新的欧美精品一区二区| 韩国高清视频一区二区三区| 美女大奶头黄色视频| 在线 av 中文字幕| 一级,二级,三级黄色视频| 性色av一级| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 国产熟女午夜一区二区三区| 欧美亚洲日本最大视频资源| 最黄视频免费看| 久久久久久久国产电影| 亚洲国产av新网站| 黑人巨大精品欧美一区二区蜜桃 | 纯流量卡能插随身wifi吗| 男人操女人黄网站| 蜜桃国产av成人99| 少妇人妻 视频| 午夜精品国产一区二区电影| 免费播放大片免费观看视频在线观看| 青春草国产在线视频| 国产黄频视频在线观看| 啦啦啦啦在线视频资源| 超色免费av| 韩国av在线不卡| 你懂的网址亚洲精品在线观看| 欧美日韩国产mv在线观看视频| 少妇高潮的动态图| 高清不卡的av网站| 亚洲中文av在线| 观看美女的网站| 国产精品国产av在线观看| 99久久中文字幕三级久久日本| 久热久热在线精品观看| 高清黄色对白视频在线免费看| 大片免费播放器 马上看| 国产不卡av网站在线观看| 中国三级夫妇交换| 狠狠精品人妻久久久久久综合| 国产极品天堂在线| 久久国内精品自在自线图片| 欧美人与善性xxx| www.熟女人妻精品国产 | 五月伊人婷婷丁香| 欧美成人午夜精品| 边亲边吃奶的免费视频| 美国免费a级毛片| 国产精品.久久久| 精品第一国产精品| 久久这里有精品视频免费| 成人国语在线视频| 亚洲图色成人| 美女主播在线视频| 婷婷成人精品国产| 免费人妻精品一区二区三区视频| 人妻少妇偷人精品九色| 一级毛片黄色毛片免费观看视频| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 中文字幕最新亚洲高清| 成年动漫av网址| 欧美精品国产亚洲| 少妇精品久久久久久久| 亚洲国产av新网站| 亚洲精品国产av蜜桃| 制服诱惑二区| 麻豆精品久久久久久蜜桃| 国产精品 国内视频| a级毛片黄视频| 亚洲第一区二区三区不卡| 日韩免费高清中文字幕av| 波野结衣二区三区在线| 日日摸夜夜添夜夜爱| 成年人免费黄色播放视频| 视频中文字幕在线观看| 国产毛片在线视频| 国产高清国产精品国产三级| 久久久精品免费免费高清| 2018国产大陆天天弄谢| 国产免费又黄又爽又色| 国产一区亚洲一区在线观看| 乱码一卡2卡4卡精品| 色婷婷av一区二区三区视频| 国产欧美日韩一区二区三区在线| 制服诱惑二区| 免费av不卡在线播放| 啦啦啦中文免费视频观看日本| 伦精品一区二区三区| 精品第一国产精品| 人妻一区二区av| 久久热在线av| 尾随美女入室| 中国国产av一级| 18+在线观看网站| 黑人猛操日本美女一级片| 久热这里只有精品99| 男女边吃奶边做爰视频| 日本91视频免费播放| 美女国产视频在线观看| 在线亚洲精品国产二区图片欧美| 国产精品无大码| 国产一级毛片在线| av播播在线观看一区| 久久国内精品自在自线图片| 大香蕉久久网| 嫩草影院入口| 国产毛片在线视频| av又黄又爽大尺度在线免费看| 极品少妇高潮喷水抽搐| 亚洲综合色惰| 一个人免费看片子| 日韩中文字幕视频在线看片| 国产xxxxx性猛交| 免费久久久久久久精品成人欧美视频 | 国产精品女同一区二区软件| 欧美激情国产日韩精品一区| 啦啦啦中文免费视频观看日本| 国产一区二区激情短视频 | 黑丝袜美女国产一区| 婷婷色综合www|