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

    六盤山華北落葉松林坡面土壤密度空間變異與合理取樣數(shù)

    2022-08-03 10:57:58鄧佳楠劉澤彬郭建斌王彥輝
    林業(yè)科學(xué)研究 2022年4期
    關(guān)鍵詞:坡面變異土層

    鄧佳楠,劉 帆,劉澤彬,郭建斌,王彥輝

    (1.北京林業(yè)大學(xué)水土保持學(xué)院,北京 100083;2.中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與自然保護(hù)研究所/國家林業(yè)和草原局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100091)

    土壤密度是土壤的基本物理性質(zhì)之一,對(duì)坡面土壤水分運(yùn)動(dòng)有重要影響[1]。受氣候(如降雨)、土壤質(zhì)地、地形和植被類型等差異的影響,土壤密度在坡面上有很大的空間異質(zhì)性。準(zhǔn)確探明坡面土壤密度空間變異及其影響因素對(duì)于深入理解坡面產(chǎn)流機(jī)制和土壤水分的坡面分配有重要意義[2]。

    土壤密度空間變異研究涉及到了區(qū)域、流域、坡面和樣地等不同的空間尺度[3-7]以及草地、灌叢和森林等不同的植被類型[8-11],但關(guān)于坡面土壤密度空間變異的研究,主要集中在草地[8]和灌叢[10]。現(xiàn)有的森林坡面研究多關(guān)注土壤密度的坡位差異[12],但通過幾個(gè)坡位的簡(jiǎn)單比較還不能完整揭示坡面微地形、立地環(huán)境等各差異下的土壤密度空間變異特征,這是因?yàn)榧词乖谔囟ㄆ挛坏耐寥烂芏纫泊嬗休^大空間變異[11]。因此,森林坡面土壤密度空間變異研究還需細(xì)化。影響土壤密度空間變異的因素在不同空間尺度存在差異,在坡面上,地形(如坡度)、植被類型和土壤質(zhì)地(結(jié)構(gòu))等是影響坡面土壤密度空間差異的重要因素[13]。但在不同地區(qū)、植被類型下,其主導(dǎo)因素也有差別,如影響土壤密度空間變異的主導(dǎo)因素在河北壩上地區(qū)草地坡面為地形[8],在四川西緣山地坡面是植被[14],在喀斯特坡面是地質(zhì)[10]。作為黃土高原的產(chǎn)流核心區(qū)-土石山林區(qū),土壤密度異質(zhì)性是影響坡面產(chǎn)流空間變化的重要因素,但坡面土壤密度空間變異的影響因素未有探索。另外,多數(shù)研究關(guān)注表層或特定深度土壤密度的空間變異[10,15-16],但受不同土層石礫含量、根系生物量等差異影響[17-18],不同土層土壤密度空間變異可能存有差異,但目前關(guān)于土壤密度空間變異隨土層深度變化還缺乏足夠探索。

    為克服坡面土壤密度空間變異的影響,往往需進(jìn)行高密度的取樣來確保坡面土壤密度的評(píng)估精度,但費(fèi)時(shí)、費(fèi)力[19]。若能量化土壤密度合理取樣數(shù)及與空間變異的關(guān)系,將顯著提高坡面土壤取樣的效率和評(píng)估精度。目前有關(guān)土壤密度合理取樣數(shù)的研究較少,劉帆等[11]雖量化了華北落葉松林表層土壤密度的合理取樣數(shù),但這僅是樣地尺度的探索,且未探討土壤密度空間變異與其合理取樣數(shù)的關(guān)系。

    華北落葉松(Larix gmelinii var.principisrupprechtii (Mayr) Pilg.)是黃土高原土石山區(qū)主要的造林樹種,對(duì)當(dāng)?shù)氐乃春B(yǎng)、保持水土有重要作用。理清華北落葉松林坡面土壤密度空間變異及合理取樣數(shù)對(duì)于準(zhǔn)確理解華北落葉松林坡面產(chǎn)流空間差異的影響機(jī)制具有重要價(jià)值。為此,本研究以六盤山半濕潤區(qū)的華北落葉松林坡面為研究對(duì)象,結(jié)合經(jīng)典統(tǒng)計(jì)學(xué)和地統(tǒng)計(jì)學(xué)方法,明確華北落葉松林坡面0~100 cm土層土壤密度的空間異質(zhì)性及影響因素,并采用Monte Carlo模擬重抽樣方法,量化各層土壤密度的合理取樣數(shù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于黃土高原土石山區(qū)的六盤山香水河小流域(106°12′10.6″~106°16′30.5″ E,35°27 ′22.5 ″~35°33 ′29.7 ″ N),流域面積為43.7 km2,海拔2 010~2 942 m,土壤類型以山地灰褐土為主,石礫含量較高;氣候?yàn)榕瘻貛О霛駶櫄夂颍昃鶜鉁?.0 ℃,年均降水量618 mm。小流域內(nèi)主要以人工林和天然次生林為主,其中華北落葉松是主要的人工林樹種,占人工林總面積的90%。

    在小流域內(nèi)選取了1個(gè)東南坡向的典型華北落葉松人工純林坡面,坡面上華北落葉松均為同齡(39 a,2020)林。海拔范圍2 259~2 478 m,坡長(zhǎng)480.6 m,平均坡度為27.8°,土層深度為1.0~1.2 m左右,土壤總孔隙度55.2%,石礫含量9.5%。林分密度為822株·hm?2,郁閉度0.73,平均樹高17.1 m,平均胸徑19.8 cm,生長(zhǎng)季平均葉面積指數(shù)2.90。林下灌木零星分布,以秦嶺小檗(Berberis circumserrata (Schneid.) Schneid.)、蒙古莢蒾(Viburnum mongolicum (Pall.) Rehd.)和西北栒子(Cotoneaster zabelii Schneid.)為主。草本覆蓋度在40%左右,以東方草莓(Fragaria orientalis Lozinsk.)和點(diǎn)葉薹草(Carex hancokiana Maxim.)為主。

    1.2 土壤密度、孔隙度和石礫含量測(cè)定

    2020年9月,參考劉春利等[20]關(guān)于坡面土壤水文物理性質(zhì)空間變異的布點(diǎn)方法,沿坡面設(shè)立3條間距為15 m的樣線,在每條樣線自下而上每隔30 m設(shè)置1個(gè)取樣點(diǎn),共48個(gè)。在每個(gè)樣點(diǎn)按 0~20、20~40、40~60、60~80和80~100 cm 5個(gè)土層,用體積為200 cm3的環(huán)刀在各層的中部豎直向下取原狀土,并記錄取樣點(diǎn)坐標(biāo)(以最左側(cè)樣線下方第一個(gè)取樣點(diǎn)為原點(diǎn))。將樣品帶回室內(nèi)后,測(cè)定土壤密度、總孔隙度等,用排水法測(cè)定石礫含量。

    1.3 空間變異分析

    變異系數(shù)(CV)是經(jīng)典統(tǒng)計(jì)學(xué)中分析空間變異程度的衡量參數(shù)。根據(jù)變異程度強(qiáng)弱,CV可劃分為3個(gè)等級(jí):弱變異(CV≤0.1)、中等變異(0.1

    半方差函數(shù)是地統(tǒng)計(jì)學(xué)中分析空間結(jié)構(gòu)的主要工具。半方差函數(shù)的主要參數(shù)為塊金值(C0)、基臺(tái)值(C0+C)和變程(a)。塊金值C0表征隨機(jī)變異程度;基臺(tái)值(C0+C)為半變異函數(shù)的最大值,該值越大,表明總的空間異質(zhì)性越高[10];塊基比[C0/(C0+C)]用于判斷變量空間自相關(guān)程度,塊基比<0.25,表明變量呈強(qiáng)空間自相關(guān)性;0.25~0.75,表明變量呈中等空間自相關(guān)性;> 0.75,表明變量呈弱空間自相關(guān)性[21]。變程(a)是變量空間自相關(guān)范圍大小的參數(shù),在變程范圍內(nèi),兩點(diǎn)空間距離越小,空間自相關(guān)性越大;當(dāng)兩點(diǎn)空間距離大于變程后,兩點(diǎn)的空間自相關(guān)性不存在。

    1.4 合理取樣數(shù)確定

    利用Monte Carlo隨機(jī)抽樣方法,確定在一定誤差下坡面土壤密度的最少取樣數(shù)。首先,分別將 0~20、20~40、40~60、60~80和80~100 cm土層的土壤密度測(cè)值作為Monte Carlo隨機(jī)抽樣的初始值;然后對(duì)每個(gè)土層48個(gè)土壤密度測(cè)值進(jìn)行不放回抽樣(抽樣數(shù)n=2~47),為了保證充足的隨機(jī)樣本,每個(gè)抽樣數(shù)下(如n=2)重復(fù)1 000次,計(jì)算不同抽樣數(shù)下土壤密度均值;最后,計(jì)算并繪制不同抽樣數(shù)土壤密度的平均值和置信區(qū)間(95%)。根據(jù)土壤密度均值及其置信區(qū)間隨抽樣數(shù)的變化,確定在10%誤差內(nèi)各土層土壤密度的最少取樣數(shù)。

    1.5 數(shù)據(jù)分析

    半方差函數(shù)分析在GS+9.0軟件中進(jìn)行,土壤密度Kriging插值采用Arc GIS 10.3軟件,Monte Carlo模擬抽樣在R 3.6.1軟件中進(jìn)行。

    2 結(jié)果分析

    2.1 坡面土壤密度基本統(tǒng)計(jì)特征

    由表1可知,隨土層加深,土壤密度逐漸增加,0~20、20~40、40~60、60~80和80~100 cm土層土壤密度平均值(變化范圍)分別為0.88(0.73~1.08)、1.04(0.78~1.55)、1.23(0.91~1.66)、1.41(1.02~1.78)和1.49(1.09~1.73)g·cm?3,變異系數(shù)為0.09、0.16、0.17、0.15和0.11,土壤密度在0~20 cm土層為弱變異,在其它土層為中等變異。由圖1可知,土壤密度空間變異系數(shù)與石礫含量的空間變異系數(shù)呈顯著的(P<0.05)線性相關(guān),表明土壤密度空間變異隨土層深度的變化主要取決于石礫含量的空間異質(zhì)性。

    圖1 土壤密度變異系數(shù)與孔隙度和石礫含量變異系數(shù)的關(guān)系Fig.1 Relationships between the coefficient of variation (CV) of soil bulk density and the CVs of soil total porosity and rock fragment

    表1 坡面土壤密度的統(tǒng)計(jì)特征Table 1 Descriptive statistics characteristics of soil bulk density on the slope

    2.2 坡面土壤密度的空間結(jié)構(gòu)特征

    2.2.1 土壤密度的半方差函數(shù)分析 由表2可知,土壤密度的半方差函數(shù)在0~20和20~40 cm土層為高斯模型,在40~60和80~100 cm土層為球狀模型,在60~80 cm土層為指數(shù)模型。各層土壤密度塊金值均為正值;土壤密度基臺(tái)值沿土層深度的變化范圍0.011~0.060,在0~20 cm土層最小,在40~60 cm土層最大。0~20、20~40、40~60、60~80和80~100 cm土層的塊基比分別為0.295、0.348、0.300、0.495和0.498,表明各層土壤密度均屬中等空間自相關(guān)性。各層土壤密度變程分別為257.70、126.20、179.30、215.42和294.50 m,表明隨土層加深,空間自相關(guān)范圍先減小后增大,在80~100 cm土層最大,在20~40 cm土層最小。

    表2 坡面不同土層土壤密度的半方差函數(shù)理論模型及其相關(guān)參數(shù)Table 2 Semi-variogram theoretical models of soil bulk density in different soil layers and their related parameters

    2.2.2 坡面土壤密度的空間分布 由圖2可知,坡面各層土壤密度有明顯的空間異質(zhì)性,自坡頂至坡底,土壤密度總體上呈“低~高~低”的變化。各層土壤密度的最大值基本在坡中,最小值在坡下(0~20 cm除外,在坡上)。0~20 cm土層土壤密度跨度范圍較小,斑塊較大,空間分布較均勻;40~60 cm土層土壤密度跨度范圍較大,條帶密集,空間分布格局最復(fù)雜。隨土層加深,土壤密度空間分布從斑塊狀逐漸變?yōu)闂l帶狀。

    圖2 坡面 0~20(a)、20~40(b)、40~60(c)、60~80(d)和80~100(e)cm 土層土壤密度的空間分布Fig.2 Spatial distribution of soil bulk density in 0~20(a),20~40(b),40~60(c),60~80(d) and 80~100(e) cm soil layer on the slope

    2.3 坡面土壤密度的最少取樣數(shù)

    基于Monte Carlo模擬重抽樣的方法確定了在95%的置信區(qū)間下,若控制坡面土壤密度在10%誤差內(nèi)的最少取樣數(shù)。由圖3可知,隨取樣數(shù)增加,土壤密度的平均值和置信區(qū)間逐漸向內(nèi)收縮。在95%置信區(qū)間下,當(dāng)取樣數(shù)分別為4、10、10、8和5個(gè)時(shí),得到的0~20、20~40、40~60、60~80和80~100 cm土層土壤密度可控制在10%誤差內(nèi)。由圖4可知,土壤密度的最少取樣數(shù)與其變異系數(shù)呈極顯著的線性相關(guān)(P<0.01),即隨土壤密度空間變異系數(shù)的增加,最少取樣數(shù)逐漸增加。

    圖3 Monte Carlo 模擬的 0~20(a)、20~40(b)、40~60(c)、60~80(d)和80~100(e)cm土層土壤密度平均值和置信區(qū)間隨取樣數(shù)的變化Fig.3 Variations of the average and confidence interval of soil bulk density in 0~20(a),20~40(b),40~60(c),60~80(d) and 80~100(e) cm soil layer simulated by Monte Carlo as the sample size

    圖4 坡面土壤密度最小取樣數(shù)隨其變異系數(shù)的變化Fig.4 Variation in the minimum sample size of soil bulk density with its coefficient of variation

    3 討論

    研究坡面100 cm土層土壤密度均值為1.21 g·cm?3,低于全國100 cm土層土壤密度均值(1.32±0.21 g·cm?3)[3],這是因森林根系的作用增加了土壤的疏松性和通氣性,使土壤密度降低[22];另外,本研究區(qū)為土石山區(qū),石礫較多,土壤孔隙度較高,導(dǎo)致土壤密度偏低。通常,土壤密度隨土層深度增加而增加,本研究也得到同樣的規(guī)律,其原因主要有兩方面:一是土壤有機(jī)質(zhì)的影響,研究表明[23],土壤密度大小與土壤有機(jī)質(zhì)含量顯著負(fù)相關(guān),土層越深,土壤有機(jī)質(zhì)含量越低,土壤密度則會(huì)越大;二是土體重量壓縮的結(jié)果,土層越深,土體越緊實(shí),土壤密度越大。這與周武等[24]得出的研究結(jié)果一致。土壤密度的變異系數(shù)在0~20 cm土層為弱變異,其他土層為中等變異,可能與不同土層石礫含量的空間異質(zhì)性有關(guān)(圖1),0~20 cm土層土壤主要來自于腐殖質(zhì)和土壤有機(jī)質(zhì),石礫含量較少,土壤較為均勻,空間變異??;其他土層石礫含量較多,且分布不均,空間變異大。

    各層土壤密度的塊金值均較小,表明隨機(jī)因素對(duì)土壤密度空間變異影響較小。本研究坡面土壤密度的基臺(tái)值與樣地尺度[11]的基臺(tái)值大致相當(dāng),表明在中小尺度上土壤密度總的空間變異相近。土壤密度的總空間異質(zhì)性在0~20 cm土層最小(基臺(tái)值最?。?0~60 cm土層最大(基臺(tái)值最大),這與經(jīng)典統(tǒng)計(jì)學(xué)的研究結(jié)果(變異系數(shù))相吻合。各層土壤密度塊金比在0.25~0.75間,表明各層土壤密度均為中等空間自相關(guān),這與周卓麗等[8]在草地坡面研究得到的不同土層土壤密度的空間自相關(guān)性存在明顯差異的結(jié)論有所不同,可能是與不同植被坡面的地形(如坡度)、立地(孔隙度、石礫含量、機(jī)械組成等)、生物因素(如根系生物量等)等空間變異不同有關(guān)[25-26]。隨土層加深,土壤密度的變程先減后增,可能是因?yàn)橥翆又胁渴芨?、石礫含量等因素的綜合作用更強(qiáng),導(dǎo)致土壤密度的空間自相關(guān)范圍小。土壤密度沿坡空間分布在不同植被類型坡面存有較大差異,如楊寬等[27]得到福建高山毛竹林上坡位和下坡位土壤密度高于中坡位;高雪松等[14]得到四川西緣山地坡面土壤密度表現(xiàn)為上坡位>中坡位>下坡位。本研究中,土壤密度從坡頂至坡腳總體呈“低~高~低”的變化,在坡中最大,這種差異可能與不同坡面土壤質(zhì)地、地形、植被類型和生長(zhǎng)特征、人為干擾等差異有關(guān)。

    本研究發(fā)現(xiàn),土壤密度最小取樣數(shù)與其變異系數(shù)呈顯著相關(guān),表明土壤密度最小取樣數(shù)主要取決于其空間變異程度。劉帆等[11]研究表明,當(dāng)林地(樣地)土壤密度空間變異系數(shù)為0.11時(shí),最少取樣數(shù)為4個(gè);本研究中,當(dāng)坡面80~100 cm土層的土壤密度空間變異系數(shù)為0.11時(shí),最少取樣數(shù)為5個(gè)。這表明在同等空間變異下,空間尺度越大,所需的最少取樣數(shù)可能越多,但該結(jié)論未來還需更多數(shù)據(jù)進(jìn)行驗(yàn)證。由于各層土壤密度空間變異不同,導(dǎo)致不同土層土壤密度最少取樣數(shù)存在較大差異,在進(jìn)行野外采樣時(shí),應(yīng)基于取樣深度及取樣深度內(nèi)土壤密度的空間變異程度制定合理的取樣策略。本研究坡面,為保證100 cm土層土壤密度的評(píng)估精度,至少需布設(shè)10個(gè)以上的取樣點(diǎn)。

    4 結(jié)論

    在本研究的華北落葉松林坡面上,在0~20 cm土層土壤密度為弱變異,在其他土層均為中等變異,不同土層土壤密度空間變異的差異與石礫含量的空間變異密切相關(guān)。土壤密度的最少取樣量取決于其空間變異程度,當(dāng)坡面各層土壤密度空間變異系數(shù)在0.09~0.17時(shí),若要控制100 cm土層土壤密度在10%誤差內(nèi)(95%的置信區(qū)間下),至少需要布設(shè)10個(gè)以上的取樣點(diǎn)。

    猜你喜歡
    坡面變異土層
    土釘噴錨在不同土層的支護(hù)應(yīng)用及效果分析
    變異危機(jī)
    變異
    土層 村與人 下
    土層——伊當(dāng)灣志
    土層 沙與土 上
    沖積扇油氣管道坡面侵蝕災(zāi)害因子分析
    超音速流越過彎曲坡面的反問題
    面板堆石壩墊層施工及坡面防護(hù)
    變異的蚊子
    欧美+日韩+精品| 91久久精品国产一区二区成人| 成人无遮挡网站| 欧美丝袜亚洲另类| 欧美潮喷喷水| 亚洲欧美精品自产自拍| 国产一区二区三区av在线| 麻豆国产97在线/欧美| 18禁裸乳无遮挡免费网站照片| 18禁动态无遮挡网站| 日韩在线高清观看一区二区三区| 蜜臀久久99精品久久宅男| 国产片特级美女逼逼视频| 亚洲av不卡在线观看| 成人鲁丝片一二三区免费| 丝袜喷水一区| 亚洲第一区二区三区不卡| 日本wwww免费看| 国产高清有码在线观看视频| 欧美成人精品欧美一级黄| 亚洲国产日韩欧美精品在线观看| 日本五十路高清| 久久久久久国产a免费观看| 国内精品宾馆在线| 日本一本二区三区精品| 美女脱内裤让男人舔精品视频| 久久精品夜夜夜夜夜久久蜜豆| 在线观看美女被高潮喷水网站| 一级黄片播放器| 亚洲av中文字字幕乱码综合| 久久精品人妻少妇| 变态另类丝袜制服| 亚洲欧美中文字幕日韩二区| 国产激情偷乱视频一区二区| 日本免费一区二区三区高清不卡| 九色成人免费人妻av| 性色avwww在线观看| 欧美成人午夜免费资源| .国产精品久久| 成人国产麻豆网| av视频在线观看入口| 亚洲婷婷狠狠爱综合网| 久热久热在线精品观看| 成人美女网站在线观看视频| 少妇熟女aⅴ在线视频| 2022亚洲国产成人精品| or卡值多少钱| 美女高潮的动态| 在线播放无遮挡| 国产三级中文精品| 舔av片在线| 少妇熟女欧美另类| 国产精品99久久久久久久久| 日韩精品有码人妻一区| 国模一区二区三区四区视频| 最近最新中文字幕大全电影3| 老司机影院毛片| 精品一区二区免费观看| 久久午夜福利片| 久久精品国产亚洲av天美| 99久久九九国产精品国产免费| 成人三级黄色视频| 亚洲aⅴ乱码一区二区在线播放| 中文字幕亚洲精品专区| 老司机影院毛片| 亚洲伊人久久精品综合 | 一本久久精品| 99热精品在线国产| 国产伦精品一区二区三区四那| 日韩在线高清观看一区二区三区| 人妻夜夜爽99麻豆av| 蜜臀久久99精品久久宅男| 美女xxoo啪啪120秒动态图| 亚洲自偷自拍三级| 亚洲五月天丁香| 好男人视频免费观看在线| 成人毛片60女人毛片免费| 日韩欧美 国产精品| av播播在线观看一区| 色综合亚洲欧美另类图片| 国产精品,欧美在线| 最近最新中文字幕免费大全7| 国产视频首页在线观看| 亚洲一区高清亚洲精品| 波野结衣二区三区在线| 久久人人爽人人爽人人片va| 99热6这里只有精品| 欧美高清成人免费视频www| 久久久午夜欧美精品| 久久久精品欧美日韩精品| 久久99热6这里只有精品| 伦精品一区二区三区| 91av网一区二区| 两个人视频免费观看高清| 久久久久久久国产电影| 精品一区二区三区人妻视频| 日韩强制内射视频| 久久久国产成人精品二区| 日本熟妇午夜| 亚洲av.av天堂| 毛片女人毛片| 看黄色毛片网站| 只有这里有精品99| 国产成人a∨麻豆精品| 国产乱人偷精品视频| 高清午夜精品一区二区三区| 亚洲内射少妇av| 亚洲欧美一区二区三区国产| 精品少妇黑人巨大在线播放 | 狂野欧美白嫩少妇大欣赏| 亚洲国产日韩欧美精品在线观看| 黄片wwwwww| 国产三级中文精品| 精品不卡国产一区二区三区| 国产精品一区二区三区四区久久| 看非洲黑人一级黄片| 丰满人妻一区二区三区视频av| 99在线人妻在线中文字幕| 18禁动态无遮挡网站| 国产日韩欧美在线精品| 免费一级毛片在线播放高清视频| 欧美+日韩+精品| 熟女电影av网| 成人三级黄色视频| 五月玫瑰六月丁香| 嫩草影院新地址| 日本欧美国产在线视频| 日本黄色片子视频| 99久久成人亚洲精品观看| 看黄色毛片网站| 夫妻性生交免费视频一级片| 免费看a级黄色片| 久久久久精品久久久久真实原创| 亚洲av福利一区| 国产在视频线精品| 中国美白少妇内射xxxbb| av在线观看视频网站免费| 日韩大片免费观看网站 | 日本黄大片高清| 色吧在线观看| 成年版毛片免费区| 久久精品人妻少妇| 亚洲av成人av| 久久久久九九精品影院| 99久久成人亚洲精品观看| 插逼视频在线观看| 99久久中文字幕三级久久日本| 一个人观看的视频www高清免费观看| 久久精品久久久久久久性| 高清毛片免费看| 国产久久久一区二区三区| 国产在线一区二区三区精 | 永久网站在线| 久久国产乱子免费精品| 高清日韩中文字幕在线| 青春草亚洲视频在线观看| 人妻少妇偷人精品九色| 高清视频免费观看一区二区 | 国语对白做爰xxxⅹ性视频网站| 免费观看精品视频网站| av在线观看视频网站免费| 69av精品久久久久久| 三级国产精品片| 秋霞伦理黄片| 久久亚洲精品不卡| 人妻系列 视频| 色综合色国产| 欧美一区二区国产精品久久精品| 亚洲国产精品国产精品| 国产乱人视频| 精品免费久久久久久久清纯| 午夜福利在线在线| 毛片一级片免费看久久久久| 亚洲欧美成人精品一区二区| 色吧在线观看| 晚上一个人看的免费电影| 欧美三级亚洲精品| 青青草视频在线视频观看| 精品99又大又爽又粗少妇毛片| 欧美日韩综合久久久久久| 床上黄色一级片| 免费看av在线观看网站| 色网站视频免费| 成人特级av手机在线观看| 欧美极品一区二区三区四区| 91久久精品电影网| 色播亚洲综合网| 啦啦啦韩国在线观看视频| 欧美日韩在线观看h| 久久6这里有精品| 国产日韩欧美在线精品| 欧美bdsm另类| 好男人在线观看高清免费视频| 久久亚洲国产成人精品v| 日韩精品有码人妻一区| 青春草国产在线视频| 一个人看的www免费观看视频| 高清av免费在线| 欧美3d第一页| 午夜福利在线观看免费完整高清在| 美女xxoo啪啪120秒动态图| av黄色大香蕉| 日韩欧美精品v在线| 插阴视频在线观看视频| 日韩成人伦理影院| 床上黄色一级片| 中文字幕免费在线视频6| 欧美精品一区二区大全| 亚洲国产精品成人综合色| 日韩视频在线欧美| 欧美日韩在线观看h| 亚洲国产欧美在线一区| 看十八女毛片水多多多| 天天躁夜夜躁狠狠久久av| 欧美bdsm另类| 青青草视频在线视频观看| 成人性生交大片免费视频hd| 免费黄网站久久成人精品| 欧美日本视频| 麻豆久久精品国产亚洲av| 久久久久久久久大av| 久久午夜福利片| 中文字幕制服av| 大话2 男鬼变身卡| 亚洲婷婷狠狠爱综合网| 午夜视频国产福利| av又黄又爽大尺度在线免费看 | 日韩,欧美,国产一区二区三区 | 麻豆成人午夜福利视频| 中文字幕亚洲精品专区| www.av在线官网国产| 男人狂女人下面高潮的视频| 欧美高清性xxxxhd video| 天天一区二区日本电影三级| 两个人的视频大全免费| 日韩欧美三级三区| 蜜桃亚洲精品一区二区三区| 亚洲精品自拍成人| 午夜免费激情av| 高清视频免费观看一区二区 | 97热精品久久久久久| 日韩欧美在线乱码| 欧美精品国产亚洲| 亚洲精品一区蜜桃| 26uuu在线亚洲综合色| 国产精品久久久久久精品电影小说 | 精品人妻视频免费看| 日本欧美国产在线视频| 国产高清不卡午夜福利| 天美传媒精品一区二区| 亚洲欧美成人精品一区二区| 国产黄a三级三级三级人| 亚洲人成网站在线播| 日韩一本色道免费dvd| 六月丁香七月| 91狼人影院| 国产高清视频在线观看网站| 尤物成人国产欧美一区二区三区| 色噜噜av男人的天堂激情| 熟女人妻精品中文字幕| 成年女人看的毛片在线观看| 69人妻影院| 日韩一区二区视频免费看| 日本wwww免费看| 久久精品久久久久久久性| 亚洲精品aⅴ在线观看| 免费搜索国产男女视频| 亚洲,欧美,日韩| 亚洲欧洲日产国产| 国产色爽女视频免费观看| 看非洲黑人一级黄片| 国产中年淑女户外野战色| 国产熟女欧美一区二区| 一个人观看的视频www高清免费观看| 尾随美女入室| 久久精品综合一区二区三区| 联通29元200g的流量卡| 国产精品一区二区三区四区久久| 国产又黄又爽又无遮挡在线| 五月伊人婷婷丁香| 热99re8久久精品国产| 夜夜看夜夜爽夜夜摸| 少妇裸体淫交视频免费看高清| 日本三级黄在线观看| 日韩高清综合在线| 日韩一区二区三区影片| 精品久久久久久久人妻蜜臀av| 纵有疾风起免费观看全集完整版 | 看十八女毛片水多多多| 日韩三级伦理在线观看| 美女高潮的动态| ponron亚洲| 亚洲综合精品二区| 欧美日韩综合久久久久久| 久久6这里有精品| 亚洲无线观看免费| 日韩强制内射视频| 欧美高清成人免费视频www| 中文精品一卡2卡3卡4更新| 免费黄网站久久成人精品| 免费看a级黄色片| 国产黄色视频一区二区在线观看 | 韩国高清视频一区二区三区| 国产女主播在线喷水免费视频网站 | 国产亚洲av嫩草精品影院| 久久鲁丝午夜福利片| 成人漫画全彩无遮挡| 日韩一区二区视频免费看| 亚洲精品色激情综合| 亚洲国产精品久久男人天堂| 久久精品久久久久久噜噜老黄 | 久久精品久久久久久噜噜老黄 | 美女cb高潮喷水在线观看| 亚洲国产精品成人久久小说| 久久久精品大字幕| 日本色播在线视频| 国产亚洲av嫩草精品影院| 国产探花在线观看一区二区| 亚洲怡红院男人天堂| 99热这里只有是精品在线观看| 插逼视频在线观看| 99九九线精品视频在线观看视频| 天堂影院成人在线观看| 亚洲精品亚洲一区二区| 美女cb高潮喷水在线观看| 日韩国内少妇激情av| av播播在线观看一区| 免费观看的影片在线观看| 天堂av国产一区二区熟女人妻| 2021天堂中文幕一二区在线观| 色播亚洲综合网| 天堂影院成人在线观看| 亚洲精品乱码久久久久久按摩| 色综合站精品国产| 成人鲁丝片一二三区免费| 欧美区成人在线视频| 国产男人的电影天堂91| 国产精品熟女久久久久浪| 久久久久免费精品人妻一区二区| 国产精品人妻久久久影院| 国产精品久久久久久久久免| 亚洲成人av在线免费| 精品久久久久久电影网 | 插逼视频在线观看| 丰满乱子伦码专区| 高清在线视频一区二区三区 | 少妇裸体淫交视频免费看高清| 亚洲国产精品国产精品| 免费观看在线日韩| 国产一区二区在线av高清观看| 中文字幕精品亚洲无线码一区| 精品酒店卫生间| 亚洲乱码一区二区免费版| 你懂的网址亚洲精品在线观看 | 欧美最新免费一区二区三区| 亚洲国产最新在线播放| 免费观看在线日韩| 禁无遮挡网站| 精品国产露脸久久av麻豆 | 少妇人妻精品综合一区二区| 蜜臀久久99精品久久宅男| 久久精品国产亚洲网站| 日韩欧美精品v在线| 欧美丝袜亚洲另类| 国产一区二区三区av在线| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| 亚洲国产成人一精品久久久| 看非洲黑人一级黄片| 中文字幕制服av| 少妇丰满av| 午夜精品国产一区二区电影 | 久久久久久久久久久免费av| 日韩欧美精品免费久久| 寂寞人妻少妇视频99o| 午夜福利视频1000在线观看| 国产精品一区二区性色av| 日本欧美国产在线视频| 少妇高潮的动态图| 一二三四中文在线观看免费高清| 国产精品一区二区三区四区免费观看| 亚洲五月天丁香| 在线a可以看的网站| 亚洲婷婷狠狠爱综合网| 免费观看在线日韩| 国产精品一区www在线观看| 日韩欧美 国产精品| 水蜜桃什么品种好| 亚洲真实伦在线观看| 中文字幕精品亚洲无线码一区| 午夜福利在线在线| 极品教师在线视频| 亚洲怡红院男人天堂| 久久精品国产鲁丝片午夜精品| 九九久久精品国产亚洲av麻豆| 天堂影院成人在线观看| 中文精品一卡2卡3卡4更新| 成人一区二区视频在线观看| 久久精品综合一区二区三区| 国产伦在线观看视频一区| 国产黄色小视频在线观看| 内射极品少妇av片p| 岛国在线免费视频观看| av天堂中文字幕网| 99久国产av精品国产电影| 两个人的视频大全免费| 久久精品国产亚洲网站| 亚洲欧美精品专区久久| 麻豆av噜噜一区二区三区| 91aial.com中文字幕在线观看| 一级av片app| 国产一区二区在线观看日韩| av黄色大香蕉| 亚洲内射少妇av| 久久久久久久久中文| 久久久精品欧美日韩精品| 久久久久久久久久成人| 亚洲高清免费不卡视频| 日本黄大片高清| 亚洲欧美日韩高清专用| 国产精品国产三级专区第一集| 婷婷色av中文字幕| 精品久久久久久电影网 | 国产精品精品国产色婷婷| 欧美+日韩+精品| 91久久精品国产一区二区成人| 国产高清不卡午夜福利| 国产淫片久久久久久久久| 日韩一本色道免费dvd| 成年免费大片在线观看| 日韩精品有码人妻一区| 国产精品麻豆人妻色哟哟久久 | 国产成人aa在线观看| 日日啪夜夜撸| 极品教师在线视频| 色吧在线观看| 波多野结衣巨乳人妻| 男的添女的下面高潮视频| 精品一区二区三区视频在线| 日本五十路高清| 欧美成人免费av一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 久久这里只有精品中国| 国产免费又黄又爽又色| 亚洲av电影在线观看一区二区三区 | 国产精品一二三区在线看| 日韩中字成人| 久久久久久久亚洲中文字幕| 一级毛片电影观看 | 久久久精品94久久精品| 亚洲欧美日韩东京热| 欧美xxxx黑人xx丫x性爽| 大话2 男鬼变身卡| 国产精品一二三区在线看| 少妇丰满av| 真实男女啪啪啪动态图| or卡值多少钱| 国产白丝娇喘喷水9色精品| 小蜜桃在线观看免费完整版高清| 亚洲国产欧美在线一区| 在线播放无遮挡| 欧美激情在线99| 国产熟女欧美一区二区| 久99久视频精品免费| 日本五十路高清| 精品99又大又爽又粗少妇毛片| 亚洲一级一片aⅴ在线观看| 免费一级毛片在线播放高清视频| 69人妻影院| 中文亚洲av片在线观看爽| 99在线视频只有这里精品首页| 日产精品乱码卡一卡2卡三| 一卡2卡三卡四卡精品乱码亚洲| 美女国产视频在线观看| 欧美潮喷喷水| 日日干狠狠操夜夜爽| 亚洲中文字幕日韩| 久久精品国产亚洲网站| 国产伦精品一区二区三区视频9| 欧美97在线视频| 乱人视频在线观看| av天堂中文字幕网| 国产成人aa在线观看| 哪个播放器可以免费观看大片| 中文在线观看免费www的网站| 麻豆国产97在线/欧美| 久久精品国产亚洲av涩爱| 一级毛片我不卡| 91久久精品电影网| 国产精品蜜桃在线观看| 爱豆传媒免费全集在线观看| 91精品一卡2卡3卡4卡| 毛片一级片免费看久久久久| 久久久久久九九精品二区国产| 男人狂女人下面高潮的视频| 免费不卡的大黄色大毛片视频在线观看 | 国产午夜精品论理片| 国产精品久久视频播放| 久久久久久九九精品二区国产| 91精品国产九色| 在线免费观看的www视频| 成人午夜精彩视频在线观看| 91久久精品国产一区二区成人| 精品一区二区三区人妻视频| 国产一级毛片在线| 亚州av有码| 色5月婷婷丁香| 午夜老司机福利剧场| 日韩一本色道免费dvd| 中文字幕亚洲精品专区| 99在线人妻在线中文字幕| 伊人久久精品亚洲午夜| 午夜精品在线福利| 少妇人妻一区二区三区视频| 国产精品麻豆人妻色哟哟久久 | 婷婷色av中文字幕| 免费不卡的大黄色大毛片视频在线观看 | www.av在线官网国产| 在现免费观看毛片| 免费观看a级毛片全部| 亚洲av成人精品一区久久| 嫩草影院精品99| 在线a可以看的网站| 中文字幕精品亚洲无线码一区| 亚洲美女视频黄频| 少妇的逼水好多| 男女视频在线观看网站免费| 亚洲中文字幕日韩| 亚洲国产色片| 国产精品一区二区性色av| 91av网一区二区| 禁无遮挡网站| 国产一级毛片七仙女欲春2| 男人和女人高潮做爰伦理| 在线a可以看的网站| 久久午夜福利片| 久久精品91蜜桃| av国产久精品久网站免费入址| 在线观看66精品国产| 美女cb高潮喷水在线观看| 亚洲中文字幕一区二区三区有码在线看| 亚洲四区av| 国产午夜福利久久久久久| 久久久色成人| 亚洲精品aⅴ在线观看| 男人和女人高潮做爰伦理| 国内精品美女久久久久久| 一本一本综合久久| 欧美一区二区精品小视频在线| 少妇人妻精品综合一区二区| 国产三级在线视频| 免费看日本二区| 久久久久精品久久久久真实原创| 久久久久久久久久久丰满| 爱豆传媒免费全集在线观看| 啦啦啦韩国在线观看视频| 我的女老师完整版在线观看| av视频在线观看入口| 午夜免费男女啪啪视频观看| 国产高清国产精品国产三级 | 国产精品不卡视频一区二区| 国产精品人妻久久久影院| 天堂网av新在线| 中文精品一卡2卡3卡4更新| 建设人人有责人人尽责人人享有的 | 日韩三级伦理在线观看| 免费一级毛片在线播放高清视频| 校园人妻丝袜中文字幕| 人体艺术视频欧美日本| 亚洲最大成人手机在线| 欧美日韩一区二区视频在线观看视频在线 | 亚洲久久久久久中文字幕| 久久婷婷人人爽人人干人人爱| 久久久久久久亚洲中文字幕| 欧美一区二区国产精品久久精品| 欧美日本视频| 少妇猛男粗大的猛烈进出视频 | 亚洲丝袜综合中文字幕| 一边亲一边摸免费视频| 日韩一区二区视频免费看| 欧美一区二区亚洲| 免费看光身美女| 亚洲高清免费不卡视频| 欧美三级亚洲精品| 国产成人aa在线观看| 欧美日韩在线观看h| 久久久久久久久久久丰满| 精品午夜福利在线看| 赤兔流量卡办理| 日本爱情动作片www.在线观看| 人人妻人人澡人人爽人人夜夜 | 搡老妇女老女人老熟妇| 一个人看的www免费观看视频| 99久久精品国产国产毛片| 亚洲精品乱码久久久久久按摩| 久久精品国产99精品国产亚洲性色| 亚洲美女视频黄频| 国产免费又黄又爽又色| 嫩草影院精品99| 女人久久www免费人成看片 | 永久网站在线| 国产精品国产三级国产专区5o | 亚洲美女搞黄在线观看| 亚洲av福利一区| 在线播放无遮挡| 免费观看的影片在线观看| 少妇丰满av| 国产伦理片在线播放av一区| 麻豆成人午夜福利视频| www.色视频.com| 久久久久久久久久成人| 熟女人妻精品中文字幕| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 啦啦啦啦在线视频资源| 国产黄a三级三级三级人|