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

    豬場肥水施用對玉米-小麥農田氨排放、氮素利用與表觀平衡的影響*

    2019-10-15 12:51:20張西群劉澤龍趙占輕張玉銘孫宏勇
    關鍵詞:肥水表觀氮素

    李 碩, 王 選, 張西群, 劉澤龍, 趙 浩, 趙占輕,張玉銘, 孫宏勇, 馬 林**

    豬場肥水施用對玉米-小麥農田氨排放、氮素利用與表觀平衡的影響*

    李 碩1, 王 選1, 張西群2, 劉澤龍1, 趙 浩1, 趙占輕1,張玉銘1, 孫宏勇1, 馬 林1**

    (1. 中國科學院遺傳與發(fā)育生物學研究所農業(yè)資源研究中心/中國科學院農業(yè)水資源重點實驗室/河北省土壤生態(tài)學重點實驗室 石家莊 050022; 2. 河北省農業(yè)機械化研究所有限公司 石家莊 050022)

    規(guī)?;i養(yǎng)殖廢棄物已成為當前重要污染來源, 為有效解決豬場廢水所引發(fā)面源污染問題, 有必要開展將其替代礦物氮肥(作為肥水)施用于農田的探索。以華北平原高度集約化玉米-小麥一年兩熟輪作體系為對象, 通過田間小區(qū)試驗, 定量研究豬場肥水施用對作物產量、氮素吸收、氮素利用效率、土壤礦質氮累積、氨揮發(fā)損失及輪作體系氮素表觀平衡的影響。試驗包括7個處理: 不施肥對照(CK)、尿素表施(CK1)、尿素注射施用(CK2)、豬場肥水替代25%尿素氮表施(25%WB)、豬場肥水替代50%尿素氮表施(50%WB)、豬場肥水替代25%尿素氮注射施用(25%WI)和豬場肥水替代50%尿素氮注射施用(50%WI)。豬場肥水作為基肥施用。結果表明, 與CK相比, 施用尿素和豬場肥水均可顯著提高玉米、小麥產量和籽粒氮吸收量, 其中25%WI最高, 50%WI次之。與尿素表施相比, 尿素注射施用、豬場肥水表施和注射施用均可明顯提高氮肥農學效率、偏生產力和表觀利用率, 且肥水注射施用最高, 肥水表施次之, 而25%WI和50%WI之間無顯著差異。與不施肥處理相比, 施用尿素和豬場肥水0~100 cm土體礦質氮殘留量顯著增加50.8%~87.9%, 其中50%WB、25%WI和50%WI無顯著差異。與尿素表施相比, 尿素注射施用、肥水表施和注射施用均可顯著降低玉米和小麥基肥期土壤氨損失總量, 降幅分別為26.5%~48.6%和11.4%~29.1%; 同時, 肥水表施和注射施用下輪作體統(tǒng)氮盈余顯著降低7.6%~16.0%, 其中25%WI降幅最高, 但與50%WI無顯著差異。綜合考慮作物產量、氮素利用和環(huán)境效應, 豬場肥水替代25%和50%尿素氮注射施用是該區(qū)玉米-小麥輪作農田豬場肥水最佳施用方式。

    玉米-小麥輪作農田;豬場肥水;注射施用;氨排放;氮素利用;氮素表觀平衡

    近年來, 我國華北平原規(guī)?;i養(yǎng)殖快速發(fā)展, 豬場廢水隨意排放不斷增加, 導致地下水硝態(tài)氮累積、氨揮發(fā)、溫室氣體排放和空氣惡臭等一系列環(huán)境問題, 嚴重制約著規(guī)?;i養(yǎng)殖的可持續(xù)發(fā)展[1-2]。豬場廢水含有豐富的有機質和作物生長所必需的氮、磷、鉀等養(yǎng)分物質, 因此, 將豬場廢水用于肥水就近施用于農田, 可有效改善土壤肥力質量, 已成為豬場廢水經濟有效的資源化利用途徑之一[3-5]。陶曉婷等[6]研究表明, 小麥(L.)越冬期施用豬場肥水60~120 m3?hm–2替代穗期施氮30~60 kg?hm–2, 可滿足小麥氮素營養(yǎng)需求。也有研究表明, 以豬場/牛場肥水替代化肥50%~75%N,可獲得等氮化肥處理相當?shù)淖魑锂a量[6-7]。然而, 在傳統(tǒng)的大水漫灌施用方式下, 豬場肥水中大量氮素會以硝態(tài)氮(NO3–-N)形式淋溶到深層土壤, 或以氨(NH3)形式排放于大氣中, 這將嚴重導致豬場肥水氮素利用效率的降低[8-9]。豬場肥水注射施用能夠降低NO3–-N淋溶和NH3排放等形態(tài)的氮損失, 進而有望提升肥水施用的氮素利用效率及其增產效應[2,10]。國內外針對于豬場肥水不同施用方式的對比研究已有很多, 然而, 這些研究主要集中于小麥、玉米(L.)、水稻(L.)或牧草等單作體系下的氨揮發(fā)、溫室氣體排放和土壤肥力質量變化[11-14], 而針對作物氮素利用效率的研究并不多[15-17]。當前, 主導華北平原的玉米-小麥一年兩熟高度集約化輪作生產體系在產量不斷提高的同時, 對耕地利用強度高, 水肥資源投入量大, 進而導致該體系肥料效益低且引發(fā)環(huán)境污染等問題[18-22]。因此, 如何在降低化肥投入的同時實現(xiàn)作物增產和養(yǎng)分高效利用已成為人們關注的熱點。以往研究大多從合理施肥、優(yōu)化作物栽培模式和養(yǎng)分管理等角度對作物生長和肥料利用率進行探討, 而該體系針對豬場肥水不同施用方式對作物產量、氮素利用和土壤氮素平衡的研究卻鮮有報道[15-17]?;诖? 本文以華北平原秸稈還田條件下玉米-小麥輪作體系為研究對象, 探討了豬場肥水不同施用方式(表施、注射)和氮素替代率(25%和50%)對作物產量、肥水中氮素利用和土壤氮素表觀平衡的影響, 以期為該地區(qū)合理進行豬場肥水施用, 提高肥水中氮素利用效率, 降低養(yǎng)殖肥水施用過程中氮素的損失提供理論依據。

    1 材料與方法

    1.1 試驗地概況

    試驗于2017年6月—2018年6月在河北省南皮縣張旗屯(38°10′N, 116°52′E)進行。該試驗區(qū)屬典型暖溫帶半濕潤季風氣候, 年平均氣溫13.4 ℃, 日照時數(shù)2 318 h, 年均降水量520 mm, 其中60%~70%的降水集中于8月。供試土壤為當?shù)氐湫偷妮p壤潮土。試驗開始時耕層(0~20 cm)土壤有機碳12.25 g?kg–1, 硝態(tài)氮13.20 mg?kg–1, 銨態(tài)氮2.64 mg?kg–1, 速效磷16.43 mg?kg–1, 速效鉀96.55 mg?kg–1。0~100 cm土體各層次(以20 cm計)的容重分別為1.30 g?cm–3、1.45 g?cm–3、1.50 g?cm–3、1.50 g?cm–3和1.50 g?cm–3。0~180 cm土體各層次(以20 cm計)礦質氮(NO3–-N和NH4+-N之和)含量分別為41.18 kg?hm–2、34.03 kg?hm–2、52.88 kg?hm–2、59.73 kg?hm–2、74.80 kg?hm–2、44.66 kg?hm–2、42.43 kg?hm–2、44.16 kg?hm–2和44.51 kg?hm–2。

    1.2 試驗設計

    試驗選取土壤肥力均一、灌溉便捷、區(qū)域代表性較強的田塊進行。實行玉米-小麥一年兩熟制, 玉米于6月上旬播種, 10月初收獲, 品種為‘正弘8號’, 行距50 cm, 密度63 000株?hm–2; 小麥于當年10月中旬播種, 翌年6月初收獲, 品種為‘小偃81’, 行距15 cm, 播種量450 kg?hm–2。兩茬作物秸稈全部還田。試驗包括7個處理(表1): 不施肥對照(CK)、尿素表施(CK1)、尿素注射施用(CK2)、豬場肥水替代25%尿素氮表施(25%WB)、豬場肥水替代50%尿素氮表施(50%WB)、豬場肥水替代25%尿素氮注射施用(25%WI)和豬場肥水替代50%尿素氮注射施用(50%WI)。其中, 注射施用采用課題組研制的畜禽肥水淺層注射施用機具(肥水注射機具)來開展。50%WB和50%WI處理豬場肥水施用量為154 m3?hm–2, 且其注射深度為5~10 cm; CK2處理50%尿素溶入154 m3?hm–2地下水中, 采用肥水注射機具施入5~10 cm土層; 其他處理以地下水補足154 m3?hm–2施用方式為表層噴灑, 以便降低施肥環(huán)節(jié)因水分投入而導致的差異。采用完全隨機區(qū)組設計, 重復3次, 小區(qū)面積72 m2(6 m×12 m)。

    表1 不同處理的玉米和小麥生長期氮素施用類型、施用量和施用方法

    玉米生長期氮、磷、鉀施用量分別為240 kg(N)?hm–2、180 kg(P2O5)?hm–2和120 kg(K2O)?hm–2, 小麥生長期氮、磷、鉀施用量分別為180 kg(N)?hm–2、120 kg(P2O5)?hm–2和90 kg(K2O)?hm–2, 其中磷、鉀肥以基肥形式一次性施用于作物播種前。氮肥分基肥和追肥兩次施用, 基追比為1∶1; 玉米追肥期為大喇叭口期, 小麥追肥期為拔節(jié)期。試驗中所用豬場肥水取自張旗屯溫氏集團養(yǎng)豬場調節(jié)池廢水, 試驗前10 d內取樣測定肥水理化性質, 其中玉米和小麥播種前肥水理化性質保持相對穩(wěn)定, 其含銨態(tài)氮545.05 mg?L–1、硝態(tài)氮46.05 mg?L–1、全氮779.33 mg?L–1、全磷39.25 mg?L–1、全鉀500.00 mg?L–1, 化學需氧量(COD)為3 610.00 mg?L–1。豬場肥水作為基肥一次性施用于農田, 養(yǎng)分不足時以礦物肥料補足。氮、磷和鉀礦物肥料分別為尿素、過磷酸鈣和硫酸鉀。玉米、小麥播種1周后進行灌溉(大水漫灌)。各試驗處理的耕作、播種、灌溉和病蟲害防治等田間操作均與當?shù)貍鹘y(tǒng)方式保持一致。

    1.3 氨揮發(fā)測定

    氨揮發(fā)田間原位測定采用間歇密閉抽氣法[23]。原理是利用真空泵減壓抽氣使罩子內土壤揮發(fā)出的氨(NH3)隨氣流通過裝有2%硼酸的洗氣瓶, 使其吸收于硼酸溶液中, 收集溶液用0.02 mol?L?1H2SO4滴定, 計算出吸收氮量。所用密閉室裝置由聚乙烯塑料(PVC)制成的密閉室(內直徑為25 cm, 高10 cm)、50 mL吸收瓶和高速真空泵組成。測定時, 在每個小區(qū)隨機選取一個樣點并將密閉室置于5 cm土壤中。通過預備試驗確定每天上午9:00—11:00和下午14:30—16:30時段測定值的平均值可基本代表全天揮發(fā)速率平均值。施肥后每天測定直至各處理與空白之間的氨揮發(fā)速率(kg?hm–2?d–1)均無明顯差異為止。通過每天氨揮發(fā)速率計算測定時期內的氨揮發(fā)累積損失(kg?hm–2)。

    1.4 樣品采集與氮素測定

    1.4.1 土壤樣品采集與測定

    試驗開始前在保護行內挖0~100 cm土壤剖面, 采用環(huán)刀法測定每20 cm土層的土壤容重, 并在試驗小區(qū)內部按“S”形路線采集混合土樣, 測定基礎養(yǎng)分含量。在玉米和小麥收獲后, 每個小區(qū)隨機選取1個樣點, 采集0~180 cm土層土樣, 層間距20 cm, 測定土壤含水量和NO3–-N和NH4+-N含量。土壤含水量采用烘干法測定; 土壤NH4+-N 和NO3–-N含量測定當天將土樣混勻后過2 mm篩, 然后用1 mol?L?1的KC1溶液浸提(水土比5∶1), 流動分析儀(OI, 美國)測定。

    1.4.2 植株樣品采集與測定

    在玉米和小麥收獲時, 每個小區(qū)隨機采集1個2 m2(2行, 2 m × 1 m)玉米樣品和1 m2(1 m × 1 m)小麥樣品。風干脫粒, 小麥樣品分籽粒和秸稈兩部分稱量其干重, 玉米樣品分籽粒、玉米穗軸和秸稈3部分稱量其干重, 然后以采樣面積折算生物量。取部分小麥和玉米樣品烘干、粉碎、混勻, 然后用濃H2SO4-H2O2消解, 采用凱氏定氮法測定植株樣品的氮含量, 用于作物吸氮量計算。

    1.5 數(shù)據計算與統(tǒng)計分析

    本研究在探討農田土壤氮平衡的過程中將種子、前茬作物秸稈、沉降、灌溉所攜入農田的氮均考慮在農田土壤氮輸入中[24]。根據氮素輸入-輸出平衡原理, 氮肥利用率和氮平衡參數(shù)的計算方法[18-19]如下:

    玉米-小麥輪作周期0~100 cm土層土壤氮素除氨揮發(fā)損失外的其他形式損失量(kg?hm–2)=施氮量+玉米播種前0~100 cm土層土壤起始礦質氮(NO3–-N和NH4+-N之和)+玉米-小麥輪作周期0~100 cm土層土壤氮素凈礦化量(kg?hm–2)+種子氮+沉降氮+灌溉氮-收獲作物攜出氮量-小麥收獲期0~100 cm土層土壤殘留礦質氮-氨揮發(fā)損失量 (1)

    氮素盈余(kg?hm?2)=玉米-小麥輪作周期0~100 cm土層土壤氮素表觀損失量+小麥收獲期0~100 cm土層土壤殘留礦質氮 (2)

    氮肥農學效率(kg?kg?1)=(施氮區(qū)產量-對照區(qū)產量)/氮肥投入量 (3)

    氮肥偏生產力(kg?kg?1)=施氮區(qū)籽粒產量/氮肥投入量 (4)

    氮肥表觀利用率(%)=(施氮區(qū)吸氮量-對照區(qū)吸氮量)/氮肥投入量×100 (5)

    土壤剖面各土層礦質氮累積量(Nmin, kg?hm–2)按照以下公式計算:

    Nmin=×b××0.1 (6)

    式中:為土層厚度(20 cm),b為土壤容重(g?cm?3),為土壤硝態(tài)氮和銨態(tài)氮含量(NO3–-N和NH4+-N, mg?kg?1), 0.1為單位換算系數(shù)。試驗前測定0~100 cm深度土壤容重(20 cm一層)分別為1.30 g?cm?3、1.45 g?cm?3、1.50 g?cm?3、1.50 g?cm?3和1.50 g?cm?3。

    試驗數(shù)據繪圖和統(tǒng)計分析分別采用Excel 2010和DPS 7.05 統(tǒng)計軟件進行, 多重比較采用Duncan法, 差異顯著性水平=0.05。

    2 結果與分析

    2.1 施用豬場肥水對作物籽粒產量的影響

    與CK相比, 施用氮肥和豬場肥水均顯著提高玉米籽粒產量(15.9%~54.1%)和小麥籽粒產量(3.9%~ 20.5%)。玉米和小麥籽粒產量均表現(xiàn)為CK

    表2 不同豬場肥水施用處理對玉米和小麥籽粒產量的影響

    相對增幅為各處理相對于CK1的增幅。同列不同小寫字母表示不同處理間差異顯著(<0.05)。The relative increase is the increase as to CK1. Different lowercase letters within the same column indicate significant differences among treatments at< 0.05.

    與CK相比, 施用氮肥和豬場肥水亦可顯著提高玉米和小麥地上部生物量產量, 其增幅分別為19.3%~55.1%和3.7%~17.5%(表3)。與CK1相比, 25%WI和50%WI處理均可顯著提高玉米地上部生物量產量, 其增幅分別為30.0%和29.9%; 25%WB、25%WI和50%WI處理均可顯著提高小麥地上部生物量產量, 其增幅分別為8.2%、12.5%和13.3%; 然而, 玉米和小麥地上部生物量產量在25%WB、50%WB、25%WI和50%WI之間無顯著差異。

    表3 不同豬場肥水施用處理對玉米和小麥地上部生物量的影響

    相對增幅: 各處理相對于CK1的增幅。同列不同小寫字母表示不同處理間差異顯著(<0.05)。The relative increase is the increase as to CK1. Different lowercase letters within the same column indicate significant differences among treatments at< 0.05.

    2.2 施用豬場肥水對作物氮素吸收和氮肥利用率的影響

    與CK相比, 施用氮肥和豬場肥水均可顯著提高作物籽粒和地上部氮素吸收量, 且玉米、小麥籽粒和地上部氮素吸收量的變化趨勢一致, 均表現(xiàn)為豬場肥水注射施用>豬場肥水表施>尿素注射施用>尿素表施>不施肥對照(表4)。其中, 25%WI處理玉米籽粒、地上部氮素吸收量和小麥籽粒、地上部氮素吸收量均最高, 與CK相比, 其增幅分別為72.6%、89.5%、25.3%和23.3%。

    在玉米和小麥生育期, 氮肥農學效率和氮肥偏生產力在各處理間變化趨勢與玉米和小麥籽粒氮素吸收量大致相同(表4)。在玉米生育期, 不同施肥處理氮肥農學效率為4.32~14.72 kg?kg–1, 其中以25%WI處理氮肥農學效率最高, 其次為50%WI和50%WB, 且3個處理間無顯著性差異; 同時, 25%WI和50%WI處理氮肥農學效率顯著高于25%WB、CK1和CK2。氮肥偏生產力為31.54~41.95 kg?kg–1, 與氮農學效率變化趨勢相似, 25%WI處理最高, 其次為50%WI、50%WB和25%WB, 且四者之間無顯著性差異; 同時25%WI和50%WI處理氮肥農學效率顯著高于CK1和CK2, 而CK1、CK2、25%WB和50%WB處理之間無顯著性差異。氮肥表觀利用率為19.02%~48.57%, 與氮肥農學效率和氮肥偏生產力變化趨勢一致, 25%WI處理最高, 其次為50%WI和50%WB, 且三者之間無顯著性差異; 同時, 25%WI、50%WI和50%WB的氮肥表觀利用率顯著高于其余各處理。

    小麥生育期, 不同施肥處理氮肥農學效率為1.22~6.37 kg?kg–1, 其中以25%WI處理氮肥農學效率最高, 其次為50%WI, 且兩個處理間無顯著性差異; 25%WI處理氮肥農學效率顯著高于其他處理。氮肥偏生產力為32.27~37.42 kg?kg–1, 同樣以25%WI最高, 50%WI次之, 且二者之間無顯著差異; 25%WI處理氮肥偏生產力顯著高于其他處理。氮肥表觀利用率為6.85%~22.36%, 與氮肥農學效率和氮肥偏生產力變化趨勢一致??傮w而言, 25%WI和50%WI可獲得最高的氮效率。

    表4 不同豬場肥水施用處理對玉米和小麥氮素吸收和氮肥利用效率的影響

    NAE: 氮肥農學效率; NPE: 氮肥偏生產力; ANR: 氮肥表觀利用率。同列不同小寫字母表示不同處理間差異顯著(<0.05)。NAE: N agronomic efficiency; NPE: N partial productivity; ANR: apparent N recovery rate. Different lowercase letters within the same column indicate significant differences among treatments at< 0.05.

    2.3 施用豬場肥水對土壤氨揮發(fā)的影響

    施用尿素和豬場肥水對玉米和小麥生育期土壤氨揮發(fā)速率的影響如圖1和圖2所示。在各測定時期, CK處理氨損失均保持較小的相對穩(wěn)定狀態(tài)。在玉米基肥期, 氨揮發(fā)峰值出現(xiàn)在施肥第1天, 且CK1處理下氨損失顯著高于其他各施肥處理, 其增幅為29.9%~64.4%, 而50%WI處理氨揮發(fā)速率最小, 其值為1.14 kg(N)?hm–2?d–1(圖1A)。各處理在第1天的氨揮發(fā)對測定期內氨損失總量貢獻率為26.1%~39.3%。隨后, 各處理氨揮發(fā)速率逐漸降低, 直至第4天達到相對穩(wěn)定狀態(tài)[~0.5 kg(N)?hm–2?d–1], 然而第5天再次出現(xiàn)一個峰值[0.55~0.92 kg(N)?hm–2?d–1], 隨后在第6天再次達到相對穩(wěn)定狀態(tài)[~0.3 kg(N)?hm–2?d–1]。與CK1相比, CK2、25%WB、50%WB、25%WI和50%WI在施肥6 d內氨損失總量分別顯著降低37.9%、34.8%、26.4%、38.6%和48.6%(表5)。

    在玉米追肥期, 氨揮發(fā)峰值同樣出現(xiàn)在第1天, 且該天氨揮發(fā)對測定期內氨損失總量貢獻率為58.2%~60.7%(圖1B)。隨后, 各處理氨揮發(fā)速率逐漸降低, 直至第6天達到相對穩(wěn)定狀態(tài)[~0.3 kg(N)?hm–2?d–1]。在第3天和第4天, 25%WI和50%WI的氨揮發(fā)速率顯著低于CK1, 其降幅分別為9.7%、13.5%(第3天)和13.5%、19.9%(第4天)。然而, 其他時間段各處理間均無顯著性差異。就玉米基肥和追肥氨損失總量而言, 25%WI和50%WI較CK1可顯著降低11.6%和16.6%, 且CK1、CK2、25%WB和50%WB之間無顯著差異(表5)。玉米基肥氨排放因子為1.46%~4.92%, 追肥氨排放因子為14.24%~14.65%。在玉米基肥和追肥氨排放總因子為7.85%~9.78%(表5)。

    在小麥基肥期, 氨揮發(fā)峰值出現(xiàn)在施肥第2天, 且CK1處理下氨損失高于其他各施肥處理, 增幅為5.9%~21.7%, 而50%WI處理氨揮發(fā)速率最小, 其值為2.76 kg(N)?hm–2?d–1(圖2A)。各處理在第2天的氨揮發(fā)對測定期內氨損失總量貢獻率為28.8%~33.4%。隨后, 各處理氨揮發(fā)速率逐漸降低, 直至第6天達到相對穩(wěn)定狀態(tài)[<0.5 kg(N)?hm–2?d–1]。與CK1相比, CK2、25%WB、50%WB、25%WI和50%WI在施肥后6 d天內氨損失總量分別顯著降低14.1%、20.6%、11.4%、16.2%和29.0%(表5)。

    在小麥追肥期, 氨揮發(fā)峰值同樣出現(xiàn)在第2天, 且該天氨揮發(fā)對測定期內氨損失總量貢獻率為32.8%~34.3%(圖2B)。隨后, 各處理氨揮發(fā)速率逐漸降低, 且各處理間無顯著差異, 直至第6天達到相對穩(wěn)定狀態(tài)[~0.5 kg(N)?hm–2?d–1]。就小麥基肥和追肥氨損失總量而言, 與CK1相比, CK2、25%WB、50%WB、25%WI和50%WI在施肥后6 d天內氨損失總量分別顯著降低9.7%、12.3%、7.4%、10.4%和16.7%(表6)。小麥基肥氨排放因子為8.08%~11.91%, 追肥氨排放因子為7.96%~8.33%, 基肥和追肥氨排放總因子為8.20~10.12%(表5)。

    圖1 不同豬場肥水施用處理下玉米基肥期(A)和追肥期(B)的氨排放速率

    施肥后同天不同小寫字母表示不同處理間差異顯著(<0.05)。Different lowercase letters at the same day after fertilization indicate significant differences among treatments at< 0.05.

    圖2 不同豬場肥水施用處理下小麥基肥期(A)和追肥期(B)的氨排放速率

    施肥后同天不同小寫字母表示不同施肥處理間差異顯著(<0.05)。Different lowercase letters at the same day after fertilization indicate significant differences among treatments at< 0.05.

    表5 玉米、小麥基肥和追肥期不同豬場肥水施用處理對氨揮發(fā)及其排放因子的影響

    同列不同小寫字母表示不同處理間差異顯著(<0.05)。Different lowercase letters within the same column indicate significant differences among treatments at< 0.05.

    2.4 施用豬場肥水對土壤NO3–-N含量的影響

    在玉米收獲期, 隨土壤深度增加, 土壤NO3–-N含量呈先降低后升高趨勢, 且在40~60 cm土層達到最小值(圖3A)。由公式(6)計算可得, 與CK相比, CK1、CK2、25%WB、50%WB、25%WI和50%WI在0~100 cm土壤剖面的硝態(tài)氮累積量顯著增加15.9%、38.3%、64.4%、89.2%、72.0%和31.7%。

    在小麥收獲期, 隨土壤深度的增加, 土壤NO3–-N呈現(xiàn)先升高后降低再升高的“S”型分布, 且在20~40 cm土層達到最大值, 在60~80 cm土層達最小值(圖3B)。由公式(6)計算可得, 與CK相比, CK1、CK2、25%WB、50%WB、25%WI和50%WI在0~100 cm土壤剖面的硝態(tài)氮累積量顯著增加53.7%、74.5%、74.4%、93.7%、81.0%和69.8%。

    圖3 不同豬場肥水施用處理下玉米(A)和小麥(B)收獲后0~180 cm土壤剖面NO3–-N分布

    2.5 豬場肥水施用下玉米-小麥輪作的土壤-作物系統(tǒng)氮素表觀平衡

    玉米-小麥輪作周期內土壤(0~100 cm)-作物系統(tǒng)氮素投入和產出見表6。系統(tǒng)氮素輸入項包括肥料(包括尿素、豬場肥水)、作物種子、前茬還田秸稈、大氣沉降、灌溉和試驗起始0~100 cm土層礦質氮累積量(起始Nmin)等形式的外源氮。土壤氮素輸出包括作物收獲攜出、試驗后0~100 cm土層礦質氮殘留(殘留Nmin)、氨揮發(fā)損失和其他形式的氮素損失。秸稈氮輸入量包括試驗開始前農田前茬小麥還田秸稈所含的氮素和玉米收獲后還田的玉米秸稈氮素。就氮輸入而言, CK處理起始Nmin占氮總輸入量的59.0%, 而其他處理起始Nmin占氮總輸入量的24.2%~25.2%; 尿素和豬場肥水施用處理肥料氮(尿素氮和肥水氮)占氮總輸入量的51.5%~53.2% (表6)。

    就氮輸出而言, 玉米和小麥兩季作物累計氮素吸收量為303.1~460.0 kg?hm–2, 且表現(xiàn)為CK

    與CK相比, 施用尿素和豬場肥水可以顯著增加玉米-小麥輪作周期內土壤氨揮發(fā)損失總量, 其增幅為4.9~6.1倍; 而與CK1相比, 其余各處理(除CK外)的土壤氨揮發(fā)損失總量均顯著降低, 其降幅為6.2%~16.6%, 且50%WI降幅最大。CK處理其他形式氮損失為-129.9 kg(N)?hm–2, 尿素和豬場肥水施用處理的其他形式氮損失為39.6~140.7 kg(N)?hm–2, 其中CK1處理最大, 而50%WB處理最小, 且在25%WB、50%WB、25%WI和50%WI之間無顯著性差異。就氮盈余而言, CK最小[32.3 kg(N)?hm–2], CK1最大[423.3 kg(N)?hm–2], 其次依次為CK2[409.2 kg(N)?hm–2]、25%WB[391.3 kg(N)?hm–2]、50%WB[376.7 kg(N)?hm–2]、50%WI[357.4 kg(N)?hm–2]和25%WI[355.3 kg(N)?hm–2], 50%WB、25%WI和50%WI處理間無顯著性差異。

    表6 不同豬場肥水施用處理下玉米-小麥輪作周期內土壤(0~100 cm)-作物系統(tǒng)氮素表觀平衡

    Nmin:礦質氮。同行不同小寫字母表示不同施肥處理間差異顯著(<0.05)。Nmin: mineral nitrogen. Different lowercase letters within the same row indicate significant differences among treatments at< 0.05.

    3 討論

    3.1 施用豬場肥水對土壤氨揮發(fā)和NO3–-N淋溶的影響

    本研究在玉米、小麥季施用尿素和豬場肥水后的氨排放因子為1.46%~15.41%, 玉米和小麥季氨排放總因子分別為7.85%~9.78%和8.20%~10.12%, 這一結果與蘇芳等[25]和董文旭等[26]在華北平原玉米-小麥輪作農田施用礦物氮肥后氨揮發(fā)損失率大致相當。在本研究中, 無論是玉米季還是小麥季, 尿素表施的氨排放均表現(xiàn)最高, 其次為尿素注射施用, 而豬場肥水施用處理可顯著降低土壤氨揮發(fā)速率及其氨揮發(fā)損失量, 其中豬場肥水替代50%尿素氮注射施用在玉米和小麥基肥期的氨排放因子僅為1.46%和8.08%, 遠低于尿素表施和尿素注射施用, 然而, 豬場肥水氮素不同替代率施用處理間無顯著差異。究其原因, 主要是尿素和豬場肥水施入土壤后發(fā)生的反應有所不同[27]: 施入土壤后的尿素在脲酶作用下很快水解成NH4HCO3, 進而轉化為NH4+-N, 從而為氨揮發(fā)提供較多的底物, 加大氨損失; 然而, 豬場肥水中活性有機物質促進土壤對肥水NH4+的吸附固定, 在一定程度上降低氨揮發(fā); 同時, 肥水注射處理將肥水施入到表層土壤(5~10 cm)中, 能夠隔絕肥水與空氣的接觸, 進一步增加土壤對NH4+的吸附固定[28]; 另外, 在硝化細菌的作用下, 施入土壤的豬場肥水氮將會轉化成NO3–-N, 從而有效降低氨揮發(fā)的發(fā)生[29]。然而, 也有研究表明向農田土壤直接施用畜禽肥水能夠提高氨揮發(fā)速率, 其原因可能為: 畜禽肥水中含有豐富NH4+-N, 且施用的農田土壤pH呈堿性, 在適宜的環(huán)境條件下肥水NH4+-N很容易排放到大氣中[30]; 同時, 畜禽肥水含有豐富活性碳源, 能夠激發(fā)表層土壤微生物活性, 促進土壤有機氮向簡單氮化合物的轉化, 進而通過脫氨基作用轉化為NH3釋放到大氣中[31]。本研究中, 除了小麥基肥期25%WB處理氨揮發(fā)量顯著低于25%WI外, 無論玉米基肥期, 還是小麥基肥期, 等N替代豬場肥水表施氨揮發(fā)量均高于豬場肥水注射施用, 且50%WI處理氨揮發(fā)量均顯著低于其他各施肥處理。因此, 在華北平原夏玉米-冬小麥一年兩熟輪作區(qū), 畜禽肥水施用過程中應避免直接撒施, 最好將肥水注射到表層土壤(5~10 cm)中, 降低土壤氨揮發(fā)損失。

    在玉米-小麥輪作農田土壤中, 土壤NO3–-N是作物吸收利用的主要氮素形態(tài), 其在土壤剖面上的分布受施氮量、灌溉和降水等因素影響[24]。研究表明, 長期過量施氮導致NO3–-N在土壤中累積, 且隨著施氮量的增加而顯著增加, 嚴重威脅地下水的安全[32]。在本研究中, 尿素或豬場肥水表施或注射(5~10 cm)施用導致NO3–-N集中在0~20 cm土層, 而富集在表層的NO3–-N在灌溉或降水等作用下均可向深層土壤遷移, 從而產生以下現(xiàn)象: 隨土壤的加深(主要80 cm以下), 土壤NO3–-N含量逐漸升高; 然而, 豬場肥水氮素不同替代率施用處理間無明顯差異。玉米和小麥收獲后NO3–-N在土壤剖面的分布趨勢并未保持一致, 其原因可能是在玉米生長中后期, 試驗所在地降水量過大, 大量NO3–-N被淋溶到深層土壤, 而小麥生長中后期降水量低, 大量NO3–-N保留在60 cm土層以上土壤中。同時, 該試驗區(qū)常年耕作深度<30 cm, 導致30~60 cm土層出現(xiàn)堅硬的犁底層, 從而導致玉米生長所需氮主要取決于60 cm以上土層。豬場肥水所含的大量有機、無機物質能夠有效協(xié)調土壤氮素供應從而降低氮素淋失; 同時豬場肥水中水、氮耦合, 能夠對土壤氮殘留和淋失起到相互促進或制約的影響[33]。在本研究中, 在玉米和小麥收獲后, 豬場肥水施用處理的土壤NO3–-N殘留量均在一定程度上高于尿素施用處理, 這與楊軍等[33]和杜會英等[34]研究結果一致。土壤膠體對NO3–-N吸附能力比較弱, 過量施氮導致大量NO3–-N累積在土壤剖面, 并在不合理灌溉或暴雨作用下向下淋溶。杜會英等[34]發(fā)現(xiàn), 在玉米-小麥輪作體系中, 20 kg(N)?hm–2的畜禽肥水氮投入導致80~100 cm土層積累大量NO3–-N(~18 kg?hm–2)。而本研究玉米-小麥輪作周期總施氮量遠高于320 kg(N)?hm–2, 因此, 在該體系下420 kg(N)?hm–2的氮投入量存在較大的環(huán)境風險。

    3.2 施用豬場肥水對作物產量、氮素利用率和氮素表觀平衡的影響

    豬場肥水中含有豐富的水溶性有機質和氮、磷、鉀等活性養(yǎng)分, 能夠被土壤微生物直接吸收利用, 利于后茬作物根際土壤微生物活性的提高和土壤氮、磷和鉀等養(yǎng)分的有效形態(tài)轉化; 同時, 土壤微生物能夠將作物來不及吸收利用的部分養(yǎng)分轉化成有機形態(tài), 當作物對養(yǎng)分需求增加時, 土壤微生物便能夠將其所固定的養(yǎng)分釋放出來, 以供作物生長所需, 進而利于作物對養(yǎng)分的吸收及其產量的提高[34-37]。豬場肥水施用的同時, 向農田土壤輸入大量水資源, 既節(jié)省水資源, 又利于種子的發(fā)芽和出苗。本研究表明, 與尿素施用處理相比, 豬場肥水施用處理, 特別是25%WI, 能夠顯著提高玉米和小麥籽粒產量。Martínez等[36]研究也表明隨著豬場廢水灌溉量的增加, 玉米產量逐漸增加, 且當豬場廢水灌溉氮帶入量為200 kg(N)?hm–2時, 玉米產量最高; 杜會英等[7]也發(fā)現(xiàn), 牛場肥水灌溉氮帶入量為240 kg(N)?hm–2時, 小麥產量最高。本研究中, 與尿素施用處理相比, 豬場肥水施用處理也能提高作物對氮的吸收利用。陶曉婷等[6]研究發(fā)現(xiàn), 施用豬場肥水灌溉能夠有效提高小麥與氮素供應水平密切相關的葉片SAPD值和植株含氮量。因此, 向農田施用豬場肥水能夠在一定程度上提高作物產量、促進氮素利用。

    氮肥農學效率、氮肥偏生產力和氮肥表觀利用率是衡量氮素效率的常用指標, 從不同方面描述作物對氮素的吸收利用[38]。Dobermann[38]認為糧食作物氮肥農學效率為10~30 kg?kg–1、氮肥偏生產力為40~70 kg?kg–1和氮肥表觀利用率為30%~50%較為適宜。在本研究中, 25%WI和50%WI處理的夏玉米生育期氮肥農學效率、氮肥偏生產力和氮肥表觀利用率介于上述適宜范圍內, 而其余各處理均低于以上范圍(除50%WB處理的氮肥農學效率和氮肥表觀利用率), 說明在玉米生育期180 kg?hm–2的施氮量下, 豬場肥水替代25%和50%尿素氮注射施用能夠獲得較好的養(yǎng)分效率; 然而, 在小麥生育期氮肥農學效率、氮肥偏生產力和氮肥表觀利用效率均低于上述適宜范圍, 因此, 在小麥生育期240 kg?hm–2的施氮量下, 無論采用何種施肥措施, 都不能取得較好的養(yǎng)分效率, 盡管豬場肥水替代25%N和50%N尿素氮注射施用的氮肥農學效率、氮肥偏生產力和氮肥表觀利用率明顯高于其他處理。豬場肥水注射施用之所以能夠較其他施肥措施提高作物氮素利用率, 是因為注射技術將肥水施入到5~10 cm土層, 能夠最大程度降低土壤氨揮發(fā)等氮損失, 進而為作物提供生長必需的水分和養(yǎng)分, 利于作物的發(fā)芽和前期生長。張福鎖等[39]報道中國玉米平均施氮量為162 kg?hm–2, 產量、氮肥農學效率、氮肥偏生產力和氮肥表觀利用效率分別為7.05 Mg?hm–2、9.8 kg?kg–1和51.6 kg?kg–1和26.1%; 小麥平均施氮量為169 kg?hm–2, 而產量、氮肥農學效率和氮肥偏生產力分別為5.72 Mg?hm–2、8.0 kg?kg–1、43.0 kg?kg–1和28.2%。在本研究中, 玉米生育期各施肥處理的玉米籽粒產量均高于7.05 Mg?hm–2, 25%WB、50%WB、25%WI和50%WI的氮肥農學效率高于9.8 kg?kg–1, 氮肥表觀利用率高于26.1%, 然而, 各處理的氮肥偏生產力均低于51.6 kg?kg–1, 小麥生育期各施肥處理的小麥籽粒產量均高于5.72 Mg?hm–2, 而各處理的氮肥農學效率均低于8.0 kg?kg–1, 氮肥偏生產力低于43.0 kg?kg–1, 而氮肥表觀利用率亦低于28.2%。其原因可能有以下兩點: 第一, 施氮量過高(玉米和小麥生育期分別為240 kg?hm–2和180 kg?hm–2); 第二, 試驗開始前土壤礦質氮累積量過高(197.8 kg?hm–2), 導致不施肥對照區(qū)作物生物量和氮素攜出量過高。因此, 在該區(qū)玉米-小麥輪作生產中, 無論采取直接施用尿素, 還是施用豬場肥水, 都應適當降低氮素投入量; 同時在施氮量一致的情況下, 采用豬場肥水替代尿素注射施用能夠獲得較高的氮素利用效率。

    分析農田土壤-作物系統(tǒng)氮素平衡狀況, 既能掌握作物對氮素施用的響應, 又能了解土壤氮素盈余狀況及肥料氮的去向, 進而為指導農業(yè)生產提供理論依據。有研究表明, 土壤礦質氮的作物有效性隨著土壤深度的增加而逐漸降低, 且0~100 cm土體中的礦質氮能夠被玉米和小麥根系直接吸收利用, 而100 cm土層以下礦質氮難以被作物吸收利用, 最終只能通過淋失進入地下水, 或通過深層反硝化作用脫氮[40]。因此, 本研究采取0~100 cm土體作為計算對象, 對玉米-小麥一年兩熟輪作體系農田土壤-作物系統(tǒng)的氮素平衡進行分析。對照處理氮盈余最低, 且氨揮發(fā)損失之外其他形式氮損失小于0。石德楊等[41]在土壤高殘留氮條件下夏玉米生育期礦化氮可達184 kg(N)?hm–2。因此, 本研究中, CK處理其他形式氮損失量可能為土壤有機氮礦化所致。在本研究尿素和豬場肥水施用各處理中, 25%WI處理的作物氮素吸收量最高, 氮盈余最低, 然而0~100 cm土體殘留礦質氮累積量仍維持在較高水平, 同時, 50%WI處理的作物籽粒產量、氮素吸收量和氮盈余均與25%WI無顯著差異, 表明豬場肥水替代25%和50%尿素氮注射施用利于作物生長和土壤氮素平衡的維持。由于本研究結果僅為1年的田間試驗結果, 對豬場肥水注射施用下氮素替代比率及作物氮素吸收利用等尚需進一步驗證。

    另外, 在灌溉條件不完善的條件下, 向農田注射施用豬場肥水能夠為種子發(fā)芽提供必須的水分和養(yǎng)分, 特別是玉米播種期。此外, 采用豬場肥水替代尿素氮施用, 既能降低生產成本, 提高經濟效益, 又能為規(guī)?;i養(yǎng)殖場解決廢水問題, 進而利于農牧的協(xié)調發(fā)展, 因此, 在不降低土壤質量的前提下, 綜合作物產量、土壤(0~100 cm)-作物系統(tǒng)氮素表觀平衡, 豬場肥水替代25%和50%尿素氮注射施用是華北平原豬場肥水施用的最佳施用方式。

    4 結論

    與不施肥處理相比, 施用尿素和豬場肥水均可顯著提高玉米、小麥產量和籽粒氮吸收量, 且豬場肥水替代25%和50%尿素氮注射施用的增幅顯著高于其他處理。與尿素表施相比, 尿素注射施用、肥水表施和注射施用在玉米和小麥基肥期的土壤氨損失總量均可顯著降低, 且豬場肥水替代50%尿素氮注射施用降幅最大; 相反, 農學利用效率、氮肥偏生產力和氮肥表觀利用率顯著提高, 且表現(xiàn)為肥水注射施用>肥水表施>尿素注射施用>尿素表施, 而豬場肥水替代25%和50%尿素氮注射施用之間無顯著差異。在玉米-小麥輪作周期內土壤(0~100 cm)-作物系統(tǒng)中, 與尿素表施相比, 尿素注射施用和豬場肥水施用, 特別是豬場肥水替代25%和50%尿素氮注射施用能夠顯著降低系統(tǒng)氮損失和氮盈余。因此, 在華北平原玉米-小麥高度集約化輪作區(qū), 豬場肥水最佳施用方式為替代25%和50%尿素氮注射施用。

    [1] BAI Z H, MA L, JIN S Q, et al. Nitrogen, phosphorus, and potassium flows through the manure management chain in China[J]. Environmental Science & Technology, 2016, 50(24): 13409–13418

    [2] WANG Y, DONG H M, ZHU Z P, et al. Mitigating greenhouse gas and ammonia emissions from swine manure management: A system analysis[J]. Environmental Science & Technology, 2017, 51(8): 4503–4511

    [3] SCHLEGEL A J, ASSEFA Y, BOND H D, et al. Changes in soil nutrients after 10 years of cattle manure and swine effluent application[J]. Soil and Tillage Research, 2017, 172: 48–58

    [4] VETSCH J A, SCHERDER E F, RUEN D C. Does liquid swine manure application timing and nitrapyrin affect corn yield and inorganic soil nitrogen?[J]. Agronomy Journal, 2017, 109(5): 2358–2370

    [5] HAO X J, ZHANG T Q, WANG Y T, et al. Soil test phosphorus and phosphorus availability of swine manures with long-term application[J]. Agronomy Journal, 2018, 110(5): 1943–1950

    [6] 陶曉婷, 朱正杰, 高威, 等. 規(guī)模化豬場處理廢水與化肥配施對小麥氮素吸收利用的影響[J]. 農業(yè)環(huán)境科學學報, 2014, 33(3): 555–561 TAO X T, ZHU Z J, GAO W, et al. Nitrogen uptake and utilization in wheat as influenced by pig slurry from large-scale pig farm[J]. Journal of Agro-Environment Science, 2014, 33(3): 555–561

    [7] 杜會英, 馮潔, 張克強, 等. 牛場肥水灌溉對冬小麥產量與氮利用效率及土壤硝態(tài)氮的影響[J]. 植物營養(yǎng)與肥料學報, 2016, 22(2): 536–541 DU H Y, FENG J, ZHANG K Q, et al. Effects of dairy effluents irrigation on yield and nitrogen use efficiency of winter wheat and soil nitrate nitrogen[J]. Journal of Plant Nutrition and Fertilizer, 2016, 22(2): 536–541

    [8] ARONSSON H, LIU J, EKRE E, et al. Effects of pig and dairy slurry application on N and P leaching from crop rotations with spring cereals and forage leys[J]. Nutrient Cycling in Agroecosystems, 2014, 98(3): 281–293

    [9] GIROTTO E, CERETTA C A, LOURENZI C R, et al. Nutrient transfers by leaching in a no-tillage system through soil treated with repeated pig slurry applications[J]. Nutrient Cycling in Agroecosystems, 2013, 95(1): 115–131

    [10] GONZATTO R, AITA C, BéLANGER G, et al. Response of no-till grain crops to pig slurry application methods and a nitrification inhibitor[J]. Agronomy Journal, 2017, 109(4): 1687–1696

    [11] AITA C, GONZATTO R, MIOLA E C C, et al. Injection of dicyandiamide-treated pig slurry reduced ammonia volatilization without enhancing soil nitrous oxide emissions from no-till corn in southern Brazil[J]. Journal of Environment Quality, 2014, 43(3): 789–800

    [12] PIETZNER B, RüCKNAGEL J, KOBLENZ B, et al. Impact of slurry strip-till and surface slurry incorporation on NH3and N2O emissions on different plot trials in central Germany[J]. Soil and Tillage Research, 2017, 169: 54–64

    [13] CONGREVES K A, GRANT B B, DUTTA B, et al. Predicting ammonia volatilization after field application of swine slurry: DNDC model development[J]. Agriculture, Ecosystems & Environment, 2016, 219: 179–189

    [14] JING Q, JéGO G, BéLANGER G, et al. Simulation of water and nitrogen balances in a perennial forage system using the STICS model[J]. Field Crops Research, 2017, 201: 10–18

    [15] HERNáNDEZ D, POLO A, PLAZA C. Long-term effects of pig slurry on barley yield and N use efficiency under semiarid Mediterranean conditions[J]. European Journal of Agronomy, 2013, 44: 78–86

    [16] MORENO-GARCíA B, GUILLéN M, QUíLEZ D. Response of paddy rice to fertilisation with pig slurry in Northeast Spain: Strategies to optimise nitrogen use efficiency[J]. Field Crops Research, 2017, 208: 44–54

    [17] PLAZA-BONILLA D, CANTERO-MARTíNEZ C, BARECHE J, et al. Do no-till and pig slurry application improve barley yield and water and nitrogen use efficiencies in rainfed Mediterranean conditions?[J]. Field Crops Research, 2017, 203: 74–85

    [18] 劉學軍, 趙紫娟, 巨曉棠, 等. 基施氮肥對冬小麥產量、氮肥利用率及氮平衡的影響[J]. 生態(tài)學報, 2002, 22(7): 1122–1128 LIU X J, ZHAO Z J, JU X T, et al. Effect of N application as basal fertilizer on grain yield of winter wheat, fertilizer N recovery and N balance[J]. Acta Ecologica Sinica, 2002, 22(7): 1122–1128

    [19] 楊曉卡, 米慧玲, 高韓鈺, 等. 不同栽培模式對冬小麥-夏玉米輪作系統(tǒng)產量、氮素累積和平衡的影響[J]. 應用生態(tài)學報, 2016, 27(6): 1935–1941 YANG X K, MI H L, GAO H Y, et al. Effects of different cultivation patterns on yield, nitrate accumulation and nitrogen balance in winter wheat and summer maize rotation system[J]. Chinese Journal of Applied Ecology, 2016, 27(6): 1935–1941

    [20] JU X T, XING G X, CHEN X P, et al. Reducing environmental risk by improving N management in intensive Chinese agricultural systems[J]. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(9): 3041–3046

    [21] PAN D, KONG F B, ZHANG N, et al. Knowledge training and the change of fertilizer use intensity: Evidence from wheat farmers in China[J]. Journal of Environmental Management, 2017, 197: 130–139

    [22] LIU Z, GAO J, GAO F, et al. Integrated agronomic practices management improve yield and nitrogen balance in double cropping of winter wheat-summer maize[J]. Field Crops Research, 2018, 221: 196–206

    [23] SUN L Y, WU Z, MA Y C, et al. Ammonia volatilization and atmospheric N deposition following straw and urea application from a rice-wheat rotation in southeastern China[J]. Atmospheric Environment, 2018, 181: 97–105

    [24] 石鵬飛, 鄭媛媛, 趙平, 等. 華北平原種養(yǎng)一體規(guī)?;r場氮素流動特征及利用效率——以河北津龍循環(huán)農業(yè)園區(qū)為例[J]. 應用生態(tài)學報, 2017, 28(4): 1281–1288 SHI P F, ZHENG Y Y, ZHAO P, et al. Nitrogen flow characteristic and use efficiency in mixed crop-livestock system in North China Plain: Hebei Jinlong circular agriculture farm as an example[J]. Chinese Journal of Applied Ecology, 2017, 28(4): 1281–1288

    [25] 蘇芳, 丁新泉, 高志嶺, 等. 華北平原冬小麥-夏玉米輪作體系氮肥的氨揮發(fā)[J]. 中國環(huán)境科學, 2007, 27(3): 409–413 SU F, DING X Q, GAO Z L, et al. Ammonia volatilization from nitrogen fertilization of winter wheat-summer maize rotation system in the North China Plain[J]. China Environmental Science, 2007, 27(3): 409–413

    [26] 董文旭, 吳電明, 胡春勝, 等. 華北山前平原農田氨揮發(fā)速率與調控研究[J]. 中國生態(tài)農業(yè)學報, 2011, 19(5): 1115–1121 DONG W X, WU D M, HU C S, et al. Ammonia volatilization and control mechanisms in the piedmont of North China plain[J]. Chinese Journal of Eco-Agriculture, 2011, 19(5): 1115–1121

    [27] 鄭鳳霞, 董樹亭, 劉鵬, 等. 長期有機無機肥配施對冬小麥籽粒產量及氨揮發(fā)損失的影響[J]. 植物營養(yǎng)與肥料學報, 2017, 23(3): 567–577 ZHENG F X, DONG S T, LIU P, et al. Effects of combined application of manure and chemical fertilizers on ammonia volatilization loss and yield of winter wheat[J]. Journal of Plant Nutrition and Fertilizer, 2017, 23(3): 567–577

    [28] THOMPSON R B, PAIN B F, REES Y J. Ammonia volatilization from cattle slurry following surface application to grassland[J]. Plant and Soil, 1990, 125(1): 119–128

    [29] 陳效民, 吳華山, 孫靜紅. 太湖地區(qū)農田土壤中銨態(tài)氮和硝態(tài)氮的時空變異[J]. 環(huán)境科學, 2006, 27(6): 1217–1222 CHEN X M, WU H S, SUN J H. Time-spatial variability of ammonium and nitrate in farmland soil of Taihu Lake region[J]. Environmental Science, 2006, 27(6): 1217–1222

    [30] MARTINES A M, NOGUEIRA M A, SANTOS C A, et al. Ammonia volatilization in soil treated with tannery sludge[J]. Bioresource Technology, 2010, 101(12): 4690–4696

    [31] SUBKE J, HAHN V, BATTIPAGLIA G, et al. Feedback interactions between needle litter decomposition and rhizosphere activity[J]. Oecologia, 2004, 139(4): 551–559

    [32] BASSO B, RITCHIE J T. Impact of compost, manure and inorganic fertilizer on nitrate leaching and yield for a 6-year maize-alfalfa rotation in Michigan[J]. Agriculture, Ecosystems & Environment, 2005, 108(4): 329–341

    [33] 楊軍, 張蕾, 張克強, 等. 豬場廢水灌溉對潮土硝態(tài)氮含量變化的影響[J]. 農業(yè)工程學報, 2009, 25(5): 35–39 YANG J, ZHANG L, ZHANG K Q, et al. Effects of irrigation with piggery wastewater on nitrate nitrogen dynamics in fluvio-aquatic soil[J]. Transactions of the CSAE, 2009, 25(5): 35–39

    [34] 杜會英, 馮潔, 郭海剛, 等. 麥季牛場肥水灌溉對冬小麥-夏玉米輪作土壤氮素平衡的影響[J]. 農業(yè)工程學報, 2015, 31(3): 159–165 DU H Y, FENG J, GUO H G, et al. Effects of dairy effluents irrigation on N balance in soil under winter wheat-summer maize rotation system[J]. Transactions of the CSAE, 2015, 31(3): 159–165

    [35] LALOR S T J, SCHR?DER J J, LANTINGA E A, et al. Nitrogen fertilizer replacement value of cattle slurry in grassland as affected by method and timing of application[J]. Journal of Environment Quality, 2011, 40(2): 362

    [36] MARTíNEZ E, MARESMA A, BIAU A, et al. Long-term effects of pig slurry combined with mineral nitrogen on maize in a Mediterranean irrigated environment[J]. Field Crops Research, 2017, 214: 341–349

    [37] 彭少兵, 黃見良, 鐘旭華, 等. 提高中國稻田氮肥利用率的研究策略[J]. 中國農業(yè)科學, 2002, 35(9): 1095–1103 PENG S B, HUANG J L, ZHONG X H, et al. Research strategy in improving fertilizer-nitrogen use efficiency of irrigated rice in China[J]. Scientia Agricultura Sinica, 2002, 35(9): 1095–1103

    [38] DOBERMANN A R. Nitrogen use efficiency —State of the art[C]//IFA International Workshop on Enhanced Efficiency Fertilizers. Frankfurt, Germany, 2005

    [39] 張福鎖, 王激清, 張衛(wèi)峰, 等. 中國主要糧食作物肥料利用率現(xiàn)狀與提高途徑[J]. 土壤學報, 2008, 45(5): 915–924 ZHANG F S, WANG J Q, ZHANG W F, et al. Nutrient use efficiencies of major cereal crops in China and measures for improvement[J]. Acta Pedologica Sinica, 2008, 45(5): 915–924

    [40] 劉學軍, 巨曉棠, 張福鎖. 減量施氮對冬小麥-夏玉米種植體系中氮利用與平衡的影響[J]. 應用生態(tài)學報, 2004, 15(3): 458–462 LIU X J, JU X T, ZHANG F S. Effect of reduced N application on N utilization and balance in winter wheat-summer maize cropping system[J]. Chinese Journal of Applied Ecology, 2004, 15(3): 458–462

    [41] 石德楊, 張海艷, 董樹亭. 土壤高殘留氮條件下施氮對夏玉米氮素平衡、利用及產量的影響[J]. 植物營養(yǎng)與肥料學報, 2013, 19(1): 37–44 SHI D Y, ZHANG H Y, DONG S T. Effects of nitrogen application on nitrogen balance and use efficiency and yield of summer maize in soil with high residual nitrogen[J]. Plant Nutrition and Fertilizer Science, 2013, 19(1): 37–44

    Effects of swine slurry application on ammonia emission, nitrogen utilization and apparent balance of a winter wheat-summer maize rotation system*

    LI Shuo1, WANG Xuan1, ZHANG Xiqun2, LIU Zelong1, ZHAO Hao1, ZHAO Zhanqing1, ZHANG Yuming1, SUN Hongyong1, MA Lin1**

    (1. Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences / Key Laboratory of Agricultural Water Resources, Chinese Academy of Sciences /Hebei Key Laboratory of Soil Ecology, Shijiazhuang 050022, China; 2. Hebei Agricultural Mechanization Institute Co. Ltd., Shijiazhuang 050022, China)

    With the development of intensive swine production in China, large amounts of slurry are produced, which causes nitrogen and phosphorus loss and serious non-point pollution. A potentially efficient way to solve these environmental problems is to substitute mineral fertilizer with swine slurry in intensive maize (L.)-wheat (L.) double-cropping rotation systems, which could promote the development of sustainable agriculture and animal husbandry. A field experiment was performed with swine slurry application in a maize-wheat double-cropping rotation system in the North China Plain. The study included the following seven treatments with three replications: zero-N control (CK); urea broadcast (CK1); urea injection (CK2); swine slurry as a substitute for 25% urea-N broadcast (25%WB), 50% urea-N broadcast (50%WB), 25% urea-N injection (25%WI), and 50% urea-N injection (50%WI). Swine slurry was applied to the soil before maize and wheat seeding. Compared with the CK, the application of urea and swine slurry significantly improved maize and wheat grain yield and N uptake, with the greatest effect seen in the 25%WI, followed by the 50%WI. Compared with the CK1, the N agronomic efficiency, N partial productivity, and apparent N recovery rate under treatments of urea injection and swine slurry broadcast and injection were significantly improved. The greatest improvement was seen in treatments of swine slurry injection, followed by treatments of swine slurry broadcast, but no significant difference was found between the 25%WI and 50%WI treatments. Compared with CK, the application of urea and swine slurry significantly improved nitrate accumulation. The increase in inorganic nitrogen ranged from 50.8%-87.9% throughout the 0-100 cm soil profile after maize harvest. No significant difference in nitrate accumulation was found among the 50%WB, 25%WI, and 50%WI treatments. Compared with CK1, swine slurry broadcast and urea and swine slurry injection significantly reduced the total amount of ammonia loss during the seeding period of maize and wheat, with reductions ranging from 26.5% to 48.6% and 11.4% to 29.1%, respectively. Meanwhile, swine slurry broadcast and injection significantly reduced N surplus from 7.6% to 16.0%, with the biggest difference found in the 25%WI treatments. However, no significant difference was found between the 25%WI and 50%WI treatments. In view of the yield, nitrogen utilization, and environmental benefits, swine slurry as a substitute for 25% and 50% urea-N injection were reasonable methods of swine slurry application in the intensive summer maize-winter wheat double-cropping rotation field.

    Maize-wheat rotation field; Swine slurry; Injection; Ammonia emission; Nitrogen utilization; Nitrogen apparent balance

    , E-mail: malin1979@sjziam.ac.cn

    Feb. 28, 2019;

    May 20, 2019

    S143.1

    2096-6237(2019)10-1502-13

    10.13930/j.cnki.cjea.190150

    馬林, 主要研究方向為農業(yè)生態(tài)學與養(yǎng)分資源管理。E-mail: malin1979@sjziam.ac.cn

    李碩, 主要研究方向為農業(yè)生態(tài)學與養(yǎng)分資源管理。E-mail: lishuo@sjziam.ac.cn

    2019-02-28

    2019-05-20

    * This study was supported by the National Key R&D Program of China (2018YFC0213300), the National Natural Science Foundation of China (31872403, 31801943), the Science and Technology Service Network Initiative of Chinese Academy of Sciences (KFJ-STS-ZDTP-053), the National Research Program for Key Issues in Air Pollution Control of China (DQGG0208), the Key Research Program of the Chinese Academy Sciences (ZDRW-ZS-2016-5),the 100-Talent Project of Chinese Academy of Sciences, the Distinguished Young Scientists Project of Natural Science Foundation of Hebei (D2017503023) and the Hebei Dairy Cattle Innovation Team of Modern Agro-industry Technology Research System (HBCT2018120206).

    * 國家重點研發(fā)計劃(2018YFC0213300)、國家自然基金項目(31872403, 31801943)、中國科學院科技服務網絡計劃(STS)項目(KFJ-STS-ZDTP-053)、大氣重污染成因與治理攻關項目(DQGG0208)、中國科學院重點部署項目(ZDRW-ZS-2016-5)、中國科學院百人計劃項目、河北省杰出青年基金項目(D2017503023)和河北省現(xiàn)代農業(yè)產業(yè)技術體系奶牛產業(yè)創(chuàng)新團隊項目(HBCT2018120206)資助

    李碩, 王選, 張西群, 劉澤龍, 趙浩, 趙占輕, 張玉銘, 孫宏勇, 馬林. 豬場肥水施用對玉米-小麥農田氨排放、氮素利用與表觀平衡的影響[J]. 中國生態(tài)農業(yè)學報(中英文), 2019, 27(10): 1502-1514

    LI S, WANG X, ZHANG X Q, LIU Z L, ZHAO H, ZHAO Z Q, ZHANG Y M, SUN H Y, MA L.Effects of swine slurry application on ammonia emission, nitrogen utilization and apparent balance of a winter wheat-summer maize rotation system[J]. Chinese Journal of Eco-Agriculture, 2019, 27(10): 1502-1514

    猜你喜歡
    肥水表觀氮素
    綠盲蝽為害與赤霞珠葡萄防御互作中的表觀響應
    河北果樹(2021年4期)2021-12-02 01:14:50
    早春水難肥,了解原因最重要!正確肥水,這三樣必須補
    當代水產(2021年4期)2021-07-20 08:11:00
    掌握這6點要素,讓肥水更高效
    當代水產(2020年4期)2020-06-16 03:23:30
    沿??O蟶養(yǎng)殖春季肥水技術
    當代水產(2020年4期)2020-06-16 03:23:26
    鋼結構表觀裂紋監(jiān)測技術對比與展望
    上海公路(2019年3期)2019-11-25 07:39:28
    例析對高中表觀遺傳學的認識
    既要環(huán)保,還要產量,魚和熊掌如何兼得——談新形式下如何肥水
    當代水產(2018年8期)2018-11-02 05:30:50
    表觀遺傳修飾在糖脂代謝中的作用
    遺傳(2014年3期)2014-02-28 20:58:52
    楸樹無性系苗期氮素分配和氮素效率差異
    基于光譜分析的玉米氮素營養(yǎng)診斷
    中国国产av一级| 一本久久精品| 丝袜美腿诱惑在线| 女人高潮潮喷娇喘18禁视频| 三上悠亚av全集在线观看| 丝瓜视频免费看黄片| 午夜日韩欧美国产| 伦理电影免费视频| 亚洲情色 制服丝袜| 99九九在线精品视频| 麻豆乱淫一区二区| 午夜视频精品福利| 国产xxxxx性猛交| 两个人免费观看高清视频| av有码第一页| 国产成人精品久久二区二区91| 肉色欧美久久久久久久蜜桃| 欧美变态另类bdsm刘玥| 老司机靠b影院| 极品少妇高潮喷水抽搐| 亚洲成人国产一区在线观看| 国产一卡二卡三卡精品| 亚洲成av片中文字幕在线观看| 欧美亚洲日本最大视频资源| 国产在线一区二区三区精| 美国免费a级毛片| 男人爽女人下面视频在线观看| 国产日韩欧美视频二区| 老熟妇仑乱视频hdxx| 十八禁高潮呻吟视频| 啦啦啦免费观看视频1| a级片在线免费高清观看视频| 一二三四社区在线视频社区8| 国产免费现黄频在线看| 亚洲七黄色美女视频| 一级片免费观看大全| 久久香蕉激情| 久久久久久久久久久久大奶| 一区二区三区激情视频| 妹子高潮喷水视频| 伦理电影免费视频| 久久青草综合色| 中文字幕精品免费在线观看视频| 国产亚洲一区二区精品| 黄网站色视频无遮挡免费观看| 亚洲国产成人一精品久久久| 精品亚洲成a人片在线观看| 美女中出高潮动态图| 捣出白浆h1v1| 久久人妻福利社区极品人妻图片| av超薄肉色丝袜交足视频| 美女主播在线视频| 国产亚洲欧美精品永久| 精品少妇内射三级| 一区二区三区激情视频| 国产一区有黄有色的免费视频| 国产淫语在线视频| 在线观看免费高清a一片| 亚洲中文日韩欧美视频| 丰满少妇做爰视频| 亚洲成人手机| 精品国内亚洲2022精品成人 | 老司机靠b影院| 国产野战对白在线观看| 午夜福利乱码中文字幕| 欧美激情极品国产一区二区三区| 亚洲欧美精品自产自拍| 久久精品亚洲av国产电影网| 最新在线观看一区二区三区| 男女国产视频网站| netflix在线观看网站| 淫妇啪啪啪对白视频 | 亚洲黑人精品在线| 欧美大码av| 亚洲精品在线美女| 90打野战视频偷拍视频| 自拍欧美九色日韩亚洲蝌蚪91| 51午夜福利影视在线观看| 亚洲伊人久久精品综合| 超碰成人久久| 黄片小视频在线播放| 国产黄频视频在线观看| 精品高清国产在线一区| 啦啦啦免费观看视频1| 日韩欧美免费精品| 亚洲美女黄色视频免费看| 人人妻,人人澡人人爽秒播| 999久久久国产精品视频| 国产成人系列免费观看| 免费在线观看完整版高清| a级片在线免费高清观看视频| 视频区欧美日本亚洲| 免费观看av网站的网址| 黄片大片在线免费观看| 一本色道久久久久久精品综合| 久久久久久久久免费视频了| 国产老妇伦熟女老妇高清| 少妇粗大呻吟视频| 精品国产一区二区三区四区第35| 十八禁网站网址无遮挡| 国产免费福利视频在线观看| 最黄视频免费看| 亚洲avbb在线观看| 亚洲午夜精品一区,二区,三区| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 国产人伦9x9x在线观看| 自线自在国产av| 黄色视频,在线免费观看| 欧美黄色淫秽网站| 一级黄色大片毛片| 国产免费现黄频在线看| 国产成人啪精品午夜网站| 国产亚洲av高清不卡| 欧美黑人欧美精品刺激| 久久久久久亚洲精品国产蜜桃av| 亚洲专区国产一区二区| 亚洲精品美女久久久久99蜜臀| av福利片在线| a级毛片黄视频| 1024香蕉在线观看| 极品人妻少妇av视频| 高清欧美精品videossex| 久久狼人影院| av福利片在线| 丝袜脚勾引网站| 国产伦理片在线播放av一区| 成人手机av| 中文字幕av电影在线播放| av天堂在线播放| 久久人妻熟女aⅴ| 国产激情久久老熟女| 秋霞在线观看毛片| 最近中文字幕2019免费版| 亚洲综合色网址| 一边摸一边做爽爽视频免费| 国产免费现黄频在线看| 婷婷丁香在线五月| 婷婷色av中文字幕| 亚洲av日韩在线播放| www.熟女人妻精品国产| 国产在线一区二区三区精| 满18在线观看网站| 久久国产精品人妻蜜桃| 国产亚洲精品第一综合不卡| 美女大奶头黄色视频| 老司机午夜福利在线观看视频 | 免费在线观看日本一区| 中文字幕精品免费在线观看视频| 亚洲精品成人av观看孕妇| 亚洲精品国产av蜜桃| 久久国产精品男人的天堂亚洲| www.999成人在线观看| 一级,二级,三级黄色视频| 日本一区二区免费在线视频| 国产日韩一区二区三区精品不卡| 在线 av 中文字幕| 巨乳人妻的诱惑在线观看| 1024香蕉在线观看| 亚洲五月色婷婷综合| 欧美av亚洲av综合av国产av| 高清视频免费观看一区二区| 亚洲一区二区三区欧美精品| 欧美日韩亚洲国产一区二区在线观看 | 叶爱在线成人免费视频播放| 国产欧美日韩一区二区精品| 两人在一起打扑克的视频| 国产黄色免费在线视频| 18禁黄网站禁片午夜丰满| 91老司机精品| 久久国产精品影院| 国产日韩欧美在线精品| 免费观看av网站的网址| 亚洲欧美日韩高清在线视频 | 大陆偷拍与自拍| 欧美另类亚洲清纯唯美| 成年av动漫网址| 99国产精品一区二区蜜桃av | 美女扒开内裤让男人捅视频| 91字幕亚洲| 亚洲精品中文字幕一二三四区 | 妹子高潮喷水视频| 99国产精品一区二区蜜桃av | 丰满人妻熟妇乱又伦精品不卡| 国产成人欧美| 日韩人妻精品一区2区三区| 日本av手机在线免费观看| 精品亚洲成a人片在线观看| 免费看十八禁软件| 人妻人人澡人人爽人人| 老汉色av国产亚洲站长工具| 夫妻午夜视频| 免费黄频网站在线观看国产| 别揉我奶头~嗯~啊~动态视频 | 亚洲专区国产一区二区| 亚洲九九香蕉| 久久久久久久精品精品| 午夜福利乱码中文字幕| 十分钟在线观看高清视频www| 亚洲专区国产一区二区| 视频区欧美日本亚洲| 亚洲国产中文字幕在线视频| 精品人妻一区二区三区麻豆| 亚洲天堂av无毛| 我的亚洲天堂| 欧美亚洲日本最大视频资源| 99re6热这里在线精品视频| 亚洲精品自拍成人| 夫妻午夜视频| 精品福利观看| 日韩 欧美 亚洲 中文字幕| 少妇裸体淫交视频免费看高清 | 天天添夜夜摸| 三级毛片av免费| 啦啦啦视频在线资源免费观看| 亚洲精品国产av蜜桃| 岛国毛片在线播放| 女警被强在线播放| 淫妇啪啪啪对白视频 | 国产精品一区二区在线观看99| 国产国语露脸激情在线看| 亚洲av成人一区二区三| 久久天堂一区二区三区四区| 日韩中文字幕视频在线看片| 日韩制服丝袜自拍偷拍| 国产av又大| 免费在线观看视频国产中文字幕亚洲 | 两人在一起打扑克的视频| 黑人巨大精品欧美一区二区mp4| 亚洲人成77777在线视频| 欧美精品av麻豆av| 美女中出高潮动态图| 操出白浆在线播放| 动漫黄色视频在线观看| 1024视频免费在线观看| 岛国在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 在线观看免费高清a一片| 国产精品国产av在线观看| 精品久久久久久电影网| 欧美 日韩 精品 国产| 久久精品国产亚洲av高清一级| 一级毛片女人18水好多| 国产伦理片在线播放av一区| 免费av中文字幕在线| 黑人巨大精品欧美一区二区蜜桃| 欧美成狂野欧美在线观看| 色婷婷久久久亚洲欧美| 最近中文字幕2019免费版| 淫妇啪啪啪对白视频 | 99热网站在线观看| 午夜精品国产一区二区电影| 国产色视频综合| 超碰成人久久| 丝袜美腿诱惑在线| 精品卡一卡二卡四卡免费| 午夜福利,免费看| 成年女人毛片免费观看观看9 | 午夜福利视频精品| 亚洲成av片中文字幕在线观看| 一级片'在线观看视频| 丝瓜视频免费看黄片| 久久免费观看电影| 黑人欧美特级aaaaaa片| 成年女人毛片免费观看观看9 | 老司机深夜福利视频在线观看 | 老司机午夜十八禁免费视频| 国产av国产精品国产| 国产成人欧美| 欧美日本中文国产一区发布| 免费在线观看黄色视频的| 亚洲av成人一区二区三| 99国产极品粉嫩在线观看| 性色av乱码一区二区三区2| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品一区二区精品视频观看| 搡老岳熟女国产| 电影成人av| 久久国产精品人妻蜜桃| 中文字幕制服av| 国产精品一区二区在线观看99| 久久久久视频综合| 午夜激情av网站| 国产91精品成人一区二区三区 | 三级毛片av免费| 欧美 亚洲 国产 日韩一| 在线观看免费视频网站a站| 日韩一卡2卡3卡4卡2021年| 国产成人欧美| 久久性视频一级片| 国产精品一区二区精品视频观看| 亚洲黑人精品在线| 美女国产高潮福利片在线看| 久久人妻熟女aⅴ| 大片免费播放器 马上看| 十八禁网站免费在线| 国产高清videossex| 少妇精品久久久久久久| 一进一出抽搐动态| 天天躁夜夜躁狠狠躁躁| 久久精品久久久久久噜噜老黄| 性色av一级| 人成视频在线观看免费观看| 欧美国产精品va在线观看不卡| 老司机影院毛片| 啦啦啦在线免费观看视频4| 欧美黄色片欧美黄色片| 青春草亚洲视频在线观看| 日韩视频在线欧美| 老司机午夜福利在线观看视频 | 午夜福利一区二区在线看| 精品欧美一区二区三区在线| 丁香六月天网| av片东京热男人的天堂| 老熟女久久久| 热99re8久久精品国产| 黄色 视频免费看| 国产成人影院久久av| 国产精品一区二区精品视频观看| 美女大奶头黄色视频| 欧美精品啪啪一区二区三区 | 亚洲色图综合在线观看| 91精品三级在线观看| 成人影院久久| 宅男免费午夜| 久久久久久久久免费视频了| 91麻豆精品激情在线观看国产 | 日韩,欧美,国产一区二区三区| 中文字幕最新亚洲高清| 欧美激情 高清一区二区三区| 免费av中文字幕在线| 男女无遮挡免费网站观看| 18禁裸乳无遮挡动漫免费视频| av天堂在线播放| 视频区图区小说| 亚洲国产成人一精品久久久| 国产国语露脸激情在线看| 每晚都被弄得嗷嗷叫到高潮| 大香蕉久久成人网| 动漫黄色视频在线观看| 后天国语完整版免费观看| 欧美成人午夜精品| 青草久久国产| 美女视频免费永久观看网站| 美女主播在线视频| 三上悠亚av全集在线观看| 欧美人与性动交α欧美精品济南到| 日韩中文字幕欧美一区二区| 一区二区三区精品91| 两性午夜刺激爽爽歪歪视频在线观看 | 人妻 亚洲 视频| 精品一区二区三区四区五区乱码| 久久精品国产a三级三级三级| 别揉我奶头~嗯~啊~动态视频 | 十分钟在线观看高清视频www| av天堂在线播放| 亚洲国产av影院在线观看| 亚洲av成人不卡在线观看播放网 | 啦啦啦免费观看视频1| 一本色道久久久久久精品综合| 精品亚洲成国产av| 各种免费的搞黄视频| 新久久久久国产一级毛片| 色婷婷av一区二区三区视频| 涩涩av久久男人的天堂| 亚洲精品乱久久久久久| 美女中出高潮动态图| 2018国产大陆天天弄谢| 欧美变态另类bdsm刘玥| 精品第一国产精品| www.999成人在线观看| 国产成人av激情在线播放| 国产精品二区激情视频| 亚洲 欧美一区二区三区| 99re6热这里在线精品视频| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产精品久久久不卡| 精品国内亚洲2022精品成人 | 韩国精品一区二区三区| 国产麻豆69| 欧美+亚洲+日韩+国产| 人人妻人人澡人人看| 成人国产一区最新在线观看| 免费不卡黄色视频| 亚洲国产毛片av蜜桃av| 亚洲欧洲日产国产| 亚洲国产精品999| 9色porny在线观看| 亚洲精品国产精品久久久不卡| 久久 成人 亚洲| 成人免费观看视频高清| 午夜91福利影院| 国产有黄有色有爽视频| 久久精品亚洲熟妇少妇任你| av网站免费在线观看视频| av免费在线观看网站| 国产黄色免费在线视频| 美女视频免费永久观看网站| 国产精品秋霞免费鲁丝片| 乱人伦中国视频| 精品久久久久久电影网| 啦啦啦免费观看视频1| 美女扒开内裤让男人捅视频| 亚洲精品一二三| 天堂8中文在线网| 99精品欧美一区二区三区四区| 青青草视频在线视频观看| 成年人免费黄色播放视频| 国产免费视频播放在线视频| 大陆偷拍与自拍| 亚洲自偷自拍图片 自拍| 国产成人精品久久二区二区91| 国产精品99久久99久久久不卡| 久久精品成人免费网站| 黄色 视频免费看| 久久久国产成人免费| 韩国精品一区二区三区| av又黄又爽大尺度在线免费看| 999久久久国产精品视频| 久久久久国内视频| 美女视频免费永久观看网站| 国产极品粉嫩免费观看在线| 亚洲精品在线美女| 日韩大码丰满熟妇| 国产高清视频在线播放一区 | 日韩欧美一区视频在线观看| 久久热在线av| 国产精品.久久久| 亚洲av成人一区二区三| 99国产精品一区二区蜜桃av | 女性被躁到高潮视频| 纵有疾风起免费观看全集完整版| 麻豆乱淫一区二区| 女警被强在线播放| 午夜影院在线不卡| 中文字幕最新亚洲高清| 国产三级黄色录像| 两性午夜刺激爽爽歪歪视频在线观看 | 美女福利国产在线| 狠狠婷婷综合久久久久久88av| 免费日韩欧美在线观看| 中国美女看黄片| 超碰97精品在线观看| 国产欧美日韩综合在线一区二区| videos熟女内射| 欧美另类亚洲清纯唯美| 国产精品二区激情视频| 亚洲欧美一区二区三区久久| 母亲3免费完整高清在线观看| 国产精品一二三区在线看| 在线观看免费午夜福利视频| 老司机福利观看| 国产真人三级小视频在线观看| 午夜久久久在线观看| 免费在线观看完整版高清| 黄色毛片三级朝国网站| 狠狠婷婷综合久久久久久88av| 水蜜桃什么品种好| tocl精华| 久久久国产欧美日韩av| 欧美国产精品一级二级三级| 热99久久久久精品小说推荐| 国产1区2区3区精品| 免费观看av网站的网址| 国产日韩欧美视频二区| 精品国产一区二区久久| 巨乳人妻的诱惑在线观看| 久久ye,这里只有精品| 久久久水蜜桃国产精品网| 丝袜美足系列| 交换朋友夫妻互换小说| 久久久精品区二区三区| 日本vs欧美在线观看视频| 久久国产精品男人的天堂亚洲| 国产精品一区二区在线不卡| 亚洲av日韩精品久久久久久密| 久久 成人 亚洲| 国产在线免费精品| 一级毛片精品| 国产成人精品在线电影| 亚洲中文日韩欧美视频| 国产成人精品在线电影| 一区二区三区乱码不卡18| 国产有黄有色有爽视频| 欧美 亚洲 国产 日韩一| 人成视频在线观看免费观看| 久热爱精品视频在线9| 午夜福利在线观看吧| 亚洲欧美日韩另类电影网站| 嫁个100分男人电影在线观看| 国产日韩欧美视频二区| 成人手机av| 精品国产超薄肉色丝袜足j| 这个男人来自地球电影免费观看| 一级a爱视频在线免费观看| a在线观看视频网站| 丰满饥渴人妻一区二区三| 制服诱惑二区| 亚洲国产欧美日韩在线播放| 男人操女人黄网站| 老汉色∧v一级毛片| 亚洲熟女精品中文字幕| 在线观看www视频免费| 又大又爽又粗| 亚洲熟女毛片儿| 免费高清在线观看视频在线观看| 老司机在亚洲福利影院| 日韩欧美免费精品| 夜夜骑夜夜射夜夜干| 日韩欧美免费精品| 午夜福利影视在线免费观看| 免费高清在线观看日韩| netflix在线观看网站| 俄罗斯特黄特色一大片| 一本久久精品| 免费观看av网站的网址| 亚洲五月色婷婷综合| 一级,二级,三级黄色视频| 侵犯人妻中文字幕一二三四区| 亚洲精品久久久久久婷婷小说| 麻豆乱淫一区二区| 大香蕉久久成人网| 999久久久精品免费观看国产| 后天国语完整版免费观看| 成人av一区二区三区在线看 | 婷婷色av中文字幕| 国产熟女午夜一区二区三区| 欧美日韩黄片免| 91国产中文字幕| 亚洲成人免费电影在线观看| 丝袜人妻中文字幕| 水蜜桃什么品种好| 久久天躁狠狠躁夜夜2o2o| 欧美国产精品va在线观看不卡| 叶爱在线成人免费视频播放| av电影中文网址| 国内毛片毛片毛片毛片毛片| 国产在视频线精品| 少妇裸体淫交视频免费看高清 | 男女下面插进去视频免费观看| 欧美 亚洲 国产 日韩一| 黑丝袜美女国产一区| 男男h啪啪无遮挡| 久久国产精品大桥未久av| netflix在线观看网站| 欧美日韩亚洲综合一区二区三区_| 精品卡一卡二卡四卡免费| 成在线人永久免费视频| 午夜91福利影院| 亚洲欧美一区二区三区黑人| 亚洲午夜精品一区,二区,三区| 1024视频免费在线观看| 激情视频va一区二区三区| 69精品国产乱码久久久| 久久人人97超碰香蕉20202| 精品久久久久久电影网| 久久久国产欧美日韩av| 欧美另类亚洲清纯唯美| 日韩欧美一区二区三区在线观看 | 成年动漫av网址| 欧美大码av| 欧美日本中文国产一区发布| 亚洲成人手机| 51午夜福利影视在线观看| 欧美日韩精品网址| 欧美日韩成人在线一区二区| 飞空精品影院首页| 精品国产一区二区久久| 美女中出高潮动态图| 91精品国产国语对白视频| 麻豆av在线久日| 极品少妇高潮喷水抽搐| 国产成人av激情在线播放| 久久99一区二区三区| 午夜精品国产一区二区电影| 男女高潮啪啪啪动态图| 欧美日韩亚洲综合一区二区三区_| 久久久久久免费高清国产稀缺| 可以免费在线观看a视频的电影网站| 国产在线视频一区二区| 九色亚洲精品在线播放| 18禁黄网站禁片午夜丰满| 在线永久观看黄色视频| 国产成+人综合+亚洲专区| 久9热在线精品视频| 久久亚洲精品不卡| 午夜日韩欧美国产| 成人手机av| 国产淫语在线视频| 王馨瑶露胸无遮挡在线观看| 大香蕉久久成人网| 日韩 欧美 亚洲 中文字幕| 91精品国产国语对白视频| 18在线观看网站| 搡老岳熟女国产| 在线观看免费高清a一片| 十八禁人妻一区二区| 国产又爽黄色视频| 热re99久久国产66热| 国产精品成人在线| 精品一品国产午夜福利视频| 亚洲av片天天在线观看| 欧美精品一区二区免费开放| 一区福利在线观看| 欧美日韩国产mv在线观看视频| 青草久久国产| h视频一区二区三区| 人妻人人澡人人爽人人| 99香蕉大伊视频| 亚洲视频免费观看视频| 免费日韩欧美在线观看| 国产成人一区二区三区免费视频网站| 欧美激情 高清一区二区三区| 亚洲一区中文字幕在线| 色婷婷av一区二区三区视频|