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

    室內(nèi)模擬水分管理對(duì)土壤溶液磷質(zhì)量濃度影響研究

    2021-09-02 08:08:26吳瑞胡正義孫藝齊耿紫琪崔巖山
    灌溉排水學(xué)報(bào) 2021年8期
    關(guān)鍵詞:土壤溶液總鐵下層

    吳瑞,胡正義,,孫藝齊,耿紫琪,崔巖山,2*

    室內(nèi)模擬水分管理對(duì)土壤溶液磷質(zhì)量濃度影響研究

    吳瑞1,胡正義1,3,孫藝齊1,耿紫琪3,崔巖山1,2*

    (1.中國(guó)科學(xué)院大學(xué) 資源與環(huán)境學(xué)院,北京 101408;2.中國(guó)科學(xué)院 生態(tài)環(huán)境研究中心,北京 100085;3.中國(guó)科學(xué)院大學(xué) 中丹學(xué)院,北京 100190)

    探究不同水分管理模式下、不同土層土壤溶液磷的環(huán)境行為。本研究采取土柱試驗(yàn)研究了3種水分管理模式下(模式1—模式3土壤含水率均恒定分別為:50%、75%、100%)84 d剖面土壤(表層0~10 cm、中層10~20 cm、下層20~30 cm)溶液磷質(zhì)量濃度變化規(guī)律及其與土壤溶液中TOC、總鐵、總鋁、總鈣之間關(guān)系。水分管理對(duì)土壤溶液磷質(zhì)量濃度有顯著的影響,并與土層有關(guān);表層土壤表現(xiàn)為:模式3(1.27 mg/L)>模式2(0.82 mg/L);中層土壤表現(xiàn)為:模式3(1.60 mg/L)>模式2(1.40 mg/L);下層土壤表現(xiàn)為:模式3(0.31 mg/L)>模式2(0.17 mg/L)>模式1(0.13 mg/L)。逐步回歸分析揭示土壤溶液總磷質(zhì)量濃度與土壤溶液TOC、總鐵、總鋁、總鈣有關(guān),并受水分管理模式和土壤層深度的影響。室內(nèi)模擬3種水分管理模式下土壤溶液總磷質(zhì)量濃度都大于地表水富營(yíng)養(yǎng)化標(biāo)準(zhǔn)值(0.1 mg/L)。退耕還濕期間土壤磷釋放可能會(huì)引起上覆水磷質(zhì)量濃度突然升高,其造成的生態(tài)環(huán)境問題值得關(guān)注。

    水分管理;土壤溶液;磷;三江平原;退耕還濕

    0 引言

    【研究意義】三江平原地區(qū)作為我國(guó)三大濕地之一[1]。歷史上,三江平原的開發(fā)曾為滿足國(guó)家糧食需求、保障國(guó)家糧食安全做出過重要貢獻(xiàn)[2]。但是,1980—2015年三江平原濕地沼澤面積減少了53.27%[3],減少的濕地主要轉(zhuǎn)化為耕地,耕地面積比從1954年15.7%增加到了2015年的56.5%[4]。三江平原濕地退化嚴(yán)重,2000—2015年濕地植被覆蓋度由91.8%減小至74%[5]。長(zhǎng)期的農(nóng)業(yè)開墾和農(nóng)業(yè)水利工程建設(shè),使三江平原地區(qū)濕地面積逐漸萎縮,旱澇災(zāi)害頻率增加,水土流失加劇。農(nóng)業(yè)堤防排水渠改變了濕地的水文條件,導(dǎo)致濕地C、N、P循環(huán)過程發(fā)生極大改變[6]。為此,三江平原自然保護(hù)區(qū)于2000年被列為國(guó)家級(jí)自然保護(hù)區(qū),隨后退耕還濕工程力度不斷加大。濕地周邊農(nóng)田停止農(nóng)業(yè)耕種后,生態(tài)補(bǔ)水是退耕還濕工程重要技術(shù)方法。有研究表明,水分管理可以顯著地影響土壤氮磷養(yǎng)分釋放[7],進(jìn)而影響植物生長(zhǎng)[8-12]及水體氮磷質(zhì)量濃度?!狙芯窟M(jìn)展】田娟等[13]研究淹水土壤磷釋放發(fā)現(xiàn)上覆水可溶磷最大值可達(dá)到0.16 mg/L,大于地表水富營(yíng)養(yǎng)化標(biāo)準(zhǔn)值0.1 mg/L[14]。有關(guān)淹水及干濕交替磷釋放機(jī)理研究也較多。崔虎等[15]通過室內(nèi)模擬干濕交替發(fā)現(xiàn)Fe-P是土壤釋放磷的主要組分,Al-P也有較少釋放。周健等[16]發(fā)現(xiàn)水庫(kù)消落帶持續(xù)淹水土壤Al-P量隨時(shí)間先降低后升高,F(xiàn)e-P量隨淹水時(shí)間先升高后降低。【切入點(diǎn)】雖然水分對(duì)土壤磷釋放研究取得了長(zhǎng)足進(jìn)展,但土壤溶液磷能更好反映磷的環(huán)境效應(yīng),下層土壤溶液磷可以影響地下水水質(zhì)、中層土壤溶液磷可以影響植物生長(zhǎng)、表層土壤溶液磷對(duì)地表水水體有重要影響。然而,不同水分模式下剖面土壤溶液中磷變化規(guī)律還缺乏研究,制約了退耕還濕過程中土壤磷釋放對(duì)水質(zhì)影響和植物磷供需的評(píng)價(jià)。【擬解決的關(guān)鍵問題】本研究以黑龍江省三江平原國(guó)家級(jí)濕地自然保護(hù)區(qū)內(nèi)的撫遠(yuǎn)市烏蘇鎮(zhèn)農(nóng)田土壤為對(duì)象,采用室內(nèi)土柱模擬實(shí)驗(yàn),研究水分管理模式對(duì)剖面土壤溶液磷變化規(guī)律,以期為退耕還濕工程設(shè)計(jì)提供參考。

    1 材料和方法

    1.1 供試土壤

    供試土壤采自黑龍江省撫遠(yuǎn)市烏蘇里鎮(zhèn)某農(nóng)田(48° 4'N,134°32 E),農(nóng)田種植的主要作物是玉米。選擇典型采樣地,挖土壤剖面,分層采集0~10、10~20、20~30 cm土壤,分別放到袋子里密封保存,帶回實(shí)驗(yàn)室內(nèi)風(fēng)干,挑去秸稈殘留,過1 mm篩。供試驗(yàn)用。留部分過1 mm篩土壤分析有效磷、pH值、機(jī)械組成。取部分過1 mm篩土壤樣品,再磨碎,過100目篩,供有機(jī)質(zhì)、總磷分析。供實(shí)驗(yàn)土壤為東北黑土;表層和中層土壤屬于粉砂質(zhì)壤土,下層土壤屬于粉砂質(zhì)黏壤土;土壤總磷量很高,但是有效磷低。土壤基本理化性質(zhì)見表1。

    表1 土壤基本理化性質(zhì)

    1.2 土柱模擬裝置及運(yùn)行

    試驗(yàn)?zāi)M裝置如圖1所示,土柱為PVC管,直徑160 mm、高450 mm。首先向每個(gè)土柱底部墊入脫脂棉與面紗布,然后放入洗凈石英砂(直徑1~2 mm)1 kg,隨后依次向每根柱子裝填下層土壤(20~30 cm)1.2 kg,中層土壤(10~20 cm)1 kg,表層土壤(0~10cm)1.4 kg土壤。水位探頭放置于觀測(cè)管中,水位控制系統(tǒng)通過膠管與觀測(cè)管相連,通過設(shè)置不同位置的水位探頭來控制土壤中的含水率。當(dāng)土柱中水位低于探頭探測(cè)的最低水位時(shí),水泵啟動(dòng),水泵通過水管自動(dòng)向土柱中注水,直到水位達(dá)到設(shè)定水位為止。土柱外表安裝有3個(gè)相距10 cm的土壤溶液采樣器[17],用于采集不同深度的土壤溶液,供分析。土柱底部設(shè)有裝置室,裝置室底部設(shè)有底座。在沖洗液收集裝置與隔離裝置連接的管路上設(shè)有開關(guān)閥。水位控制系統(tǒng)包括液位探頭、水泵、繼電器和接觸器;液位探頭通過電纜與繼電器連接。

    1.試驗(yàn)?zāi)M柱;2.淋洗液收集裝置;3.觀測(cè)管;4.橡膠軟管;5.土培室;6.設(shè)備室;7.底座;8.隔離裝置;9.開關(guān)閥;10.液位探頭;11.水泵;12.水箱;13.土壤溶液采樣器

    1.3 試驗(yàn)處理與采樣

    設(shè)3種水分管理模式,模式1為土壤質(zhì)量含水率恒定50%、模式2為土壤質(zhì)量含水率恒定75%、模式3為土壤質(zhì)量含水率恒定100%。預(yù)實(shí)驗(yàn)測(cè)定實(shí)現(xiàn)土柱含水率50%、75%、100%分別對(duì)應(yīng)水位-37、-22、-5 cm(以土柱外壁最高點(diǎn)為相對(duì)高程0點(diǎn),以豎直向上為正方向)。所以將水位探頭分別設(shè)置位距離土柱外壁最高點(diǎn)37、22、5 cm處(最低水位為以探頭中心點(diǎn)為原點(diǎn)-2.5 cm、最高水位為+2.5 cm)。每個(gè)處理設(shè)3個(gè)重復(fù)。試驗(yàn)期間室溫為23℃。土壤溶液采樣時(shí)間前期5 d左右1次,28 d以后大約14 d為1次。

    1.4 分析與測(cè)試方法

    用膜孔為0.45 μm土壤溶液取樣器(Rhizon 19.60.22 F)收集土壤溶液[18]。土壤經(jīng)雙氧水消化,激光粒度儀(Mastersizer 2000)測(cè)定土壤粒徑[19];ICP-OES(Optima 5300 DV )測(cè)定土壤溶液總磷、總鋁、總鐵、總鈣[20];TOC分析儀(Multi N/C 3100)測(cè)得土壤溶液中TOC[20];土壤溶液pH值用pH計(jì)(上海梅特勒-托利多公司SG2)測(cè)定[20];土壤溶液采用便攜數(shù)字化多參數(shù)分析儀[20](HQ40d HACH);鉬銻抗分光光度法測(cè)定土壤溶液中無機(jī)磷[20];土壤pH值采用pH計(jì)(上海梅特勒-托利多公司SG2)測(cè)定(土、水比為1∶2.5)[19]。用HCl(0.025 mol/L)+NH4F(0.03 mol/L),浸提土壤有效磷,鉬銻抗比色法測(cè)定[19];土壤有機(jī)質(zhì)采用高溫外熱重鉻酸鉀氧化法-滴定法測(cè)定[19];土壤經(jīng)過微波消解儀(美國(guó)CEM公司Mars 6)(硝酸、雙氧水、氫氟酸比為5∶3∶2)消解,通過ICP-OES(Optima 5300 DV )測(cè)定總磷、總鐵[19];土壤無機(jī)磷組分采用酸性土壤無機(jī)磷形態(tài)的分級(jí)測(cè)定[19]。

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

    采用AutoCAD 2012 (Autodesk) 繪制土柱裝置示意圖。采用Excel 2019(Microsoft)繪制參數(shù)動(dòng)態(tài)變化圖。所有測(cè)定指標(biāo)應(yīng)用SPSS 22(IBM, USA)軟件中Wilcxon檢驗(yàn)方法分析處理間差異顯著性水平,并應(yīng)用該軟件分析土壤溶液總磷與總鐵、總鋁、總鈣、TOC之間相關(guān)性。用SPSS 22(IBM, USA)多元逐步回歸擬合土壤溶液總磷與土壤溶液總鐵、總鋁、總鈣、TOC之間關(guān)系。

    2 結(jié)果與分析

    2.1 不同水分管理模式下土壤溶液pH值和Eh的變化

    模式1表層和中層土壤沒有收集到土壤溶液,下層土壤收集到土壤溶液,模式2和模式3各層土壤都收集到土壤溶液。3種水分管理模式下土柱中土壤溶液pH值、變化見圖2。3種水分管理模式下,表層、中層土壤溶液pH值表現(xiàn)為:模式3>模式2;下層表現(xiàn)為:模式3>模式2>模式1。各土層土壤溶液pH值:模式3表現(xiàn)為:表層>中層>下層,模式2表現(xiàn)為:中層>下層>表層。

    3種水分管理模式下,表層和中層土壤溶液表現(xiàn)為,模式3<模式2;下層土壤溶液表現(xiàn)為,模式3<模式2<模式1;各土層土壤溶液:模式3表現(xiàn)為:下層>中層>表層,模式2表現(xiàn)為:下層>表層>中層,表層和中層土壤中,培養(yǎng)到20 d以后差別更加明顯,而下層土壤直到40 d以后才明顯。

    在培養(yǎng)期間土壤溶液pH值與與水分管理模式和土壤層次有關(guān)且呈相反的變化規(guī)律,pH值除了模式3表層外在20 d以后呈緩慢上升趨勢(shì),則相反呈下降趨勢(shì)。模式2表層pH值與在20 d前沒有穩(wěn)定,20 d后pH值下降第70 d開始再次上升,則波動(dòng)上升。

    圖2 3種水分管理模式下土柱中土壤溶液pH值和Eh動(dòng)態(tài)變化

    2.2 不同水分管理模式下土壤溶液TOC、總鐵、總磷和總鈣質(zhì)量濃度的變化

    3種水分管理模式下土柱中土壤溶液TOC、總鐵、總鋁和總鈣質(zhì)量濃度的變化見圖3。由圖3可見,水分管理模式對(duì)土壤溶液TOC、總鐵、總鋁和總鈣質(zhì)量濃度有顯著影響,并與土層和時(shí)間有關(guān)。表層和中層土壤溶液TOC、總鐵、總鋁和總鈣質(zhì)量濃度表現(xiàn)為:模式3>模式2;而下層土壤溶液TOC、總鐵和總鈣質(zhì)量濃度表現(xiàn)為:模式3>模式2>模式1,下層土壤溶液總鋁質(zhì)量濃度表現(xiàn)為:模式3>模式1>模式2。模式2和模式3下,各土層土壤溶液TOC、總鐵、總鋁和總鈣質(zhì)量濃度表現(xiàn)為:中層>表層>下層。培養(yǎng)期間,土壤TOC、總鐵、總鋁和總鈣質(zhì)量濃度動(dòng)態(tài)變化規(guī)律與水分管理模式、土層次有關(guān)。隨培養(yǎng)時(shí)間延長(zhǎng),表層土壤溶液鋁質(zhì)量濃度逐漸增加,而表層土壤溶液TOC、總鐵、總鈣質(zhì)量濃度先升高,在40~54 d達(dá)到最大值,隨后逐漸下降(圖3)。隨培養(yǎng)時(shí)間延長(zhǎng),下層土壤溶液TOC、總鐵、總鈣質(zhì)量濃度逐漸升高,而下層土壤溶液總鋁質(zhì)量濃度先逐漸升高,在16~21 d達(dá)到最大值,隨后逐漸下降。隨培養(yǎng)時(shí)間延長(zhǎng),中層土壤溶液總鐵、鈣濃質(zhì)量濃度逐漸升高,達(dá)到最大值(總鐵在40~54 d,總鈣在21~38 d)后,再逐漸下降。隨培養(yǎng)時(shí)間延長(zhǎng),模式2中層土壤溶液TOC和總鋁質(zhì)量濃度逐漸升高,達(dá)到最大值(TOC在21 d,總鋁在11 d)后,再逐漸下降;模式2中層土壤溶液TOC質(zhì)量濃度逐漸下降,而模式3中層土壤溶液總鋁質(zhì)量濃度趨于逐漸升高。

    圖3 3種水分管理模式下土壤溶液TOC、鐵、鋁、鈣動(dòng)態(tài)變化

    表2 3種水分管理模式下土壤溶液理化性質(zhì)

    注 表中數(shù)值以整個(gè)培養(yǎng)階段平均值形式表示,同列不同小寫字母表示不同處理之間差異顯著(<0.05)。

    2.3 不同水分條件下土壤溶液無機(jī)磷和總磷質(zhì)量濃度的變化

    3種水分管理模式下土柱中土壤溶液無機(jī)磷、總磷質(zhì)量濃度的變化于圖4。由圖4可見,水分管理模式對(duì)土壤溶液無機(jī)磷和總磷質(zhì)量濃度有顯著影響,并與土層和時(shí)間有關(guān)。而3層土壤溶液無機(jī)磷質(zhì)量濃度差異不顯著。各土層土壤溶液總磷、無機(jī)磷質(zhì)量濃度表現(xiàn)為:中層>表層>下層。表層、中層和下層土壤溶液總磷質(zhì)量濃度表現(xiàn)為:模式3>模式2[表層:模式3(1.27 mg/L)>模式2(0.82 mg/L);中層:模式3(1.60 mg/L)>模式2(1.40 mg/L);下層:模式3(0.31 mg/L)>模式2(0.17 mg/L)>模式1(0.13 mg/L)]。

    圖4 3種水分管理模式下土壤溶液無機(jī)磷和總磷動(dòng)態(tài)變化

    培養(yǎng)期間模式1—模式3無機(jī)磷變化均呈“Z”型雙峰變化(圖4)。隨著培養(yǎng)時(shí)間延長(zhǎng)模式2—模式3表層和模式1下層在第21天出現(xiàn)最大值隨后迅速下降,下降至最小值后又緩慢升高。而模式2—模式3中層和下層則在第28 d達(dá)到最大值。隨后同樣迅速下降至最小值后又緩慢上升。

    在培養(yǎng)期間,模式2和模式3表層和中層土壤溶液總磷質(zhì)量濃度呈倒“V”型單峰變化,而下層土壤溶液總磷質(zhì)量濃度呈“Z”型雙峰變化(圖4)。隨培養(yǎng)時(shí)間延長(zhǎng),模式2和模式3表層土壤溶液總磷質(zhì)量逐漸升高,在40~50 d達(dá)到最大值,隨后逐漸下降,中層土壤溶液總磷質(zhì)量濃度也有類似規(guī)律,但是最大值出現(xiàn)在20 d;隨時(shí)間延長(zhǎng),下層土壤溶液總磷質(zhì)量濃度升高,在10 d達(dá)到最大值,隨后下降,在40 d達(dá)到最小值,然后再繼續(xù)升高。

    2.4 土壤溶液總磷與溶液總鐵、總鋁、總鈣和TOC之間關(guān)系

    表3為土壤溶液中總磷質(zhì)量濃度與土壤溶液中其他物質(zhì)(總鐵、總鋁、總Ca和TOC)之間的相關(guān)系數(shù),表4為土壤溶液中總磷質(zhì)量濃度()與土壤溶液總鐵(1)、總鋁(2)、鈣(3)和TOC(4)之間的逐步回歸方程。由表3、表4可知,模式1下層土壤溶液總磷質(zhì)量濃度與下層土壤溶液總鐵、總鈣、TOC質(zhì)量濃度之間呈正相關(guān)(<0.01),而與土壤溶液總鋁質(zhì)量濃度沒有相關(guān)性;模式1下層土壤溶液總磷質(zhì)量濃度與下層土壤溶液總鐵、總鋁、總鈣、TOC質(zhì)量濃度逐步回歸分析,只有TOC包含在方程中。模式2表層土壤溶液總磷質(zhì)量濃度與總鋁、總鈣、TOC正極顯著相關(guān),而與土壤溶液總鐵正顯著相關(guān);表層土壤溶液總磷質(zhì)量濃度與表層土壤溶液總鐵、總鋁、總鈣、TOC質(zhì)量濃度之間逐步回歸方程中僅保留總鋁和TOC。模式2中層土壤溶液總磷質(zhì)量濃度與中層土壤溶液總鐵、總鋁、總鈣質(zhì)量濃度極顯著正相關(guān),與土壤溶液TOC質(zhì)量濃度沒有相關(guān)性;中層土壤溶液總磷質(zhì)量濃度與表層土壤溶液總鐵、總鋁、總鈣、TOC質(zhì)量濃度之間逐步回歸方程中僅保留總鋁。模式2下層土壤溶液總磷質(zhì)量濃度與下層土壤溶液總鐵、總鈣、TOC質(zhì)量濃度極顯著正相關(guān),而與土壤溶液總鋁質(zhì)量濃度負(fù)相關(guān)達(dá)到顯著水平;下層土壤溶液總磷質(zhì)量濃度與下層土壤溶液總鐵、總鋁、總鈣、TOC質(zhì)量濃度之間逐步回歸方程中僅保留總鈣和TOC。

    模式3表層和中層土壤溶液總磷與土壤溶液總鐵、總鋁、總鈣、TOC之間呈正相關(guān)(<0.01)。模式3表層土壤溶液總磷質(zhì)量濃度與表層土壤溶液總鐵、總鋁、總鈣、TOC之間逐步回歸方程中僅保留總鐵;模式3中層土壤溶液總磷質(zhì)量濃度與中層土壤溶液總鐵、總鋁、總鈣、TOC質(zhì)量濃度之間逐步回歸方程中僅保留總鈣。

    模式3下層土壤溶液總磷與下層土壤溶液總鐵、總鈣、TOC正相關(guān),而與總鋁相關(guān)性沒有達(dá)到顯著性水平。模式3下層土壤溶液總磷質(zhì)量濃度與下層土壤溶液總鐵、總鋁、總鈣、TOC之間逐步回歸方程中僅保留總鈣和TOC。

    總體來講,前40 d,土壤溶液總磷質(zhì)量濃度與TOC、總鐵、總鋁和總鈣相關(guān)性要強(qiáng)于整個(gè)時(shí)期。

    表3 土壤溶液中總磷與土壤溶液中總鐵、總鋁、總Ca和TOC的相關(guān)系數(shù)

    注 *和**分別表示相關(guān)性達(dá)到0.05和0.01顯著性水平。

    表4 土壤溶液中總磷質(zhì)量濃度(y)與土壤溶液總鐵(x1)、 總鋁(x2)、鈣(x3)和TOC(x4)之間逐步回歸方程

    3 討論

    在整個(gè)培養(yǎng)期間內(nèi),所有土層土壤溶液總磷質(zhì)量濃度都是隨土壤水分增加而升高(圖4)??赡苡幸韵聨追矫嬖?。水分升高促進(jìn)土壤有機(jī)磷降解、釋放:隨水分升高,土壤微生物活動(dòng)增強(qiáng)[21],進(jìn)而加速土壤有機(jī)磷分解;TOC可以與磷酸根競(jìng)爭(zhēng)吸附土壤上的吸附位點(diǎn),促進(jìn)可溶性磷釋放[22]。TOC也可以與磷產(chǎn)生絡(luò)合增溶溶解現(xiàn)象,進(jìn)一步促進(jìn)土壤磷的釋放[23-24];隨土壤水分升高,土壤溶液TOC質(zhì)量濃度也增加(圖2),土壤溶液TOC與土壤溶液總磷相關(guān)性顯著(表3),3種水分管理模式下層土壤溶液總磷與總鐵、總鋁、總鈣和TOC逐步方程中都含TOC(表4),這些證實(shí)水分升高促進(jìn)土壤有機(jī)質(zhì)分解,間接證實(shí)水分升高促進(jìn)土壤有機(jī)磷分解;水分升高驅(qū)動(dòng)了土壤無機(jī)磷釋放。隨著土壤水分升高,導(dǎo)致土壤還原性增強(qiáng),土壤Fe3+被還原為Fe2+[25],從而導(dǎo)致土壤中Fe-P釋放,也破壞了土壤Al-P和Ca-P結(jié)構(gòu)。隨土壤水分升高,土壤溶液趨于下降(圖2)。3種水分管理模式中土壤溶液總磷與總鐵、總鋁、總鈣和TOC逐步方程中含總鐵、總鋁、總鈣因子(表4),指示水分驅(qū)動(dòng)土壤無機(jī)磷組分釋放。土壤TOC維持土壤溶液鐵穩(wěn)定性[25]。土壤可溶有機(jī)質(zhì)與土壤釋放到溶液鐵磷、鋁磷形成絡(luò)合物維持了土壤溶液高質(zhì)量濃度磷存在。研究證實(shí),低分子有機(jī)物(TOC)可以與鐵磷或者鋁磷形成三元復(fù)合體[22, 24, 26-27]。水分誘導(dǎo)土壤pH值上升,促進(jìn)土壤磷溶解釋放。隨著土壤水分升高,土壤溶液pH值從4.45上升到6.5左右,這與蔡振國(guó)等[25]試驗(yàn)結(jié)果相同,淹水處理高磷土壤pH值從4.85增加到6.5,pH值上升進(jìn)而促進(jìn)酸性土壤磷溶解釋放。酸性土壤中磷有效性一般較低,pH值升高到6.5~6.8,土壤磷固定作用下降,有利于磷釋放,提高土壤磷有效性[28];土壤淹水引起土壤死亡微生物分解導(dǎo)致的磷釋放。土壤淹水滲透壓上升,會(huì)導(dǎo)致死亡微生物細(xì)胞內(nèi)的磷釋放[29]。王子琬等[29]Blackwell等[30]證實(shí)干濕交替過程中釋放的總?cè)芙饬字辽?8%來自裂解的微生物細(xì)胞。

    研究結(jié)果發(fā)現(xiàn),表層、中層、下層土壤溶液磷質(zhì)量濃度均值分別達(dá)到0.82~1.27、1.40~1.60 mg/L和0.13~0.31 mg/L(圖4),大于地表水富營(yíng)養(yǎng)化標(biāo)準(zhǔn)值(0.1 mg/L)[14]。因此,退耕還濕土壤磷釋放進(jìn)入土壤溶液,隨后通過向上覆水?dāng)U散,提高湖泊水體磷濃度,其影響湖泊水體磷質(zhì)量濃度最關(guān)鍵時(shí)期在退耕還濕40 d前后,因?yàn)榇藭r(shí)土壤磷釋放導(dǎo)致土壤水磷質(zhì)量濃度最大(圖4)。因此,退耕還濕土壤磷釋放引起的磷釋放造成生態(tài)環(huán)境問題值得關(guān)注。隨水分升高,土壤溶液磷質(zhì)量濃度增加(圖4)。因此,退耕還濕過程中,采取逐步灌水措施有利于降低灌水對(duì)水體磷沖擊。這些研究有待野外觀察驗(yàn)。

    4 結(jié)論

    土壤水分顯著影響土壤溶液磷質(zhì)量濃度,且隨著水分增加而增加。3種水分管理模式下土壤溶液總磷質(zhì)量濃度都大于地表水富營(yíng)養(yǎng)化標(biāo)準(zhǔn)值(0.1 mg/L)[14],尤其是模式3土壤淹水飽和條件下。因此,退耕還濕期間土壤磷釋放,特別是土壤淹水飽和情況下,可能會(huì)引起上覆水磷質(zhì)量濃度突然升高,其造成的生態(tài)環(huán)境問題值得關(guān)注。

    [1] 李紅艷, 王維峰. 淺淡三江平原退耕還濕工程的深遠(yuǎn)意義[J]. 現(xiàn)代化農(nóng)業(yè), 2005(2): 1-3.

    LI Hongyan, WANG Weifeng. Talking about profound significance of returning arable land back to wetland in Sanjiang plain[J]. Modernizing Agriculture, 2005(2): 1-3.

    [2] 張文琦, 宋戈. 三江平原濕地農(nóng)田化進(jìn)程中土地利用變化及影響因素分析[J]. 北方園藝, 2019, 442(19): 150-158.

    ZHANG Wenqi, SONG Ge. Analysis of land use change and its influencing factors in the process of Wetland Farmland in Sanjiang Plain [J]. Northern Horticulture, 2019, 442(19): 150-158.

    [3] 王延吉, 神祥金, 呂憲國(guó). 1980—2015年?yáng)|北沼澤濕地景觀格局及氣候變化特征[J]. 地球與環(huán)境, 2020, 48(3): 348-357.

    WANG Yanji, SHEN Xiangjin, LYU Xianguo. Change characteristics of landscape pattern and climate in marsh areas of northeast China during 1980-2015[J]. Earth and Environment, 2020, 48(3): 348-357.

    [4] 楊春霞, 鄭華, 歐陽(yáng)志云. 三江平原土地利用變化、效應(yīng)與驅(qū)動(dòng)力[J]. 環(huán)境保護(hù)科學(xué), 2020, 46(5): 99-104.

    YANG Chunxia, ZHENG Hua, OUYANG Zhiyun. Changes, effect and driving forces of land use in Sanjiang plain[J]. Environmental Protection Science, 2020, 46(5): 99-104.

    [5] 何興元, 賈明明, 王宗明, 等. 基于遙感的三江平原濕地保護(hù)工程成效初步評(píng)估[J]. 中國(guó)科學(xué)院院刊, 2017, 32(1): 3-10.

    HE Xingyuan, JIA Mingming, WANG Zongming, et al. Remote sensing based evaluation of effectiveness of wetland protection project in Sanjiang plain, China[J]. Bulletin of Chinese Academy of Sciences, 2017, 32(1): 3-10.

    [6] YANG J S, LIU J S, HU X J, et al. Changes of soil organic carbon, nitrogen and phosphorus concentrations under different land uses in marshes of Sanjiang Plain[J]. Acta Ecologica Sinica, 2013, 33(6): 332-337.

    [7] 高居娟. 水分對(duì)若爾蓋濕地不同微生境土壤溫室氣體排放的影響[D]. 北京: 北京林業(yè)大學(xué), 2016..

    GAO Jujuan. Effects of moisture on greenhouse gas emission of different microhabitats soil in zoige wetland[D]. Beijing: Beijing Forestry University, 2016.

    [8] 童夢(mèng)瑩, 黃家權(quán), 李長(zhǎng)江. 淹水脅迫對(duì)櫻桃番茄苗期形態(tài)特征及葉綠素?zé)晒馓匦缘挠绊?[J]. 灌溉排水學(xué)報(bào), 2019, 38(11): 8-15.

    TONG Mengying, HUANG Jiaquan, LI Changjiang. The effect of flooding stress on the morphological characteristics and chlorophyll fluorescence characteristics of cherry tomato seedling [J]. Journal of Irrigation and Drainage, 2019, 38(11): 8-15.

    [9] 晏軍, 王偉義, 費(fèi)月躍,等. 鹽土條件下不同品系藜麥生理特性和產(chǎn)量指標(biāo)對(duì)花期漬水脅迫的響應(yīng) [J]. 灌溉排水學(xué)報(bào), 2021, 40(2): 47-54.

    YAN Jun, WANG Weiyi, FEI Yueyue, et al. Response of physiological characteristics and yield indexes of quinoa in different lines to waterlogging stress in flowering period [J]. Journal of Irrigation and Drainage, 2021, 40 (2): 47-54.

    [10] 周青云, 李夢(mèng)初, 漆棟良, 等. 拔節(jié)期淹水條件下施氮量對(duì)春玉米生理特性的影響 [J]. 灌溉排水學(xué)報(bào), 2020, 39(S2): 40-44.

    ZHOU Qingyun, LI Mengchu, QI Dongliang, et al. Effects of nitrogen application rate on physiological characteristics of spring maize under waterlogging at jointing stage [J]. Journal of Irrigation and Drainage, 2020, 39 (S2): 40-44.

    [11] 張新燕, 王浩翔, 牛文全. 不同水氮供應(yīng)模式對(duì)設(shè)施番茄生長(zhǎng)及產(chǎn)量的影響 [J]. 灌溉排水學(xué)報(bào), 2020, 39(11): 55-60.

    ZHANG Xinyan, WANG Haoxiang, NIU Wenquan. The effect of different water and nitrogen supply modes on the growth and yield of tomato facilities [J]. Journal of Irrigation and Drainage, 2020, 39 (11): 55-60.

    [12] 賈詠霖, 屈忠義, 丁艷宏, 等. 不同灌溉方式下施用生物炭對(duì)土壤水鹽運(yùn)移規(guī)律及玉米水分利用效率的影響 [J]. 灌溉排水學(xué)報(bào), 2020, 39(8): 44-51.

    JIA Yonglin, QU Zhongyi, DING Yanhong, et al. Effects of biochar application on soil water and salt movement and water use efficiency of Maize under different irrigation methods [J]. Journal of Irrigation and Drainage, 2020, 39(8): 44-51.

    [13] 田娟, 劉凌, 丁海山, 等. 淹水土壤土-水界面磷素遷移轉(zhuǎn)化研究[J]. 環(huán)境科學(xué), 2008, 29(7): 1 818-1 823.

    TIAN Juan, LIU Ling, DING Haishan, et al. Mobilization and transformation of phosphorus from water-soil interface of flooded soil[J]. Environmental Science, 2008, 29(7): 1 818-1 823.

    [14] 濕地生態(tài)風(fēng)險(xiǎn)評(píng)估技術(shù)規(guī)范. GB/T 27647—2011[S].

    Technical specification for wetland ecological risk assessment. GB/T 27647—2011[S].

    [15] 崔虎, 王莉霞, 歐洋, 等. 濕地生態(tài)系統(tǒng)磷遷移轉(zhuǎn)化機(jī)制研究進(jìn)展[J]. 水生態(tài)學(xué)雜志, 2020, 41(2): 105-112.

    CUI Hu, WANG Lixia, OU Yang, et al. Research progress on phosphorus migration and transformation in wetland ecosystem[J]. Journal of Hydroecology, 2020, 41(2): 105-112.

    [16] 周健, 李春輝, 張志永, 等. 淹水落干下三峽水庫(kù)消落帶土壤無機(jī)磷形態(tài)轉(zhuǎn)化特征[J]. 環(huán)境科學(xué), 2018, 39(1): 130-136.

    ZHOU Jian, LI Chunhui, ZHANG Zhiyong, et al. Effects of flooding and drying on the transformation of soil inorganic phosphorus in the water-level-fluctuating zone of the Three Gorges reservoir, China[J]. Environmental Science, 2018, 39(1): 130-136.

    [17] LAING G, MEERS E, DEWISPELAERE M, et al. Effect of water table level on metal mobility at different depths in wetland soils of the Scheldt estuary (Belgium)[J]. Water, Air, and Soil Pollution, 2009, 202(1/2/3/4): 353-367.

    [18] CLARK J M, HEINEMEYER A, MARTIN P, et al. Processes controlling DOC in pore water during simulated drought cycles in six different UK peats[J]. Biogeochemistry, 2012, 109(1/2/3): 253-270.

    [19] 魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京: 中國(guó)農(nóng)業(yè)科技出版社, 2000.

    LU Rukun. Analytical methods of soil agrochemistry[M]. Beijing: China Agricultural Science and Technology Press, 2000.

    [20] 魏復(fù)盛. 水和廢水監(jiān)測(cè)分析方法[M]. 第4版. 北京:中國(guó)環(huán)境科學(xué)出版社, 2002.

    WEI Fusheng. Methods for monitoring and analyzing water and wastewater [M]. 4th Edition. Beijing: China Environmental Science Press, 2002.

    [21] 陳昊. 水分狀況對(duì)內(nèi)蒙古典型草原土壤微生物量及群落結(jié)構(gòu)的影響[D]. 北京: 中國(guó)農(nóng)業(yè)大學(xué), 2018.

    CHEN Hao. Impact of water status on soil microbial biomass and community structure in typical grasslands of Inner Mongolia[D]. Beijing: China Agricultural University, 2018.

    [22] HUTCHISON K J, HESTERBERG D. Dissolution of phosphate in a phosphorus-enriched ultisol as affected by microbial reduction[J]. Journal of Environmental Quality, 2004, 33(5): 1 793-1 802.

    [23] JOHNSON S E, LOEPPERT R H. Role of organic acids in phosphate mobilization from iron oxide[J]. Soil Science Society of America Journal, 2006, 70(1): 222-234.

    [24] WILLETT I R. Causes and prediction of changes in extractable phosphorus during flooding[J]. Soil Research, 1989, 27(1): 45.

    [25] 蔡振國(guó). 淹水——落干條件下紅壤中磷與鐵的形態(tài)轉(zhuǎn)化及其耦合機(jī)制[D]. 武漢: 華中農(nóng)業(yè)大學(xué), 2019.

    CAI Zhenguo. Speciation transformation and coupling mechanism of phosphorus and iron in red soil under flooding-drying condition[D]. Wuhan: Huazhong Agricultural University, 2019.

    [26] 田娟, 劉凌, 董貴明, 等. 淹水土壤磷釋放機(jī)理研究進(jìn)展[J]. 土壤通報(bào), 2008, 39(2): 426-430.

    TIAN Juan, LIU Ling, DONG Guiming, et al. Study progress of phosphorus release mechanics in flooded soils[J]. Chinese Journal of Soil Science, 2008, 39(2): 426-430.

    [27] 陳凌玉, 劉飛, 趙雙嬌, 等. 干濕交替條件下鐵氧化對(duì)水稻土CO2排放的影響[J]. 河南農(nóng)業(yè)科學(xué), 2017, 46(8): 72-76.

    CHEN Lingyu, LIU Fei, ZHAO Shuangjiao, et al. Effects of iron oxidation on CO2emission in paddy soils under dry-wet alternation conditions[J]. Journal of Henan Agricultural Sciences, 2017, 46(8): 72-76.

    [28] 黃昌勇. 土壤學(xué)[M]. 北京: 中國(guó)農(nóng)業(yè)出版社, 2000..

    HUANG Changyong. Soil science[M]. Beijing: Chinese Agriculture Press, 2000.

    [29] 王子琬, 梁新強(qiáng). 土壤干濕交替對(duì)磷素釋放的影響機(jī)制[J]. 環(huán)境生態(tài)學(xué), 2020, 2(5): 54-58.

    WANG Ziwan, LIANG Xinqiang. Effects of alternate drying- rewetting process on migration and transformation of soil phosphorus[J]. Environmental Ecology, 2020, 2(5): 54-58.

    [30] BLACKWELL M S A, BROOKES P C, DE LA FUENTE-MARTINEZ N, et al. Phosphorus solubilization and potential transfer to surface waters from the soil microbial biomass following drying-rewetting and freezing-thawing[J]. Advances in Agronomy, 2010, 106: 1-35.

    The Effects of Soil Water Content on Distribution of Mobile Phosphorus Concentration in Soil

    WU Rui1, HU Zhengyi1,3, SUN Yiqi1,GENG Ziqi3,CUI Yanshang1,2*

    (1. School of Resources and Environment, University of Chinese Academy of Sciences, Beijing 101408, China;2. Research Center for Eco-environmental Sciences , CAS, Beijing 100085, China;3. Sino-danish College, University of Chinese Academy of Sciences, Beijing 100190, China)

    】Concentration of mobile phosphorus (P) in soil affects not only plant growth but also water quality because of its potential leaching. Most studies on soil P focused on topsoil while overlooking subsoil which is equally important.The objective of this paper is to study the dynamics of mobile P at different soil depths and its driving forces in attempts to help improve P uptake by plants while in the meantime reducing its leaching.【】Column experiments filled with soil taken from Sanjiang plain in China were conducted with average volumetric soil water content controlled at 50% (M1), 75% (M2) and 100% (M3) respectively. In each treatment, we measured the changes in mobile phosphorus in the depths of 0~10 cm (topsoil) 10~20 cm (middle soil) and 20~30 cm (subsoil). We also analyzed its relationship with total organic carbon (TOC), total Fe, total Al and total Ca in soil solution; the experiments lasted 84 days.【】When soil water content was 50%, soil solution was not extractable from the topsoil and the middle soil, and only was limited water extracted from the subsoil. When soil water content exceeded 75%, soil solution was more mobile for extraction in which the P concentration in the middle soil was found to be the highest, 1.3~1.7 times that in the topsoil and 5.2~8.2 times that in the subsoil. The impact of soil water on P concentration varied with soil depth and was ranked in the following order: In the topsoil, M3 (1.27 mg/L) > M2 (0.82 mg/L), in the middle soil M3 (1.60 mg/L) > M2 (1.40 mg/L), and in the subsoil M3 (0.31 mg/L) > M2 (0.17 mg/L) >M1 (0.13 mg/L). A stepwise regression analysis indicated that P concentration in the soil solution was related to TOC, total Fe, total Al, and total Ca in soil solution, with the correlations varying with soil water and soil depth. These suggested that improving soil water management could shift dissolution and/or decomposition of P in different forms in the soil, and can thus bed used to manage soil P. Regardless of soil water content, P concentration in the soil was much higher than the critical P concentration, 0.1 mg/L, which could trigger eutrophication in surface water.【】Increasing soil water could result in a quick release of P from soil to the overlying surface water, triggering eutrophication as a result. Care thus must be taken when temporally returning dry farmland to wetland.

    water management; soil solution; phosphorus; Sanjiang plain; returning farmland to wetland

    X53

    A

    10.13522/j.cnki.ggps.2021065

    1672 – 3317(2021)08 - 0073 - 08

    吳瑞, 胡正義, 孫藝齊, 等.室內(nèi)模擬水分管理對(duì)土壤溶液磷質(zhì)量濃度影響研究[J]. 灌溉排水學(xué)報(bào), 2021, 40(8): 73-80.

    WU Rui, HU Zhengyi, SUN Yiqi, et al. The Effects of Soil Water Content on Distribution of Mobile Phosphorus Concentration in Soil [J]. Journal of Irrigation and Drainage, 2021, 40(8): 73-80.

    2021-02-13

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFC0500408)

    吳瑞(1994-),男,遼寧沈陽(yáng)人。碩士研究生,主要研究土壤環(huán)境化學(xué)。E-mail: wurui18@mails.ucas.ac.cn

    崔巖山(1972-),男。教授,博士生導(dǎo)師,研究方向土壤污染物環(huán)境行為及其修復(fù)技術(shù)。E-mail: cuiyanshan@.ucas.ac.cn

    責(zé)任編輯:趙宇龍

    猜你喜歡
    土壤溶液總鐵下層
    洋沙泡水庫(kù)底泥總鐵遷移影響因素實(shí)驗(yàn)研究
    氮添加對(duì)亞熱帶常綠闊葉林土壤溶液化學(xué)特性的影響
    互花米草入侵對(duì)膠州灣河口濕地土壤總鐵分布的影響
    生活飲用水相關(guān)指標(biāo)與細(xì)菌總數(shù)的相關(guān)性研究
    臨桂區(qū)土壤溶液的水化學(xué)特征及其溶蝕能力淺析
    一類多個(gè)下層的雙層規(guī)劃問題
    積雪
    陜西橫山羅圪臺(tái)村元代壁畫墓發(fā)掘簡(jiǎn)報(bào)
    考古與文物(2016年5期)2016-12-21 06:28:48
    土壤增溫對(duì)杉木幼林不同深度土壤溶液N03-濃度的影響
    土壤增溫對(duì)杉木幼林不同深度土壤溶液DOM的影響
    午夜久久久在线观看| 亚洲少妇的诱惑av| 亚洲精品美女久久久久99蜜臀 | 午夜精品国产一区二区电影| 热re99久久精品国产66热6| 精品第一国产精品| 99re6热这里在线精品视频| 亚洲久久久国产精品| 久热这里只有精品99| 亚洲国产毛片av蜜桃av| 午夜福利一区二区在线看| 亚洲精品美女久久av网站| 一区福利在线观看| 亚洲第一av免费看| 女性生殖器流出的白浆| 欧美日韩综合久久久久久| 老鸭窝网址在线观看| 两性夫妻黄色片| 日本av免费视频播放| 爱豆传媒免费全集在线观看| 爱豆传媒免费全集在线观看| 国产日韩一区二区三区精品不卡| 中文乱码字字幕精品一区二区三区| a 毛片基地| 婷婷色av中文字幕| 国产亚洲精品第一综合不卡| 国产淫语在线视频| 十八禁网站网址无遮挡| 51午夜福利影视在线观看| 国产精品久久久久久人妻精品电影 | 天天躁夜夜躁狠狠久久av| 国产亚洲欧美在线一区二区| 成人亚洲欧美一区二区av| 天天操日日干夜夜撸| 欧美精品一区二区免费开放| 少妇人妻久久综合中文| 国产视频首页在线观看| 超碰成人久久| 啦啦啦中文免费视频观看日本| 午夜91福利影院| 亚洲自偷自拍图片 自拍| 国产精品久久久久久精品电影小说| 国产精品久久久久成人av| 欧美精品亚洲一区二区| 欧美精品啪啪一区二区三区 | 久久人人爽av亚洲精品天堂| av电影中文网址| 少妇 在线观看| 男人舔女人的私密视频| 热99久久久久精品小说推荐| 亚洲成色77777| av网站免费在线观看视频| 亚洲色图 男人天堂 中文字幕| 少妇的丰满在线观看| 国产精品一区二区在线不卡| 国产免费视频播放在线视频| 精品少妇一区二区三区视频日本电影| 日韩电影二区| 免费高清在线观看日韩| 国产精品国产三级专区第一集| 操美女的视频在线观看| 免费在线观看日本一区| 国产片特级美女逼逼视频| 女人被躁到高潮嗷嗷叫费观| 久久狼人影院| 午夜福利在线免费观看网站| 两性夫妻黄色片| 一区在线观看完整版| 婷婷丁香在线五月| 亚洲国产最新在线播放| 蜜桃国产av成人99| 亚洲国产看品久久| 午夜福利免费观看在线| 久久精品亚洲av国产电影网| 美女扒开内裤让男人捅视频| 天天躁日日躁夜夜躁夜夜| 国产视频一区二区在线看| av天堂在线播放| 国产av一区二区精品久久| 女人久久www免费人成看片| 满18在线观看网站| 日韩一卡2卡3卡4卡2021年| 男女高潮啪啪啪动态图| 人人澡人人妻人| 99国产精品99久久久久| 国产高清视频在线播放一区 | 尾随美女入室| 久久精品国产亚洲av高清一级| 又大又爽又粗| 亚洲成国产人片在线观看| 色视频在线一区二区三区| 无限看片的www在线观看| 激情视频va一区二区三区| 丁香六月天网| 日韩熟女老妇一区二区性免费视频| 大型av网站在线播放| 国产伦理片在线播放av一区| 女性生殖器流出的白浆| 成人手机av| 日韩熟女老妇一区二区性免费视频| 久久精品国产亚洲av高清一级| 91国产中文字幕| 少妇人妻久久综合中文| 国产成人啪精品午夜网站| www日本在线高清视频| 久久 成人 亚洲| 丰满少妇做爰视频| 欧美日韩黄片免| 国产野战对白在线观看| 国产亚洲欧美在线一区二区| 大香蕉久久网| e午夜精品久久久久久久| 可以免费在线观看a视频的电影网站| 欧美日韩视频高清一区二区三区二| 亚洲欧美成人综合另类久久久| 久久人妻福利社区极品人妻图片 | 国产伦人伦偷精品视频| 久久性视频一级片| 91精品国产国语对白视频| 亚洲av综合色区一区| 国产成人一区二区在线| 嫁个100分男人电影在线观看 | 看十八女毛片水多多多| 国产色视频综合| 国产成人啪精品午夜网站| 九色亚洲精品在线播放| 久久亚洲国产成人精品v| 国产一区二区三区av在线| 精品第一国产精品| 亚洲伊人色综图| 欧美精品高潮呻吟av久久| 午夜福利在线免费观看网站| 操出白浆在线播放| 男女国产视频网站| 美女国产高潮福利片在线看| 男人添女人高潮全过程视频| 1024视频免费在线观看| 韩国精品一区二区三区| 国产免费福利视频在线观看| 天堂8中文在线网| 亚洲国产欧美日韩在线播放| 国产精品av久久久久免费| 各种免费的搞黄视频| 亚洲国产精品成人久久小说| 国产一卡二卡三卡精品| 菩萨蛮人人尽说江南好唐韦庄| 国语对白做爰xxxⅹ性视频网站| 精品人妻1区二区| 日本a在线网址| 亚洲国产中文字幕在线视频| 久久久久久久国产电影| 国产成人精品久久久久久| 搡老乐熟女国产| 婷婷色av中文字幕| 在线看a的网站| 国产视频首页在线观看| 久久久国产欧美日韩av| 制服诱惑二区| 咕卡用的链子| 热99久久久久精品小说推荐| 99热网站在线观看| 久久久久国产精品人妻一区二区| 黄片小视频在线播放| 中文字幕人妻熟女乱码| 亚洲一区中文字幕在线| 女性生殖器流出的白浆| av片东京热男人的天堂| 韩国高清视频一区二区三区| 中文字幕人妻丝袜制服| 中文精品一卡2卡3卡4更新| 久久国产精品影院| 黄片小视频在线播放| 人人妻人人澡人人看| 国产男人的电影天堂91| 亚洲三区欧美一区| 中文字幕另类日韩欧美亚洲嫩草| www.精华液| 日日摸夜夜添夜夜爱| 三上悠亚av全集在线观看| 久久狼人影院| 精品久久久久久电影网| 国语对白做爰xxxⅹ性视频网站| 99久久99久久久精品蜜桃| 国产高清不卡午夜福利| 午夜老司机福利片| 亚洲欧美日韩高清在线视频 | 日本午夜av视频| 亚洲国产精品一区三区| 午夜福利视频在线观看免费| 亚洲av日韩在线播放| 七月丁香在线播放| 妹子高潮喷水视频| 亚洲久久久国产精品| 首页视频小说图片口味搜索 | 欧美老熟妇乱子伦牲交| 国产成人影院久久av| av又黄又爽大尺度在线免费看| 亚洲精品第二区| 国产精品久久久人人做人人爽| av福利片在线| 美女福利国产在线| 18禁国产床啪视频网站| 精品国产一区二区久久| 国产精品熟女久久久久浪| 国产成人一区二区在线| 男的添女的下面高潮视频| 亚洲国产欧美一区二区综合| 国产精品一区二区在线不卡| 亚洲精品国产av蜜桃| 久久久久久久久久久久大奶| 国产精品一区二区免费欧美 | 飞空精品影院首页| 国产精品九九99| 2018国产大陆天天弄谢| 精品少妇一区二区三区视频日本电影| 一级毛片黄色毛片免费观看视频| 亚洲欧洲精品一区二区精品久久久| 2021少妇久久久久久久久久久| 国产成人系列免费观看| 大陆偷拍与自拍| 男女高潮啪啪啪动态图| 欧美成人精品欧美一级黄| 天堂中文最新版在线下载| 亚洲综合色网址| 久久久久久亚洲精品国产蜜桃av| 欧美日韩精品网址| 青春草视频在线免费观看| 在线观看免费午夜福利视频| 亚洲国产欧美网| 男女无遮挡免费网站观看| 老鸭窝网址在线观看| 男女国产视频网站| 亚洲国产毛片av蜜桃av| 亚洲黑人精品在线| 国产欧美日韩综合在线一区二区| 欧美久久黑人一区二区| 欧美精品高潮呻吟av久久| 日本午夜av视频| 亚洲中文日韩欧美视频| 午夜免费男女啪啪视频观看| 99热国产这里只有精品6| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 99国产精品免费福利视频| 国产精品人妻久久久影院| 在线观看免费高清a一片| 韩国高清视频一区二区三区| 高清黄色对白视频在线免费看| 欧美中文综合在线视频| 免费在线观看视频国产中文字幕亚洲 | 国产高清视频在线播放一区 | 人妻一区二区av| 久久久久视频综合| 99九九在线精品视频| 悠悠久久av| 99国产精品免费福利视频| 精品国产乱码久久久久久小说| 亚洲欧美中文字幕日韩二区| 久久av网站| 大香蕉久久成人网| 午夜福利视频精品| 在线观看一区二区三区激情| 精品福利永久在线观看| 脱女人内裤的视频| 国产欧美日韩一区二区三 | 色网站视频免费| 十八禁人妻一区二区| 亚洲成人手机| av欧美777| 极品人妻少妇av视频| 高清欧美精品videossex| 在线精品无人区一区二区三| 18禁裸乳无遮挡动漫免费视频| 老司机靠b影院| 精品国产乱码久久久久久男人| 国产成人91sexporn| 另类精品久久| 亚洲成色77777| 精品久久久久久电影网| 亚洲欧美激情在线| 五月开心婷婷网| 亚洲欧洲精品一区二区精品久久久| 免费看av在线观看网站| 国产福利在线免费观看视频| 国产精品国产av在线观看| 亚洲人成电影观看| 欧美中文综合在线视频| 国产精品香港三级国产av潘金莲 | 精品一区二区三区av网在线观看 | 色94色欧美一区二区| 欧美xxⅹ黑人| 青青草视频在线视频观看| 伊人久久大香线蕉亚洲五| 男女下面插进去视频免费观看| 激情五月婷婷亚洲| 国产在线观看jvid| 国产色视频综合| 美女视频免费永久观看网站| 久久人妻福利社区极品人妻图片 | 久久精品国产亚洲av涩爱| 欧美中文综合在线视频| 日韩视频在线欧美| bbb黄色大片| 夜夜骑夜夜射夜夜干| 人人妻人人澡人人爽人人夜夜| 在线av久久热| 亚洲精品一卡2卡三卡4卡5卡 | 十八禁高潮呻吟视频| 精品免费久久久久久久清纯 | 亚洲成人手机| 91精品三级在线观看| 狂野欧美激情性xxxx| 我要看黄色一级片免费的| 人人妻人人澡人人爽人人夜夜| 大片电影免费在线观看免费| 国产一级毛片在线| 欧美日韩成人在线一区二区| 咕卡用的链子| 亚洲,欧美精品.| 最新在线观看一区二区三区 | videos熟女内射| 最近中文字幕2019免费版| 黄色毛片三级朝国网站| h视频一区二区三区| 精品国产乱码久久久久久男人| 热re99久久国产66热| 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 国产成人免费无遮挡视频| 欧美日韩亚洲国产一区二区在线观看 | 久热爱精品视频在线9| 黄色一级大片看看| 亚洲 欧美一区二区三区| 在线观看一区二区三区激情| xxx大片免费视频| 免费看av在线观看网站| 色网站视频免费| 亚洲国产精品一区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲国产欧美在线一区| 免费少妇av软件| 亚洲天堂av无毛| 天堂俺去俺来也www色官网| 精品免费久久久久久久清纯 | 国产野战对白在线观看| 麻豆av在线久日| 精品卡一卡二卡四卡免费| 欧美变态另类bdsm刘玥| 亚洲国产av影院在线观看| 最近手机中文字幕大全| 捣出白浆h1v1| 精品少妇久久久久久888优播| 日本91视频免费播放| 免费在线观看黄色视频的| 亚洲欧美精品综合一区二区三区| 高清视频免费观看一区二区| 我的亚洲天堂| 久久亚洲精品不卡| 最黄视频免费看| 美女视频免费永久观看网站| 免费看十八禁软件| 在现免费观看毛片| 国产欧美亚洲国产| 欧美日韩一级在线毛片| 亚洲精品第二区| 狠狠婷婷综合久久久久久88av| 久久午夜综合久久蜜桃| 热99国产精品久久久久久7| 你懂的网址亚洲精品在线观看| 久久国产精品人妻蜜桃| 国产免费又黄又爽又色| 亚洲av欧美aⅴ国产| 亚洲中文日韩欧美视频| 国产不卡av网站在线观看| 色综合欧美亚洲国产小说| 妹子高潮喷水视频| 日本午夜av视频| 老司机在亚洲福利影院| 国产精品一区二区在线观看99| 王馨瑶露胸无遮挡在线观看| 国产精品国产av在线观看| 亚洲av男天堂| 亚洲三区欧美一区| a级毛片在线看网站| 一区二区三区四区激情视频| 男女床上黄色一级片免费看| av网站在线播放免费| 亚洲专区中文字幕在线| 免费女性裸体啪啪无遮挡网站| 99国产精品一区二区三区| 啦啦啦 在线观看视频| 午夜老司机福利片| 国产日韩欧美亚洲二区| 波多野结衣一区麻豆| 国产人伦9x9x在线观看| 久久av网站| av又黄又爽大尺度在线免费看| 欧美97在线视频| 国产成人一区二区三区免费视频网站 | 久久 成人 亚洲| 一区二区日韩欧美中文字幕| 黑人猛操日本美女一级片| 女人高潮潮喷娇喘18禁视频| 国产一区二区三区av在线| 午夜视频精品福利| 91国产中文字幕| 精品国产国语对白av| 中文字幕另类日韩欧美亚洲嫩草| 午夜福利视频精品| 熟女av电影| 精品一区在线观看国产| 欧美在线黄色| 亚洲国产欧美在线一区| 一个人免费看片子| 亚洲国产av新网站| 十八禁高潮呻吟视频| av福利片在线| 亚洲少妇的诱惑av| 国产男女内射视频| 极品少妇高潮喷水抽搐| 人人妻人人澡人人看| 性色av一级| 亚洲av日韩在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 黄网站色视频无遮挡免费观看| 欧美精品人与动牲交sv欧美| 国产精品久久久久成人av| 精品人妻熟女毛片av久久网站| 在线观看免费高清a一片| 午夜福利,免费看| 中文精品一卡2卡3卡4更新| 精品一品国产午夜福利视频| 午夜激情久久久久久久| 成人国产一区最新在线观看 | 丝袜在线中文字幕| 国产精品 国内视频| 亚洲自偷自拍图片 自拍| 国产福利在线免费观看视频| 精品第一国产精品| 久久久精品94久久精品| 91精品伊人久久大香线蕉| 脱女人内裤的视频| 亚洲国产欧美日韩在线播放| 欧美精品亚洲一区二区| 搡老岳熟女国产| 亚洲国产最新在线播放| 伊人久久大香线蕉亚洲五| 男女之事视频高清在线观看 | 精品少妇一区二区三区视频日本电影| 人妻一区二区av| 女人精品久久久久毛片| 亚洲欧美中文字幕日韩二区| 亚洲精品国产av蜜桃| 亚洲九九香蕉| 国产成人精品久久久久久| 首页视频小说图片口味搜索 | 亚洲第一av免费看| 亚洲五月婷婷丁香| 午夜激情久久久久久久| 亚洲国产精品国产精品| 亚洲精品日韩在线中文字幕| 亚洲精品国产区一区二| 欧美在线一区亚洲| 丰满少妇做爰视频| 国产不卡av网站在线观看| 少妇人妻久久综合中文| 999精品在线视频| 久久精品成人免费网站| 超碰成人久久| 1024香蕉在线观看| av国产精品久久久久影院| 日韩欧美一区视频在线观看| 免费在线观看黄色视频的| 亚洲三区欧美一区| 嫁个100分男人电影在线观看 | 香蕉丝袜av| 国产av国产精品国产| 新久久久久国产一级毛片| 午夜av观看不卡| 欧美日韩av久久| 一级毛片黄色毛片免费观看视频| 日本欧美国产在线视频| 精品免费久久久久久久清纯 | 国产欧美日韩一区二区三区在线| 丝袜人妻中文字幕| 日韩 亚洲 欧美在线| 欧美日韩综合久久久久久| 亚洲精品日本国产第一区| 99国产综合亚洲精品| 青青草视频在线视频观看| 老司机影院毛片| 久久99热这里只频精品6学生| 亚洲伊人色综图| 最黄视频免费看| 久久人人爽人人片av| 日韩一本色道免费dvd| 五月天丁香电影| 高潮久久久久久久久久久不卡| 丁香六月天网| √禁漫天堂资源中文www| 国产精品二区激情视频| 国产精品三级大全| 自线自在国产av| 亚洲黑人精品在线| 免费在线观看视频国产中文字幕亚洲 | 色精品久久人妻99蜜桃| 日韩免费高清中文字幕av| 午夜福利免费观看在线| 久久精品久久精品一区二区三区| 精品亚洲乱码少妇综合久久| 国产男女内射视频| 精品一区二区三卡| 成人三级做爰电影| 午夜福利视频在线观看免费| 欧美少妇被猛烈插入视频| videos熟女内射| 99久久精品国产亚洲精品| 韩国精品一区二区三区| 国产极品粉嫩免费观看在线| 久久久精品区二区三区| 国产黄频视频在线观看| 精品国产超薄肉色丝袜足j| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 亚洲中文av在线| 制服人妻中文乱码| 波野结衣二区三区在线| 亚洲欧美日韩高清在线视频 | 高清视频免费观看一区二区| 国产免费一区二区三区四区乱码| 国产高清不卡午夜福利| 欧美中文综合在线视频| 亚洲国产日韩一区二区| 一区二区三区四区激情视频| 成年人免费黄色播放视频| 欧美 亚洲 国产 日韩一| 免费观看人在逋| 成年人午夜在线观看视频| 久久99一区二区三区| 日韩免费高清中文字幕av| 亚洲精品久久久久久婷婷小说| 妹子高潮喷水视频| 欧美在线一区亚洲| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 成人18禁高潮啪啪吃奶动态图| 欧美变态另类bdsm刘玥| 欧美精品av麻豆av| 黄片小视频在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 巨乳人妻的诱惑在线观看| 青草久久国产| 日本91视频免费播放| av又黄又爽大尺度在线免费看| 丁香六月天网| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 黄色毛片三级朝国网站| 国产成人精品在线电影| 日日摸夜夜添夜夜爱| 男女下面插进去视频免费观看| 97人妻天天添夜夜摸| 国产爽快片一区二区三区| 欧美亚洲 丝袜 人妻 在线| 啦啦啦啦在线视频资源| 18禁黄网站禁片午夜丰满| 亚洲伊人色综图| 麻豆国产av国片精品| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 精品人妻一区二区三区麻豆| 国产精品一区二区免费欧美 | 精品国产一区二区三区四区第35| 久久久久视频综合| 国产精品一二三区在线看| 十分钟在线观看高清视频www| 日韩 欧美 亚洲 中文字幕| 大片电影免费在线观看免费| 国产深夜福利视频在线观看| 少妇人妻久久综合中文| 黑人巨大精品欧美一区二区蜜桃| 亚洲av日韩在线播放| a级毛片黄视频| 亚洲国产欧美一区二区综合| www日本在线高清视频| av网站免费在线观看视频| 我要看黄色一级片免费的| 国产精品国产三级专区第一集| 一边摸一边做爽爽视频免费| 免费看十八禁软件| 欧美av亚洲av综合av国产av| 国产激情久久老熟女| 国产xxxxx性猛交| 91精品伊人久久大香线蕉| 精品国产乱码久久久久久男人| 老汉色∧v一级毛片| av福利片在线| 极品少妇高潮喷水抽搐| 国产高清videossex| 亚洲第一青青草原| 久久久久国产一级毛片高清牌| 巨乳人妻的诱惑在线观看| 欧美国产精品一级二级三级| 亚洲伊人色综图| 精品久久蜜臀av无| 国产又爽黄色视频| 99热网站在线观看| 国产xxxxx性猛交| 婷婷丁香在线五月| 九草在线视频观看| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 免费在线观看日本一区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av欧美aⅴ国产|