• <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系列| 亚洲国产精品久久男人天堂| 99久国产av精品国产电影| 国产黄a三级三级三级人| 久久国产乱子免费精品| 成年女人毛片免费观看观看9| 久久综合国产亚洲精品| 免费高清视频大片| 精品乱码久久久久久99久播| 嫩草影视91久久| av福利片在线观看| 免费无遮挡裸体视频| 国产伦精品一区二区三区四那| 午夜福利在线观看吧| 天天躁日日操中文字幕| 日韩av在线大香蕉| 搡老妇女老女人老熟妇| 久久久精品欧美日韩精品| 校园人妻丝袜中文字幕| 深夜精品福利| 波多野结衣高清作品| 色吧在线观看| 午夜激情欧美在线| 国产伦一二天堂av在线观看| 免费搜索国产男女视频| av在线天堂中文字幕| 九九热线精品视视频播放| 精品国产三级普通话版| 亚洲美女黄片视频| 搡老岳熟女国产| 午夜福利18| 久久精品国产自在天天线| 国产一区二区三区av在线 | 一区二区三区免费毛片| 99在线人妻在线中文字幕| 天堂√8在线中文| 97超视频在线观看视频| 18禁黄网站禁片免费观看直播| 大型黄色视频在线免费观看| 久久国产乱子免费精品| a级毛片a级免费在线| 亚洲精品影视一区二区三区av| 精品久久国产蜜桃| 深夜a级毛片| 丰满的人妻完整版| 久久亚洲精品不卡| 国产精品久久久久久亚洲av鲁大| 日本色播在线视频| 波多野结衣高清作品| 男人的好看免费观看在线视频| 中出人妻视频一区二区| 中文资源天堂在线| 少妇丰满av| 综合色丁香网| 欧美成人免费av一区二区三区| 不卡视频在线观看欧美| 亚洲成人中文字幕在线播放| 性插视频无遮挡在线免费观看| 亚洲av不卡在线观看| 特级一级黄色大片| 人人妻,人人澡人人爽秒播| 狂野欧美激情性xxxx在线观看| 精品人妻熟女av久视频| 成人欧美大片| .国产精品久久| 亚洲美女黄片视频| 两个人的视频大全免费| 国产乱人偷精品视频| 真实男女啪啪啪动态图| 欧美日本视频| 三级国产精品欧美在线观看| 国产一区二区激情短视频| 毛片一级片免费看久久久久| 精品久久久久久成人av| 三级经典国产精品| 99热只有精品国产| 狂野欧美白嫩少妇大欣赏| 精品久久久久久久人妻蜜臀av| av天堂中文字幕网| 久久人人爽人人爽人人片va| 久久久色成人| 成人欧美大片| 国产 一区精品| 少妇裸体淫交视频免费看高清| 免费看av在线观看网站| 人人妻人人澡人人爽人人夜夜 | 国产成人福利小说| 日韩成人av中文字幕在线观看 | 日韩欧美在线乱码| 春色校园在线视频观看| 最近的中文字幕免费完整| 久久草成人影院| 亚洲人成网站高清观看| 中文资源天堂在线| 国产黄片美女视频| 悠悠久久av| 日产精品乱码卡一卡2卡三| 九九久久精品国产亚洲av麻豆| 卡戴珊不雅视频在线播放| 国产麻豆成人av免费视频| 中文字幕av在线有码专区| 国产午夜福利久久久久久| 日韩制服骚丝袜av| 18禁在线无遮挡免费观看视频 | av国产免费在线观看| 我要搜黄色片| 午夜久久久久精精品| av福利片在线观看| 亚洲av免费在线观看| 一边摸一边抽搐一进一小说| 亚洲丝袜综合中文字幕| 日本五十路高清| 成年女人永久免费观看视频| 国产视频内射| 一个人看视频在线观看www免费| 亚洲成人精品中文字幕电影| 91在线精品国自产拍蜜月| 成年免费大片在线观看| 国内精品一区二区在线观看| 日韩欧美一区二区三区在线观看| 18禁裸乳无遮挡免费网站照片| 直男gayav资源| 精品不卡国产一区二区三区| 久久精品影院6| 亚洲人与动物交配视频| 又黄又爽又免费观看的视频| 在线免费观看的www视频| 18禁裸乳无遮挡免费网站照片| 亚洲av成人精品一区久久| 亚洲一区二区三区色噜噜| 亚洲熟妇中文字幕五十中出| 干丝袜人妻中文字幕| 成人永久免费在线观看视频| avwww免费| 夜夜爽天天搞| 精品久久久噜噜| 欧美区成人在线视频| 在线免费观看的www视频| 成人性生交大片免费视频hd| 久久精品国产亚洲av涩爱 | av在线播放精品| 女的被弄到高潮叫床怎么办| 日本黄色视频三级网站网址| 99九九线精品视频在线观看视频| 女人被狂操c到高潮| a级一级毛片免费在线观看| 我的老师免费观看完整版| 国国产精品蜜臀av免费| 成熟少妇高潮喷水视频| 永久网站在线| 天美传媒精品一区二区| 内地一区二区视频在线| 色哟哟哟哟哟哟| 亚洲最大成人av| 久久精品夜色国产| 波野结衣二区三区在线| 色综合站精品国产| 欧美国产日韩亚洲一区| 久久欧美精品欧美久久欧美| 天天躁日日操中文字幕| 久久久久久久久中文| 日本色播在线视频| 听说在线观看完整版免费高清| 亚洲五月天丁香| 精品一区二区免费观看| 变态另类丝袜制服| 丰满的人妻完整版| 十八禁网站免费在线| 国产精品免费一区二区三区在线| 午夜激情福利司机影院| 黄色日韩在线| 久久精品国产清高在天天线| 国产精品久久久久久av不卡| 狂野欧美白嫩少妇大欣赏| 国产一区二区在线av高清观看| 国产精品女同一区二区软件| 丝袜美腿在线中文| 亚洲精品影视一区二区三区av| 成人鲁丝片一二三区免费| 成人欧美大片| 夜夜看夜夜爽夜夜摸| 男人舔女人下体高潮全视频| 免费大片18禁| 亚洲成a人片在线一区二区| 国内精品久久久久精免费| 伦精品一区二区三区| 大型黄色视频在线免费观看| 日韩欧美三级三区| 九色成人免费人妻av| 深夜a级毛片| 我要看日韩黄色一级片| 国产中年淑女户外野战色| 国产乱人偷精品视频| 国产精品伦人一区二区| 97热精品久久久久久| 给我免费播放毛片高清在线观看| 老司机福利观看| 高清午夜精品一区二区三区 | 女人十人毛片免费观看3o分钟| 亚洲在线观看片| 精品人妻视频免费看| 一本久久中文字幕| 国产精品亚洲美女久久久| 丰满人妻一区二区三区视频av| 麻豆av噜噜一区二区三区| 看十八女毛片水多多多| 色吧在线观看| 美女xxoo啪啪120秒动态图| 久久精品影院6| 亚洲四区av| 中出人妻视频一区二区| a级一级毛片免费在线观看| 九九在线视频观看精品| 午夜精品一区二区三区免费看| 亚洲中文日韩欧美视频| 亚洲欧美精品自产自拍| 99久久中文字幕三级久久日本| 久久天躁狠狠躁夜夜2o2o| 欧美另类亚洲清纯唯美| 国产三级中文精品| 国产高清激情床上av| 99久久精品国产国产毛片| 亚洲av中文字字幕乱码综合| 真人做人爱边吃奶动态| 一级a爱片免费观看的视频| 1024手机看黄色片| 少妇的逼水好多| 99热网站在线观看| 国产亚洲欧美98| 欧美激情久久久久久爽电影| 色哟哟·www| 禁无遮挡网站| 国产精品电影一区二区三区| 久久中文看片网| 给我免费播放毛片高清在线观看| 久久精品国产鲁丝片午夜精品| 成人亚洲精品av一区二区| 日本黄色片子视频| 夜夜爽天天搞| 亚洲最大成人手机在线| 国产精品久久视频播放| 男人和女人高潮做爰伦理| 99riav亚洲国产免费| 老熟妇仑乱视频hdxx| 亚洲av中文av极速乱| 日韩精品有码人妻一区| 亚洲精品粉嫩美女一区| 日日干狠狠操夜夜爽| 最近视频中文字幕2019在线8| 亚洲精品久久国产高清桃花| 毛片女人毛片| 综合色丁香网| 亚洲自偷自拍三级| 久久久国产成人免费| 免费观看的影片在线观看| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 神马国产精品三级电影在线观看| av免费在线看不卡| 神马国产精品三级电影在线观看| 国产亚洲精品久久久久久毛片| 成人高潮视频无遮挡免费网站| 伊人久久精品亚洲午夜| 你懂的网址亚洲精品在线观看 | ponron亚洲| 日韩欧美三级三区| 秋霞在线观看毛片| 美女高潮的动态| 最新中文字幕久久久久| 最近的中文字幕免费完整| 男女之事视频高清在线观看| 能在线免费观看的黄片| 日韩成人av中文字幕在线观看 | 国产高清三级在线| av.在线天堂| 99久久无色码亚洲精品果冻| 人妻丰满熟妇av一区二区三区| 99久久精品热视频| 在线免费观看不下载黄p国产| 亚洲第一区二区三区不卡| 亚洲18禁久久av| 中国美女看黄片| 成人精品一区二区免费| 久久婷婷人人爽人人干人人爱| 桃色一区二区三区在线观看| 成人特级黄色片久久久久久久| 两个人视频免费观看高清| 中文在线观看免费www的网站| 美女大奶头视频| 一个人观看的视频www高清免费观看| 国产乱人偷精品视频| 99久久久亚洲精品蜜臀av| 一区福利在线观看| 久久久精品欧美日韩精品| 中出人妻视频一区二区| 日韩欧美在线乱码| 国产成人freesex在线 | 老熟妇仑乱视频hdxx| 久久久久久九九精品二区国产| 村上凉子中文字幕在线| 天堂动漫精品| 精品一区二区三区视频在线| 欧美最新免费一区二区三区| 日韩三级伦理在线观看| 国产一区二区亚洲精品在线观看| 少妇人妻一区二区三区视频| 免费看a级黄色片| 天堂影院成人在线观看| 亚洲aⅴ乱码一区二区在线播放| 51国产日韩欧美| 少妇高潮的动态图| 99久久中文字幕三级久久日本| 精品人妻熟女av久视频| 国内精品美女久久久久久| 国产精品女同一区二区软件| 日韩欧美在线乱码| 国产成年人精品一区二区| 精品久久久久久久久久久久久| 亚洲第一区二区三区不卡| 久久精品影院6| 亚洲精品456在线播放app| 免费av不卡在线播放| 淫秽高清视频在线观看| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 性欧美人与动物交配| 日韩人妻高清精品专区| 18禁裸乳无遮挡免费网站照片| 网址你懂的国产日韩在线| 国产男人的电影天堂91| 淫妇啪啪啪对白视频| 99国产极品粉嫩在线观看| 最近2019中文字幕mv第一页| 天天一区二区日本电影三级| 精品一区二区免费观看| 一卡2卡三卡四卡精品乱码亚洲| 久久久久免费精品人妻一区二区| 丰满乱子伦码专区| 欧美xxxx性猛交bbbb| 91久久精品电影网| 两个人视频免费观看高清| 日韩国内少妇激情av| 色综合亚洲欧美另类图片| 亚洲av第一区精品v没综合| 看黄色毛片网站| 亚洲精品亚洲一区二区| 熟女人妻精品中文字幕| 国产高清不卡午夜福利| 成人无遮挡网站| 麻豆一二三区av精品| 我的女老师完整版在线观看| 中文字幕免费在线视频6| 亚洲国产精品sss在线观看| 最近手机中文字幕大全| av.在线天堂| 偷拍熟女少妇极品色| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月| 国产精品久久电影中文字幕| 国产高清视频在线播放一区| 狂野欧美白嫩少妇大欣赏| 欧美区成人在线视频| 中国美女看黄片| 免费看光身美女| 亚洲七黄色美女视频| 精品一区二区三区av网在线观看| 免费高清视频大片| 人妻久久中文字幕网| 少妇丰满av| 久久人人爽人人片av| 亚洲精品在线观看二区| 久久午夜福利片| 伦理电影大哥的女人| 国国产精品蜜臀av免费| 91在线观看av| 日韩精品有码人妻一区| 天堂影院成人在线观看| 草草在线视频免费看| 国内久久婷婷六月综合欲色啪| 亚洲不卡免费看| 久久久久久九九精品二区国产| 美女内射精品一级片tv| 99视频精品全部免费 在线| 久久久精品欧美日韩精品| 搡老熟女国产l中国老女人| 久久精品影院6| 日韩一区二区视频免费看| 成人国产麻豆网| 国内久久婷婷六月综合欲色啪| 国产精品综合久久久久久久免费| 99热6这里只有精品| 老熟妇乱子伦视频在线观看| 亚洲高清免费不卡视频| 免费人成在线观看视频色| 变态另类丝袜制服| 国产中年淑女户外野战色| 日韩av不卡免费在线播放| 99久久精品国产国产毛片| 免费观看在线日韩| 亚洲欧美清纯卡通| 日本成人三级电影网站| 欧美日韩乱码在线| 日韩精品中文字幕看吧| 国产一区二区三区av在线 | av卡一久久| 韩国av在线不卡| 中文字幕久久专区| 国产精品久久久久久亚洲av鲁大| 一区二区三区免费毛片| 亚洲国产精品久久男人天堂| 亚洲精品一卡2卡三卡4卡5卡| 成年女人看的毛片在线观看| 日韩强制内射视频| 一本精品99久久精品77| 久久精品国产自在天天线| 欧美国产日韩亚洲一区| 人妻夜夜爽99麻豆av| 国产精品一区www在线观看| 欧美日韩综合久久久久久| 性色avwww在线观看| 国产亚洲精品av在线| 69人妻影院| 真人做人爱边吃奶动态| 一级毛片久久久久久久久女| 欧美潮喷喷水| 97在线视频观看| 男女那种视频在线观看| 色哟哟·www| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 床上黄色一级片| 国产精品久久电影中文字幕| 亚洲18禁久久av| 亚洲欧美精品综合久久99| 午夜福利高清视频| 午夜福利在线观看吧| 成人特级黄色片久久久久久久| 国产黄色小视频在线观看| h日本视频在线播放| 亚洲18禁久久av| 久久草成人影院| 久久热精品热| 91在线精品国自产拍蜜月| 国产高清三级在线| 国产午夜福利久久久久久| 一个人看视频在线观看www免费| 亚洲国产高清在线一区二区三| av在线天堂中文字幕| 久久精品国产99精品国产亚洲性色| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久精品电影| 床上黄色一级片| 在线免费观看不下载黄p国产| 免费av观看视频| 有码 亚洲区| 久久午夜亚洲精品久久| 久久久精品94久久精品| 久久久久久久久久久丰满| 午夜福利在线在线| 久久久国产成人免费| av黄色大香蕉| 国产精品久久电影中文字幕| 亚洲经典国产精华液单| 欧美bdsm另类| 我要看日韩黄色一级片| 国产亚洲av嫩草精品影院| 亚洲av成人精品一区久久| 97碰自拍视频| 少妇熟女欧美另类| 亚洲丝袜综合中文字幕| 国产精品女同一区二区软件| 精品久久久久久久人妻蜜臀av| 国产黄色小视频在线观看| 波多野结衣巨乳人妻| 亚洲七黄色美女视频| 欧美潮喷喷水| 日日干狠狠操夜夜爽| 不卡视频在线观看欧美| 免费观看的影片在线观看| 九九在线视频观看精品| 国产中年淑女户外野战色| 女同久久另类99精品国产91| 欧美色视频一区免费| 小说图片视频综合网站| 中文字幕精品亚洲无线码一区| av天堂中文字幕网| 国产淫片久久久久久久久| 网址你懂的国产日韩在线| 国产亚洲av嫩草精品影院| 精品不卡国产一区二区三区| 能在线免费观看的黄片| 亚洲精品日韩av片在线观看| 精品久久久久久成人av| 床上黄色一级片| 中文字幕av在线有码专区| 亚洲aⅴ乱码一区二区在线播放| 久久亚洲国产成人精品v| 卡戴珊不雅视频在线播放| 看黄色毛片网站| 久久韩国三级中文字幕| 久久久久免费精品人妻一区二区| 亚洲精品影视一区二区三区av| 91久久精品电影网| 好男人在线观看高清免费视频| 色在线成人网| 最近2019中文字幕mv第一页| 久久精品国产自在天天线| 一区二区三区免费毛片| 国内少妇人妻偷人精品xxx网站| 99久久精品国产国产毛片| 国产男靠女视频免费网站| 老司机影院成人| 人人妻,人人澡人人爽秒播| 18禁在线无遮挡免费观看视频 | 少妇被粗大猛烈的视频| 国产大屁股一区二区在线视频| 午夜福利18| 最近最新中文字幕大全电影3| 亚洲美女黄片视频| 中文字幕免费在线视频6| 精品国产三级普通话版| 高清毛片免费看| 亚洲国产欧美人成| 啦啦啦观看免费观看视频高清| 欧美性感艳星| 欧美最黄视频在线播放免费| 国产精品女同一区二区软件| 日日摸夜夜添夜夜添小说| 久久久精品欧美日韩精品| 国产欧美日韩精品一区二区| 亚洲七黄色美女视频| av.在线天堂| 给我免费播放毛片高清在线观看| 国产午夜精品久久久久久一区二区三区 | 99riav亚洲国产免费| 国产三级中文精品| 精品无人区乱码1区二区| 丰满的人妻完整版| 精品99又大又爽又粗少妇毛片| 精品乱码久久久久久99久播| 久久亚洲国产成人精品v| 少妇裸体淫交视频免费看高清| 国产午夜福利久久久久久| 午夜免费激情av| 无遮挡黄片免费观看| 在线免费观看不下载黄p国产| 男女下面进入的视频免费午夜| 麻豆精品久久久久久蜜桃| av中文乱码字幕在线| 国产乱人偷精品视频| 久久精品影院6| 婷婷亚洲欧美| 欧美在线一区亚洲| 91久久精品电影网| 午夜精品国产一区二区电影 | av女优亚洲男人天堂| 亚洲精品456在线播放app| 国产精品不卡视频一区二区| 婷婷精品国产亚洲av在线| 男插女下体视频免费在线播放| 真实男女啪啪啪动态图| 亚洲av中文av极速乱| 国产精品嫩草影院av在线观看| 国产成人一区二区在线| 国产精品av视频在线免费观看| 成年av动漫网址| 日本五十路高清| 久久99热6这里只有精品| 亚洲精品亚洲一区二区| 日本撒尿小便嘘嘘汇集6| 国产午夜精品久久久久久一区二区三区 | 久久久色成人| 久久精品91蜜桃| 男人狂女人下面高潮的视频| 亚洲美女黄片视频| 青春草视频在线免费观看| 久久这里只有精品中国| 啦啦啦韩国在线观看视频| 日日撸夜夜添| 色5月婷婷丁香| 国产精品国产三级国产av玫瑰| 内地一区二区视频在线| 人人妻,人人澡人人爽秒播| 蜜桃久久精品国产亚洲av| 国产国拍精品亚洲av在线观看| 亚洲中文日韩欧美视频| 男人的好看免费观看在线视频| 成人亚洲欧美一区二区av| aaaaa片日本免费| 国产精品国产三级国产av玫瑰| 小蜜桃在线观看免费完整版高清| 毛片一级片免费看久久久久| 国产成人精品久久久久久| 成熟少妇高潮喷水视频| 国产三级中文精品| 色av中文字幕| 亚洲成人中文字幕在线播放| 亚洲久久久久久中文字幕| 成人漫画全彩无遮挡| 看片在线看免费视频| 久久久久九九精品影院| 亚洲真实伦在线观看| 国产91av在线免费观看| 亚洲性夜色夜夜综合| 久久久精品大字幕| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 免费av毛片视频| 国产av一区在线观看免费| 天天躁夜夜躁狠狠久久av| 91午夜精品亚洲一区二区三区| 免费大片18禁| 免费人成在线观看视频色|