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

    地表覆蓋對紫色土坡耕地壤中流流動和硝態(tài)氮運(yùn)移過程的影響

    2018-10-29 08:18:52龍?zhí)煊?/span>劉金輝劉祥章曹懷亮田大浪
    水土保持研究 2018年6期
    關(guān)鍵詞:結(jié)皮耕作層硝態(tài)

    龍?zhí)煊? 劉金輝, 劉祥章, 曹懷亮, 田大浪

    (1.重慶大學(xué) 三峽庫區(qū)生態(tài)環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室, 重慶 400045;2.重慶大學(xué) 低碳綠色建筑國際聯(lián)合研究中心, 重慶 400045; 3.重慶大學(xué) 土木工程學(xué)院, 重慶 400045)

    地表覆蓋主要影響降雨過程中水分的下滲,從而影響地表徑流與壤中流的分布。地表覆蓋的影響因素有很多,包括表土結(jié)皮、稻草覆蓋、植被覆蓋和礫石覆蓋等。三峽庫區(qū)的土壤類型以紫色土為主,占庫區(qū)總面積的近40%,總耕地面積高達(dá)75%以上,且以坡耕地為主[1]。紫色土坡耕地土層淺薄,土壤田間持水力低、土壤疏松,土壤水分極易蓄滿,持續(xù)降雨時,水分下滲至土壤—基巖界面,因紫色泥頁巖滲透性弱,水分難以垂直入滲而沿坡向匯集形成壤中流,壤中流十分發(fā)育。在紫色土坡耕地中,氮素主要隨壤中流流失。與平地中氮素隨壤中流垂向一維運(yùn)移不同,紫色土坡耕地氮素隨壤中流的運(yùn)移是二維、甚至是三維的,因此,針對三峽庫區(qū)紫色土坡耕地壤中流流動和養(yǎng)分流失機(jī)理的特點(diǎn)開展相關(guān)研究具有重要的理論和實(shí)際意義。

    通常土壤在降雨或灌溉的作用下或頻繁的耕作狀況下都會結(jié)皮,土壤表層形成一層致密層[2-5]。而紫色土即使無降雨作用,仍能形成穩(wěn)定的結(jié)皮層并影響水分入滲[6-7]。地表稻草覆蓋具有能減少土壤侵蝕、抑制降雨時土壤結(jié)皮等優(yōu)點(diǎn),逐漸在紫色土坡耕地中采用。迄今為止,對表土結(jié)皮和稻草覆蓋對紫色土坡耕地壤中流動和養(yǎng)分運(yùn)移影響研究的學(xué)者和發(fā)表的論文為數(shù)不多[8-10],并且所采用的方法均是把坡耕地作為“黑箱”的試驗(yàn)研究,有關(guān)土壤中水分流動和養(yǎng)分運(yùn)移過程的數(shù)值模擬研究未見報(bào)道。為此,本研究采用數(shù)值模擬的方法,通過建立土壤水分和養(yǎng)分運(yùn)移的三維數(shù)學(xué)模型,借助HYDRUS軟件,針對表土結(jié)皮和稻草覆蓋這兩種較為常見的地表覆蓋狀況,對降雨情況下紫色土坡耕地水分和硝態(tài)氮隨壤中流的運(yùn)移過程進(jìn)行模擬,分析紫色土坡耕地水分和硝態(tài)氮的流失特性。

    1 材料與方法

    1.1 數(shù)學(xué)模型與模擬區(qū)域

    1.1.1 土壤水分和養(yǎng)分運(yùn)移控制方程組 不計(jì)空氣、溫度等對水分運(yùn)移的作用,忽略水分運(yùn)移對土壤結(jié)構(gòu)和性質(zhì)的影響,采用Richards方程描述土壤中的水分運(yùn)移過程,即

    (1)

    采用對流—彌散方程對硝態(tài)氮運(yùn)移過程進(jìn)行模擬,即

    (2)

    式中:c為硝態(tài)氮濃度(mg/cm3);Dxi為xi方向硝態(tài)氮的彌散系數(shù)(cm2/min);ui為xi方向上的水流速度(cm/min);s為源項(xiàng),包括根系吸收、土壤吸附與解析等,本文中忽略不計(jì)。

    1.1.2 模擬區(qū)域 模擬的坡耕地的幾何尺寸如圖1所示。坡耕地長600 cm、寬300 cm、高50 cm,坡度分別為5°,10°,15°,20°和25°。土壤耕作不可避免地使土壤分層,因此,為了使模型更符合實(shí)際土層情況,把坡耕地分為耕作層和非耕作層。地表結(jié)皮和稻草覆蓋層分別作為在耕作層上的單獨(dú)一層。耕作層高度為20 cm,非耕作層高度為30 cm,地表結(jié)皮層或稻草覆蓋層的高度均為1 cm[8]。

    圖1 模擬區(qū)域和邊界條件示意

    1.1.3 模型參數(shù) 土壤水分特征曲線和非飽和土壤水力傳導(dǎo)率采用van Geneuchten 模型來表示。對結(jié)皮土壤的水力特性參數(shù)采用Chahinian N的確定方法和相關(guān)研究結(jié)果[11-12];稻草覆蓋層的水力特性參數(shù)的確定采用虎膽·吐馬爾白等的研究結(jié)果[13];耕作層和非耕作層的水力特性參數(shù)采用HYDRUS自帶的神經(jīng)元網(wǎng)絡(luò)反演法,由已有的試驗(yàn)數(shù)據(jù)反演優(yōu)化得出[14-15](表1)。

    表1 土壤水力特性參數(shù)

    溶質(zhì)運(yùn)移參數(shù)同樣采用優(yōu)化反演法確定。耕作層和非耕作層土壤體積密度為1.5 g/cm3,縱向彌散度為0.5 cm,橫向彌散度為0.01 cm,硝態(tài)氮在水中的分子擴(kuò)散系數(shù)為0.001 5 cm2/min,在土壤空氣中的分子擴(kuò)散系數(shù)為0。

    1.1.4 區(qū)域離散與網(wǎng)格劃分 借助HYDRUS軟件對數(shù)學(xué)模型進(jìn)行數(shù)值計(jì)算。采用等間距三角形網(wǎng)格,經(jīng)網(wǎng)格無關(guān)性檢驗(yàn),表土結(jié)皮情況下X方向設(shè)置61個節(jié)點(diǎn),Y方向31個,Z方向22個,共生成41 620個節(jié)點(diǎn)。第一層的高度為1 cm,第二層為1.5 cm,其他均為2.5 cm。稻草覆蓋情況下,X方向設(shè)置61個節(jié)點(diǎn),Y方向31個,Z方向21個,共生成39 711個節(jié)點(diǎn)。第一層的高度為1 cm,其他均為2.5 cm。模擬均采用隱式模式。

    1.2 定解條件

    降雨條件為降雨強(qiáng)度為30 mm/h,降雨歷時為3 h,該降雨條件屬三峽庫區(qū)典型大雨。分別針對坡度為5°,10°,15°,20°和25°,模擬坡耕地中水分和硝態(tài)氮濃度的時空分布,模擬時間均為4 h。不同表土性質(zhì)下的土壤中氮素均勻分布、含量相同,為150 kg/hm2。根據(jù)相關(guān)研究[16],稻草覆蓋能夠完全抑制降雨時土壤結(jié)皮的發(fā)生,因此本文不考慮表土結(jié)皮和稻草覆蓋的交互作用。為對照分析,在相同降雨特性和坡度的情況下,對不結(jié)皮、無稻草覆蓋時的坡耕地中水分和硝態(tài)氮濃度的時空分布也進(jìn)行了模擬。

    初始條件給定初始壓力水頭,設(shè)置地下水位線為-200 cm。施肥方式采用全層施肥,即肥料與表層0—5 cm的土壤均勻混合,由此得到初始氮素濃度0.3 mg/cm3。模型邊界條件如圖1所示,土壤表面為降雨條件下的大氣邊界,坡下為自由出流邊界,其他為零通量邊界。

    1.3 模型驗(yàn)證

    1.3.1 模型驗(yàn)證的基本條件 為率定模型參數(shù)和驗(yàn)證模型模擬結(jié)果,構(gòu)建如圖2所示的試驗(yàn)裝置并進(jìn)行相關(guān)試驗(yàn)[17]。如圖2 所示,土槽尺寸為1.0 m×0.5 m×0.5 m(長×寬×高),坡度可在0~30°內(nèi)任意調(diào)節(jié)。在裝置內(nèi)填裝紫色土,土層分為兩層,其中耕作層厚度20 cm,非耕作層厚度30 cm。試驗(yàn)坡度為15°,降雨強(qiáng)度為0.89 mm/min,降雨歷時62 min,試驗(yàn)重復(fù)進(jìn)行3次,取其平均,模擬施肥選用KNO3,不考慮根系影響。

    驗(yàn)證模型基本條件與室內(nèi)試驗(yàn)保持一致,尺寸設(shè)為100 cm×500 cm×500 cm(長×寬×高)。模型邊界條件為:ABCD面為大氣邊界條件,BFGC面為自由出流邊界條件,其他邊界條件為0通量邊界(圖1)。在BFGC面距地表20 cm處設(shè)一個觀察點(diǎn),用以檢測該點(diǎn)的硝態(tài)氮濃度。

    1.3.2 模型可靠性分析 把模型的模擬值與坡度15°,降雨強(qiáng)度為0.89 mm/min時的實(shí)測值進(jìn)行對比,該實(shí)測值沒有參與參數(shù)率定,結(jié)果如圖3、圖4所示。對觀察點(diǎn)硝態(tài)氮濃度,累積硝態(tài)氮流失量和實(shí)測值進(jìn)行兩配對樣本T 檢驗(yàn)驗(yàn)證模型的可靠性。如表2所示,結(jié)果表明:土壤硝態(tài)氮濃度和累積硝態(tài)氮出流量配對T 檢驗(yàn)的顯著性水平p值均不在置信區(qū)間(α=0.05),說明土壤硝態(tài)氮濃度和累積硝態(tài)氮出流量的模擬值與實(shí)測值均無顯著差異,模擬結(jié)果合理,參數(shù)較為可靠,構(gòu)建的模型可用于紫色土坡耕地的模擬與分析。

    圖2 試驗(yàn)裝置示意

    圖3 流失硝態(tài)氮濃度實(shí)測值與模擬值對比

    圖4 硝態(tài)氮出流量實(shí)測值與模擬值對比

    均值標(biāo)準(zhǔn)誤差T檢驗(yàn)概率P相關(guān)系數(shù)硝態(tài)氮濃度0.930.100.991累積硝態(tài)氮量0.00210.170.994

    2 結(jié)果與分析

    2.1 地表覆蓋對紫色土壤中流的影響

    2.1.1 不同地表覆蓋下的水分入滲 壤中流的產(chǎn)生和發(fā)展受土壤入滲的影響,分析土壤水分入滲規(guī)律有助于對壤中流的理解。圖5為三峽庫區(qū)典型坡度15°情況下,不同地表覆蓋下土壤入滲率隨時間的變化過程。從圖5中可以看到,在降雨初期,3種不同表土性質(zhì)下的土壤地表入滲率均為0.05 cm/min,等于降雨強(qiáng)度,此時降雨全部入滲。降雨大約25 min后,結(jié)皮的表土入滲率開始下降。不結(jié)皮、無稻草覆蓋對照情況下的土壤表土入滲率開始減小的時間約為60 min,比結(jié)皮情況下大為延后,土壤結(jié)皮使土壤入滲能力減弱。稻草覆蓋入滲率下降的時間更晚,約在降雨70 min后。在整個降雨期內(nèi),表土結(jié)皮時土壤表層平均入滲率為0.027 cm/min,無結(jié)皮無覆蓋時為0.040 cm/min,稻草覆蓋時為0.042 cm/min。表土結(jié)皮在很大程度上降低了土壤的入滲率,而稻草覆蓋略增加了入滲率。

    圖6為坡度15°時不同地表覆蓋下土壤濕潤鋒隨時間的變化。從圖6中可以看到,在降雨初期,3種表土情況下濕潤鋒運(yùn)移率(濕潤鋒/時間)均比較快,稻草覆蓋下濕潤鋒運(yùn)移率略低于對照情況,這是由于土壤表層有稻草層對降雨進(jìn)行截留和緩沖。經(jīng)過緩沖以后,稻草覆蓋下濕潤鋒運(yùn)移率基本與對照情況相當(dāng)。結(jié)皮下濕潤鋒運(yùn)移率明顯低于其他兩種情況,說明結(jié)皮阻礙了水分下滲;降雨約15 min后,3種表土狀況下的濕潤鋒均到達(dá)耕作層底部(20 cm處),因非耕作層滲透系數(shù)較小,此后濕潤鋒運(yùn)移率均明顯減慢。在降雨結(jié)束后,稻草覆蓋和對照土壤濕潤鋒均運(yùn)移到45 cm處,而結(jié)皮濕潤鋒僅運(yùn)移到30 cm處,表土結(jié)皮導(dǎo)致入滲量減少。

    圖5 不同地表覆蓋下表土水分入滲速率變化

    圖6 不同地表覆蓋下濕潤鋒隨時間變化

    2.1.2 不同地表覆蓋下壤中流的產(chǎn)生和發(fā)展 壤中流是指降雨、灌溉過程中及其隨后一段時間內(nèi),水分在土壤內(nèi)部的運(yùn)動,包括水分在土壤內(nèi)的下滲和側(cè)流。側(cè)流是坡地徑流的重要組成部分,對土壤水分和養(yǎng)分流失等具有重要影響[18]。圖7為坡度15°時3種不同地表狀況下土壤中水分隨時間的運(yùn)移過程。由圖可以看出,降雨初期,耕作層中水分的運(yùn)移基本上為垂直入滲。隨著降雨的持續(xù),耕作層逐漸從上到下在水分飽和區(qū)產(chǎn)生側(cè)流。在約90 min后,稻草覆蓋下側(cè)流基本上擴(kuò)展到整個耕作層。而在表土結(jié)皮和對照條件下在降雨進(jìn)行到190 min時,側(cè)流才出現(xiàn)在整個耕作層。從圖中3種不同表土性質(zhì)下的速度矢量的方向可以看出,當(dāng)水分進(jìn)入非耕作層后,在耕作層和非耕作層流速方向明顯不同,在耕作層的水分為側(cè)流,在非耕作層的水分主要為垂直入滲,這種情況持續(xù)到降雨結(jié)束。模擬結(jié)果表明,與對照情況相比,稻草覆蓋增加了水分入滲,使較多水分入滲到非耕作層,而表土結(jié)皮的情況正好相反。說明坡耕地中水分的流失過程主要發(fā)生在耕作層。

    圖7 表土結(jié)皮、稻草覆蓋、對照條件下各典型時刻水分運(yùn)移速度矢量

    2.1.3 表土性質(zhì)對壤中流出流量的影響 圖8為3種不同地表覆蓋下不同坡度在降雨過程中形成的累積的壤中流出流量。不同坡度下,結(jié)皮組的壤中流累積出流量為12 997~14 141 cm3,稻草覆蓋組的壤中流累積出流量為24 100~25 200 cm3,對照組的壤中流累積出流量為21 400~23 200 cm3。稻草覆蓋的壤中流出流量>對照>表土結(jié)皮,稻草覆蓋會增加壤中流出流量,而表土結(jié)皮會明顯降低壤中流出流量。稻草覆蓋時,稻草層的孔隙率大,滲透能力強(qiáng),降雨能夠更多地滲入土壤,導(dǎo)致更多的水分以壤中流的方式流出。與此相反,表土結(jié)皮時,土壤滲透性差,更多的降雨隨地表徑流流失,減少了入滲,從而減少壤中流流量。在5種坡度下,坡度由小到大,結(jié)皮情況下壤中流出流量分別比未結(jié)皮降低41.6%,43.0%,36.9%,40.5%和34.8%。坡度不同,表土結(jié)皮對壤中流的影響也不相同,在25°時土壤結(jié)皮對壤中流影響最小,10°時土壤結(jié)皮對壤中流影響最大。5種坡度下,稻草覆蓋比對照情況下的壤中流出流量分別增加了8.6%,10.1%,10.7%,11.9%和12.6%。說明坡度越大,稻草覆蓋對壤中流的影響越明顯,由于稻草覆蓋削弱了坡度對地表徑流的增加效果,使得坡度增加更有利于壤中流的產(chǎn)生和出流。

    圖8 不同地表覆蓋下壤中流累積出流量

    為了進(jìn)一步地研究不同表土性質(zhì)下壤中流出流的動態(tài)變化,圖9列出了坡度15°下,3種不同表土性質(zhì)下出流邊界壤中流速率變化對比。從圖中看出,表土結(jié)皮和稻草覆蓋下壤中流變化趨勢仍都遵循紫色土壤中流流失規(guī)律:分為降雨開始時的出流強(qiáng)度增加階段,穩(wěn)定或平穩(wěn)增加階段和降雨結(jié)束后的衰減階段。稻草覆蓋下,壤中流穩(wěn)定強(qiáng)度約為200 cm3/min,對比無稻草覆蓋下增加了約11.1%;而表土結(jié)皮下壤中流穩(wěn)定強(qiáng)度為95 cm3/min,對比無結(jié)皮情況下紫色土飽和穩(wěn)定強(qiáng)度約降低47.2%。在整個模擬期內(nèi),稻草覆蓋明顯增加壤中流強(qiáng)度而表土結(jié)皮均明顯降低壤中流強(qiáng)度,可能原因是稻草覆蓋和表土結(jié)皮后,土壤表層入滲能力變化,導(dǎo)致土壤內(nèi)部水壓力發(fā)生改變。

    2.2 表土性質(zhì)對紫色土硝態(tài)氮隨壤中流流失的影響

    圖10為3種表土性質(zhì)下硝態(tài)氮隨壤中流的流失量。在5個坡度下,表土結(jié)皮時硝態(tài)氮流失量分別減少了14.2%,15.9%,8.6%,34.4%和18.9%,稻草覆蓋時硝態(tài)氮流失量分別增加了8.6%,15.1%,15.4%,13.5%和11.1%。與表土結(jié)皮和稻草覆蓋對壤中流影響對比分析可以發(fā)現(xiàn),表土結(jié)皮下硝態(tài)氮隨壤中流流失的減小幅度明顯小于表土結(jié)皮對水分的影響,而稻草覆蓋下硝態(tài)氮隨壤中流流失的增加幅度卻略大于稻草覆蓋對水分的影響。與壤中流水分流失類似,不同坡度下,表土結(jié)皮對硝態(tài)氮流失量的影響程度也不相同,在坡度為15°時,表土結(jié)皮對硝態(tài)氮流失影響小,坡度為20°時影響最大,這與表土結(jié)皮對壤中流影響不同步。稻草覆蓋對硝態(tài)氮在15°時影響最大。

    圖9 不同地表覆蓋下邊界壤中流強(qiáng)度變化

    圖10 不同地表覆蓋下硝態(tài)氮隨壤中流流失量

    為了進(jìn)一步研究表土結(jié)皮下出流硝態(tài)氮濃度變化,選取典型坡度15°時不同表土性質(zhì)下出流硝態(tài)氮濃度隨時間變化進(jìn)行分析。如圖11所示,在降雨剛開始時,表土結(jié)皮下硝態(tài)氮濃度比未結(jié)皮時增長速率慢,到達(dá)濃度峰值的時間約推遲6 min,但表土結(jié)皮下硝態(tài)氮濃度峰值卻較未結(jié)皮時高,且?guī)缀跽麄€模擬期內(nèi)表土結(jié)皮時硝態(tài)氮濃度均高于土壤未結(jié)皮時;對比稻草覆蓋和對照情況的硝態(tài)氮濃度也可以發(fā)現(xiàn),稻草覆蓋下硝態(tài)氮濃度高于對照情況。這種現(xiàn)象的原因可能是土壤結(jié)皮和稻草覆蓋時,土壤內(nèi)部的水分入滲速率比對照情況時較慢,水分在土壤內(nèi)部停留時間較長,土壤中硝態(tài)氮溶解更為充分,單位水分溶解的硝態(tài)氮更多,這也解釋了土壤結(jié)皮對硝態(tài)氮隨壤中流流失的影響小于對水分流失的影響,而稻草覆蓋下硝態(tài)氮隨壤中流流失的增加幅度大于稻草覆蓋對水分的影響。

    圖11 不同地表覆蓋下出流硝態(tài)氮濃度隨時間變化

    3 結(jié) 論

    (1) 紫色土表土結(jié)皮會減少土壤表面水分的入滲,壤中流出流量降低34.8%~43.0%,壤中流穩(wěn)定出流強(qiáng)度約降低42.7%,硝態(tài)氮隨壤中流流失量減少8.6%~34.4%,但是流失的硝態(tài)氮濃度有所增加。稻草覆蓋下,壤中流出流量增加8.6%~12.6%,壤中流出流強(qiáng)度約增加11.1%,硝態(tài)氮隨壤中流流失量增加8.6%~15.4%,硝態(tài)氮流失濃度也有所增加。

    (2) 坡度對表土結(jié)皮下水分和硝態(tài)氮運(yùn)移的影響比較復(fù)雜,坡度不同,表土結(jié)皮對紫色土壤中流以及硝態(tài)氮流失的影響程度也不相同。坡度越大,稻草覆蓋對壤中流的影響越大,而對硝態(tài)氮的影響在15°時最大。

    (3) 對于紫色土坡耕地,在采用稻草覆蓋增加入滲的情況下,應(yīng)采取相應(yīng)的措施保水保肥。

    猜你喜歡
    結(jié)皮耕作層硝態(tài)
    耕地耕作層資源保護(hù)利用的意義、問題和對策
    NaCl與Na2SO4風(fēng)沙土鹽結(jié)皮力學(xué)特性比較分析?
    黃土高原不同發(fā)育階段生物結(jié)皮的導(dǎo)水和持水特征
    黃土高原典型生物結(jié)皮對坡面產(chǎn)流產(chǎn)沙過程的影響
    耕作層剝離再利用實(shí)踐與政策差距淺談
    耕作層剝離再利用實(shí)踐與政策差距分析
    土壤風(fēng)蝕環(huán)境下鹽結(jié)皮研究及展望?
    低C/N比污水反硝化過程中亞硝態(tài)氮累積特性研究
    硝態(tài)氮供應(yīng)下植物側(cè)根生長發(fā)育的響應(yīng)機(jī)制
    關(guān)于耕作層土壤剝離用于土壤培肥的必要條件探討
    极品教师在线视频| 午夜福利高清视频| 春色校园在线视频观看| 国产免费一级a男人的天堂| 国产视频内射| 男女下面进入的视频免费午夜| 日本在线视频免费播放| 亚洲av.av天堂| 午夜精品一区二区三区免费看| 国产精品爽爽va在线观看网站| 黄色一级大片看看| 午夜精品在线福利| 无遮挡黄片免费观看| 久久久久精品国产欧美久久久| 国国产精品蜜臀av免费| 成人国产综合亚洲| 1024手机看黄色片| 搡老岳熟女国产| 亚洲人成网站高清观看| 乱人视频在线观看| 网址你懂的国产日韩在线| 久久久色成人| 国产一区二区亚洲精品在线观看| 少妇高潮的动态图| 精品人妻熟女av久视频| 黄片wwwwww| 美女被艹到高潮喷水动态| 丝袜美腿在线中文| 国产爱豆传媒在线观看| 国产久久久一区二区三区| 亚洲国产欧美人成| 欧美区成人在线视频| 麻豆精品久久久久久蜜桃| 香蕉av资源在线| 国产午夜精品论理片| 久久久色成人| 久久亚洲精品不卡| 亚洲中文日韩欧美视频| 亚洲精华国产精华液的使用体验 | 简卡轻食公司| 国产不卡一卡二| 亚洲av中文字字幕乱码综合| 亚洲国产精品成人综合色| 久久热精品热| 级片在线观看| 少妇裸体淫交视频免费看高清| 欧美三级亚洲精品| 特大巨黑吊av在线直播| 男女边吃奶边做爰视频| 国产老妇女一区| 国产精品伦人一区二区| 亚洲国产欧美人成| 夜夜看夜夜爽夜夜摸| 精品日产1卡2卡| a级毛片免费高清观看在线播放| av福利片在线观看| 美女高潮的动态| 变态另类丝袜制服| 久久这里只有精品中国| 国产精品综合久久久久久久免费| 亚洲不卡免费看| 亚洲人成伊人成综合网2020| 亚洲成人久久爱视频| 三级毛片av免费| 91精品国产九色| 国产在视频线在精品| 真实男女啪啪啪动态图| 国产 一区 欧美 日韩| 免费人成在线观看视频色| 中文资源天堂在线| av在线亚洲专区| 精品福利观看| 内射极品少妇av片p| 免费观看在线日韩| 婷婷精品国产亚洲av在线| 国产蜜桃级精品一区二区三区| 在线观看66精品国产| 嫩草影视91久久| 午夜免费成人在线视频| 一进一出好大好爽视频| 精品人妻视频免费看| 亚洲内射少妇av| 观看美女的网站| 亚洲七黄色美女视频| 国产精品自产拍在线观看55亚洲| 女生性感内裤真人,穿戴方法视频| 色吧在线观看| 久久久久久伊人网av| 一级毛片久久久久久久久女| 午夜福利在线观看吧| 亚洲av免费高清在线观看| 亚洲欧美日韩高清专用| 一个人免费在线观看电影| 免费观看的影片在线观看| 非洲黑人性xxxx精品又粗又长| 精品午夜福利视频在线观看一区| 啦啦啦观看免费观看视频高清| 久久精品国产亚洲av涩爱 | 国模一区二区三区四区视频| 国产av一区在线观看免费| 变态另类成人亚洲欧美熟女| 久久久久久久午夜电影| 国产成人福利小说| 在线国产一区二区在线| 老熟妇乱子伦视频在线观看| 毛片女人毛片| 熟妇人妻久久中文字幕3abv| 熟女电影av网| 久久久精品大字幕| 麻豆精品久久久久久蜜桃| 嫩草影院新地址| 日本欧美国产在线视频| 尤物成人国产欧美一区二区三区| 亚洲专区国产一区二区| 国产精品国产三级国产av玫瑰| 亚洲黑人精品在线| 精品人妻一区二区三区麻豆 | 自拍偷自拍亚洲精品老妇| 中文字幕熟女人妻在线| 又爽又黄无遮挡网站| 成人永久免费在线观看视频| 国产伦精品一区二区三区视频9| 久久亚洲真实| 老师上课跳d突然被开到最大视频| 黄色日韩在线| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| a级一级毛片免费在线观看| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 亚洲成a人片在线一区二区| 在线观看舔阴道视频| 久久人人精品亚洲av| 俄罗斯特黄特色一大片| 99热这里只有精品一区| 日本黄大片高清| 国产国拍精品亚洲av在线观看| 国产淫片久久久久久久久| 国产探花在线观看一区二区| 亚洲精品一区av在线观看| 国产精品女同一区二区软件 | 亚洲美女视频黄频| 桃红色精品国产亚洲av| 午夜视频国产福利| 欧美日韩精品成人综合77777| 国产伦人伦偷精品视频| 97超级碰碰碰精品色视频在线观看| 欧美xxxx性猛交bbbb| 欧美成人性av电影在线观看| 伦精品一区二区三区| 在线播放无遮挡| 男人狂女人下面高潮的视频| 亚洲熟妇中文字幕五十中出| 男人舔奶头视频| 亚洲成av人片在线播放无| 可以在线观看毛片的网站| 韩国av在线不卡| 亚洲av免费高清在线观看| 毛片一级片免费看久久久久 | 欧美xxxx性猛交bbbb| 亚洲最大成人手机在线| 男女视频在线观看网站免费| 男人舔奶头视频| 国国产精品蜜臀av免费| 极品教师在线视频| 午夜福利在线观看免费完整高清在 | 别揉我奶头 嗯啊视频| 小说图片视频综合网站| 久久精品综合一区二区三区| 色av中文字幕| 国产极品精品免费视频能看的| 亚洲国产欧美人成| 日韩欧美在线二视频| 欧美日韩黄片免| 噜噜噜噜噜久久久久久91| 成人鲁丝片一二三区免费| 精品福利观看| 女人被狂操c到高潮| 亚洲欧美激情综合另类| 看免费成人av毛片| 免费高清视频大片| 99久国产av精品| 在线观看66精品国产| 亚洲精品色激情综合| 久99久视频精品免费| 熟妇人妻久久中文字幕3abv| 国产男人的电影天堂91| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 久久久久久大精品| 直男gayav资源| 国产免费一级a男人的天堂| 亚洲最大成人av| 精品午夜福利视频在线观看一区| 亚洲一区二区三区色噜噜| 九九热线精品视视频播放| 老司机午夜福利在线观看视频| 免费av不卡在线播放| 国产精品一区二区三区四区久久| 极品教师在线视频| 亚洲专区中文字幕在线| 嫩草影院入口| 亚洲精华国产精华液的使用体验 | 亚洲国产日韩欧美精品在线观看| 日本三级黄在线观看| 简卡轻食公司| 午夜激情福利司机影院| 99久久无色码亚洲精品果冻| 日本精品一区二区三区蜜桃| 国产主播在线观看一区二区| 午夜a级毛片| 亚洲内射少妇av| 美女高潮的动态| 在线天堂最新版资源| 久久精品综合一区二区三区| 日韩欧美国产在线观看| 亚洲av中文av极速乱 | 久久久精品大字幕| 国产老妇女一区| 人人妻人人澡欧美一区二区| 丰满乱子伦码专区| 日本 欧美在线| 人妻久久中文字幕网| 亚洲电影在线观看av| 精品久久久久久久久久久久久| 91在线观看av| 国产久久久一区二区三区| 美女 人体艺术 gogo| 少妇高潮的动态图| 久久99热这里只有精品18| 亚洲最大成人手机在线| 超碰av人人做人人爽久久| 能在线免费观看的黄片| a级一级毛片免费在线观看| 国产免费一级a男人的天堂| 少妇的逼好多水| 18禁黄网站禁片午夜丰满| 伦精品一区二区三区| 日韩欧美精品v在线| 麻豆国产97在线/欧美| 欧美成人一区二区免费高清观看| 成人二区视频| 亚洲精华国产精华液的使用体验 | 男插女下体视频免费在线播放| 一级av片app| 亚洲国产精品合色在线| 99国产极品粉嫩在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久久久久末码| 亚洲男人的天堂狠狠| 网址你懂的国产日韩在线| www日本黄色视频网| 国产成人福利小说| 能在线免费观看的黄片| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 国产亚洲欧美98| 亚洲第一区二区三区不卡| 欧美成人性av电影在线观看| 伦精品一区二区三区| 哪里可以看免费的av片| 老女人水多毛片| 国产一区二区三区av在线 | 成人无遮挡网站| 国产精品野战在线观看| 亚洲成人精品中文字幕电影| 深爱激情五月婷婷| 欧美3d第一页| 免费av不卡在线播放| 国产黄a三级三级三级人| 国产精品永久免费网站| 婷婷色综合大香蕉| 国产亚洲欧美98| 免费av不卡在线播放| 免费电影在线观看免费观看| 久久久久国内视频| 国产高清三级在线| 变态另类成人亚洲欧美熟女| 女生性感内裤真人,穿戴方法视频| 麻豆精品久久久久久蜜桃| 欧美日韩国产亚洲二区| 国内少妇人妻偷人精品xxx网站| 亚洲狠狠婷婷综合久久图片| 特大巨黑吊av在线直播| 麻豆久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 中文字幕精品亚洲无线码一区| 午夜精品一区二区三区免费看| 免费在线观看影片大全网站| 欧美性猛交╳xxx乱大交人| 真实男女啪啪啪动态图| 精品一区二区免费观看| 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| 在线免费十八禁| 三级男女做爰猛烈吃奶摸视频| 欧美一区二区亚洲| 久久久久免费精品人妻一区二区| 久久精品综合一区二区三区| av专区在线播放| 97超视频在线观看视频| 九色国产91popny在线| 午夜福利欧美成人| 欧美极品一区二区三区四区| av国产免费在线观看| 国产日本99.免费观看| 国内精品美女久久久久久| 欧美黑人巨大hd| 欧美激情国产日韩精品一区| 黄色视频,在线免费观看| 日日摸夜夜添夜夜添小说| 国产91精品成人一区二区三区| 中出人妻视频一区二区| 神马国产精品三级电影在线观看| 麻豆国产av国片精品| 在线a可以看的网站| 日韩av在线大香蕉| 亚洲国产精品合色在线| 69av精品久久久久久| 国产探花极品一区二区| 免费大片18禁| 搡老熟女国产l中国老女人| 精品久久久久久久久av| 九色国产91popny在线| 日韩欧美一区二区三区在线观看| 蜜桃久久精品国产亚洲av| 美女免费视频网站| 一级av片app| 人人妻人人澡欧美一区二区| 欧美3d第一页| 亚洲成人中文字幕在线播放| 久久草成人影院| 亚洲自拍偷在线| 美女 人体艺术 gogo| ponron亚洲| 禁无遮挡网站| 亚洲电影在线观看av| www.www免费av| av天堂中文字幕网| a在线观看视频网站| 精品福利观看| 久久久久国内视频| 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 国产视频一区二区在线看| 啪啪无遮挡十八禁网站| 精品福利观看| 久久久久国内视频| 日本-黄色视频高清免费观看| 欧美在线一区亚洲| 又黄又爽又免费观看的视频| 搡老妇女老女人老熟妇| 久久久久久久久久成人| 岛国在线免费视频观看| 国产毛片a区久久久久| 国产伦人伦偷精品视频| 成人国产综合亚洲| 欧美成人一区二区免费高清观看| 全区人妻精品视频| 1024手机看黄色片| 日韩一本色道免费dvd| 亚洲国产色片| 国产老妇女一区| 国内毛片毛片毛片毛片毛片| 欧美国产日韩亚洲一区| ponron亚洲| 亚洲欧美日韩卡通动漫| 淫妇啪啪啪对白视频| 免费看a级黄色片| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 欧美+亚洲+日韩+国产| 99国产精品一区二区蜜桃av| 最近最新中文字幕大全电影3| 国产精品国产高清国产av| 成人国产一区最新在线观看| 亚洲熟妇中文字幕五十中出| 日本撒尿小便嘘嘘汇集6| 搡老熟女国产l中国老女人| 少妇被粗大猛烈的视频| 久久久久久久精品吃奶| 日日摸夜夜添夜夜添小说| 国产高清激情床上av| 国产日本99.免费观看| 麻豆成人av在线观看| 男女边吃奶边做爰视频| 亚洲av电影不卡..在线观看| 在线观看免费视频日本深夜| 久久久精品欧美日韩精品| 国产又黄又爽又无遮挡在线| 2021天堂中文幕一二区在线观| 午夜福利在线观看吧| 成人一区二区视频在线观看| 12—13女人毛片做爰片一| 亚洲最大成人av| 午夜久久久久精精品| 中文字幕av成人在线电影| 狠狠狠狠99中文字幕| 国产成人a区在线观看| 久久热精品热| 春色校园在线视频观看| 91在线观看av| 91麻豆精品激情在线观看国产| 亚洲人成伊人成综合网2020| 亚洲最大成人中文| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 久久婷婷人人爽人人干人人爱| 日本三级黄在线观看| 久久中文看片网| 男女那种视频在线观看| 久久久久久久久久黄片| 国产在线男女| 日本a在线网址| 日韩精品中文字幕看吧| 午夜福利在线观看吧| 欧美三级亚洲精品| 人妻少妇偷人精品九色| 村上凉子中文字幕在线| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 又爽又黄a免费视频| 国产精品,欧美在线| 大又大粗又爽又黄少妇毛片口| 99在线人妻在线中文字幕| 国产色婷婷99| 国产精品精品国产色婷婷| 欧美+日韩+精品| 99久久无色码亚洲精品果冻| 亚洲第一区二区三区不卡| 人妻少妇偷人精品九色| 无遮挡黄片免费观看| 精品久久久久久成人av| 色吧在线观看| 成人精品一区二区免费| 欧美成人性av电影在线观看| 国产亚洲精品久久久久久毛片| 久9热在线精品视频| 亚洲专区中文字幕在线| 蜜桃亚洲精品一区二区三区| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 亚洲国产欧美人成| 亚洲最大成人中文| 午夜福利视频1000在线观看| 大型黄色视频在线免费观看| 如何舔出高潮| 亚洲av.av天堂| 亚洲在线观看片| 日韩强制内射视频| 国内精品一区二区在线观看| 国产精华一区二区三区| 黄片wwwwww| 男人舔奶头视频| 国产美女午夜福利| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩卡通动漫| 中文字幕高清在线视频| 热99在线观看视频| 成人亚洲精品av一区二区| 成人国产一区最新在线观看| 国产人妻一区二区三区在| 亚洲av二区三区四区| av视频在线观看入口| av在线亚洲专区| 亚洲精品国产成人久久av| 亚洲中文日韩欧美视频| 免费高清视频大片| 国产色婷婷99| 免费电影在线观看免费观看| 高清在线国产一区| 久久国产乱子免费精品| 久久99热6这里只有精品| 久久热精品热| 午夜老司机福利剧场| 亚洲国产日韩欧美精品在线观看| 国产欧美日韩精品亚洲av| 亚洲第一电影网av| www日本黄色视频网| 中出人妻视频一区二区| 亚洲国产精品sss在线观看| 无人区码免费观看不卡| 久久精品91蜜桃| 99久久精品国产国产毛片| 亚洲国产高清在线一区二区三| 性色avwww在线观看| 免费无遮挡裸体视频| 日本免费a在线| 亚洲va日本ⅴa欧美va伊人久久| 校园人妻丝袜中文字幕| 国产成人一区二区在线| a级毛片免费高清观看在线播放| 老师上课跳d突然被开到最大视频| 午夜福利欧美成人| 人妻久久中文字幕网| 最好的美女福利视频网| 国产人妻一区二区三区在| 国产一级毛片七仙女欲春2| 中国美女看黄片| 久久久久国产精品人妻aⅴ院| 免费大片18禁| 欧美日本视频| 久久精品国产亚洲av香蕉五月| 校园春色视频在线观看| 国产一区二区三区av在线 | 亚洲精品一卡2卡三卡4卡5卡| 俄罗斯特黄特色一大片| 日日夜夜操网爽| 欧美黑人欧美精品刺激| 尤物成人国产欧美一区二区三区| 91av网一区二区| 午夜激情欧美在线| 精品人妻偷拍中文字幕| 午夜精品在线福利| 亚洲中文字幕日韩| 国产亚洲精品久久久com| 麻豆成人av在线观看| 欧美日韩乱码在线| 韩国av一区二区三区四区| АⅤ资源中文在线天堂| 天堂网av新在线| 久久久久久久久久黄片| 亚洲美女搞黄在线观看 | 人人妻,人人澡人人爽秒播| 男插女下体视频免费在线播放| 国语自产精品视频在线第100页| 免费看美女性在线毛片视频| 成人一区二区视频在线观看| 亚洲自偷自拍三级| 日韩欧美国产在线观看| 丰满的人妻完整版| 亚洲人成网站在线播| 五月伊人婷婷丁香| 亚洲av不卡在线观看| av在线蜜桃| 午夜福利视频1000在线观看| 黄色日韩在线| 两个人视频免费观看高清| 在线播放无遮挡| 久久久久久大精品| 窝窝影院91人妻| 亚洲18禁久久av| 久久精品国产亚洲av天美| 午夜a级毛片| 成人av一区二区三区在线看| 毛片女人毛片| 人妻夜夜爽99麻豆av| 给我免费播放毛片高清在线观看| 亚洲电影在线观看av| 性色avwww在线观看| 国产高清不卡午夜福利| 夜夜爽天天搞| av福利片在线观看| 色精品久久人妻99蜜桃| 白带黄色成豆腐渣| 色精品久久人妻99蜜桃| 在线播放无遮挡| 变态另类丝袜制服| 久久久久久久亚洲中文字幕| 午夜福利欧美成人| 亚洲成a人片在线一区二区| 日韩欧美国产一区二区入口| 桃红色精品国产亚洲av| 欧美不卡视频在线免费观看| 69av精品久久久久久| 韩国av一区二区三区四区| 亚洲精品粉嫩美女一区| 免费观看精品视频网站| 香蕉av资源在线| 中文字幕免费在线视频6| 3wmmmm亚洲av在线观看| 午夜激情福利司机影院| 三级国产精品欧美在线观看| 男人的好看免费观看在线视频| 国产视频一区二区在线看| 高清毛片免费观看视频网站| 国产视频一区二区在线看| 亚洲av免费高清在线观看| 精品一区二区三区人妻视频| 国产人妻一区二区三区在| 色哟哟哟哟哟哟| 午夜爱爱视频在线播放| 91午夜精品亚洲一区二区三区 | 国产精品1区2区在线观看.| www.www免费av| 不卡一级毛片| 国产亚洲精品久久久com| 欧美日韩黄片免| 午夜爱爱视频在线播放| 日韩精品青青久久久久久| xxxwww97欧美| 日韩中字成人| 欧美国产日韩亚洲一区| 高清日韩中文字幕在线| av在线老鸭窝| 国产av不卡久久| 精品久久久久久久人妻蜜臀av| 精品久久久久久久久av| 看黄色毛片网站| 九九热线精品视视频播放| 亚洲无线观看免费| 岛国在线免费视频观看| 在线观看免费视频日本深夜| 日本一二三区视频观看| 又粗又爽又猛毛片免费看| 久久99热6这里只有精品| 国产不卡一卡二| 亚洲人成网站高清观看| 露出奶头的视频| 少妇熟女aⅴ在线视频| videossex国产| 国产一区二区三区视频了| 国产欧美日韩一区二区精品| 国产男靠女视频免费网站| 亚洲av电影不卡..在线观看|