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

    生物炭對不同坡度坡耕地土壤水動力學參數(shù)的影響

    2019-04-01 12:55:52魏永霞
    農(nóng)業(yè)機械學報 2019年3期
    關鍵詞:導水率水率非飽和

    魏永霞 王 鶴 吳 昱 劉 慧

    (1.東北農(nóng)業(yè)大學水利與土木工程學院, 哈爾濱 150030;2.東北農(nóng)業(yè)大學農(nóng)業(yè)農(nóng)村部農(nóng)業(yè)水資源高效利用重點實驗室, 哈爾濱 150030;3.黑龍江農(nóng)墾勘測設計研究院, 哈爾濱 150090; 4.東北林業(yè)大學林學院, 哈爾濱 150040;5.東北農(nóng)業(yè)大學理學院, 哈爾濱 150030)

    0 引言

    生物炭是指在限氧或無氧的環(huán)境條件下,通過高溫(一般小于700℃)裂解將小薪柴、農(nóng)作物秸稈、雜草等生物質經(jīng)碳化而形成的一種穩(wěn)定難溶、高度芳香化、碳含量極其豐富(高達60%)的固態(tài)產(chǎn)物[1]。早有研究表明,將生物炭施入土壤,可降低土壤容重、提高總孔隙率,這種結構改變對土壤持水能力也產(chǎn)生一定影響,如提高土壤含水率和降水的入滲量等[2],特別是土壤中有效水含量的增加更利于作物的生長[3-4];生物炭的吸濕能力比土壤高1~2個數(shù)量級[5],使含有生物炭的土壤比普通土壤田間持水率高約18%[6]。通常認為,生物炭對土壤物理結構等性狀的改良以及對土壤水分的影響均與生物炭自身具有多孔結構和吸附能力密不可分[7-8]。隨著生物炭對土壤性質改良、土壤水分調控等研究的進一步深入,部分研究逐漸著眼于生物炭對土壤水動力學參數(shù)的影響,如田丹等[9]在研究中得出,添加0.05、0.10、0.15 g/g的秸稈木炭和花生殼炭均可有效減小土壤水分擴散率,增加粉砂壤土的土壤持水性。劉志凱[10]在研究中也發(fā)現(xiàn),當土壤體積含水率小于0.4 cm3/cm3時,生物炭抑制土壤水分擴散率,當土壤體積含水率大于0.4 cm3/cm3時,生物炭則提高土壤水分擴散率,且生物炭可增加土壤的非飽和導水率,提高其導水能力。

    位于我國松花江流域和遼河流域的東北黑土區(qū)是全球三大黑土區(qū)之一,東北黑土區(qū)土壤主要以黑土、黑鈣土和草甸黑土為主,黑土土壤肥沃、高產(chǎn)且結構性好,具有東北地區(qū)主要農(nóng)作物生長所需的適宜土壤條件[11]。但隨著黑土區(qū)水土流失現(xiàn)象日益嚴重,黑土耕層逐步變薄,使土壤中有機質含量下降,土壤一系列理化性狀也逐漸惡化,使土壤容重增大、孔隙度減小,進而大大降低了土壤持水性能[12]。東北黑土區(qū)耕地中有60%為2°~5°坡耕地[13],多年以來人們過度開墾和掠奪式經(jīng)營導致水土流失加劇,使黑土區(qū)農(nóng)作物實際平均產(chǎn)量僅為潛在產(chǎn)量的48.49%[14],東北黑土區(qū)作為我國糧食主要生產(chǎn)基地,改善其土壤性質、促進其土壤水分向作物所需方向運動,以保證糧食安全已刻不容緩。生物炭作為土壤改良劑、土壤保水劑、導水劑等,將會改善此現(xiàn)象[15-16]。

    目前針對生物炭對土壤水動力學參數(shù)影響的研究多集中在無坡度土壤的一年試驗中,對于不同坡度土壤水動力學參數(shù)在施加生物炭后兩年及以上的研究較少。本文對不同坡度坡耕地黑土中添加適量生物炭,研究分析坡度變化和添加生物炭后兩年內對土壤水動力學參數(shù)的影響,并使用HYDRUS-1D軟件模擬土壤水分運動情況,以期為高效利用東北黑土區(qū)坡耕地農(nóng)業(yè)水土資源提供相應理論依據(jù)與方法。

    1 材料與方法

    1.1 試驗區(qū)概況

    試驗區(qū)位于黑龍江省水利科學研究院綜合試驗基地(45°43′09″N,126°36′35″E),總面積55 hm2,多年平均氣溫3.1℃,無霜期130~140 d,多年平均水面蒸發(fā)量796 mm,耕地土壤多以黑土為主,入滲能力弱。年平均降水量介于400~650 mm,多集中在7—9月,且歷時較短,僅7—9月的降水量就占全年降水總量的70%以上。供試土壤主要為黑土,速效氮(N)質量比為154.4 mg/kg,速效磷(P2O5)質量比為40.1 mg/kg,速效鉀(K2O)質量比為376.8 mg/kg,pH值7.27。0~80 cm土層平均田間持水率(質量分數(shù))為29.4%,土壤干容重為1.22 g/cm3。大豆是該地區(qū)主要農(nóng)作物之一。

    1.2 試驗材料

    供試生物炭購于遼寧金和福開發(fā)有限公司,使用玉米秸稈在450℃無氧條件下燒制而成,其基本理化性質如表1所示。供試大豆品種為黑河3號。

    表1 生物炭的基本性質Tab.1 Basic feature of biochar

    1.3 試驗設計

    試驗在黑龍江省水利科學研究院綜合試驗基地的徑流小區(qū)內進行,小區(qū)規(guī)格為2 m×5 m,坡度選擇1.5°、3°、5°共3種在東北黑土區(qū)比較有代表性的坡度。每種坡度的徑流小區(qū)分別設置施加生物炭和不施加生物炭處理,小區(qū)編號分別為T1.5、T3、T5和CK1.5、CK3、CK5,共計6個小區(qū)。根據(jù)前期研究成果[16],生物炭施用量選擇施用當年效果較好的75 kg/hm2,施用時間為2016年,2017年不再施加。生物炭施用方法為播種前將生物炭粉均勻鋪撒于土壤表面,人工攪拌至與0~20 cm表層土壤充分均勻混合后靜止待用。

    1.4 測定參數(shù)及方法

    各參數(shù)測定按坡上、坡中、坡下3點取樣,取其平均值,每個樣點3次重復。

    1.4.1土壤水分常數(shù)

    采用環(huán)刀法和干燥法測定田間持水率和飽和含水率。使用環(huán)刀在徑流小區(qū)各位置點取深度15~20 cm處土壤作為試驗土樣,將土樣帶回室內處理后,采用干燥法測量土壤含水率。

    1.4.2土壤水分特征曲線

    土壤水分特征曲線(Soil water characteristic curve, SWCC)采用離心機法進行測定。使用CR21-GIII型高速離心機配套環(huán)刀取樣帶回,環(huán)刀土樣處理后放入盛有蒸餾水的容器中浸泡12 h,水面低于環(huán)刀上緣2 mm,土壤達飽和后取出環(huán)刀離心備用。離心機轉速設定為310、540、690、980、1 700、2 190、3 100、5 370、6 210、6 930、8 210、9 300 r/min,相對應平衡時間為10、12、16、26、42、49、58、73、77、80、85、89 min,每次離心結束后將離心的水分擦去,然后將環(huán)刀與土樣一同稱量并記錄數(shù)據(jù)。

    目前應用最為廣泛的推求土壤水分特征曲線的模型是Van Genuchten模型,表達式為

    (1)

    為了更加明顯地顯示體積含水率與負壓水頭間的關系,將式(1)轉換為

    (2)

    (3)

    (4)

    式中θ——體積含水率,cm3/cm3

    h——負壓水頭,cm

    θr——殘余含水率,cm3/cm3

    θs——飽和含水率,cm3/cm3

    α、n、m——經(jīng)驗擬合參數(shù)(或曲線性狀參數(shù)),其中α為進氣值倒數(shù)

    hb——進氣吸力,cm

    1.4.3土壤飽和導水率

    土壤飽和導水率(Ks)采用定水頭法測定。所用設備為自制馬氏瓶和內徑10 cm、高80 cm的有機玻璃筒(圖1)。進行試驗裝土時,先在有機玻璃筒壁涂滿凡士林后,將土柱內壁邊緣土壤壓實,確保無邊緣效應。試驗過程中,馬氏瓶作為供水系統(tǒng)裝置將供水水頭控制在200 cm。在定水頭作用下放置一段時間,待出水管管口有水溢出且保持均勻出水量時,則認為土樣已達飽和,此時開始試驗,各處理每隔60 min測量一次出流量,并記錄數(shù)據(jù)。

    圖1 土壤飽和導水率試驗裝置示意圖Fig.1 Schematic of experimental device for soil saturated hydraulic conductivity1.平臺 2.進水口 3.馬氏瓶 4.進氣管 5.土柱 6.滴水管7.止水閥 8.排水孔 9.濾層

    土壤飽和導水率計算公式為

    (5)

    式中Ks——土壤飽和導水率,cm/min

    V——一定時間下出流水量,cm3

    L——土樣長度,cm

    H——水頭差,cm

    A——土樣截面積,cm2

    t——滲流時間,min

    1.4.4土壤非飽和導水率

    土壤非飽和導水率直接測量十分困難且復雜,故采用間接公式推求方法來獲取。根據(jù)Van Genuchten模型及土壤飽和導水率可計算出土壤非飽和導水率理論值

    (6)

    式中K(h)——土壤非飽和導水率,cm/min

    1.4.5土壤非飽和水分擴散率

    土壤非飽和水分擴散率采用水平土柱法測定。所用裝備示意圖如圖2所示,自制帶通氣管的馬氏瓶與直徑5 cm、長60 cm的圓柱透明有機玻璃管,在裝土過程中,要保證各處理填裝過程都與對照組具有統(tǒng)一的擊實次數(shù)和擊實壓力。各處理土壤均分層裝入試驗裝置的透明有機玻璃管中,填裝過程也要保證各處理土壤的均一性。供水系統(tǒng)依然由馬氏瓶供水,待濕潤峰到達整個土柱的2/3(40 cm)時,關閉供水閥并記錄試驗結束時間,同時按節(jié)取出土壤,用干燥法測定土壤含水率,記錄試驗數(shù)據(jù)。

    圖2 土壤水分擴散率試驗裝置示意圖Fig.2 Schematic of experimental device for soil water diffusion rate1.馬氏瓶 2.儲水室 3.濾料室 4.過濾板 5.土柱

    水平入滲是測定土壤水分擴散率非穩(wěn)定流的一種方法,此法利用水平土柱吸滲試驗數(shù)據(jù),再結合解析法計算出土壤水分擴散率。采用Boltzmann變換,對一維水平流動微分方程求解,得出

    (7)

    將式(7)轉換為差分的形式

    (8)

    式中D(θ)——土壤非飽和水分擴散率,cm2/min

    θa——初始土壤含水率,cm3/cm3

    ξ——Boltzmann變換參數(shù)

    1.4.6比水容量

    對Van Genuchten模型式(2)兩邊同時求導即可得到土壤的比水容量

    (9)

    式中C(h)——土壤比水容量,cm3/cm4

    1.5 數(shù)據(jù)處理方法

    采用Excel 2010、Origin 9.1進行基本數(shù)據(jù)處理和繪圖,利用SPSS 19.0進行擬合度檢驗和顯著性分析,顯著性水平為0.05、0.01,用HYDRUS-1D軟件對土壤水分運動情況進行擬合。

    2 結果與分析

    2.1 生物炭對坡耕地土壤水分常數(shù)的影響

    2016、2017年各處理飽和含水率和田間持水率如表2所示。無施加生物炭條件下隨地形坡度增加飽和含水率與田間持水率減小,與已有結果一致[17]。施加生物炭當年飽和含水率與田間持水率均有所增加,且隨地形坡度增加二者變幅增強,坡度為1.5°、3°、5°時田間持水率增加率分別為4.89%、7.35%、10.98%,飽和含水率增加率分別是10.91%、15.10%、16.46%。施加生物炭次年飽和含水率與田間持水率也均有增加,且隨地形坡度增加二者變幅增強,坡度為1.5°、3°、5°時田間持水率增加率分別為6.69%、10.03%、11.94%,飽和含水率增加率分別為11.91%、15.48%、16.18%。

    表2 各處理的土壤水分常數(shù)Tab.2 Soil moisture constants of each treatment cm3/cm3

    施加生物炭增大了兩土壤水分常數(shù),這與生物炭自身具有較低容重、較高的吸水能力密不可分,生物炭降低土壤容重,增加土壤孔隙度,使土壤含水率增加、土壤儲水空間變大,進而提高土壤持水能力。隨地形坡度增加兩土壤水分常數(shù)增幅增大,這是因為地形坡度越大,地表徑流和地下水越容易流失,土壤水土流失情況越加嚴重,生物炭的強吸附性可使更多流失的水分留于土壤中,對土壤持水性改善能力更大。但2017年數(shù)值均較2016年有小幅減少,這與施加生物炭后的年限有關。

    為進一步分析坡度和是否施加生物炭兩因素對各土壤水分常數(shù)影響的顯著程度,分別建立2016年飽和含水率(y11)、田間持水率(y21)和2017年飽和含水率(y12)、田間持水率(y22)關于坡度和是否施用生物炭的回歸方程,其中坡度因子(x1)為定量變量,施用生物炭(x2)為定性變量,x2取0表示不施用生物炭,取1表示施用生物炭,得出回歸方程分別為

    y11=0.430-0.002x1+0.6x2
    (R2=0.976,F(xiàn)=62.25,P=0.004)

    (10)

    y21=0.329-0.001x1+0.25x2
    (R2=0.905,F(xiàn)=14.23,P=0.029)

    (11)

    y12=0.425-0.003x1+0.061x2
    (R2=0.985,F(xiàn)=100.26,P=0.002)

    (12)

    y22=0.319-0.001x1+0.03x2
    (R2=0.955,F(xiàn)=31.78,P=0.010)

    (13)

    4個方程的R2均大于0.9,P值均小于0.05,表明回歸方程擬合效果較好。從兩個變量的顯著性看,生物炭施用當年,坡度對飽和含水率和田間持水率的影響均不顯著(P分別為0.653、0.823),而是否施用生物炭對飽和含水率影響極顯著(P=0.002),對田間持水率的影響也達顯著程度(P=0.013);生物炭施用次年,坡度對飽和含水率和田間持水率的影響均不顯著(P分別為0.381、0.841),而是否施用生物炭對飽和含水率和田間持水率影響均為極顯著(P分別為0.001、0.004)。說明是否施用生物炭對飽和含水率和田間持水率的影響顯著,坡度因子則不顯著。

    2.2 生物炭對坡耕地土壤持水性的影響

    2.2.1土壤水分特征曲線的擬合參數(shù)

    土壤水分特征曲線可較好地反映土壤持水性能,采用HYDRUS-1D軟件中Van Genuchten模型對各處理土壤水分特征曲線進行擬合,得到2016、2017年擬合參數(shù)及決定系數(shù),如表3所示,各處理決定系數(shù)R2均大于0.99,表明HYDRUS-1D軟件中Van Genuchten模型對各處理均適用。同時從表3可看出,生物炭對參數(shù)θs和α有增大的作用,對參數(shù)θr和n有減小的作用,且參數(shù)θs、θr和n的變幅與坡度呈正相關,參數(shù)α的變幅與坡度呈負相關。

    表3 各處理Van Genuchten模型的擬合參數(shù)Tab.3 Fitting parameters of Van Genuchten model for each treatment

    2.2.2不同坡度施加生物炭對土壤水分特征曲線的影響

    圖3 施用生物炭不同坡度的土壤水分特征曲線Fig.3 Soil water characteristic curves of different slopes with application of biological carbon

    土壤水分特征曲線如圖3所示,在低吸力段(h<500 cm)時,曲線坡度陡峭而密集;當吸力高于500 cm時,曲線變得平緩稀疏。這是由于在低吸力段,土壤排水主要是在大孔隙中進行,盡管吸力變化不大,但足夠引起大孔隙土壤內含水率發(fā)生明顯變化,故SWCC呈“陡直”狀;當吸力在500~9 000 cm時,土壤中只有較小的孔隙能保留水分,土壤對水分的吸持力較大,土壤含水率隨吸力的變化不明顯,故在此吸力范圍內SWCC呈“平緩”狀。兩年中T1.5、T3、T5的SWCC均位于CK1.5、CK3、CK5曲線上方,且坡度越大二者距離越遠,這是因為生物炭增加了土壤孔隙度,已有研究[18]表明,施加75 kg/hm2生物炭主要增加土壤的中、小孔隙相對比例,土壤中、小孔隙的水分不易流出,使有效水含量增加,進而增強土壤持水性,故施加生物炭的SWCC位于對照曲線上方,且地形坡度最大的土壤其θr減小幅度最大,有效水含量增加最多,土壤持水性最優(yōu),由于θr變幅與地形坡度呈正相關關系,土壤持水性與地形坡度亦呈正相關關系,故是否施加生物炭二者SWCC間距與地形坡度呈正相關關系。

    表4為不同坡度坡耕地施加生物炭處理較對照處理土壤體積含水率隨吸力的變化,2016年施用生物炭土壤含水率的總變化量在1.5°、3°、5°坡度分別增加0.282 2、0.292 5、0.404 0 cm3/cm3,總變化率為8.16%、8.69%、14.21%;2017年1.5°、3°、5°坡度的土壤含水率總變化量各增大0.268 7、0.281 9、0.397 5 cm3/cm3,總變化率為7.77%、8.37%、13.98%。2017年較2016年相比施用生物炭后土壤體積含水率總增加量略有減少,土壤持水性能改善略弱,并且隨地形坡度增大總變化率幅度也略有減少,說明施用生物炭次年作用效果較施用當年有所減弱,但由于生物炭降解較慢,對土壤持水性能仍有增加的作用。

    2.2.3生物炭對坡耕地土壤比水容量的影響

    土壤比水容量是指單位基質勢的變化引起的含水率變化,是衡量土壤水分對植物的有效性、反映土壤供水性能的一個重要指標,一般土壤水在吸力h=1 000 cm時的比水容量可以表征土壤的供水能力。各處理h=1 000 cm時土壤比水容量見表5,2016、2017年T1.5、T3、T5的比水容量均較CK1.5、CK3、CK5有所增加,2016年1.5°、3°、5°坡度各自增加了1.160 44×10-3、1.282 537×10-3、1.830 207×10-3cm3/cm4,2017年分別增加1.215 2×10-5、2.000 1×10-5、4.9×10-5cm3/cm4,兩年增長率均隨著坡度的增加而增大,表明生物炭可使土壤比水容量增加,且比水容量變幅隨地形坡度增加而增大。5°小區(qū)比水容量增長率最大與生物炭對其田間持水率和殘余含水率變幅最大有關,坡度最大其田間持水率增幅最大,殘余含水率減幅最大,使其土壤有效含水量增加最大,供給作物的水分增加量也就最多,進而供水能力增強最明顯,比水容量增長率最大。2017年比水容量較2016年小,說明施用生物炭次年與當年相比供水能力的增加程度會減弱。

    表4 不同坡度坡耕地施加生物炭處理較對照處理土壤體積含水率隨吸力的變化Tab.4 Changes of soil volumetric water content with suction in different slopes sloping farmland by biological carbon treatment

    表5 不同坡度土壤吸力1 000 cm時的比水容量Tab.5 Water capacity of soil with different slopes under suction of 1 000 cm cm3/cm4

    2.3 生物炭對坡耕地土壤非飽和導水率的影響

    土壤非飽和導水率代表土壤水分在水頭壓力差作用下流動的性能,受土壤容重、含水率、孔隙分布特征等因素的共同影響[19],當容重小、孔隙度大、有機質養(yǎng)分含量高時,K(h)往往相對較大。由圖4可知,當土壤含水率小于0.30 cm3/cm3時,隨著土壤含水率的增加非飽和導水率增加緩慢,增幅并不明顯;當土壤含水率大于0.30 cm3/cm3時,非飽和導水率隨θ的增加迅速增大。施用生物炭當年與次年T1.5、T3、T5的K(h)曲線均位于CK1.5、CK3、CK5上方,且隨地形坡度增大K(h)曲線略向左上方移動。施加生物炭使非飽和導水率增大主要是由于生物炭使土壤容重降低,孔隙度增大,增加了土壤對養(yǎng)分的吸附和有機質含量,進而使非飽和導水率增加,說明生物炭可以提高土壤在非飽和狀態(tài)下的導水能力,但坡度對土壤非飽和導水率的影響并不十分明顯。

    施加生物炭當年和次年T1.5、T3、T5的K(h)值較CK1.5、CK3、CK5均明顯增大,2016年1.5°、3°、5°坡度分別增加172.72%、239.61%、214.27%,2017年1.5°、3°、5°坡度分別增加143.03%、161.9%、164.04%。坡度相同時,2017年較2016年K(h)值均有所減小,且2017年較2016年施加生物炭后K(h)增量也減小。與SWCC相同,施加生物炭次年對土壤持水性、導水性仍有效但效果有所減弱,故施用生物炭K(h)依然增加但增加幅度減小,但坡度對土壤非飽和導水率的影響并無明顯規(guī)律。

    2.4 生物炭對坡耕地土壤非飽和水分擴散率的影響

    土壤非飽和水分擴散率D(θ)不僅反映出土壤孔隙狀況、孔隙分布以及導水性能,同時對土壤水分運動狀況也有一定影響[20]。土壤水分擴散率受土壤質地、土壤密度、孔隙度、有機質含量和地面坡度等諸多因素的影響[21],在相同介質、不同環(huán)境下表現(xiàn)為不相同性。由水平土柱試驗所得數(shù)據(jù)計算得出Boltzmann的變換參數(shù)ξ,采用式(8)計算各處理2016、2017年的土壤非飽和水分擴散率,如圖5所示,各處理D(θ)隨著θ的增大迅速增大。采用經(jīng)驗公式D(θ)=aebθ擬合2016、2017年土壤非飽和水分擴散率與土壤含水率的關系,擬合結果如表6所示。各處理R2均大于0.96,擬合效果極為理想,表明土壤非飽和水分擴散率與土壤含水率符合經(jīng)驗公式呈顯著的指數(shù)型曲線變化。

    圖5 施用生物炭不同坡度的土壤非飽和水分擴散率Fig.5 Water diffusion rate of unsaturated soil on different slopes with application of biochar

    當θ小于等于0.42 cm3/cm3時,T1.5、T3、T5的D(θ)曲線均位于CK1.5、CK3、CK5曲線下方,說明生物炭此時抑制了土壤水分擴散;當θ大于0.42 cm3/cm3時,T1.5、T3、T5的D(θ)數(shù)值迅速增大,曲線向上越過CK1.5、CK3、CK5曲線,說明該階段生物炭的存在使土壤水分水平擴散能力有所提高,在降雨時促進水分移動,提高降雨時土壤入滲速率,對夏季多發(fā)暴雨的黑土區(qū)尤為重要。同時隨地形坡度增大生物炭的作用效果增大,由于坡度增加水土流失嚴重,導致土壤有機質大量流失,生物炭可起到滯留劑作用,且坡度大滯留有機質多,故D(θ)增加也就越多。根據(jù)徑流小區(qū)實際情況θ多處于0.20~0.35 cm3/cm3范圍內,即生物炭抑制土壤水分的水平運動,使更多水分留于土壤中進而促進了作物生長發(fā)育。2016年T1.5、T3、T5的D(θ)增加幅度相對2017年更大一些,即施用生物炭次年對土壤水分擴散率的影響比施用生物炭當年弱化一些。

    表6 土壤非飽和水分擴散率與土壤含水率的擬合方程參數(shù)Tab.6 Fitting equation parameters of water diffusion rate and soil water content

    3 討論

    土壤水分特征曲線、土壤非飽和導水率、土壤非飽和水分擴散率以及比水容量等是土壤水分運動基本方程中的重要參數(shù),反映了土壤孔隙狀況、導水能力、持水性能及供水能力,受到土壤質地、結構、坡度等多重因素的共同影響。本研究在東北黑土區(qū)坡耕地開展了施加適量生物炭的對比試驗,研究了坡度和施用生物炭當年、次年對上述土壤水動力學參數(shù)的影響。

    施加生物炭后的兩年內飽和含水率和田間持水率均增大,這與生物炭自身特性密不可分,施加生物炭改良土壤性質,使土壤含水率增加,持水能力加強,與已有研究[22-24]結果一致;并且坡度越大,效果越顯著,原因是坡度越大水土流失越嚴重,施用生物炭后提高了土壤持水能力,進而吸附更多流失的水分。由2016、2017年飽和含水率和田間持水率關于坡度和是否施用生物炭的回歸方程可知:是否施用生物炭對飽和含水率和田間持水率的影響均較坡度因子顯著。施用生物炭次年的2017年飽和含水率和田間持水率都比2016當年增加量少,這是由于施用生物炭后年限越長生物炭產(chǎn)生效果越弱。

    施用生物炭后兩年內對SWCC參數(shù)均有影響,生物炭對參數(shù)θs和α有增大的作用,對參數(shù)θr和n有減小的作用,與王丹丹等[26]研究結果吻合,且參數(shù)θs、θr和n的變幅與地形坡度呈正相關關系,參數(shù)α的變幅與地形坡度呈負相關關系。兩年內T1.5、T3、T5的SWCC均位于對照組上方,土壤持水能力增強[27],隨坡度增大兩曲線距離越大,且均是T5曲線與CK5曲線距離最大,土壤有效水含量的變化是導致這一現(xiàn)象的主要原因;但2017年各曲線間間距與2016年相比略小,這與生物炭降解緩慢有關。

    作物以相同的能量吸水,在不同基質勢下從各種土壤中所吸收的水量因比水容量不同而形成很大的差別,比水容量越高,土壤供水能力越強,持水量越大,土壤有效水含量越高,作物吸水量也就越大,在本研究中施用生物炭當年和次年都使比水容量增大,對土壤水分調控、作物生長均有可觀效果,同田丹等[9]對生物炭研究結果一致。

    本研究中,施加生物炭當年和次年T1.5、T3、T5小區(qū)的K(h)較CK1.5、CK3、CK5均增大,即土壤在非飽和狀態(tài)下的導水能力提高,劉志凱[10]、王睿垠等[28]得到過相似結論,但K(h)與坡度變化關系不明顯。生物炭增大K(h)的原因:土壤水分總向土水勢降低的方向移動,水分在吸力的作用下從較小的孔隙流向大孔隙進而流走,施加生物炭的土壤孔隙度較大,水分移動更為容易,與已有結論一致[29]。且施加生物炭次年較當年K(h)增加量減小,仍與施用年限和生物炭微降解相關。

    由施加生物炭當年與次年對不同坡度土壤水分擴散率的影響可知,當θ≤0.42 cm3/cm3時,施加生物炭抑制土壤水分擴散;當θ>0.42 cm3/cm3時,添加生物炭促進土壤水分擴散,東北黑土區(qū)土壤含水率多集中在0.20~0.35 cm3/cm3范圍內,即生物炭抑制本地區(qū)土壤水分擴散,已有研究得到類似結論[9],且與潘英華等[30]研究結果一致。坡度對土壤非飽和水分擴散率也有影響,在東北黑土區(qū)生物炭對土壤水分擴散的抑制程度隨地形坡度增大而變大,坡度變大土壤中可移動水量增多,故施加生物炭后對水分流失控制更大,更加有利于作物生長發(fā)育,且D(θ)-θ曲線符合指數(shù)模型,R2均大于0.96,擬合效果較為理想。施用生物炭次年土壤水分擴散率抑制效果比施用生物炭當年略差,可能是生物炭施用次年效果減弱的原因。

    本文僅分析了施加生物炭后兩年內的情況,對于施加生物炭多年后的效應有待增加試驗年限進一步探究。

    4 結論

    (1)施用生物炭當年與次年均可增大土壤的田間持水率和飽和含水率,并且地形坡度越大效果越顯著,其中施用生物炭次年較當年田間持水率和飽和含水率均略有減小。施用生物炭對飽和含水率和田間持水率的影響均較坡度因子更為顯著。

    (2)施用生物炭當年與次年對SWCC參數(shù)和土壤水吸力在1 000 cm時的比水容量均有影響,生物炭可增大參數(shù)θs、α和比水容量,減小參數(shù)θr和n,且參數(shù)θs、θr、n及比水容量的變幅與坡度呈正相關關系,參數(shù)α的變幅與坡度呈負相關關系。在各個土壤吸力條件下,施用生物炭當年和次年均使土壤含水率增加、土壤持水性增強,生物炭施用次年效果較施用當年減弱。

    (3)坡度對K(h)影響無明顯規(guī)律性,但施加生物炭可使K(h)增大,使非飽和土壤的導水性增強,并且施加生物炭當年的效果較次年更為顯著,2016年K(h)增長率以3°坡度施加適量生物炭的處理最高,為239.61%,2017年K(h)增長率以5°坡度施加適量生物炭的處理最高,為164.04%。

    (4)對于土壤含水率多集中在0.20~0.35 cm3/cm3范圍內的黑土區(qū)土壤,施加生物炭可降低D(θ),抑制土壤水分的水平運動,隨地形坡度增加抑制效果增強,并且施用生物炭次年較施用當年對土壤水分擴散率的影響小。

    猜你喜歡
    導水率水率非飽和
    煤骨料充填材料析水率試驗研究
    陜北山地蘋果園土壤飽和導水率和植物導水率特征
    基于線性源法與圖像處理的土壤飽和導水率快速測量方法
    非飽和原狀黃土結構強度的試驗研究
    工程與建設(2019年1期)2019-09-03 01:12:24
    海城市土壤飽和導水率空間分布特征
    非飽和多孔介質應力滲流耦合分析研究
    非飽和土基坑剛性擋墻抗傾覆設計與參數(shù)分析
    非飽和地基土蠕變特性試驗研究
    持水率測量結果影響因素試驗與分析
    草海流域3種優(yōu)勢樹種凋落物葉分解歷程中的水文特征
    欧美日韩精品网址| netflix在线观看网站| 国产精品秋霞免费鲁丝片| 一个人免费在线观看的高清视频| 亚洲成人免费电影在线观看| 成人18禁在线播放| 亚洲熟女精品中文字幕| 亚洲人成77777在线视频| 亚洲av成人不卡在线观看播放网| 欧美国产精品va在线观看不卡| 两个人看的免费小视频| 欧美日韩一级在线毛片| 久久精品成人免费网站| 国产一区二区激情短视频| 精品亚洲成国产av| 亚洲精品成人av观看孕妇| 成人影院久久| 欧美人与性动交α欧美精品济南到| 1024香蕉在线观看| 波多野结衣av一区二区av| 黑人欧美特级aaaaaa片| 国产精品成人在线| 丝袜美腿诱惑在线| 久久久久精品人妻al黑| 精品国产一区二区三区四区第35| 动漫黄色视频在线观看| 欧美日韩国产mv在线观看视频| 欧美国产精品一级二级三级| 如日韩欧美国产精品一区二区三区| 在线看a的网站| 热99久久久久精品小说推荐| 亚洲少妇的诱惑av| 欧美精品亚洲一区二区| av国产精品久久久久影院| 一级片'在线观看视频| 亚洲精华国产精华精| 亚洲久久久国产精品| 亚洲少妇的诱惑av| 动漫黄色视频在线观看| 欧美日韩av久久| 9色porny在线观看| 夜夜夜夜夜久久久久| 精品国产乱码久久久久久小说| 久久性视频一级片| 亚洲欧美一区二区三区黑人| 五月开心婷婷网| 1024香蕉在线观看| 啦啦啦免费观看视频1| 蜜桃国产av成人99| 亚洲自偷自拍图片 自拍| 久久99热这里只频精品6学生| 在线播放国产精品三级| 亚洲免费av在线视频| 国产极品粉嫩免费观看在线| 一边摸一边抽搐一进一小说 | 国产亚洲精品一区二区www | 两个人看的免费小视频| 精品久久蜜臀av无| 无人区码免费观看不卡 | 两性午夜刺激爽爽歪歪视频在线观看 | a级片在线免费高清观看视频| 久久精品国产a三级三级三级| 中文字幕制服av| 中文字幕制服av| 在线观看人妻少妇| 国产在线观看jvid| 男女无遮挡免费网站观看| 男女无遮挡免费网站观看| 日韩欧美一区二区三区在线观看 | 一本久久精品| 青青草视频在线视频观看| 久久午夜综合久久蜜桃| 高清在线国产一区| 色综合欧美亚洲国产小说| 日本wwww免费看| 男女高潮啪啪啪动态图| 丁香六月欧美| 99精品在免费线老司机午夜| 亚洲色图综合在线观看| 伦理电影免费视频| 99热网站在线观看| 国产黄频视频在线观看| av超薄肉色丝袜交足视频| 久久人妻福利社区极品人妻图片| 天堂动漫精品| 亚洲第一青青草原| 午夜福利乱码中文字幕| 欧美老熟妇乱子伦牲交| 人人妻人人爽人人添夜夜欢视频| 1024香蕉在线观看| 99re6热这里在线精品视频| 蜜桃在线观看..| 一区二区日韩欧美中文字幕| 久久精品熟女亚洲av麻豆精品| 精品一区二区三卡| 国产在线观看jvid| 日韩成人在线观看一区二区三区| 人人妻,人人澡人人爽秒播| 悠悠久久av| 美女视频免费永久观看网站| 国产在视频线精品| 高潮久久久久久久久久久不卡| 国产亚洲精品一区二区www | 两性夫妻黄色片| 久久久久国内视频| 黄色丝袜av网址大全| av福利片在线| 久久精品亚洲av国产电影网| 免费看十八禁软件| av网站免费在线观看视频| 69av精品久久久久久 | 久久久久久亚洲精品国产蜜桃av| 99热国产这里只有精品6| 波多野结衣一区麻豆| 免费在线观看黄色视频的| 超色免费av| a级毛片在线看网站| 高清黄色对白视频在线免费看| 极品教师在线免费播放| 久久天躁狠狠躁夜夜2o2o| 极品教师在线免费播放| 精品福利永久在线观看| 免费观看av网站的网址| 母亲3免费完整高清在线观看| 夜夜爽天天搞| 波多野结衣一区麻豆| 91国产中文字幕| 亚洲欧美一区二区三区黑人| 久久人妻熟女aⅴ| 99国产综合亚洲精品| 国产成人精品无人区| 午夜福利免费观看在线| 欧美精品一区二区大全| 交换朋友夫妻互换小说| 99国产综合亚洲精品| 久久亚洲真实| 又大又爽又粗| 亚洲精品自拍成人| 青草久久国产| 欧美日韩一级在线毛片| 性少妇av在线| 自拍欧美九色日韩亚洲蝌蚪91| 成人三级做爰电影| 国产亚洲午夜精品一区二区久久| www.精华液| 午夜成年电影在线免费观看| 日本撒尿小便嘘嘘汇集6| 国产亚洲欧美精品永久| 两个人看的免费小视频| 久久精品国产综合久久久| 亚洲中文字幕日韩| 精品亚洲成国产av| 国产欧美日韩一区二区三| 国产精品自产拍在线观看55亚洲 | 99国产精品一区二区蜜桃av | 欧美日韩一级在线毛片| 人妻久久中文字幕网| 亚洲欧美日韩高清在线视频 | 国产精品自产拍在线观看55亚洲 | 午夜福利,免费看| 午夜福利,免费看| 90打野战视频偷拍视频| 叶爱在线成人免费视频播放| 叶爱在线成人免费视频播放| 亚洲精品自拍成人| 人人妻人人澡人人看| 一本—道久久a久久精品蜜桃钙片| 久久久精品94久久精品| 午夜日韩欧美国产| 我要看黄色一级片免费的| 精品福利观看| 久久婷婷成人综合色麻豆| 少妇粗大呻吟视频| 国产精品影院久久| 日本欧美视频一区| 中文字幕人妻丝袜一区二区| 精品国产一区二区三区久久久樱花| 满18在线观看网站| 99riav亚洲国产免费| 99久久国产精品久久久| 久久久国产一区二区| 淫妇啪啪啪对白视频| 女人爽到高潮嗷嗷叫在线视频| 9热在线视频观看99| 久久久精品区二区三区| 大型黄色视频在线免费观看| 国产精品香港三级国产av潘金莲| 中文字幕人妻丝袜一区二区| 国产av精品麻豆| a级毛片在线看网站| 欧美日韩亚洲国产一区二区在线观看 | 搡老熟女国产l中国老女人| 香蕉久久夜色| 热99国产精品久久久久久7| 女人被躁到高潮嗷嗷叫费观| 亚洲国产看品久久| 日本精品一区二区三区蜜桃| 免费观看人在逋| 欧美乱妇无乱码| av片东京热男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三区av网在线观看 | 欧美日韩av久久| 国产无遮挡羞羞视频在线观看| 国产av又大| 99香蕉大伊视频| 国产精品亚洲av一区麻豆| 成人三级做爰电影| 黄片播放在线免费| 亚洲av成人一区二区三| 国产精品久久久久久精品电影小说| 热re99久久精品国产66热6| 丝袜喷水一区| 19禁男女啪啪无遮挡网站| 天天操日日干夜夜撸| 日本一区二区免费在线视频| 成人18禁在线播放| 我要看黄色一级片免费的| 女人被躁到高潮嗷嗷叫费观| 国产高清视频在线播放一区| 亚洲伊人色综图| 亚洲精品久久成人aⅴ小说| 日本黄色日本黄色录像| 亚洲五月婷婷丁香| 国产亚洲欧美精品永久| 99九九在线精品视频| 黄色 视频免费看| 他把我摸到了高潮在线观看 | 91老司机精品| 美女国产高潮福利片在线看| 女人被躁到高潮嗷嗷叫费观| 精品亚洲成a人片在线观看| 女人高潮潮喷娇喘18禁视频| 1024香蕉在线观看| 狠狠婷婷综合久久久久久88av| 香蕉国产在线看| 成人三级做爰电影| 天天躁日日躁夜夜躁夜夜| 在线观看66精品国产| 新久久久久国产一级毛片| 欧美日本中文国产一区发布| 高清毛片免费观看视频网站 | 国产男女内射视频| 青青草视频在线视频观看| 久久精品国产综合久久久| 黄色视频在线播放观看不卡| 老熟女久久久| 久久中文看片网| 国产免费福利视频在线观看| 欧美日韩成人在线一区二区| 一本—道久久a久久精品蜜桃钙片| 久久久久久久久免费视频了| 欧美性长视频在线观看| 亚洲国产成人一精品久久久| 一区福利在线观看| 中文字幕人妻丝袜一区二区| 久久久久久久久久久久大奶| 99精品久久久久人妻精品| 国产男女内射视频| 色综合婷婷激情| 亚洲欧美日韩另类电影网站| 久久人人97超碰香蕉20202| 另类亚洲欧美激情| 成在线人永久免费视频| 国产成人影院久久av| 操美女的视频在线观看| 成人影院久久| 亚洲自偷自拍图片 自拍| 韩国精品一区二区三区| 久久中文看片网| 精品免费久久久久久久清纯 | 国产福利在线免费观看视频| 99精品久久久久人妻精品| 一级毛片电影观看| 国产在视频线精品| 悠悠久久av| 日韩欧美一区二区三区在线观看 | 另类精品久久| 免费少妇av软件| 男女床上黄色一级片免费看| 久久99一区二区三区| 国产成人av教育| 色尼玛亚洲综合影院| 在线观看舔阴道视频| 亚洲精品在线美女| 搡老岳熟女国产| 国产亚洲欧美在线一区二区| 18在线观看网站| 午夜激情久久久久久久| 午夜91福利影院| 视频在线观看一区二区三区| 狠狠狠狠99中文字幕| 免费在线观看黄色视频的| 91大片在线观看| 在线 av 中文字幕| 亚洲伊人色综图| 亚洲熟妇熟女久久| 丝袜人妻中文字幕| 日韩中文字幕视频在线看片| 中文字幕高清在线视频| 伊人久久大香线蕉亚洲五| 午夜福利一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 69精品国产乱码久久久| 国产高清videossex| 热99国产精品久久久久久7| 亚洲 国产 在线| 美女高潮喷水抽搐中文字幕| 高清黄色对白视频在线免费看| 日韩一卡2卡3卡4卡2021年| 脱女人内裤的视频| 午夜两性在线视频| 两个人免费观看高清视频| 最近最新免费中文字幕在线| 日本撒尿小便嘘嘘汇集6| 久久天堂一区二区三区四区| 丝袜美腿诱惑在线| 国产麻豆69| 一级毛片精品| 日本黄色视频三级网站网址 | 亚洲av成人不卡在线观看播放网| 十八禁网站网址无遮挡| 欧美日韩成人在线一区二区| 男女高潮啪啪啪动态图| 国产精品久久久久久精品电影小说| 日本欧美视频一区| 性色av乱码一区二区三区2| 99热网站在线观看| 9191精品国产免费久久| 日韩免费av在线播放| 捣出白浆h1v1| 国产av又大| 国产成人av激情在线播放| 亚洲五月色婷婷综合| 成年动漫av网址| 纵有疾风起免费观看全集完整版| 欧美黄色淫秽网站| 亚洲午夜理论影院| 日韩欧美国产一区二区入口| 亚洲精品粉嫩美女一区| 一区在线观看完整版| 青青草视频在线视频观看| 亚洲国产欧美日韩在线播放| 国产精品国产av在线观看| 啪啪无遮挡十八禁网站| 色综合欧美亚洲国产小说| 又黄又粗又硬又大视频| 久久久国产欧美日韩av| 久久中文字幕人妻熟女| e午夜精品久久久久久久| 成人精品一区二区免费| 国产深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 久久热在线av| av电影中文网址| 午夜福利视频在线观看免费| 亚洲av电影在线进入| 日本撒尿小便嘘嘘汇集6| 色综合婷婷激情| 国产男靠女视频免费网站| 国产一区二区三区在线臀色熟女 | 不卡一级毛片| 99精品欧美一区二区三区四区| 超色免费av| 日本av手机在线免费观看| 精品国产超薄肉色丝袜足j| 美女高潮喷水抽搐中文字幕| 成人黄色视频免费在线看| 亚洲精品久久午夜乱码| 午夜免费成人在线视频| 中文欧美无线码| 国产99久久九九免费精品| 亚洲国产成人一精品久久久| 国产精品.久久久| 国产国语露脸激情在线看| 国产成人精品久久二区二区免费| 精品人妻1区二区| 黄频高清免费视频| 午夜激情久久久久久久| 97人妻天天添夜夜摸| 两个人免费观看高清视频| 久久免费观看电影| 国产在线精品亚洲第一网站| 亚洲精品国产区一区二| 岛国毛片在线播放| 国产国语露脸激情在线看| 午夜福利视频精品| 亚洲少妇的诱惑av| 欧美激情高清一区二区三区| 国产成人av教育| 亚洲中文日韩欧美视频| 人人妻,人人澡人人爽秒播| 精品国产乱码久久久久久男人| 色综合欧美亚洲国产小说| 一边摸一边做爽爽视频免费| 欧美黄色片欧美黄色片| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜免费鲁丝| 成人国产一区最新在线观看| 精品国产乱码久久久久久男人| 亚洲第一欧美日韩一区二区三区 | 午夜视频精品福利| 天天添夜夜摸| 久久精品国产99精品国产亚洲性色 | a在线观看视频网站| 欧美日韩国产mv在线观看视频| 波多野结衣av一区二区av| 国产成人精品无人区| 丰满人妻熟妇乱又伦精品不卡| 国产成人欧美| 亚洲av日韩精品久久久久久密| 女人精品久久久久毛片| 色婷婷久久久亚洲欧美| 在线永久观看黄色视频| 如日韩欧美国产精品一区二区三区| 国产成人免费无遮挡视频| 男女午夜视频在线观看| 久久久久久久久久久久大奶| 我的亚洲天堂| 老熟女久久久| 国产成人精品久久二区二区91| 久久久久精品国产欧美久久久| 飞空精品影院首页| 99国产综合亚洲精品| 一级片'在线观看视频| 亚洲五月色婷婷综合| 国产午夜精品久久久久久| 日本撒尿小便嘘嘘汇集6| 国产精品国产高清国产av | 99久久人妻综合| av超薄肉色丝袜交足视频| 50天的宝宝边吃奶边哭怎么回事| 日韩制服丝袜自拍偷拍| 狠狠狠狠99中文字幕| 黄色a级毛片大全视频| 久久国产精品男人的天堂亚洲| 69精品国产乱码久久久| 欧美日韩av久久| 91成年电影在线观看| 色尼玛亚洲综合影院| 精品人妻1区二区| 天天操日日干夜夜撸| 真人做人爱边吃奶动态| 在线永久观看黄色视频| 99热网站在线观看| 亚洲精品在线美女| 男女免费视频国产| 成人精品一区二区免费| 国产91精品成人一区二区三区 | 乱人伦中国视频| 91国产中文字幕| 超碰成人久久| 日本黄色日本黄色录像| 欧美另类亚洲清纯唯美| 亚洲专区中文字幕在线| av超薄肉色丝袜交足视频| 男男h啪啪无遮挡| 国产免费av片在线观看野外av| 亚洲人成伊人成综合网2020| 在线观看免费视频日本深夜| 俄罗斯特黄特色一大片| 久久久久久久久免费视频了| 久久免费观看电影| 久久久国产成人免费| 成年人午夜在线观看视频| av有码第一页| 国产人伦9x9x在线观看| videos熟女内射| 欧美乱妇无乱码| www日本在线高清视频| 久久久水蜜桃国产精品网| 国产精品亚洲av一区麻豆| 超碰97精品在线观看| 少妇裸体淫交视频免费看高清 | 黄片大片在线免费观看| 咕卡用的链子| kizo精华| 精品久久久久久电影网| 一级黄色大片毛片| 俄罗斯特黄特色一大片| 国产一区二区 视频在线| 少妇精品久久久久久久| 大香蕉久久成人网| 我要看黄色一级片免费的| 99re6热这里在线精品视频| 日韩中文字幕视频在线看片| 欧美日韩黄片免| 久久99一区二区三区| 黄色毛片三级朝国网站| 一区福利在线观看| 精品国产一区二区久久| 黄色成人免费大全| 午夜福利免费观看在线| 正在播放国产对白刺激| 亚洲美女黄片视频| 国产欧美日韩一区二区精品| 亚洲一区中文字幕在线| 国产三级黄色录像| 国产高清videossex| 91成人精品电影| tocl精华| 9热在线视频观看99| 一个人免费在线观看的高清视频| 正在播放国产对白刺激| 国产亚洲精品第一综合不卡| 香蕉国产在线看| 国产精品国产高清国产av | 蜜桃在线观看..| 国产精品亚洲av一区麻豆| 黄色 视频免费看| 亚洲av欧美aⅴ国产| 99九九在线精品视频| 欧美黄色片欧美黄色片| 婷婷丁香在线五月| 中文字幕制服av| 十八禁人妻一区二区| a级毛片黄视频| 国产av精品麻豆| 国精品久久久久久国模美| 80岁老熟妇乱子伦牲交| 十分钟在线观看高清视频www| 国产高清视频在线播放一区| 这个男人来自地球电影免费观看| 高清欧美精品videossex| 五月开心婷婷网| 久久久久久久久免费视频了| h视频一区二区三区| 国产97色在线日韩免费| 欧美日韩亚洲高清精品| 18在线观看网站| 成年人午夜在线观看视频| 精品久久久久久久毛片微露脸| 青青草视频在线视频观看| 欧美黑人精品巨大| 18禁裸乳无遮挡动漫免费视频| 老熟妇仑乱视频hdxx| 高清毛片免费观看视频网站 | 午夜福利视频在线观看免费| 热99re8久久精品国产| 欧美日韩亚洲高清精品| 免费人妻精品一区二区三区视频| 一边摸一边做爽爽视频免费| 欧美亚洲 丝袜 人妻 在线| 黄色成人免费大全| 亚洲av欧美aⅴ国产| 久久久久网色| 国产福利在线免费观看视频| 老司机深夜福利视频在线观看| 国产一区二区三区在线臀色熟女 | 久久精品aⅴ一区二区三区四区| 99久久精品国产亚洲精品| 午夜精品国产一区二区电影| 亚洲五月婷婷丁香| 成人特级黄色片久久久久久久 | 精品午夜福利视频在线观看一区 | 天天添夜夜摸| 热re99久久国产66热| 亚洲人成77777在线视频| 色播在线永久视频| 成人亚洲精品一区在线观看| 久久国产亚洲av麻豆专区| 十八禁网站免费在线| 国产一区二区 视频在线| 亚洲,欧美精品.| 欧美精品av麻豆av| 国产一卡二卡三卡精品| www.熟女人妻精品国产| 在线永久观看黄色视频| 久久久精品94久久精品| 一进一出好大好爽视频| 成人影院久久| 午夜福利,免费看| 看免费av毛片| 亚洲,欧美精品.| 我的亚洲天堂| 十八禁网站网址无遮挡| 看免费av毛片| 国产午夜精品久久久久久| 久久国产亚洲av麻豆专区| 日韩熟女老妇一区二区性免费视频| 中文字幕av电影在线播放| 国产一区二区三区综合在线观看| 久久人妻av系列| 亚洲av美国av| 国产日韩欧美视频二区| 99国产综合亚洲精品| 香蕉久久夜色| 中文字幕人妻丝袜制服| 高清欧美精品videossex| 国产伦人伦偷精品视频| 老司机亚洲免费影院| av网站免费在线观看视频| 一边摸一边抽搐一进一出视频| 男女下面插进去视频免费观看| 激情在线观看视频在线高清 | 国产成人精品无人区| 亚洲五月色婷婷综合| 亚洲五月婷婷丁香| 老司机靠b影院| 国产一区二区在线观看av| 国产精品麻豆人妻色哟哟久久| 美女扒开内裤让男人捅视频| 成年动漫av网址| 成人影院久久| 99在线人妻在线中文字幕 | 免费一级毛片在线播放高清视频 | 欧美精品一区二区大全| 又紧又爽又黄一区二区| 欧美成狂野欧美在线观看| 少妇粗大呻吟视频| 亚洲国产av新网站| 亚洲精品在线观看二区| 日本撒尿小便嘘嘘汇集6| 少妇精品久久久久久久| 高清av免费在线|