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

    潮土磷素累積流失風(fēng)險及環(huán)境閾值

    2020-12-25 07:44:14張乃明張淑香
    關(guān)鍵詞:潮土磷素徑流

    劉 娟,張乃明,張淑香

    潮土磷素累積流失風(fēng)險及環(huán)境閾值

    劉 娟1,3,張乃明2,3※,張淑香4

    (1. 云南農(nóng)業(yè)大學(xué)植物保護(hù)學(xué)院,昆明 650201;2. 云南農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,昆明 650201;3. 云南省土壤培肥與污染修復(fù)工程實(shí)驗(yàn)室,昆明 650201;4. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081)

    潮土是中國分布比較廣、施肥強(qiáng)度大的典型耕作土壤,潮土中磷素累積與流失對區(qū)域水環(huán)境的污染風(fēng)險不容忽視。該研究在潮土面積最大的河南省采集磷素水平不同的典型潮土作為供試土壤,采用人工模擬降雨及土柱模擬試驗(yàn)方法,通過測定土壤中Olsen-P和溶解態(tài)活性磷CaCl2-P含量以及徑流或淋濾液中各形態(tài)磷濃度,研究了潮土中磷素隨地表徑流和下滲流失特征,并通過分段線性模型對潮土的磷素環(huán)境閾值進(jìn)行擬合。結(jié)果表明:1)不同形態(tài)磷在潮土土壤剖面中均有一定程度的累積,土壤Olsen-P和CaCl2-P含量表現(xiàn)為高磷最大,中磷次之,低磷最小,而磷吸持指數(shù)值表現(xiàn)為低磷最大,中磷次之,高磷最小。從磷素的剖面分布來看,低磷和中磷水平潮土Olsen-P和CaCl2-P含量隨著土壤深度的增加而降低,而高磷水平的潮土Olsen-P和CaCl2-P含量在20~40 cm土層含量最高。2)不同磷水平潮土徑流中總磷(Total Phosphorus,TP)、可溶性總磷(Total Dissolved phosphorus,TDP)和顆粒磷(Particulate Phosphorus,PP)濃度和流失量大小表現(xiàn)為高磷最高,中磷和低磷水平土壤次之,潮土徑流流失以PP為主。3)低磷和中磷水平潮土淋濾液中的各形態(tài)磷濃度和流失量隨著土層深度的增加而降低,而在高磷水平的潮土淋濾液中,20~40 cm土層淋濾液中磷濃度和流失量要顯著高于其他土層,在整個土壤剖面磷素濃度隨著土層深度的增加呈現(xiàn)先上升后下降的趨勢,潮土淋濾流失以TDP為主,其中,高磷和低磷水平潮土以可溶性有機(jī)磷占主導(dǎo),而中磷水平潮土以鉬酸鹽反應(yīng)磷占主導(dǎo)。4)通過分段回歸模型將不同含磷水平潮土的水溶性磷與土壤中Olsen-P含量進(jìn)行擬合,得出潮土土壤磷素環(huán)境閾值為24.65 mg/kg,研究還表明徑流和滲漏液中TP濃度與土壤CaCl2-P含量呈顯著正相關(guān),因此可通過測定CaCl2-P來預(yù)測并判斷土壤磷素流失風(fēng)險。

    磷;土壤;徑流;潮土;環(huán)境閾值;流失風(fēng)險

    0 引 言

    潮土是中國北方重要耕地土壤,面積約1 267萬hm2,占中國耕地面積的9.39%[1]。潮土區(qū)最普遍的種植制度為夏玉米-冬小麥一年兩熟制[2],這種復(fù)種指數(shù)高、利用強(qiáng)度大的種植模式對潮土本身養(yǎng)分消耗量很大[3]。其次,由于潮土區(qū)雨量分布不均,易受旱、澇、和次生鹽漬化威脅,同時由于黃河多次泛濫和改道沉積的作用,導(dǎo)致潮土肥力水平低[4]。在各項(xiàng)肥力指標(biāo)中,土壤磷素是潮土生產(chǎn)力最重要的限制因子。雖然潮土全磷含量較高,但有效磷在大部分地區(qū)含量較低,為了使小麥和玉米獲得高產(chǎn),需要在農(nóng)業(yè)生產(chǎn)中大量施用磷肥,而高頻率、過量的施肥使得潮土區(qū)土壤形成了巨大的磷庫,潮土本身屬于輕質(zhì)屬性,且潮土區(qū)屬于暖溫帶半濕潤型季風(fēng)氣候,雨季多集中于7—9月[5],磷素累積引起的環(huán)境問題不容忽視。

    土壤中磷素遷移到水體中主要包括3種途徑,即地表徑流、侵蝕和淋濾,而地表徑流和淋濾通常被認(rèn)為是農(nóng)田土壤中磷素遷移到水環(huán)境最主要的途徑[6-7]。由于農(nóng)田土壤磷素流失為主的農(nóng)田非點(diǎn)源污染已成為水環(huán)境污染的一個重要來源[8-9],使得磷素非點(diǎn)源污染研究的重要性更加凸顯[10-11]。過去在黃淮海地區(qū)的潮土中,雖然土壤全磷含量較高,但有效磷在大部分地區(qū)呈虧缺狀態(tài),在第二次土壤普查時潮土有效磷甚至不足5 mg/kg,磷素流失的風(fēng)險小,因此人們對潮土的研究主要集中在潮土不同肥力土壤的基礎(chǔ)地力、長期施肥對潮土磷素利用與累積、長期不同施肥模式潮土肥力演變規(guī)律及其持續(xù)利用以及夏玉米-冬小麥輪作系統(tǒng)下,玉米和小麥兩大主要作物的Olsen-P農(nóng)學(xué)閾值研究等方面[5,12-13]。但是隨著施磷量提高以及磷肥不斷投入,土壤中的磷素含量逐漸由虧轉(zhuǎn)為盈,導(dǎo)致土壤磷素大量累積[14],潮土區(qū)土壤有效磷平均值已經(jīng)上升到了27.9 mg/kg[15],且在潮土區(qū)小麥-玉米輪作體系下,小麥和玉米Olsen-P農(nóng)學(xué)閾值分別為13.1和7.5 mg/kg[12],這就表明潮土區(qū)大多數(shù)土壤中的Olsen-P含量已經(jīng)遠(yuǎn)遠(yuǎn)地超過作物生長的需要,繼續(xù)施肥不僅不能提高作物產(chǎn)量,還會增加土壤磷素隨徑流和淋濾到水環(huán)境中的風(fēng)險,對環(huán)境造成潛在威脅[16]。

    磷素在土壤中累積引起磷素流失的報(bào)道很多,目前關(guān)于農(nóng)田磷素流失的研究主要是以地表徑流流失規(guī)律的探究為主,在一些富磷的砂土或結(jié)構(gòu)發(fā)育良好的黏土中也有磷淋濾的報(bào)道[17]。已有研究表明,土壤中磷的流失不僅受到土壤有效磷含量的影響,還受到土壤理化性質(zhì),例如土壤有機(jī)質(zhì)、pH、孔隙度、陽離子交換量,土壤容重、交換性鈣、鎂、土壤質(zhì)地以及土壤中鐵、鋁和碳酸鈣含量的影響[18]。其中,土壤磷素水平是影響耕地土壤磷素流失的重要原因之一,當(dāng)土壤有效磷含量超過其環(huán)境閾值時,土壤徑流和滲漏液中磷含量將明顯增加,從而對水環(huán)境構(gòu)成威脅[19]。關(guān)于坡耕地紅壤、紫色土等磷素遷移與流失規(guī)律以及太湖流域水稻土徑流和滲濾液中磷素組分的研究已經(jīng)有大量報(bào)道[17],而針對潮土在土壤中累積引起的徑流和淋濾遷移和土壤磷素環(huán)境閾值以及徑流或淋濾液中全磷含量與土壤CaCl2-P之間的關(guān)系方面的研究還很少。本文選取潮土作為研究對象,通過人工模擬降雨試驗(yàn)和土柱淋濾試驗(yàn)方法,探討不同磷水平潮土徑流和淋濾液磷素流失狀況,并在此基礎(chǔ)上提出土壤磷素環(huán)境閾值,為潮土區(qū)水環(huán)境保護(hù)和土壤磷素管理提供重要的科學(xué)依據(jù)。

    1 材料與方法

    1.1 土壤樣品采集

    土壤樣品取自河南省鄭州市農(nóng)科院試驗(yàn)地(34o47′25"N,113o 40′42" E),該區(qū)域?qū)儆跍貛駶?、半濕潤季風(fēng)氣候,冬春干旱少雨、夏季炎熱濕潤,年平均降水量665 mm,雨季多集中在7—9月。分析該區(qū)域降雨資料發(fā)現(xiàn),潮土地區(qū)最大降雨強(qiáng)度為78 mm/h[20]。該區(qū)域地形較為平坦,坡度在5°~15°之間,試驗(yàn)所用土壤為潮土(淡半水成土),成土母質(zhì)為黃河沖積物,土壤質(zhì)地為砂質(zhì)黏土,于2017年4月,按照“隨機(jī)”多點(diǎn)混合的原則,采集耕層(0~20 cm)及剖面土壤樣品,土壤剖面按層次采樣,自下而上分層采取,每個層次分為0~20、>20~40、>40~60和>60~80 cm共4個層面。

    根據(jù)耕層(0~20 cm)土壤磷素含量狀況,將采集的潮土大致分為低、中、高3種磷素水平的土壤。低磷水平潮土連續(xù)27 a每年平均施磷肥(普通過磷酸鈣,P2O5含量12.05%,下同)75 kg/hm2,中磷水平潮土每年平均施磷肥175 kg/hm2,而高磷水平潮土每年平均施磷肥265 kg/hm2;磷肥施用量是根據(jù)潮土區(qū)施肥現(xiàn)狀進(jìn)行設(shè)置,在研究區(qū)域內(nèi)具有典型性和代表性。將采集回來的土壤樣品,剔除土樣中的石塊、粗有機(jī)物和侵入體等,風(fēng)干、磨碎,過2?mm孔徑篩,分層充分混勻后備用。供試土壤的主要基本性質(zhì)如表1所示。

    表1 供試土壤0~20 cm土層的理化性質(zhì)

    1.2 徑流試驗(yàn)設(shè)計(jì)

    徑流試驗(yàn)于2017年6—7月在云南省農(nóng)業(yè)節(jié)水工程技術(shù)研究中心模擬降雨大廳進(jìn)行,設(shè)備采用NLJY?10型室內(nèi)人工模擬降雨控制系統(tǒng)(南林電子生產(chǎn)),噴頭高度為16 m,噴頭類型為散噴型,有效降雨面積約40 m2。該人工模擬降雨控制系統(tǒng)的精度主要是通過雨滴分布(小、中、大雨滴所占比例數(shù))及雨滴降落終速來控制,采用色斑法測定雨滴直徑大小,采用高速攝影法校驗(yàn)雨滴降落終速;在進(jìn)行徑流試驗(yàn)以前,進(jìn)行降雨均勻度的測定及降雨強(qiáng)度的標(biāo)定,其主要是通過計(jì)算降雨期間地面設(shè)置的雨量筒(8個)觀測到的降雨量的平均值實(shí)現(xiàn),率定后降雨均勻系數(shù)大于95%,雨滴降落終速可達(dá)到自然雨滴的99%,降雨特性接近于天然降雨。

    試驗(yàn)采用70 cm(長)×28 cm(寬)×22 cm(高)的土槽作為盛土容器,然后向土槽中填裝等量0~20 cm耕層土,為了確保填裝土壤的均勻,采用分層填裝的方法即每隔5 cm填裝一層,在填裝下層之前,抓毛上一層已填進(jìn)土壤的表層,以防止各土層之間出現(xiàn)分層現(xiàn)象,同時,通過分層填裝與壓實(shí)相結(jié)合的方法控制土壤容重,確保填裝好的土壤容重保持在1.3 g/m3左右?;诔蓖羺^(qū)坡耕地現(xiàn)狀及降雨發(fā)生頻率,本試驗(yàn)坡度設(shè)計(jì)為8°,降雨強(qiáng)度設(shè)計(jì)為30 mm/h,每7 d降雨1次,每次降雨32 min,每4 min用徑流桶收集1次徑流,并測定徑流體積,降雨試驗(yàn)結(jié)束后取各時段采集的徑流樣進(jìn)行測定,每個磷素水平設(shè)置3個重復(fù)。

    1.3 淋濾試驗(yàn)設(shè)計(jì)

    淋濾試驗(yàn)于2017年8—10月在云南農(nóng)業(yè)大學(xué)化學(xué)樓進(jìn)行,將PVC管(長80 cm,直徑10 cm)截取為4段,每段20 cm,用尼龍網(wǎng)包裹每段PVC管的底端,然后分別將3種含磷水平的潮土分層填充到土柱中,為了可以方便地分離和恢復(fù)完整的土柱,該研究采用PVC接口(直徑為11 cm)將每段土柱連接成高80 cm的土柱桶。土柱桶底部填充約2 cm厚的石英砂,石英砂和土壤之間由1張48m的紗網(wǎng)分隔,然后用PVC封蓋(直徑為11 cm)蓋住土柱底部。為了能夠方便地收集淋濾液,在每個封蓋上鉆1個小孔,小孔直徑約0.5 cm,用橡膠塞將排水管固定在封蓋上收集淋濾液,最后將土柱直立地放在架子上。淋濾試驗(yàn)開始以前,向土壤中加入超純水達(dá)到飽和持水量后靜置24 h,淋濾試驗(yàn)連續(xù)進(jìn)行42 d,每7 d灌水1次,每次灌水約300 mL,在灌水后的24 h內(nèi)收集淋濾液,測定淋濾液體積。為避免土壤受到干擾,在土壤表面放置1張無磷濾紙,同時,濾紙還可使水樣勻速下滲,避免產(chǎn)生邊際效應(yīng),每個磷素水平潮土設(shè)置3個重復(fù)。

    1.4 分析項(xiàng)目與方法

    溶液中總磷(Total Phosphorus,TP)濃度采用過硫酸鉀?鉬藍(lán)比色法(GB 11893—89)測定;可溶性總磷(Total Dissolved Phosphorus,TDP)濃度采用0.45m濾膜過濾,過硫酸鉀?鉬藍(lán)比色法(GB 11893—89)測定;鉬酸鹽反應(yīng)磷(Molybdate Reactive Phosphorus,MRP)采用0.45m濾膜過濾,鉬藍(lán)比色法測定;顆粒磷(Particulate Phosphorus,PP)為TP與TDP的差值;可溶性有機(jī)磷(Dissolved Organic Phosphorus,DOP)為TDP與MRP的差值;

    土壤pH值采用酸度計(jì)(NY-T 1121.2—2006)測定;有機(jī)質(zhì)采用重鉻酸鉀氧化-容量法(NY-T 1121.6—2006)測定;堿解氮采用堿解擴(kuò)散法測定;有效磷(Olsen-P)采用鉬銻抗比色法(HJ704-2014)測定;速效鉀采用1 mol/L的中性NH4OAc(pH值為7)溶液浸提,火焰光度法測定;全磷采用堿熔?鉬銻抗比色法(HJ632-2011)測定;溶解態(tài)活性磷(CaCl2-P)采用0.01 mol/L氯化鈣溶液浸提(土水質(zhì)量比1∶5),鉬銻抗比色法測定;土壤磷吸持指數(shù)(Phosphate Sorption Index,PSI):在水土質(zhì)量比為10∶1的條件下,按每克土中加入1.5 mg磷酸鹽,平衡后測得土壤吸磷量與平衡溶液磷濃度,采用計(jì)算式(1)[21]計(jì)算。

    PSI=100/lg(1)

    式(1)中為土壤吸磷量,mg/g;為平衡溶液中磷濃度,mol/L。

    1.5 數(shù)據(jù)統(tǒng)計(jì)與分析方法

    采用Excel 2007和Sigmaplot 10.0對數(shù)據(jù)進(jìn)行分析,通過軟件SPSS 19.0對所有數(shù)據(jù)進(jìn)行單因素方差分析和顯著性檢驗(yàn),檢驗(yàn)方法采用ANOVA法,采用OriginPro 9.0進(jìn)行作圖,所有結(jié)果數(shù)據(jù)均采用3個重復(fù)試驗(yàn)所獲數(shù)據(jù)的平均值。

    為了計(jì)算場降雨徑流磷素流失量,采用公式

    式中為磷素流失量,mg/m2;為取樣次數(shù);c為各時段采集徑流樣的磷素濃度,mg/L;v為各時段采集徑流樣的體積,L;S為土槽的面積,m2。

    采用分段線性模型(Split-line Model)計(jì)算潮土磷素流失的環(huán)境閾值。根據(jù)Olsen-P與CaCl2-P之間的關(guān)系圖,采用數(shù)學(xué)方法,以不偏離突變點(diǎn)為基礎(chǔ),分別擬合低和高Olsen-P含量與CaCl2-P之間的分段線性回歸方程,2個方程的交點(diǎn)即為土壤磷素環(huán)境閾值[17]。

    2 結(jié)果與分析

    2.1 潮土中有效磷和易溶解態(tài)磷及吸持指數(shù)

    不同磷水平潮土中Olsen-P累積狀況如圖1a所示,由圖可知,3個磷素水平潮土在土壤剖面中Olsen-P含量(質(zhì)量分?jǐn)?shù),下同)分別為3.0~8.5、11.0~52.5和8.2~33.4 mg/kg,不同磷水平潮土Olsen-P含量在各土層中表現(xiàn)出一致規(guī)律,即高磷水平最高,中磷水平次之,而低磷水平最低,3個磷素水平之間呈顯著性差異。從土壤剖面分布來看,低磷和中磷土壤中Olsen-P含量隨著土壤剖面深度的增加呈降低趨勢,而高磷水平土壤在20~40 cm土層最高。全國第二次土壤普查養(yǎng)分分級標(biāo)準(zhǔn)將土壤中Olsen-P含量分為6個等級,將3個磷素水平潮土表層土壤(0~20 cm)中的Olsen-P含量與全國第二次土壤普查的養(yǎng)分分級標(biāo)準(zhǔn)進(jìn)行比較,高磷水平潮土達(dá)到一級標(biāo)準(zhǔn),土壤Olsen-P極高;中磷水平潮土達(dá)到二級標(biāo)準(zhǔn),土壤Olsen-P很高;低磷水平潮土達(dá)到四級標(biāo)準(zhǔn),土壤Olsen-P中等,這說明潮土表層均出現(xiàn)了一定程度的累積。

    注:不同大寫字母表示同一土層不同磷素水平間P累積差異顯著(P<0.05),不同小寫字母表示同一磷素水平不同土層間P累積差異顯著(P<0.05)。

    不同磷水平潮土中CaCl2-P累積狀況如圖1b所示,由圖可知,土壤剖面中CaCl2-P含量范圍為1.01~4.20 mg/kg。不同磷水平潮土CaCl2-P含量在各土層中存在差異,在0~20 cm土層,高磷和中磷水平土壤CaCl2-P含量顯著高于低磷水平土壤,而高磷和中磷水平土壤之間差異不顯著;在20~40 cm土層,3個磷素水平土壤之間CaCl2-P含量高磷最高,中磷次之,低磷最低,且各磷水平之間呈顯著性差異;而在40~60 cm和60~80 cm土層,3個磷素水平土壤之間CaCl2-P含量未達(dá)到顯著性差異。從土壤剖面分布來看,低磷土壤CaCl2-P含量雖然隨著土壤剖面深度的增加呈降低趨勢,但各層次土壤中的CaCl2-P含量差異性不顯著。中磷土壤在整個土壤剖面CaCl2-P隨著土壤剖面深度的增加,CaCl2-P含量降低,在40~60 cm和60~80 cm土層中CaCl2-P含量差異不顯著。高磷土壤中CaCl2-P質(zhì)量分?jǐn)?shù)范圍為1.15~4.20 mg/kg,CaCl2-P含量在20~40 cm土層最高,而在40~60 cm和60~80 cm土層CaCl2-P含量未達(dá)到顯著性差異,在整個剖面上CaCl2-P含量呈現(xiàn)先上升,再下降的趨勢。

    土壤PSI反映了土壤磷吸持或釋放潛力,一般用于表征土壤固相磷向液相釋放的可能性。對不同磷水平潮土的PSI值進(jìn)行分析計(jì)算得出,各層次土壤的PSI值的范圍為4.38~22.42,其中,低、中、高磷水平潮土在土壤剖面中的PSI值分別為4.66~22.42、6.22~10.58和4.75~9.30,在整個土壤剖面中,低磷水平土壤磷吸持指數(shù)最大,中磷水平土壤次之,高磷水平土壤磷吸持指數(shù)最小,PSI越大,說明土壤吸持磷的能力越強(qiáng),而PSI越小,說明土壤中的磷越容易發(fā)生遷移。從土壤剖面來看,不同磷水平潮土PSI值隨著土壤剖面深度的增加表現(xiàn)出不同的趨勢,高磷和低磷水平潮土在整個土壤剖面隨著土壤深度的增加,PSI值降低,而中磷水平潮土在整個土壤剖面變化幅度整體較小,大致呈先下降再上升,再下降的趨勢(圖2)。

    圖2 不同磷水平潮土各土層磷吸持指數(shù)(PSI)

    2.2 不同磷水平潮土徑流中各形態(tài)磷濃度及流失量

    不同磷水平潮土徑流中各形態(tài)磷濃度如圖3a所示,表層徑流中TP濃度為0.146~0.416 mg/L,可溶性總磷TDP濃度為0.032~0.086 mg/L,顆粒磷PP濃度為0.109~0.330 mg/L,不同磷水平潮土TP、TDP和PP濃度大小表現(xiàn)為高磷水平潮土最高,中磷和低磷水平潮土次之。通過方差分析得出,在第1場降雨和第2場降雨中,高磷水平潮土徑流中TP濃度要顯著高于低磷和中磷水平,而低磷和中磷水平之間TP濃度未達(dá)到顯著性差異;在第3場和第4場降雨中,3種磷水平潮土徑流中TP濃度均達(dá)到了顯著性差異。從徑流中各形態(tài)磷占TP的比例來看,不同磷水平潮土TDP占TP的比例為9.76%~34.02%,而PP占TP的比例為65.98%~90.24%,PP所占比例明顯高于TDP,說明3種磷水平潮土徑流流失都以顆粒態(tài)為主。

    不同磷水平潮土徑流中各形態(tài)磷流失量如圖3b所示,3種磷水平潮土徑流中TP流失量為2.05~6.05 mg/m2,TDP流失量為0.31~1.29 mg/m2,PP流失量為1.48~5.16 mg/m2,不同磷水平潮土徑流TP、TDP和PP流失量大小表現(xiàn)為高磷水平土壤最大,中磷和低磷水平土壤次之,3種磷素水平潮土徑流中TP流失量的大小順序與土壤表層有效磷含量具有一致性,這說明土壤表層磷素累積狀況可能會影響徑流磷素流失量的大小。從各形態(tài)磷占總磷的比例來看,不同磷水平潮土TDP流失量占TP流失量的10.31%~37.46%,PP流失量占TP流失量的62.54%~89.69%,由此可見,潮土徑流中土壤表層磷素流失以顆粒態(tài)為主。

    注:Ⅰ~Ⅳ代表不同降雨場次;不同小寫字母表示同一降雨場次不同磷水平之間潮土徑流液中磷素濃度或流失量差異顯著(P<0.05)。

    2.3 不同磷水平潮土淋濾液中各形態(tài)磷濃度及流失量

    不同磷水平潮土淋濾液中各形態(tài)磷濃度如圖4a所示,由圖可知,不同磷水平潮土淋濾液中各形態(tài)磷濃度差異大,TP、TDP、PP、MRP和DOP濃度分別為0.114~0.995、0.084~0.660、0.031~0.335、0.027~0.378和0.050~0.321 mg/L,潮土淋濾流失以TDP為主,在TDP中,高磷和低磷水平土壤以DOP為主,而中磷水平土壤以MRP為主。從淋濾液中各形態(tài)磷的剖面分布來看,低磷和中磷水平潮土淋濾液中的各形態(tài)磷濃度隨著剖面深度的增加磷素濃度降低,在表層(0~20 cm)淋濾液中的磷素濃度最高;而在高磷水平的潮土淋濾液中,各形態(tài)磷濃度在20~40 cm土層最高。3種磷素水平潮土底層(60~80 cm)土壤淋濾液中的磷素濃度普遍較低,但3種磷水平土壤淋濾液中TP濃度均超過了水體富營養(yǎng)化閾值(0.02 mg/L)[22]。

    注:不同大寫字母表示同一土壤剖面不同磷水平之間差異顯著(P<0.05);不同小寫字母代表同一磷素水平不同剖面之間差異顯著(P<0.05);

    不同磷水平潮土淋濾液中各形態(tài)磷流失量如圖4b所示,由圖可知,不同磷水平潮土淋濾液中TP、TDP、PP、MRP和DOP流失量分別為1.97~14.84、1.29~12.48、0.41~2.99、0.42~7.60和0.51~4.92 mg/m2。潮土淋濾液中以TDP為主,在TDP中,不同含磷水平土壤占主導(dǎo)的磷形式又存在差異,其中,高磷和中磷水平潮土淋濾流失以MRP為主,而低磷水平潮土淋濾流失以DOP為主。從淋濾液中各形態(tài)磷的剖面分布來看,低磷和中磷水平潮土淋濾液中的各形態(tài)磷流失量隨著剖面深度的增加而降低,而在高磷水平的潮土淋濾液中,20~40 cm土層的磷流失量要顯著高于其他土層。

    2.4 土壤磷素的環(huán)境閾值擬合

    為了更好地預(yù)測磷的損失潛力,利用分段線性模型(Split-line Model),將不同含磷水平潮土的Olsen-P含量與土壤中CaCl2-P進(jìn)行回歸擬合(如圖5),得出潮土Olsen-P含量與CaCl2-P之間并不是1條直線關(guān)系,而是2線斜率明顯不同的直線,其相關(guān)方程如下:

    式中x為Olsen-P含量,mg/kg;y為土壤中CaCl2-P含量,mg/kg。

    通過計(jì)算得出,研究區(qū)域內(nèi)潮土土壤磷素環(huán)境閾值為24.65 mg/kg。當(dāng)土壤Olsen-P含量小于或等于土壤磷素環(huán)境閾值時,CaCl2-P含量隨Olsen-P增加而增加的幅度小,而當(dāng)土壤中Olsen-P的含量大于土壤磷素環(huán)境閾值時,CaCl2-P含量隨Olsen-P增加而增加的幅度大。

    進(jìn)一步分析潮土徑流液或淋濾液中TP濃度與土壤CaCl2-P之間的關(guān)系,發(fā)現(xiàn)線性方程可以較好地?cái)M合潮土徑流或淋濾液中TP濃度與土壤CaCl2-P含量之間的關(guān)系,通過線性回歸分析得出徑流液和淋濾液中TP濃度與土壤CaCl2-P關(guān)系如圖6所示。

    徑流中TP濃度與土壤CaCl2-P含量的相關(guān)方程為

    =0.06+0.11 (2=0.80,<0.01,樣本量為12) (4)

    淋濾液中TP濃度與土壤CaCl2-P含量的相關(guān)方程為

    =0.21+0.03(2=0.79,<0.01,樣本量為36) (5)

    式中為土壤CaCl2-P含量,mg/kg;為土壤中TP濃度,mg/L。

    由擬合方程發(fā)現(xiàn),土壤中CaCl2-P含量與徑流或淋濾液TP濃度存在正相關(guān)關(guān)系,隨著土壤中CaCl2-P含量增加,徑流或淋濾液TP濃度也隨之增加。通過測定土壤中CaCl2-P來預(yù)測徑流或淋濾液中TP濃度,可以評估土壤磷素流失對水體的造成的污染。

    圖6 潮土徑流液和淋濾液中總磷濃度(TP)與CaCl2-P之間的關(guān)系

    3 討 論

    3.1 不同磷水平潮土磷素累積流失特征

    潮土區(qū)是中國重要的小麥-玉米主產(chǎn)區(qū),這種復(fù)種指數(shù)高、利用強(qiáng)度大的種植模式對潮土本身養(yǎng)分消耗量很大,需要在農(nóng)業(yè)生產(chǎn)中大量施用磷肥,而高頻率、過量的施肥使得潮土區(qū)土壤形成了巨大的磷庫。在本研究中發(fā)現(xiàn),潮土表層土壤(0~20 cm)已經(jīng)出現(xiàn)了一定程度的磷素累積,將3個磷素水平潮土表層土壤(0~20 cm)中的Olsen-P含量與全國第二次土壤普查的養(yǎng)分分級標(biāo)準(zhǔn)進(jìn)行比較,高磷水平潮土達(dá)到一級標(biāo)準(zhǔn),土壤Olsen-P極高;中磷水平潮土達(dá)到二級標(biāo)準(zhǔn),土壤Olsen-P很高;低磷水平潮土達(dá)到四級標(biāo)準(zhǔn),土壤Olsen-P中等,這說明潮土表層均出現(xiàn)了一定程度的累積。當(dāng)土壤中磷素累積超過環(huán)境閾值,土壤中磷素可以通過地表徑流、侵蝕以及淋濾等途徑遷移到水體環(huán)境中,進(jìn)而加劇水環(huán)境污染負(fù)荷。

    在本研究中,3種磷水平表層徑流中TP濃度為0.146~0.416 mg/L,TDP濃度為0.032~0.086 mg/L,PP濃度為0.109~0.330 mg/L,不同磷水平潮土TP、TDP和PP濃度大小表現(xiàn)為高磷水平潮土最高,中磷和低磷水平潮土次之,這與崔力拓等[23]先前通過天然降雨徑流試驗(yàn)得到的褐土地表徑流中不同形態(tài)磷的含量隨土壤磷素水平的提高也相應(yīng)地增加的結(jié)論一致,并且在其他土壤類型如黃壤、紫色土等也發(fā)現(xiàn)了類似的觀察結(jié)果[24-25],這說明徑流液中磷素濃度與土壤本身磷素累積狀況有密切關(guān)系。本實(shí)驗(yàn)室還通過人工模擬降雨試驗(yàn),比較了黑土、潮土和紅壤中磷素隨地表徑流流失的特征,發(fā)現(xiàn)黑土、潮土和紅壤表層徑流中磷素濃度和磷素流失量均表現(xiàn)為潮土最大,紅壤次之,黑土最小[26]。磷素在潮土中的擴(kuò)散性要強(qiáng)于其他2種土壤,主要是因?yàn)槌蓖了趨^(qū)域?yàn)橹袊饕霓r(nóng)耕區(qū),磷肥的逐年投入及施磷量的提高,加劇土壤中磷素累積狀況,潮土為砂質(zhì)黏土,土壤中黏粒含量較黑土和紅壤少,所以當(dāng)徑流發(fā)生時,潮土比其他2種土壤更容易發(fā)生遷移。從各形態(tài)的比例來看,本研究得出顆粒態(tài)磷是潮土徑流流失的主要形態(tài),這與秦華等[27]研究得出徑流作用下PP濃度占總磷濃度80%以上以及左繼超等[28]研究得出的地表徑流磷素輸出以PP為主的結(jié)論相同,在降雨的過程中,徑流中磷素的流失不僅與土壤本身含磷狀況密切相關(guān),還可能與坡度、降雨強(qiáng)度、徑流量等因素相關(guān),所以,在季風(fēng)區(qū),應(yīng)結(jié)合土壤本身性質(zhì)以及該區(qū)域降雨條件,對養(yǎng)分流失進(jìn)行管理。

    在本研究中,不同磷水平潮土淋濾液中TP濃度為0.114~0.995 mg/L,TDP濃度為0.084~0.660 mg/L,PP濃度為0.031~0.335 mg/L,MRP濃度為0.027~0.378 mg/L,DOP濃度為0.050~0.321 mg/L。低磷和中磷水平潮土淋濾液中的各形態(tài)磷含量隨著剖面深度的增加,磷素濃度降低,而高磷水平潮土淋濾液中磷素濃度隨著土層深度的增加,呈現(xiàn)先上升后下降的趨勢,這種現(xiàn)象的出現(xiàn)與土壤中Olsen-P含量的高低密切相關(guān),在高磷土壤中,20~40 cm土層徑流液中的磷素濃度要高于土壤表層,主要是因?yàn)楸旧硎┝姿捷^高,作物吸收不了過多的磷,在降雨的作用下,由磷素向下淋溶所致。本試驗(yàn)室還通過土柱淋溶試驗(yàn)[1],比較黑土、潮土、紅壤和水稻土淋溶流失差異,發(fā)現(xiàn)不同類型土壤TP流失量以潮土最大,水稻土和黑土次之,紅壤最小,這主要是由于長期耕作條件下,潮土本身磷含量累積過高,以及潮土砂質(zhì)黏土的土壤質(zhì)地,黏粒含量比其他幾種土壤少造成的。本研究還發(fā)現(xiàn),不同磷水平潮土淋濾流失以TDP為主,這與項(xiàng)大力等[29]在塿土以及陸欣欣等[30]的得出的淋濾流失以TDP為主的結(jié)論一致,這主要是因?yàn)榱転V與徑流磷素流失不同,淋濾流失受外部因素(坡度、降雨強(qiáng)度、徑流量等)的影響要小于土壤本身性質(zhì)的影響,土壤中磷素含量、土壤質(zhì)地等在潮土淋濾過程中發(fā)揮著極其重要的作用。

    3.2 潮土土壤磷素環(huán)境閾值的科學(xué)性及應(yīng)用

    不同類型土壤磷素環(huán)境閾值差異很大,鐘曉英等[31]通過土壤培養(yǎng)試驗(yàn),得出中國11種土壤類型23個土壤的土壤磷素環(huán)境閾值在29.96~156.78 mg/kg之間變化。聶敏等[32]通過室內(nèi)模擬試驗(yàn),發(fā)現(xiàn)中國16種可變電荷土壤(紅壤、黃壤等)的土壤磷素環(huán)境閾值在56~123 mg/kg之間變化,同時,趙小蓉等[33]還通過研究土壤磷素環(huán)境閾值與土壤理化性質(zhì)的關(guān)系時發(fā)現(xiàn),土壤磷素環(huán)境閾值與土壤pH呈拋物線關(guān)系,當(dāng)pH值大于6時,土壤磷素環(huán)境閾值隨著pH的升高而減小。在本研究中,利用分段回歸模型(Split-line Model),將不同含磷水平潮土的水溶性磷與土壤中Olsen-P含量進(jìn)行回歸擬合,分析得出潮土的土壤磷素環(huán)境閾值為24.65 mg/kg,潮土本身為石灰性土壤,土壤pH值在8.01~8.25之間,所以與其他土壤類型相比,潮土土壤磷素環(huán)境閾值可能偏低,而潮土中土壤有機(jī)質(zhì)、土壤質(zhì)地、無定形氧化鐵、鋁含量、碳酸鈣等是否對土壤磷素環(huán)境閾值造成影響還需要進(jìn)一步深入研究。

    一般而言,當(dāng)水體中的含磷量超過了水體富營養(yǎng)化閾值(0.02 mg/L)就可能造成水體中藻類的大量繁殖。在本研究中發(fā)現(xiàn),當(dāng)潮土Olsen-P含量超過土壤磷素環(huán)境閾值(即:24.65 mg/kg)時,從土體中排出的水進(jìn)入河流或湖泊中,極可能對水體造成污染。因此,在潮土區(qū)域,應(yīng)該根據(jù)潮土目前的磷素累積狀況,適當(dāng)減少施磷量,在維持作物產(chǎn)量的同時,降低P對水環(huán)境造成污染的風(fēng)險,在保證作物不減產(chǎn)的條件下,適當(dāng)減少施磷量還可以在一定程度提高土壤中已累積磷素的再利用[34],這對于農(nóng)業(yè)生產(chǎn)與生態(tài)環(huán)境保護(hù)都具有重要的現(xiàn)實(shí)意義。本研究中,還建立了徑流或淋濾液中總磷與CaCl2-P之間的關(guān)系,發(fā)現(xiàn)潮土徑流或淋濾液TP濃度與土壤中CaCl2-P含量存在正相關(guān)關(guān)系,徑流或淋濾液TP濃度會隨著土壤CaCl2-P含量的增大迅速增加,因此,可以通過測定土壤中CaCl2-P的含量來預(yù)測徑流和淋濾液中TP濃度,從而對土壤磷素遷移對水環(huán)境造成的污染進(jìn)行評價。

    4 結(jié) 論

    1)土壤剖面中的Olsen-P含量(質(zhì)量分?jǐn)?shù))為3.0~33.4 mg/kg,溶解態(tài)活性磷(CaCl2-P)含量為1.01~4.20 mg/kg,土壤磷吸持指數(shù)(PSI)值為4.38~22.42,Olsen-P和溶解態(tài)活性磷(CaCl2-P)含量表現(xiàn)為高磷最高,中磷次之,低磷最低,而PSI值表現(xiàn)為低磷最高,中磷次之,高磷最低。

    2)不同磷水平潮土徑流中總磷、可溶性總磷和顆粒磷濃度和流失量大小為高磷最高,中磷和低磷水平土壤次之,潮土徑流流失以顆粒態(tài)為主。

    3)低磷和中磷水平潮土淋濾液中的各形態(tài)磷濃度和流失量隨著土層深度的增加而降低,而在高磷水平的潮土淋濾液中,20~40 cm土層淋濾液中磷濃度和流失量要顯著高于其他土層,潮土淋濾流失以可溶態(tài)為主。

    4)通過分段回歸模型將不同含磷水平潮土的水溶性磷與土壤中Olsen-P含量進(jìn)行擬合,得出潮土土壤磷素環(huán)境閾值為24.65 mg/kg。同時還得出,徑流液和淋濾液中總磷濃度與土壤中CaCl2-P含量呈顯著正相關(guān),可通過測定CaCl2-P來預(yù)測從土體排出水的磷濃度。

    [1] 劉娟,包立,張乃明,等. 我國4種土壤磷素淋溶流失特征[J]. 水土保持學(xué)報(bào),2018,32(5):67-73. Liu Juan, Bao Li, Zhang Naiming, et al. Characteristics of phosphorus leaching losses in four soils in China[J]. Journal of soil and Water conservation, 2018, 32(5): 67-73. (in Chinese with English abstract)

    [2] 黃欣欣,廖文華,劉建玲,等. 長期秸稈還田對潮土土壤各形態(tài)磷的影響[J]. 土壤學(xué)報(bào),2016,53(5):779-789. Huang Xinxin, Liao Wenhua, Liu Jianling, et al. Effects of long-term straw return on various fractions of phosphorus in Fluvo-aquic soil[J]. Acta Pedologica Sinica, 2016, 53(5): 779-789. (in Chinese with English abstract)

    [3] 信秀麗,欽繩武,張佳寶,等. 長期不同施肥下潮土磷素的演變特征[J]. 植物營養(yǎng)與肥料學(xué)報(bào),2015,21(6):1514-1520. Xiu Xiuli, Qin Shengwu, Zhang Jiabao, et al. Dynamics of phosphorus in Fluvo-aquic soil under long-term fertilization[J]. Journal of P 1ant Nutrition and Fertilizer, 2015, 21(6): 1514-1520. (in Chinese with English abstract)

    [4] 申艷,段英華,黃紹敏,等. 潮土 CaCl2-P 含量對磷肥施用的響應(yīng)及其淋失風(fēng)險分析[J]. 植物營養(yǎng)與肥料學(xué)報(bào),2018,24(6):1689-1696. Shen Yan, Duan Yinghua, Huang Shaomin, et al. Response of CaCl2-P to phosphorus fertilization and leaching risk in fluvo-aquic soils[J]. Journal of Plant Nutrition and Fertilizers, 2018, 24(6): 1689-1696. (in Chinese with English abstract)

    [5] 黃紹敏,郭斗斗,張水清. 長期施用有機(jī)肥和過磷酸鈣對潮土有效磷積累與淋溶的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào),2011,22(1):97-102. Huang Shaomin, Guo Doudou, Zhang Shuiqing. Effects of long-term application of organic fertilizer and superphosphate on accumulation and leaching of Olsen-P in fluvo-aquic soil[J]. Chinese Journal of Applied Ecology, 2011, 22(1): 97-102. (in Chinese with English abstract)

    [6] Koopmans G F, McDowell R W, Chardon R W, et al. Soil phosphorus quantility-intensity relationships to predict increased soil phosphorus loss to overland and subsurface flow[J]. Chemosphere, 2002, 48: 697-687.

    [7] Simmonds B, Mcdowell R W, Condron L M, et al. Can phosphorus fertilizers sparingly soluble in water decrease phosphorus leaching loss from an acid peat soil[J/OL]. Soil Use and Management, 2016. doi:10.1111/sum.12274

    [8] Wang G, Wu B, Zhang L, et al. Role of soil erodibility in affecting available nitrogen and phosphorus losses under simulated rainfall[J]. Journal of Hydrology, 2014, 514: 180-191.

    [9] 吳漢卿,萬煒,單艷軍,等. 基于磷指數(shù)模型的海河流域農(nóng)田磷流失環(huán)境風(fēng)險評價[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(14):17-27. Wu Hanqing, Wan Wei, Shan Yanjun, et al. Environmental risk assessment of phosphorus loss from farmland based on phosphorus index model in the Haihe River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 17-27. (in Chinese with English abstract)

    [10] Kleinman P J A, Church C, Saporito L S, et al. Phosphorus leaching from agricultural soils of the delmarva peninsula, USA[J]. Journal of Environment Quality, 2015, 44(2): 524-534.

    [11] Jarvie H P, Johnson L T, Sharpley A N, et al. Increased soluble phosphorus loads to lake erie: Unintended consequences of conservation practices[J]. Journal of Environmental Quality, 2017, 46(1): 123-132.

    [12] 張水清,黃紹敏,婁翼來,等. 華北潮土區(qū)不同肥力水平土壤基礎(chǔ)地力研究[J]. 中國農(nóng)學(xué)通報(bào),2016,32(20):97-100. Zhang Shuiqing, Huang Shaomin, Lou Yilai, et al. Basic soil productivity study on different fertility levels in fluvo-aquic soil region, North China[J]. Chinese Agricultural Science Bulletin, 2016, 32(20): 97-100. (in Chinese with English abstract)

    [13] 郭斗斗,黃紹敏,張水清,等. 潮土小麥和玉米 Olsen-P農(nóng)學(xué)閾值及其差異分析[J]. 植物營養(yǎng)與肥料學(xué)報(bào),2017,23(5):1184–1190. Guo Doudou, Huang Shaomin, Zhang Shuiqing, et al. Threshold values o f soil Olsen-P for maize and wheat in fluvo-aquicsoil[J]. Journal of Plant Nutrition and Fertilizer, 2017, 23(5): 1184–1190. (in Chinese with English abstract)

    [14] 萬煒,師紀(jì)博,劉忠,等. 棲霞市蘋果園氮磷養(yǎng)分平衡及環(huán)境風(fēng)險評價[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(4):211-219. Wan Wei, Shi Jibo, Liu Zhong, et al. Nitrogen and phosphorus nutrient balance and environmental risk assessment of apple orchard in Qixia city[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 211-219. (in Chinese with English abstract)

    [15] 李少叢,萬紅友,王興科,等. 河南省潮土、砂姜黑土基本性質(zhì)變化分析[J]. 土壤,2014,46(5):920-926. Li Shaocong, Wan Hongyou, Wang Xingke, et al. Basic properties variation of fluvo-aquic soil and lime calcic black soil in Henan province[J]. Soils, 2014, 46(5): 920-926. (in Chinese with English abstract)

    [16] 王帥兵,宋婭麗,王克勤,等. 不同雨型下反坡臺階減少紅壤坡耕地氮磷流失的效果[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(13):160-169. Wang Shuaibing, Song Yali, Wang Keqin, et al. Effects of reverse-slope terrace on nitrogen and phosphorus loss in sloping farmland of red loam under different rainfall patterns[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(13): 160-169. (in Chinese with English abstract)

    [17] Li YY, Gao R, Yang R, et al. Using a simple soil column method to evaluate soil phosphorus leaching risk[J]. Clean-Soil, Air, Water, 2013, 41(11): 1100-1107.

    [18] Heckrath G, Brookes P C, Poulton P R, et al. Phosphorus leaching from containing different phosphorus concentrations in the Broadbalk experiment[J]. J Environ Qual, 1995, 24: 904-910.

    [19] Hesketh N, Brookes P C. Development of an Indicator for risk of phosphorus leaching[J]. J Environ Qual, 2000, 29: 105-110.

    [20] 杜明成,王振龍,姜翠玲,等. 淮北平原黃潮土多雨強(qiáng)變坡度產(chǎn)流產(chǎn)沙規(guī)律試驗(yàn)?zāi)M[J]. 水土保持學(xué)報(bào),2018,32(6):34-39. Du Mingcheng, Wang Zhenlong, Jiang Cuiling, et al. Simulation of runoff and sediment production regularity of different rainfall intensity and changeable slope gradients in the Yellow fluvo-aquic soil of the Huaibei plain[J]. Journal of Ecology and Rural Environment, 2018, 32(6):34-39. (in Chinese with English abstract)

    [21] 徐紅嬌,包立,張乃明,等. 滇池流域不同利用方式紅壤滲濾液的磷素形態(tài)變化[J]. 水土保持學(xué)報(bào),2015,29(3):267-271. Xu Hongjiao, Bao Li, Zhang Naiming, et al. Changes of phosphorus fractions in leachate from red soil under different land use types in Dianchi Basin[J]. Journal of Soil and Water Conservation, 2015, 29(3): 267-271. (in Chinese with English abstract)

    [22] OECD. Eutrophication of waters[M]//Monitoring, Assessment and Control. Paris, France,OECD: 1982.

    [23] 崔力拓,李志偉. 洋河流域緩坡地土壤磷素徑流輸出特征[J]. 水土保持學(xué)報(bào),2006,37(3):10-12. Cui Lituo, Li Zhiwei. Characteristics of soil phosphorus export by surface runoff from gentle slope in Yanghe watershed[J]. Journal of Soil and Water Conservation, 2006, 37(3): 10-12. (in Chinese with English abstract)

    [24] 劉方,何騰兵,錢曉剛,等. 不同利用方式下黃壤旱坡地磷素狀況及環(huán)境影響分析[J]. 土壤與環(huán)境,2002,11(3):232-236. Liu Fang, He Tengbing, Qian Xiaogang, et al. The status of phosphorus and its environmental effects analysis under various land use patterns in Yellow soil of hilly areas[J]. Soi1 and Environmenta1 Sciences, 2002, 11(3): 232-236. (in Chinese with English abstract)

    [25] 李學(xué)平,孫燕,石孝均. 紫色土稻田磷素淋失特征及其對地下水的影響[J]. 環(huán)境科學(xué)學(xué)報(bào),2008,28(9):1832-1838. Li Xueping, Sun Yan, Shi Xiaojun. Characteristics of phosphorus leaching and its impact on ground water in purple paddy soil[J]. Acta Scientiae Circumstantiae, 2008, 28(9): 1832 -1838. (in Chinese with English abstract)

    [26] 劉娟,張淑香,寧東衛(wèi),等. 3種耕作土壤磷隨地表徑流流失的特征及影響因素[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2019,35(10):1346-1352. Liu Juan, Zhang Shuxiang, Ning Dongwei, et al. Characteristics of phosphorus runoff losses and influencing facters in three cultivated soils[J]. Journal of Ecology and Rural Environment, 2019, 35(10): 1346-1352. (in Chinese with English abstract)

    [27] 秦華,李曄,李波,等. 人工模擬降雨條件下石灰土養(yǎng)分流失規(guī)律[J]. 水土保持學(xué)報(bào),2016,30(1):1-4. Qin Hua, Li Ye, Li Bo, et al. Nutrient loss of limeston soil under artificial simulated rainfall[J]. Journal of Ecology and Rural Environment, 2016, 30(1): 1-4. (in Chinese with English abstract)

    [28] 左繼超,鄭海金,奚同行,等. 自然降雨條件下紅壤坡地磷素隨徑流垂向分層輸出特征[J]. 環(huán)境科學(xué),2017,38(10):4178-4186. Zuo Jichao, Zheng Haijin, Xi Tonghang, et al. Characteristics of phosphorus output through runoff on a red soil slope under natural rainfall conditions[J]. Environmental Science, 2017, 38(10): 4178-4186. (in Chinese with English abstract)

    [29] 項(xiàng)大力,楊學(xué)云,孫本華,等. 灌溉水平對塿土磷素淋失的影響[J]. 植物營養(yǎng)與肥料學(xué)報(bào),2010,16(1):112-117. Xiang Dali, Yang Xueyun, Sun Benhua, et al. Impacts of irrigation regimes on phosphorus leaching in manural loessial soil[J]. Plant Nutrition and Fertilizer Science, 2010, 16(1): 112-117. (in Chinese with English abstract)

    [30] 陸欣欣,岳玉波,趙崢,等. 不同施肥處理稻田系統(tǒng)磷素輸移特征研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2014,22(4):394-400. Lu Xinxin, Yue Yubo, Zhao Zheng, et a1. Phosphorus loss and migration characteristics in paddy fields under different fertilization treatments[J]. Chinese Journal of Eco-Agriculture, 2014, 22(4): 394-400. (in Chinese with English abstract)

    [31] 鐘曉英,趙小蓉,鮑華軍,等. 我國23個土壤磷素淋失風(fēng)險評估Ⅰ. 淋失臨界值[J]. 生態(tài)學(xué)報(bào),2004,24(10):2275-2280. Zhong Xiaoying, Zhao Xiaorong, Bao Huajun, et al. The evaluation of phosphorus leaching risk of 23 Chinese soilsⅠ. Leaching criterion[J]. Acta Ecologica Sinica, 2004, 24(10): 2275-2280. (in Chinese with English abstract)

    [32] 聶敏,肖和艾,廖敦秀,等. 亞熱帶可變電荷土壤磷素淋失臨界點(diǎn)及其與土壤特性的關(guān)系[J]. 環(huán)境科學(xué)學(xué)報(bào),2013,33(2):579-586. Nie Min, Xiao Heai, Liao Dunxiu, et al. Phosphorus leaching change point of subtropical variable-charge soils and its relations with soil properties[J]. Acta Scientiae Circumstantiae, 2013, 33(2): 579-586. (in Chinese with English abstract)

    [33] 趙小蓉,鐘曉英,李貴桐,等. 我國23個土壤磷素淋失風(fēng)險評估Ⅱ. 淋失臨界值與土壤理化性質(zhì)和磷吸附特性的關(guān)系[J]. 生態(tài)學(xué)報(bào),2005,26(9):3011-3017. Zhao Xiaorong, Zhong Xiaoying, Li Guitong, et al. The evaluation of phosphorus leaching risk of 23 Chinese soilsⅡ. The relationships between soil properties, P absorption characteristics and the leaching criterion[J]. Acta Ecologica Sinica, 2005, 26(9): 3011-3017. (in Chinese with English abstract)

    [34] 張福鎖,崔振嶺,陳新平. 最佳養(yǎng)分管理技術(shù)單列[M]. 北京:中國農(nóng)業(yè)大學(xué)出版社,2010:5-7.

    Accumulative loss risk of phosphorus and its environmental threshold in fluvo-aquic soil

    Liu Juan1,3, Zhang Naiming2,3※, Zhang Shuxiang4

    (1.,,650201,; 2.,,650201,; 3.,650201,; 4.,,100081,)

    Fluvo-aquic soil is a typical cultivated soil with the widely distribution and high intensity of fertilization in China. However, the accumulation and loss of phosphorus in fluvo-aquic soil cannot be ignored as the ecological changes. In this study, the typical fluvo-aquic soils were collected as the test soils with different phosphorus levels (high, medium, and low) in Henan province, China, in order to explore the characteristics of phosphorus loss from fluvo-aquic soil induced by surface runoff and leaching. Simulation experiments of an artificial rainfall and a soil column were conducted to measure the content of Olsen-P and soluble active phosphorus (CaCl2-P) in soils, and the concentrations of different forms of phosphorus in runoff or leaching. In a split-line model, the relationship between soil Olsen-P content and Total Phosphorus (TP) in solution was used to evaluate P risk from black soil in surface runoff and leaching. The results showed that: 1) The contents of Olsen-P and CaCl2-P in fluvo-aquic soil profiles were in the order of high, medium and low phosphorus level, whereas, the soil phosphorus sorption index value were in the order of low, medium and high phosphorus level. In the profile distribution of phosphorus, the content of Olsen-P and CaCl2-P in fluvo-aquic soil with low and medium phosphorus level decreased with increasing soil profile, while the content of Olsen-P and CaCl2-P in fluvo-aquic soil with high phosphorus level was the highest in the 20-40 cm soil layer. 2) The concentration and losses of TP, Total Dissolved Phosphorus (TDP) and particulate phosphorus (PP) in the runoff from fluvo-aquic soils with different phosphorus levels were in the order of high > medium > low phosphorus level, and the PP dominated in fluvo-aquic soil runoff. 3) The concentration and loss of various forms of phosphorus in the leachate of low- and medium-phosphorus level from fluvo-aquic soil decreased with the increasing soil depth. However, in the leachate of high phosphorus level, the concentration and loss of phosphorus in the 20-40 cm soil profile were significantly higher than those of in other soil profiles, where the phosphorus concentration in the whole profile increased first and then decreased with the increasing soil profile, whereas, the profile of high phosphorus level presented the distribution pattern that upper and lower was low, and middle was high. The leaching loss of fluvo-aquic soil was dominated by TDP. The high- and low-phosphorus fluvo-aquic soil was dominated by Dissolved Organic Phosphorus (DOP), while the medium-phosphorus fluvo-aquic soil was dominated by Molybdate Reactive Phosphorus (MRP). 4) The relationship between CaCl2-P and Olsen-P were characterized using split-line models, identifying the thresholds at 24.65 Olsen-P mg/kg for the fluvo-aquic soil. It was also pointed out that the TP concentration in runoff and leachate was positively correlated with the CaCl2-P content in soil, and that CaCl2-P can be used to predict the phosphorus concentration of water discharged from soil. The outcomes derived from these experimental conditions regarding the determination of P loss in runoff by the artificially simulated rainfall experiment. The P loss in leaching by the soil column method and P environmental thresholds also need further verification in the field. This finding can help to assess phosphorus loss from fluvo-aquic soil under different phosphorus levels, and further to provide a scientific basis for water environment protection and soil phosphorus management in fluvo-aquic soil areas.

    phosphorus; soils; runoff; fluvo-aquic soil; environmental threshold; loss risk

    劉娟,張乃明,張淑香. 潮土磷素累積流失風(fēng)險及環(huán)境閾值[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(20):8-16.doi:10.11975/j.issn.1002-6819.2020.20.002 http://www.tcsae.org

    Liu Juan, Zhang Naiming, Zhang Shuxiang. Accumulative loss risk of phosphorus and its environmental threshold in fluvo-aquic soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(20): 8-16. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.20.002 http://www.tcsae.org

    2020-06-01

    2020-10-10

    國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFD0201208-2);云南省對外科技合作計(jì)劃院士工作站(2015IC022)

    劉娟,博士生,主要從事農(nóng)用化學(xué)物質(zhì)與環(huán)境方面的研究。Email:15587214232@163.com

    張乃明,博士,教授,主要從事土壤培肥與污染修復(fù)方面的研究。

    10.11975/j.issn.1002-6819.2020.20.002

    S158.5;X522

    A

    1002-6819(2020)-20-0008-09

    猜你喜歡
    潮土磷素徑流
    磷素添加對土壤水分一維垂直入滲特性的影響
    不同土壤類型小麥測土配方施肥對當(dāng)季土壤供氮的影響
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    探秘“大徑流”
    攻克“大徑流”
    江埡水庫降雨徑流相關(guān)圖的建立
    作物高效利用土壤磷素的研究進(jìn)展
    土壤與作物(2013年4期)2013-03-11 18:47:35
    長期施鉀和秸稈還田對河北潮土區(qū)作物產(chǎn)量和土壤鉀素狀況的影響
    磷素營養(yǎng)對大豆磷素吸收及產(chǎn)量的影響
    典型潮土N2O排放的DNDC模型田間驗(yàn)證研究
    a级片在线免费高清观看视频| 国产精品久久久久久精品电影小说| 免费在线观看成人毛片| 国产又色又爽无遮挡免| 男的添女的下面高潮视频| 成人国产麻豆网| 亚洲精品国产色婷婷电影| 2021少妇久久久久久久久久久| 久久久久久久国产电影| 亚洲经典国产精华液单| 草草在线视频免费看| 日本av手机在线免费观看| a 毛片基地| 国产高清三级在线| 午夜福利网站1000一区二区三区| xxx大片免费视频| 精品卡一卡二卡四卡免费| a级毛色黄片| 久热这里只有精品99| 在线看a的网站| 97超视频在线观看视频| 一级毛片我不卡| 免费黄网站久久成人精品| 国产在线视频一区二区| 99热国产这里只有精品6| 午夜免费鲁丝| 成年av动漫网址| 国产一区二区三区综合在线观看 | 少妇的逼水好多| 在线看a的网站| 日韩中文字幕视频在线看片| a 毛片基地| 亚洲丝袜综合中文字幕| 汤姆久久久久久久影院中文字幕| 一级毛片久久久久久久久女| 久久精品久久久久久久性| 国产伦精品一区二区三区四那| 我的老师免费观看完整版| 国产一区亚洲一区在线观看| 欧美日韩av久久| 男人添女人高潮全过程视频| 日韩一本色道免费dvd| 高清av免费在线| 女性被躁到高潮视频| tube8黄色片| 国产熟女欧美一区二区| 国产一区亚洲一区在线观看| 精品一区二区三区视频在线| 97超碰精品成人国产| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 建设人人有责人人尽责人人享有的| 午夜视频国产福利| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 丝袜脚勾引网站| 99久久精品一区二区三区| 男的添女的下面高潮视频| 一级a做视频免费观看| 亚洲国产精品成人久久小说| 久久毛片免费看一区二区三区| 亚洲激情五月婷婷啪啪| 日韩成人av中文字幕在线观看| 亚洲av综合色区一区| 亚洲,欧美,日韩| 亚洲美女黄色视频免费看| 久久综合国产亚洲精品| 超碰97精品在线观看| 高清视频免费观看一区二区| 国产欧美日韩一区二区三区在线 | 欧美最新免费一区二区三区| av免费在线看不卡| 欧美精品亚洲一区二区| kizo精华| 免费观看a级毛片全部| a级毛片免费高清观看在线播放| 国产精品99久久99久久久不卡 | 亚洲激情五月婷婷啪啪| 免费观看性生交大片5| 国产精品伦人一区二区| 久久久久久久久久人人人人人人| 自拍欧美九色日韩亚洲蝌蚪91 | 免费观看a级毛片全部| 亚洲,一卡二卡三卡| 亚洲国产成人一精品久久久| 成年女人在线观看亚洲视频| 久久热精品热| 国产精品国产三级国产av玫瑰| 一级二级三级毛片免费看| 久久精品久久久久久久性| 国产一区二区三区综合在线观看 | 欧美精品一区二区大全| 精品一区在线观看国产| 欧美丝袜亚洲另类| 人体艺术视频欧美日本| 大码成人一级视频| 麻豆精品久久久久久蜜桃| 午夜福利,免费看| 国产 一区精品| 91在线精品国自产拍蜜月| 中文字幕制服av| 伦理电影大哥的女人| 国产精品秋霞免费鲁丝片| 精品久久久精品久久久| 秋霞伦理黄片| 日韩伦理黄色片| 视频中文字幕在线观看| 亚洲不卡免费看| 高清午夜精品一区二区三区| 亚洲精品国产成人久久av| 妹子高潮喷水视频| 九九在线视频观看精品| 精品一区二区免费观看| 丰满乱子伦码专区| 国产精品成人在线| 久久女婷五月综合色啪小说| 亚洲美女黄色视频免费看| 在线观看美女被高潮喷水网站| 春色校园在线视频观看| 蜜桃久久精品国产亚洲av| 少妇的逼水好多| 精品国产一区二区三区久久久樱花| 久久久久久久久久成人| 两个人免费观看高清视频 | av免费在线看不卡| 校园人妻丝袜中文字幕| 久久久久国产精品人妻一区二区| 纯流量卡能插随身wifi吗| 人妻 亚洲 视频| 久久99一区二区三区| 欧美日韩视频精品一区| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品国产成人久久av| 91久久精品国产一区二区成人| 国产黄色视频一区二区在线观看| 成年av动漫网址| 午夜福利,免费看| 日韩av不卡免费在线播放| 亚洲av男天堂| 亚洲精品久久午夜乱码| 777米奇影视久久| 深夜a级毛片| 国产亚洲欧美精品永久| 国产免费一区二区三区四区乱码| 国产精品99久久99久久久不卡 | 国产男人的电影天堂91| 亚洲国产精品一区三区| 99热这里只有是精品在线观看| 91精品一卡2卡3卡4卡| 国产精品无大码| 亚洲情色 制服丝袜| 在线观看www视频免费| 国产毛片在线视频| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 亚洲国产色片| 中文字幕久久专区| 国产成人精品久久久久久| 3wmmmm亚洲av在线观看| 人人妻人人添人人爽欧美一区卜| 精品99又大又爽又粗少妇毛片| 性高湖久久久久久久久免费观看| 草草在线视频免费看| 欧美日韩国产mv在线观看视频| 久热久热在线精品观看| 欧美老熟妇乱子伦牲交| 在线观看免费日韩欧美大片 | 中文字幕av电影在线播放| 国产无遮挡羞羞视频在线观看| 日韩视频在线欧美| 一级,二级,三级黄色视频| 国产精品女同一区二区软件| 天堂俺去俺来也www色官网| 国产成人午夜福利电影在线观看| 日韩av在线免费看完整版不卡| 青青草视频在线视频观看| 亚洲精品中文字幕在线视频 | 哪个播放器可以免费观看大片| 久久青草综合色| 日韩一区二区三区影片| 赤兔流量卡办理| 天堂中文最新版在线下载| 又粗又硬又长又爽又黄的视频| 男人爽女人下面视频在线观看| 少妇精品久久久久久久| 免费在线观看成人毛片| 日韩 亚洲 欧美在线| 久久精品国产a三级三级三级| 国产毛片在线视频| 免费人妻精品一区二区三区视频| 亚洲国产欧美在线一区| 美女内射精品一级片tv| 欧美成人午夜免费资源| 色网站视频免费| 777米奇影视久久| 女性被躁到高潮视频| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 色婷婷av一区二区三区视频| 天美传媒精品一区二区| 国产精品一二三区在线看| 男女无遮挡免费网站观看| 一个人看视频在线观看www免费| 国产亚洲5aaaaa淫片| 日韩强制内射视频| 秋霞在线观看毛片| 少妇被粗大的猛进出69影院 | 亚洲国产精品999| 国产精品秋霞免费鲁丝片| 男人和女人高潮做爰伦理| 免费看光身美女| av又黄又爽大尺度在线免费看| 国产亚洲一区二区精品| 夜夜看夜夜爽夜夜摸| 久久综合国产亚洲精品| 少妇精品久久久久久久| 人人妻人人澡人人看| 国产成人aa在线观看| 欧美最新免费一区二区三区| 天美传媒精品一区二区| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 在线免费观看不下载黄p国产| 黄色日韩在线| 伦理电影免费视频| 女的被弄到高潮叫床怎么办| 久久久久久久久久成人| 在线免费观看不下载黄p国产| 精品少妇内射三级| 另类精品久久| 色5月婷婷丁香| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 香蕉精品网在线| 国产69精品久久久久777片| 久久这里有精品视频免费| 日本wwww免费看| 国产乱来视频区| 肉色欧美久久久久久久蜜桃| 热re99久久国产66热| 欧美高清成人免费视频www| 我的女老师完整版在线观看| 狠狠精品人妻久久久久久综合| 午夜福利网站1000一区二区三区| 午夜激情久久久久久久| 日韩中文字幕视频在线看片| 欧美区成人在线视频| av免费观看日本| 久久午夜综合久久蜜桃| 日本欧美国产在线视频| 亚洲不卡免费看| 国产成人精品福利久久| 菩萨蛮人人尽说江南好唐韦庄| 春色校园在线视频观看| 日韩视频在线欧美| 日韩一本色道免费dvd| 中国国产av一级| 99热这里只有是精品50| 七月丁香在线播放| 自线自在国产av| 汤姆久久久久久久影院中文字幕| 免费观看的影片在线观看| 国产在线视频一区二区| 亚洲情色 制服丝袜| 成人综合一区亚洲| 黑人高潮一二区| 欧美日韩亚洲高清精品| 黄片无遮挡物在线观看| 秋霞在线观看毛片| 日韩av不卡免费在线播放| 一本大道久久a久久精品| 岛国毛片在线播放| 国产综合精华液| 男女无遮挡免费网站观看| 精品一区二区免费观看| 九九爱精品视频在线观看| 一区二区三区精品91| 亚洲av免费高清在线观看| 成人亚洲精品一区在线观看| av在线播放精品| 97超碰精品成人国产| 国产亚洲91精品色在线| 又爽又黄a免费视频| 欧美日韩视频精品一区| 热re99久久精品国产66热6| 中文天堂在线官网| 老司机影院成人| 亚洲经典国产精华液单| 91精品伊人久久大香线蕉| 亚洲精品一二三| 国产成人精品婷婷| 精品久久久噜噜| www.色视频.com| 特大巨黑吊av在线直播| 成人黄色视频免费在线看| 亚洲精品第二区| 三级国产精品欧美在线观看| 中文字幕精品免费在线观看视频 | 纯流量卡能插随身wifi吗| 日产精品乱码卡一卡2卡三| 五月天丁香电影| 免费av不卡在线播放| 日本av手机在线免费观看| 欧美高清成人免费视频www| 九色成人免费人妻av| 涩涩av久久男人的天堂| 极品人妻少妇av视频| 女的被弄到高潮叫床怎么办| 日韩中字成人| 18禁动态无遮挡网站| 国产精品福利在线免费观看| 啦啦啦啦在线视频资源| 午夜av观看不卡| 99热6这里只有精品| 永久免费av网站大全| 春色校园在线视频观看| 成人亚洲欧美一区二区av| 91久久精品国产一区二区成人| 一个人看视频在线观看www免费| 中国国产av一级| 80岁老熟妇乱子伦牲交| 久久人人爽人人爽人人片va| 国产淫语在线视频| 亚洲欧洲日产国产| 久久狼人影院| 有码 亚洲区| 免费少妇av软件| 一本色道久久久久久精品综合| 亚洲国产精品国产精品| 日日撸夜夜添| 亚洲丝袜综合中文字幕| 日日啪夜夜撸| 91在线精品国自产拍蜜月| av黄色大香蕉| 亚洲av中文av极速乱| 午夜福利视频精品| 在线观看免费视频网站a站| 日韩不卡一区二区三区视频在线| 三上悠亚av全集在线观看 | 中文字幕免费在线视频6| 亚洲天堂av无毛| 婷婷色av中文字幕| 最黄视频免费看| 美女视频免费永久观看网站| 国产一区二区在线观看日韩| 亚洲不卡免费看| 亚洲欧洲国产日韩| 日日啪夜夜爽| 中国三级夫妇交换| 国产日韩欧美在线精品| 国产精品免费大片| 久久毛片免费看一区二区三区| 多毛熟女@视频| 亚洲丝袜综合中文字幕| 精品一区在线观看国产| 亚洲国产精品999| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 狂野欧美激情性xxxx在线观看| 两个人免费观看高清视频 | 少妇高潮的动态图| 亚洲高清免费不卡视频| 一级毛片黄色毛片免费观看视频| 一本一本综合久久| 天天躁夜夜躁狠狠久久av| 精品一区二区三区视频在线| 制服丝袜香蕉在线| 国产美女午夜福利| 国产亚洲5aaaaa淫片| 99热这里只有是精品在线观看| 只有这里有精品99| 欧美日韩亚洲高清精品| 国产一区二区在线观看av| 精品亚洲乱码少妇综合久久| 色婷婷av一区二区三区视频| 久久久久久久亚洲中文字幕| 久久人人爽人人爽人人片va| 人人妻人人看人人澡| 成人影院久久| 在线精品无人区一区二区三| 亚洲欧美日韩卡通动漫| 成人亚洲欧美一区二区av| 一本久久精品| 两个人的视频大全免费| 99热全是精品| 美女主播在线视频| 十八禁高潮呻吟视频 | 18+在线观看网站| 成人18禁高潮啪啪吃奶动态图 | 成人亚洲欧美一区二区av| 少妇猛男粗大的猛烈进出视频| 精品久久久噜噜| 在线免费观看不下载黄p国产| 人妻 亚洲 视频| 深夜a级毛片| 黄色视频在线播放观看不卡| 国产精品99久久久久久久久| 少妇的逼好多水| 极品少妇高潮喷水抽搐| 精华霜和精华液先用哪个| 亚洲精品国产av蜜桃| 人人妻人人澡人人看| 国产乱人偷精品视频| 亚洲精品乱久久久久久| a级一级毛片免费在线观看| 国产视频首页在线观看| 免费看光身美女| 少妇猛男粗大的猛烈进出视频| 国产精品一区二区在线观看99| 欧美日韩国产mv在线观看视频| 免费观看在线日韩| 欧美xxxx性猛交bbbb| 在线观看免费日韩欧美大片 | 国产免费又黄又爽又色| 成人特级av手机在线观看| 色5月婷婷丁香| 久久97久久精品| 亚洲经典国产精华液单| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 国产成人aa在线观看| 免费人妻精品一区二区三区视频| 久久久久久久大尺度免费视频| 搡女人真爽免费视频火全软件| 国产日韩欧美亚洲二区| 色哟哟·www| 日韩在线高清观看一区二区三区| 免费观看性生交大片5| 欧美最新免费一区二区三区| 成年人午夜在线观看视频| 亚洲精品视频女| 另类精品久久| 色婷婷av一区二区三区视频| 久久狼人影院| 免费大片18禁| 国产高清不卡午夜福利| 欧美日韩av久久| 国产精品偷伦视频观看了| 久久久久久久国产电影| av又黄又爽大尺度在线免费看| 两个人免费观看高清视频 | 一级毛片黄色毛片免费观看视频| 熟女人妻精品中文字幕| 亚洲美女黄色视频免费看| 亚洲美女搞黄在线观看| 国产成人精品婷婷| 午夜福利,免费看| 少妇猛男粗大的猛烈进出视频| 日本免费在线观看一区| 秋霞在线观看毛片| 少妇人妻 视频| 亚洲av成人精品一区久久| 成人特级av手机在线观看| 十八禁网站网址无遮挡 | 日韩人妻高清精品专区| 极品教师在线视频| 我的女老师完整版在线观看| 少妇的逼水好多| 国产高清有码在线观看视频| 亚洲电影在线观看av| 黄色视频在线播放观看不卡| 日本色播在线视频| tube8黄色片| 伦理电影免费视频| 看十八女毛片水多多多| 国产精品熟女久久久久浪| 最新的欧美精品一区二区| 午夜福利影视在线免费观看| 两个人的视频大全免费| 欧美3d第一页| 成人国产av品久久久| 国产精品一区二区在线观看99| 热99国产精品久久久久久7| freevideosex欧美| 九九爱精品视频在线观看| 国产成人免费观看mmmm| 18+在线观看网站| 91精品国产国语对白视频| 欧美xxⅹ黑人| 另类精品久久| 最近最新中文字幕免费大全7| 国产伦精品一区二区三区四那| 嫩草影院新地址| 国产精品久久久久久久电影| 在线免费观看不下载黄p国产| 日韩伦理黄色片| av在线老鸭窝| 18禁在线播放成人免费| 插逼视频在线观看| 一级黄片播放器| 男女边摸边吃奶| 97精品久久久久久久久久精品| 国产精品一二三区在线看| 在线看a的网站| 亚洲精品日本国产第一区| 永久免费av网站大全| 嫩草影院入口| 夫妻性生交免费视频一级片| 美女cb高潮喷水在线观看| 欧美xxⅹ黑人| 久久久久网色| 91久久精品电影网| 99热6这里只有精品| 国产午夜精品久久久久久一区二区三区| 妹子高潮喷水视频| 3wmmmm亚洲av在线观看| 乱人伦中国视频| 十分钟在线观看高清视频www | 亚洲国产日韩一区二区| 日日摸夜夜添夜夜添av毛片| 男女边吃奶边做爰视频| 久热这里只有精品99| 国产男人的电影天堂91| 日本欧美视频一区| 日韩中文字幕视频在线看片| 久久狼人影院| 91久久精品国产一区二区三区| 国产精品国产av在线观看| 又粗又硬又长又爽又黄的视频| 国产在线男女| 欧美另类一区| 国产伦理片在线播放av一区| 精品少妇久久久久久888优播| 少妇高潮的动态图| av在线播放精品| 欧美人与善性xxx| 亚洲欧美日韩卡通动漫| 人人妻人人澡人人爽人人夜夜| 国产在视频线精品| 国产白丝娇喘喷水9色精品| 最后的刺客免费高清国语| 99视频精品全部免费 在线| 一本大道久久a久久精品| 纵有疾风起免费观看全集完整版| 简卡轻食公司| 免费高清在线观看视频在线观看| 国产69精品久久久久777片| 免费人妻精品一区二区三区视频| 免费看av在线观看网站| 国产熟女午夜一区二区三区 | 男人添女人高潮全过程视频| 黑人猛操日本美女一级片| 国产毛片在线视频| 日日爽夜夜爽网站| 一级毛片黄色毛片免费观看视频| 男人爽女人下面视频在线观看| freevideosex欧美| 18禁动态无遮挡网站| 一级毛片 在线播放| 日韩电影二区| 久久久久人妻精品一区果冻| 精品少妇黑人巨大在线播放| 在线观看免费高清a一片| 女人久久www免费人成看片| 国产免费视频播放在线视频| 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 这个男人来自地球电影免费观看 | 中国国产av一级| 国产精品人妻久久久久久| 久久人人爽av亚洲精品天堂| 亚洲在久久综合| av卡一久久| 中文字幕人妻丝袜制服| 乱系列少妇在线播放| 免费人妻精品一区二区三区视频| 国产精品久久久久久精品古装| 最黄视频免费看| av黄色大香蕉| 日韩中文字幕视频在线看片| 久热这里只有精品99| 在线观看美女被高潮喷水网站| 性高湖久久久久久久久免费观看| 国语对白做爰xxxⅹ性视频网站| 午夜福利网站1000一区二区三区| 国产精品99久久99久久久不卡 | 夜夜看夜夜爽夜夜摸| 水蜜桃什么品种好| 欧美日韩在线观看h| 亚洲精品乱码久久久v下载方式| 国产成人91sexporn| 久久婷婷青草| av国产久精品久网站免费入址| 国产精品福利在线免费观看| 成人综合一区亚洲| 99久久综合免费| 国产精品成人在线| 久久精品熟女亚洲av麻豆精品| 美女内射精品一级片tv| 亚洲欧美精品专区久久| 亚洲三级黄色毛片| 精品少妇内射三级| 边亲边吃奶的免费视频| 日韩一区二区三区影片| 一级毛片aaaaaa免费看小| 亚洲伊人久久精品综合| 极品教师在线视频| 亚洲av成人精品一二三区| 国产一区有黄有色的免费视频| 亚洲精品国产成人久久av| 人妻制服诱惑在线中文字幕| 99九九在线精品视频 | 少妇裸体淫交视频免费看高清| 伦理电影大哥的女人| 大香蕉97超碰在线| 婷婷色麻豆天堂久久| 在线亚洲精品国产二区图片欧美 | 精品久久国产蜜桃| 精品少妇黑人巨大在线播放| 成人无遮挡网站| 日本黄大片高清| 国产日韩欧美视频二区| 搡女人真爽免费视频火全软件| 中文字幕久久专区|