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

    滲透指數(shù)在黃土丘陵溝壑區(qū)徑流輸沙與覆被格局關(guān)系研究中的適用性

    2017-03-23 03:48:50
    中國(guó)農(nóng)村水利水電 2017年11期
    關(guān)鍵詞:輸沙覆蓋度柵格

    劉 宇

    (中國(guó)科學(xué)院地理科學(xué)與資源研究所 中國(guó)科學(xué)院生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測(cè)與模擬重點(diǎn)實(shí)驗(yàn)室,北京 100101)

    植被覆蓋格局是流域水土流失的重要控制因素。從徑流小區(qū)、坡面到流域的實(shí)驗(yàn)研究表明,植被覆蓋的空間分布格局是侵蝕產(chǎn)沙、泥沙輸移、地表產(chǎn)流、徑流量的調(diào)控因素[1]。因此,定量刻畫植被格局的水土保持效應(yīng)是水土流失研究的熱點(diǎn)議題之一[2]。在斑塊-基底-廊道景觀范式下發(fā)展起來(lái)的斑塊密度、斑塊性狀、景觀破碎度、景觀連通度等景觀指數(shù)被大量應(yīng)用到反映給定土地覆被/利用格局下的水土流失[2]。然而,這些景觀指數(shù)中僅有部分指數(shù)能反映水土流失的某些方面,大多數(shù)指數(shù)與水土流失過程表征變量之間沒有明確的統(tǒng)計(jì)關(guān)系[3],無(wú)法對(duì)水土流失做出準(zhǔn)確的預(yù)測(cè),因而難以用它們的變化預(yù)測(cè)給定覆被格局下的水土流失發(fā)展趨勢(shì)。水土流失過程的方向性與景觀指數(shù)各向同性的差別等原因是景觀指數(shù)不能準(zhǔn)確反映水土流失狀況的重要原因。針對(duì)景觀指數(shù)在水土流失研究中的局限性,一些研究開始從水土流失過程機(jī)理的理解入手,構(gòu)建針對(duì)水土流失過程的格局指數(shù),將覆被格局及其變化與水土流失聯(lián)系起來(lái)[4,5]。在水力侵蝕區(qū),水沙輸移沿地形梯度方向進(jìn)行,即從上坡向下坡,從上游向下游。植被斑塊和地形洼地等是泥沙輸移的主要阻礙。在半干旱區(qū)草地生態(tài)系統(tǒng)實(shí)驗(yàn)觀測(cè)的基礎(chǔ)上,Ludwig等[6]結(jié)合植被覆蓋度的空間分布和地形格局,提出了滲透指數(shù)(LI)來(lái)定量不同植被覆蓋格局的泥沙和徑流截留功能。從構(gòu)建原理上看LI指數(shù)考慮了水土流失的方向性及其與植被覆蓋度之間的非線性關(guān)系,是針對(duì)水土流失過程而發(fā)展起來(lái)的指數(shù),在澳大利亞半干旱區(qū)草原生態(tài)系統(tǒng)中得到驗(yàn)證[6]。由于其是針對(duì)澳大利亞半干旱區(qū)草地生態(tài)系統(tǒng)而建立,對(duì)林、灌木、草等多植被類型鑲嵌的覆被格局和在其他區(qū)域的適用性還需進(jìn)一步驗(yàn)證。黃土高原是我國(guó)水土流失最為嚴(yán)重的地區(qū),黃土丘陵溝壑區(qū)尤為突出[7]。過去數(shù)十年的植被恢復(fù)形成林、灌、草、農(nóng)作物鑲嵌的植被覆蓋格局??紤]到LI指數(shù)耦合覆被格局與水沙輸移的特點(diǎn)和其最初發(fā)展的環(huán)境背景,本研究旨在流域尺度探討LI指數(shù)在黃土高原丘陵溝壑區(qū)的適用性。

    1 研究區(qū)概況

    本研究選擇了黃河右岸黃土丘陵溝壑區(qū)6個(gè)有長(zhǎng)期水文監(jiān)測(cè)的流域?yàn)檠芯繀^(qū),即皇甫川(皇甫川站)、孤山河(高石崖站)、窟野河(溫家川站)、無(wú)定河(丁家溝站)、北洛河(狀頭站)、禿尾河(高家川站)(圖1)。這些流域地處黃土丘陵溝壑區(qū),地形破碎,溝壑密度大,土壤質(zhì)地疏松,侵蝕強(qiáng)烈,是黃河泥沙的重要來(lái)源區(qū)。本區(qū)屬于暖溫帶半干旱向半濕潤(rùn)氣候過渡區(qū),多年平均降水量372.9~502.0 mm。所選各流域多年平均降水量在360~500 mm之間,蒸發(fā)量在325~452 mm之間(表1)。各流域面積及其控制站點(diǎn)如表1所示。

    表1 研究流域1982-2008平均水文特征Tab.1 Average hydrological characteristic ofthe study catchments during 1982-2008

    圖1 研究流域地理位置示意圖Fig.1 Geographical location of watersheds focused on in this study

    長(zhǎng)期以來(lái)人類活動(dòng)對(duì)植被的改變是黃土丘陵溝壑土壤侵蝕速率變化的主要因素之一。黃土丘陵溝壑區(qū)地形陡峭,坡耕地比重大,是黃土高原退耕環(huán)林還草等植被恢復(fù)的主要區(qū)域。自發(fā)撂荒、退耕還林還草工程和禁牧措施是區(qū)域植被恢復(fù)的主要驅(qū)動(dòng)。在這些行動(dòng)的推動(dòng)下,過去30年來(lái),尤其是在1998年大規(guī)模退耕還林(還草)工程實(shí)施后,流域植被得到較好的恢復(fù),覆蓋度顯著提高。農(nóng)、林、草的空間鑲嵌構(gòu)成流域的基本土地覆蓋格局。

    2 滲透指數(shù)計(jì)算方法

    2.1 植被覆蓋度參數(shù)提取

    2000-2008 年MODIS 數(shù)據(jù)為來(lái)源于美國(guó)國(guó)家宇航局(NASA)EOS 數(shù)據(jù)中心發(fā)布的16 NDVI 合成250 m分辨率數(shù)據(jù)產(chǎn)品。采用像元二分法[8]提取所選流域2000-2008 年植被覆蓋度,計(jì)算公式如式(1)。

    (1)

    式中:FCi為像元i的植被覆蓋度;NDVIi為像元i的NDVI值;NDVIv、NDVIs分別為純植被和純土壤的NDVI值。根據(jù)公式(1)獲得的植被覆蓋度數(shù)據(jù),采用最大值合成法,生成2000-2008年每年的植被覆蓋度。

    2.2 流域滲透指數(shù)計(jì)算

    滲透指數(shù)(LI)是在半干旱草地生態(tài)系統(tǒng)實(shí)驗(yàn)研究的基礎(chǔ)上提出的,基于植被對(duì)侵蝕產(chǎn)沙的抑制和對(duì)泥沙輸送的阻礙作用構(gòu)建。其計(jì)算式如式(2)。

    (2)

    式中:Lcal為流域現(xiàn)實(shí)植被覆蓋(覆蓋度)下的泥沙留滯能力;Lmin和Lmax分別為最小和最大留滯能力,這里分別設(shè)置為植被覆蓋度為100%時(shí)的Lcal,植被覆蓋度為0時(shí)的Lcal;k表征LI與植被覆蓋之間關(guān)系曲線的形狀和陡度。

    由式(2)可知,LI∈[0,1]。LI值越小則土壤流失風(fēng)險(xiǎn)越小,值越大則土壤流失風(fēng)險(xiǎn)越大。Lcal采用八流向累積算法根據(jù)式(3)計(jì)算。

    (3)

    式中:n=1,2,3,4,5,6,7,8,分別代表與柵格(i,j)相鄰的8個(gè)柵格;Sn是柵格n流向柵格(i,j)的比例,由柵格間的坡降確定[9]。由于地形坡度(Si,j)是徑流、泥沙的輸出的主要影響因素,式(3)中將坡度因子Sin(Si,j)作為一個(gè)參數(shù)。當(dāng)把Lij替換為L(zhǎng)ij_min或Lij_max,所得Pi,j即分別為Pij_min、Pij_max。

    (4)

    式中:Pout,j為與目標(biāo)水體鄰接且流向目標(biāo)水體的柵格i逐步累積滲透值,即Lcal為所有鄰接且水流方向指向目標(biāo)水體的柵格累積滲透值之和;m為與目標(biāo)水體鄰接的柵格總數(shù)。

    Lmin、Lmax的計(jì)算與公式(4)相似,將Pout_j用與河道相鄰且水流方向指向河道的Pij_min、Pij_max分別替代即得。

    植被覆蓋度是影響侵蝕產(chǎn)沙、輸沙的重要因子之一。潛在土壤流失風(fēng)險(xiǎn)與植被覆蓋度呈負(fù)指數(shù)關(guān)系[6]。在公式(3)中,Li,j用植被覆蓋度構(gòu)造的公式(5)求算,它表示植被覆蓋度為ci,j(%)的柵格的土壤流失風(fēng)險(xiǎn)。

    Lij=e-b×cij

    (5)

    式中:Li,j為潛在土壤流失風(fēng)險(xiǎn);b為目標(biāo)流域土壤流失隨植被覆蓋度增加而下降的曲線陡度,與流域土壤、地形的特征有關(guān),可利用侵蝕產(chǎn)沙和植被覆蓋度的回歸擬合推算。根據(jù)公式(2)~式(5),按照如圖2所示的流程,利用Arcgis 10.0水文分析功能和Arc/info腳本編程相結(jié)合完成流域LI指數(shù)的計(jì)算。

    圖2 基于ArcGIS的滲透指數(shù)計(jì)算流程Fig.2 Diagram of leakiness index calculation using ArcGIS

    圖2中,F(xiàn)C為植被覆蓋度,%;SP為坡度因子=Sin(slope);FD為水流方向;P_min、P、P_max分別為最小、現(xiàn)實(shí)、最大流失風(fēng)險(xiǎn)權(quán)重下的累積匯流量;Flow(FD,w),w為柵格權(quán)重,這里等于Lij、Lij_min或Lij_max與坡度因子的乘積,累積匯流量計(jì)算用Arcgis desktop 10.0水文分析模塊Flowaccumulation工具完成;河網(wǎng)的提取采用匯流面積閾值法,即根據(jù)實(shí)際河流分布確定產(chǎn)生河道的最小匯流面積閾值,面積大于閾值的即為河網(wǎng);提取與河道相鄰且水流方向指向河道的柵格的P、P_min、P_max時(shí),可在包括多個(gè)流域的整個(gè)研究區(qū)尺度上一起提取,然后用包括多個(gè)流域的流域邊界文件(polygon)分別與提取出的P、P_min、P_max做Zonal Statistics (用Zonal Statistics as table工具),獲得的表格中的sum項(xiàng)即分別為要求的Lcal、L_min、L_max;根據(jù)公式(2)即可得到各個(gè)流域的滲透指數(shù)。

    2.3 數(shù)據(jù)來(lái)源和分析方法

    本研究利用的泥沙數(shù)據(jù)為6個(gè)黃河一級(jí)支流水文站2000-2008年的泥沙監(jiān)測(cè)數(shù)據(jù)。地形數(shù)據(jù)為STRM90m分辨率柵格數(shù)據(jù)。植被指數(shù)(NDVI)、地形數(shù)據(jù)皆來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http:∥www.gscloud.cn)。根據(jù)2.2節(jié)的方法,計(jì)算出每個(gè)流域2000-2008年逐年的滲透指數(shù)。利用流域控制水文站的徑流、輸沙數(shù)據(jù)和滲透指數(shù)之間的非線性回歸關(guān)系,刻畫流域徑流、輸沙與植被格局的統(tǒng)計(jì)關(guān)系。

    3 結(jié)果分析

    3.1 流域植被覆蓋度、滲透指數(shù)、徑流和輸沙變化

    基于遙感的植被覆蓋度監(jiān)測(cè)顯示,2000-2008年,6個(gè)流域平均植被覆蓋度都呈增加趨勢(shì)。最南邊的北洛河降雨較多,平均植被覆蓋度最高,整個(gè)時(shí)期都在70%以上,在時(shí)段末期增加到約85%,遠(yuǎn)高于其他流域。北部的無(wú)定河、皇甫川、窟野河等流域植被覆蓋度較低,此間從30%~40%增加到50%左右(圖3)。

    圖3 2000-2008年流域植被覆蓋度變化Fig.3 Vegetation coverage dynamic during 2000-2008

    植被覆蓋度增加的同時(shí),反映流域植被覆蓋格局水土保持能力的流域滲透指數(shù)也總體上呈逐漸下降的趨勢(shì)(圖4)。除北洛河流域外,其余各流域在2004年出現(xiàn)一個(gè)由降轉(zhuǎn)升的階段,至2006年后轉(zhuǎn)為降低趨勢(shì)。各流域之間的比較表明,北洛河具有最低的滲透指數(shù)值,而無(wú)定河最高。

    圖4 各流域滲透指數(shù)時(shí)間變化Fig.4 Temporal dynamic of leakiness index of selected watersheds

    黃土高原丘陵溝壑區(qū)降雨集中,高強(qiáng)度的暴雨貢獻(xiàn)了徑流、輸沙的絕大部分[10]。降雨年際變異大,進(jìn)而導(dǎo)致流域徑流、輸沙的巨大變化幅度(圖5)。單位降雨輸沙量的變異系數(shù)除窟野河流域外,都在80%以上,最高為北洛河,達(dá)到90.21%;徑流系數(shù)的變異系數(shù)除皇甫川流域外(變異系數(shù)為66.54%),都低于40%,最低為禿尾河,為20.31%(圖6)。

    圖5 流域單位降雨輸沙模數(shù)年際變化Fig.5 Annual variation of sediment yield of unit precipitation in studied watersheds

    圖6 流域徑流系數(shù)年際變化Fig.6 Annual variation of runoff coefficient for each watershed

    3.2 滲透指數(shù)表征流域產(chǎn)流和輸沙能力的有效性

    從結(jié)構(gòu)上看,滲透指數(shù)利用了植被覆蓋度作為參數(shù)構(gòu)建的指標(biāo),與植被覆蓋度之間是負(fù)指數(shù)關(guān)系。利用流域平均植被覆蓋度與滲透指數(shù)的進(jìn)行回歸擬合,得到如圖7所示的結(jié)果,證實(shí)了植被覆蓋度與滲透指數(shù)之間的負(fù)指數(shù)關(guān)系。滲透指數(shù)通過逐柵格的累積算法,不僅僅考慮植被覆蓋度,還將植被空間分布格局的作用隱含進(jìn)來(lái),因而既能反映植被覆蓋度對(duì)流域徑流和泥沙的影響,也能反映植被覆蓋度的空間分布格局變化的水土流失效應(yīng)。

    采用逐流域?qū)?000-2008年滲透指數(shù)與單位降雨輸沙模數(shù)進(jìn)行回歸擬合,結(jié)果顯示二者間存在顯著的指數(shù)關(guān)系(圖8)。滲透指數(shù)對(duì)單位降雨量輸沙模數(shù)的解釋度因流域而異,最低的皇甫川流域,僅為27.36%,其余都在40%以上,最高為孤山河流域,達(dá)到64.58%。這些結(jié)果表明,LI指數(shù)能在一定程度上準(zhǔn)確反映流域植被覆蓋格局變化的輸沙效應(yīng)。

    圖7 滲透指數(shù)與植被覆蓋度之間的關(guān)系Fig.7 Statistical relationships between leakiness index and vegetation coverage

    圖8 流域單位降雨輸沙模數(shù)與滲透指數(shù)的回歸關(guān)系Fig.8 Regression relationship between sediment yield of unit rainfall and leakiness index

    流域徑流反映流域?qū)⒔涤贽D(zhuǎn)化為徑流的能力,受植被覆蓋度及其分布格局的影響。LI指數(shù)與流域徑流系數(shù)間除窟野河(溫家川站)呈顯著的指數(shù)關(guān)系外,其余流域均無(wú)顯著關(guān)系(圖9)。表明LI指數(shù)在表征植被格局變化的徑流效應(yīng)方面存在局限。其原因是滲透指數(shù)基于徑流、泥沙在地表的傳輸過程建立,并未考慮地下徑流的作用,而流域徑流中有相當(dāng)一部分由地下徑流貢獻(xiàn)[11]。加之植被覆蓋對(duì)泥沙傳輸?shù)淖铚饔靡獜?qiáng)于對(duì)徑流的阻滯作用,使得滲透指數(shù)與流域徑流系數(shù)的回歸關(guān)系并不如與單位降雨輸沙量的回歸關(guān)系好。

    在黃土丘陵溝壑區(qū),魚鱗坑、梯田、淤地壩等水土保持措施廣泛分布,而在LI指數(shù)中因數(shù)據(jù)的缺乏而沒有將這些工程性因素的影響考慮在內(nèi),這限制了LI指數(shù)對(duì)給定植被覆蓋格局下流域輸沙能力、特別是流域產(chǎn)流能力的表達(dá)。這也是指數(shù)應(yīng)用效果的重要原因。此外,計(jì)算LI采用的地形數(shù)的分辨率會(huì)影響LI的大小。一般而言,在流域尺度上,低分辨率的地形數(shù)據(jù)導(dǎo)致高估LI。盡管如此,滲透指數(shù)通過逐柵格的累積算法,一定程度上反映了水沙沿地表的傳輸過程,因而是將水沙傳輸與植被覆蓋格局聯(lián)系起來(lái)的有用工具。

    4 結(jié) 語(yǔ)

    本研究結(jié)果顯示LI指數(shù)與流域輸沙模數(shù)、單位降雨輸沙模數(shù)具有很好的回歸關(guān)系,大致呈指數(shù)關(guān)系。這表明LI指數(shù)能準(zhǔn)確表征流域尺度植被覆蓋格局變化的輸沙效應(yīng)。雖然徑流與LI指數(shù)呈正相關(guān),但回歸關(guān)系并不顯著。由此看,LI指數(shù)難以反映該區(qū)域植被格局對(duì)流域徑流影響。LI指數(shù)在表征植被格局對(duì)流域徑流和輸沙的影響方面的差別可歸因于指數(shù)本身的性質(zhì)與流域徑流、輸沙影響因素的多樣性。今后需要將植被結(jié)構(gòu)、微地形改造工程等信息納入進(jìn)來(lái),以增強(qiáng)指數(shù)的有效性。LI指數(shù)是基于水沙在地表的輸移構(gòu)建。而流域徑流量中有相當(dāng)部分是地下徑流。此外,植被對(duì)地表徑流的阻滯作用并不像對(duì)泥沙的阻滯那樣強(qiáng)烈。這些原因使得LI指數(shù)更適合表征流域植被覆蓋格局對(duì)輸沙的影響,而對(duì)植被覆蓋格局對(duì)徑流的影響指示效果較差。由于滲透指數(shù)具有一定的過程意義,可作為評(píng)價(jià)植被配置格局土壤保持、甄選流域植被格局設(shè)計(jì)方案的工具。需要注意的是,該指數(shù)僅利用了植被覆蓋度、地形兩個(gè)方面的參數(shù),是以區(qū)域土壤、降雨等其他因素?zé)o空間異質(zhì)性為前提的。因此,在應(yīng)用時(shí)應(yīng)注意這一點(diǎn)。在可能的條件下,構(gòu)造表征其他因素影響土壤侵蝕的指標(biāo)作為參數(shù)與LI融合時(shí)最佳的選擇。

    圖9 流域徑流系數(shù)與滲透指數(shù)的回歸關(guān)系Fig.9 Regression curves between watershed runoff coefficients and leakiness index

    [1] Puigdefábregas J. The role of vegetation patterns in structuring runoff and sediment fluxes in drylands[J]. Earth Surface Processes and Landforms,2005,30(2):133-147.

    [2] Ouyang Wei, Andrew K S, Hao Fanghua, et al. Soil erosion dynamics response to landscape pattern[J]. Science of the Total Environment, 2010,408:1 358-1 366.

    [3] 劉 宇, 呂一河, 傅伯杰. 景觀格局-土壤侵蝕研究中景觀指數(shù)的意義解釋及局限性[J]. 生態(tài)學(xué)報(bào), 2011,30(1):267-275.

    [4] Borselli L, Cassi P, Torri D. Prolegomena to sediment and flow connectivity in the landscape: A GIS and field numerical assessment[J]. Catena, 2008,(75):268-277.

    [5] Imeson A C, Prinsen H A M. Vegetation patterns as biological indicators for identifying runoff and sediment source and sink areas for semi-arid landscapes in Spain[J]. Agriculture, Ecosystems & Environment, 2004,104:333-342.

    [6] Ludwig J A, Bastin G N, Chewings V H, et al. Leakiness: A new index for monitoring the health of arid and semiarid landscapes using remotely sensed vegetation cover and elevation data[J]. Ecological Indicators, 2007,(7):442-454.

    [7] Chen Liding, Wei Wei, Fu Bojie, et al. Soil and water conservation on the Loess Plateau in China: review and perspective[J]. Progress in Physical Geography, 2007,31:389-403.

    [8] 趙英時(shí). 遙感應(yīng)用分析原理與方法[M]. 北京:科學(xué)出版社,2003:391-392.

    [9] 陳 軍,周啟星. 數(shù)字地形分析[M].北京:科學(xué)出版社,2006:158-161.

    [10] 汪麗娜, 王 勇, 高 鵬,等. 黃土高原粗泥沙集中來(lái)源區(qū)水沙變化特征及趨勢(shì)性分析[J]. 水土保持通報(bào), 2008,28(2):11-16.

    [11] 錢云平, 蔣秀華, 金雙彥,等. 黃河中游黃土高原區(qū)河川基流特點(diǎn)及變化分析[J]. 地球科學(xué)與環(huán)境學(xué)報(bào),2004,26(2):88-91.

    猜你喜歡
    輸沙覆蓋度柵格
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測(cè)
    基于NDVI的晉州市植被覆蓋信息提取
    基于鄰域柵格篩選的點(diǎn)云邊緣點(diǎn)提取方法*
    低覆蓋度CO分子在Ni(110)面的吸附研究
    不同剖面形狀的柵格壁對(duì)柵格翼氣動(dòng)特性的影響
    基于CVT排布的非周期柵格密度加權(quán)陣設(shè)計(jì)
    湖南省四水流域森林、徑流、輸沙變化特性
    基于分布式水文模型的流域輸沙過程模擬
    塔克拉瑪干沙漠腹地塔中地區(qū)的風(fēng)動(dòng)力環(huán)境
    推移質(zhì)輸沙率計(jì)算公式分析
    亚洲色图综合在线观看| 国产又色又爽无遮挡免| 熟女少妇亚洲综合色aaa.| 丝袜人妻中文字幕| 亚洲国产av影院在线观看| 亚洲av男天堂| 亚洲国产精品一区三区| 国产精品熟女久久久久浪| 国产无遮挡羞羞视频在线观看| 精品国产超薄肉色丝袜足j| 国产乱人偷精品视频| kizo精华| 免费观看人在逋| 波多野结衣一区麻豆| www.精华液| 啦啦啦在线免费观看视频4| 又黄又粗又硬又大视频| 国产一级毛片在线| 国产成人一区二区在线| 国产黄频视频在线观看| 精品第一国产精品| 亚洲精品成人av观看孕妇| 国产成人系列免费观看| 国产av精品麻豆| 69精品国产乱码久久久| 天堂8中文在线网| 欧美日韩av久久| 亚洲精品中文字幕在线视频| 国产一卡二卡三卡精品 | 一区在线观看完整版| 国产激情久久老熟女| 婷婷色综合www| 久久人人97超碰香蕉20202| 亚洲成人国产一区在线观看 | 成人三级做爰电影| 亚洲精品av麻豆狂野| 国产精品一二三区在线看| 一区福利在线观看| 欧美精品一区二区大全| www.精华液| av不卡在线播放| 亚洲伊人色综图| 国产探花极品一区二区| 香蕉国产在线看| 国产片内射在线| 9热在线视频观看99| 午夜福利视频在线观看免费| 国产av精品麻豆| 日本一区二区免费在线视频| 日韩成人av中文字幕在线观看| 老司机在亚洲福利影院| 在线观看免费视频网站a站| 欧美xxⅹ黑人| 美女扒开内裤让男人捅视频| 97人妻天天添夜夜摸| 亚洲欧美色中文字幕在线| 成人亚洲欧美一区二区av| 操出白浆在线播放| 91精品国产国语对白视频| 婷婷色av中文字幕| 日韩,欧美,国产一区二区三区| 只有这里有精品99| 亚洲情色 制服丝袜| 97人妻天天添夜夜摸| 在线观看一区二区三区激情| 交换朋友夫妻互换小说| 美女扒开内裤让男人捅视频| 一二三四在线观看免费中文在| 久久国产精品男人的天堂亚洲| 亚洲成人av在线免费| 亚洲av成人精品一二三区| 免费观看性生交大片5| 麻豆乱淫一区二区| 99re6热这里在线精品视频| 亚洲精品久久久久久婷婷小说| 夫妻午夜视频| 久久人人爽人人片av| 巨乳人妻的诱惑在线观看| 国产有黄有色有爽视频| 久久女婷五月综合色啪小说| 亚洲美女黄色视频免费看| 亚洲精品国产色婷婷电影| netflix在线观看网站| 久久久久久久大尺度免费视频| 男女床上黄色一级片免费看| 男人操女人黄网站| 亚洲,一卡二卡三卡| av不卡在线播放| 亚洲熟女精品中文字幕| 成人漫画全彩无遮挡| 亚洲图色成人| 国产精品一区二区在线不卡| 女人久久www免费人成看片| 又黄又粗又硬又大视频| 伊人久久大香线蕉亚洲五| 国产一区二区三区综合在线观看| 亚洲欧美成人综合另类久久久| 国产99久久九九免费精品| 婷婷成人精品国产| www.av在线官网国产| 最黄视频免费看| 最近中文字幕高清免费大全6| 国产精品一区二区在线不卡| videos熟女内射| 亚洲美女黄色视频免费看| 新久久久久国产一级毛片| 一本色道久久久久久精品综合| 久久久国产精品麻豆| 老司机靠b影院| 精品卡一卡二卡四卡免费| 亚洲精品乱久久久久久| 麻豆av在线久日| 在线天堂中文资源库| 国产又爽黄色视频| 下体分泌物呈黄色| 97在线人人人人妻| 久久久久久久久久久久大奶| 久久午夜综合久久蜜桃| 老司机靠b影院| 伦理电影大哥的女人| 成人国产麻豆网| 啦啦啦视频在线资源免费观看| 伊人久久大香线蕉亚洲五| 黑丝袜美女国产一区| 十分钟在线观看高清视频www| 亚洲四区av| 精品一区二区三卡| 老鸭窝网址在线观看| 激情五月婷婷亚洲| 亚洲成人一二三区av| 在线天堂中文资源库| 亚洲av成人精品一二三区| 一区二区三区乱码不卡18| 无遮挡黄片免费观看| 天堂中文最新版在线下载| 大片电影免费在线观看免费| 亚洲五月色婷婷综合| 国产一卡二卡三卡精品 | 日本猛色少妇xxxxx猛交久久| 国产 精品1| 18禁国产床啪视频网站| 欧美日韩福利视频一区二区| 日韩免费高清中文字幕av| 一个人免费看片子| 麻豆av在线久日| 亚洲一卡2卡3卡4卡5卡精品中文| 日本欧美视频一区| 欧美久久黑人一区二区| 久久午夜综合久久蜜桃| av一本久久久久| 啦啦啦视频在线资源免费观看| a级毛片黄视频| 精品国产国语对白av| 人妻人人澡人人爽人人| 女人高潮潮喷娇喘18禁视频| 老鸭窝网址在线观看| 精品少妇久久久久久888优播| 亚洲av成人不卡在线观看播放网 | 久久狼人影院| 亚洲熟女精品中文字幕| 又黄又粗又硬又大视频| 十八禁网站网址无遮挡| 2018国产大陆天天弄谢| 熟妇人妻不卡中文字幕| 最近最新中文字幕大全免费视频 | 亚洲熟女毛片儿| 国产精品.久久久| 在线观看三级黄色| 亚洲国产成人一精品久久久| 中文天堂在线官网| av卡一久久| 久久狼人影院| 国产极品粉嫩免费观看在线| 最近2019中文字幕mv第一页| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图| 飞空精品影院首页| 日韩人妻精品一区2区三区| 看免费av毛片| 欧美变态另类bdsm刘玥| 精品一区二区免费观看| 啦啦啦啦在线视频资源| 精品免费久久久久久久清纯 | 999精品在线视频| 最近中文字幕2019免费版| 国产片内射在线| 丝瓜视频免费看黄片| 深夜精品福利| 亚洲欧洲国产日韩| 最近2019中文字幕mv第一页| 欧美激情高清一区二区三区 | 国产成人精品在线电影| av天堂久久9| 国产乱人偷精品视频| 亚洲七黄色美女视频| 国产97色在线日韩免费| 青青草视频在线视频观看| 少妇猛男粗大的猛烈进出视频| av卡一久久| 少妇 在线观看| 亚洲精品一二三| 最黄视频免费看| 不卡av一区二区三区| 黄片小视频在线播放| 精品久久久久久电影网| 久久久精品区二区三区| 国产有黄有色有爽视频| 人人澡人人妻人| 免费黄频网站在线观看国产| 大香蕉久久成人网| 青春草亚洲视频在线观看| 男男h啪啪无遮挡| 老司机靠b影院| 岛国毛片在线播放| 久久精品久久久久久噜噜老黄| 一本久久精品| 男人舔女人的私密视频| 尾随美女入室| 街头女战士在线观看网站| 国产人伦9x9x在线观看| 亚洲欧美一区二区三区国产| 日本91视频免费播放| 人人妻人人爽人人添夜夜欢视频| 男的添女的下面高潮视频| 国产亚洲欧美精品永久| 中文字幕精品免费在线观看视频| 建设人人有责人人尽责人人享有的| 午夜免费观看性视频| 精品国产一区二区久久| 亚洲精品国产区一区二| 男女免费视频国产| 99九九在线精品视频| 免费看不卡的av| 桃花免费在线播放| 久久久久久人妻| 麻豆av在线久日| 精品久久久精品久久久| 国产精品.久久久| 欧美精品av麻豆av| 精品一区二区三卡| 一区二区av电影网| 国产精品久久久久久久久免| 婷婷色麻豆天堂久久| 成人影院久久| 欧美人与性动交α欧美软件| 国产爽快片一区二区三区| 亚洲人成网站在线观看播放| 一级毛片黄色毛片免费观看视频| 欧美日韩福利视频一区二区| 丰满乱子伦码专区| 国产精品久久久久久精品电影小说| 亚洲伊人色综图| 亚洲av电影在线进入| 国产成人免费观看mmmm| av卡一久久| 伊人久久国产一区二区| 免费高清在线观看日韩| 少妇被粗大猛烈的视频| 精品少妇久久久久久888优播| 亚洲成人一二三区av| 亚洲成av片中文字幕在线观看| 超色免费av| 国产成人啪精品午夜网站| netflix在线观看网站| av线在线观看网站| 久久鲁丝午夜福利片| 亚洲精品久久久久久婷婷小说| 亚洲av电影在线观看一区二区三区| 女性被躁到高潮视频| 捣出白浆h1v1| www.熟女人妻精品国产| 成年动漫av网址| 日本vs欧美在线观看视频| 麻豆精品久久久久久蜜桃| 国产福利在线免费观看视频| 国语对白做爰xxxⅹ性视频网站| 国产探花极品一区二区| 久久久久久久久免费视频了| 国产一卡二卡三卡精品 | 免费看av在线观看网站| 亚洲激情五月婷婷啪啪| av女优亚洲男人天堂| 最近中文字幕2019免费版| 亚洲,一卡二卡三卡| a级片在线免费高清观看视频| 精品久久蜜臀av无| 国产乱人偷精品视频| 男人操女人黄网站| 精品少妇一区二区三区视频日本电影 | 老司机影院毛片| 亚洲精品久久午夜乱码| 一本一本久久a久久精品综合妖精| 日韩av免费高清视频| 精品亚洲成a人片在线观看| 哪个播放器可以免费观看大片| 在线免费观看不下载黄p国产| 51午夜福利影视在线观看| 色播在线永久视频| 亚洲在久久综合| 亚洲精品国产色婷婷电影| 亚洲国产看品久久| 成人三级做爰电影| 精品国产一区二区三区四区第35| 国产精品一国产av| 欧美国产精品va在线观看不卡| 国产爽快片一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产黄频视频在线观看| av国产精品久久久久影院| 欧美黑人精品巨大| 欧美国产精品一级二级三级| 免费少妇av软件| 99热网站在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲国产av影院在线观看| 亚洲国产成人一精品久久久| 操出白浆在线播放| 中文字幕亚洲精品专区| 日韩精品有码人妻一区| 久久久久久久国产电影| 国产成人精品无人区| 美女扒开内裤让男人捅视频| av国产久精品久网站免费入址| 国产爽快片一区二区三区| 欧美激情极品国产一区二区三区| 国产激情久久老熟女| 搡老岳熟女国产| 中国三级夫妇交换| 搡老岳熟女国产| 一本—道久久a久久精品蜜桃钙片| 久久精品国产a三级三级三级| 69精品国产乱码久久久| 纯流量卡能插随身wifi吗| 两性夫妻黄色片| 建设人人有责人人尽责人人享有的| 亚洲精品自拍成人| 精品少妇一区二区三区视频日本电影 | 欧美精品高潮呻吟av久久| 国产欧美日韩综合在线一区二区| 国产老妇伦熟女老妇高清| 久久久欧美国产精品| 亚洲国产精品一区二区三区在线| av一本久久久久| 麻豆乱淫一区二区| 精品人妻熟女毛片av久久网站| 成人手机av| 免费观看性生交大片5| 久久人人97超碰香蕉20202| 成年美女黄网站色视频大全免费| 国产日韩一区二区三区精品不卡| 男女午夜视频在线观看| 欧美日韩成人在线一区二区| 亚洲自偷自拍图片 自拍| 欧美黄色片欧美黄色片| 精品久久久久久电影网| av不卡在线播放| 欧美人与善性xxx| 日韩制服骚丝袜av| 国产老妇伦熟女老妇高清| 80岁老熟妇乱子伦牲交| 亚洲欧美成人综合另类久久久| 欧美精品亚洲一区二区| 九草在线视频观看| 亚洲欧洲日产国产| 伊人久久大香线蕉亚洲五| 亚洲精品自拍成人| 午夜91福利影院| 亚洲国产精品999| 可以免费在线观看a视频的电影网站 | 新久久久久国产一级毛片| 久久99精品国语久久久| 美国免费a级毛片| 欧美人与善性xxx| 大码成人一级视频| 男人舔女人的私密视频| 丝袜人妻中文字幕| 精品少妇内射三级| av线在线观看网站| 狠狠婷婷综合久久久久久88av| 国产高清国产精品国产三级| 久久久久久久国产电影| 狠狠婷婷综合久久久久久88av| 国产视频首页在线观看| 赤兔流量卡办理| 免费av中文字幕在线| 国产午夜精品一二区理论片| 青春草亚洲视频在线观看| 国产成人一区二区在线| 看免费av毛片| 老司机影院毛片| 人妻人人澡人人爽人人| 高清黄色对白视频在线免费看| 成年人午夜在线观看视频| 欧美精品人与动牲交sv欧美| 精品少妇一区二区三区视频日本电影 | 欧美日韩亚洲国产一区二区在线观看 | 自线自在国产av| 久久久精品国产亚洲av高清涩受| 久久久国产欧美日韩av| 午夜久久久在线观看| 另类亚洲欧美激情| 人人妻人人添人人爽欧美一区卜| 性高湖久久久久久久久免费观看| 一级毛片我不卡| 免费观看人在逋| 中文乱码字字幕精品一区二区三区| 久久久久久久久久久久大奶| 久久人人爽人人片av| 黄片无遮挡物在线观看| 伊人久久大香线蕉亚洲五| 日日撸夜夜添| 精品视频人人做人人爽| 成人18禁高潮啪啪吃奶动态图| 国产精品二区激情视频| 国产免费视频播放在线视频| 丝袜美腿诱惑在线| 国产精品久久久久久精品电影小说| 色精品久久人妻99蜜桃| 热re99久久精品国产66热6| 最黄视频免费看| 人人妻,人人澡人人爽秒播 | 综合色丁香网| 色94色欧美一区二区| 少妇猛男粗大的猛烈进出视频| 免费在线观看完整版高清| 亚洲国产av新网站| 亚洲欧美激情在线| 波野结衣二区三区在线| 国产精品欧美亚洲77777| 中文字幕制服av| 一区二区av电影网| 美女午夜性视频免费| 免费av中文字幕在线| 国产成人一区二区在线| 亚洲婷婷狠狠爱综合网| 在线看a的网站| 免费黄色在线免费观看| 最黄视频免费看| 男女午夜视频在线观看| 国产精品 欧美亚洲| 亚洲男人天堂网一区| 在线观看免费高清a一片| 欧美日韩综合久久久久久| 女性被躁到高潮视频| 又大又爽又粗| 18禁观看日本| 日韩中文字幕视频在线看片| 少妇被粗大的猛进出69影院| 久久97久久精品| 99国产综合亚洲精品| 亚洲欧美中文字幕日韩二区| 亚洲熟女精品中文字幕| 丁香六月天网| 欧美国产精品一级二级三级| 亚洲国产日韩一区二区| 亚洲精品国产av成人精品| 男女免费视频国产| 精品国产超薄肉色丝袜足j| 日韩不卡一区二区三区视频在线| 国产成人系列免费观看| 成人18禁高潮啪啪吃奶动态图| 岛国毛片在线播放| 免费日韩欧美在线观看| 超色免费av| 九九爱精品视频在线观看| xxx大片免费视频| 宅男免费午夜| 男男h啪啪无遮挡| svipshipincom国产片| 国产一卡二卡三卡精品 | 免费看不卡的av| 中文字幕av电影在线播放| 最近中文字幕2019免费版| 777久久人妻少妇嫩草av网站| 亚洲精品,欧美精品| 亚洲精品国产av蜜桃| 成人国语在线视频| 亚洲精华国产精华液的使用体验| 少妇猛男粗大的猛烈进出视频| 制服丝袜香蕉在线| 另类精品久久| 最新在线观看一区二区三区 | 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 制服诱惑二区| 国产一区有黄有色的免费视频| 久久ye,这里只有精品| 午夜精品国产一区二区电影| 久久久久国产一级毛片高清牌| 亚洲欧美成人精品一区二区| 亚洲国产欧美网| 伊人亚洲综合成人网| 国产成人欧美| 一本—道久久a久久精品蜜桃钙片| 久久毛片免费看一区二区三区| 啦啦啦在线免费观看视频4| 天堂中文最新版在线下载| 黄网站色视频无遮挡免费观看| 久久这里只有精品19| 秋霞伦理黄片| 在线观看免费高清a一片| 男女下面插进去视频免费观看| 下体分泌物呈黄色| 亚洲七黄色美女视频| 亚洲婷婷狠狠爱综合网| 热99久久久久精品小说推荐| 在线观看三级黄色| 成年人午夜在线观看视频| 女人久久www免费人成看片| 亚洲国产中文字幕在线视频| 国产亚洲av片在线观看秒播厂| 极品人妻少妇av视频| 老司机靠b影院| 国产精品无大码| 亚洲欧美日韩另类电影网站| av网站免费在线观看视频| 成年人午夜在线观看视频| 黑丝袜美女国产一区| 亚洲综合精品二区| 伦理电影大哥的女人| 丁香六月欧美| 看免费成人av毛片| 免费观看性生交大片5| 亚洲精品美女久久久久99蜜臀 | 黄色怎么调成土黄色| 亚洲国产精品999| 青春草亚洲视频在线观看| 午夜福利,免费看| 丰满迷人的少妇在线观看| 亚洲人成网站在线观看播放| av网站在线播放免费| tube8黄色片| 一级毛片 在线播放| 亚洲av成人不卡在线观看播放网 | 久久久欧美国产精品| 国产精品欧美亚洲77777| 免费少妇av软件| 成人三级做爰电影| 久久狼人影院| 777米奇影视久久| 久久久国产欧美日韩av| 日本wwww免费看| 国产在线一区二区三区精| 看非洲黑人一级黄片| 天天躁夜夜躁狠狠躁躁| 国产精品秋霞免费鲁丝片| 国产乱人偷精品视频| 最新在线观看一区二区三区 | 老司机靠b影院| 亚洲精品乱久久久久久| 捣出白浆h1v1| 欧美日韩精品网址| 国产一区二区在线观看av| 宅男免费午夜| 老司机深夜福利视频在线观看 | 1024视频免费在线观看| 亚洲少妇的诱惑av| 黄片播放在线免费| 亚洲一码二码三码区别大吗| 18禁观看日本| 国产人伦9x9x在线观看| 免费观看av网站的网址| 国产精品蜜桃在线观看| 男女边吃奶边做爰视频| 欧美日韩综合久久久久久| 精品国产超薄肉色丝袜足j| 777米奇影视久久| 亚洲中文av在线| 国产黄色视频一区二区在线观看| 国产成人91sexporn| 国产成人一区二区在线| 中文乱码字字幕精品一区二区三区| 国产欧美日韩综合在线一区二区| 欧美变态另类bdsm刘玥| 一本久久精品| 69精品国产乱码久久久| 国产av一区二区精品久久| 一区二区三区乱码不卡18| 亚洲熟女精品中文字幕| 秋霞伦理黄片| 一区二区三区乱码不卡18| 成人黄色视频免费在线看| 亚洲精品久久成人aⅴ小说| 午夜91福利影院| 欧美日韩视频精品一区| 午夜福利乱码中文字幕| 一级毛片我不卡| 亚洲天堂av无毛| 日本欧美视频一区| 中文字幕高清在线视频| 久久性视频一级片| 亚洲国产看品久久| 精品国产国语对白av| 国产精品一国产av| 一区二区日韩欧美中文字幕| netflix在线观看网站| 亚洲色图 男人天堂 中文字幕| avwww免费| 亚洲欧美一区二区三区久久| 成人黄色视频免费在线看| 久久久精品免费免费高清| 亚洲精品国产av成人精品| 99精品久久久久人妻精品| 蜜桃在线观看..| 热99国产精品久久久久久7| 99热国产这里只有精品6| kizo精华| 男人爽女人下面视频在线观看| 久久这里只有精品19| 中文字幕人妻丝袜一区二区 | 国产极品粉嫩免费观看在线| a级片在线免费高清观看视频| 国产片特级美女逼逼视频| 菩萨蛮人人尽说江南好唐韦庄| 人妻一区二区av|