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

    青藏高原地區(qū)草地綠度對物候變化的響應(yīng)

    2018-03-28 10:40馬雪紅程恒
    南方農(nóng)業(yè)·下旬 2018年12期
    關(guān)鍵詞:植被指數(shù)物候青藏高原

    馬雪紅 程恒

    摘 要 采用遙感植被指數(shù)數(shù)據(jù)評估植被綠度最大值、生長季開始時間、生長季結(jié)束時間三者與植被綠度年度變率之間的定量關(guān)系。結(jié)果表明:青藏高原植被指數(shù)的變化趨勢呈東部地區(qū)增加、西部地區(qū)減少,且越往東植被指數(shù)增加幅度越大,越往西減少幅度越大,整體呈上升趨勢;青藏高原植被指數(shù)最大值整體呈上升態(tài)勢,但在東南地區(qū)散布著許多植被綠度最大值呈增加趨勢的小區(qū)域;青藏高原地區(qū)植被生長季開始時間整體呈推后態(tài)勢,變化趨勢呈東南地區(qū)提前,西部、北部和西北地區(qū)推遲的空間格局;青藏高原地區(qū)植被生長季結(jié)束時間整體呈推遲態(tài)勢,其中東部、南部地區(qū)在推遲,西北地區(qū)在提前的變化趨勢;生長季開始時間與年NDVI指數(shù)呈負(fù)相關(guān),從東部向中部呈現(xiàn)出相關(guān)性逐漸增高的態(tài)勢;生長季結(jié)束時間SOS與年NDVI指數(shù)呈正相關(guān)。

    關(guān)鍵詞 青藏高原;物候;植被生理;植被指數(shù)

    中圖分類號:S812 文獻(xiàn)標(biāo)志碼:B DOI:10.19415/j.cnki.1673-890x.2018.36.060

    1 研究背景

    物候是指各類植物受自然環(huán)境(氣候、水文和土壤等)影響而呈現(xiàn)以年為周期的萌芽、展葉、開花和落葉等自然現(xiàn)象的現(xiàn)象[1]。決定植被活躍光合作用的時間和持續(xù)時間的植物物候?qū)W,在陸地生態(tài)系統(tǒng)年度碳吸收中具有重要作用。此外,植被物候的變化可能會影響各種生態(tài)系統(tǒng)特性,并通過改變陸地-大氣能量交換來反饋區(qū)域氣候。例如,春季綠化和植被秋季休眠的時間極大地影響植被生長、植被競爭和地表反照率,可能會進(jìn)一步影響局部溫度。因此,過去幾十年全球變化研究中,植被物候及其相關(guān)環(huán)境因子的空間和時間模式受到了越來越多的關(guān)注。遙感技術(shù)(如航空遙感、計算機(jī)遙感)和計算機(jī)技術(shù)(如圖形處理、大數(shù)據(jù)模擬、機(jī)器學(xué)習(xí))的應(yīng)用,使物候研究在數(shù)據(jù)采集、物候規(guī)律探索和物候綜合應(yīng)用方面取得了較大發(fā)展和重大突破。

    青藏高原的溫度、降水量和太陽輻射等氣候因子的變化,直接或間接地影響了青藏高原的植被生長,而植被生長又反過來影響碳循環(huán)、微氣候、土地粗糙度、地表反照率以及放牧及農(nóng)作物估產(chǎn)等經(jīng)濟(jì)活動[2]。新世紀(jì)以來,青藏高原綠度變化很大。張鐿鋰指出[3],青藏高原總體生長季植被覆蓋度呈增加趨勢,將青藏高原正在變綠的現(xiàn)象歸因于相關(guān)區(qū)域的氣候。此外,自然保護(hù)區(qū)的建設(shè)對增加中東部綠度和減緩西部綠度下降也有突出成效。卓嘎[4]將青藏高原植被覆蓋度變化進(jìn)行分區(qū)討論,認(rèn)為高原植被覆蓋既有增長又有退化趨勢。劉振元[5]則認(rèn)為,降水量、氣溫和太陽輻射是導(dǎo)致植被指數(shù)NDVI發(fā)生變化的原因。有研究表明[6],溫度正向影響NDVI,降水反向影響NDVI,得出溫度是三江源地區(qū)植被生長的主導(dǎo)因子。同時,氣候因子與植被指數(shù)相關(guān)性分析中,西部優(yōu)于中東部,海拔不同則響應(yīng)程度也不同??梢钥闯?,有關(guān)青藏高原的植被覆蓋度變化的研究很多,大部分研究認(rèn)可青藏高原植被覆蓋度整體上呈增加的趨勢。但是,對于植被覆蓋度變化的歸因分析中,多數(shù)研究將原因歸結(jié)于溫度、降水和輻射通量等氣候因素和退耕還林、畜牧等人為因素。所以,基于全球變化的現(xiàn)狀和生物圈物質(zhì)循環(huán)系統(tǒng)的重要性,研究青藏高原物候具有重要價值。

    因為所用數(shù)據(jù)和研究方法不同,不同研究中青藏高原物候變化結(jié)果大有不同。比如,張鐿鋰等人[3]的研究表明,自2000年以來,青藏高原總體生長季植被覆蓋度逐年增加,即青藏高原正在變綠。于海英指出[7],由于青藏高原冬季氣溫提高,一定程度上導(dǎo)致植物無法春化,從而推遲了返青期。田柳茜認(rèn)為[8],凍融過程、降水等因素導(dǎo)致了青藏高原植被覆蓋度的變化。學(xué)界對青藏高原綠度如何變化、為什么變化的成因分析很多,且多數(shù)研究都是氣候引起的植物物候變化和生理變化,而植物物候與生理的相對重要性對每年的NDVI變化影響卻不清楚。因此,大量研究皆將青藏高原綠度變化原因歸結(jié)于青藏高原地區(qū)的氣候變化和人為活動,卻幾乎沒有將其歸因于植被本身生理變化。基于此,研究青藏高原植被綠度變化與其生理之間的關(guān)系,填補(bǔ)此類研究的缺失。

    通過文獻(xiàn)調(diào)研,選擇GIMMS NDVI數(shù)據(jù)進(jìn)行青藏高原菜地綠度與物候間關(guān)系的研究。在數(shù)據(jù)降噪方面選擇奇異譜濾波方法,物候因素的提取選擇動態(tài)閾值法,即利用動態(tài)閾值法提取出生長季開始時間、生長季結(jié)束時間和植被綠度最大值三個物候因子。利用MATLAB軟件對三個因子進(jìn)行簡單相關(guān)分析,檢驗三大因子的相互獨立性,再對所提的三個因子與植被綠度變化做定性的偏相關(guān)分析,從而獲得了三大因子與植被綠度變化的相關(guān)性結(jié)果。

    2 數(shù)據(jù)與方法

    2.1 研究區(qū)介紹

    青藏高原位于中國的西北部,處于26°~39°N、73°~104°E。青藏高原西邊的起始點處于帕米爾高原,東部的終結(jié)點位于橫斷山,北部的邊界為昆侖山、阿爾金山,最南部抵達(dá)喜馬拉雅山[9]。青藏高原總面積2.5×106? km2,占全國總面積的26%,區(qū)域內(nèi)包含有西藏、青海、新疆、甘肅、四川和云南等地。

    受地形和氣候的影響,青藏高原生態(tài)系統(tǒng)獨特,植被類型復(fù)雜多樣,包括有山地森林、高寒草甸和高寒草原等,并且呈現(xiàn)出自東南向西北,植被類型從森林過渡到草原的水平地帶性分布規(guī)律。整體上看,高原區(qū)域草地占地面積最大,面積高達(dá)青藏高原區(qū)域總面積的1/2,牧草資源豐富。青藏高原作為中國的一個重要天然牧場,區(qū)域內(nèi)天然草地類型多樣,包括有高山草甸草地類、地灌叢草地類等八類。

    2.2 數(shù)據(jù)準(zhǔn)備

    2.2.1 NDVI數(shù)據(jù)

    本研究選擇歸一化植被指數(shù)進(jìn)行植被物候變化趨勢的分析。NDVI與植被的綠度和植被覆蓋都具有很好的相關(guān)性[10],被廣泛應(yīng)用于植被生長過程和趨勢變化研究中。歸一化植被指數(shù)數(shù)據(jù)有GIMMS NDVI、MODIS NDVI等遙感數(shù)據(jù)。在杜加強(qiáng)等人[11]的研究中,兩者捕獲植被物候數(shù)據(jù)精確度的能力基本一致。GIMMS NDVI具有時間橫跨范圍廣、數(shù)據(jù)研究區(qū)范圍廣、時間空間可比性強(qiáng)、具有較強(qiáng)的植被動態(tài)變化表現(xiàn)能力等特點。監(jiān)測植被的變化趨勢需要時間序列數(shù)據(jù)的時間周期長,而MODIS NDVI、SPOT VGT NDVI都是從2000年開始進(jìn)行數(shù)據(jù)獲取的,數(shù)據(jù)的時間跨度較短,沒有海量數(shù)據(jù)。歐盟支持生產(chǎn)的VEGETATION傳感器在1998年3月安裝在SPOT-4航天器上,升空后進(jìn)行數(shù)據(jù)獲取工作。因此,AVHRR NDVI數(shù)據(jù)集是連接植被生長周期歷史數(shù)據(jù)和當(dāng)代遙感衛(wèi)星數(shù)據(jù)的連接點。所以,這是在GIMMS NDVI數(shù)據(jù)存在不足的情況下依舊選取GIMMS NDVI的原因。

    2.2.2 土地覆蓋類型數(shù)據(jù)

    本研究選擇ESA CCI土地覆蓋產(chǎn)品。歐空局地面覆蓋數(shù)據(jù)涉及1992—2015年的24年,空間分辨率為300 m。這些地圖以聯(lián)合國土地覆蓋分類系統(tǒng)為基礎(chǔ),描述了37個原始LC類別中地球的陸地表面。這個獨特的長期土地覆蓋時間序列通過五種不同觀測系統(tǒng)的全球日常地表反射率實現(xiàn)。因為研究需要研究植被物候因子的時序變化之間的相關(guān)性與貢獻(xiàn),所以應(yīng)該保證在研究的時間區(qū)間內(nèi)該研究區(qū)的土地覆蓋類型一直是草地。所以,利用MATLAB將這24年間土地覆蓋類型沒有發(fā)生改變的地區(qū)標(biāo)為有效區(qū)域。

    2.3 平滑降噪NDVI數(shù)據(jù)

    雖然NDVI數(shù)據(jù)經(jīng)過了大氣、輻射、幾何校正等去掉數(shù)據(jù)噪聲因素的操作,產(chǎn)生了減弱NDVI噪聲因素的影響,但由于傳感器性的技術(shù)限制和使用限制,以及天氣等隨機(jī)因素的干擾,NDVI數(shù)據(jù)還存在一些殘余云及云霾等噪聲。這種隨機(jī)噪聲使得反映植被生長過程變化的時間序列曲線呈鋸齒狀分布,不利于進(jìn)行各種信息提取和趨勢分析。因此,對數(shù)據(jù)進(jìn)行平滑降噪處理必不可少,這里選用適合青藏高原地區(qū)的SSA濾波法。

    SSA(奇異譜分析濾波模型)是一種用于時間序列分析的增長分析方法,已實際應(yīng)用于多個科學(xué)領(lǐng)域。標(biāo)準(zhǔn)的SSA算法包括嵌入、SVD分解(本質(zhì)上為非平穩(wěn)時間序列信號處理的一種遞歸過程,具有改進(jìn)分析信號的時頻特性功能)、分組和對角平均(即重構(gòu))。綜上所述,SSA能夠?qū)r間序列分解,并對其進(jìn)行降噪平滑處理。

    2.4 研究區(qū)植被物候點的提取方法

    許多研究者采用的遙感植被指數(shù)數(shù)據(jù)分辨率各不相同,提出了多種以遙感數(shù)據(jù)提取物候信息的相關(guān)算法,如閾值法、主成分分析法、曲線擬合法、最累積頻率法、延遲滑動平均法和大變化斜率法等。當(dāng)代,獲取植被類型物候時期的方法多種多樣,全球范圍內(nèi)尚且沒有一個統(tǒng)一的方法,再加上青藏高原的植被類型屬于高寒植被類型,能采用的方法繁多。常清等人[12]選擇動態(tài)閾值法、最大變化斜率法及l(fā)ogistic曲線擬合法對青藏高原物候進(jìn)行物候期提取,發(fā)現(xiàn)動態(tài)閾值法比其他兩種方法速度更快,數(shù)據(jù)結(jié)果更加細(xì)致和精確?;诖?,研究選擇動態(tài)閾值法[13]:當(dāng)NDVI的數(shù)據(jù)值增長達(dá)到當(dāng)年NDVI數(shù)據(jù)值波動振幅的20%,即可定義為生長季的起始時間;而NDVI的數(shù)據(jù)值降低達(dá)到當(dāng)年NDVI數(shù)據(jù)值波動振幅的20%,即可定義生長季的結(jié)束時間;而將大于振幅95%的點全部挑選出來求平均值,定為NDVI最大值。因為該方法采用的動態(tài)閾值法的閾值是一種動態(tài)的給定像元和給定年份的植被指數(shù)值與當(dāng)年振幅的比值,可以有效避免由于因地域水熱狀況、植被類型和土壤類型等不同而造成的相互干擾。因此,動態(tài)閾值在時間和空間上的表現(xiàn)都優(yōu)于固定閾值。

    3 結(jié)果與分析

    3.1 青藏高原物候的綠度變化趨勢

    為了剖析青藏高原植被物候1982—2015年間的時空變化趨勢,將研究區(qū)內(nèi)每一個柵格中34年的物候數(shù)據(jù)都進(jìn)行線性回歸,再用每個柵格線性回歸系數(shù)表示物候參數(shù)變化幅度的大小。需要說明的是,年NDVI值和NDVI最大值的線性回歸系數(shù)的正負(fù)表示NDVI的增大和縮小,生長季開始/結(jié)束時間的線性回歸系數(shù)正負(fù)表示物候的推遲和提前。

    青藏高原在34年間的植被綠度變化率在空間上的分布狀況。青藏高原植被綠度呈現(xiàn)出年均不斷增長的趨勢,面積占青藏高原總面積的56.93%;植被綠度呈減少趨勢的面積占青藏高總面積的42.34%;青藏高原地區(qū)植被綠度在整個研究區(qū)上呈現(xiàn)出數(shù)據(jù)值不斷增長的趨勢。從空間分布分析可得,青藏高原東部地區(qū)、南部地區(qū)的植被綠度呈增加趨勢,變化率處于0~0.6。在地圖上可以明顯看出,東部三江源地區(qū)的增加趨勢最顯著,變化數(shù)值率處于0.2~0.6;青藏高原西部地區(qū)植被綠度呈減少趨勢,變化率在-0.2~0,且西部地區(qū)有部分散點的變化率達(dá)到了-1~-0.6。綜合上述研究數(shù)據(jù)和整個青藏高原研究區(qū)的數(shù)據(jù)分布看,變化趨勢呈東部增加、西部減少、向東增加的趨勢。

    3.2 生長季開始時間時空變化趨勢

    青藏高原1982—2015年植被生長季開始時間變化率的空間分布狀況。青藏高原生長季開始時間的數(shù)據(jù)值為負(fù)數(shù)的面積占青藏高原總面積的31.36%,植被物候終止時期的數(shù)據(jù)值為正數(shù)的面積占青藏高總面積的64.17%,因此青藏高原地區(qū)植被物候起始時間整體上在延后。在空間分布上,青藏高原東南部的三江源地區(qū)的植被生長季開始時間呈提前趨勢,變化率處于

    -0.5~0 d/a;青藏高原西部地區(qū)植被開始時間呈推遲趨勢,變化率在0~1 d/年,且西北地區(qū)植被開始時間推遲趨勢最明顯,變化率達(dá)到了0.5~-1.0 d/年。所以,青藏高原植被開始時間的變化趨勢整體上呈現(xiàn)的空間格局為:東南地區(qū)在提前,西部、北部和西北地區(qū)在推遲。

    3.3 生長季結(jié)束時間時空變化趨勢

    青藏高原1982—2015年植被生長季結(jié)束時間變化率的空間分布狀況。青藏高原生長季結(jié)束時間的數(shù)據(jù)值為負(fù)數(shù)的面積占青藏高原總面積的42.43%,生長季結(jié)束時間的數(shù)據(jù)值為正數(shù)的面積占青藏高總面積的56.12%,因此青藏高原地區(qū)生長季結(jié)束時間整體上呈推后趨勢。在空間分布上,青藏高原東部地區(qū)、南部地區(qū)的植被生長季結(jié)束時間的數(shù)據(jù)值為正數(shù),變化率處于0~7 d/年;西北角長條形區(qū)域也呈物候延后2 d左右的態(tài)勢,變化率處于2.6~7.1 d/年,且東部地區(qū)的物候終止期的數(shù)據(jù)值為正數(shù)的態(tài)勢相當(dāng)明顯;青藏高原西部地區(qū)植被數(shù)據(jù)值的變化率在-0.5~0 d/年,且西北部地區(qū)有部分散點的變化率達(dá)到了-2~-7.0 d/年。所以,青藏高原生長季結(jié)束時間的變化趨勢整體呈現(xiàn)的空間格局為:東部地區(qū)、南部地區(qū)在推遲,西北地區(qū)在提前。

    4 結(jié)論與展望

    本文立足于生物地球化學(xué)調(diào)控機(jī)制,利用遙感植被物候參數(shù),研究青藏高原地區(qū)草地綠度變化的生物成因,明晰高寒草地生態(tài)系統(tǒng)植被活性的驅(qū)動因子及其定量貢獻(xiàn)。結(jié)果表明,青藏高原植被指數(shù)東部地區(qū)增加、西部地區(qū)減少,且越往東植被指數(shù)增加幅度越大,越往西減少幅度越大,整體呈上升趨勢;青藏高原植被指數(shù)最大值整體呈上升態(tài)勢,但在東南地區(qū)散布著許多植被綠度最大值呈增加趨勢的小區(qū)域;青藏高原地區(qū)植被生長季開始時間的空間格局為東南地區(qū)在提前,西部、北部和西北地區(qū)在推遲,整體呈推后態(tài)勢;青藏高原地區(qū)植被生長季結(jié)束時間的變化趨勢為東部、南部地區(qū)在推遲,西北地區(qū)在提前,整體呈推遲態(tài)勢;生長季開始時間與年NDVI指數(shù)呈負(fù)相關(guān),從東部向中部呈現(xiàn)出相關(guān)性逐漸增高的態(tài)勢;生長季結(jié)束時間SOS與年NDVI指數(shù)呈正相關(guān),西部敏感性高于東部,高緯度地區(qū)的敏感性高于低緯度地區(qū)。

    未來研究中,可從以下方面加以改進(jìn):1)可以選擇多源遙感數(shù)據(jù),同時選用GIMMS NDVI和MODIS NDVI兩種數(shù)據(jù),對兩組結(jié)果進(jìn)行比對,得出的結(jié)論更貼合實際;2)物候提取方法可以多選擇幾種,然后進(jìn)行數(shù)據(jù)比對,選擇提取效果更好的數(shù)據(jù);3)可以考慮物候因子對氣候因素變化的響應(yīng),探究物候因子年際變化的原因。

    參考文獻(xiàn):

    [1] 竺可楨,宛敏渭.物候?qū)W[M].長沙:湖南教育出版社,1999.

    [2] Luo C, Xu G, Chao Z,et al.Effect ofwarming and grazing on litter mass loss and temperature sensitivity of litter and dung mass loss onthe Tibetan plateau[J].Global Change Biology,2010,16(5):1606-1617.

    [3] 張鐿鋰,李蘭暉,丁明軍,等.新世紀(jì)以來青藏高原綠度變化及動因[J].自然雜志,2017,39(3):173-178.

    [4] 卓嘎.青藏高原植被覆蓋時空變化及其對氣候因子的響應(yīng)[J].生態(tài)學(xué)報,2017,38(9):220-230.

    [5] 劉振元,張杰,陳立.青藏高原植被指數(shù)最新變化特征及其與氣候因子的關(guān)系[J].氣候與環(huán)境研究,2017,22(3):289-300.

    [6] 陳瓊,周強(qiáng),張海峰,等.三江源地區(qū)基于植被生長季的NDVI對氣候因子響應(yīng)的差異性研究[J].生態(tài)環(huán)境學(xué)報,2010,26(6):1284-1289.

    [7] 于海英.青藏高原植被生長季變化及對氣候的響應(yīng)研究[D].昆明:昆明植物所,2010.

    [8] 田柳茜.基于遙感技術(shù)的青藏高原返青期的研究[D].咸陽:西北農(nóng)林科技大學(xué),2015.

    [9] 張鐿鋰,李炳元,鄭度.論青藏高原范圍與面積[J].地理研究,2002(1):1-8.

    [10] Zhang G, Zhang Y, Dong J,et al.Green-up dates in the Tibetan Plateau have continuously advanced from 1982 to 2011[J].Proc Natl AcadSci USA,2013,110(11):4309-4314.

    [11] 杜加強(qiáng),舒儉民,王躍輝,等.青藏高原MODIS NDVI與GIMMS NDVI的對比[J].應(yīng)用生態(tài)學(xué)報,2014,25(2):533-544.

    [12] 常清,王思遠(yuǎn),孫云曉,等.青藏高原典型植被生長季遙感模型提取分析[J].地球信息科學(xué)學(xué)報,2014,16(5):815-823.

    [13] 張晗,任志遠(yuǎn).多種時序NDVI重建方法比較與應(yīng)用分析[J].中國農(nóng)業(yè)科學(xué),2014,47(15):2998-3008.

    (責(zé)任編輯:劉昀)

    猜你喜歡
    植被指數(shù)物候青藏高原
    耳釘
    氣候變化對民和杏發(fā)育期影響分析
    the HIGH WAY
    淺談植被指數(shù)的分類與應(yīng)用
    古代物候知識和物候歷
    山東省冬小麥2001—2010年生長季變化分析
    青藏高原筑“天路”
    欧美一区二区国产精品久久精品| 国产av一区在线观看免费| 韩国av一区二区三区四区| 亚洲第一区二区三区不卡| 国产一区二区在线观看日韩| 岛国在线免费视频观看| 色av中文字幕| 成人高潮视频无遮挡免费网站| 91字幕亚洲| 亚洲,欧美精品.| 日本免费a在线| 给我免费播放毛片高清在线观看| 日本在线视频免费播放| 欧美日韩综合久久久久久 | ponron亚洲| 亚洲av免费在线观看| 夜夜爽天天搞| 丰满人妻熟妇乱又伦精品不卡| 免费看光身美女| 欧美日本亚洲视频在线播放| 亚洲av日韩精品久久久久久密| 精品欧美国产一区二区三| 精品久久久久久久久久久久久| 国产伦一二天堂av在线观看| 久久国产乱子伦精品免费另类| 国产精品爽爽va在线观看网站| 俺也久久电影网| 日本黄色片子视频| 色尼玛亚洲综合影院| 搞女人的毛片| 99视频精品全部免费 在线| 无遮挡黄片免费观看| 露出奶头的视频| 亚洲人成电影免费在线| 久久久久久久久中文| 91久久精品电影网| 国产一区二区亚洲精品在线观看| 少妇人妻精品综合一区二区 | 久久99热这里只有精品18| 免费大片18禁| 免费看日本二区| 搡老熟女国产l中国老女人| 91午夜精品亚洲一区二区三区 | 一级作爱视频免费观看| 亚洲中文日韩欧美视频| 2021天堂中文幕一二区在线观| 九九热线精品视视频播放| 亚洲人与动物交配视频| 亚洲国产日韩欧美精品在线观看| 午夜福利18| 我要看日韩黄色一级片| 美女高潮的动态| 91狼人影院| 国产精品影院久久| 欧美日韩乱码在线| 人妻久久中文字幕网| 欧美精品啪啪一区二区三区| 九九久久精品国产亚洲av麻豆| 久久亚洲真实| 国产 一区 欧美 日韩| 小说图片视频综合网站| 国产免费av片在线观看野外av| 最近最新免费中文字幕在线| 又黄又爽又刺激的免费视频.| 午夜福利成人在线免费观看| 黄色一级大片看看| 一a级毛片在线观看| 久久6这里有精品| 极品教师在线视频| 亚洲专区国产一区二区| 99热只有精品国产| 国产麻豆成人av免费视频| 精品人妻一区二区三区麻豆 | 又爽又黄无遮挡网站| 成人鲁丝片一二三区免费| 亚洲男人的天堂狠狠| 亚洲第一电影网av| 小蜜桃在线观看免费完整版高清| 中文字幕av在线有码专区| 在线观看av片永久免费下载| 一级毛片久久久久久久久女| 成熟少妇高潮喷水视频| 国产精品国产高清国产av| 一区福利在线观看| 午夜免费激情av| 国产色爽女视频免费观看| 午夜福利在线观看免费完整高清在 | 国产成人av教育| 真人做人爱边吃奶动态| 亚洲精品日韩av片在线观看| 日韩欧美精品免费久久 | 欧美性猛交黑人性爽| 美女 人体艺术 gogo| 国产亚洲精品综合一区在线观看| av国产免费在线观看| 欧美激情在线99| 一卡2卡三卡四卡精品乱码亚洲| 成人国产综合亚洲| 精品乱码久久久久久99久播| 亚洲av免费高清在线观看| 国产探花在线观看一区二区| 久久久久亚洲av毛片大全| 俄罗斯特黄特色一大片| 久久伊人香网站| 精品不卡国产一区二区三区| 色哟哟·www| 亚洲人与动物交配视频| 两性午夜刺激爽爽歪歪视频在线观看| 精品国产三级普通话版| 午夜福利视频1000在线观看| 99热精品在线国产| 少妇人妻一区二区三区视频| 成人无遮挡网站| 日本一本二区三区精品| 精品人妻偷拍中文字幕| 亚洲av一区综合| 国产v大片淫在线免费观看| 国产毛片a区久久久久| 一进一出好大好爽视频| 国产欧美日韩一区二区精品| 亚洲av中文字字幕乱码综合| 亚洲精品在线观看二区| 亚洲国产精品成人综合色| 高清在线国产一区| 成人毛片a级毛片在线播放| a级毛片免费高清观看在线播放| 亚洲人成伊人成综合网2020| 又黄又爽又免费观看的视频| 国产高清视频在线观看网站| 性欧美人与动物交配| 在线十欧美十亚洲十日本专区| 一区福利在线观看| 日韩有码中文字幕| 一本综合久久免费| 日韩欧美一区二区三区在线观看| 能在线免费观看的黄片| 欧美精品啪啪一区二区三区| 一二三四社区在线视频社区8| 国产美女午夜福利| 黄色配什么色好看| 欧美精品国产亚洲| 女人被狂操c到高潮| 黄色日韩在线| 精品久久国产蜜桃| 欧美极品一区二区三区四区| 我的老师免费观看完整版| 久久久久久大精品| 亚洲人成网站在线播| h日本视频在线播放| 久久久久久久久久黄片| 国产日本99.免费观看| 少妇熟女aⅴ在线视频| 久久久成人免费电影| 麻豆成人午夜福利视频| 少妇高潮的动态图| 一二三四社区在线视频社区8| 97超级碰碰碰精品色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av免费在线观看| x7x7x7水蜜桃| 在线观看免费视频日本深夜| 国产精品自产拍在线观看55亚洲| 又爽又黄a免费视频| 免费在线观看影片大全网站| 免费人成视频x8x8入口观看| 成人美女网站在线观看视频| www日本黄色视频网| 国产精华一区二区三区| 国产成+人综合+亚洲专区| 草草在线视频免费看| 亚洲欧美清纯卡通| 一本综合久久免费| 哪里可以看免费的av片| 成熟少妇高潮喷水视频| 亚洲 国产 在线| 天美传媒精品一区二区| 成人一区二区视频在线观看| 一a级毛片在线观看| 在现免费观看毛片| 特级一级黄色大片| 欧美精品国产亚洲| 亚洲男人的天堂狠狠| 12—13女人毛片做爰片一| 美女xxoo啪啪120秒动态图 | 久久精品综合一区二区三区| 尤物成人国产欧美一区二区三区| 亚洲精品在线观看二区| 嫩草影院入口| 亚洲午夜理论影院| 99热只有精品国产| 中文字幕免费在线视频6| 日韩中文字幕欧美一区二区| 亚洲av熟女| 亚洲在线自拍视频| av欧美777| 午夜免费激情av| 国产熟女xx| 91麻豆精品激情在线观看国产| 热99re8久久精品国产| 又黄又爽又免费观看的视频| 久久久精品欧美日韩精品| 亚洲avbb在线观看| 精品人妻1区二区| 亚洲五月天丁香| 国产欧美日韩精品一区二区| 国产精品1区2区在线观看.| АⅤ资源中文在线天堂| 国产精品永久免费网站| 亚洲国产精品sss在线观看| 国产亚洲av嫩草精品影院| 在线播放无遮挡| 亚洲综合色惰| 女人被狂操c到高潮| 精品一区二区三区视频在线观看免费| 亚洲av熟女| 99热精品在线国产| 国产精品免费一区二区三区在线| 亚洲精品456在线播放app | 亚洲av成人av| 色尼玛亚洲综合影院| 欧美+亚洲+日韩+国产| 亚洲精品影视一区二区三区av| 可以在线观看毛片的网站| 精品国产三级普通话版| 欧美中文日本在线观看视频| ponron亚洲| 亚洲五月天丁香| 亚洲午夜理论影院| 亚洲av美国av| 久久草成人影院| 不卡一级毛片| 岛国在线免费视频观看| 中文字幕高清在线视频| 欧美黄色片欧美黄色片| 国产69精品久久久久777片| 精品一区二区三区视频在线观看免费| 亚洲精华国产精华精| 少妇的逼水好多| 久久婷婷人人爽人人干人人爱| 精品欧美国产一区二区三| 欧美日韩瑟瑟在线播放| 国产精品三级大全| 99久久成人亚洲精品观看| 亚洲精品在线观看二区| 简卡轻食公司| 三级男女做爰猛烈吃奶摸视频| or卡值多少钱| 免费观看精品视频网站| 亚洲激情在线av| 桃色一区二区三区在线观看| 午夜福利高清视频| 欧美乱妇无乱码| 欧美激情久久久久久爽电影| 免费看美女性在线毛片视频| 真人做人爱边吃奶动态| 亚洲熟妇中文字幕五十中出| 亚洲美女搞黄在线观看 | 日本免费一区二区三区高清不卡| 97热精品久久久久久| 好看av亚洲va欧美ⅴa在| 国产高清视频在线观看网站| 国产又黄又爽又无遮挡在线| 国产蜜桃级精品一区二区三区| 国产成人av教育| 一区二区三区激情视频| 18美女黄网站色大片免费观看| 91狼人影院| 少妇的逼好多水| av黄色大香蕉| 亚洲精品日韩av片在线观看| 夜夜夜夜夜久久久久| 人妻夜夜爽99麻豆av| 色精品久久人妻99蜜桃| 黄色配什么色好看| 51午夜福利影视在线观看| 看免费av毛片| 我的女老师完整版在线观看| 亚洲人成网站在线播| 最好的美女福利视频网| 亚洲av免费高清在线观看| h日本视频在线播放| 亚洲av一区综合| 麻豆一二三区av精品| 精品乱码久久久久久99久播| 日本撒尿小便嘘嘘汇集6| 好看av亚洲va欧美ⅴa在| 亚洲自拍偷在线| 91在线精品国自产拍蜜月| 成熟少妇高潮喷水视频| 欧美日韩瑟瑟在线播放| 中文字幕人成人乱码亚洲影| 欧美激情在线99| www.熟女人妻精品国产| 亚洲国产高清在线一区二区三| 99在线人妻在线中文字幕| 桃色一区二区三区在线观看| h日本视频在线播放| 亚洲av成人av| 欧美黄色淫秽网站| 精品不卡国产一区二区三区| 国产黄a三级三级三级人| 美女 人体艺术 gogo| 757午夜福利合集在线观看| 色吧在线观看| 少妇人妻精品综合一区二区 | 亚洲精品一卡2卡三卡4卡5卡| 国产精品精品国产色婷婷| 国产欧美日韩精品一区二区| 亚洲精品粉嫩美女一区| 国产色婷婷99| 一二三四社区在线视频社区8| 午夜久久久久精精品| 久久亚洲精品不卡| 久久伊人香网站| 亚洲国产精品999在线| 51国产日韩欧美| 99视频精品全部免费 在线| 丝袜美腿在线中文| 国产乱人伦免费视频| 一区二区三区四区激情视频 | www.www免费av| 看十八女毛片水多多多| 毛片一级片免费看久久久久 | 动漫黄色视频在线观看| 久久性视频一级片| 亚洲 欧美 日韩 在线 免费| 亚洲成av人片免费观看| 99视频精品全部免费 在线| 毛片一级片免费看久久久久 | 麻豆久久精品国产亚洲av| 一本久久中文字幕| 成年女人看的毛片在线观看| 日韩高清综合在线| 欧美区成人在线视频| 亚洲黑人精品在线| 欧美最新免费一区二区三区 | 国产亚洲精品av在线| 久久精品综合一区二区三区| 赤兔流量卡办理| 亚洲国产欧美人成| 色综合婷婷激情| 成人特级av手机在线观看| 可以在线观看毛片的网站| 国产欧美日韩精品亚洲av| 99热这里只有是精品50| 亚洲av免费在线观看| 国产大屁股一区二区在线视频| 一区福利在线观看| www.熟女人妻精品国产| 国产精品免费一区二区三区在线| 国产成人a区在线观看| 亚洲综合色惰| 国产高清三级在线| 免费高清视频大片| 99精品在免费线老司机午夜| 色哟哟·www| 国产欧美日韩精品一区二区| 久久精品综合一区二区三区| 国产精品不卡视频一区二区 | 波野结衣二区三区在线| 欧美成人一区二区免费高清观看| 99视频精品全部免费 在线| 久久久久国内视频| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 在线a可以看的网站| 夜夜爽天天搞| 中文亚洲av片在线观看爽| 午夜老司机福利剧场| 精品久久久久久,| 亚洲欧美激情综合另类| av天堂在线播放| 在线免费观看不下载黄p国产 | 欧美性猛交╳xxx乱大交人| 国产亚洲精品综合一区在线观看| 亚洲av一区综合| 内射极品少妇av片p| 深爱激情五月婷婷| 成人精品一区二区免费| 99国产综合亚洲精品| 好看av亚洲va欧美ⅴa在| 我要搜黄色片| 99久久九九国产精品国产免费| 丝袜美腿在线中文| 欧美区成人在线视频| 桃红色精品国产亚洲av| 99久久久亚洲精品蜜臀av| 久99久视频精品免费| 91久久精品电影网| xxxwww97欧美| 午夜福利欧美成人| 精品人妻一区二区三区麻豆 | 我的老师免费观看完整版| 久久久国产成人免费| 国内少妇人妻偷人精品xxx网站| 久久午夜亚洲精品久久| 国内少妇人妻偷人精品xxx网站| 51国产日韩欧美| 中文字幕av成人在线电影| 精品无人区乱码1区二区| 神马国产精品三级电影在线观看| 欧美成狂野欧美在线观看| 人妻制服诱惑在线中文字幕| 啪啪无遮挡十八禁网站| 久久久久国产精品人妻aⅴ院| 欧美日韩中文字幕国产精品一区二区三区| 国内精品久久久久精免费| 哪里可以看免费的av片| 在线播放无遮挡| 性色av乱码一区二区三区2| 日韩精品中文字幕看吧| 精品久久久久久久末码| 色哟哟哟哟哟哟| 日韩 亚洲 欧美在线| 久久精品人妻少妇| 老司机午夜福利在线观看视频| av国产免费在线观看| 欧美高清性xxxxhd video| 一进一出好大好爽视频| 脱女人内裤的视频| av在线天堂中文字幕| 日韩国内少妇激情av| 免费搜索国产男女视频| 日韩中文字幕欧美一区二区| 国产成人aa在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲欧美精品综合久久99| 亚洲美女视频黄频| 欧美一级a爱片免费观看看| 麻豆国产97在线/欧美| 在线播放无遮挡| 日日干狠狠操夜夜爽| 看十八女毛片水多多多| 欧美成狂野欧美在线观看| 日韩人妻高清精品专区| 国产亚洲欧美98| 久久精品国产自在天天线| 久久久久精品国产欧美久久久| 亚洲,欧美精品.| 美女被艹到高潮喷水动态| 最近视频中文字幕2019在线8| 91狼人影院| 免费人成视频x8x8入口观看| 国产午夜福利久久久久久| 国产人妻一区二区三区在| 精品人妻1区二区| 亚洲av不卡在线观看| 久久久久久久久大av| 精品国内亚洲2022精品成人| 免费看美女性在线毛片视频| av福利片在线观看| 精品熟女少妇八av免费久了| 最好的美女福利视频网| 别揉我奶头 嗯啊视频| 精品乱码久久久久久99久播| 中文字幕高清在线视频| 欧美激情在线99| 亚洲中文字幕日韩| 一级黄片播放器| eeuss影院久久| 一边摸一边抽搐一进一小说| 亚洲国产高清在线一区二区三| 精品人妻熟女av久视频| 18禁黄网站禁片免费观看直播| 蜜桃久久精品国产亚洲av| 深夜精品福利| 欧美高清成人免费视频www| 一区二区三区四区激情视频 | 99久久精品国产亚洲精品| 色噜噜av男人的天堂激情| 老女人水多毛片| 国产亚洲欧美在线一区二区| 在线十欧美十亚洲十日本专区| 午夜老司机福利剧场| 亚洲专区中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲18禁久久av| 亚洲国产日韩欧美精品在线观看| 嫩草影院入口| 日韩人妻高清精品专区| 欧美日韩亚洲国产一区二区在线观看| 五月玫瑰六月丁香| 亚洲,欧美,日韩| 一边摸一边抽搐一进一小说| 91午夜精品亚洲一区二区三区 | 免费无遮挡裸体视频| 在现免费观看毛片| 亚洲欧美日韩无卡精品| 尤物成人国产欧美一区二区三区| 丰满乱子伦码专区| 精品福利观看| 国产成人aa在线观看| 欧美在线黄色| 亚洲三级黄色毛片| 超碰av人人做人人爽久久| 在线a可以看的网站| 日韩 亚洲 欧美在线| 欧美bdsm另类| 成人一区二区视频在线观看| 91在线精品国自产拍蜜月| 动漫黄色视频在线观看| 国产69精品久久久久777片| 国产欧美日韩一区二区三| 波多野结衣巨乳人妻| 91久久精品国产一区二区成人| 色哟哟哟哟哟哟| 又爽又黄无遮挡网站| 丰满人妻一区二区三区视频av| 欧美高清成人免费视频www| 老鸭窝网址在线观看| 久久久久久九九精品二区国产| 亚洲精品456在线播放app | 麻豆国产av国片精品| 精品久久久久久,| 九色国产91popny在线| 极品教师在线免费播放| 99视频精品全部免费 在线| 亚洲人成伊人成综合网2020| 日韩欧美一区二区三区在线观看| 在线国产一区二区在线| 国模一区二区三区四区视频| 亚洲av成人不卡在线观看播放网| 国内久久婷婷六月综合欲色啪| 色综合站精品国产| 国产白丝娇喘喷水9色精品| 亚洲人成电影免费在线| 国产主播在线观看一区二区| 天堂网av新在线| 精品免费久久久久久久清纯| 色哟哟哟哟哟哟| 婷婷丁香在线五月| 欧美精品国产亚洲| 久久久久国产精品人妻aⅴ院| 国产淫片久久久久久久久 | 日本撒尿小便嘘嘘汇集6| 免费av毛片视频| 赤兔流量卡办理| 久久天躁狠狠躁夜夜2o2o| 在线观看av片永久免费下载| 日本在线视频免费播放| 亚洲狠狠婷婷综合久久图片| 三级男女做爰猛烈吃奶摸视频| 一级黄片播放器| 久久久久久久久大av| 韩国av一区二区三区四区| 日韩欧美精品免费久久 | 级片在线观看| 我要看日韩黄色一级片| 成人特级黄色片久久久久久久| 麻豆av噜噜一区二区三区| 在线免费观看的www视频| 特级一级黄色大片| 精品一区二区三区视频在线| av欧美777| 中文字幕av在线有码专区| 国模一区二区三区四区视频| 2021天堂中文幕一二区在线观| 国产高清三级在线| 91在线精品国自产拍蜜月| 成人毛片a级毛片在线播放| 内地一区二区视频在线| 一个人免费在线观看的高清视频| 九九在线视频观看精品| 国产成人欧美在线观看| 亚洲18禁久久av| 国产毛片a区久久久久| 国产麻豆成人av免费视频| www.色视频.com| 9191精品国产免费久久| 给我免费播放毛片高清在线观看| 五月伊人婷婷丁香| 国产伦精品一区二区三区四那| 精品久久久久久久末码| 久久人人精品亚洲av| 最近最新中文字幕大全电影3| 乱人视频在线观看| 日本黄大片高清| 内地一区二区视频在线| 我要看日韩黄色一级片| 精品一区二区三区视频在线观看免费| 老司机福利观看| 亚洲无线在线观看| 国产老妇女一区| 一级av片app| 成年女人永久免费观看视频| 一级黄色大片毛片| 免费看日本二区| 国产一区二区亚洲精品在线观看| 国产高潮美女av| 国产单亲对白刺激| 麻豆国产97在线/欧美| 麻豆成人av在线观看| 久久精品国产99精品国产亚洲性色| 波多野结衣高清作品| 成年版毛片免费区| 非洲黑人性xxxx精品又粗又长| 欧美午夜高清在线| 亚洲一区二区三区不卡视频| 欧美色视频一区免费| 久久婷婷人人爽人人干人人爱| 欧美日韩国产亚洲二区| 一本一本综合久久| 男女做爰动态图高潮gif福利片| 国产亚洲精品综合一区在线观看| 欧美乱妇无乱码| 国产高清激情床上av| netflix在线观看网站| 精品欧美国产一区二区三| 婷婷六月久久综合丁香| 国产精品亚洲美女久久久| 国产人妻一区二区三区在| 国产精品一及| 国产精品国产高清国产av| 99久久精品国产亚洲精品| 99在线人妻在线中文字幕|