• <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)診斷
    中国美白少妇内射xxxbb| 你懂的网址亚洲精品在线观看| 夜夜看夜夜爽夜夜摸| 亚洲精品亚洲一区二区| 亚洲精品乱久久久久久| 国产精品久久久久久av不卡| 亚洲av免费高清在线观看| 干丝袜人妻中文字幕| av专区在线播放| 久久久久国产网址| 亚洲精品成人久久久久久| 日韩三级伦理在线观看| 久久久久久国产a免费观看| 一本久久精品| videos熟女内射| 少妇熟女aⅴ在线视频| 婷婷色麻豆天堂久久| 日韩三级伦理在线观看| 毛片女人毛片| 亚洲国产成人一精品久久久| 不卡视频在线观看欧美| 国产黄片美女视频| 免费少妇av软件| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av天美| av在线老鸭窝| 精品久久久噜噜| 欧美一级a爱片免费观看看| 亚洲国产日韩欧美精品在线观看| 在线播放无遮挡| 亚洲最大成人中文| 国产白丝娇喘喷水9色精品| 欧美日韩综合久久久久久| 2022亚洲国产成人精品| 亚洲成人中文字幕在线播放| 91在线精品国自产拍蜜月| 非洲黑人性xxxx精品又粗又长| 精品久久久噜噜| 肉色欧美久久久久久久蜜桃 | 自拍偷自拍亚洲精品老妇| 国产男人的电影天堂91| 亚洲av免费高清在线观看| 天堂影院成人在线观看| 九九在线视频观看精品| 黑人高潮一二区| 嫩草影院入口| 亚洲国产成人一精品久久久| 伊人久久精品亚洲午夜| 日本午夜av视频| 激情 狠狠 欧美| 国产美女午夜福利| 午夜视频国产福利| 午夜福利高清视频| 亚洲最大成人av| 国产老妇伦熟女老妇高清| 国产在视频线在精品| 看十八女毛片水多多多| 日韩制服骚丝袜av| 最后的刺客免费高清国语| 亚洲图色成人| 欧美区成人在线视频| 男女下面进入的视频免费午夜| 国产高清国产精品国产三级 | 美女内射精品一级片tv| 国产乱人偷精品视频| 岛国毛片在线播放| 亚洲怡红院男人天堂| 亚洲欧美日韩卡通动漫| 男人狂女人下面高潮的视频| 搞女人的毛片| 26uuu在线亚洲综合色| 国产精品女同一区二区软件| 国产v大片淫在线免费观看| 久久亚洲国产成人精品v| 91午夜精品亚洲一区二区三区| 亚洲精品色激情综合| videossex国产| 熟妇人妻不卡中文字幕| 三级经典国产精品| 成人国产麻豆网| 两个人视频免费观看高清| 91精品国产九色| 久久亚洲国产成人精品v| 2022亚洲国产成人精品| av在线播放精品| 亚洲国产成人一精品久久久| 欧美区成人在线视频| 亚洲av男天堂| 亚洲av免费高清在线观看| 国产成人a区在线观看| 在线天堂最新版资源| 免费av毛片视频| 美女内射精品一级片tv| 国产成人精品一,二区| 夜夜看夜夜爽夜夜摸| 日韩欧美三级三区| 精品人妻熟女av久视频| 18+在线观看网站| 日韩av在线免费看完整版不卡| 亚洲国产欧美人成| 色视频www国产| 在线播放无遮挡| 蜜臀久久99精品久久宅男| 综合色av麻豆| 国产精品蜜桃在线观看| av福利片在线观看| 久久人人爽人人爽人人片va| 日韩一区二区三区影片| 国产毛片a区久久久久| eeuss影院久久| 一区二区三区四区激情视频| 国产亚洲av片在线观看秒播厂 | 亚洲av一区综合| 三级毛片av免费| 免费看不卡的av| 亚洲成人精品中文字幕电影| 精品久久久精品久久久| 成年女人在线观看亚洲视频 | 亚洲精品国产成人久久av| 99热全是精品| 国产黄片美女视频| 国产成人福利小说| 成人二区视频| 成人亚洲精品一区在线观看 | 日本-黄色视频高清免费观看| 国产欧美另类精品又又久久亚洲欧美| 久久精品夜色国产| 免费大片黄手机在线观看| 久久精品夜夜夜夜夜久久蜜豆| 欧美日韩国产mv在线观看视频 | 极品少妇高潮喷水抽搐| 国产精品综合久久久久久久免费| 在线观看美女被高潮喷水网站| 久久韩国三级中文字幕| 亚洲人成网站高清观看| 欧美精品一区二区大全| 综合色丁香网| 中文字幕免费在线视频6| 亚洲av二区三区四区| 精品一区二区免费观看| 国产精品一区www在线观看| 欧美日本视频| 亚洲伊人久久精品综合| 熟女人妻精品中文字幕| 九草在线视频观看| 精品久久久久久久久av| 九九爱精品视频在线观看| 国产老妇女一区| 日日撸夜夜添| 少妇人妻精品综合一区二区| 国产一区二区三区综合在线观看 | 亚洲欧美日韩无卡精品| 性色avwww在线观看| 精品国内亚洲2022精品成人| 亚洲乱码一区二区免费版| 日韩成人av中文字幕在线观看| 亚洲精品中文字幕在线视频 | 国产高清国产精品国产三级 | 青春草国产在线视频| 国产成人一区二区在线| 亚洲av成人精品一二三区| 国模一区二区三区四区视频| 免费观看精品视频网站| 精华霜和精华液先用哪个| 亚洲精品影视一区二区三区av| 极品少妇高潮喷水抽搐| 中文字幕av成人在线电影| 亚洲av电影在线观看一区二区三区 | 插阴视频在线观看视频| 国产在线男女| 亚洲真实伦在线观看| 亚洲av一区综合| 日本三级黄在线观看| 亚洲欧美一区二区三区黑人 | 不卡视频在线观看欧美| 午夜激情久久久久久久| 如何舔出高潮| 亚洲av日韩在线播放| 国产视频首页在线观看| 亚洲精品一区蜜桃| 国产亚洲5aaaaa淫片| 免费观看性生交大片5| 高清毛片免费看| 日本-黄色视频高清免费观看| 男女国产视频网站| av线在线观看网站| 91午夜精品亚洲一区二区三区| 亚洲av一区综合| 亚洲国产精品国产精品| 亚洲av中文av极速乱| 国产免费视频播放在线视频 | 欧美日韩视频高清一区二区三区二| 午夜精品国产一区二区电影 | 久久99精品国语久久久| 一区二区三区四区激情视频| 可以在线观看毛片的网站| 蜜臀久久99精品久久宅男| 在线观看免费高清a一片| 免费黄频网站在线观看国产| 美女大奶头视频| 欧美成人一区二区免费高清观看| 亚洲精华国产精华液的使用体验| 亚洲国产精品专区欧美| 国产乱来视频区| 色综合色国产| 久久久久久国产a免费观看| 国产精品伦人一区二区| 国产女主播在线喷水免费视频网站 | 亚洲天堂国产精品一区在线| av播播在线观看一区| 少妇猛男粗大的猛烈进出视频 | 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 99久国产av精品| 熟女电影av网| 女人十人毛片免费观看3o分钟| 老师上课跳d突然被开到最大视频| 久久热精品热| 亚洲精品影视一区二区三区av| 啦啦啦韩国在线观看视频| 自拍偷自拍亚洲精品老妇| 成人无遮挡网站| 赤兔流量卡办理| 日韩中字成人| 亚洲av免费高清在线观看| 人妻一区二区av| 久久久成人免费电影| 乱系列少妇在线播放| 日本熟妇午夜| 精品久久久久久久久av| 天堂影院成人在线观看| 国产亚洲5aaaaa淫片| av国产久精品久网站免费入址| 成人一区二区视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 非洲黑人性xxxx精品又粗又长| 网址你懂的国产日韩在线| 一级片'在线观看视频| 亚洲一区高清亚洲精品| 人体艺术视频欧美日本| 国产高清国产精品国产三级 | 国产精品伦人一区二区| 久久热精品热| 禁无遮挡网站| 黄色一级大片看看| 91久久精品国产一区二区三区| 亚洲人成网站在线播| 欧美精品一区二区大全| 天美传媒精品一区二区| 国产激情偷乱视频一区二区| 国产成人免费观看mmmm| 乱系列少妇在线播放| 精品少妇黑人巨大在线播放| 有码 亚洲区| 中文字幕亚洲精品专区| 亚洲真实伦在线观看| 久久久久九九精品影院| 亚洲在久久综合| 中文精品一卡2卡3卡4更新| 欧美丝袜亚洲另类| h日本视频在线播放| 国产亚洲午夜精品一区二区久久 | 黄色一级大片看看| 欧美性感艳星| 青春草视频在线免费观看| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产| 国产黄片美女视频| 校园人妻丝袜中文字幕| 成年女人在线观看亚洲视频 | 亚洲精品影视一区二区三区av| 午夜免费观看性视频| 久久精品国产自在天天线| 国产美女午夜福利| 亚洲欧美一区二区三区黑人 | 国产女主播在线喷水免费视频网站 | 美女被艹到高潮喷水动态| 欧美bdsm另类| 青春草国产在线视频| 免费观看av网站的网址| 久久韩国三级中文字幕| 伦精品一区二区三区| 日韩欧美三级三区| 精品久久久久久久久久久久久| 只有这里有精品99| 精品欧美国产一区二区三| 看免费成人av毛片| 亚洲av一区综合| 嘟嘟电影网在线观看| 午夜福利网站1000一区二区三区| 久久午夜福利片| 人妻一区二区av| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影小说 | 国产麻豆成人av免费视频| 插逼视频在线观看| 日韩av在线大香蕉| 亚洲最大成人中文| 青春草亚洲视频在线观看| 丝袜喷水一区| 国产成人91sexporn| 欧美丝袜亚洲另类| 日韩精品青青久久久久久| 免费观看的影片在线观看| 伊人久久精品亚洲午夜| 在线观看美女被高潮喷水网站| 一级毛片aaaaaa免费看小| 搡女人真爽免费视频火全软件| 99热6这里只有精品| 午夜福利在线观看免费完整高清在| 日本爱情动作片www.在线观看| 成年女人在线观看亚洲视频 | 97热精品久久久久久| 国产伦一二天堂av在线观看| 乱人视频在线观看| 亚洲在线观看片| 日韩精品有码人妻一区| 极品教师在线视频| 日韩人妻高清精品专区| 一本一本综合久久| 国产av码专区亚洲av| 嫩草影院新地址| 午夜视频国产福利| 亚洲欧洲日产国产| 国产免费又黄又爽又色| 一级毛片 在线播放| 蜜桃亚洲精品一区二区三区| 国产av在哪里看| 成人性生交大片免费视频hd| 成年女人看的毛片在线观看| 91精品伊人久久大香线蕉| 免费av观看视频| 亚洲av福利一区| 内地一区二区视频在线| 男人舔女人下体高潮全视频| 国产单亲对白刺激| 免费高清在线观看视频在线观看| 人妻少妇偷人精品九色| 全区人妻精品视频| 日韩欧美精品免费久久| 女人久久www免费人成看片| 亚洲欧洲日产国产| 欧美97在线视频| 少妇丰满av| 精品一区二区三区视频在线| 国产熟女欧美一区二区| 一区二区三区免费毛片| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 精品少妇黑人巨大在线播放| 一个人免费在线观看电影| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 精品亚洲乱码少妇综合久久| 最新中文字幕久久久久| 中文字幕免费在线视频6| 免费观看a级毛片全部| 建设人人有责人人尽责人人享有的 | 日本黄色片子视频| 亚洲精品亚洲一区二区| 亚洲精品第二区| 亚洲av日韩在线播放| 在线观看免费高清a一片| 国产亚洲av片在线观看秒播厂 | 最近2019中文字幕mv第一页| 纵有疾风起免费观看全集完整版 | 亚洲欧美日韩东京热| 欧美一区二区亚洲| 中国美白少妇内射xxxbb| 男女国产视频网站| 久久久欧美国产精品| 2018国产大陆天天弄谢| 亚洲三级黄色毛片| 国产男人的电影天堂91| 欧美一级a爱片免费观看看| 亚洲一级一片aⅴ在线观看| 蜜桃久久精品国产亚洲av| 91av网一区二区| a级毛色黄片| 日产精品乱码卡一卡2卡三| 成人av在线播放网站| 亚洲第一区二区三区不卡| av国产久精品久网站免费入址| 大香蕉久久网| 日本av手机在线免费观看| 特级一级黄色大片| 深爱激情五月婷婷| 日日摸夜夜添夜夜添av毛片| 国产亚洲午夜精品一区二区久久 | av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 女人久久www免费人成看片| 日本爱情动作片www.在线观看| 伦精品一区二区三区| 激情 狠狠 欧美| 波多野结衣巨乳人妻| 爱豆传媒免费全集在线观看| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区久久| 午夜精品一区二区三区免费看| 中文精品一卡2卡3卡4更新| 黄色欧美视频在线观看| 亚洲欧洲日产国产| 日韩在线高清观看一区二区三区| 免费看美女性在线毛片视频| 特级一级黄色大片| av在线播放精品| 日韩,欧美,国产一区二区三区| 日本爱情动作片www.在线观看| 国产一区二区三区综合在线观看 | 日韩大片免费观看网站| 99热这里只有是精品50| 久久久国产一区二区| 免费少妇av软件| 国产成人精品婷婷| 成年女人在线观看亚洲视频 | 51国产日韩欧美| 国产成人精品久久久久久| 最近手机中文字幕大全| 国产免费一级a男人的天堂| 午夜精品在线福利| 亚洲欧美一区二区三区黑人 | 2021天堂中文幕一二区在线观| 国产黄色视频一区二区在线观看| 免费高清在线观看视频在线观看| 成人二区视频| 日韩不卡一区二区三区视频在线| 丰满乱子伦码专区| 91精品国产九色| 国产乱人视频| 国产av国产精品国产| 天堂网av新在线| 欧美bdsm另类| 午夜激情欧美在线| 亚洲电影在线观看av| 天美传媒精品一区二区| 免费观看精品视频网站| 婷婷六月久久综合丁香| 好男人在线观看高清免费视频| 成人高潮视频无遮挡免费网站| 国产精品三级大全| 免费看日本二区| 免费观看在线日韩| 搡老妇女老女人老熟妇| 国产探花在线观看一区二区| 国内揄拍国产精品人妻在线| 久久亚洲国产成人精品v| 成人鲁丝片一二三区免费| 久久精品国产亚洲av涩爱| 国产成年人精品一区二区| 少妇裸体淫交视频免费看高清| 欧美高清成人免费视频www| 超碰av人人做人人爽久久| 又黄又爽又刺激的免费视频.| 日韩三级伦理在线观看| 国产在视频线在精品| 少妇的逼好多水| 乱码一卡2卡4卡精品| 中文精品一卡2卡3卡4更新| 国产在视频线在精品| 国产国拍精品亚洲av在线观看| 国产在视频线在精品| 亚洲欧洲国产日韩| 男女国产视频网站| 少妇的逼水好多| 青青草视频在线视频观看| 十八禁国产超污无遮挡网站| 久久久久久伊人网av| 成人欧美大片| 国产一区亚洲一区在线观看| 亚洲精品一二三| 国产成人精品婷婷| 久久久久免费精品人妻一区二区| 激情五月婷婷亚洲| 不卡视频在线观看欧美| 日本av手机在线免费观看| 国产老妇女一区| 一二三四中文在线观看免费高清| av福利片在线观看| 午夜福利视频精品| 国产高潮美女av| 99久国产av精品国产电影| 国产极品天堂在线| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 国产老妇伦熟女老妇高清| 白带黄色成豆腐渣| 国产亚洲av嫩草精品影院| 国产老妇女一区| 国产麻豆成人av免费视频| 亚洲第一区二区三区不卡| 久久精品国产亚洲网站| 神马国产精品三级电影在线观看| 久久久久久久大尺度免费视频| 久久99精品国语久久久| 亚洲精品一区蜜桃| 色综合亚洲欧美另类图片| 国产伦精品一区二区三区四那| 久久久久久久久久黄片| 中国美白少妇内射xxxbb| av一本久久久久| 伦精品一区二区三区| 最近视频中文字幕2019在线8| 国产黄片美女视频| 午夜福利视频精品| 国产爱豆传媒在线观看| 人妻少妇偷人精品九色| 美女主播在线视频| 久久草成人影院| 日韩三级伦理在线观看| 国产日韩欧美在线精品| 国产综合精华液| 免费看不卡的av| 亚洲高清免费不卡视频| 国产成人精品福利久久| 男人爽女人下面视频在线观看| 亚洲欧美精品自产自拍| 超碰97精品在线观看| av播播在线观看一区| 国产综合懂色| 99久久中文字幕三级久久日本| 国产又色又爽无遮挡免| 国产一区二区亚洲精品在线观看| 伊人久久精品亚洲午夜| 校园人妻丝袜中文字幕| 身体一侧抽搐| 日本欧美国产在线视频| 亚洲成人中文字幕在线播放| 久久久久久久大尺度免费视频| 亚洲18禁久久av| 国产日韩欧美在线精品| 免费观看性生交大片5| 内射极品少妇av片p| 能在线免费看毛片的网站| 亚洲成色77777| 狂野欧美白嫩少妇大欣赏| 边亲边吃奶的免费视频| 亚洲人成网站在线播| 亚洲av成人精品一二三区| 精品99又大又爽又粗少妇毛片| 精华霜和精华液先用哪个| 又爽又黄无遮挡网站| 日本午夜av视频| 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 黄色配什么色好看| 成人午夜高清在线视频| 亚洲精品aⅴ在线观看| 久久久久网色| av国产免费在线观看| 韩国高清视频一区二区三区| 亚洲乱码一区二区免费版| 久久精品国产亚洲网站| 国产精品国产三级国产av玫瑰| 精品久久久久久久久久久久久| 中国国产av一级| 欧美日韩一区二区视频在线观看视频在线 | 亚洲美女视频黄频| 国产成人午夜福利电影在线观看| 最后的刺客免费高清国语| 自拍偷自拍亚洲精品老妇| 国产精品麻豆人妻色哟哟久久 | 肉色欧美久久久久久久蜜桃 | 黄色一级大片看看| 午夜精品一区二区三区免费看| 一个人免费在线观看电影| 我的老师免费观看完整版| 亚洲人成网站在线播| 免费av毛片视频| 亚洲精品乱码久久久v下载方式| 成人特级av手机在线观看| 精品不卡国产一区二区三区| 三级国产精品片| 免费电影在线观看免费观看| 国产午夜精品论理片| 日韩亚洲欧美综合| 婷婷色综合大香蕉| 午夜视频国产福利| 插阴视频在线观看视频| 亚洲熟女精品中文字幕| 欧美97在线视频| av卡一久久| 91久久精品电影网| 久久精品国产亚洲av天美| av在线天堂中文字幕| 日本黄大片高清| 国产精品蜜桃在线观看| 97在线视频观看| av福利片在线观看| 成人鲁丝片一二三区免费| 少妇的逼好多水| 国产男女超爽视频在线观看| 人妻系列 视频| 777米奇影视久久| 三级经典国产精品| 天堂俺去俺来也www色官网 | 国产精品一二三区在线看| 国产 亚洲一区二区三区 | 97超碰精品成人国产| 国语对白做爰xxxⅹ性视频网站| 国产真实伦视频高清在线观看| 日韩视频在线欧美| 美女内射精品一级片tv| 亚洲自拍偷在线| 日韩一区二区视频免费看| 亚洲高清免费不卡视频| 黄片无遮挡物在线观看| 国产色婷婷99| 在线观看免费高清a一片| 自拍偷自拍亚洲精品老妇| 联通29元200g的流量卡| 国模一区二区三区四区视频| 91av网一区二区|