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

    基于被動微波遙感技術(shù)的玉米冠層葉面積指數(shù)反演

    2013-09-26 02:25:38馬紅章劉素美朱曉波孫根云柳欽火
    自然資源遙感 2013年3期
    關(guān)鍵詞:發(fā)射率冠層樣地

    馬紅章,劉素美,朱曉波,孫根云,孫 林,柳欽火

    (1.中國石油大學(xué)(華東)理學(xué)院,青島 266580;2.中國資源衛(wèi)星應(yīng)用中心,北京 100094;3.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,青島 266580;4.山東科技大學(xué)測繪科學(xué)學(xué)院,青島 266590;5.中國科學(xué)院遙感應(yīng)用研究所,北京 100101)

    0 引言

    植被在自然界物質(zhì)能量交換和全球氣候變化中扮演著重要角色[1]。葉面積指數(shù)(leaf area index,LAI)是陸地生態(tài)系統(tǒng)中重要的特征參量,控制著植被的生物物理過程[2];植被的蒸騰作用、光合作用以及地表凈初級生產(chǎn)力均與LAI密切相關(guān)[3-4]。研究表明,LAI可以很好地反映農(nóng)作物的生長狀況,在農(nóng)作物估產(chǎn)上具有重要應(yīng)用價值[5-6]。利用遙感數(shù)據(jù)反演大尺度LAI對研究陸地生態(tài)系統(tǒng)以及全球氣候變化具有重要意義。目前,大尺度LAI衛(wèi)星產(chǎn)品數(shù)據(jù)大多采用光學(xué)多光譜傳感器數(shù)據(jù)計算得到[7-9],但光學(xué)遙感受大氣和光照條件的影響較大,且對LAI的敏感性隨著LAI的增大而下降。被動微波遙感具有較大的觀測視場和全天時、全天候的工作能力,且微波對介質(zhì)具有一定的穿透性,能提供植被冠層的垂直廓線信息。目前,基于被動微波遙感技術(shù)的LAI反演研究還較少。為此,本文基于Matrix-Doubling(MD)模型對不同生長期玉米冠層微波輻射特性數(shù)據(jù)的模擬,研究了多波段被動微波亮溫數(shù)據(jù)在玉米冠層LAI反演方面的可行性與應(yīng)用潛力。

    1 試驗(yàn)概況與數(shù)據(jù)源

    1.1 試驗(yàn)概況

    圖1 試驗(yàn)區(qū)及玉米樣地位置示意圖Fig.1 Position of test site and sampling points

    試驗(yàn)區(qū)(圖1)位于河北省懷來縣官廳水庫的南岸,以實(shí)驗(yàn)站駐地(N40°20’56”,E115°47’04”)為中心,方圓近百km的區(qū)域。2010年5—10月,中國科學(xué)院遙感應(yīng)用研究所和北京師范大學(xué)等多家單位在本試驗(yàn)區(qū)開展了全波段主、被動聯(lián)合觀測遙感試驗(yàn),獲得了大量數(shù)據(jù),建立了不同生長期玉米冠層參數(shù)測量數(shù)據(jù)集。本研究應(yīng)用了數(shù)據(jù)集中不同生長期玉米結(jié)構(gòu)參數(shù)測量數(shù)據(jù)、LAI測量數(shù)據(jù)以及多玉米樣地4波段雙極化的微波輻射亮溫觀測數(shù)據(jù)。

    1.2 試驗(yàn)數(shù)據(jù)

    在懷來試驗(yàn)區(qū)1.5 km尺度上共測量了9塊玉米地,樣地編號為M1—M9(圖1),除M2樣地外,均采用平行壟向的觀測方式。對M2樣地采用平行壟和與壟成45°兩種方位測量,分別用編號M2-1和M2-2區(qū)分。冠層微波輻射亮溫數(shù)據(jù)由4頻段8通道車載微波輻射計測量得到。輻射計參數(shù)如表1所示;觀測平臺如圖2所示。

    表1 車載微波輻射計的基本參數(shù)Tab.1 Basic parameters of microwave radiometer

    圖2 觀測儀器與觀測平臺Fig.2 Observation instrument and observation platform

    為消除樣地玉米長勢不均對測量數(shù)據(jù)的影響,每個樣地均進(jìn)行多位置點(diǎn)的測量,取平均值作為該樣地觀測值。

    與樣地配套的觀測參數(shù)包括土壤0~3 cm深度內(nèi)的平均溫度(針式溫度計測量,精度優(yōu)于0.5 K)、觀測視場地表平均溫度(熱紅外熱像儀測量,精度優(yōu)于0.2 K)、植被含水量(烘干法測量,精度優(yōu)于1%)及植被葉片面積(收獲法LAI3000測量,精度優(yōu)于10 cm2)等等。在微波輻射觀測試驗(yàn)中,因?yàn)槊總€樣地測量過程持續(xù)時間短(不超過15 min),并且在多云天氣下進(jìn)行的,因此,玉米冠層溫度與土壤溫度較為接近。所以,在本研究中,輻射計觀測視場的地表微波等效溫度用視場平均溫度和土壤溫度的均值代替,各參數(shù)如表2所示。

    表2 樣地配套參數(shù)Tab.2 Sampling point supporting parameters

    2 研究方法

    2.1 研究思路

    如圖3所示,以不同生長期玉米冠層參數(shù)測量數(shù)據(jù)集作為MD微波輻射傳輸模型的輸入文件,通過回歸方法得到不同微波波段及不同極化方式下的玉米冠層微波發(fā)射率與透射率之間的關(guān)系,在此基礎(chǔ)上構(gòu)建參數(shù)化的微波輻射模型;應(yīng)用LAI和土壤水分參數(shù)驅(qū)動微波輻射傳輸模型,通過迭代計算玉米場景微波輻射模擬亮溫值與觀測亮溫值的匹配絕對誤差值,尋找最佳匹配的驅(qū)動參數(shù)值,使M值達(dá)到最小,實(shí)現(xiàn)對LAI的反演。

    圖3 LAI反演技術(shù)路線Fig.3 Technical route for LAI inversion

    2.2 模型參數(shù)的輸入

    在本試驗(yàn)期間共砍伐不同生長期玉米近100株。測量了株高、桿高、桿直徑、葉片面積、葉片及桿含水量等參數(shù);通過對玉米苗期到抽穗期的冠層結(jié)構(gòu)參數(shù)測量數(shù)據(jù)的回歸分析,得到冠層LAI與株高、株高與桿高、桿高與桿直徑的統(tǒng)計關(guān)系。另外,通過對大量實(shí)測數(shù)據(jù)的分析發(fā)現(xiàn),在非極端天氣條件下,不同生長期的玉米的含水量比葉片含水量平均高出12.5%,玉米葉片的含水量基本在70% ~80%范圍內(nèi)變化,觀測樣地土壤含水量在4% ~10%范圍內(nèi)變化。為了刻畫與觀測樣地更接近的玉米植被場景,研究中把以上觀測數(shù)據(jù)作為先驗(yàn)知識,微波輻射傳輸模型的輸入?yún)?shù)文件設(shè)置如表3所示。

    表3 輸入?yún)?shù)設(shè)置Fig.3 Set of model parameters

    2.3 模型及模擬數(shù)據(jù)庫

    研究表明,對C波段以上頻率的微波,植被內(nèi)部的多次散射以及植被層與土壤邊界之間的多次散射對研究相關(guān)問題變得越來越重要[10-11]。植被層應(yīng)被看作離散的散射體集合,用電磁波理論在更寬的頻率范圍內(nèi)更準(zhǔn)確地描述其散射體的散射和衰減特性[12-13]?;诠饩€跟蹤原理的MD模型算法最早應(yīng)用于計算大氣中非相干多次散射,后被引入到植被覆蓋地表的輻射計算中。在MD模型中[14]考慮了植被層內(nèi)部的多次散射,把植被冠層分割成N個無限薄子層,穿過任一薄子層的散射認(rèn)為單次的;同時,為計算任意方向上的散射強(qiáng)度,在薄子層入射和散射方向上用高斯-勒讓德(Gaussian-Legendre)方法分別劃分出M個方向,對于某個入射方向來說,相鄰2個薄子層Δτ1和Δτ2之間的多次散射過程如圖4所示。

    圖4 多次散射示意圖Fig.4 Diagram of multiple scattering

    前向透過矩陣S和后向散射矩陣T的計算式為

    通過反復(fù)迭代計算即可得到整個植被層任意方向上的散射矩陣,進(jìn)而可以獲取植被整層透過率與發(fā)射率。

    利用MD模型模擬玉米出苗到抽穗期的微波輻射特征。選擇玉米冠層場景5 000個,在每個場景下模擬觀測角度為55°時 4 個波段(6.925 GHz,10.7 GHz,18.7 GHz和 36.5GHz)雙極化(H,V 極化)的植被層發(fā)射率和透過率。

    2.4 模型參數(shù)化

    分析模擬數(shù)據(jù)發(fā)現(xiàn),在相同極化條件下,冠層發(fā)射率與透過率在4個波段上均保持良好線性關(guān)系,圖5給出了10.65 GHz和18.7 GHz波段在V和H極化方式時的冠層發(fā)射率與透過率間的關(guān)系。

    圖5 冠層發(fā)射率與透過率關(guān)系Fig.5 Relations between canopy emissivity and transmissivity

    其他波段植被冠層透過率與發(fā)射率之間的線性關(guān)系模型為

    式中:γp與分別為p極化條件下植被冠層微波透過率與發(fā)射率;下腳標(biāo)p表示極化方式(V極化或H極化);a與b為模型系數(shù)(表4)。

    表4 冠層微波透過率與發(fā)射率模型系數(shù)Tab.4 Model coefficient between transmittance and emissivity

    另外,6.925 GHz,18.7 GHz和36.5GHz波段的冠層發(fā)射率均與10.7 GHz的同極化冠層發(fā)射率具有較高的一致性關(guān)系,10.7 GHz冠層H極化發(fā)射率與V極化發(fā)射率同樣具有較好的一致性關(guān)系,因此,其他3個波段的冠層發(fā)射率均可由10.7 GHz冠層V極化發(fā)射率來推算。圖6給出了6.925 GHz和18.7 GHz波段的冠層發(fā)射率與10.7 GHz波段冠層發(fā)射率間的關(guān)系。

    圖6 波段間冠層發(fā)射率的關(guān)系Fig.6 Relations of bands canopy emissivity

    由10.7 GHz冠層發(fā)射率計算其他波段冠層發(fā)射率的表達(dá)式為

    表5列出了不同波段冠層發(fā)射率轉(zhuǎn)化模型系數(shù)。

    表5 冠層發(fā)射率轉(zhuǎn)化模型系數(shù)Fig.5 Model coefficient for canopy emissivity transform

    由式(2)(3),可通過10.7 GHz垂直極化玉米冠層發(fā)射率計算 6.925 GHz,18.7 GHz和 36.5 GHz波段冠層的雙極化發(fā)射率、透過率以及10.7 GHz玉米冠層水平極化發(fā)射率。

    在已知土壤含水量和粗糙度參數(shù)的條件下,根據(jù)QP模型[15]可計算波段10.7 GHz的土壤發(fā)射率,而其他3波段土壤發(fā)射率可用10.7GHz的同極化發(fā)射率來表示。即

    表6 土壤發(fā)射率模型系數(shù)Tab.6 Coefficient of soil emissivity model

    3 LAI的反演

    根據(jù)0階微波輻射傳輸方程,玉米覆蓋地表的微波輻射亮溫可用式(6)表示,即

    式中:γp為植被冠層微波透過率;和分別為土壤和植被冠層發(fā)射率;Tv和Ts分別為植被冠層和土壤溫度。

    由式(1)—(5)和QP模型知,只要給定10.7 GHz波段V極化冠層發(fā)射率、土壤體積含水量以及土壤層和冠層溫度,則通過式(6)可以得到4個波段雙極化的微波輻射亮溫。通過調(diào)整10.7 GHz波段冠層V極化發(fā)射率,模擬得到對應(yīng)場景下的4頻段、雙極化的微波輻射亮溫;通過模擬亮溫與實(shí)測亮溫的比較,可確定10.7 GHz波段冠層V極化發(fā)射率的最佳值,而10.7 GHz波段植被層V極化發(fā)射率與植被LAI之間的關(guān)系可由玉米冠層微波輻射特性模擬數(shù)據(jù)庫得到,即

    本文采用 6.925 GHz,10.7 GHz,18.7 GHz 和36.5 GHz 4個波段的玉米樣地微波亮溫觀測數(shù)據(jù),利用循環(huán)迭代最小二乘法,搜索10.7 GHz波段V極化發(fā)射率與土壤水分的最佳值,再結(jié)合式(7)計算樣地LAI。通過計算得到的LAI與測量值相關(guān)性高達(dá)0.91,如圖7所示。

    圖7 LAI反演結(jié)果Fig.7 Inversion result of LAI

    反演值與實(shí)測值均方根誤差RMSE=0.173 1。由此可見,多波段被動微波輻射亮溫數(shù)據(jù)在農(nóng)作物冠層LAI反演方面具有一定可行性和應(yīng)用潛力。

    4 結(jié)論與討論

    本文利用玉米結(jié)構(gòu)參數(shù)實(shí)測數(shù)據(jù)生成Matrix-Doubling模型的輸入文件,通過建立玉米植被冠層多波段雙極化輻射特性數(shù)據(jù)庫,回歸分析得到了玉米冠層各波段之間的發(fā)射率關(guān)系以及各波段發(fā)射率與透過率的關(guān)系。通過微波輻射傳輸參數(shù)化模型,利用循環(huán)迭代算法對河北懷來10個玉米樣地進(jìn)行了LAI的反演,LAI反演值與實(shí)測LAI相關(guān)性大于0.9,表明多波段被動微波數(shù)據(jù)在玉米LAI反演方面具有一定的應(yīng)用潛力。

    LAI反演誤差產(chǎn)生原因是植被覆蓋地表微波輻射亮溫模擬和地表參數(shù)反演中的幾個突出問題,具體分析如下:

    1)本文利用部分玉米結(jié)構(gòu)參數(shù)測量數(shù)據(jù)回歸得到模型輸入文件,是對“一般”種植和生長狀態(tài)玉米的概括,具有局地適用性。

    2)反演中沒有考慮大氣微波下行輻射及背景微波輻射的影響,是造成LAI反演值與實(shí)測值間具有一定偏差的原因之一。

    3)在0階微波輻射傳輸模型的應(yīng)用中,為考慮植被冠層對高頻微波的散射影響,文中以半球各向到觀測方向的透過率之和作為植被層觀測方向的透過率,雖然這樣能部分彌補(bǔ)植被下行輻射被地表散射到觀測方向部分以及非觀測方向土壤微波發(fā)射被植被層散射到觀測方向的部分,但是,土壤表面的極化差異以及各向的非朗伯特性等因素對微波輻射亮溫的影響仍然不可忽略。

    另外,隨著微波頻率的增加,植被層對微波信號衰減增強(qiáng),本文中應(yīng)用的是C波段6.925 GHz以上頻率的微波,所以地表微波亮溫主要為植被層信息,而土壤信號較弱,導(dǎo)致玉米樣地土壤水分的反演結(jié)果(文中并沒有給出相應(yīng)反演結(jié)果)不佳。對有植被覆蓋地表而言,由于不同頻段微波對冠層的穿透能力不同,低頻微波對土壤水分信息會更加敏感,所以在植被覆蓋地表的土壤水分反演中,高低頻被動微波的聯(lián)合應(yīng)用會更具優(yōu)勢。

    [1]Kuusk A.Monitoring of vegetation parameters on large areas by the inversion of a canopy reflectance model[J].International Journal of Remote Sensing,1998,19(15):2893-2905.

    [2]Chen J M,Cihlar J.Retrieving leaf area index of boreal conifer forests using landsat TM images[J].Remote Sensing of Environment,1996,55(43):153-162.

    [3]Bonan G.Importance of leaf area index and forest type when estimating photosynthesis in boreal forests[J].Remote Sensing of Environment,1993,43(3):303-314.

    [4]Pierce L L,Running S W.Rapid estimation of coniferous forest leaf area index using a portable integrating radiometer[J].Ecology,1988,69(18):1762-1767.

    [5]Holben B N,Tucer C J.Spectral assessment of soybean leaf area and leaf biomass[J].Photogramnetric Engineering and Remote Sensing,1980,46(5):651-656.

    [6]Rasmussen M S.Operational yield forecast using AVHRR NDVI data:Reduction of environmental and inter- annual variability[J].International Journal of Remote Sensing,1997,18(5):1059-1077.

    [7]Myneni R B,Nemani R R,Running S W.Estimation of global leaf area index and absorbed PAR using radiative transfer models[J].IEEE Transactions on Geoscience and Remote Sensing,1997,35(66):1380-1393.

    [8]Verhoef W.Light scattering by leaf layers with application to canopy reflectance modeling,the SAIL model[J].Remote Sensing of Environment,1984,16(2):125-141.

    [9]Li X W,Strahler A H.Geometric-optical model of a conifer forest canopy[J].IEEE Transactions on Geoscience and Remote Sensing,1985,23(5):705-720.

    [10]Njoku J.Soil moisture retrieval from AMSR- E[J].IEEE Transactions on Geoscience and Remote Sensing,2003,41(2):215-229.

    [11]Ulaby F T,Moore R K,F(xiàn)ung A K.Microwave remote sensing:Active and passive[M].Dedham,MA:Artech House,1986.

    [12]Wigneron J P,Calvet J C,Kerr Y.Microwave emission of vegetation:A sensitivity to leaf characteristics[J].IEEE Trans Geosci Remote Sensing,1993,23(31):716-726.

    [13]Press W H,Teukolsky S A.Numerical recipes in C:The art of scientific computing[M]2nd edition.London:Cambridge University Press,1992.

    [14]Karam M A.A physical model for microwave radiometry of vegetation[J].IEEE Transactions on Geoscience and Remote Sensing,1997,35(4):1045-1058.

    [15]Chen K S,Jackson T J,Oneill P E.A parameterized surface reflectivity model and estimation of bare surface soil moisture with L-band radiometer[J].IEEE Trans Geosci Remote Sens,2002,40(12):2674-2686.

    [16]Chen K S,Wu T D,Tsang L.Emission of rough surfaces calculated by the integral equation method with a comparison to a three-dimensional moment method simulations[J].IEEE Trans Geosci Remote Sens,2003,41(1):90-101.

    猜你喜歡
    發(fā)射率冠層樣地
    基于低空遙感的果樹冠層信息提取方法研究
    基于激光雷達(dá)的樹形靶標(biāo)冠層葉面積探測模型研究
    額爾古納市興安落葉松中齡林植被碳儲量研究
    安徽省淮南森林冠層輻射傳輸過程的特征
    昆明市主要綠化樹種閾值測定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    氧氣A(O,O)波段氣輝體發(fā)射率和臨邊輻射強(qiáng)度模擬與分析
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    低溫狀態(tài)下的材料法向發(fā)射率測量
    塔克拉瑪干沙漠地表發(fā)射率及分布變化特征
    色哟哟哟哟哟哟| 丝袜美腿在线中文| 国产av麻豆久久久久久久| 爱豆传媒免费全集在线观看| 亚洲图色成人| 成人三级黄色视频| 少妇熟女欧美另类| 日韩av不卡免费在线播放| 麻豆成人午夜福利视频| av.在线天堂| 97热精品久久久久久| 国产老妇女一区| 一边摸一边抽搐一进一小说| 变态另类成人亚洲欧美熟女| videossex国产| 高清午夜精品一区二区三区 | 欧美一级a爱片免费观看看| 国产精品免费一区二区三区在线| 尤物成人国产欧美一区二区三区| 热99在线观看视频| 别揉我奶头 嗯啊视频| 国产伦精品一区二区三区视频9| 一级av片app| 最新中文字幕久久久久| 高清毛片免费看| 亚洲综合色惰| 午夜激情福利司机影院| 亚洲精品乱码久久久久久按摩| 狂野欧美白嫩少妇大欣赏| 亚洲激情五月婷婷啪啪| a级毛色黄片| 国产 一区精品| 大香蕉久久网| 国产白丝娇喘喷水9色精品| 精品久久久久久久久久久久久| 老司机影院成人| 免费看美女性在线毛片视频| 欧美精品国产亚洲| 波多野结衣巨乳人妻| 日本爱情动作片www.在线观看| 久久久久久久久中文| 草草在线视频免费看| 国产黄色视频一区二区在线观看 | 欧美日本视频| av在线播放精品| 欧美日韩一区二区视频在线观看视频在线 | 国产成人a区在线观看| 中文字幕精品亚洲无线码一区| 69av精品久久久久久| 免费看美女性在线毛片视频| 熟女电影av网| 韩国av在线不卡| 草草在线视频免费看| av专区在线播放| 天美传媒精品一区二区| 精品不卡国产一区二区三区| 日韩亚洲欧美综合| 国产精品无大码| 国产综合懂色| 色综合站精品国产| 亚洲在久久综合| 一本久久中文字幕| 91午夜精品亚洲一区二区三区| 免费看日本二区| 99久久中文字幕三级久久日本| 久久欧美精品欧美久久欧美| 欧美极品一区二区三区四区| 国产精品免费一区二区三区在线| 18+在线观看网站| 久久韩国三级中文字幕| 亚洲一区二区三区色噜噜| 成年免费大片在线观看| 久久亚洲精品不卡| 欧美变态另类bdsm刘玥| 91久久精品国产一区二区成人| 国产女主播在线喷水免费视频网站 | 成人av在线播放网站| 国产不卡一卡二| 亚洲丝袜综合中文字幕| 两个人视频免费观看高清| 国产高潮美女av| 中文字幕制服av| 99在线人妻在线中文字幕| 不卡视频在线观看欧美| 亚洲av免费高清在线观看| 97人妻精品一区二区三区麻豆| 真实男女啪啪啪动态图| 偷拍熟女少妇极品色| 中文亚洲av片在线观看爽| 一级av片app| 国产黄a三级三级三级人| 少妇熟女aⅴ在线视频| 精品一区二区三区视频在线| 色5月婷婷丁香| 一边摸一边抽搐一进一小说| 日本与韩国留学比较| 国产高清三级在线| 免费观看人在逋| 亚洲中文字幕日韩| 性欧美人与动物交配| 国产高清视频在线观看网站| 欧美人与善性xxx| 1024手机看黄色片| 在线观看一区二区三区| 一级毛片aaaaaa免费看小| 亚洲国产精品成人综合色| 国产三级中文精品| av专区在线播放| 岛国毛片在线播放| 久久久久久久久久久免费av| 国产欧美日韩精品一区二区| 亚洲aⅴ乱码一区二区在线播放| 国内少妇人妻偷人精品xxx网站| 国产精品久久电影中文字幕| 我要看日韩黄色一级片| 日本-黄色视频高清免费观看| 99久国产av精品| a级毛片免费高清观看在线播放| 在线观看免费视频日本深夜| 欧美3d第一页| 亚洲国产欧美人成| 在线天堂最新版资源| 精品日产1卡2卡| 国产一区二区三区在线臀色熟女| 夜夜看夜夜爽夜夜摸| 人妻久久中文字幕网| 青春草亚洲视频在线观看| 亚洲一区二区三区色噜噜| 99久久成人亚洲精品观看| 日韩av在线大香蕉| 国产精品一区www在线观看| 男插女下体视频免费在线播放| 午夜福利在线观看吧| 国产高清有码在线观看视频| 深爱激情五月婷婷| 97人妻精品一区二区三区麻豆| 精品少妇黑人巨大在线播放 | 亚洲精品影视一区二区三区av| 色视频www国产| 国产黄a三级三级三级人| 亚洲一区二区三区色噜噜| 日本一本二区三区精品| 欧美成人一区二区免费高清观看| 最新中文字幕久久久久| 狂野欧美激情性xxxx在线观看| 亚洲无线在线观看| 国产一区二区亚洲精品在线观看| 国产成年人精品一区二区| 亚洲综合色惰| 国产成人91sexporn| 久久久久久国产a免费观看| 欧美+亚洲+日韩+国产| 国产精品不卡视频一区二区| 在线免费观看不下载黄p国产| 国产精品三级大全| 久久久午夜欧美精品| 亚洲国产日韩欧美精品在线观看| 欧美+亚洲+日韩+国产| 狂野欧美白嫩少妇大欣赏| 91久久精品国产一区二区三区| 神马国产精品三级电影在线观看| 麻豆一二三区av精品| 51国产日韩欧美| 91久久精品国产一区二区三区| 又粗又爽又猛毛片免费看| 亚洲丝袜综合中文字幕| 狂野欧美白嫩少妇大欣赏| 91久久精品国产一区二区三区| 欧美又色又爽又黄视频| 国产成人午夜福利电影在线观看| 最近的中文字幕免费完整| 激情 狠狠 欧美| 激情 狠狠 欧美| 亚洲真实伦在线观看| 尾随美女入室| 免费看日本二区| 亚洲aⅴ乱码一区二区在线播放| 午夜a级毛片| 国产熟女欧美一区二区| 老师上课跳d突然被开到最大视频| 国模一区二区三区四区视频| 美女黄网站色视频| 成年女人看的毛片在线观看| 久久久久久久亚洲中文字幕| 国产亚洲av嫩草精品影院| 床上黄色一级片| 日韩人妻高清精品专区| 亚洲精华国产精华液的使用体验 | 网址你懂的国产日韩在线| 久久精品国产亚洲av天美| 国模一区二区三区四区视频| 亚洲精品色激情综合| 成人毛片a级毛片在线播放| 看黄色毛片网站| 国产69精品久久久久777片| 搡老妇女老女人老熟妇| 99久久精品热视频| 国产免费一级a男人的天堂| 国产激情偷乱视频一区二区| 国产视频内射| 亚洲av熟女| 久久久久久伊人网av| 91久久精品电影网| 国产伦理片在线播放av一区 | 大又大粗又爽又黄少妇毛片口| 国产中年淑女户外野战色| 免费观看人在逋| 丰满的人妻完整版| 一进一出抽搐gif免费好疼| 国产精品不卡视频一区二区| 中文字幕免费在线视频6| 国内精品宾馆在线| 嫩草影院精品99| 男人舔女人下体高潮全视频| 欧美日韩一区二区视频在线观看视频在线 | 国产精品国产高清国产av| 真实男女啪啪啪动态图| 精品无人区乱码1区二区| 99国产精品一区二区蜜桃av| 高清日韩中文字幕在线| 99久久久亚洲精品蜜臀av| 久久久久久久久久成人| 一本精品99久久精品77| 51国产日韩欧美| 99久久无色码亚洲精品果冻| 国产高清视频在线观看网站| 最好的美女福利视频网| av天堂中文字幕网| av专区在线播放| 成人亚洲精品av一区二区| 好男人在线观看高清免费视频| 少妇裸体淫交视频免费看高清| 丰满的人妻完整版| 亚洲精品乱码久久久久久按摩| 免费无遮挡裸体视频| 一区二区三区高清视频在线| 国产高清激情床上av| av又黄又爽大尺度在线免费看 | 一进一出抽搐gif免费好疼| 亚洲欧美日韩无卡精品| 国产精品一区二区性色av| 麻豆成人av视频| 欧美一级a爱片免费观看看| 午夜福利成人在线免费观看| 精品一区二区三区视频在线| 日日干狠狠操夜夜爽| 又粗又硬又长又爽又黄的视频 | 国产男人的电影天堂91| 成人国产麻豆网| 国产亚洲5aaaaa淫片| 亚洲欧美精品自产自拍| 日韩人妻高清精品专区| 狠狠狠狠99中文字幕| 国产91av在线免费观看| 日本成人三级电影网站| 国产男人的电影天堂91| 久99久视频精品免费| 欧美日本视频| 久久国产乱子免费精品| 能在线免费看毛片的网站| 又爽又黄a免费视频| 我的老师免费观看完整版| 乱人视频在线观看| 男人狂女人下面高潮的视频| 白带黄色成豆腐渣| 99热网站在线观看| 好男人在线观看高清免费视频| 久久久a久久爽久久v久久| 久久精品91蜜桃| 欧美一区二区国产精品久久精品| 热99re8久久精品国产| 看片在线看免费视频| 三级国产精品欧美在线观看| 国产精品av视频在线免费观看| 一级毛片我不卡| 亚洲国产精品国产精品| 日本成人三级电影网站| 国产精品一区二区三区四区久久| 好男人视频免费观看在线| 国产成人a区在线观看| ponron亚洲| av免费在线看不卡| 精品久久久久久久末码| 十八禁国产超污无遮挡网站| 22中文网久久字幕| 97在线视频观看| 老司机影院成人| 最近的中文字幕免费完整| 99久久人妻综合| .国产精品久久| 丰满人妻一区二区三区视频av| 日韩高清综合在线| 国产成人a区在线观看| 男的添女的下面高潮视频| 噜噜噜噜噜久久久久久91| 青春草国产在线视频 | 中文欧美无线码| av在线播放精品| 欧美成人免费av一区二区三区| 久久亚洲国产成人精品v| 自拍偷自拍亚洲精品老妇| 日本爱情动作片www.在线观看| 国产精品伦人一区二区| 小说图片视频综合网站| 亚洲精品色激情综合| 久久久久性生活片| 赤兔流量卡办理| avwww免费| 18禁裸乳无遮挡免费网站照片| 精品久久久噜噜| 2022亚洲国产成人精品| 国产黄色视频一区二区在线观看 | 国产一区二区三区av在线 | 26uuu在线亚洲综合色| 免费在线观看成人毛片| 日产精品乱码卡一卡2卡三| 国产精品一区www在线观看| 国产极品精品免费视频能看的| 毛片一级片免费看久久久久| 欧美人与善性xxx| 中国国产av一级| 波多野结衣巨乳人妻| 长腿黑丝高跟| 一级毛片aaaaaa免费看小| 亚洲18禁久久av| 在线播放国产精品三级| 99久久精品热视频| 国产在视频线在精品| 国产不卡一卡二| 亚洲精品456在线播放app| 最好的美女福利视频网| 色视频www国产| 亚洲精品久久久久久婷婷小说 | 91精品国产九色| 只有这里有精品99| 搡老妇女老女人老熟妇| 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| 成人一区二区视频在线观看| 男女边吃奶边做爰视频| 看黄色毛片网站| 色综合站精品国产| 欧美高清成人免费视频www| 天堂√8在线中文| 99热这里只有是精品50| 国产成人精品婷婷| 国产伦精品一区二区三区四那| 久久精品影院6| 日韩视频在线欧美| 又粗又硬又长又爽又黄的视频 | 国产蜜桃级精品一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 黄片无遮挡物在线观看| 欧美丝袜亚洲另类| 精品人妻偷拍中文字幕| 国产精品美女特级片免费视频播放器| 国产爱豆传媒在线观看| 秋霞在线观看毛片| 亚洲成人精品中文字幕电影| 欧美成人精品欧美一级黄| 少妇人妻一区二区三区视频| 成熟少妇高潮喷水视频| 亚洲美女视频黄频| 国产大屁股一区二区在线视频| 日日撸夜夜添| 午夜精品一区二区三区免费看| 国产私拍福利视频在线观看| 99热这里只有是精品50| 九色成人免费人妻av| 国产精品电影一区二区三区| 综合色av麻豆| 日韩成人伦理影院| 在线观看美女被高潮喷水网站| 亚洲国产日韩欧美精品在线观看| 国产黄色小视频在线观看| 不卡视频在线观看欧美| 男女视频在线观看网站免费| 插逼视频在线观看| 麻豆国产av国片精品| 天堂中文最新版在线下载 | 欧美最新免费一区二区三区| 最近手机中文字幕大全| 亚洲一区二区三区色噜噜| 精品欧美国产一区二区三| 天堂中文最新版在线下载 | 人体艺术视频欧美日本| 亚洲精品国产av成人精品| 日本-黄色视频高清免费观看| 国内精品宾馆在线| 日韩一本色道免费dvd| 在线免费观看不下载黄p国产| 午夜福利在线在线| 蜜桃亚洲精品一区二区三区| 久久精品国产亚洲av香蕉五月| 日本三级黄在线观看| 99热网站在线观看| 国产探花在线观看一区二区| 亚洲最大成人av| av在线亚洲专区| 亚洲av一区综合| 久久久久国产网址| 欧美成人a在线观看| 成人亚洲欧美一区二区av| 午夜免费激情av| 精品一区二区三区人妻视频| 午夜激情福利司机影院| 日本色播在线视频| 尤物成人国产欧美一区二区三区| 欧美精品国产亚洲| 97超碰精品成人国产| 国产成人精品久久久久久| 午夜免费激情av| 色哟哟哟哟哟哟| 乱系列少妇在线播放| 日日啪夜夜撸| 中文字幕制服av| 日韩制服骚丝袜av| 国产精品麻豆人妻色哟哟久久 | 3wmmmm亚洲av在线观看| 级片在线观看| 有码 亚洲区| 欧美变态另类bdsm刘玥| 日本免费a在线| 白带黄色成豆腐渣| 国内精品宾馆在线| 日日啪夜夜撸| 欧美另类亚洲清纯唯美| 99久久中文字幕三级久久日本| 亚洲在线观看片| 亚洲精品色激情综合| 中文字幕久久专区| 日韩强制内射视频| 国产精品国产高清国产av| 中国美女看黄片| 你懂的网址亚洲精品在线观看 | 免费人成在线观看视频色| 综合色av麻豆| 国产成人freesex在线| 亚洲高清免费不卡视频| 1000部很黄的大片| 麻豆av噜噜一区二区三区| 美女大奶头视频| 国产精品永久免费网站| 亚洲国产欧洲综合997久久,| 国产精品不卡视频一区二区| 一个人观看的视频www高清免费观看| 人体艺术视频欧美日本| 国产成人aa在线观看| 日本与韩国留学比较| 精品99又大又爽又粗少妇毛片| 91麻豆精品激情在线观看国产| 啦啦啦观看免费观看视频高清| a级毛片a级免费在线| 免费看光身美女| 久久韩国三级中文字幕| 色哟哟哟哟哟哟| 人妻系列 视频| 老司机福利观看| 亚洲真实伦在线观看| 午夜精品国产一区二区电影 | 国产久久久一区二区三区| 美女黄网站色视频| 午夜免费男女啪啪视频观看| 69av精品久久久久久| 国产精品久久久久久久电影| 国产中年淑女户外野战色| 免费观看a级毛片全部| 国产在视频线在精品| 国产黄a三级三级三级人| 欧美日韩国产亚洲二区| 天天躁日日操中文字幕| 日本黄大片高清| 美女 人体艺术 gogo| 3wmmmm亚洲av在线观看| 久久久精品大字幕| 我要看日韩黄色一级片| .国产精品久久| 在线免费十八禁| 2021天堂中文幕一二区在线观| 成熟少妇高潮喷水视频| 国产男人的电影天堂91| 国产精品99久久久久久久久| 一本久久精品| 亚洲精品乱码久久久久久按摩| 特大巨黑吊av在线直播| 亚洲国产欧美在线一区| 欧美人与善性xxx| 一级毛片我不卡| 欧美+亚洲+日韩+国产| 中文字幕制服av| 精品午夜福利在线看| 久久亚洲国产成人精品v| 亚洲高清免费不卡视频| 男女啪啪激烈高潮av片| 国产精品蜜桃在线观看 | 亚洲精品乱码久久久v下载方式| 亚洲熟妇中文字幕五十中出| 国产免费男女视频| 国产亚洲精品久久久com| av天堂在线播放| 99热这里只有精品一区| 又粗又硬又长又爽又黄的视频 | 色综合站精品国产| 成年版毛片免费区| 成人av在线播放网站| av国产免费在线观看| 久久久a久久爽久久v久久| 中文字幕熟女人妻在线| 国产黄片美女视频| 成人性生交大片免费视频hd| 欧美bdsm另类| 亚洲丝袜综合中文字幕| 91久久精品国产一区二区成人| 1000部很黄的大片| 亚洲av不卡在线观看| 亚洲av成人av| 国产精品久久久久久精品电影小说 | 亚洲第一区二区三区不卡| 中文字幕精品亚洲无线码一区| 国产一区二区在线观看日韩| 联通29元200g的流量卡| 国产黄片视频在线免费观看| 国产精品一区二区性色av| 久久精品国产99精品国产亚洲性色| 国产精品国产三级国产av玫瑰| 中文资源天堂在线| 男插女下体视频免费在线播放| 国产三级在线视频| 一本精品99久久精品77| 99九九线精品视频在线观看视频| 全区人妻精品视频| 亚洲一区高清亚洲精品| 亚洲自偷自拍三级| 亚洲最大成人中文| 老司机影院成人| 欧美日韩精品成人综合77777| 日韩大尺度精品在线看网址| 校园人妻丝袜中文字幕| 人妻夜夜爽99麻豆av| 亚洲av不卡在线观看| 亚洲精品粉嫩美女一区| 国产黄片美女视频| 免费在线观看成人毛片| 国产精品免费一区二区三区在线| 伦精品一区二区三区| 国内精品宾馆在线| 欧美日韩精品成人综合77777| 欧美最黄视频在线播放免费| 亚洲aⅴ乱码一区二区在线播放| 91久久精品电影网| 日本三级黄在线观看| or卡值多少钱| 亚洲国产高清在线一区二区三| 国产视频内射| 亚洲自偷自拍三级| 黄色配什么色好看| 国产真实乱freesex| 看非洲黑人一级黄片| 国产精品电影一区二区三区| av卡一久久| 99国产极品粉嫩在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲自拍偷在线| 亚洲在线观看片| 久久6这里有精品| 中文字幕av成人在线电影| 在线播放国产精品三级| 免费大片18禁| 悠悠久久av| 高清毛片免费看| 青春草亚洲视频在线观看| 亚洲第一区二区三区不卡| 特大巨黑吊av在线直播| 久久精品综合一区二区三区| 欧美区成人在线视频| 成人二区视频| 有码 亚洲区| 一级黄色大片毛片| 女人被狂操c到高潮| 哪里可以看免费的av片| 少妇高潮的动态图| 99久久精品国产国产毛片| 精品一区二区三区人妻视频| 久久午夜亚洲精品久久| 成人性生交大片免费视频hd| 国产女主播在线喷水免费视频网站 | 非洲黑人性xxxx精品又粗又长| 自拍偷自拍亚洲精品老妇| 国产单亲对白刺激| 亚洲精品成人久久久久久| 国产大屁股一区二区在线视频| 欧美性猛交黑人性爽| 色综合色国产| 欧美成人一区二区免费高清观看| 一个人看视频在线观看www免费| 色播亚洲综合网| 天天一区二区日本电影三级| 男插女下体视频免费在线播放| 男的添女的下面高潮视频| 性欧美人与动物交配| av免费观看日本| 一本精品99久久精品77| 麻豆国产97在线/欧美| 国产精品久久电影中文字幕| 国产伦精品一区二区三区视频9| 欧美极品一区二区三区四区| 丰满的人妻完整版| 欧美最黄视频在线播放免费| 欧美日韩乱码在线| 欧美三级亚洲精品| 美女高潮的动态| 观看免费一级毛片| 丰满乱子伦码专区| 老熟妇乱子伦视频在线观看| 亚洲va在线va天堂va国产|