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

    植被指數(shù)對藏北高寒草甸干旱的敏感性分析

    2018-03-05 08:51:30周艷蓮居為民戴聲佩包永康
    江蘇農(nóng)業(yè)科學 2018年1期
    關鍵詞:時間尺度草甸植被指數(shù)

    李 晗,周艷蓮,居為民,戴聲佩,宋 練,包永康

    (1.南京大學國際地球系統(tǒng)科學研究所,江蘇南京 210046;2.南京大學地理與海洋科學學院,江蘇南京 210046;3.江蘇省地理信息技術(shù)重點實驗室,江蘇南京 210046)

    當前人類活動的不斷加劇會對全球大氣循環(huán)的格局造成越來越嚴重的影響。許多大氣循環(huán)模型預測未來全球氣候變化的特征,結(jié)果是極端降水事件和極端干旱事件發(fā)生的頻率會顯著增加[1-4]。植被生長季干旱嚴重會影響植被指數(shù)、地表溫度、土壤濕度、反射率、蒸散等地表生物物理參數(shù)[5]。

    干旱指數(shù)是反映氣候干旱程度的指標,許多學者利用氣象資料、水文資料、土壤資料或衛(wèi)星遙感資料等來建立不同的干旱指數(shù),從而進行干旱監(jiān)測預警。資料顯示常用的干旱指數(shù)有50多種,例如標準化降水指數(shù)(SPI)[6]、Palmer干早指數(shù)(PDSI)[7]、Z指數(shù)[8]、綜合氣象干旱指數(shù)(CI)[9]等。這些指數(shù)在一定程度上為全球或者地區(qū)的干旱監(jiān)測提供了幫助,但都不可避免地需要降水數(shù)據(jù)的地面采集[10-11]。

    與常規(guī)地面干旱監(jiān)測相比,利用遙感手段進行干旱監(jiān)測具有時空分辨率相對較高、傳感器覆蓋面廣、數(shù)據(jù)獲取較為容易、消耗較少人力物力等特點。隨著遙感對地觀測能力的提升,大量研究者利用遙感數(shù)據(jù)不同波段組合獲取植被指數(shù),通過簡化復雜的干旱現(xiàn)象評估干旱的嚴重程度。遙感數(shù)據(jù)可以獲取與植物生長有關的植被指數(shù),如歸一化差值植被指數(shù)(NDVI)和增強型植被指數(shù)(EVI)等。NDVI可以反映年尺度和季節(jié)尺度上植被生長活動的變化,由于它進行了比值和歸一化處理,在一定程度上可以消除如太陽高度角、衛(wèi)星觀測角、地形、云、陰影和大氣條件有關的輻照度變化等的部分影響[12]。EVI利用藍光波段和紅光波段通過大氣氣溶膠時的差別,補償了殘留氣溶膠對紅光的吸收,綜合采用抗大氣植被指數(shù)和抗土壤植被指數(shù),克服了土壤背景的影響和NDVI在植被高覆蓋區(qū)易飽和、植被低覆蓋區(qū)受土壤植被影響較大、對大氣衰減去除不徹底等缺點[13]。Song等利用中國地區(qū)1982—2005年GIMMS NDVI數(shù)據(jù)集和降水數(shù)據(jù)分析認為,在所有氣象要素中,每月NDVI的最大值與降水的相關性最強,尤其是在干旱與半干旱地區(qū)其相關性表現(xiàn)最為明顯[14]。趙偉等采用NDVI時間序列諧波分析法(HANTS)對中分辨率成像光譜儀MODIS產(chǎn)品的EVI數(shù)據(jù)進行去云處理,重構(gòu)EVI數(shù)據(jù)并結(jié)合地面的氣象數(shù)據(jù),分析表明處理后的EVI時間序列數(shù)據(jù)能夠較好地反映干旱對地表植被的影響[12]。但是,NDVI、EVI等均是比較保守的指數(shù),在葉片受到外界環(huán)境脅迫而導致水分含量發(fā)生變化時,NDVI與EVI自身不一定能及時地響應[15-16]。因此,近年來一些與水分相關的植被指數(shù)被越來越多地用于植被水分含量的估計與干旱監(jiān)測。無論是物理輻射傳輸模型還是試驗統(tǒng)計均表明,植物的水分含量在很大程度上影響著葉片在短波紅外波段下的反射率,由于這個區(qū)域是葉片水分光譜的吸收范圍,所以植被在短波紅外的反射率與它們?nèi)~片的水分含量是呈負相關的[16]。研究表明,將短波紅外波段反射率和近紅外波段反射率一起組合建立如歸一化差值水體指數(shù)(NDWI)和地表水分指數(shù)(LSWI)等植被指數(shù)可以對植被水分含量進行監(jiān)測[17]。一些國外學者研究發(fā)現(xiàn),與植被水分有關的植被指數(shù)(如LSWI)在植被干旱監(jiān)測上的敏感性優(yōu)于與植被生長狀況有關的植被指數(shù)[18-19]。Bajgain等以美國季風區(qū)高桿草原為研究區(qū),分析認為當干旱發(fā)生時,LSWI在年際尺度上的下降幅度優(yōu)于NDVI與EVI,當土壤水分含量發(fā)生變化時,LSWI更能對變化做出及時的響應[20]。但目前在國內(nèi),尤其是關于半干旱區(qū)高原草甸不同類型植被指數(shù)對干旱的敏感性研究尚少,不同時間尺度、不同類型植被指數(shù)對干旱追蹤能力的分析也較欠缺,因此利用不同類型植被指數(shù)對半干旱區(qū)高原草甸干旱的敏感性進行系統(tǒng)的分析,就顯得很有必要。

    高寒草甸是青藏高原廣泛分布的典型植被類型之一,面積約為1.2×106km2,是亞洲中部高寒地區(qū)乃至世界高寒地區(qū)極具代表性的植被[21-22]。水分作為高寒草甸這種干旱、半干旱區(qū)草原植物生長發(fā)育的限制性資源,其缺失會直接導致區(qū)域內(nèi)地生物量的減少,而草原生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)中對干旱非常敏感的系統(tǒng)[1]。因此,對高寒草甸的干旱進行分析與評估是至關重要的。目前由于人類活動的加劇,藏北牧區(qū)超載放牧現(xiàn)象十分嚴重,極端天氣事件時有發(fā)生,草地利用壓力很大,因此對類似于當雄草原這種典型高寒草甸植被的干旱監(jiān)測研究顯得很有必要[21]。本研究以高寒草甸干旱的監(jiān)測與評估為目的,以MODIS的反射率數(shù)據(jù)為遙感數(shù)據(jù)源,計算NDVI、EVI、LSWI等植被指數(shù),結(jié)合同時相內(nèi)觀測的降水與土壤濕度數(shù)據(jù),分析在不同時間尺度上,不同類型植被指數(shù)對干旱的敏感性,篩選出適宜高寒草甸干旱監(jiān)測與評估的最佳遙感指數(shù)。

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

    1.1 研究地概況

    西藏自治區(qū)拉薩市當雄縣草原站,地理坐標為 30°29′50.21″N,91°03′58.90″E,地處念青唐古拉山的南緣。當雄草原植被屬于典型的高寒嵩草草甸植被,建群種有高山嵩草、矮嵩草、絲穎針茅,伴生種有多種密叢生嵩草,本類草甸由于超載過牧而退化,針茅類入侵較多,常形成草原草甸,群落覆蓋度約80%。氣候為高原季風氣候,年平均氣溫1.3 ℃,最冷月1月最低平均氣溫-10.4 ℃,最熱月7月10.7 ℃。多年平均降水量450 mm,其中85%集中在6—8月。冰凍期較長,持續(xù)時間為11月至次年1月。年太陽總輻射 7 527.6 MJ/m2,光合有效輻射3 213.3 MJ/m2。土壤類型為高寒草甸土,土壤厚度為30~50 cm。

    1.2 降水與土壤濕度數(shù)據(jù)

    采用2個不同深度(5、10 cm)的土壤水分含量(SWC)來表征土壤的濕潤程度。同時,為了與MODIS地表反射率數(shù)據(jù)匹配,對站點獲取的日平均降水和土壤水分含量數(shù)據(jù)分別作求和與平均處理,獲得時間尺度為8 d的數(shù)據(jù),進一步計算得到月平均降水、年平均降水和土壤水分含量。

    1.3 遙感數(shù)據(jù)源及數(shù)據(jù)處理

    MODIS數(shù)據(jù)可以同時獲得可見光至近紅外的36個波段,常被用于植被遙感研究。本研究通過美國航空航天宇航局的戈達德航天中心數(shù)據(jù)網(wǎng)站(LAADS,https://ladsweb.nascom.nasa.gov/)下載覆蓋當雄通量觀測站點每8 d的合成地表反射率產(chǎn)品MOD09A1,時間序列為2004年至2012年,空間分辨率為500 m,該產(chǎn)品主要包括MODIS傳感器1~7波段的地表反射率數(shù)據(jù)和各個波段的質(zhì)量控制數(shù)據(jù),同時還包括反射率數(shù)據(jù)狀態(tài)質(zhì)量保證數(shù)據(jù)。利用MODIS數(shù)據(jù)的處理工具(MRT)將投影方式由Sinusoidal投影轉(zhuǎn)換為UTM/WGS 84投影,并通過最近鄰法重采樣,獲取站點中心單像元紅光波段(620~670 nm)、近紅外波段1(841~876 nm)、藍光波段(459~479 nm)以及短波紅外波段1(1 628~1 652 nm)反射率數(shù)據(jù)。根據(jù)公式(1~3)分別計算不同時間尺度的NDVI、EVI和LSWI。

    (1)

    (2)

    (3)

    式中:ρX為MOIDS數(shù)據(jù)X波段的反射率;公式(2)中C1和C2是氣溶膠等大氣影響物的校正系數(shù),L為背景(土壤)調(diào)整系數(shù),G為常量,取值分別為C1=6,C2=7.5,L=1,G=1。

    由于LSWI與天氣條件有關,無法像葉面積指數(shù)(LAI)和光合作用有效輻射比率(FPAR)一樣進行數(shù)據(jù)平滑處理[21]。本研究結(jié)合質(zhì)量控制數(shù)據(jù),采用兩步差補方法處理時間序列中受云影響、質(zhì)量較差的植被指數(shù)數(shù)據(jù)[22]。NDVI、EVI和LSWI均未作平滑處理。

    2 結(jié)果與分析

    2.1 年時間尺度敏感性分析

    圖1表明,2004—2012年,拉薩當雄高原草甸站點植被生長季的年平均降水量為428 mm(±126 mm)。年平均降水量大于2004—2012年平均降水量與標準差之和(554 mm)的年份被認為是濕潤年,而年平均降水量小于2004—2012年平均降水量與標準差之差(302 mm)的年份被認為是干旱年。拉薩當雄高原草甸站點2006、2009年的降水量分別只有241、291 mm,因此為干旱年。而2008年的降水量為648 mm,因此為濕潤年。干旱年份(2006、2009年)植被生長季(5—9月)的不同深度(5、10 cm)土壤水分含量均<9%,而濕潤年份(2008年)土壤水分含量為12.6%。

    圖2表明,干旱年份(2006、2009年)植被生長季平均NDVI(NDVIavg)、EVI(EVIavg)和LSWI(LSWIavg)在分析的時間序列中均較低(NDVIavg<0.28、EVIavg<0.18、LSWIavg<0),而在濕潤年份(2008年)植被生長季NDVI、EVI和LSWI均達到最大值(NDVIavg=0.41、EVIavg=0.28、LSWIavg=0.07)。

    干旱年份(2006、2009年)植被生長季植被指數(shù)最大值NDVI(NDVImax)、EVI(EVImax)和LSWI(LSWImax)均取得最小值(NDVImax<0.44,EVImax<0.29,LSWImax<0.07),在濕潤年份(2008年)植被生長季植被指數(shù)只有最大值NDVI和LSWI達到峰值(NDVImax=0.64,LSWImax=0.20),而EVImax僅為0.48,遠遠小于2011年(EVImax=0.67)。

    植被指數(shù)均值和最大值的變化幅度可以用來表征植被指數(shù)對干旱的相對敏感性。變化幅度用當年植被指數(shù)的均值或最大值減去多年相應的平均值表示。變化幅度百分比用當年植被指數(shù)的均值或最大值減去多年相應的平均值后的結(jié)果再除以多年相應的平均值表示。如圖3所示,2006年LSWIavg與2004—2012年平均值相比,減少了0.052(316%),減少幅度最大。2009年的結(jié)果與2006年相似,LSWIavg、EVIavg、NDVIavg

    分別比2004—2012年平均值減少0.037(227%)、0.068(16%)、0.054(29%)。植被生長季LSWImax與NDVImax、EVImax相比也表現(xiàn)出同樣的敏感性,2006年LSWImax減少0.13(96%),均大于EVImax減少的0.14(35%)與NDVImax減少的0.11(22%)。同樣的,2009年LSWImax減少0.07(52%),均大于EVImax減少的0.13(32%)與NDVImax減少的0.08(16%)。

    濕潤年份(2008年)LSWI持續(xù)大于0的時間最大,而2006年僅有1個LSWI大于0。由于氣溫低,海拔高,4月中旬生長季才剛剛開始。經(jīng)過春季生長之后,植被生長的速率很大一部分取決于土壤水分含量是否充足。若土壤水分含量不充足,植被的生長就會受到抑制。相反,干旱年份LSWI小于0的天數(shù)多于濕潤年份與正常年份,而濕潤年份LSWI大于0.1的時間多于干旱年份與正常年份。如圖4所示,LSWI小于0的天數(shù)隨降水量的增加整體是下降的。2006年LSWI小于0的天數(shù)為120 d,大于2007年的72 d和2008年的 40 d,而2006年LSWI大于0.1的天數(shù)僅為8 d,而2008年為88 d,大于2007年的56 d。

    2.2 月時間尺度敏感性分析

    在分析以植被指數(shù)為基礎的干旱指標對降水的敏感性時,易忽略植被指數(shù)和降水的同時相季節(jié)變化,可能對其相關性產(chǎn)生影響。本研究為消除季節(jié)變化引起的降水與植被指數(shù)的正弦波動,對月時間尺度的降水數(shù)據(jù)和植被指數(shù)數(shù)據(jù)均作了距平處理。具體做法是用當月降水的總累積量和植被指數(shù)的平均值減去多年的月平均值。圖5顯示,當前月份LSWI與降水距平的相關性最高,r2達到0.32(P<0.001),其次是NDVI,r2為0.16(P<0.1),EVI的r2最低,僅為0.05。而在1個月滯后的情況下,NDVI與EVI的r2分別為0.20(P<0.01)和0.24(P<0.01),與當前月份降水與植被指數(shù)距平相關性相比有一定的提高,而LSWI和降水距平的r2僅為0.17(P<0.1),有一定的下降。

    3 結(jié)論與討論

    相比利用氣象、水文、土壤等資料建立干旱指數(shù)用于干旱監(jiān)測等方法而言,以植被指數(shù)為核心的遙感干旱監(jiān)測指數(shù)具有很多優(yōu)勢。其成本低、時空分辨率較高、對降水能夠及時響應,而且它對地表植被的葉綠素含量、水分含量均比較敏感,優(yōu)勢明顯。植被在不同波段的光譜特征能反映指數(shù)的不同要素或者某種特征的差異,植被指數(shù)利用多個光譜波段的反射率經(jīng)分析運算(線性或非線性組合方式),產(chǎn)生對植被長勢、水分含量等有一定指示意義的數(shù)值。這是利用植被指數(shù)進行干旱監(jiān)測與評估的理論基礎[20]。Chandrasekar等以印度南岸的城市——安得拉邦為研究區(qū),以整個地區(qū)MODIS反射率數(shù)據(jù)為基礎,結(jié)合氣象數(shù)據(jù)分析認為,LSWI在低降水地區(qū)與降水的相關性更好,優(yōu)于NDVI;而在高降水地區(qū),LSWI與降水的相關性要差于NDVI,且都存在滯后現(xiàn)象[23]。但是該研究沒有考慮因季節(jié)變化引起的降水和植被指數(shù)在年內(nèi)的正弦波動,這種季節(jié)特性在一定程度上影響了結(jié)果的可靠性。本研究利用NDVI、EVI、LSWI這3種指數(shù)對高寒草甸進行干旱監(jiān)測,同時為消除季節(jié)變化引起的降水與植被指數(shù)的正弦波動,對降水和植被指數(shù)在月尺度上進行距平處理,結(jié)果顯示月時間尺度上LSWI對降水的敏感性優(yōu)于NDVI與EVI,且LSWI不存在滯后現(xiàn)象。這與草原生態(tài)系統(tǒng)下墊面比較均一、受背景反射率影響較小等原因有關。2006年夏季站點周圍發(fā)生了嚴重的干旱現(xiàn)象,LSWI在8、9月已經(jīng)顯示出明顯的下降趨勢,而NDVI和EVI直到10月才開始下降。進一步說明LSWI對降水的敏感性在月尺度上較NDVI、EVI更好,且不存在滯后現(xiàn)象(圖6)。

    Bajgain以美國季風區(qū)的2個高桿草原站點為研究區(qū)、以遙感獲取的LSWI為基礎對干旱進行的分級(LSWI>0.1;00.1(0)表示正常。海北高寒草甸位于半干旱區(qū),正常年份均值處于0

    以植被指數(shù)為核心的遙感干旱監(jiān)測指數(shù)雖然具有很多優(yōu)勢,但在其大范圍應用之前應該對不同的植被類型進行遙感干旱檢測研究,以確定合適的研究指標和研究手段。本研究的目的就是探討藏北典型高寒草甸干旱檢測的最佳遙感指數(shù),以期更精確地對高寒草甸干旱事件的發(fā)生進行有效地監(jiān)測與評估。當然利用植被指數(shù)進行干旱監(jiān)測與評估也有諸多限制,例如遙感觀測資料的時間尺度,在天時間尺度或者更高時間尺度上的實時監(jiān)測目前來講仍存在極大挑戰(zhàn),不過隨著區(qū)域性小衛(wèi)星的陸續(xù)發(fā)射升空,也為干旱的實時遙感監(jiān)測提供了可能。

    對干旱指數(shù)的評價應該從多方面進行,綜合考慮。本研究評價干旱指數(shù)的效果時,采用被廣泛使用的相關系數(shù)r2和顯著性水平P值作為主要評價指標,也定性地分析了在干旱發(fā)生當月,最早作出反應的植被指數(shù),并驗證了以LSWI為基礎的干旱指標在半干旱區(qū)的適應性問題,進一步說明LSWI對草原干旱事件的響應優(yōu)于NDVI與EVI。

    本研究利用MODIS反射率數(shù)據(jù)計算NDVI、EVI和LSWI,結(jié)合降水與土壤濕度數(shù)據(jù)對藏北高寒草甸干旱進行監(jiān)測,干旱事件發(fā)生時植被生長受到了很大的限制,3種植被指數(shù)在干旱年份(2006年)減小幅度也基本相同,LSWI略為敏感。月尺度上LSWI與降水距平的相關性最高,r2達到0.32,其次是NDVI與EVI,且LSWI不存在滯后現(xiàn)象。LSWI<0的天數(shù)隨降水量的增加呈明顯的減少趨勢。利用LSWI對干旱進行分級(LSWI>0.1;0

    [1]張 彬,朱建軍,劉華民,等.極端降水和極端干旱事件對草原生態(tài)系統(tǒng)的影響[J].植物生態(tài)學報,2014,38(9):1008-1018.

    [2]International Panel on Climate Change.Climate change 2007:impacts,adaption and vulnerability[R].New York:Cambridge University Press,2007.

    [3]Easterling D R,Meehl G A,Parmesan C,et al.Climate extremes:observations,modeling,and impacts[J].Science,2000,289(5487):2068-2074.

    [4]Meehl G A,Tebaldi C.More intense,more frequent,and longer lasting heat waves in the 21st century[J].Science,2004,305(5686):994-997.

    [5]Ghulam A,Li Z L,Qin Q,et al.Exploration of the spectral space basedon vegetation index and lbedo for surface drought estimation[J].Journal of Applied Remote Sensing,2007,1(1):13529-13529.

    [6]王勁松,李耀輝,王潤元,等.我國氣象干旱研究進展評述[J].干旱氣象,2012,30(4):497-508.

    [7]楊小利,楊興國,馬鵬里,等.PDSI在甘肅中東部地區(qū)的修正和應用[J].地球科學進展,2005,20(9):1022-1028.

    [8]王勁松,黃玉霞,馮建英,等.徑流量Z指數(shù)與Palmer指數(shù)對河西干旱的監(jiān)測[J].應用氣象學報,2009,20(4):471-477.

    [9]李樹巖,劉榮花,師麗魁,等.基于CI指數(shù)的河南省近40年干旱特征分析[J].干旱氣象,2009,27(2):97-102.

    [10]王春林,郭 晶,薛麗芳,等.改進的綜合氣象干旱指數(shù)CInew及其適用性分析[J].中國農(nóng)業(yè)氣象,2011,32(4):621-626.

    [11]張紅衛(wèi),陳懷亮,申雙和,等.基于表層水分含量指數(shù)(SWCI)的土壤干旱遙感監(jiān)測[J].遙感技術(shù)與應用,2008,23(6):624-628.

    [12]趙 偉,李召良.利用MODIS/EVI時間序列數(shù)據(jù)分析干旱對植被的影響[J].地理科學進展,2007,26(6):40-46.

    [13]彭代亮,黃敬峰,王秀珍.基于MODIS-EVI區(qū)域植被季節(jié)變化與氣象因子的關系[J].應用生態(tài)學報,2007,18(5):983-989.

    [14]Song Y,Ma M G.A statistical analysis of the relationship between climatic factors and the normalized difference vegetation index in China[J].International Journal of Remote Sensing,2011,32(14):3947-3965.

    [15]劉小磊,覃志豪.NDWI與NDVI指數(shù)在區(qū)域干旱監(jiān)測中的比較分析——以2003年江西夏季干旱為例[J].遙感技術(shù)與應用,2007,22(5):608-612.

    [16]劉立文,張吳平,段永紅,等.TVDI模型的農(nóng)業(yè)旱情時空變化遙感應用[J].生態(tài)學報,2014,34(13):3704-3711.

    [17]Chen D Y,Huang J F,Jackson T J.Vegetation water content estimation for corn and soybeans using spectral indices derived from MODIS near- and short-wave infrared bands[J].Remote Sensing of Enviroment,2005,98(2):225-236.

    [18]張友水,謝元禮.MODIS影像的NDVI和LSWI植被水分含量估算[J].地理科學,2008,28(1):72-76.

    [19]劉元亮,李 艷,吳劍亮.基于LSWI和NDVI時間序列的水田信息提取研究[J].地理與地理信息科學,2015,31(3):32-37.

    [20]Bajgain R,Xiao X M,Wagle P,et al.Sensitivity analysis of vegetation indices to drought over two tallgrass prairie sites[J].ISPRS Journal of Photogrammetry and Remote Sensing,2015(108):151-160.

    [21]周宇庭,付 剛,沈振西,等.藏北典型高寒草甸地上生物量的遙感估算模型[J].草業(yè)學報,2013,22(1):120-129.

    [22]Zhang Y,Song C H,Sun G,et al.Understanding moisture stress on light use efficiency across terrestrial ecosystems based on global flux and remote sensing data[J].Journal of Geophysical Research:Biogeosciences,2015,120(10):2053-2066.

    [23]Chandrasekar K,Sesha S,Roy P S,et al.Land surface water index (LSWI) response to rainfall and NDVI using the MODIS vegetation index product[J].International Journal of Remote Sensing,2010,31(15):3987-4005.

    猜你喜歡
    時間尺度草甸植被指數(shù)
    高山草甸
    時間尺度上非完整系統(tǒng)的Noether準對稱性與守恒量
    The world's narrowest river
    山地草甸
    時間尺度上Lagrange 系統(tǒng)的Hojman 守恒量1)
    力學學報(2021年10期)2021-12-02 02:32:04
    交直流混合微電網(wǎng)多時間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    大連市暴雨多時間尺度研究分析
    91午夜精品亚洲一区二区三区| 亚洲精品第二区| 丰满乱子伦码专区| 亚洲精品久久久久久婷婷小说| 男女边摸边吃奶| 在线观看av片永久免费下载| 高清黄色对白视频在线免费看 | 中文字幕av电影在线播放| 在线观看美女被高潮喷水网站| 18禁在线无遮挡免费观看视频| 另类精品久久| 精品视频人人做人人爽| 大话2 男鬼变身卡| 精品久久国产蜜桃| 国产精品国产三级国产专区5o| 精品亚洲成国产av| 亚洲欧美一区二区三区黑人 | 少妇被粗大的猛进出69影院 | 免费人成在线观看视频色| 精品国产露脸久久av麻豆| 五月伊人婷婷丁香| 噜噜噜噜噜久久久久久91| 下体分泌物呈黄色| 另类亚洲欧美激情| 亚洲精品,欧美精品| 黑人猛操日本美女一级片| 男人和女人高潮做爰伦理| 精品人妻一区二区三区麻豆| 熟女电影av网| 中文精品一卡2卡3卡4更新| 亚洲va在线va天堂va国产| 丝瓜视频免费看黄片| 人妻制服诱惑在线中文字幕| 人人妻人人澡人人看| 日本爱情动作片www.在线观看| 少妇 在线观看| 青春草视频在线免费观看| 中文字幕人妻丝袜制服| 国产成人精品婷婷| 国产美女午夜福利| 日韩成人伦理影院| 欧美日韩视频精品一区| 日日爽夜夜爽网站| 精品亚洲乱码少妇综合久久| 欧美三级亚洲精品| 国产男女内射视频| a 毛片基地| 国产色婷婷99| 人妻夜夜爽99麻豆av| 日韩大片免费观看网站| 五月天丁香电影| 丝袜喷水一区| av.在线天堂| 亚洲一级一片aⅴ在线观看| 性色avwww在线观看| 国产在线一区二区三区精| 性高湖久久久久久久久免费观看| 色婷婷av一区二区三区视频| 国产在视频线精品| 久久狼人影院| 中文天堂在线官网| 亚洲av电影在线观看一区二区三区| 狂野欧美白嫩少妇大欣赏| h视频一区二区三区| 国产免费一级a男人的天堂| 免费不卡的大黄色大毛片视频在线观看| 一区二区三区免费毛片| 国产色爽女视频免费观看| 亚洲成人av在线免费| 国产伦精品一区二区三区视频9| 久久精品久久精品一区二区三区| 久久久久久久国产电影| 久久国产精品大桥未久av | 又黄又爽又刺激的免费视频.| 欧美区成人在线视频| 亚洲欧美成人综合另类久久久| 高清不卡的av网站| 精品亚洲成国产av| 日韩av在线免费看完整版不卡| 97在线视频观看| 国产美女午夜福利| 国产一区二区三区综合在线观看 | 亚洲精品乱码久久久久久按摩| 欧美少妇被猛烈插入视频| 国产黄色视频一区二区在线观看| 人妻 亚洲 视频| 久久久亚洲精品成人影院| 国产精品三级大全| 国产精品伦人一区二区| 大陆偷拍与自拍| 99热网站在线观看| 久久久久精品性色| 亚洲在久久综合| av免费观看日本| 亚洲,一卡二卡三卡| 亚洲人与动物交配视频| 涩涩av久久男人的天堂| 国产在线视频一区二区| 制服丝袜香蕉在线| 免费人妻精品一区二区三区视频| 桃花免费在线播放| 99九九线精品视频在线观看视频| 丰满少妇做爰视频| 麻豆乱淫一区二区| 国产色婷婷99| 一级,二级,三级黄色视频| 国产精品一区二区在线不卡| 国产91av在线免费观看| 熟女电影av网| 久久午夜综合久久蜜桃| 自线自在国产av| 丰满迷人的少妇在线观看| 婷婷色综合www| 国产爽快片一区二区三区| 夜夜爽夜夜爽视频| a级一级毛片免费在线观看| 黄色配什么色好看| 国产精品久久久久久久电影| 日本色播在线视频| 国产亚洲午夜精品一区二区久久| 日韩一区二区三区影片| 少妇高潮的动态图| 精品久久久久久久久av| 久热久热在线精品观看| 全区人妻精品视频| 国产日韩一区二区三区精品不卡 | 少妇 在线观看| 成年女人在线观看亚洲视频| 插逼视频在线观看| 91午夜精品亚洲一区二区三区| 国产欧美日韩综合在线一区二区 | 成年人午夜在线观看视频| 欧美精品亚洲一区二区| 十八禁高潮呻吟视频 | 一本大道久久a久久精品| 久久精品久久久久久久性| 亚洲三级黄色毛片| 久久99蜜桃精品久久| 免费观看无遮挡的男女| 在线观看一区二区三区激情| 国产黄频视频在线观看| 一本久久精品| 亚洲性久久影院| 最新中文字幕久久久久| a 毛片基地| 欧美+日韩+精品| 成人国产麻豆网| 99久久精品国产国产毛片| 亚洲三级黄色毛片| 男女免费视频国产| 有码 亚洲区| 国产精品一区www在线观看| 日韩成人伦理影院| 一本色道久久久久久精品综合| 最近中文字幕高清免费大全6| 麻豆乱淫一区二区| 黄色日韩在线| 国产亚洲5aaaaa淫片| 中文精品一卡2卡3卡4更新| 免费观看的影片在线观看| 午夜精品国产一区二区电影| 日韩视频在线欧美| 日韩熟女老妇一区二区性免费视频| 最新的欧美精品一区二区| 一级a做视频免费观看| 高清欧美精品videossex| 91成人精品电影| 中文字幕人妻丝袜制服| 精品少妇久久久久久888优播| 国产深夜福利视频在线观看| 在线观看av片永久免费下载| 美女内射精品一级片tv| 成年女人在线观看亚洲视频| 26uuu在线亚洲综合色| 插阴视频在线观看视频| 不卡视频在线观看欧美| 啦啦啦在线观看免费高清www| 久久人人爽人人片av| 日本黄色日本黄色录像| 九色成人免费人妻av| 韩国高清视频一区二区三区| 久久人妻熟女aⅴ| 亚洲国产精品一区三区| 国产精品麻豆人妻色哟哟久久| 男女边摸边吃奶| 欧美国产精品一级二级三级 | 国产在线免费精品| 天天躁夜夜躁狠狠久久av| av福利片在线| 成人特级av手机在线观看| 国产精品熟女久久久久浪| 又爽又黄a免费视频| h日本视频在线播放| 草草在线视频免费看| 欧美精品亚洲一区二区| 国产在线免费精品| 在线观看国产h片| 狠狠精品人妻久久久久久综合| 日本-黄色视频高清免费观看| 狂野欧美激情性xxxx在线观看| 亚洲精品国产av蜜桃| 全区人妻精品视频| 免费观看无遮挡的男女| 亚洲欧美成人精品一区二区| 黑人高潮一二区| 99热这里只有是精品在线观看| 久久国产精品大桥未久av | 高清视频免费观看一区二区| 有码 亚洲区| 成年av动漫网址| 久久久久久久久久久久大奶| 午夜老司机福利剧场| av视频免费观看在线观看| 亚洲天堂av无毛| 色婷婷久久久亚洲欧美| 黄色配什么色好看| 亚洲自偷自拍三级| 最后的刺客免费高清国语| 色视频www国产| 国产伦精品一区二区三区四那| 精品视频人人做人人爽| 亚洲成人av在线免费| 少妇丰满av| 丰满乱子伦码专区| 国产成人a∨麻豆精品| 如何舔出高潮| 国产成人一区二区在线| 亚洲精品视频女| 国产欧美另类精品又又久久亚洲欧美| 免费观看的影片在线观看| 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www | 免费观看在线日韩| 亚洲精品久久午夜乱码| 久久毛片免费看一区二区三区| 一区二区av电影网| 超碰97精品在线观看| 日本-黄色视频高清免费观看| 午夜免费鲁丝| 赤兔流量卡办理| 亚洲av男天堂| 中文字幕免费在线视频6| 成人毛片60女人毛片免费| 成人18禁高潮啪啪吃奶动态图 | 成人午夜精彩视频在线观看| 一级毛片 在线播放| 乱码一卡2卡4卡精品| 久久久午夜欧美精品| 国产免费一区二区三区四区乱码| 国产欧美亚洲国产| 亚洲av电影在线观看一区二区三区| 一个人免费看片子| 亚洲av成人精品一二三区| 51国产日韩欧美| 国产黄片美女视频| 综合色丁香网| 蜜桃在线观看..| 最近最新中文字幕免费大全7| 麻豆乱淫一区二区| 日本黄大片高清| 综合色丁香网| 久久婷婷青草| 我的女老师完整版在线观看| 内地一区二区视频在线| 日日啪夜夜撸| 国产免费视频播放在线视频| videos熟女内射| 国产免费福利视频在线观看| 亚洲精品日韩av片在线观看| 视频中文字幕在线观看| 美女视频免费永久观看网站| 伊人久久精品亚洲午夜| 一区在线观看完整版| 国产精品久久久久久av不卡| 亚洲国产精品国产精品| 男女免费视频国产| 久久久久人妻精品一区果冻| h视频一区二区三区| 亚洲av成人精品一二三区| 国产精品国产三级国产专区5o| 永久网站在线| 午夜福利影视在线免费观看| 黑丝袜美女国产一区| 搡老乐熟女国产| 中文字幕人妻熟人妻熟丝袜美| 深夜a级毛片| 人妻少妇偷人精品九色| 日本黄大片高清| 一区二区三区免费毛片| 女人久久www免费人成看片| 中国国产av一级| 亚洲三级黄色毛片| 亚洲av福利一区| 秋霞在线观看毛片| 亚洲久久久国产精品| 亚洲精品久久久久久婷婷小说| 少妇的逼好多水| 人妻系列 视频| 日日摸夜夜添夜夜添av毛片| 免费人妻精品一区二区三区视频| 国精品久久久久久国模美| tube8黄色片| 久久久亚洲精品成人影院| 亚洲情色 制服丝袜| 热re99久久精品国产66热6| 日韩成人伦理影院| 国产精品无大码| 黄色毛片三级朝国网站 | 日韩中字成人| 亚洲欧美一区二区三区国产| 国内少妇人妻偷人精品xxx网站| 亚洲欧美清纯卡通| 精品国产一区二区三区久久久樱花| 亚洲,一卡二卡三卡| 亚洲精品一二三| 色哟哟·www| 中文在线观看免费www的网站| 26uuu在线亚洲综合色| 韩国高清视频一区二区三区| 国产精品福利在线免费观看| 性色av一级| 亚洲激情五月婷婷啪啪| 午夜激情福利司机影院| 日日啪夜夜撸| 国产亚洲一区二区精品| 国产黄频视频在线观看| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| 久久精品久久久久久久性| 天天躁夜夜躁狠狠久久av| 国产老妇伦熟女老妇高清| 肉色欧美久久久久久久蜜桃| 大香蕉97超碰在线| 欧美三级亚洲精品| 在现免费观看毛片| 纵有疾风起免费观看全集完整版| 人妻少妇偷人精品九色| 91午夜精品亚洲一区二区三区| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 成人无遮挡网站| 三级经典国产精品| 王馨瑶露胸无遮挡在线观看| 亚洲成人av在线免费| 热re99久久国产66热| 黄色一级大片看看| 嫩草影院入口| 国产男女超爽视频在线观看| 精品国产国语对白av| 成人免费观看视频高清| 91精品伊人久久大香线蕉| 十八禁网站网址无遮挡 | h日本视频在线播放| 啦啦啦视频在线资源免费观看| 亚洲性久久影院| 一级二级三级毛片免费看| 国产av一区二区精品久久| 久久久国产一区二区| 成人综合一区亚洲| 九色成人免费人妻av| 美女中出高潮动态图| 啦啦啦在线观看免费高清www| 久久久久久久亚洲中文字幕| 国产亚洲91精品色在线| 免费人成在线观看视频色| 一级片'在线观看视频| 国产精品麻豆人妻色哟哟久久| 日韩大片免费观看网站| 中文字幕久久专区| 自拍欧美九色日韩亚洲蝌蚪91 | 我要看日韩黄色一级片| 人妻一区二区av| 一级二级三级毛片免费看| 亚洲国产精品999| 搡老乐熟女国产| 久久亚洲国产成人精品v| 成人漫画全彩无遮挡| 国产精品不卡视频一区二区| 日韩成人av中文字幕在线观看| 中文欧美无线码| 18禁在线播放成人免费| 男女无遮挡免费网站观看| 日韩伦理黄色片| 国产美女午夜福利| 中文字幕人妻丝袜制服| 久久6这里有精品| 久久这里有精品视频免费| av免费观看日本| 欧美日韩av久久| 人人妻人人澡人人爽人人夜夜| 精品一区在线观看国产| 久久久久久久亚洲中文字幕| 人人妻人人澡人人看| 国产中年淑女户外野战色| 美女脱内裤让男人舔精品视频| 亚洲精品国产成人久久av| 中文字幕人妻熟人妻熟丝袜美| 只有这里有精品99| 麻豆成人午夜福利视频| 99精国产麻豆久久婷婷| 美女xxoo啪啪120秒动态图| 少妇精品久久久久久久| 亚洲av成人精品一二三区| 99热6这里只有精品| 久久久久国产网址| av天堂久久9| 欧美一级a爱片免费观看看| 久久国产乱子免费精品| 国产黄色视频一区二区在线观看| 亚洲欧美一区二区三区黑人 | 精品国产露脸久久av麻豆| 美女cb高潮喷水在线观看| 日本-黄色视频高清免费观看| 曰老女人黄片| 久久热精品热| 麻豆精品久久久久久蜜桃| 亚洲欧美精品专区久久| 久久久久久久国产电影| 色吧在线观看| 美女中出高潮动态图| 97超视频在线观看视频| 中文字幕亚洲精品专区| 免费高清在线观看视频在线观看| 中文乱码字字幕精品一区二区三区| 亚洲国产日韩一区二区| 免费久久久久久久精品成人欧美视频 | 亚洲欧美中文字幕日韩二区| 国产精品一区二区在线观看99| 成人18禁高潮啪啪吃奶动态图 | 不卡视频在线观看欧美| 嫩草影院入口| 大话2 男鬼变身卡| 日韩 亚洲 欧美在线| 男的添女的下面高潮视频| 国产亚洲精品久久久com| 日韩熟女老妇一区二区性免费视频| 色婷婷av一区二区三区视频| 国产视频内射| 亚洲国产精品一区二区三区在线| 国产成人a∨麻豆精品| 91精品一卡2卡3卡4卡| 国产男人的电影天堂91| 亚洲精品一区蜜桃| 免费观看a级毛片全部| 色婷婷久久久亚洲欧美| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久成人av| 制服丝袜香蕉在线| 国产日韩欧美视频二区| av福利片在线| 又黄又爽又刺激的免费视频.| 黑人高潮一二区| 亚洲人成网站在线观看播放| 亚洲精品久久午夜乱码| 九九在线视频观看精品| 曰老女人黄片| 国产 一区精品| 免费黄网站久久成人精品| 亚洲国产精品专区欧美| 在线观看免费视频网站a站| 噜噜噜噜噜久久久久久91| 国产在线视频一区二区| 日韩免费高清中文字幕av| 日韩欧美 国产精品| 国产精品一二三区在线看| 伦理电影免费视频| 欧美 亚洲 国产 日韩一| 女性生殖器流出的白浆| 免费观看性生交大片5| 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 日韩视频在线欧美| 成人国产av品久久久| a级片在线免费高清观看视频| 国产亚洲欧美精品永久| 成人美女网站在线观看视频| 我的女老师完整版在线观看| av在线app专区| 国产真实伦视频高清在线观看| 97精品久久久久久久久久精品| 看十八女毛片水多多多| 国产亚洲91精品色在线| 在线观看www视频免费| 国产白丝娇喘喷水9色精品| 亚洲自偷自拍三级| 狂野欧美激情性bbbbbb| 日韩成人伦理影院| 26uuu在线亚洲综合色| 免费观看的影片在线观看| 成人美女网站在线观看视频| 国产黄频视频在线观看| 久久精品久久精品一区二区三区| 久久久久久久久久久丰满| 国产精品一区二区性色av| 91久久精品国产一区二区成人| 亚洲在久久综合| 人人澡人人妻人| 国产欧美日韩精品一区二区| 爱豆传媒免费全集在线观看| 新久久久久国产一级毛片| 免费播放大片免费观看视频在线观看| 在线观看www视频免费| 国产黄色免费在线视频| 少妇人妻精品综合一区二区| 成人毛片a级毛片在线播放| 国内精品宾馆在线| 精品熟女少妇av免费看| 最后的刺客免费高清国语| 人人妻人人澡人人爽人人夜夜| 黄色毛片三级朝国网站 | 国产精品女同一区二区软件| 亚洲无线观看免费| 天堂8中文在线网| 国产av一区二区精品久久| 久久久精品94久久精品| 亚洲人成网站在线观看播放| 国产日韩一区二区三区精品不卡 | 极品少妇高潮喷水抽搐| 高清午夜精品一区二区三区| 久久热精品热| av国产久精品久网站免费入址| 在线播放无遮挡| 亚洲内射少妇av| 国产成人精品一,二区| 女人精品久久久久毛片| 亚洲精品乱码久久久v下载方式| 亚洲精品乱码久久久久久按摩| 99热6这里只有精品| av网站免费在线观看视频| 丰满少妇做爰视频| av.在线天堂| 国产亚洲精品久久久com| 日本黄大片高清| 久久国内精品自在自线图片| 寂寞人妻少妇视频99o| 精品国产国语对白av| 国产女主播在线喷水免费视频网站| 亚洲国产精品成人久久小说| a级毛片在线看网站| 简卡轻食公司| 三级国产精品欧美在线观看| 在线观看免费高清a一片| 亚洲成人av在线免费| 五月天丁香电影| 肉色欧美久久久久久久蜜桃| 在线观看免费高清a一片| 中文字幕人妻丝袜制服| 乱系列少妇在线播放| 亚洲精品视频女| 超碰97精品在线观看| 草草在线视频免费看| 99热这里只有精品一区| 一区二区av电影网| 一个人免费看片子| 欧美国产精品一级二级三级 | 视频区图区小说| 亚洲精品456在线播放app| 亚洲伊人久久精品综合| 亚洲三级黄色毛片| 久久久国产一区二区| 久久久精品免费免费高清| 国产亚洲91精品色在线| 免费黄网站久久成人精品| 亚洲av日韩在线播放| 亚洲欧美中文字幕日韩二区| 2022亚洲国产成人精品| 久久久久久久久久久丰满| av在线观看视频网站免费| 国产精品伦人一区二区| 少妇被粗大猛烈的视频| 亚洲国产精品国产精品| 亚洲综合色惰| 久久精品国产亚洲av天美| freevideosex欧美| 久久久久久人妻| 久久精品国产自在天天线| 午夜精品国产一区二区电影| 人人妻人人澡人人看| 精品亚洲成国产av| 欧美精品国产亚洲| 国产成人精品一,二区| 男女国产视频网站| 看免费成人av毛片| 免费大片18禁| 黑丝袜美女国产一区| 狂野欧美激情性bbbbbb| 欧美日韩精品成人综合77777| 久久久a久久爽久久v久久| 国产视频首页在线观看| 免费大片黄手机在线观看| 午夜福利影视在线免费观看| 亚洲一区二区三区欧美精品| 久久久久久人妻| 高清av免费在线| 日本黄色日本黄色录像| 免费观看的影片在线观看| 久久99蜜桃精品久久| 99久久人妻综合| 国产男女内射视频| 黄色一级大片看看| 女性被躁到高潮视频| freevideosex欧美| 2021少妇久久久久久久久久久| 免费人妻精品一区二区三区视频| 国产在线一区二区三区精| 国国产精品蜜臀av免费| 大香蕉久久网| av女优亚洲男人天堂| 亚洲熟女精品中文字幕| 丁香六月天网| 中国国产av一级| 免费看光身美女| 亚洲国产最新在线播放| 国产av码专区亚洲av| 简卡轻食公司|