• <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
    黄片播放在线免费| 国产欧美日韩综合在线一区二区| 欧美日韩视频精品一区| 欧美日韩一区二区视频在线观看视频在线| 黄色配什么色好看| 99热这里只有是精品在线观看| 久久影院123| 国产福利在线免费观看视频| 欧美日本中文国产一区发布| 午夜福利,免费看| 深夜精品福利| 国产精品.久久久| 十八禁网站网址无遮挡| 一级爰片在线观看| 日韩制服丝袜自拍偷拍| 国产免费视频播放在线视频| 精品午夜福利在线看| 欧美精品国产亚洲| 亚洲成国产人片在线观看| 最新中文字幕久久久久| 欧美激情 高清一区二区三区| 久久这里有精品视频免费| 国产一区二区三区综合在线观看 | 国产欧美另类精品又又久久亚洲欧美| 中文字幕av电影在线播放| 51国产日韩欧美| 国产精品久久久av美女十八| 精品人妻熟女毛片av久久网站| 久久久欧美国产精品| 亚洲精品美女久久久久99蜜臀 | av.在线天堂| 亚洲中文av在线| 欧美激情国产日韩精品一区| 国产熟女欧美一区二区| 国产有黄有色有爽视频| 精品国产露脸久久av麻豆| 一级毛片黄色毛片免费观看视频| 精品熟女少妇av免费看| 国产有黄有色有爽视频| 精品国产一区二区三区久久久樱花| 亚洲第一区二区三区不卡| 欧美精品高潮呻吟av久久| 国产成人午夜福利电影在线观看| 欧美日韩国产mv在线观看视频| 久久人人爽av亚洲精品天堂| 卡戴珊不雅视频在线播放| 国产亚洲最大av| 91在线精品国自产拍蜜月| 国产精品久久久久成人av| 国产精品久久久久久精品电影小说| 午夜精品国产一区二区电影| 成年动漫av网址| 精品福利永久在线观看| 男人舔女人的私密视频| 欧美人与善性xxx| 成人手机av| 亚洲国产精品一区三区| 久久久久久久大尺度免费视频| 亚洲国产日韩一区二区| 日韩制服骚丝袜av| 国产精品嫩草影院av在线观看| 免费黄频网站在线观看国产| 99视频精品全部免费 在线| 亚洲伊人久久精品综合| 丰满少妇做爰视频| 91成人精品电影| 美女脱内裤让男人舔精品视频| 国产一区二区在线观看av| 久久韩国三级中文字幕| 免费黄频网站在线观看国产| 22中文网久久字幕| 777米奇影视久久| 久久国内精品自在自线图片| 80岁老熟妇乱子伦牲交| 国产精品无大码| 天堂8中文在线网| 超碰97精品在线观看| 精品人妻熟女毛片av久久网站| 五月伊人婷婷丁香| 亚洲综合精品二区| 丝袜喷水一区| 男女边摸边吃奶| 国产片内射在线| 三级国产精品片| 伦精品一区二区三区| 国产 一区精品| 九色成人免费人妻av| 亚洲欧美日韩卡通动漫| 成人无遮挡网站| 精品国产国语对白av| 欧美日韩国产mv在线观看视频| 午夜福利影视在线免费观看| 亚洲精品国产av成人精品| 国内精品宾馆在线| 久久国产亚洲av麻豆专区| 女性被躁到高潮视频| 亚洲国产欧美日韩在线播放| 男人添女人高潮全过程视频| 在线天堂中文资源库| 亚洲中文av在线| 九九在线视频观看精品| 精品福利永久在线观看| 婷婷色麻豆天堂久久| 69精品国产乱码久久久| xxx大片免费视频| 成人18禁高潮啪啪吃奶动态图| 国产综合精华液| 91精品国产国语对白视频| av不卡在线播放| 久久久久精品性色| 人妻少妇偷人精品九色| 亚洲欧美一区二区三区国产| 亚洲精品456在线播放app| 在线 av 中文字幕| a级毛色黄片| 黑人巨大精品欧美一区二区蜜桃 | 69精品国产乱码久久久| 亚洲国产欧美在线一区| av天堂久久9| 在线观看国产h片| 美女内射精品一级片tv| 欧美精品亚洲一区二区| 性色av一级| kizo精华| 欧美xxⅹ黑人| 22中文网久久字幕| 亚洲精品乱码久久久久久按摩| 99热国产这里只有精品6| 精品一区在线观看国产| 国产精品熟女久久久久浪| 在线观看一区二区三区激情| 夫妻性生交免费视频一级片| 90打野战视频偷拍视频| 欧美xxⅹ黑人| 日本午夜av视频| xxx大片免费视频| 亚洲精品久久午夜乱码| 欧美精品国产亚洲| 尾随美女入室| a级片在线免费高清观看视频| 黄色一级大片看看| 亚洲第一区二区三区不卡| 欧美精品人与动牲交sv欧美| 精品一区二区免费观看| 多毛熟女@视频| 久久久久久久精品精品| 99久国产av精品国产电影| 天天躁夜夜躁狠狠躁躁| 看免费成人av毛片| 22中文网久久字幕| 久久精品国产自在天天线| 中文精品一卡2卡3卡4更新| 亚洲欧美日韩卡通动漫| 亚洲国产最新在线播放| 99热全是精品| av黄色大香蕉| 两个人免费观看高清视频| 女性生殖器流出的白浆| 中文天堂在线官网| 亚洲国产精品专区欧美| 亚洲婷婷狠狠爱综合网| 国产高清三级在线| 男女下面插进去视频免费观看 | 国产免费福利视频在线观看| 人成视频在线观看免费观看| 国产一区二区在线观看av| 一级a做视频免费观看| 久久 成人 亚洲| 赤兔流量卡办理| 久久97久久精品| 国产成人精品在线电影| 国产黄色视频一区二区在线观看| 日日啪夜夜爽| 女性生殖器流出的白浆| 欧美亚洲 丝袜 人妻 在线| 免费久久久久久久精品成人欧美视频 | 亚洲欧美中文字幕日韩二区| 女性被躁到高潮视频| 日日啪夜夜爽| 大陆偷拍与自拍| 日韩精品有码人妻一区| 99九九在线精品视频| 香蕉精品网在线| 97超碰精品成人国产| 最近最新中文字幕大全免费视频 | 亚洲国产av影院在线观看| 亚洲成国产人片在线观看| 久久国产亚洲av麻豆专区| 最后的刺客免费高清国语| 欧美成人午夜免费资源| 免费观看在线日韩| 黑人猛操日本美女一级片| 蜜桃国产av成人99| 久久女婷五月综合色啪小说| 日日撸夜夜添| 一级毛片电影观看| 一区二区三区乱码不卡18| 国产在线一区二区三区精| 精品人妻熟女毛片av久久网站| 亚洲国产精品一区二区三区在线| 亚洲欧美成人综合另类久久久| 午夜福利视频精品| 国产成人欧美| 免费少妇av软件| 考比视频在线观看| 欧美最新免费一区二区三区| 满18在线观看网站| 一本色道久久久久久精品综合| 午夜91福利影院| 老女人水多毛片| 性高湖久久久久久久久免费观看| 伊人亚洲综合成人网| 草草在线视频免费看| 精品久久国产蜜桃| 一级毛片我不卡| 欧美亚洲日本最大视频资源| 一区二区日韩欧美中文字幕 | 亚洲国产av新网站| 中文精品一卡2卡3卡4更新| 久久国产亚洲av麻豆专区| 久久久久国产精品人妻一区二区| 精品国产一区二区三区四区第35| 精品熟女少妇av免费看| 久久久国产一区二区| 一本久久精品| 乱码一卡2卡4卡精品| 在线天堂最新版资源| 一区二区av电影网| 国产亚洲最大av| 侵犯人妻中文字幕一二三四区| 亚洲性久久影院| 观看av在线不卡| 中文字幕最新亚洲高清| 美国免费a级毛片| 在线观看人妻少妇| 欧美最新免费一区二区三区| 日韩中字成人| 少妇的丰满在线观看| 26uuu在线亚洲综合色| 欧美国产精品一级二级三级| 婷婷色综合www| 老熟女久久久| www.av在线官网国产| 少妇被粗大猛烈的视频| 亚洲精品久久午夜乱码| 丝袜脚勾引网站| 精品亚洲成国产av| 国产精品久久久久久久电影| 亚洲欧洲国产日韩| 卡戴珊不雅视频在线播放| 另类精品久久| 欧美日本中文国产一区发布| 亚洲天堂av无毛| av天堂久久9| 亚洲欧美日韩另类电影网站| 老女人水多毛片| 国产女主播在线喷水免费视频网站| 日韩大片免费观看网站| av卡一久久| 欧美成人午夜精品| 久久毛片免费看一区二区三区| 日本与韩国留学比较| 色哟哟·www| 欧美国产精品va在线观看不卡| 精品视频人人做人人爽| 视频在线观看一区二区三区| 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说| 人妻一区二区av| 26uuu在线亚洲综合色| 精品第一国产精品| 丝袜在线中文字幕| 久热这里只有精品99| 国产福利在线免费观看视频| √禁漫天堂资源中文www| 国产在视频线精品| 乱人伦中国视频| 18+在线观看网站| 狠狠婷婷综合久久久久久88av| 看非洲黑人一级黄片| 亚洲婷婷狠狠爱综合网| 新久久久久国产一级毛片| 永久网站在线| 男人添女人高潮全过程视频| 欧美精品一区二区大全| 久久热在线av| 飞空精品影院首页| 九草在线视频观看| 成人黄色视频免费在线看| 亚洲高清免费不卡视频| av在线app专区| 免费观看av网站的网址| 国产精品三级大全| 久久人人爽人人片av| 十八禁网站网址无遮挡| 又粗又硬又长又爽又黄的视频| 日韩 亚洲 欧美在线| 赤兔流量卡办理| 免费人妻精品一区二区三区视频| 国产男女内射视频| 欧美日韩视频高清一区二区三区二| 下体分泌物呈黄色| 国产在线一区二区三区精| 狠狠婷婷综合久久久久久88av| 五月玫瑰六月丁香| 建设人人有责人人尽责人人享有的| 色5月婷婷丁香| 国产国语露脸激情在线看| 18禁在线无遮挡免费观看视频| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 我的女老师完整版在线观看| 中文乱码字字幕精品一区二区三区| 亚洲三级黄色毛片| 亚洲国产精品一区二区三区在线| 国产激情久久老熟女| 日韩一区二区三区影片| 成年人免费黄色播放视频| 亚洲av.av天堂| 丝袜脚勾引网站| 亚洲第一av免费看| 视频在线观看一区二区三区| 黄片播放在线免费| 韩国精品一区二区三区 | 最近中文字幕高清免费大全6| 少妇的丰满在线观看| 日本av手机在线免费观看| 久久久久久人妻| 中文字幕人妻丝袜制服| 18禁裸乳无遮挡动漫免费视频| 国产 精品1| 精品少妇黑人巨大在线播放| 亚洲国产精品一区三区| 午夜福利,免费看| 国产亚洲最大av| 人妻 亚洲 视频| 大片电影免费在线观看免费| 青春草视频在线免费观看| 少妇的丰满在线观看| 亚洲国产精品国产精品| 一级毛片 在线播放| 高清毛片免费看| 精品国产乱码久久久久久小说| 亚洲成人一二三区av| 亚洲精品,欧美精品| 久久久久国产精品人妻一区二区| 久久久久久久久久久久大奶| 久久久久网色| 日韩人妻精品一区2区三区| videosex国产| 色哟哟·www| 亚洲三级黄色毛片| 国产精品三级大全| 国产又色又爽无遮挡免| 久久热在线av| 十八禁网站网址无遮挡| 中文字幕av电影在线播放| 成人午夜精彩视频在线观看| 精品国产一区二区久久| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 高清在线视频一区二区三区| 最近最新中文字幕免费大全7| 韩国精品一区二区三区 | 色视频在线一区二区三区| 色5月婷婷丁香| 国产亚洲欧美精品永久| 亚洲图色成人| 五月伊人婷婷丁香| 如何舔出高潮| 搡女人真爽免费视频火全软件| 亚洲国产成人一精品久久久| 日本wwww免费看| 老熟女久久久| 各种免费的搞黄视频| 日本91视频免费播放| 99久久中文字幕三级久久日本| av在线app专区| 七月丁香在线播放| 国产乱人偷精品视频| 精品一区二区免费观看| av在线老鸭窝| 又粗又硬又长又爽又黄的视频| 亚洲三级黄色毛片| www.熟女人妻精品国产 | 91精品国产国语对白视频| 国产精品.久久久| 免费观看无遮挡的男女| 成人亚洲欧美一区二区av| 亚洲国产最新在线播放| 中文字幕人妻丝袜制服| 亚洲国产精品成人久久小说| 亚洲精品乱码久久久久久按摩| 久久久久国产精品人妻一区二区| 日本与韩国留学比较| 亚洲av综合色区一区| 一区在线观看完整版| 女的被弄到高潮叫床怎么办| 极品人妻少妇av视频| 高清在线视频一区二区三区| 亚洲国产欧美在线一区| 少妇被粗大的猛进出69影院 | 久久久a久久爽久久v久久| 桃花免费在线播放| 免费少妇av软件| 午夜激情久久久久久久| 久久久精品免费免费高清| 国产亚洲最大av| 免费在线观看完整版高清| 天美传媒精品一区二区| 成年人免费黄色播放视频| 久久韩国三级中文字幕| 九九在线视频观看精品| 美女国产高潮福利片在线看| 日韩av不卡免费在线播放| a级毛色黄片| 一区在线观看完整版| 男人舔女人的私密视频| 日韩大片免费观看网站| 久久精品aⅴ一区二区三区四区 | 精品少妇内射三级| 久久人妻熟女aⅴ| 精品一区二区三区视频在线| 亚洲一区二区三区欧美精品| 亚洲精品456在线播放app| 又大又黄又爽视频免费| 国产精品麻豆人妻色哟哟久久| xxx大片免费视频| 人妻 亚洲 视频| 亚洲国产精品一区三区| 国产一区二区激情短视频 | 国产亚洲精品第一综合不卡 | 夜夜爽夜夜爽视频| 免费在线观看黄色视频的| 侵犯人妻中文字幕一二三四区| 中文字幕另类日韩欧美亚洲嫩草| 成人影院久久| 成年女人在线观看亚洲视频| 国产熟女午夜一区二区三区| 老司机亚洲免费影院| 黄色配什么色好看| 在线观看国产h片| 国产又爽黄色视频| 极品人妻少妇av视频| 久久久久久久大尺度免费视频| 亚洲av国产av综合av卡| 老司机影院成人| 久久久久久久久久人人人人人人| 亚洲情色 制服丝袜| 精品视频人人做人人爽| 中国国产av一级| 看免费成人av毛片| 婷婷色综合www| 精品人妻偷拍中文字幕| 国产成人欧美| 99久久中文字幕三级久久日本| 青春草视频在线免费观看| 国产高清国产精品国产三级| www.av在线官网国产| 国语对白做爰xxxⅹ性视频网站| 最近最新中文字幕大全免费视频 | 午夜av观看不卡| 国产深夜福利视频在线观看| 亚洲精品中文字幕在线视频| 成人国产麻豆网| 日本91视频免费播放| √禁漫天堂资源中文www| 久久久久久久久久成人| 亚洲国产最新在线播放| 天天躁夜夜躁狠狠久久av| 激情视频va一区二区三区| 99热国产这里只有精品6| tube8黄色片| √禁漫天堂资源中文www| 亚洲成av片中文字幕在线观看 | 亚洲精品日本国产第一区| 丝袜人妻中文字幕| 精品99又大又爽又粗少妇毛片| 爱豆传媒免费全集在线观看| 宅男免费午夜| 好男人视频免费观看在线| 亚洲欧美精品自产自拍| 激情视频va一区二区三区| 国产精品秋霞免费鲁丝片| 精品一区二区三卡| 午夜免费鲁丝| 免费观看av网站的网址| 久久99热6这里只有精品| 一区二区av电影网| 久久精品熟女亚洲av麻豆精品| 成人漫画全彩无遮挡| av有码第一页| 亚洲国产色片| 免费观看a级毛片全部| 国产成人午夜福利电影在线观看| 有码 亚洲区| 亚洲精品美女久久久久99蜜臀 | 国产女主播在线喷水免费视频网站| 欧美人与性动交α欧美软件 | 亚洲国产最新在线播放| 国产av码专区亚洲av| 性色avwww在线观看| 国产av精品麻豆| 咕卡用的链子| 精品国产露脸久久av麻豆| 人妻人人澡人人爽人人| 久久人妻熟女aⅴ| av不卡在线播放| 最近最新中文字幕大全免费视频 | 日韩精品免费视频一区二区三区 | 久久精品久久久久久久性| 国产黄频视频在线观看| 国产成人欧美| 日韩精品有码人妻一区| 十八禁高潮呻吟视频| 日日撸夜夜添| 日本欧美视频一区| 大片电影免费在线观看免费| 两个人免费观看高清视频| 国产视频首页在线观看| 高清视频免费观看一区二区| 亚洲色图综合在线观看| 国产免费一级a男人的天堂| 亚洲av综合色区一区| 爱豆传媒免费全集在线观看| 亚洲久久久国产精品| 春色校园在线视频观看| 中文欧美无线码| 涩涩av久久男人的天堂| 丰满迷人的少妇在线观看| 老司机影院毛片| 一级黄片播放器| 色婷婷av一区二区三区视频| 久久亚洲国产成人精品v| 精品久久久精品久久久| 最新的欧美精品一区二区| 亚洲av在线观看美女高潮| 91久久精品国产一区二区三区| 国产精品一区二区在线观看99| 欧美xxxx性猛交bbbb| 精品国产一区二区久久| 国产不卡av网站在线观看| 欧美日韩国产mv在线观看视频| 最新的欧美精品一区二区| videos熟女内射| 9色porny在线观看| 一级毛片黄色毛片免费观看视频| 另类精品久久| 国产成人精品婷婷| 五月开心婷婷网| 最近的中文字幕免费完整| 少妇人妻 视频| 久久久欧美国产精品| 丰满少妇做爰视频| 欧美精品一区二区免费开放| 国产免费又黄又爽又色| 亚洲,一卡二卡三卡| 高清黄色对白视频在线免费看| 久久久久久人人人人人| 国产激情久久老熟女| 人人妻人人澡人人爽人人夜夜| 在线免费观看不下载黄p国产| 亚洲综合色惰| 欧美激情国产日韩精品一区| 男人爽女人下面视频在线观看| 免费播放大片免费观看视频在线观看| 久久精品熟女亚洲av麻豆精品| 国产成人a∨麻豆精品| 这个男人来自地球电影免费观看 | 久久久久久人人人人人| 97在线视频观看| 精品国产一区二区三区久久久樱花| 欧美丝袜亚洲另类| 久久人人爽人人爽人人片va| www.av在线官网国产| 最近的中文字幕免费完整| av片东京热男人的天堂| 久久久久久久久久成人| 曰老女人黄片| 精品一区二区免费观看| 欧美+日韩+精品| 男人爽女人下面视频在线观看| 一区二区av电影网| 色5月婷婷丁香| 中文乱码字字幕精品一区二区三区| 人人妻人人澡人人看| kizo精华| 国产精品嫩草影院av在线观看| 国产国语露脸激情在线看| 国产精品不卡视频一区二区| 2021少妇久久久久久久久久久| 精品国产乱码久久久久久小说| 久久毛片免费看一区二区三区| 99久久中文字幕三级久久日本| 人人妻人人爽人人添夜夜欢视频| 男人爽女人下面视频在线观看| 久久这里只有精品19| 人人妻人人爽人人添夜夜欢视频| 亚洲国产色片| 99热网站在线观看| 免费看光身美女| 午夜激情av网站| 成人国产av品久久久| 晚上一个人看的免费电影| 国国产精品蜜臀av免费| 人体艺术视频欧美日本| 丝袜在线中文字幕| www.熟女人妻精品国产 | 久久精品久久久久久久性| 日韩av免费高清视频| 国产成人精品在线电影| 另类精品久久| 黑人猛操日本美女一级片| 狠狠婷婷综合久久久久久88av| 精品99又大又爽又粗少妇毛片|