• <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)勢樹種凋落物葉分解歷程中的水文特征
    91久久精品国产一区二区成人| 嫩草影院新地址| 在线免费观看的www视频| 日韩欧美精品v在线| 一区二区三区高清视频在线| 欧美日韩一区二区视频在线观看视频在线 | 赤兔流量卡办理| 一区二区三区高清视频在线| 日韩国内少妇激情av| 伦精品一区二区三区| 国产在线男女| 成年女人在线观看亚洲视频 | 国产一区有黄有色的免费视频 | 国产白丝娇喘喷水9色精品| 熟妇人妻久久中文字幕3abv| 青春草视频在线免费观看| 国产综合懂色| 九九爱精品视频在线观看| 晚上一个人看的免费电影| 欧美 日韩 精品 国产| 亚洲av不卡在线观看| 欧美成人精品欧美一级黄| 一级av片app| 又爽又黄无遮挡网站| 成年免费大片在线观看| 亚洲精品第二区| 在线免费观看不下载黄p国产| av女优亚洲男人天堂| 亚洲综合精品二区| 18禁在线无遮挡免费观看视频| 成人午夜高清在线视频| 日本爱情动作片www.在线观看| 国产亚洲最大av| 亚洲最大成人av| 色综合站精品国产| 久久久精品94久久精品| 国产男人的电影天堂91| 高清日韩中文字幕在线| 777米奇影视久久| 久久人人爽人人爽人人片va| 国产欧美另类精品又又久久亚洲欧美| 国产欧美另类精品又又久久亚洲欧美| 亚洲18禁久久av| 麻豆国产97在线/欧美| 久久精品综合一区二区三区| 欧美精品国产亚洲| 黑人高潮一二区| xxx大片免费视频| 一级a做视频免费观看| 丰满乱子伦码专区| 日韩av免费高清视频| 日韩三级伦理在线观看| 亚洲av男天堂| 国产成人免费观看mmmm| 一级片'在线观看视频| 国产真实伦视频高清在线观看| 寂寞人妻少妇视频99o| 国产精品久久久久久久电影| 国产男女超爽视频在线观看| 国产精品蜜桃在线观看| 免费无遮挡裸体视频| 日韩成人av中文字幕在线观看| 亚洲内射少妇av| 舔av片在线| 天堂√8在线中文| 69人妻影院| 男人和女人高潮做爰伦理| 在线 av 中文字幕| 欧美丝袜亚洲另类| 最近的中文字幕免费完整| 高清在线视频一区二区三区| 少妇熟女aⅴ在线视频| 2021天堂中文幕一二区在线观| 2021天堂中文幕一二区在线观| 丰满人妻一区二区三区视频av| 女人十人毛片免费观看3o分钟| 免费黄色在线免费观看| 可以在线观看毛片的网站| 亚洲第一区二区三区不卡| 黄片wwwwww| 国产精品精品国产色婷婷| 非洲黑人性xxxx精品又粗又长| 成人一区二区视频在线观看| 日本免费在线观看一区| 国产成年人精品一区二区| 亚洲精品国产av成人精品| av黄色大香蕉| 我的老师免费观看完整版| 亚洲av免费在线观看| 亚洲精品国产av成人精品| 最近最新中文字幕大全电影3| 夫妻性生交免费视频一级片| 日本wwww免费看| 日韩成人av中文字幕在线观看| 韩国高清视频一区二区三区| 亚洲精华国产精华液的使用体验| 国产黄色免费在线视频| 久久99热这里只频精品6学生| 国产一级毛片在线| 亚洲性久久影院| 国产精品福利在线免费观看| 美女xxoo啪啪120秒动态图| 久久久久性生活片| 免费大片黄手机在线观看| 美女主播在线视频| 一级毛片我不卡| 老师上课跳d突然被开到最大视频| 国产老妇女一区| 99热网站在线观看| 国产精品久久久久久精品电影小说 | 精品一区二区免费观看| 国内精品宾馆在线| 亚洲欧美成人精品一区二区| 自拍偷自拍亚洲精品老妇| 国产精品女同一区二区软件| 女人被狂操c到高潮| 熟妇人妻久久中文字幕3abv| 婷婷色综合大香蕉| 国精品久久久久久国模美| 亚洲人成网站在线播| 网址你懂的国产日韩在线| 草草在线视频免费看| 久久午夜福利片| 国产亚洲av片在线观看秒播厂 | 超碰av人人做人人爽久久| 亚洲精品456在线播放app| 久久鲁丝午夜福利片| 精品久久久久久久久亚洲| 女人十人毛片免费观看3o分钟| 九九爱精品视频在线观看| .国产精品久久| av免费在线看不卡| 亚洲人与动物交配视频| 成人高潮视频无遮挡免费网站| 国产不卡一卡二| 美女内射精品一级片tv| 国产欧美日韩精品一区二区| 亚洲av中文av极速乱| 国产高清国产精品国产三级 | 成人性生交大片免费视频hd| 日韩一区二区视频免费看| 亚洲av国产av综合av卡| 国产免费视频播放在线视频 | 一级a做视频免费观看| 亚洲精品第二区| 国产成人一区二区在线| 亚洲av男天堂| 久久人人爽人人片av| 日韩 亚洲 欧美在线| 一级毛片我不卡| 青青草视频在线视频观看| 色哟哟·www| 美女cb高潮喷水在线观看| 99热这里只有是精品50| 国产黄频视频在线观看| 亚洲综合精品二区| 可以在线观看毛片的网站| 国产亚洲5aaaaa淫片| 啦啦啦中文免费视频观看日本| 亚洲美女搞黄在线观看| 在线观看av片永久免费下载| 听说在线观看完整版免费高清| 一级av片app| 亚洲av二区三区四区| 97超视频在线观看视频| 久久99热6这里只有精品| 亚洲欧美日韩东京热| 搡女人真爽免费视频火全软件| 午夜免费激情av| 美女被艹到高潮喷水动态| 国产精品国产三级专区第一集| 久久99热这里只频精品6学生| 亚洲精品一区蜜桃| 国产黄色免费在线视频| 日韩欧美精品免费久久| 亚洲精华国产精华液的使用体验| 国产一区二区在线观看日韩| 插阴视频在线观看视频| 亚洲无线观看免费| 亚州av有码| 久久精品夜色国产| 在线a可以看的网站| 国产一区亚洲一区在线观看| 最近2019中文字幕mv第一页| 欧美+日韩+精品| 草草在线视频免费看| 老司机影院成人| 尾随美女入室| 亚洲精品乱久久久久久| 干丝袜人妻中文字幕| 亚洲精品一二三| 日本-黄色视频高清免费观看| 国产毛片a区久久久久| 国产精品一区www在线观看| 亚洲一区高清亚洲精品| 中文字幕亚洲精品专区| av专区在线播放| 成人亚洲精品一区在线观看 | 日韩av免费高清视频| 最近最新中文字幕免费大全7| 国产综合精华液| 97精品久久久久久久久久精品| 国产成人a区在线观看| 色5月婷婷丁香| 一个人看的www免费观看视频| 在线观看一区二区三区| 我要看日韩黄色一级片| 丝袜美腿在线中文| 国产精品久久视频播放| 国产一区二区亚洲精品在线观看| 免费大片18禁| 久久精品夜色国产| 精品国内亚洲2022精品成人| 午夜福利在线观看免费完整高清在| 国内揄拍国产精品人妻在线| 两个人视频免费观看高清| 麻豆av噜噜一区二区三区| 久久久久久久久中文| 国产男女超爽视频在线观看| 美女大奶头视频| 美女xxoo啪啪120秒动态图| 麻豆乱淫一区二区| 成年女人在线观看亚洲视频 | 午夜福利网站1000一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 一二三四中文在线观看免费高清| 又爽又黄a免费视频| 婷婷色麻豆天堂久久| 丰满少妇做爰视频| 黄色一级大片看看| 麻豆国产97在线/欧美| 亚洲人成网站在线播| 久久久亚洲精品成人影院| 一个人观看的视频www高清免费观看| 欧美潮喷喷水| 婷婷六月久久综合丁香| 少妇熟女欧美另类| 高清午夜精品一区二区三区| 免费大片18禁| 久久久欧美国产精品| 国产精品久久久久久精品电影| 免费在线观看成人毛片| av在线蜜桃| 卡戴珊不雅视频在线播放| 久久这里有精品视频免费| 热99在线观看视频| 亚洲精品456在线播放app| 午夜福利网站1000一区二区三区| 国产成人a∨麻豆精品| 久久久久久久久中文| 色吧在线观看| 狠狠精品人妻久久久久久综合| 成人毛片a级毛片在线播放| 免费观看无遮挡的男女| 看黄色毛片网站| 亚洲国产精品专区欧美| 亚洲精品第二区| 一级毛片黄色毛片免费观看视频| 国产高清有码在线观看视频| 亚洲三级黄色毛片| 菩萨蛮人人尽说江南好唐韦庄| 国产乱来视频区| 肉色欧美久久久久久久蜜桃 | 亚洲一区高清亚洲精品| 国产黄色免费在线视频| 能在线免费观看的黄片| 亚洲欧美精品自产自拍| 欧美日韩国产mv在线观看视频 | 可以在线观看毛片的网站| 日韩制服骚丝袜av| 中文字幕亚洲精品专区| 亚洲成人精品中文字幕电影| 午夜爱爱视频在线播放| av福利片在线观看| 日韩av在线大香蕉| 熟妇人妻久久中文字幕3abv| 久久久久久国产a免费观看| 伊人久久国产一区二区| 久久这里只有精品中国| 亚洲人成网站在线播| 免费看a级黄色片| 男人爽女人下面视频在线观看| 国产黄片美女视频| 国精品久久久久久国模美| 欧美高清成人免费视频www| 日韩制服骚丝袜av| 在线 av 中文字幕| 天堂影院成人在线观看| 亚洲精品第二区| 国产伦精品一区二区三区四那| 亚洲精品乱码久久久久久按摩| 久久综合国产亚洲精品| 久久这里有精品视频免费| 99久久中文字幕三级久久日本| 午夜福利成人在线免费观看| av在线观看视频网站免费| 亚洲国产色片| 日本黄大片高清| 六月丁香七月| 国产白丝娇喘喷水9色精品| 美女xxoo啪啪120秒动态图| 欧美潮喷喷水| 五月玫瑰六月丁香| 伊人久久国产一区二区| eeuss影院久久| 国产亚洲一区二区精品| 国产精品国产三级专区第一集| 中文乱码字字幕精品一区二区三区 | 最新中文字幕久久久久| 女的被弄到高潮叫床怎么办| 日产精品乱码卡一卡2卡三| 好男人在线观看高清免费视频| 免费av不卡在线播放| 嫩草影院新地址| 亚洲精华国产精华液的使用体验| 男女边吃奶边做爰视频| 黄色欧美视频在线观看| 我的女老师完整版在线观看| 日本三级黄在线观看| av国产免费在线观看| 久久人人爽人人片av| 久久久久久九九精品二区国产| 午夜免费激情av| 人人妻人人澡人人爽人人夜夜 | 国产又色又爽无遮挡免| 成人特级av手机在线观看| 国产淫语在线视频| 精品久久久久久久久av| 亚洲丝袜综合中文字幕| 久久久久久国产a免费观看| 欧美97在线视频| 国产成人免费观看mmmm| 中文字幕人妻熟人妻熟丝袜美| 久久精品夜夜夜夜夜久久蜜豆| 蜜桃亚洲精品一区二区三区| 国产黄频视频在线观看| 亚洲一区高清亚洲精品| 亚洲乱码一区二区免费版| 国产黄色视频一区二区在线观看| 中文天堂在线官网| 欧美xxxx性猛交bbbb| 男的添女的下面高潮视频| 精品欧美国产一区二区三| 欧美日韩一区二区视频在线观看视频在线 | 男女边吃奶边做爰视频| 午夜视频国产福利| 亚洲精品国产成人久久av| 午夜免费观看性视频| 精品午夜福利在线看| 十八禁国产超污无遮挡网站| 免费av毛片视频| 亚洲精品乱码久久久v下载方式| 亚洲欧美日韩东京热| 嫩草影院精品99| 欧美变态另类bdsm刘玥| 日韩强制内射视频| 久久亚洲国产成人精品v| 久久久久久久午夜电影| 午夜视频国产福利| 亚洲av电影在线观看一区二区三区 | 精品少妇黑人巨大在线播放| 搡老乐熟女国产| 人人妻人人看人人澡| 亚洲精品国产成人久久av| 久久久精品免费免费高清| 有码 亚洲区| 久久精品熟女亚洲av麻豆精品 | 伊人久久国产一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩视频在线欧美| 国产国拍精品亚洲av在线观看| 亚洲人成网站高清观看| 又粗又硬又长又爽又黄的视频| 日本一本二区三区精品| 精品午夜福利在线看| 国产成人精品一,二区| 亚州av有码| 一级爰片在线观看| 2018国产大陆天天弄谢| 99久久精品一区二区三区| 精品久久国产蜜桃| 国产精品不卡视频一区二区| 我要看日韩黄色一级片| 国产探花极品一区二区| 两个人视频免费观看高清| 免费大片黄手机在线观看| 69av精品久久久久久| 精品酒店卫生间| av卡一久久| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 免费观看在线日韩| 国产精品美女特级片免费视频播放器| ponron亚洲| 乱系列少妇在线播放| 午夜福利高清视频| 韩国av在线不卡| 亚洲真实伦在线观看| 麻豆av噜噜一区二区三区| 亚洲欧美日韩无卡精品| 熟妇人妻不卡中文字幕| 一级毛片 在线播放| 国产精品麻豆人妻色哟哟久久 | 欧美人与善性xxx| 观看美女的网站| 少妇熟女aⅴ在线视频| 少妇猛男粗大的猛烈进出视频 | 久久久久久久久大av| 国产在视频线在精品| 国产成人免费观看mmmm| 日韩电影二区| 嫩草影院精品99| 特级一级黄色大片| 男女视频在线观看网站免费| 国产高清不卡午夜福利| 在线 av 中文字幕| 汤姆久久久久久久影院中文字幕 | 中文字幕av成人在线电影| 久久精品夜色国产| 国产 一区精品| 欧美日韩亚洲高清精品| 亚洲av成人精品一二三区| 国产在视频线在精品| 亚洲精品色激情综合| 成年女人在线观看亚洲视频 | 国产精品一二三区在线看| 美女大奶头视频| 九九爱精品视频在线观看| 一级毛片 在线播放| 成人综合一区亚洲| 成人特级av手机在线观看| 亚洲无线观看免费| 少妇猛男粗大的猛烈进出视频 | 99久久中文字幕三级久久日本| 亚洲在线观看片| 国国产精品蜜臀av免费| 亚洲综合色惰| 最近最新中文字幕大全电影3| 国产精品人妻久久久影院| 十八禁国产超污无遮挡网站| 色综合站精品国产| 尾随美女入室| 亚洲高清免费不卡视频| av在线蜜桃| 免费在线观看成人毛片| 在线观看一区二区三区| 亚洲人成网站高清观看| 亚洲最大成人av| 搞女人的毛片| 亚洲婷婷狠狠爱综合网| 老司机影院毛片| 久久6这里有精品| 尾随美女入室| 97在线视频观看| 亚洲在线观看片| 免费电影在线观看免费观看| 久久久久久久久久久丰满| 亚洲国产精品sss在线观看| 国产中年淑女户外野战色| 国产午夜福利久久久久久| 久久韩国三级中文字幕| 黄色配什么色好看| 男女国产视频网站| 国产免费又黄又爽又色| 久久久久久久久大av| 又爽又黄无遮挡网站| 22中文网久久字幕| 国产亚洲精品久久久com| 久久久久国产网址| av一本久久久久| 精品国内亚洲2022精品成人| 我的女老师完整版在线观看| 精品久久久久久成人av| 国产午夜精品久久久久久一区二区三区| 免费黄色在线免费观看| 永久免费av网站大全| 一级毛片我不卡| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频 | 最近中文字幕2019免费版| 嫩草影院新地址| 亚洲无线观看免费| or卡值多少钱| 免费av观看视频| 国产黄色小视频在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美中文字幕日韩二区| 街头女战士在线观看网站| 免费黄网站久久成人精品| 国产成年人精品一区二区| 亚洲精品一二三| 男女视频在线观看网站免费| 中文天堂在线官网| 中文字幕制服av| 久久久久久久久久黄片| 欧美日本视频| 18禁在线无遮挡免费观看视频| 成人毛片60女人毛片免费| 久久精品久久久久久久性| www.色视频.com| 天堂av国产一区二区熟女人妻| 亚洲av一区综合| 成人亚洲精品一区在线观看 | 夫妻午夜视频| 久久鲁丝午夜福利片| 99久国产av精品国产电影| 国产av在哪里看| 亚洲精品一区蜜桃| 国产亚洲最大av| 亚洲色图av天堂| 亚洲欧美中文字幕日韩二区| 看非洲黑人一级黄片| 国产成人精品一,二区| 偷拍熟女少妇极品色| 国产精品1区2区在线观看.| 国产精品一区www在线观看| 日韩中字成人| 国产淫语在线视频| 成人午夜高清在线视频| 边亲边吃奶的免费视频| 午夜福利成人在线免费观看| 九色成人免费人妻av| 国产一区二区三区综合在线观看 | 国产高清有码在线观看视频| 亚洲av日韩在线播放| 亚洲成人一二三区av| 久久久久久久久久黄片| 少妇被粗大猛烈的视频| 日韩电影二区| 丝袜美腿在线中文| 22中文网久久字幕| 九九爱精品视频在线观看| 男女边吃奶边做爰视频| h日本视频在线播放| 人妻少妇偷人精品九色| 在线观看一区二区三区| 国产激情偷乱视频一区二区| 欧美zozozo另类| 永久免费av网站大全| 亚洲国产成人一精品久久久| 国产成人福利小说| 久久久久久久久久久免费av| 三级毛片av免费| 国产一区二区在线观看日韩| 成人欧美大片| 国产精品久久久久久av不卡| 午夜日本视频在线| 久久精品国产鲁丝片午夜精品| 国产精品爽爽va在线观看网站| 国产精品一区二区在线观看99 | 人妻夜夜爽99麻豆av| 女的被弄到高潮叫床怎么办| 色哟哟·www| 人人妻人人澡人人爽人人夜夜 | 三级国产精品欧美在线观看| 欧美成人精品欧美一级黄| 国产男人的电影天堂91| 不卡视频在线观看欧美| 国产老妇伦熟女老妇高清| 美女主播在线视频| a级一级毛片免费在线观看| 亚洲精品影视一区二区三区av| 日韩av不卡免费在线播放| 久久精品国产鲁丝片午夜精品| 激情 狠狠 欧美| 夜夜看夜夜爽夜夜摸| 免费看av在线观看网站| 国产精品一区二区三区四区久久| 亚洲欧美日韩东京热| 人妻少妇偷人精品九色| 久久99热这里只有精品18| 欧美日本视频| 在现免费观看毛片| 日本午夜av视频| 成年女人在线观看亚洲视频 | 又爽又黄a免费视频| 偷拍熟女少妇极品色| 精品一区二区三区人妻视频| 亚洲av电影不卡..在线观看| 国产精品日韩av在线免费观看| 网址你懂的国产日韩在线| 久久精品国产亚洲网站| 国产一区亚洲一区在线观看| 国产成人精品福利久久| 秋霞在线观看毛片| 国产精品人妻久久久影院| 少妇高潮的动态图| 久久人人爽人人片av| 97超视频在线观看视频| 九草在线视频观看| 非洲黑人性xxxx精品又粗又长| 国产精品国产三级国产av玫瑰| 亚洲欧美日韩东京热| 国产黄色视频一区二区在线观看| 国产成人免费观看mmmm| 国产一区二区三区综合在线观看 | 亚洲四区av| 黄片wwwwww| 亚洲精华国产精华液的使用体验| 国产亚洲精品久久久com| 欧美bdsm另类| 亚洲国产日韩欧美精品在线观看| 国产精品嫩草影院av在线观看| 精品久久久久久久久亚洲| 久久久久久久久中文| 国产亚洲精品久久久com| 久久精品久久久久久噜噜老黄| 亚洲精品一二三| 久久久久网色| 久久99热这里只有精品18| 特级一级黄色大片| 国产黄色视频一区二区在线观看| 丰满人妻一区二区三区视频av| 欧美激情久久久久久爽电影| 久久久久性生活片|