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

    2007—2014年錫林浩特草原返青期變化的遙感監(jiān)測

    2016-03-12 15:22陳玫琪郭秋成徐媛張靖
    天津農(nóng)業(yè)科學(xué) 2016年3期

    陳玫琪++郭秋成++徐媛++張靖

    摘 要:草原返青期是指在溫度和水分條件達到適宜狀態(tài)時,牧草開始萌發(fā)、變綠、生長的時間。通過監(jiān)測草原返青期不僅可以預(yù)測年度草原植被長勢及產(chǎn)量,還可為草原管理部門合理安排季節(jié)畜牧業(yè)生產(chǎn)的實施提供科學(xué)依據(jù)。本研究選取內(nèi)蒙古錫林浩特市草原區(qū)作為典型的研究區(qū),采用非對稱高斯函數(shù)對2007—2014年該區(qū)NDVI時間序列濾波,采用動態(tài)閾值法對該地區(qū)返青期的變化和3—4月份平均氣溫和累計降水量的分析。結(jié)果顯示:2007—2014年錫林浩特草原區(qū)各返青期總體為提前趨勢,但受到地理區(qū)位的影響,東南部區(qū)域表現(xiàn)為提前,而西部返青期延后;相對于3—4月份累計降水量,返青期的變化更易受到該時間段平均氣溫的影響,溫度越高越有利于返青期的提前。

    關(guān)鍵詞:返青期;遙感監(jiān)測;NDVI

    中圖分類號:S127 文獻標(biāo)識碼:A DOI 編碼:10.3969/j.issn.1006-6500.2016.03.015

    草地生態(tài)系統(tǒng)是我國陸地上分布最廣的生態(tài)系統(tǒng),也是極為豐富的草地資源[1]。從資源利用角度來講它屬于生產(chǎn)資料,是畜牧業(yè)的基礎(chǔ);作為地帶性植被,它具有不可替代的生態(tài)功能,對維持地球生命支持系統(tǒng)的穩(wěn)定具有重要作用。內(nèi)蒙古自治區(qū)位于我國北疆,草地面積占全國草地面積的21.7%,橫跨“三北”、靠近京津,是北方尤其是京津地區(qū)的天然生態(tài)屏障[2]。但是近年來由于氣候干旱、降雨減少等自然因素及超載過牧等人為因素,草地生產(chǎn)力不斷下降、草原生態(tài)環(huán)境惡化現(xiàn)象日益嚴(yán)重、超過2/3的草地出現(xiàn)不同程度退化,造成以草原畜牧業(yè)為生的牧民的收入降低,牧民與草原生態(tài)環(huán)境關(guān)系失衡,引發(fā)出一系列社會、經(jīng)濟、生態(tài)問題。

    草原返青期是指在溫度和水分條件達到適宜狀態(tài)時,牧草開始萌發(fā)、變綠、生長的時間。目視觀察法是最常用的監(jiān)測方法,該方法簡單易行,但觀測點數(shù)量較多時需耗費大量人力和物力;且在大尺度研究區(qū)域,受到交通、地形地貌等環(huán)境因素的限制,只能做到“點”的觀察和監(jiān)測。隨著遙感技術(shù)的發(fā)展,特別是遙感物候?qū)W的出現(xiàn)[3],其應(yīng)用也擴展到了植被物候監(jiān)測方面 [4-5],結(jié)合少量的觀測驗證點便可以用相對較低的費用,提供大區(qū)域評估數(shù)據(jù),且能快速、有效監(jiān)測植被返青期[6]。通過監(jiān)測草原返青狀況及其動態(tài)變化,不僅有助于了解草原植被返青動態(tài)規(guī)律、植被生長發(fā)育狀況,預(yù)測本年度天然草原植被長勢及產(chǎn)量,更能為草原管理部門合理安排季節(jié)畜牧業(yè)生產(chǎn)的實施提供科學(xué)依據(jù)。

    為確保內(nèi)蒙古草地生態(tài)環(huán)境穩(wěn)定,為了提供決策所需的重要參數(shù),開展草地返青期遙感監(jiān)測技術(shù)研究是十分必要和迫切的。本研究選取內(nèi)蒙古錫林浩特市草原區(qū)作為典型的研究區(qū),采用非對稱高斯濾波歸一化植被指數(shù)(NDVI)重構(gòu)法,以及動態(tài)閾值返青期提取方法,對該地區(qū)草原區(qū)的返青期進行監(jiān)測,以期為該地區(qū)草地畜牧生產(chǎn)提供決策支持。

    1 研究區(qū)概況

    錫林浩特草原位于內(nèi)蒙古高原中部,地理位置北緯43°02′~44°52′,東經(jīng)115°13′~117°06′,南北長208 km,東西長143 km。地勢南高北低,南部為低山丘陵,北部為平緩的波狀平原,平均海拔高度988.5 m。屬于中溫帶半干旱大陸性氣候。降水量294.9 mm,降水分布不均勻,且主要集中于7、8月份。一年中氣溫最高月份是7月,1月份氣溫最低,無霜期110 d。年平均氣溫為1~2 ℃,年均氣溫差距較大??衫貌輬雒娣e1.4×105 hm2,地跨草甸草原、典型草原和沙丘沙地草原,草原內(nèi)的主要植被類型為濕地、森林、草原,伴隨有沙地灌叢、沼澤、草甸與河漫灘等。以大針茅(Stipagrandis)、羊草(Leymuschinense)和克氏針茅(S. krylovii)為主的叢生禾—根莖禾草草原,重點分布于河流域中部高平原。

    2 研究方法

    2.1 數(shù)據(jù)來源及處理

    研究表明,歸一化植被指數(shù)(NDVI)與綠色植物長勢有很好的相關(guān)性,可反映其生長狀況,因此,可以使用NDVI來提取草原植被的返青期。研究所采用的遙感數(shù)據(jù)為EOS/Terra衛(wèi)星的MODIS產(chǎn)品數(shù)據(jù)中的16 d合成產(chǎn)品MOD13Q1,該產(chǎn)品使用MVC方法合成,空間分辨率為250 m,來源于NASA網(wǎng)站(下載地址:http://ladsweb.nascom.nasa.gov/data/search.html),時間為2007—2014年。應(yīng)用ArcGIS 10.0軟件對該數(shù)據(jù)進行格式轉(zhuǎn)換、坐標(biāo)轉(zhuǎn)換、圖像裁剪等預(yù)處理。

    2.2 返青期提取方法

    2.2.1 NDVI時間序列重構(gòu) 由于受到云、氣溶膠和傳感器精度變化等因素的影響[7],在進行遙感數(shù)據(jù)預(yù)處理前,常采用最大合成法(MVC)去除時間序列數(shù)據(jù)中的局部異常點,以降低噪聲[3]。但仍存在“突降”或“針尖狀突峰”等奇異值,影響到對時間序列分析的精度。消除這些奇異值對NDVI曲線的影響,需要對NDVI時間序列進行重構(gòu)(或濾波)然后再判斷返青期的起始日期[8]。非對稱性高斯函數(shù)(Asymmetric Gaussian,AG)擬合法是主要的數(shù)據(jù)重構(gòu)方法之一[8]。

    非對稱高斯函數(shù)擬合方法是一個由半局部到整體的擬合方法,首先獲得NDVI時間序列曲線中谷、峰值的原始數(shù)據(jù),然后使用非對稱性高斯函數(shù)將曲線左右兩部分擬合(局部擬合),最后通過局部擬合構(gòu)造整體擬合函數(shù)。其中,局部擬合公式為:

    式(1)與(2)中,g(t;a1,… ,a5)為高斯函數(shù);t為時間變量;a1 為對應(yīng)t的峰值或谷值的位置參數(shù);a2 為決定函數(shù)曲線右邊寬度的參數(shù);a3 為決定函數(shù)曲線右邊平度的參數(shù);a4為決定函數(shù)曲線左邊寬度的參數(shù);a5 為決定函數(shù)曲線左邊平度的參數(shù);c1為基線的線性參數(shù);c2為振幅的線性參數(shù)。

    若局部擬合函數(shù)上曲線凸起部分的效果不明顯,需將曲線分為左、右邊谷值間隔區(qū)和中部峰值間隔區(qū),通過局部擬合構(gòu)造整體擬合函數(shù)。整體擬合函數(shù)為

    (3)

    式(3)中[tL,tR] 為整個NDVI 時間序列數(shù)據(jù)變化區(qū)間;fL(t) 為區(qū)間[tL,tR]內(nèi)左邊谷值所對應(yīng)的局部函數(shù);fC(t) 為區(qū)間[tL,tR]內(nèi)中間峰值所對應(yīng)的局部函數(shù);fR(t) 為區(qū)間[tL,tR]內(nèi)右邊谷值所對應(yīng)的局部函數(shù);α(t)和β(t)為介于[0,1]之間的剪切系數(shù)。

    2.2.2 返青期信息提取 NDVI時間序列圖像經(jīng)過濾波后,便可以進行草原返青期信息提取。常用的提取方法有:動態(tài)閾值法與NDVI閾值法。(1)閾值法即根據(jù)預(yù)先定義的NDVI參考值(如NDVI>0.3)確定返青期[6]。該方法操作簡單,但由于不同區(qū)域和不同植被類型的NDVI值存在顯著差異,確定適合大尺度物候期提取的單一NDVI閾值是不可行的。(2)動態(tài)閾值法,其改進之處在于通過考慮NDVI季節(jié)變化幅度,實現(xiàn)在像元尺度上閾值的動態(tài)設(shè)定,從而在一定程度上消除了土壤背景值和植被類型的影響[9]。目前,該方法已經(jīng)應(yīng)用于全球很多區(qū)域,取得了較好的應(yīng)用結(jié)果。本研究根據(jù)一些學(xué)者的研究結(jié)果,將返青期開始期定為NDVI年振幅的20%[9]。

    3 結(jié)果與分析

    3.1 錫林郭勒草原返青期變化

    根據(jù)錫林浩特草原返青期的遙感監(jiān)測可以看出(圖1和2),在2007—2014年間,錫林浩特草原返青期變化情況較大:2007年錫林浩特草原大部分區(qū)域返青期主要集中在3月份,返青期區(qū)域分布于該地區(qū)的東南部;2008年則推遲到5月上下旬,較2007年推遲近1個月,返青區(qū)域呈東南部先返青、西北部后返青的狀態(tài);2009—2011年返青期情況比較相似,從3月上旬草地開始返青(位于東南部),一直持續(xù)到5月下旬;2012年草原返青主要集中在4月上旬到5月上旬之間;2013年草原返青期較2014年提前1個月,出現(xiàn)與其他年份不同的現(xiàn)象,西北部分提前返青;而2014年返青期推遲到5月上旬,幾乎整個中部和北部地區(qū)都推遲返青??傮w來說,該地區(qū)牧草返青的日期大約為4月初至5月初。

    3.2 錫林郭勒草原氣候變化分析

    降水量變化和溫度變化可以影響返青期的提前和延后。因此,筆者計算了3—4月份的累積降水量和平均溫度情況(圖3)。錫林浩特草原3—4月份平均氣溫均呈不顯著升高趨勢,年際波動相對較大(變異系數(shù)為40.5%,年平均氣溫7.2 ℃),最低值出現(xiàn)在2010年(2.4 ℃),最高值出現(xiàn)在2014年(11.6 ℃)。3—4月份累計降水量變化趨勢從2007年到2014年呈減少狀態(tài)(平均值為18.1 mm),但是年際波動較大(變異系數(shù)73.0%),降水量最大值出現(xiàn)在2009年(47.0 mm),最低值出現(xiàn)在2011年(3.6 mm)。

    4 討 論

    4.1 氣候變化與返青期變化的相關(guān)性

    錫林浩特草原返青期總體呈現(xiàn)出提前趨勢,這與一些學(xué)者的研究結(jié)果相似[10]。牧草返青既受水分條件的影響,也受熱量條件的影響,其中降水和氣溫影響因素尤為重要[11]。在當(dāng)年3—4月份溫度和同期降水、氣溫條件較好年份,牧草的返青期有所提前,而在雨水較少和氣溫較低的條件下,可能導(dǎo)致牧草的返青期推遲。根據(jù)返青期與3—4月份降水和氣溫情況的分析來看,影響返青期的因素較多,雖然氣溫和降水只能解釋一些年份的返青期提前和延后,但仍然表現(xiàn)出以下規(guī)律:當(dāng)氣溫和降水同步時,有利于返青期的提前,如3—4月份平均氣溫和累計降水較高可以解釋2007、2009、2013年返青期提前的現(xiàn)象,2010年返青期相對于前幾年延后;當(dāng)降水和溫度不同步時,溫度的增加有利于返青期提前,如2011和2014年返青期提前,2012年返青期延后(當(dāng)年氣溫較低)。春季氣溫的高低通過控制融雪和土壤解凍過程對植物返青期的早晚起到關(guān)鍵的作用[12]。但也出現(xiàn)反常的年份,如2008年,雨熱同期,但是返青期卻大大延后,可能與其他因素有關(guān),如放牧影響、上一年冬季的低溫天氣等。

    4.2 政策實施與影響返青期變化的因素

    草原是我國最大、資源極為豐富的生態(tài)系統(tǒng),但是由于長期超負荷的使用草場以及牧民保護意識較弱等因素,草場利用率逐漸降低,草原退化問題越來越嚴(yán)重。這不單涉及草原生態(tài)平衡,還限制了草原區(qū)經(jīng)濟的快速發(fā)展,使牧民收入增長緩慢。為解決草原退化問題,從2011年開始中央財政每年發(fā)放資金134億元,對草原生態(tài)問題相當(dāng)嚴(yán)重、草原退化現(xiàn)象極其惡劣、不適合放牧的草場實行禁止放牧的保護補助獎勵機制。在禁止放牧區(qū)域以外的草場,對通過審核達到規(guī)定載畜量后的草場實行草原生態(tài)保護補助獎勵體制。增加牧民補貼,加強牧民對草原的保護意識,使其認識到草原問題不僅與經(jīng)濟問題關(guān)聯(lián),更與生存環(huán)境問題息息相關(guān)。

    由圖2可看出,2011年草原區(qū)實施“草原保護補助機制”后,錫林郭勒草原返青期整體提早,主要集中于3月下旬以及4月上旬,全區(qū)基本在5月上旬就已經(jīng)全面返青。返青期的總體提早與草原區(qū)實施“草原保護補助獎勵機制”有一定的關(guān)系,因為在草原植被階段較為敏感的時候,實行禁牧、休牧等方法,秋季有利于植物存儲較多的營養(yǎng)物質(zhì),保證來年的返青,對草原植被的生長發(fā)育較為有利。

    4.3 本研究不足之處

    本研究的誤差來源主要來自數(shù)據(jù)、濾波方法和動態(tài)閾值3個方面。首先,使用的MODIS NDVI時間分辨率較粗(16 d),為此研究草原植被返青期監(jiān)測只能獲得旬的水平(±16 d);定位觀測的方法,返青期時間間隔可以做到每天記錄。改進遙感數(shù)據(jù)的時間分辨率會改變這種不足,如MODIS NDVI可以處理為8 d的數(shù)據(jù)。其次,因為一些像元的NDVI時間序列曲線起伏變化較大,由于濾波函數(shù)的擬合關(guān)系,引起了擬合值與實際值有較大的偏差,從而降低了擬合精度[13],可以應(yīng)用多次設(shè)定濾波參數(shù)的方法,提高函數(shù)的擬合精度。本研究采用非對稱高斯濾波(AG)歸一化植被指數(shù)(NDVI)重構(gòu)法,濾波函數(shù)可能是誤差來源之一,濾波函數(shù)的實質(zhì)是做時間域上的均值,對時間序列NDVI進行濾波時會消弱最大值、增加最小值。相對于其他濾波函數(shù),非對稱高斯濾波可以提供更好的時間域平滑功能,減少信息的損耗。再次,動態(tài)閾值也是誤差的來源之一,本研究采用國內(nèi)外常用的閾值20%,但是有研究發(fā)現(xiàn)不同的草原類型,動態(tài)閾值范圍不同(15%~30%)。但是由于不同時期的NDVI值由于牲畜采食或刈割等原因造成其變化差異較大(降低NDVI的峰值,濾波結(jié)果使提取的返青期提前),所以在人類干擾比較強的地區(qū)使用動態(tài)閾值,可能會誤判返青期,對于放牧壓力大的區(qū)域,會誤判為返青期提前,而草地長勢比較好的區(qū)域,會誤判為返青期延后。但是,隨著遙感物候?qū)W監(jiān)測手段的成熟,并結(jié)合實地調(diào)查,將會大大提升數(shù)據(jù)的準(zhǔn)確性,推動遙感物候?qū)W的應(yīng)用。

    5 結(jié) 論

    本研究選取內(nèi)蒙古錫林浩特市草原區(qū)作為典型的研究區(qū),采用非對稱高斯函數(shù)對2007—2014年該區(qū)NDVI時間序列濾波,采用動態(tài)閾值法對該地區(qū)返青期的變化和3—4月份平均氣溫和累計降水量進行分析,結(jié)論如下:

    (1)錫林浩特草原區(qū)返青期總體為提前趨勢,但地域差異明顯,表現(xiàn)為東南部區(qū)域提前,而西部返青期延后。

    (2)該區(qū)返青期主要受3—4月平均氣溫的影響,在氣溫恒定的情況下,3—4月份累計降水量的增加有利于返青期的提前。

    參考文獻:

    [1] 李博. 中國北方草地退化及其防治對策[J]. 中國農(nóng)業(yè)科學(xué), 1997, 30(6): 1-9.

    [2] 謝高地, 張釔鋰, 魯春霞, 等. 中國自然草地生態(tài)系統(tǒng)服務(wù)價值[J]. 自然資源學(xué)報, 2001, 16(1): 47-53.

    [3] 陳效逑, 王林海. 遙感物候?qū)W研究進展[J]. 地理科學(xué)進展, 2009(1): 33-40.

    [4] SCHWARTZ M, REED B, WHITE M. Assessing satellite-derived start-of-season measures in the conterminous USA[J]. International Journal of Climatology, 2002, 22(14): 1793-1805.

    [5] ZHANG X, FRIEDL M, SCHAAF C. Global vegetation phenology from moderate resolution imaging spectroradiometer (MODIS): evaluation of global patterns and comparison with in situ measurements[J]. Journal of Geophysical Research, 2006, 111(G4): G04017.

    [6] KANG S, RUNNING S W, LIM J H, et al. A regional phenology model for detecting onset of greenness in temperate mixed forests, Korea: an application of MODIS leaf area index[J]. Remote Sensing of Environmental Management, 2003, 86(2):32-42.

    [7] 張峰, 吳炳方. 利用時序植被指數(shù)監(jiān)測作物物候的方法研究[J]. 農(nóng)業(yè)工程學(xué)報, 2004, 20(1): 155-159.

    [8] BRADLEY B A, JACOB R W. A curve fitting procedure to derive inter -annual phonologies from time series of noisy satellite NDVI data[J]. Remote Sensing of Environment, 2007, 10(6):137-145.

    [9] JONSSON P, EKLUNDH L. Seasonality extraction by function fitting to time-series of satellite sensor data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(8): 1824-1832.

    [10] 顧潤源, 周偉燦, 白美蘭, 等. 氣候變化對內(nèi)蒙古草原典型植物物候的影響[J]. 生態(tài)學(xué)報, 2012, (3): 767-776.

    [11] 白美蘭, 郝潤全, 李喜倉,等. 草原區(qū)草本植物物候期對氣候變化的響應(yīng)[C]//第28屆中國氣象學(xué)會年會:S4應(yīng)對氣候變化,發(fā)展低碳經(jīng)濟. 廈門: 中國氣象學(xué)會, 2011.

    [12] 陳效逑, 李倞. 內(nèi)蒙古草原羊草物候與氣象因子的關(guān)系[J].生態(tài)學(xué)報, 2009, (10): 5280-5290.

    [13] 王宏, 李曉兵, 余弘婧. 基于NOAA/AVHRR NDVI監(jiān)測中國北方典型草原的生長季及變化[J]. 植物生態(tài)學(xué)報, 2006(3): 365-374.

    丝袜喷水一区| 日韩一卡2卡3卡4卡2021年| 国产精品欧美亚洲77777| 丰满少妇做爰视频| 亚洲专区中文字幕在线| 国产精品九九99| 免费看十八禁软件| 国产精品一区二区在线不卡| 男人爽女人下面视频在线观看| 日本91视频免费播放| 99re6热这里在线精品视频| 国产成人影院久久av| 大码成人一级视频| 免费在线观看视频国产中文字幕亚洲 | 在线精品无人区一区二区三| 成年人黄色毛片网站| 久久ye,这里只有精品| 女性被躁到高潮视频| 女人高潮潮喷娇喘18禁视频| 精品国产乱码久久久久久男人| 一区二区av电影网| 天天躁日日躁夜夜躁夜夜| 成人国产一区最新在线观看| 国产av精品麻豆| 国产国语露脸激情在线看| 国产99久久九九免费精品| 韩国高清视频一区二区三区| 久久久久精品人妻al黑| 一进一出抽搐动态| 欧美日韩国产mv在线观看视频| 男人舔女人的私密视频| 成人亚洲精品一区在线观看| 最新在线观看一区二区三区| 一本久久精品| 免费少妇av软件| 天堂中文最新版在线下载| 久久久久久亚洲精品国产蜜桃av| 免费日韩欧美在线观看| 国产淫语在线视频| 久久久久国内视频| 国产亚洲午夜精品一区二区久久| av免费在线观看网站| 巨乳人妻的诱惑在线观看| 男女之事视频高清在线观看| 亚洲精品国产一区二区精华液| 岛国在线观看网站| 两个人免费观看高清视频| 天堂中文最新版在线下载| 爱豆传媒免费全集在线观看| 丝袜脚勾引网站| 精品国产一区二区三区四区第35| 国产福利在线免费观看视频| 国产人伦9x9x在线观看| 黄色 视频免费看| 日本黄色日本黄色录像| 黄色片一级片一级黄色片| 久久国产亚洲av麻豆专区| 久久久久久久大尺度免费视频| 99国产精品免费福利视频| 如日韩欧美国产精品一区二区三区| 国产精品欧美亚洲77777| 国产精品免费视频内射| 一二三四社区在线视频社区8| 国产av又大| 50天的宝宝边吃奶边哭怎么回事| 自线自在国产av| 久久女婷五月综合色啪小说| 可以免费在线观看a视频的电影网站| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 精品少妇黑人巨大在线播放| 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 久久国产精品男人的天堂亚洲| 久久99一区二区三区| 一级片'在线观看视频| 色老头精品视频在线观看| 精品一区二区三卡| 黄片大片在线免费观看| 91精品伊人久久大香线蕉| 嫩草影视91久久| 老司机影院毛片| 国产成人免费观看mmmm| 国产精品免费视频内射| 黄色视频不卡| 成人国语在线视频| 久久ye,这里只有精品| 黄色毛片三级朝国网站| tube8黄色片| 国产成人精品久久二区二区91| 国产精品免费大片| 亚洲专区国产一区二区| 亚洲男人天堂网一区| 成人黄色视频免费在线看| 精品国内亚洲2022精品成人 | 精品卡一卡二卡四卡免费| 国产日韩欧美视频二区| 91成年电影在线观看| 免费不卡黄色视频| 欧美日韩av久久| 日韩欧美一区视频在线观看| 亚洲专区中文字幕在线| 国产一区二区三区在线臀色熟女 | 亚洲欧美激情在线| 色婷婷av一区二区三区视频| 精品一区二区三区四区五区乱码| 精品一区在线观看国产| 搡老乐熟女国产| 国产av国产精品国产| 五月天丁香电影| 悠悠久久av| 日本91视频免费播放| 日韩欧美一区二区三区在线观看 | videosex国产| 国产精品 国内视频| 波多野结衣一区麻豆| 可以免费在线观看a视频的电影网站| 制服人妻中文乱码| 天天躁日日躁夜夜躁夜夜| 婷婷成人精品国产| 在线观看免费视频网站a站| 久久精品人人爽人人爽视色| 桃红色精品国产亚洲av| 中文精品一卡2卡3卡4更新| 日日摸夜夜添夜夜添小说| 黄网站色视频无遮挡免费观看| 少妇精品久久久久久久| 岛国毛片在线播放| 一区二区三区精品91| 午夜久久久在线观看| 女人爽到高潮嗷嗷叫在线视频| 日本欧美视频一区| 麻豆av在线久日| 色视频在线一区二区三区| 久久人人97超碰香蕉20202| 狠狠婷婷综合久久久久久88av| 国产男女内射视频| 亚洲成人免费电影在线观看| 精品久久久久久久毛片微露脸 | 伊人久久大香线蕉亚洲五| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| 在线观看人妻少妇| 黄片大片在线免费观看| 一本大道久久a久久精品| 两个人免费观看高清视频| 制服诱惑二区| 女性被躁到高潮视频| 宅男免费午夜| av片东京热男人的天堂| 亚洲午夜精品一区,二区,三区| 久久中文看片网| 国产国语露脸激情在线看| 欧美 亚洲 国产 日韩一| 午夜免费观看性视频| 久久九九热精品免费| 91精品伊人久久大香线蕉| 久久国产精品男人的天堂亚洲| 又大又爽又粗| 国产亚洲精品一区二区www | av欧美777| 国产av又大| 人妻人人澡人人爽人人| 一区福利在线观看| 欧美精品一区二区免费开放| 正在播放国产对白刺激| 欧美激情久久久久久爽电影 | 日韩大码丰满熟妇| √禁漫天堂资源中文www| 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 国产在线视频一区二区| 别揉我奶头~嗯~啊~动态视频 | 精品人妻熟女毛片av久久网站| 欧美精品人与动牲交sv欧美| 亚洲国产中文字幕在线视频| 曰老女人黄片| 男女午夜视频在线观看| 欧美日本中文国产一区发布| 桃红色精品国产亚洲av| 丰满人妻熟妇乱又伦精品不卡| a 毛片基地| 国产在线观看jvid| 精品人妻1区二区| 成人免费观看视频高清| 99热国产这里只有精品6| 啦啦啦在线免费观看视频4| 亚洲精品国产av蜜桃| 亚洲成av片中文字幕在线观看| 欧美日韩一级在线毛片| 黄色毛片三级朝国网站| 日本五十路高清| 香蕉丝袜av| 在线av久久热| 热re99久久国产66热| 18禁观看日本| 女性生殖器流出的白浆| 国产一区二区三区av在线| 夫妻午夜视频| 黑人巨大精品欧美一区二区mp4| 亚洲视频免费观看视频| 中亚洲国语对白在线视频| 亚洲综合色网址| 狂野欧美激情性bbbbbb| 欧美激情 高清一区二区三区| 欧美日韩视频精品一区| 看免费av毛片| 十八禁高潮呻吟视频| 国产伦人伦偷精品视频| av网站在线播放免费| 亚洲av片天天在线观看| 国产精品麻豆人妻色哟哟久久| 高清黄色对白视频在线免费看| 精品一区二区三区四区五区乱码| 我要看黄色一级片免费的| av欧美777| 欧美日韩视频精品一区| 亚洲av成人一区二区三| 国产在线观看jvid| 伊人亚洲综合成人网| 久久ye,这里只有精品| 三上悠亚av全集在线观看| 久久国产精品大桥未久av| 日韩大片免费观看网站| 水蜜桃什么品种好| 国产欧美日韩一区二区三 | 最近中文字幕2019免费版| tocl精华| 欧美精品av麻豆av| 蜜桃国产av成人99| 久久中文字幕一级| 成人三级做爰电影| 亚洲专区国产一区二区| av天堂在线播放| 90打野战视频偷拍视频| 一本—道久久a久久精品蜜桃钙片| 欧美人与性动交α欧美软件| 国产又爽黄色视频| 纵有疾风起免费观看全集完整版| 国产高清videossex| 日韩免费高清中文字幕av| 国产xxxxx性猛交| 满18在线观看网站| 91大片在线观看| 国产精品一二三区在线看| 99国产精品免费福利视频| 国产极品粉嫩免费观看在线| 一本综合久久免费| 国产xxxxx性猛交| 亚洲精品自拍成人| 精品少妇久久久久久888优播| 亚洲专区中文字幕在线| 一个人免费在线观看的高清视频 | 亚洲精品成人av观看孕妇| 五月开心婷婷网| 国产激情久久老熟女| 国产精品久久久久久精品电影小说| 欧美老熟妇乱子伦牲交| 午夜精品久久久久久毛片777| 在线观看舔阴道视频| 好男人电影高清在线观看| 俄罗斯特黄特色一大片| 欧美一级毛片孕妇| 欧美成人午夜精品| 亚洲成人国产一区在线观看| 久久精品成人免费网站| 中国国产av一级| 超碰97精品在线观看| 国产一区二区激情短视频 | 人人妻,人人澡人人爽秒播| 久久精品国产亚洲av高清一级| 多毛熟女@视频| 久久久精品区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 手机成人av网站| 亚洲第一青青草原| 亚洲国产av新网站| 无限看片的www在线观看| 午夜免费观看性视频| 香蕉国产在线看| 王馨瑶露胸无遮挡在线观看| 精品视频人人做人人爽| 99精品久久久久人妻精品| 亚洲精品成人av观看孕妇| 亚洲七黄色美女视频| 亚洲国产欧美网| 国产一区二区三区综合在线观看| 亚洲情色 制服丝袜| 精品视频人人做人人爽| 午夜激情久久久久久久| 欧美在线黄色| 国产男人的电影天堂91| 久久国产精品男人的天堂亚洲| 菩萨蛮人人尽说江南好唐韦庄| 亚洲熟女毛片儿| 国产免费现黄频在线看| 亚洲精品国产一区二区精华液| 亚洲一区中文字幕在线| 国产伦理片在线播放av一区| 国产极品粉嫩免费观看在线| 免费在线观看视频国产中文字幕亚洲 | 国产成人精品久久二区二区91| 亚洲国产精品一区三区| 久久人人97超碰香蕉20202| 久久精品国产亚洲av高清一级| 成人18禁高潮啪啪吃奶动态图| 丝瓜视频免费看黄片| 国产成人精品久久二区二区91| 亚洲情色 制服丝袜| 亚洲伊人色综图| 国产成人精品无人区| 国产成人av激情在线播放| 国产亚洲精品第一综合不卡| 黄片大片在线免费观看| 高清欧美精品videossex| 亚洲国产日韩一区二区| 亚洲精品乱久久久久久| 男人舔女人的私密视频| 人人妻人人澡人人看| 精品国产乱码久久久久久小说| 亚洲欧美日韩另类电影网站| 中亚洲国语对白在线视频| 在线 av 中文字幕| 高清黄色对白视频在线免费看| 美国免费a级毛片| 日韩熟女老妇一区二区性免费视频| 久久久久久久精品精品| 国产在线视频一区二区| 啦啦啦视频在线资源免费观看| a级片在线免费高清观看视频| 一区二区av电影网| 纯流量卡能插随身wifi吗| 中国美女看黄片| 亚洲国产精品成人久久小说| 一级片免费观看大全| a级毛片黄视频| 欧美激情高清一区二区三区| 免费一级毛片在线播放高清视频 | 黄色片一级片一级黄色片| 欧美日韩黄片免| 在线观看免费日韩欧美大片| 亚洲精品在线美女| 久9热在线精品视频| 99热全是精品| 香蕉丝袜av| 人人澡人人妻人| 人人妻,人人澡人人爽秒播| 天天操日日干夜夜撸| 一级毛片精品| 欧美97在线视频| 午夜精品久久久久久毛片777| 欧美性长视频在线观看| 一二三四社区在线视频社区8| 成人国语在线视频| 不卡一级毛片| 51午夜福利影视在线观看| 亚洲国产精品一区三区| 91成人精品电影| 久久精品亚洲av国产电影网| 久久精品国产综合久久久| 另类亚洲欧美激情| 精品一区二区三区四区五区乱码| 精品少妇一区二区三区视频日本电影| 视频区欧美日本亚洲| 中文字幕最新亚洲高清| 欧美精品高潮呻吟av久久| 日韩有码中文字幕| 国产深夜福利视频在线观看| 国产区一区二久久| 国产xxxxx性猛交| 久久久国产一区二区| 日韩电影二区| www.熟女人妻精品国产| 久久中文看片网| 久久热在线av| 亚洲精品美女久久久久99蜜臀| 两个人看的免费小视频| 啦啦啦 在线观看视频| 亚洲精品av麻豆狂野| 亚洲精品国产精品久久久不卡| 国产av国产精品国产| 国产三级黄色录像| 男男h啪啪无遮挡| 久久精品人人爽人人爽视色| 午夜福利影视在线免费观看| 亚洲国产精品一区二区三区在线| 日韩 亚洲 欧美在线| 麻豆av在线久日| 麻豆乱淫一区二区| 亚洲人成电影观看| 女人被躁到高潮嗷嗷叫费观| 三级毛片av免费| 国产在线一区二区三区精| av片东京热男人的天堂| 日本欧美视频一区| 久久久久精品国产欧美久久久 | 大片电影免费在线观看免费| 少妇被粗大的猛进出69影院| 国产精品一区二区免费欧美 | 国产精品久久久久成人av| 亚洲av日韩在线播放| 国产精品久久久久久精品电影小说| 一级片免费观看大全| 黑人巨大精品欧美一区二区mp4| 欧美日韩成人在线一区二区| 久久久久久人人人人人| 国产一区二区在线观看av| 国产精品久久久久久精品古装| 午夜福利在线观看吧| 中亚洲国语对白在线视频| 欧美一级毛片孕妇| 老汉色∧v一级毛片| 一级片免费观看大全| 中文字幕精品免费在线观看视频| 日韩,欧美,国产一区二区三区| 欧美国产精品一级二级三级| 国产野战对白在线观看| 91国产中文字幕| 水蜜桃什么品种好| 日韩欧美免费精品| 91老司机精品| 欧美在线黄色| 亚洲人成电影免费在线| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 久久av网站| 日本猛色少妇xxxxx猛交久久| 99国产极品粉嫩在线观看| 成人黄色视频免费在线看| 久久久久国产精品人妻一区二区| 狠狠婷婷综合久久久久久88av| 国产成人a∨麻豆精品| 美女高潮喷水抽搐中文字幕| 久久久国产一区二区| 久久久国产欧美日韩av| 欧美日韩一级在线毛片| 性色av一级| 一级,二级,三级黄色视频| 国产精品熟女久久久久浪| 成在线人永久免费视频| 精品国内亚洲2022精品成人 | 欧美精品一区二区免费开放| 亚洲专区国产一区二区| 欧美午夜高清在线| 精品人妻熟女毛片av久久网站| 国产亚洲午夜精品一区二区久久| 午夜日韩欧美国产| av线在线观看网站| 波多野结衣一区麻豆| 97人妻天天添夜夜摸| 男女之事视频高清在线观看| 亚洲成人手机| 国产精品久久久人人做人人爽| 亚洲成国产人片在线观看| www.自偷自拍.com| 国产黄频视频在线观看| 亚洲色图综合在线观看| videos熟女内射| 操出白浆在线播放| 动漫黄色视频在线观看| 久久精品成人免费网站| 亚洲va日本ⅴa欧美va伊人久久 | 久久天躁狠狠躁夜夜2o2o| 色老头精品视频在线观看| 最新的欧美精品一区二区| 飞空精品影院首页| 爱豆传媒免费全集在线观看| 国产一卡二卡三卡精品| 免费在线观看日本一区| 免费观看人在逋| 97精品久久久久久久久久精品| 亚洲熟女精品中文字幕| 亚洲成人免费av在线播放| 一级毛片电影观看| 亚洲中文av在线| 免费日韩欧美在线观看| 香蕉国产在线看| 首页视频小说图片口味搜索| 国产精品国产av在线观看| 精品少妇久久久久久888优播| 最近中文字幕2019免费版| 水蜜桃什么品种好| 操出白浆在线播放| 美女中出高潮动态图| 久久综合国产亚洲精品| 成在线人永久免费视频| 一级黄色大片毛片| 精品国产一区二区久久| 久久久精品免费免费高清| 50天的宝宝边吃奶边哭怎么回事| 香蕉丝袜av| 国产亚洲av高清不卡| 99香蕉大伊视频| 国产成人系列免费观看| 成人18禁高潮啪啪吃奶动态图| 国产成人精品在线电影| 国产精品免费视频内射| 人人妻人人澡人人看| 青春草视频在线免费观看| 黑人猛操日本美女一级片| 欧美黑人欧美精品刺激| 国产在视频线精品| 免费人妻精品一区二区三区视频| av天堂在线播放| 欧美少妇被猛烈插入视频| 亚洲国产日韩一区二区| 精品国产国语对白av| 亚洲av欧美aⅴ国产| 12—13女人毛片做爰片一| 中文字幕色久视频| 热99久久久久精品小说推荐| 窝窝影院91人妻| 国产亚洲精品久久久久5区| 欧美人与性动交α欧美软件| 亚洲一码二码三码区别大吗| 啦啦啦啦在线视频资源| 亚洲成国产人片在线观看| 欧美97在线视频| www.精华液| 欧美老熟妇乱子伦牲交| 国产在视频线精品| 久久九九热精品免费| 亚洲成人免费av在线播放| 国产日韩一区二区三区精品不卡| 国产成人欧美在线观看 | av天堂在线播放| 伦理电影免费视频| 黑丝袜美女国产一区| 亚洲精品第二区| 极品人妻少妇av视频| 国产欧美日韩一区二区三 | 日韩三级视频一区二区三区| a级毛片黄视频| 两性夫妻黄色片| 纯流量卡能插随身wifi吗| 久久久精品国产亚洲av高清涩受| 啦啦啦啦在线视频资源| 精品亚洲成a人片在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲美女黄色视频免费看| 久久人人爽人人片av| 波多野结衣一区麻豆| 97在线人人人人妻| 高清av免费在线| 国产欧美日韩一区二区三 | 激情视频va一区二区三区| 欧美成狂野欧美在线观看| 成年av动漫网址| 首页视频小说图片口味搜索| 欧美97在线视频| 欧美久久黑人一区二区| 岛国毛片在线播放| 欧美人与性动交α欧美软件| 亚洲激情五月婷婷啪啪| 国产高清视频在线播放一区 | 国产日韩欧美亚洲二区| 首页视频小说图片口味搜索| 国产精品国产av在线观看| 熟女少妇亚洲综合色aaa.| 久久毛片免费看一区二区三区| 欧美久久黑人一区二区| 国产在线免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 另类亚洲欧美激情| 精品少妇久久久久久888优播| 日韩视频一区二区在线观看| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 妹子高潮喷水视频| 精品熟女少妇八av免费久了| 国产野战对白在线观看| 高清av免费在线| 18禁国产床啪视频网站| 国产1区2区3区精品| 成年美女黄网站色视频大全免费| 久热爱精品视频在线9| 别揉我奶头~嗯~啊~动态视频 | 老司机福利观看| 久久中文看片网| 午夜福利视频精品| 美女脱内裤让男人舔精品视频| 99国产精品一区二区三区| 免费高清在线观看视频在线观看| 女人被躁到高潮嗷嗷叫费观| e午夜精品久久久久久久| 国产成人系列免费观看| 在线av久久热| 欧美乱码精品一区二区三区| 国产有黄有色有爽视频| 咕卡用的链子| 亚洲av欧美aⅴ国产| 亚洲五月色婷婷综合| 别揉我奶头~嗯~啊~动态视频 | 国产成人精品无人区| 国产视频一区二区在线看| 成人亚洲精品一区在线观看| 首页视频小说图片口味搜索| 欧美在线黄色| 丰满迷人的少妇在线观看| 大香蕉久久网| 天天躁日日躁夜夜躁夜夜| 电影成人av| 欧美日韩黄片免| 夜夜夜夜夜久久久久| avwww免费| 黄色怎么调成土黄色| 国产精品av久久久久免费| 精品福利观看| 美女高潮到喷水免费观看| 国产一区二区在线观看av| 国产又爽黄色视频| 9色porny在线观看| 91av网站免费观看| 9热在线视频观看99|