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

    基于遙感的海南島橡膠林零平面位移和粗糙度估算

    2016-05-30 08:42:29劉少軍易雪張京紅蔡大鑫田光輝李偉光
    熱帶作物學(xué)報(bào) 2016年10期
    關(guān)鍵詞:橡膠林海南島粗糙度

    劉少軍 易雪 張京紅 蔡大鑫 田光輝 李偉光

    摘 要 在橡膠風(fēng)害評估中,橡膠林零平面位移和粗糙度是建立定量評估模型的重要?jiǎng)恿W(xué)參數(shù)。以MODIS遙感數(shù)據(jù)為信息源,計(jì)算橡膠林的葉面積指數(shù),并利用形態(tài)學(xué)Raupach方法,估算橡膠林歸一化零平面位移和粗糙度。該方法充分發(fā)揮了遙感大面積獲取數(shù)據(jù)的優(yōu)勢,同時(shí)解決了傳統(tǒng)氣象學(xué)方法獲取橡膠林零平面位移和粗糙度需要大量長期的野外觀測資料的限制。結(jié)果表明:橡膠林葉面積指數(shù)與歸一化零平面位移和粗糙度分別存在正相關(guān)和負(fù)相關(guān);估算結(jié)果精度能滿足橡膠風(fēng)害評估模型參數(shù)的需求。研究結(jié)果可為開展大范圍內(nèi)橡膠林風(fēng)害評估提供技術(shù)支撐。

    關(guān)鍵詞 橡膠林;零平面位移;粗糙度;海南島

    中圖分類號 P49 文獻(xiàn)標(biāo)識碼 A

    Abstract Zero-plane displacement and roughness length of rubber are crucial parameters in the typhoon impact assessment model on rubber plant. Leaf area index was derived based on the MODIS remote sensing data. Zero-plane displacement and roughness length of rubber were derived by the morphology model. This method not only made full use of the advantage of remote sensing, but also solved the traditional meteorology method to get zero-plane displacement and roughness depending on a long-term field observation data. The results showed that leaf area indices were positively and negatively correlated with d/h and z0/h. The accuracy of result could meet the need of rubber wind damage assessment model parameters. The results provided technical support for carrying out large scale rubber wind damage assessment.

    Key words Rubber plant;Zero-Plane Displacement;Roughness Length;Hainan Island

    doi 10.3969/j.issn.1000-2561.2016.10.029

    零平面位移和粗糙度是描述下墊面空氣動力學(xué)特征的重要參數(shù)[1]。零平面位移為下墊面開始吸收動量的高度,地表粗糙度是下墊面風(fēng)速為零的高度[2],但在不同天氣背景下, 零平面位移和粗糙度長度并不是常數(shù)[3]。在海南橡膠林臺風(fēng)災(zāi)害動力學(xué)評估模型中,橡膠樹斷、倒條件的判識,必須解決大范圍、不同季節(jié)、不同地形條件下橡膠林零平面位移和粗糙度。由于早期的研究對零平面位移和粗糙度的理解不確切,一般通過冠層高度和葉面積指數(shù)等分別建立零平面位移和粗糙度的經(jīng)驗(yàn)公式進(jìn)行估算[4-7]。零平面位移和粗糙度是風(fēng)速和大氣層結(jié)的函數(shù),而不僅僅是由冠層幾何高度來確定[3,8]。目前,零平面位移和粗糙度的計(jì)算一般存在氣象學(xué)和形態(tài)學(xué)兩種方法[9]。氣象學(xué)方法常采用三維超聲風(fēng)溫儀,觀測多層的風(fēng)速觀測資料,然后利用風(fēng)速廓線公式通過迭代求解[1,10]。如張學(xué)仕等[1]利用應(yīng)用溫度方差方法計(jì)算了江蘇句容下蜀次生櫟林空氣動力學(xué)參數(shù);鐘中等[3]根據(jù)空氣質(zhì)量守恒原理, 對長白山天然闊葉紅松林冠層空氣動力學(xué)參數(shù)進(jìn)行了數(shù)值計(jì)算;趙曉松等[10]利用廓線法中的牛頓迭代法計(jì)算了長白山闊葉紅松林的零平面位移和粗糙度;Lo[11]根據(jù)冠層上下多層風(fēng)、溫觀測資料計(jì)算高大森林冠層零平面位移和粗糙度的方法。而形體學(xué)方法主要是采用粗糙元高度、粗糙元頂部面積、粗糙元迎風(fēng)面積等計(jì)算零平面位移和粗糙度[12]。如楊阿強(qiáng)等[9]采用粗糙元形態(tài)學(xué)方法,計(jì)算中國東部不同季節(jié)空氣動力學(xué)零平面位移和粗糙度;Macdonald等[13]根據(jù)粗糙元迎風(fēng)面積和頂部面積的方法計(jì)算零平面位移和粗糙度。

    形態(tài)學(xué)方法比氣象學(xué)方法更具有優(yōu)勢。氣象學(xué)方法物理意義更為確切,不同參數(shù)具有不同的物理含義,但計(jì)算過程中需要在穩(wěn)定中性大氣條件下準(zhǔn)確測量氣象數(shù)據(jù),推導(dǎo)橡膠林零平面位移和粗糙度對于大范圍長時(shí)間序列的計(jì)算需要大量長期的野外觀測,成本較高;而形態(tài)學(xué)方法是根據(jù)物理基礎(chǔ)經(jīng)過簡化綜合經(jīng)驗(yàn)系數(shù)得到的,精度需要不斷提高,但其方法可以采用遙感數(shù)據(jù)作為數(shù)據(jù)源,而且可以進(jìn)行大范圍內(nèi)參數(shù)的反演[9]。目前,關(guān)于采用遙感和形態(tài)學(xué)方法估算橡膠林零平面位移和粗糙度的研究未見報(bào)道。因此,本研究采用遙感數(shù)據(jù)源,嘗試?yán)眯螒B(tài)學(xué)模型的方法,估算海南島橡膠林零平面位移和粗糙度,為快速大范圍獲取橡膠林動力參數(shù)提供了新的技術(shù)途徑。

    1 材料與方法

    1.1 材料

    海南島MODIS遙感數(shù)據(jù)、TM遙感數(shù)據(jù)提取海南天然橡膠空間分布圖(具體方法參考文獻(xiàn)[14-15])、2011年海南天然橡膠普查數(shù)據(jù)、橡膠生產(chǎn)分布以農(nóng)場為單位進(jìn)行統(tǒng)計(jì),數(shù)據(jù)來源于文獻(xiàn)[16]。省界和縣界行政區(qū)劃圖來源于國家基礎(chǔ)地理信息網(wǎng)站提供的1 ∶ 400萬基礎(chǔ)地理信息數(shù)據(jù)(http://ngcc.sbsm.gov.cn/)。

    1.2 方法

    在形態(tài)學(xué)方法中,一般認(rèn)為Raupach[17]方法精度較高[18],因此選擇Raupach方法進(jìn)行海南島橡膠林零平面位移和粗糙度估算。通過提取橡膠林葉面積指數(shù)、相關(guān)參數(shù)計(jì)算,代入Raupach 的形態(tài)學(xué)模型,計(jì)算海南島橡膠林零平面位移和粗糙度。具體技術(shù)流程見圖1。

    1.2.1 橡膠葉面積指數(shù)提取 葉面積指數(shù)(Leaf Area Index,LAI)又叫葉面積系數(shù),是指單位土地面積上植物葉片總面積占土地面積的倍數(shù)。即:葉面積指數(shù)=葉片總面積/土地面積。葉面積指數(shù)(LAI)作為植被冠層的重要結(jié)構(gòu)參數(shù)之一,是許多生態(tài)系統(tǒng)生產(chǎn)力模型和全球氣候、水文、生物地球化學(xué)和生態(tài)學(xué)模型的關(guān)鍵輸入?yún)?shù)。葉面積指數(shù)的計(jì)算方法有:直接方法(葉面積的測定、描形稱重法、儀器測定法);間接方法(點(diǎn)接觸法、消光系數(shù)法、經(jīng)驗(yàn)公式法、遙感方法)。本文通過橡膠林葉面積指數(shù)(LAI)和歸一化植被指數(shù)(NDVI)統(tǒng)計(jì)關(guān)系,計(jì)算得到如下回歸方程,進(jìn)行天然橡膠葉面積指數(shù)的反演。

    LAI=7.82NDVI+0.781 (1)

    其中,歸一化植被指數(shù)的計(jì)算公式如下:

    式中,NDVI表示歸一化植被指數(shù),Rnir表示近紅外波段反照率,Rred表示可見光波段(紅光)反照率。

    1.2.2 零平面位移和粗糙度估算 Raupach的形態(tài)學(xué)公式表示如下[17]:

    d為零平面位移,ZO為粗糙度,h為粗糙元高度,λ為單位地面粗糙元迎風(fēng)面積,ψh為粗糙層影響函數(shù),b為冠層寬度,Cd為經(jīng)驗(yàn)系數(shù)。β為粗糙元阻力系數(shù)與地表阻力系數(shù)之比,γ為距地面高度z 處的水平風(fēng)速與摩擦速度之比。其中單位地面粗糙元迎風(fēng)面積λ可以為:

    A表示為單位面積上植被對風(fēng)的阻擋面積,值與葉面積指數(shù)很接近[19],因此采用遙感反演的葉面積指數(shù)LAI表示單位面積上植被對風(fēng)的阻擋面積。

    模型中其他參數(shù)參考相關(guān)文獻(xiàn)獲取[4-8,10-13]。

    2 結(jié)果與分析

    本次研究選擇了2015年3月中旬云量相對較少的MODIS遙感數(shù)據(jù),并結(jié)合海南島橡膠分布[14-16],分別提取了海南島橡膠林的葉面積指數(shù)、歸一化零平面位移、歸一化粗糙度分布圖。

    結(jié)果顯示:海南島橡膠林的葉面積指數(shù)在1.0~8.6之間,高值區(qū)分布在海南島的中部和東南部區(qū)域,低值區(qū)域分布在海南島的北部和西部區(qū)域(圖2),遙感葉面積指數(shù)的差異體現(xiàn)了此時(shí)橡膠冠層物候的差異。其主要原因是該時(shí)段海南島橡膠普遍進(jìn)入第一蓬葉展葉期,但由于不同區(qū)域的溫度、降水、光照、風(fēng)速等的差異,對此時(shí)橡膠葉的生長有著不同的影響,一般以東南和南部地區(qū)最早出現(xiàn),最后是西北部和中部地區(qū)[20]。

    海南橡膠林歸一化零平面位移在0.5~0.82之間,歸一化零平面位移最大值為0.82,最小值為0.5,平均值為0.75,標(biāo)準(zhǔn)差0.046。高值區(qū)分布在海南島的中部和東南部區(qū)域,低值區(qū)域分布在海南島的北部和西部區(qū)域(圖3)。由上可知,橡膠林的歸一化平均零平面位移與葉面積指數(shù)變化趨勢相一致,高值和低值區(qū)分布基本一致。通過ArcGIS軟件的Spatial analyst tools-Multivariate-Band collection statistics功能模塊,計(jì)算兩者的相關(guān)性為0.85。

    海南橡膠林歸一化粗糙度在0.056~0.10之間,歸一化粗糙度最大值為0.10,最小值為0.056,平均值為0.073,標(biāo)準(zhǔn)差0.012。低值區(qū)分布在海南島的中部區(qū)域,高值區(qū)域分布在海南島的北部和西南部區(qū)域(圖4)。歸一化粗糙度的變化趨勢與橡膠林葉面積指數(shù)高值和低值區(qū)的分布相反。通過ArcGIS軟件計(jì)算兩者的相關(guān)性為-0.84。其主要原因可能是當(dāng)橡膠林葉面積指數(shù)高值時(shí),橡膠樹葉片相對繁茂,使氣流從作物頂部通過,從而減小了氣流同葉片的摩擦降低了空氣動力學(xué)粗糙度。

    3 討論與結(jié)論

    本文在借鑒前人研究方法的基礎(chǔ)上,結(jié)合資料可獲得性及滿足大范圍獲取橡膠林零平面位移和粗糙度的需求,嘗試?yán)肕ODIS遙感數(shù)據(jù)和Raupach的形態(tài)學(xué)模型,采用的計(jì)算方法是建立在嚴(yán)格的理論推導(dǎo)基礎(chǔ)上,同時(shí)兼顧經(jīng)濟(jì)性原則,希望突破以往有關(guān)零平面位移和粗糙度取值的經(jīng)驗(yàn)性約束,為橡膠林動力學(xué)風(fēng)害評估模型提供動力參數(shù)。

    橡膠林零平面位移和粗糙度實(shí)際上是隨大氣穩(wěn)定情況發(fā)生改變的,可能同一天都會相差很大,不同大氣層結(jié)狀態(tài),零平面位移和粗糙度均會有所差異,同時(shí)橡膠林平面位移和粗糙度隨季節(jié)會發(fā)生變化。遙感反演的海南橡膠林d/h和z0/h與橡膠林葉面積指數(shù)分別存在正相關(guān)和負(fù)相關(guān)的關(guān)系。這與Shaw等[21]、趙曉松等[10]利用氣象觀測塔實(shí)測數(shù)據(jù)計(jì)算的規(guī)律一致。

    本研究由于缺乏與遙感數(shù)據(jù)相對應(yīng)時(shí)間段的通量塔觀測資料與之進(jìn)行對比驗(yàn)證,導(dǎo)致無法進(jìn)行遙感反演精度的直接驗(yàn)證。但橡膠樹屬于闊葉林,相關(guān)參數(shù)的驗(yàn)證僅與前人[9-10,22-24]的森林的結(jié)果進(jìn)行對比分析,根據(jù)2015年3月15日MODIS數(shù)據(jù)反演的海南島橡膠林歸一化零平面位移和粗糙度均在同一量級范圍內(nèi),反演結(jié)果整體上反映了橡膠林歸一化零平面位移和粗糙度在空間分布的差異,能滿足橡膠林風(fēng)害模型的參數(shù)要求,但具體的誤差仍需要與后期建立橡膠林觀測數(shù)據(jù)進(jìn)行對比分析,不斷改進(jìn)Raupach的形態(tài)學(xué)模型相關(guān)參數(shù),以提高遙感反演橡膠林零平面位移和粗糙度的精度。

    由于本文僅停留在利用Raupach方法,探索性地估算橡膠林零平面位移和粗糙度,而且相關(guān)參數(shù)主要參考了經(jīng)驗(yàn)值;同時(shí)由于橡膠林的面積的來源主要參考文獻(xiàn)[14-16]的成果,沒能及時(shí)更新,難免存在誤差和不足。

    參考文獻(xiàn)

    [1] 張學(xué)仕, 蔣 琰, 薛建輝, 等.江蘇句容下蜀次生櫟林的空氣動力學(xué)參數(shù)研究[J]. 南京林業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版, 2010, 34(6): 61-65.

    [2] 蔡 福, 周廣勝, 明惠青, 等. 玉米農(nóng)田空氣動力學(xué)參數(shù)動態(tài)及其與影響因子的關(guān)系[J]. 生態(tài)學(xué)報(bào), 2013, 33(17): 5 339-5 352.

    [3]鐘 中, 韓士杰. 長白山闊葉紅松林冠層空氣動力學(xué)參數(shù)的計(jì)算[J]. 南京大學(xué)學(xué)報(bào)(自然科學(xué)), 2002, 38(4): 565-571.

    [4] Szeicz G, Endrodi G, Tajchman S. Aerodynamic and surface factors in evaporation[J]. Water Resources Research, 1969, 5(2): 380-394 .

    [5] Stanhill G. A Simple instrument for field measurement of turbulent diffusion flux[J]. Journal of Applied Meteorology, 1969(8): 509-513 .

    [6] Zoumakis N M. Determination of the mean wind speed and momentum diffusivity profiles above tall vegetation and forest canopies using a mass conservation assumption[J]. Journal of Applied Meteorology, 1994, 33: 295-302.

    [7] Raupach M R. Simplified expressions forvegetation roughness length and zeroplane displacement as functions of canopy height and area index[J]. Boundary-Layer Meteorology, 1994, 71(1-2): 211-216.

    [8] Lesnik G E.Estimation of displacement layer thickness and roughness parameter from measurements above a forest[J]. Izvestiya Atmospheric And Oceanic Physics, 1973, 9(1): 88-90. [9] 楊阿強(qiáng), 孫國清, 盧立新, 等.基于MODIS資料的中國東部時(shí)間序列空氣動力學(xué)粗糙度和零平面位移高度估算[J]. 氣象科學(xué), 2011, 31(4): 516-524.

    [10] 趙曉松, 關(guān)德新, 吳家兵, 等.長白山闊葉紅松林的零平面位移和粗糙度[J]. 生態(tài)學(xué)雜志, 2004, 23(5): 84-88.

    [11] Lo A K. Determination of zero-plane displacement and roughness length of a forest canopy using profiles of limited height[J]. Boundary-Layer Meteorology, 1995, 75: 81-402.

    [12] RaupachM R, Antonia R A, Rajagopalan S. Rough-wall turbulent boundary layers[J]. Applied Mechanics Reviews, 1991, 44(1): 1-25.

    [13] Macdonald R W, Griffiths R F, Hall D J. An improved method for estimation of surface roughness of obstacle arrays[J]. Atmospheric Environment, 1998, 32(11): 1 857-1 864.

    [14] 張京紅, 陶忠良, 劉少軍, 等.基于TM影像的海南島橡膠種植面積信息提取[J]. 2010, 31(4): 661-665.

    [15]劉少軍, 張京紅, 何政偉. 基于面向?qū)ο蟮南鹉z分布面積估算研究[J]. 廣東農(nóng)業(yè)科學(xué), 2010, 37(1): 168-170.

    [16] 唐群鋒, 郭澎濤, 劉志崴, 等. 基于FMT-AHP的海南農(nóng)墾花崗巖類多雨區(qū)橡膠園地力評價(jià)[J]. 生態(tài)學(xué)報(bào), 2014, 34(15): 4 435-4 445.

    [17] Jasinski M F, Borak J, Crago R. Bulk surface momentum parameters for satellite-derived vegetation elds[J]. Agricultural and Forest Meteorology, 2005, 133: 55-68.

    [18] Grimmond C S B, Oke T R.Aerodynamic properties of urban areas derived from analysis of surface form[J]. Journal of Applied Meteorology, 1999, 38: 1 262-1 292.

    [19] ZENG Xubin, Shaikh Muhammad, Dai Yongjiu, et al. Coupling of the common land model to the NCAR community climate model[J]. Journal of Climate, 2001, 15: 1 832-1 854.

    [20] 陳小敏, 陳匯林, 李偉光, 等. 海南島天然橡膠林春季物候期的遙感監(jiān)測[J]. 中國農(nóng)業(yè)氣象, 2016, 37(1): 111-116.

    [21] Shaw R H, Pereira A R. Aerodynamic roughness of a plant canopy: a numerical experiment[J]. Agricultural Meteorology, 1982, 26: 51-65.

    [22] 連士華. 橡膠樹風(fēng)害成因間題的探討[J]. 熱帶作物學(xué)報(bào), 1984, 5(1): 59-72.

    [23] 劉和平, 劉樹華, 朱廷曜, 等.森林冠層空氣動力學(xué)參數(shù)的確定[J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 1997, 33(4): 522-528.

    [24] Borak Jordan S, Jasinski Michael F, Crago Richard D. Time seriesvegetation aerodynamic roughness fields estimated from MODISobservations[J]. Agricultural and Forest Meteorology, 2005, 135: 252-268.

    猜你喜歡
    橡膠林海南島粗糙度
    海南橡膠林生態(tài)系統(tǒng)凈碳交換物候特征
    橡膠林
    基于無人機(jī)影像的巖體結(jié)構(gòu)面粗糙度獲取
    甘肅科技(2020年20期)2020-04-13 00:30:18
    冷沖模磨削表面粗糙度的加工試驗(yàn)與應(yīng)用
    模具制造(2019年4期)2019-06-24 03:36:48
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    鋼材銹蝕率與表面三維粗糙度參數(shù)的關(guān)系
    西雙版納橡膠林土壤有機(jī)碳分布特征研究
    是海南省還是海南島?
    關(guān)于海南島戰(zhàn)役作戰(zhàn)方針的考察
    軍事歷史(1994年6期)1994-08-15 08:56:38
    回憶海南島渡海登陸戰(zhàn)役
    軍事歷史(1987年2期)1987-08-20 06:10:48
    久9热在线精品视频| av在线蜜桃| 老熟妇乱子伦视频在线观看| 久久久久免费精品人妻一区二区| 成人高潮视频无遮挡免费网站| 成人国产综合亚洲| 久久国产精品影院| 午夜福利在线观看免费完整高清在 | 亚洲国产精品合色在线| 国产野战对白在线观看| 日韩欧美免费精品| 国产精品久久久久久久久免 | 99热只有精品国产| 1024手机看黄色片| 欧美最黄视频在线播放免费| 身体一侧抽搐| 久久精品国产清高在天天线| 美女cb高潮喷水在线观看| 欧美国产日韩亚洲一区| 小说图片视频综合网站| 国产精品一及| 亚洲av免费高清在线观看| 欧美3d第一页| 在线十欧美十亚洲十日本专区| 免费看a级黄色片| 久久久久久大精品| 久久久久性生活片| 国产伦在线观看视频一区| 97碰自拍视频| 黄色视频,在线免费观看| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 国产精品久久久久久亚洲av鲁大| 久久精品影院6| 免费大片18禁| 欧美另类亚洲清纯唯美| 午夜两性在线视频| 日本一本二区三区精品| 给我免费播放毛片高清在线观看| 日韩有码中文字幕| 99久国产av精品| 亚洲美女搞黄在线观看 | 国产精品免费一区二区三区在线| 亚洲综合色惰| 亚洲内射少妇av| 欧美一区二区亚洲| 男女视频在线观看网站免费| 久久久国产成人免费| 十八禁国产超污无遮挡网站| 国产高潮美女av| 波多野结衣高清作品| 久久九九热精品免费| 国产色爽女视频免费观看| 国产精品98久久久久久宅男小说| 国产精品免费一区二区三区在线| 99在线视频只有这里精品首页| 成人美女网站在线观看视频| 亚洲 国产 在线| 色哟哟·www| 成人av在线播放网站| 啦啦啦观看免费观看视频高清| 我的老师免费观看完整版| 亚洲18禁久久av| 日本一二三区视频观看| 男女那种视频在线观看| 91午夜精品亚洲一区二区三区 | 日韩有码中文字幕| 日日摸夜夜添夜夜添av毛片 | 中文字幕熟女人妻在线| 99久国产av精品| 国产毛片a区久久久久| 精品久久久久久久久亚洲 | 老司机福利观看| 久久久久国产精品人妻aⅴ院| 搡老岳熟女国产| а√天堂www在线а√下载| 99久国产av精品| 成人国产综合亚洲| 天堂动漫精品| 热99re8久久精品国产| 人人妻人人看人人澡| 久久欧美精品欧美久久欧美| 嫩草影院精品99| 少妇人妻精品综合一区二区 | 国产av不卡久久| 1024手机看黄色片| 色综合婷婷激情| 欧美日本亚洲视频在线播放| 欧美激情在线99| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 免费大片18禁| 别揉我奶头 嗯啊视频| 亚洲专区国产一区二区| 亚洲精品色激情综合| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 亚洲美女黄片视频| 成人午夜高清在线视频| 国产精品一区二区免费欧美| 亚洲精品在线观看二区| 性色av乱码一区二区三区2| 亚洲成a人片在线一区二区| 欧美成人一区二区免费高清观看| 日韩欧美在线二视频| 日韩有码中文字幕| 高清在线国产一区| 搡老熟女国产l中国老女人| 一进一出抽搐gif免费好疼| 亚洲国产精品sss在线观看| 国产一区二区激情短视频| 精品久久久久久久人妻蜜臀av| 少妇熟女aⅴ在线视频| 天堂av国产一区二区熟女人妻| 久久久久久久精品吃奶| 欧美日韩乱码在线| 两个人视频免费观看高清| 最后的刺客免费高清国语| 老司机午夜十八禁免费视频| 淫妇啪啪啪对白视频| 亚洲国产精品合色在线| 日本五十路高清| 噜噜噜噜噜久久久久久91| 日韩亚洲欧美综合| 国产大屁股一区二区在线视频| 国产精品国产高清国产av| 午夜影院日韩av| АⅤ资源中文在线天堂| 97热精品久久久久久| 蜜桃亚洲精品一区二区三区| 天堂av国产一区二区熟女人妻| 欧美日韩黄片免| 黄色女人牲交| 免费大片18禁| 欧美乱妇无乱码| 国产又黄又爽又无遮挡在线| 大型黄色视频在线免费观看| 91在线精品国自产拍蜜月| 久久伊人香网站| 在线a可以看的网站| 波多野结衣高清作品| 非洲黑人性xxxx精品又粗又长| 国产探花极品一区二区| 日韩欧美一区二区三区在线观看| 午夜激情欧美在线| 又黄又爽又刺激的免费视频.| 中文资源天堂在线| 久久国产精品影院| 高潮久久久久久久久久久不卡| 国产av麻豆久久久久久久| 亚洲,欧美精品.| 露出奶头的视频| 一级作爱视频免费观看| a级一级毛片免费在线观看| 婷婷精品国产亚洲av在线| а√天堂www在线а√下载| 日韩欧美国产在线观看| 成人国产一区最新在线观看| 嫩草影院入口| 日日摸夜夜添夜夜添小说| 日本与韩国留学比较| 亚洲欧美日韩卡通动漫| 成人一区二区视频在线观看| 亚洲欧美成人综合另类久久久| 婷婷色综合www| 午夜视频国产福利| 菩萨蛮人人尽说江南好唐韦庄| 在现免费观看毛片| 五月开心婷婷网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 在线观看一区二区三区激情| 亚洲欧美日韩无卡精品| 国产淫片久久久久久久久| 尾随美女入室| 国产亚洲91精品色在线| 狠狠精品人妻久久久久久综合| 欧美xxⅹ黑人| 免费观看性生交大片5| tube8黄色片| 国产亚洲5aaaaa淫片| 国产免费视频播放在线视频| 熟女av电影| 亚洲欧洲国产日韩| 人妻系列 视频| 人人妻人人澡人人爽人人夜夜| 麻豆国产97在线/欧美| 免费看av在线观看网站| 午夜福利网站1000一区二区三区| 国产精品久久久久久av不卡| 91精品一卡2卡3卡4卡| 精品一区二区免费观看| 韩国av在线不卡| 国产精品.久久久| 精品久久久久久久久av| 久热久热在线精品观看| 国产永久视频网站| 国产免费一区二区三区四区乱码| 国产精品99久久久久久久久| 久久久久久国产a免费观看| 极品教师在线视频| 日韩一本色道免费dvd| 欧美zozozo另类| 日韩视频在线欧美| 深夜a级毛片| 99久久精品国产国产毛片| 亚洲精华国产精华液的使用体验| 亚洲av福利一区| 美女视频免费永久观看网站| 人妻夜夜爽99麻豆av| 精品少妇久久久久久888优播| 老女人水多毛片| 2021天堂中文幕一二区在线观| av国产久精品久网站免费入址| a级一级毛片免费在线观看| 在线观看一区二区三区| 免费黄色在线免费观看| 你懂的网址亚洲精品在线观看| 舔av片在线| 国产乱人偷精品视频| 亚洲一区二区三区欧美精品 | 日韩av不卡免费在线播放| 一二三四中文在线观看免费高清| 最近最新中文字幕大全电影3| 欧美区成人在线视频| 欧美+日韩+精品| 99热全是精品| 干丝袜人妻中文字幕| 欧美97在线视频| 制服丝袜香蕉在线| 亚洲自偷自拍三级| 麻豆成人午夜福利视频| 国产淫语在线视频| 国产黄色免费在线视频| 久久久精品欧美日韩精品| 天堂俺去俺来也www色官网| 国产精品不卡视频一区二区| 一区二区三区精品91| 青青草视频在线视频观看| 又爽又黄无遮挡网站| 天天一区二区日本电影三级| 免费观看a级毛片全部| 午夜福利视频精品| 久久精品国产亚洲av天美| 久久久亚洲精品成人影院| 国模一区二区三区四区视频| 黄色一级大片看看| 欧美国产精品一级二级三级 | 大香蕉97超碰在线| av天堂中文字幕网| 色播亚洲综合网| 国产午夜精品一二区理论片| 一级片'在线观看视频| 亚洲av不卡在线观看| 在线观看av片永久免费下载| 老司机影院成人| 黄色日韩在线| 91在线精品国自产拍蜜月| 成人午夜精彩视频在线观看| 99久久中文字幕三级久久日本| 国产一级毛片在线| 免费看a级黄色片| 男人爽女人下面视频在线观看| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| 国产精品麻豆人妻色哟哟久久| 看非洲黑人一级黄片| 免费黄色在线免费观看| 国产精品一区二区性色av| 欧美zozozo另类| 国产精品不卡视频一区二区| 亚洲内射少妇av| 涩涩av久久男人的天堂| 日韩制服骚丝袜av| 久久久久久久大尺度免费视频| 亚洲精品日韩在线中文字幕| 99久国产av精品国产电影| 久久久久久久久久久免费av| 欧美激情在线99| 熟女av电影| 国产亚洲91精品色在线| 一级爰片在线观看| 精品久久久久久电影网| 日本欧美国产在线视频| 日日啪夜夜爽| 黄色怎么调成土黄色| 男人舔奶头视频| 搡女人真爽免费视频火全软件| 一本久久精品| 男人添女人高潮全过程视频| 久久热精品热| 国产综合精华液| 久久久久久久国产电影| 人体艺术视频欧美日本| 国产亚洲午夜精品一区二区久久 | 精品国产一区二区三区久久久樱花 | kizo精华| 精品国产乱码久久久久久小说| 最后的刺客免费高清国语| 黄色配什么色好看| av国产精品久久久久影院| 国产又色又爽无遮挡免| 免费不卡的大黄色大毛片视频在线观看| 中文字幕制服av| 街头女战士在线观看网站| 精品人妻熟女av久视频| 午夜福利高清视频| 下体分泌物呈黄色| 久久久久久九九精品二区国产| 久久国内精品自在自线图片| 一级毛片黄色毛片免费观看视频| 一本久久精品| 肉色欧美久久久久久久蜜桃 | 亚洲精品aⅴ在线观看| 日韩免费高清中文字幕av| 黄色一级大片看看| 国产av国产精品国产| 亚洲无线观看免费| 久久人人爽av亚洲精品天堂 | 神马国产精品三级电影在线观看| 成年人午夜在线观看视频| 在线天堂最新版资源| 亚洲高清免费不卡视频| 国产男女超爽视频在线观看| 久久午夜福利片| 少妇高潮的动态图| 女人十人毛片免费观看3o分钟| 人妻一区二区av| 97在线视频观看| 精品一区二区三卡| 精品久久久久久久久av| 久热这里只有精品99| 自拍欧美九色日韩亚洲蝌蚪91 | 18禁动态无遮挡网站| 中文欧美无线码| 亚洲欧美日韩另类电影网站 | 日韩伦理黄色片| 婷婷色av中文字幕| 成人黄色视频免费在线看| 免费在线观看成人毛片| 一区二区三区精品91| 亚洲精品第二区| 男的添女的下面高潮视频| 男人爽女人下面视频在线观看| 永久网站在线| 中文乱码字字幕精品一区二区三区| 九九在线视频观看精品| 婷婷色综合www| 99视频精品全部免费 在线| 好男人视频免费观看在线| 国产成人精品福利久久| 性插视频无遮挡在线免费观看| 大香蕉久久网| 大片免费播放器 马上看| 网址你懂的国产日韩在线| 91久久精品国产一区二区成人| 久久鲁丝午夜福利片| 天美传媒精品一区二区| 岛国毛片在线播放| 18禁在线无遮挡免费观看视频| 欧美激情久久久久久爽电影| 日日摸夜夜添夜夜爱| 黄色一级大片看看| 好男人视频免费观看在线| 免费av不卡在线播放| 好男人视频免费观看在线| 波野结衣二区三区在线| 插阴视频在线观看视频| 在线天堂最新版资源| 美女cb高潮喷水在线观看| 肉色欧美久久久久久久蜜桃 | 免费看a级黄色片| 青春草国产在线视频| 免费看a级黄色片| 日本黄色片子视频| 欧美日韩视频精品一区| 熟妇人妻不卡中文字幕| 国产成人免费观看mmmm| 夜夜爽夜夜爽视频| 精品久久国产蜜桃| 亚洲精品一二三| 亚洲最大成人av| 日韩在线高清观看一区二区三区| 观看美女的网站| 亚洲,欧美,日韩| 成人特级av手机在线观看| 欧美高清性xxxxhd video| 国产精品国产三级国产专区5o| 国产精品久久久久久精品电影小说 | 一级毛片久久久久久久久女| 乱码一卡2卡4卡精品| 国产午夜精品久久久久久一区二区三区| 美女cb高潮喷水在线观看| 内射极品少妇av片p| 九色成人免费人妻av| 别揉我奶头 嗯啊视频| av在线观看视频网站免费| 美女cb高潮喷水在线观看| 久久鲁丝午夜福利片| a级毛片免费高清观看在线播放| 精品人妻一区二区三区麻豆| av播播在线观看一区| 国产老妇伦熟女老妇高清| 成人亚洲欧美一区二区av| 成人无遮挡网站| 肉色欧美久久久久久久蜜桃 | 精品人妻视频免费看| 国产精品蜜桃在线观看| 亚洲国产精品专区欧美| 黄色欧美视频在线观看| 综合色av麻豆| 欧美激情国产日韩精品一区| 天天躁日日操中文字幕| 日本-黄色视频高清免费观看| 毛片女人毛片| 日韩一本色道免费dvd| 建设人人有责人人尽责人人享有的 | 国产免费一级a男人的天堂| 人妻少妇偷人精品九色| 一区二区三区精品91| 国产v大片淫在线免费观看| 又大又黄又爽视频免费| videos熟女内射| 亚洲自偷自拍三级| 国产又色又爽无遮挡免| av在线播放精品| 成人免费观看视频高清| 亚洲综合精品二区| 国产高潮美女av| 久久精品国产亚洲网站| 久久精品国产自在天天线| 日本三级黄在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲三级黄色毛片| 青春草亚洲视频在线观看| 伊人久久精品亚洲午夜| 性色avwww在线观看| 人妻夜夜爽99麻豆av| 亚洲精品成人av观看孕妇| 日本与韩国留学比较| 国产欧美日韩一区二区三区在线 | 国产伦精品一区二区三区四那| 舔av片在线| 成人高潮视频无遮挡免费网站| 久久午夜福利片| 亚洲第一区二区三区不卡| 熟女人妻精品中文字幕| 亚洲欧美日韩无卡精品| 97在线人人人人妻| 一级a做视频免费观看| 亚洲精华国产精华液的使用体验| 免费观看无遮挡的男女| 欧美xxⅹ黑人| 久久99热6这里只有精品| 国产免费视频播放在线视频| 美女国产视频在线观看| 日韩免费高清中文字幕av| 99久久九九国产精品国产免费| 舔av片在线| 联通29元200g的流量卡| 日韩制服骚丝袜av| 国产欧美日韩一区二区三区在线 | 中文字幕亚洲精品专区| 亚洲精品自拍成人| 国产片特级美女逼逼视频| 亚洲国产高清在线一区二区三| 亚洲国产精品成人综合色| 精品少妇黑人巨大在线播放| 舔av片在线| 免费av观看视频| 97在线视频观看| 欧美潮喷喷水| 中文精品一卡2卡3卡4更新| 在线观看人妻少妇| 国产成人91sexporn| 蜜桃久久精品国产亚洲av| 97人妻精品一区二区三区麻豆| 别揉我奶头 嗯啊视频| 午夜激情久久久久久久| 特级一级黄色大片| 亚洲国产最新在线播放| 亚洲精品视频女| 嘟嘟电影网在线观看| 久久久久久久久大av| 国产淫语在线视频| 91在线精品国自产拍蜜月| 高清午夜精品一区二区三区| 亚洲人与动物交配视频| 免费播放大片免费观看视频在线观看| 亚洲三级黄色毛片| 校园人妻丝袜中文字幕| 成年女人看的毛片在线观看| 色视频www国产| 亚洲av成人精品一二三区| 男人添女人高潮全过程视频| 黄色怎么调成土黄色| 国产精品久久久久久久电影| 免费在线观看成人毛片| 日韩制服骚丝袜av| 另类亚洲欧美激情| 毛片一级片免费看久久久久| 交换朋友夫妻互换小说| 亚洲av不卡在线观看| 免费不卡的大黄色大毛片视频在线观看| 18+在线观看网站| 天天躁日日操中文字幕| 日韩一本色道免费dvd| 欧美极品一区二区三区四区| 少妇的逼好多水| 天天一区二区日本电影三级| 99久国产av精品国产电影| av播播在线观看一区| 老师上课跳d突然被开到最大视频| 亚洲精品国产成人久久av| 成年版毛片免费区| 熟女人妻精品中文字幕| 日韩不卡一区二区三区视频在线| av在线app专区| 我的女老师完整版在线观看| 热re99久久精品国产66热6| 亚洲在线观看片| 久久久久久久久大av| 极品少妇高潮喷水抽搐| 国产精品国产三级国产av玫瑰| 可以在线观看毛片的网站| 我要看日韩黄色一级片| 国产欧美亚洲国产| 天天躁日日操中文字幕| 最近中文字幕2019免费版| 精品一区二区三区视频在线| 汤姆久久久久久久影院中文字幕| 国产一级毛片在线| 亚洲国产av新网站| 国产 精品1| 自拍偷自拍亚洲精品老妇| 亚洲精品456在线播放app| 亚洲av成人精品一二三区| 涩涩av久久男人的天堂| 国产一区亚洲一区在线观看| 一级毛片aaaaaa免费看小| 搡老乐熟女国产| 大香蕉久久网| 亚洲欧美清纯卡通| 在线观看av片永久免费下载| 99热全是精品| 久久精品国产自在天天线| 九九爱精品视频在线观看| 免费av毛片视频| 国产免费又黄又爽又色| 一级a做视频免费观看| 久久亚洲国产成人精品v| 校园人妻丝袜中文字幕| 一个人看的www免费观看视频| 男男h啪啪无遮挡| 久久久久久久久久久免费av| 91久久精品国产一区二区成人| 一个人看的www免费观看视频| videossex国产| 1000部很黄的大片| a级毛色黄片| 美女主播在线视频| 尤物成人国产欧美一区二区三区| 一级毛片 在线播放| av专区在线播放| 午夜亚洲福利在线播放| 国产黄片美女视频| 国产亚洲av嫩草精品影院| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 久久99蜜桃精品久久| 国产v大片淫在线免费观看| 男女下面进入的视频免费午夜| 婷婷色av中文字幕| 久久精品人妻少妇| 嘟嘟电影网在线观看| 国产精品秋霞免费鲁丝片| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 秋霞伦理黄片| 在线观看美女被高潮喷水网站| 国产探花在线观看一区二区| 国产亚洲午夜精品一区二区久久 | 免费看a级黄色片| 蜜桃亚洲精品一区二区三区| 黄片wwwwww| 涩涩av久久男人的天堂| 一级片'在线观看视频| 免费av观看视频| 女人十人毛片免费观看3o分钟| 国产高清不卡午夜福利| 久久久久网色| 只有这里有精品99| 晚上一个人看的免费电影| 少妇人妻久久综合中文| 国产 一区 欧美 日韩| 又粗又硬又长又爽又黄的视频| 秋霞伦理黄片| 欧美日韩视频高清一区二区三区二| 看十八女毛片水多多多| 禁无遮挡网站| 男女边摸边吃奶| 香蕉精品网在线| 国产精品国产三级国产av玫瑰| 波多野结衣巨乳人妻| 综合色丁香网| 最近中文字幕2019免费版| 在线 av 中文字幕| 中文欧美无线码| 直男gayav资源| 啦啦啦在线观看免费高清www| 欧美老熟妇乱子伦牲交| 91精品伊人久久大香线蕉| 身体一侧抽搐| 啦啦啦啦在线视频资源| 亚洲国产av新网站| 99久久精品一区二区三区|