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

    基于氮收支平衡的河套灌區(qū)春小麥農(nóng)田灌溉和施氮策略

    2023-01-16 09:42:12李根東陳志君張雪晨黃冠華

    李 越,李根東,陳志君,張雪晨,黃冠華※

    ·農(nóng)業(yè)水土工程·

    基于氮收支平衡的河套灌區(qū)春小麥農(nóng)田灌溉和施氮策略

    李 越1,2,李根東3,陳志君1,4,張雪晨1,2,黃冠華1,2※

    (1. 中國(guó)農(nóng)業(yè)大學(xué),中國(guó)-以色列國(guó)際農(nóng)業(yè)研究培訓(xùn)中心,北京 100083;2. 中國(guó)農(nóng)業(yè)大學(xué),中國(guó)農(nóng)業(yè)水問(wèn)題研究中心,北京 100083;3.內(nèi)蒙古河套灌區(qū)水利發(fā)展中心,巴彥淖爾 015000;4. 沈陽(yáng)農(nóng)業(yè)大學(xué)水利學(xué)院,沈陽(yáng) 110866)

    針對(duì)中國(guó)黃河中上游河套灌區(qū)不合理灌溉和施肥造成的土壤氮素流失嚴(yán)重及氮收支不平衡等問(wèn)題,該研究于2019-2021年開(kāi)展田間試驗(yàn),探討不同灌溉和施肥策略對(duì)土壤氮損失、作物氮吸收及氮收支的影響。試驗(yàn)設(shè)置了3個(gè)灌溉水平(高水I1:450 mm,中水I2:315 mm,低水I3:180 mm)和2個(gè)施氮水平(高氮N1:340 kg/hm2,低氮N3:170 kg/hm2),此外,2020和2021年在中等灌溉水平I2下補(bǔ)充了中等施氮水平(250 kg/hm2,N2),對(duì)不同處理的土壤氮損失、作物氮吸收及氮收支等指標(biāo)進(jìn)行了對(duì)比分析。結(jié)果表明,肥料氮是農(nóng)田氮輸入的主要來(lái)源,其次是灌溉水、大氣沉降和非豆科作物固定。作物吸氮占土壤氮輸出的比例最大,其次是NO3--N淋失、NH3揮發(fā)和N2O排放。對(duì)于氮輸入而言,其值隨著灌水量和施氮量的減少而降低。對(duì)于土壤氮輸出而言,減少灌水量和施氮量可顯著降低土壤總氮損失量,但過(guò)低的灌水量和施氮量將導(dǎo)致小麥吸氮量的降低。傳統(tǒng)的N1施氮處理可導(dǎo)致土壤氮素盈余,而施氮量降低50%的N3處理則導(dǎo)致土壤氮素大量虧缺。對(duì)照處理(I1N1)的土壤氮損失量最高,該處理氮損失占土壤氮輸出的比例高達(dá)25%~35%,其中NO3--N淋失和NH3揮發(fā)占總氮損失的95%以上。與對(duì)照處理相比,I2N2處理可減少21%~29%的氮損失,且作物吸氮幾乎未受到影響。同時(shí),該處理土壤氮素處于輕度虧缺狀態(tài),其虧缺量為28~50 kg/hm2,占總施氮量的11%~20%。若在收獲后將4~8 t/hm2的小麥秸稈還田,則可保持麥田土壤的氮收支平衡。因此,通過(guò)改善灌溉和施肥策略并配合適當(dāng)?shù)慕斩掃€田可以有效緩解河套灌區(qū)春小麥農(nóng)田的氮損失且實(shí)現(xiàn)土壤氮平衡,該研究可為干旱半干旱地區(qū)春小麥農(nóng)田可持續(xù)生產(chǎn)和氮污染物減排提供科學(xué)依據(jù)。

    氮;地下水;溫室氣體;水肥管理;排放;污染物;環(huán)境

    0 引 言

    全球范圍內(nèi),灌溉和施肥已經(jīng)成為干旱半干旱地區(qū)糧食生產(chǎn)的重要農(nóng)田管理措施。在過(guò)去的幾十年里,大量的灌溉水和氮肥被用于農(nóng)業(yè)生產(chǎn),盡管在很大程度上提高了作物產(chǎn)量,但也產(chǎn)生了大量的氮污染物,如NH3、N2O和NO3--N等[1-2]。這些氮污染物的排放帶來(lái)了溫室效應(yīng)、臭氧層破壞、地下水硝酸鹽污染等一系列不利的環(huán)境和氣候問(wèn)題[3]。然而,如果氮肥供應(yīng)量過(guò)低則不能滿足作物對(duì)氮素的需求,從而導(dǎo)致土壤氮虧缺和土壤肥力下降,最終導(dǎo)致農(nóng)田系統(tǒng)功能退化[4]。因此,平衡農(nóng)田系統(tǒng)中土壤氮收支對(duì)于降低環(huán)境污染風(fēng)險(xiǎn)、維持作物產(chǎn)量和土壤肥力至關(guān)重要。

    灌溉可以增加農(nóng)田系統(tǒng)的氮輸入,但它對(duì)不同氮輸出途徑的影響不同[4]。許多研究發(fā)現(xiàn),過(guò)度灌溉會(huì)導(dǎo)致較高的NO3--N淋失,而減少灌水量是降低土壤NO3--N淋失的有效途徑[5]。與土壤NO3--N淋失相反,由于NH3分子的親水性,其揮發(fā)量通常隨著灌水量的減少而增加。例如,李銀坤等[6]研究表明,減少灌水量加劇了農(nóng)田土壤NH3揮發(fā)。此外,灌溉對(duì)N2O排放和作物氮吸收的影響在不同條件下有所不同[7]。雷宏軍等[8]研究發(fā)現(xiàn),在中國(guó)北方干旱半干旱地區(qū),土壤N2O排放隨著灌水量的增加而增加。但Lam等[7]研究發(fā)現(xiàn),在半干旱地區(qū),補(bǔ)充灌溉反而顯著降低了土壤N2O排放。對(duì)于作物氮吸收而言,Djaman等[9]研究發(fā)現(xiàn),在0~180 mm的灌水量范圍內(nèi),作物吸氮量隨著灌水量的增加而增加。但楊曉亞等[10]的研究表明,作物吸氮量在灌溉120和180 mm之間無(wú)顯著差異。可以看出,灌溉對(duì)農(nóng)田土壤氮輸出的綜合影響尚不完全清楚。

    氮肥是農(nóng)田土壤氮輸入的主要途徑,同時(shí)也與土壤各項(xiàng)氮輸出密切相關(guān)。在合理的范圍內(nèi)增加施氮量可以通過(guò)提高土壤速效氮含量進(jìn)而促進(jìn)作物對(duì)氮素的吸收[11]。例如,Shi等[12]研究發(fā)現(xiàn),在小麥農(nóng)田系統(tǒng)中,將施氮量從0增加至300 kg/hm2可以顯著提高作物的氮素吸收量。然而,過(guò)度的氮輸入對(duì)于作物吸氮量增加效果不顯著,反而會(huì)導(dǎo)致大量氮素在土壤中殘留,這增加了土壤氮損失的風(fēng)險(xiǎn)[13]。Zhang等[4]的研究發(fā)現(xiàn),在小麥農(nóng)田系統(tǒng)中,與施氮420 kg/hm2相比,施氮600 kg/hm2使土壤氣態(tài)氮損失(NH3和N2O)顯著增加了約43%。Liu等[13]的研究發(fā)現(xiàn),在小麥-玉米種植系統(tǒng)中,土壤NO3--N淋失量與氮肥投入量呈正相關(guān)關(guān)系,將施氮量從240 kg/hm2增加到360 kg/hm2使作物生育期內(nèi)的NO3--N淋失量提高了約42%。

    氣候干旱的河套灌區(qū)是中國(guó)重要的優(yōu)質(zhì)春小麥生產(chǎn)基地之一,在保障中國(guó)糧食安全方面發(fā)揮著重要作用。在該地區(qū),農(nóng)民通常進(jìn)行過(guò)度灌溉和施肥來(lái)保證小麥高產(chǎn)。相關(guān)研究表明,春小麥生育期內(nèi)的灌水量可達(dá)300 mm,部分地區(qū)甚至超過(guò)了450 mm,并且施氮量高達(dá)350 kg/hm2[14]。過(guò)度的氮肥投入無(wú)疑會(huì)導(dǎo)致土壤氮素大量盈余,顯著提高了NH3揮發(fā)、N2O排放和NO3--N淋失的風(fēng)險(xiǎn)。在大量灌溉水的驅(qū)動(dòng)下,超過(guò)18%的肥料氮以NH3揮發(fā)和N2O排放的方式進(jìn)入大氣[15],約30%的肥料氮以NO3--N的形式被淋洗進(jìn)入地下水[16]。如此長(zhǎng)期大量的氮污染物排放嚴(yán)重威脅到中國(guó)黃河中上游地區(qū)的地下水和大氣環(huán)境。因此,在保持作物氮素吸收和土壤肥力的前提下,改善該地區(qū)傳統(tǒng)的灌溉和施肥策略,降低農(nóng)田系統(tǒng)的氮損失勢(shì)在必行。然而,關(guān)于中國(guó)黃河中上游河套灌區(qū)不同灌溉和施肥條件下農(nóng)田氮損失及土壤氮收支的研究仍較為罕見(jiàn)。因此,為解決農(nóng)田土壤氮流失造成的一系列環(huán)境和土壤問(wèn)題,本研究探討灌溉施肥對(duì)農(nóng)田土壤氮收支的影響,在此基礎(chǔ)上考慮土壤氮損失減排、作物氮吸收和土壤肥力保持,提出適宜該地區(qū)春小麥可持續(xù)生產(chǎn)的灌溉和施肥策略,以期為試驗(yàn)區(qū)春小麥農(nóng)田可持續(xù)生產(chǎn)和氮污染物減排提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    本研究于2019—2021年在內(nèi)蒙古自治區(qū)巴彥淖爾市中國(guó)農(nóng)業(yè)大學(xué)河套實(shí)驗(yàn)站(41°09′N;107°39′E)開(kāi)展,實(shí)驗(yàn)站具體位置見(jiàn)圖1。該地區(qū)平均海拔1 031 m,屬半干旱溫帶大陸性氣候,日照充足,晝夜溫差大,多年平均氣溫6.8 ℃。降水量少,蒸發(fā)量大,地下水埋深淺。多年平均降水量?jī)H160~240 mm,蒸發(fā)量可達(dá)2 200~2 400 mm,作物生長(zhǎng)季節(jié)內(nèi)的地下水埋深一般在1.2~2.7 m之間變化[17]。2019—2021年春小麥生育期內(nèi)逐日降水量和平均氣溫見(jiàn)圖2。試驗(yàn)地0~70 cm土壤為粉質(zhì)壤土,70~100 cm為壤土。0~100 cm土壤平均有機(jī)質(zhì)含量為14.33 g/kg,平均容重為1.44 g/cm3,詳細(xì)的土壤理化性質(zhì)見(jiàn)表1。

    注:I1,高水450 mm;I2,中水315 mm;I3,低水180 mm。N1,高氮340 kg·hm-2;N2,中氮250 kg·hm-2;N3,低氮170 kg·hm-2。下同。

    圖2 2019—2021年春小麥生育期內(nèi)逐日降水量和平均空氣溫度

    表1 試驗(yàn)田土壤的主要物理和化學(xué)性狀指標(biāo)

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

    田間試驗(yàn)的春小麥種植品種為永良4號(hào),2019和2020年在3月15日播種,7月20日收獲,由于耕作原因,2021年在4月6日播種,7月20日收獲。種植行距15 cm,每穴間隔5 cm,種植密度6.03×106株/hm2。在該地區(qū),春小麥生育期內(nèi)灌水量最高可達(dá)450 mm,施氮量可達(dá)340 kg/hm2[5,14]。因此,本研究以當(dāng)?shù)毓喔群褪┑繛閷?duì)照,參考前人推薦的灌水和施氮量[5,16],設(shè)置了不同的節(jié)水和減氮處理。具體如下:試驗(yàn)設(shè)置3個(gè)灌溉水平(高水I1:450 mm,中水I2:315 mm和低水I3:180 mm)和2個(gè)施氮水平(高氮N1:340 kg/hm2和低氮N3:170 kg/hm2)。共6種灌溉施肥組合:I1N1,I2N1,I3N1,I1N3,I2N3和I3N3。此外,2020和2021年在中等灌溉水平I2下還考慮了中等施氮水平N2:250 kg/hm2。每個(gè)處理3個(gè)重復(fù),2019年共18個(gè)小區(qū),2020和2021年共21個(gè)小區(qū),每個(gè)小區(qū)面積175 m2,長(zhǎng)×寬為25 m×7 m。試驗(yàn)田四周預(yù)留寬10 m的保護(hù)行以避免當(dāng)?shù)剞r(nóng)民灌溉和施肥的影響。此外,每個(gè)小區(qū)四周設(shè)置田埂高×寬為0.4 m× 0.6 m,以盡量減少相鄰小區(qū)的水平干擾。在所有處理中,I1N1處理為對(duì)照,與當(dāng)?shù)剞r(nóng)民傳統(tǒng)的灌溉和施肥措施一致。

    灌溉方式為畦灌,每個(gè)處理的灌水量用精度為0.001 m3的LXS-25E型旋翼式水表(寧波甬升水表實(shí)業(yè)有限公司)控制。2019—2021年均在分蘗期,拔節(jié)期和抽穗期進(jìn)行3次等量灌溉(2019和2020年灌溉時(shí)間依次為5月10日,5月21日和6月18日,2021年灌溉時(shí)間依次為5月17日,5月28日和6月25日)。對(duì)于I1,I2和I3灌溉處理,每次灌溉的灌水定額分別為150、105和60 mm。為了確保每個(gè)小區(qū)的灌溉均勻性,在每季小麥播種前采用激光平地機(jī)對(duì)試驗(yàn)田土地進(jìn)行平整,同時(shí),在灌水過(guò)程中,每個(gè)小區(qū)縱向均勻布置了3個(gè)出水口。每季的基肥隨小麥播種同時(shí)施入,各處理基肥為相同用量的磷肥(以P2O5計(jì),95 kg/hm2),鉀肥(以K2O計(jì),95 kg/hm2)和氮肥(以N計(jì),50 kg/hm2)。根據(jù)每個(gè)處理設(shè)計(jì)的總施氮量,施入基肥后剩余氮肥的60%在小麥分蘗期追施,40%在拔節(jié)期追施。在這2個(gè)生育期灌溉前將肥料均勻撒施在土壤表面(2019和2020年追肥時(shí)間為5月10日和5月21日,2021年追肥時(shí)間為5月17日和5月28日)。磷肥、鉀肥和氮肥選用當(dāng)?shù)爻S玫牧姿岫@、硫酸鉀和尿素。肥料類型、施肥方式、施肥時(shí)間與當(dāng)?shù)卮盒←溕a(chǎn)一致。具體的灌溉和施肥制度見(jiàn)表2。

    表2 春小麥農(nóng)田的灌溉和施肥制度

    1.3 樣品采集及確定

    土壤氮素淋失的測(cè)量方法與Li等[18]相同。在每個(gè)小區(qū)中間埋設(shè)一個(gè)高100 cm,直徑40 cm的圓柱體和40 cm高的錐形容器組成的鋼桶(滲漏桶),與土壤表面相平。100 cm高的圓柱體內(nèi)填入與田間相同容重的土壤,底部的錐形容器清空用于收集滲濾液。每次灌溉和降雨后48 h內(nèi)用真空泵抽取滲濾液,直至完全清空,并據(jù)此測(cè)量滲濾液的體積。同時(shí),使用紫外分光光度計(jì)(UV-3100,中國(guó))測(cè)量滲濾液中硝態(tài)氮的濃度。滲漏桶中收集的土壤溶液僅用于定量描述隨著灌溉和降雨而引起的氮素淋失量,暫未考慮地下水的毛管補(bǔ)償作用。

    使用通風(fēng)法測(cè)量每個(gè)小區(qū)的土壤NH3揮發(fā)量[19]。將直徑16 cm,高25 cm的PMMA采樣管插入土壤5 cm深,并將兩層浸泡過(guò)磷酸甘油的海綿置于通氣管內(nèi)。上層海綿用于消除大氣中NH3的影響,下層海綿用于吸收土壤中揮發(fā)的NH3。在試驗(yàn)過(guò)程中平均每10 d收集和更換一次海綿。灌溉和施肥后,采樣間隔縮短至2 d。采樣后立即將含有NH4+-N樣品的海綿帶到實(shí)驗(yàn)室并浸入500 mL濃度為1.0 mol/L的KCl溶液中。將浸泡的樣品海綿和溶液用體積為1 L聚乙烯瓶密封并在往復(fù)式震蕩器上以200 rpm/min的速度振動(dòng)1 h。使用紫外分光光度計(jì)測(cè)量每個(gè)樣品提取液中的NH4+-N濃度。

    使用靜態(tài)箱法收集每個(gè)小區(qū)土壤排放的N2O。靜態(tài)箱系統(tǒng)由一個(gè)固定在土壤中20 cm深的底座(50 cm× 50 cm×25 cm)和一個(gè)可移動(dòng)的頂箱(50 cm×50 cm× 100 cm)組成。底座為不銹鋼材質(zhì),頂箱為聚甲基丙烯酸甲酯材質(zhì)。靜態(tài)箱底座和頂箱連接處采用水封的方式保證氣密性。在頂箱上部安裝一個(gè)小風(fēng)扇以混合靜態(tài)箱內(nèi)的氣體。靜態(tài)箱內(nèi)安裝溫度計(jì)用于測(cè)量?jī)?nèi)部空氣溫度。每次從09:00至10:00,使用50 mL聚丙烯注射器從靜態(tài)箱內(nèi)連續(xù)收集4個(gè)氣體樣品,采樣間隔為15 min,并將氣體儲(chǔ)存在氣袋中。使用氣相色譜儀(GC-2014C,日本)在48 h內(nèi)對(duì)氣體樣品進(jìn)行分析。

    在小麥?zhǔn)斋@后,在0~100 cm土層采集土壤樣品用于測(cè)定土壤NO3--N含量。0~20 cm土層的采樣間隔為10 cm,>20~100 cm土層的采樣間隔為20 cm。將采集的土樣在2 mol/L的KCl溶液里浸提,土壤與KCl溶液的比例為1∶5。然后用紫外分光光度計(jì)測(cè)量浸提液中的NO3--N濃度。詳細(xì)的樣品制備和計(jì)算過(guò)程在Li等[16]的研究中進(jìn)行了描述。

    作物氮素吸收量采用干燒法和元素分析儀測(cè)定。收獲時(shí),每個(gè)小區(qū)中選取1 m2的小麥植株,并從地表收割以測(cè)量地上部分干生物量。采集的植物樣品首先在105 ℃下烘1 h殺青,然后再以75 ℃的溫度將它們烘至恒質(zhì)量。記錄此時(shí)的質(zhì)量為地上部分干生物量。然后,選擇10株烘干的植物樣品,用超微粉碎機(jī)(XL-30C,中國(guó))粉碎至0.5 mm。稱取50.0 mg粉碎的植物樣品,用錫紙包裹,然后放入元素分析儀(Vario Macro CN,德國(guó))的自動(dòng)進(jìn)樣盤中以測(cè)量小麥植株中的氮含量。

    1.4 數(shù)據(jù)計(jì)算及統(tǒng)計(jì)分析

    本研究中氮收支根據(jù)Howarth等[20]推薦的模型建立。農(nóng)田氮輸入考慮了4種途徑,包括來(lái)自肥料的氮素、來(lái)自灌溉水的氮素、大氣氮沉降(干+濕)和非豆科作物固定。氮輸出途徑包括作物氮吸收、NH3揮發(fā)、N2O排放和NO3--N淋失。氮輸入和輸出之間的差值為0~100 cm作物主根區(qū)土壤總氮儲(chǔ)量的變化(反映土壤氮素的盈余或虧損)[5]。

    來(lái)自灌溉水的氮素計(jì)算如下:

    式中irrigation為灌溉水中的總氮含量,kg/hm2;depth為灌水量,mm;NI為灌溉水中總氮的濃度,mg/L。本研究中,NI參考馮兆忠等[21]推薦的值,約為22 mg/L。

    大氣氮沉降包括干沉降和濕沉降。干沉降是從大氣輸送到水面和陸地表面的氣態(tài)氮和固態(tài)顆粒狀氮[22]。濕沉降為降水中的氮[23]。在中國(guó)北方地區(qū),干沉降量每年約為50 kg/hm2[24],濕沉降范圍為每年19~29 kg/hm2[4]。本研究將年大氣氮沉降值(干+濕)設(shè)定為每年75 kg/hm2。根據(jù)小麥生育期持續(xù)時(shí)間估算,2019和2020年春小麥生育期內(nèi)大氣氮沉降值約為26 kg/hm2,2021年約為22 kg/hm2。在本研究區(qū),由于春小麥生育期內(nèi)降水量及降水中氮含量均相對(duì)較低,因此忽略了年際間降水變化而造成的濕沉降差異。

    對(duì)于非豆科作物,非生物固氮的范圍為4.5~20 kg/hm2[25]。根據(jù)Liu等[26]的研究,本研究將非豆科作物氮固定值設(shè)為20 kg/hm2。

    小麥生育期內(nèi)氮素吸收量計(jì)算如下:

    式中uptake為小麥吸氮量,kg/hm2;b為收獲時(shí)地上部分干生物量,kg/hm2;TNr為元素分析儀測(cè)定的植物氮含量。

    小麥生育期內(nèi)的累積NH3揮發(fā)量計(jì)算如下:

    式中NH3為累積NH3揮發(fā)量,kg/hm2;avg為小麥生育期內(nèi)的平均NH3揮發(fā)通量,mg/(m2?h);R為第次采樣時(shí)的NH3揮發(fā)通量,mg/(m2?h);為2次相鄰采樣的間隔天數(shù),d;為小麥生育期的總天數(shù),d;M為第次采集的樣品海綿中NH3-N的含量(mg),等于提取液中NH4+-N的含量;為海綿的橫截面積,m2。

    小麥生育期內(nèi)的累積N2O排放量計(jì)算如下:

    式中N2O為累積N2O排放量,kg/hm2;avg為小麥生育期內(nèi)的平均N2O排放通量,mg/(m2?h);F為第次采樣時(shí)的氣體排放通量,mg/(m2?h);為小麥生育期的總天數(shù),d;d為2次相鄰采樣的間隔天數(shù),d;dc/d為第次采樣時(shí)靜態(tài)箱內(nèi)N2O濃度變化速率,10?6/min;為被測(cè)氣體的摩爾質(zhì)量,g/mol;Pa為第次采樣周期內(nèi)的平均大氣壓力,Pa;為第次采樣時(shí)靜態(tài)箱內(nèi)的平均空氣溫度,K;0為標(biāo)準(zhǔn)條件下氣體的體積,m3;0和0為標(biāo)準(zhǔn)條件下的壓力(Pa)和絕對(duì)溫度(K);為靜態(tài)箱高度,cm。

    小麥生育期內(nèi)的NO3--N淋失量計(jì)算如下:

    式中l(wèi)eaching為生育期內(nèi)累積的NO3--N淋失量,kg/hm2;Cl為每次采樣時(shí)浸出0~100 cm土層的滲濾液中NO3--N的濃度,mg/L;V為每次采樣時(shí)測(cè)量的滲濾液體積,L;1為滲漏桶的橫截面積,m2。

    土壤氮素的盈余(或虧損)計(jì)算如下:

    式中sd為土壤中盈余(或虧缺)的氮素(正值表示氮盈余,負(fù)值表示氮虧缺),kg/hm2;input為土壤總氮輸入量,kg/hm2;output為土壤總氮輸出量,kg/hm2。

    每項(xiàng)氮收支占總氮收支的比例計(jì)算如下:

    式中rate為每項(xiàng)氮輸入(或輸出)占總氮輸入(或輸出)的比例,%;N為每項(xiàng)氮輸入(或輸出)的值,kg/hm2;Total是總氮輸出(或輸出)的值,kg/hm2。

    本研究使用SPSS 20.0軟件(SPSS Inc., Chicago, USA)對(duì)數(shù)據(jù)結(jié)果進(jìn)行方差分析(Analysis of Variance, ANOVA)。用最小顯著性差異(Least Significant Difference, LSD)檢驗(yàn)各試驗(yàn)處理間的差異,顯著性水平0.05。

    2 結(jié)果與分析

    2.1 農(nóng)田土壤氮輸入

    2019-2021年春小麥生育期內(nèi)土壤氮收支見(jiàn)圖3??梢钥闯?,土壤總氮輸入量在252~485 kg/hm2范圍內(nèi)。肥料氮是小麥農(nóng)田土壤氮輸入的主要來(lái)源,其次是來(lái)自灌溉水中的氮、大氣氮沉降和非豆科作物固氮(表3)。肥料氮、來(lái)自灌溉水中的氮、大氣氮沉降和非豆科作物固氮分別約占總氮輸入的60%~81%、9%~32%、5%~10%和4%~8%。該結(jié)果表明,在農(nóng)業(yè)生產(chǎn)中進(jìn)行施氮時(shí),除了考慮肥料氮外,由于灌溉水而引入的氮素也不容忽視。因此,農(nóng)田系統(tǒng)的總氮輸入量主要受施氮和灌溉的影響,土壤總氮輸入量隨著施氮量和灌水量的減少而顯著降低。與對(duì)照處理(I1N1)相比,I2N2和I3N3處理使土壤總氮輸入量分別降低了25%和48%。

    注:不同小寫字母表示在0.05的水平上差異顯著。下同。

    表3 2019—2021年春小麥生育期內(nèi)各項(xiàng)氮輸入占總氮輸入的比例

    2.2 農(nóng)田土壤氮輸出

    2019-2021年不同處理土壤氮輸出量見(jiàn)圖3。2019-2021年小麥生育期內(nèi)土壤總氮輸出量分別在338~431、335~457和300~436 kg/hm2范圍內(nèi)。灌溉和施肥均對(duì)總氮輸出量有較明顯的影響。減少灌水量和施氮量均可以顯著降低農(nóng)田土壤的總氮輸出量。與對(duì)照處理相比,I2N2和I3N3處理使總氮輸出量分別降低了9%~11%和20%~23%。在各項(xiàng)氮輸入途徑中,小麥氮素吸收量高達(dá)177~315 kg/hm2。與I1和I2灌溉處理相比,I3灌溉處理使小麥氮吸收量分別降低了10%~17%和5%~14%(圖4),但其值在I1和I2灌溉處理間沒(méi)有顯著差異。此外,小麥吸氮量隨著施氮量的減少呈降低趨勢(shì),但其值在N1和N2施氮處理間也沒(méi)有顯著差異。與N3施氮處理相比,N1和N2施氮處理使小麥吸氮量提高了約23%。

    對(duì)于土壤氮損失而言,隨著灌水量的減少土壤NO3--N淋失量和N2O排放量顯著降低(圖4)。與I1灌溉處理相比,I2和I3灌溉處理使NO3--N淋失量分別降低了6%~15%和8%~24%,同時(shí)使N2O排放量分別降低了17%~27%和31%~42%。與N2O排放和NO3--N淋失相反,土壤NH3揮發(fā)量隨灌水量減少呈增加趨勢(shì)(圖4d~4f)。但總體而言,土壤總氮損失量隨著灌水量減少呈降低趨勢(shì)。與I1灌溉處理相比,I2和I3灌溉處理使土壤總氮損失分別降低了3%~11%和4%~17%。此外,減少施氮量可以同時(shí)降低土壤NO3--N淋失、NH3揮發(fā)和N2O排放。因此,土壤總氮損失量隨著施氮量的減少顯著降低。與N1施氮處理相比,N2和N3施氮處理使土壤NO3--N淋失、NH3揮發(fā)和N2O排放和總氮損失分別降低了4%~22%和11%~26%、20%~21%和18%~31%、11%~16%和24%~27%、13%~21%和7%~27%??紤]灌溉和施氮對(duì)土壤氮損失的綜合影響,與對(duì)照處理相比,I2N2和I3N3處理可使總氮損失分別降低21%~29%和16%~34%。上述這些結(jié)果表明,將灌水量從450 mm減少至315 mm同時(shí)將施氮量從340 kg/hm2減少至250 kg/hm2,小麥吸氮幾乎不受影響,但顯著降低了土壤總氮損失和總氮輸出量。

    小麥氮吸收占土壤氮輸出的比例最大,其次是NO3--N淋失、NH3揮發(fā)和N2O排放(表4)。在2019-2021年小麥氮吸收、NO3--N淋失、NH3揮發(fā)和N2O排放分別占總氮輸出的59%~78%、18%~34%、3%~7%和0.5%~0.9%。盡管小麥吸氮量占總氮輸出的比例在不同灌溉處理間的變化規(guī)律不明顯,但整體上看來(lái),其值仍然隨灌水量的減少呈降低趨勢(shì)。土壤氮損失占總氮輸出的22%~41%,其中NO3--N淋失和NH3揮發(fā)貢獻(xiàn)了土壤氮損失的95%以上。盡管氮損失占土壤總氮輸出的比例在不同灌溉和施肥處理間的變化無(wú)明顯規(guī)律,但與其他處理相比,I2N2處理的氮損失在土壤總氮輸出中所占比例處在較低水平。

    圖4 2019—2021年小麥生育期內(nèi)的氮吸收和土壤氮損失

    表4 2019—2021年春小麥生育期各項(xiàng)氮輸出占總氮輸出的比例

    2.3 土壤剖面NO3--N含量

    2019-2021年春小麥?zhǔn)斋@后0~100 cm根區(qū)土壤NO3--N含量及剖面分布見(jiàn)圖5??梢钥闯觯斋@后的各層土壤NO3--N含量均在0~35 mg/kg范圍內(nèi)變化。灌溉和施肥對(duì)收獲后土壤NO3--N含量及剖面分布有明顯的影響。春小麥?zhǔn)斋@后,I1灌溉處理的土壤NO3--N均集中在80~100 cm土層,而I2和I3灌溉處理的土壤NO3--N主要集中在20~60 cm土層。并且,I1灌溉處理0~60 cm土壤的NO3--N含量明顯低于I2和I3灌溉處理。這一結(jié)果表明,過(guò)高的灌水量會(huì)加速土壤NO3--N的垂向遷移,進(jìn)而降低作物主要根系區(qū)的土壤NO3--N含量。而增加施氮量則顯著提高了0~100 cm土壤的NO3--N含量。例如,在2020和2021年同等灌溉條件下(I2),N1和N2施氮處理0~100 cm土壤的平均NO3--N含量分別比N3施氮處理提高了39%~45%和18%~43%。相反,該結(jié)果也表明減少施氮量則會(huì)顯著降低土壤NO3--N含量??紤]灌溉和施肥對(duì)收獲后土壤NO3--N的綜合影響,I2N1,I3N1和I2N2處理的土壤NO3--N含量處理較高水平,而對(duì)照處理的土壤NO3--N含量?jī)H與I2N2處理相當(dāng),且明顯低于I3N3處理。這些結(jié)果表明,較高的灌水量和較低的施氮量均不利于土壤速效氮累積和土壤肥力保持,適當(dāng)?shù)慕档凸嗨亢驮黾邮┑渴谴盒←溵r(nóng)田節(jié)水和土壤肥力保持的有效途徑。

    2.4 土壤氮盈虧

    土壤氮盈虧反映了生育期末土壤總氮含量的增加或降低。從圖3中可以看出,土壤氮素的盈余或虧損被施氮量決定。在2019—2021年春小麥生育期內(nèi),所有高施氮(340 kg/hm2)處理均表現(xiàn)為土壤氮盈余,所有低施氮(170 kg/hm2)處理均表現(xiàn)為氮虧缺(圖3)。在2019、2020和2021年,高施氮處理各灌溉水平下的土壤氮盈余量分別為46~54、16~47和46~80 kg/hm2,而施氮量降低50%的低施氮處理各灌溉水平下的土壤虧缺量分別為53~101、65~91和48~92 kg/hm2。而在2020和2021年灌溉315 mm(I2)條件下,施氮250 kg/hm2(N2)時(shí)土壤氮素接近平衡狀態(tài)(僅輕度虧缺28~50 kg/hm2,占總施氮量的11%~20%),并且該處理在幾乎不降低小麥氮吸收的前提下有效減少了土壤氮損失。盡管在該灌溉條件下高施氮量處理(I2N1)的氮盈余量也相對(duì)較低,但這種情況由于氮肥用量過(guò)高,因此伴隨著大量的氮損失。這些結(jié)果表明,當(dāng)?shù)貍鹘y(tǒng)的施氮量340 kg/hm2對(duì)于作物生產(chǎn)是過(guò)度的,并且在過(guò)度灌溉的條件下導(dǎo)致了較高的土壤氮損失。過(guò)度的減少施氮量則會(huì)導(dǎo)致土壤氮虧缺,降低土壤肥力。適當(dāng)?shù)臏p少灌水量和施氮量可以有效降低農(nóng)田氮損失帶來(lái)的環(huán)境污染風(fēng)險(xiǎn),且不影響作物生產(chǎn)和土壤肥力。

    圖5 2019—2021年小麥?zhǔn)斋@后0~100 cm土壤剖面NO3--N含量

    3 討 論

    3.1 灌溉施肥對(duì)小麥氮吸收的影響

    灌溉可以通過(guò)改變根區(qū)土壤水分條件來(lái)影響作物根系的活性和生長(zhǎng),從而影響其對(duì)土壤養(yǎng)分的吸收[16,27]。本研究發(fā)現(xiàn),小麥吸氮量隨著灌水量的增加呈增加趨勢(shì),但I(xiàn)1和I2灌溉處理間的小麥吸氮量沒(méi)有顯著差異,即過(guò)高的灌水量對(duì)小麥氮素吸收的改善效果不顯著。造成上述結(jié)果的主要原因有2種解釋,一方面可能是過(guò)度灌溉(I1)導(dǎo)致土壤含水量超過(guò)了田間持水量,因此,將灌水量從I1減少到I2并沒(méi)有明顯降低土壤有效水分含量,對(duì)作物根系活性及其氮吸收的影響較小。而較低的灌水量I3顯著了降低土壤有效水分,導(dǎo)致小麥根系遭受一定程度的水分脅迫,從而限制了小麥對(duì)氮素的吸收[16]。另一方面,過(guò)高的灌水量伴隨著大量的土壤速效氮流失,如NO3--N淋失,N2O排放等,因此不會(huì)明顯增加小麥吸氮量[4]。這些結(jié)果與楊曉亞等[10]的研究結(jié)果一致,他們發(fā)現(xiàn)過(guò)度灌溉沒(méi)有顯著促進(jìn)作物生長(zhǎng),因此,適當(dāng)?shù)慕档凸嗨縿t不會(huì)影響作物對(duì)土壤養(yǎng)分的吸收。

    施用尿素氮肥可以通過(guò)改善根區(qū)土壤速效氮含量來(lái)影響作物生長(zhǎng)和氮素吸收[12-13]。在本研究中,與N3施氮處理相比,N2和N1施氮處理顯著增加了小麥吸氮量,但其值在N1和N2施氮處理間的無(wú)顯著差異。造成這一結(jié)果的原因可能是增加尿素氮肥的施用量提高了作物主根區(qū)的土壤速效氮(NO3--N和NH4+-N)含量,從而使更多的速效氮被小麥植株吸收[16]。而過(guò)量的氮素供應(yīng)可能使土壤氮素含量超過(guò)了植物的最大氮素需求,因此,小麥吸氮量不再隨施氮量的增加而增加[8]。這一結(jié)果許多以前的研究一致[12,18],他們發(fā)現(xiàn)在合理的范圍內(nèi)增加施氮量有效促進(jìn)了作物根系氮吸收,但持續(xù)增加施氮量對(duì)植物氮吸收無(wú)顯著影響。

    3.2 灌溉施肥對(duì)土壤氮損失的影響

    適當(dāng)減少灌水量是降低土壤氮素?fù)p失的有效途徑之一。在本研究中,隨著灌水量的減少,土壤NO3--N淋失量和N2O排放顯著降低。對(duì)于NO3--N淋失量的降低,產(chǎn)生這種結(jié)果的原因是較低的灌水量使上層土壤的壓力水頭較低,難以驅(qū)動(dòng)大量的土壤速效氮淋出0~100 cm根區(qū)土壤[18]。因此,本研究發(fā)現(xiàn),低灌水量處理的土壤NO3--N遷移深度明顯低于高灌水量處理。對(duì)于N2O排放而言,相關(guān)研究表明[28],土壤水分通過(guò)調(diào)節(jié)土壤微生物的有效氧進(jìn)而影響硝化過(guò)程和反硝化過(guò)程,最終影響N2O排放。土壤N2O排放的最佳充水孔隙空間(Water Filled Pore Space,WFPS)較高,一般在60%~70%之間[28]。減少灌水量會(huì)降低土壤含水率,進(jìn)而可能導(dǎo)致土壤WFPS低于適宜N2O排放的最佳WFPS,從而顯著降低了N2O排放。上述結(jié)果與之前許多研究的結(jié)果一致[4,8]。然而,減少灌水量會(huì)增加土壤NH3揮發(fā),是因?yàn)镹H3分子具有親水性,較高的灌溉水可以溶解更多的NH3[6]。Xu等[29]的研究也發(fā)現(xiàn),較低的灌水量會(huì)增加土壤NH3揮發(fā)的風(fēng)險(xiǎn)??傮w而言,將灌水量從I1減少至I2和I3,NO3--N淋失量分別降低了5~17和6~31 kg/hm2,N2O排放量分別降低了0.5~0.9和1~1.6 kg/hm2,而NH3揮發(fā)量?jī)H增加了1.3~3.5和2.9~8.1 kg/hm2。因此,土壤總氮損失隨著灌水量的減少而降低,并且由于灌水量減少而降低的土壤氮損失主要?dú)w因于NO3--N淋失和N2O排放的降低。此外,由于降低灌水量減少了生育期內(nèi)的土壤氮損失,本研究發(fā)現(xiàn)I2和I3灌溉處理收獲后的土壤NO3--N含量高于I1灌溉處理。

    許多研究表明,減少氮肥施用量可以通過(guò)降低作物根區(qū)土壤氮素累積量直接減少土壤氮損失[16,18]。本研究發(fā)現(xiàn),與高施氮處理相比,低施氮處理會(huì)顯著降低收獲后土壤NO3--N累積量。此外,將施氮量從N1減少至N2和N3分別使NO3--N淋失量降低了4%~22%和11%~26%。該結(jié)果與Shi等[12]早期研究結(jié)果一致,他們發(fā)現(xiàn),與高施氮處理相比,低施氮處理顯著降低了土壤速效氮累積量,同時(shí)降低了氮素淋失量。Li等[18]的研究也表明,土壤中較低的氮素積累量可以降低NO3--N淋失進(jìn)入地下水的風(fēng)險(xiǎn)。此外,土壤N2O排放和NH3揮發(fā)均隨施氮量的減少而顯著降低。這一發(fā)現(xiàn)主要是因?yàn)闇p少施氮量直接降低了土壤中硝化過(guò)程和反硝化過(guò)程發(fā)生的基質(zhì)以及NH3揮發(fā)的底物,即NH4+-N和NO3--N。Zhang等[4]也獲得了類似的結(jié)果,他們的研究表明土壤氣態(tài)氮損失量與施氮量顯著正相關(guān)關(guān)系。

    3.3 灌溉施肥對(duì)土壤氮素盈虧的影響

    在農(nóng)田系統(tǒng)中,氮輸入和輸出的不平衡可能導(dǎo)致土壤氮素盈余或虧缺[4,12]。較高的土壤氮素盈余會(huì)增加氮損失的風(fēng)險(xiǎn),而嚴(yán)重的土壤氮素虧缺則會(huì)降低土壤肥力,影響作物生產(chǎn)[13]。在本研究中,N1施氮處理的土壤氮素均處于盈余狀態(tài),可能是因?yàn)?40 kg/hm2的施氮量超過(guò)了作物對(duì)氮素的最大需求,進(jìn)而導(dǎo)致大量氮素殘留在土壤中。本研究也發(fā)現(xiàn),小麥的最高氮素吸收量約為315 kg/hm2,遠(yuǎn)低于當(dāng)?shù)貍鹘y(tǒng)的施氮量。此外,N3施氮處理的土壤氮素均處于虧缺狀態(tài)。顯然是因?yàn)?70 kg/hm2的施氮量不能滿足作物的氮素需求量,并且其生長(zhǎng)過(guò)程中還不可避免地產(chǎn)生了一定程度的氮損失。與N1施氮處理相比,N3施氮處理減少了50%的施氮量,但總氮輸出量?jī)H降低了5%~20%,這表明土壤中累積的氮素被消耗,即N3施氮處理的氮素供應(yīng)量不足。這一結(jié)果與Abera等[30]研究一致,他們發(fā)現(xiàn)當(dāng)肥料氮的施用量減少50%時(shí),總氮損失僅降低了7%~17%。

    3.4 土壤氮平衡分析

    研究表明,估算土壤氮輸入和輸出之間的平衡是評(píng)估農(nóng)田氮肥管理合理性的有效方法[4,20]。但實(shí)際上,在較高的灌水量條件下,高施氮量通常伴隨著高氮輸出,如植株氮吸收、NO3--N的淋失和氣態(tài)氮損失,這種情況也可能接近“土壤氮平衡”。顯然,這種伴隨著高氮輸出的氮平衡不利于農(nóng)業(yè)生產(chǎn)的可持續(xù)發(fā)展,因?yàn)樗鼘?dǎo)致了較高的土壤氮污染物排放,增加了環(huán)境污染和氣候變暖的風(fēng)險(xiǎn)。例如,在中國(guó)西北地區(qū)的春小麥農(nóng)田中,當(dāng)?shù)剞r(nóng)民的灌溉施肥措施,即I1N1處理的土壤氮盈余量相對(duì)較低(圖3),但該措施在所有處理中導(dǎo)致了最高的土壤損失(圖4)。在該處理下,氮損失占總氮輸出的比例高達(dá)25%~35%(表4)。實(shí)際上,在正常農(nóng)業(yè)生產(chǎn)過(guò)程中完全達(dá)到氮輸入和輸出間的平衡幾乎是不可能的,因此,只能在有效降低土壤氮損失,且不影響作物生長(zhǎng)的情況下將土壤氮素的盈余量或虧缺量控制在相對(duì)較低的水平。本研究發(fā)現(xiàn),將灌水量從I1減少至I2并將施氮量從N1減少至N2可以在不影響小麥氮吸收的前提下降低21%~29%的氮損失,并且此時(shí)的土壤氮素處于輕度虧缺狀態(tài),僅虧缺約28~50 kg/hm2,占總施氮量的11%~20%。

    從長(zhǎng)期角度,收獲后的小麥秸稈可作為一種綠色的氮肥來(lái)源用于修復(fù)土壤氮虧缺[31]。以往研究表明,秸稈對(duì)土壤肥力的影響在長(zhǎng)期還田后逐漸體現(xiàn),但不能直接地增加土壤速效氮含量[32]。Christensen等[33]的研究也表明收獲后的小麥秸稈長(zhǎng)期還田可以增加土壤總氮含量。因此,將收獲后的秸稈部分還田可逐漸緩解土壤氮虧缺,但不會(huì)像施入尿素一樣迅速增加土壤NH4+-N和NO3--N含量,這在一定程度上避免了氮損失發(fā)生,尤其是NO3--N淋失。因此,在本研究中,為了緩解I2N2處理的土壤氮虧缺,根據(jù)該處理土壤氮虧缺量及小麥秸稈的氮含量進(jìn)行估算,如果在每季作物收獲后將4~8 t/hm2小麥秸稈還田則保持土壤氮收支平衡。該推薦值與韓新忠等[34]通過(guò)田間試驗(yàn)得出得結(jié)果一致,他們的研究發(fā)現(xiàn)小麥秸稈還田3~7.5 t/hm2對(duì)土壤肥力補(bǔ)充的綜合效果最佳??偠灾?,畦灌315 mm、施氮250 kg/hm2并在收獲后以4~8 t/hm2的小麥秸稈還田的管理模式能夠在維持作物正常生產(chǎn)的前提下有效緩解氮損失帶來(lái)的環(huán)境污染風(fēng)險(xiǎn)和土壤氮收支不平衡問(wèn)題。

    4 結(jié) 論

    本研究于2019-2021年開(kāi)展了為期3 a的田間試驗(yàn),探討了不同灌溉施肥措施對(duì)黃河中上游河套灌區(qū)春小麥農(nóng)田土壤氮損失、作物氮吸收及氮收支的影響,得出了以下結(jié)論:

    1)肥料氮、來(lái)自灌溉水中的氮、大氣氮沉降和非豆科作物固氮分別占土壤總氮輸入的60%~81%、9%~32%、5%~10%和4%~8%。減少灌水量和施氮量可降低農(nóng)田土壤的總氮輸入量。

    2)小麥氮吸收、NO3--N淋失、NH3揮發(fā)和N2O排放分別占總氮輸出的59%~78%、18%~34%、3%~7%和0.5%~0.9%。減少灌水量和施氮量可顯著降低土壤總氮損失,但過(guò)低灌水量和施氮量會(huì)導(dǎo)致小麥吸氮量降低。

    3)當(dāng)?shù)貍鹘y(tǒng)的氮肥施用量N1可導(dǎo)致土壤氮素盈余,而施氮量降低50%的N3處理則導(dǎo)致土壤氮素大量虧缺。對(duì)照處理(I1N1)的土壤氮素?fù)p失量最高。與對(duì)照處理相比,I2N2處理可減少21%~29%的氮損失,且作物氮吸收幾乎未受到影響。

    4)I2N2處理的土壤氮素僅輕度虧缺28~50 kg/hm2,占總施氮量的11%~20%。根據(jù)該處理的土壤氮虧缺量和小麥秸稈含氮量估算,若在收獲后將4~8 t/hm2的小麥秸稈還田,則可保持小麥農(nóng)田土壤的氮收支平衡。

    5)綜合考慮氮污染物減排,作物生產(chǎn)及土壤氮收支平衡在河套灌區(qū)和條件類似地區(qū)的春小麥農(nóng)田中建議采用畦灌315 mm、施氮250 kg/hm2且收獲后4~8 t/hm2小麥秸稈還田的管理模式,該研究可為干旱半干旱地區(qū)春小麥農(nóng)田綠色可持續(xù)發(fā)展提供科學(xué)依據(jù)。

    [1] Dobbie K E, Smith K A. The effects of temperature, water-filled pore space and land use on N2O emissions from an imperfectly drained gleysol[J]. European Journal of Soil Science, 2001, 52: 667-673.

    [2] 劉青麗,李志宏,陳順輝,等. 稻草還田對(duì)煙田追肥氣態(tài)氮損失及相關(guān)微生物的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(22):246-253.

    Liu Qingli, Li Zhihong, Chen Shunhui, et al. Effects of rice-straw returning on gaseous nitrogen loss and microorganisms in tobacco field after topdressing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(22): 246-253. (in Chinese with English abstract)

    [3] Drecht G V, Bouwman A F, Knoop J M, et al. Global modeling of the fate of nitrogen from point and nonpoint sources in soils, groundwater, and surface water[J]. Global Biogeochemical Cycles, 2003, 17(4): 2601-2620.

    [4] Zhang X, Xiao G, Bol R, et al. Influences of irrigation and fertilization on soil N cycle and losses from wheat-maize cropping system in northern China[J]. Environmental Pollution, 2021, 278: 116852.

    [5] 歐陽(yáng)威,郭波波,張璇,等. 北方典型灌區(qū)不同灌期農(nóng)田系統(tǒng)中氮素遷移特征分析[J]. 中國(guó)環(huán)境科學(xué),2013,33(1):123-131.

    Eouyang Wei, Guo Bobo, Zhang Xuan, et al. Transfer characteristics of soil nitrogen in northern typical irrigation area under different irrigation periods[J]. China Environmental Science, 2013, 33(1): 123-131. (in Chinese with English abstract)

    [6] 李銀坤,武雪萍,梅旭榮,等. 常規(guī)灌溉條件下施氮對(duì)溫室土壤氨揮發(fā)的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2011,27(7):23-30.

    Li Yinkun, Wu Xueping, Mei Xurong, et al. Effects of nitrogen application on ammonia volatilization in greenhouse soil under condition of conventional irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(7): 23-30. (in Chinese with English abstract)

    [7] Lam S K, Chen D, Norton R, et al. Influence of elevated atmospheric carbon dioxide and supplementary irrigation on greenhouse gas emissions from a spring wheat crop in southern Australia[J]. Journal of Agricultural Science, 2012, 151(2): 201-208.

    [8] 雷宏軍,楊宏光,劉歡,等. 水肥氣耦合滴灌番茄地土壤N2O排放特征及影響因素分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(11):95-104.

    Lei Hongjun, Yang Honguang, Liu Huan, et al. Characteristics and influencing factors of N2O emission from greenhouse tomato field soil under water-fertilizer-air coupling drip irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(11): 95-104. (in Chinese with English abstract)

    [9] Djaman K, Irmak S, Martin D L, et al. Plant nutrient uptake and soil nutrient dynamics under full and limited irrigation and rainfed maize production[J]. Agronomy Journal, 2013, 105(2): 527-538.

    [10] 楊曉亞,于振文,許振柱. 灌水量和灌水時(shí)期對(duì)小麥耗水特性和氮素積累分配的影響[J]. 生態(tài)學(xué)報(bào),2009,29(2):846-853.

    Yang Xiaoya, Yu Zhenwen, Xu Zhenzhu. Effects of irrigation regimes on water consumption characteristics and nitrogen accumulation and allocation in wheat[J]. Acta Ecologica Sinica, 2009, 29(2): 846-853. (in Chinese with English abstract)

    [11] 谷少委,高劍民,鄧忠,等. 畦灌與施肥時(shí)機(jī)對(duì)土壤硝態(tài)氮分布和冬小麥產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(9):134-142.

    Gu Shaowei, Gao Jianmin, Deng Zhong, et al. Effects of border irrigation and fertilization timing on soil nitrate nitrogen distribution and winter wheat yield[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(9): 134-142. (in Chinese with English abstract)

    [12] Shi Z, Li D, Jing Q, et al. Effects of nitrogen applications on soil nitrogen balance and nitrogen utilization of winter wheat in a rice-wheat rotation[J]. Field Crops Research, 2012, 127: 241-247.

    [13] Liu X, Ju X, Zhang F, et al. Nitrogen dynamics and budgets in a winter wheat–maize cropping system in the North China Plain[J]. Field Crops Research, 2003, 83: 111-124.

    [14] 董文娟,張永平,謝岷,等. 河套灌區(qū)春小麥高產(chǎn)栽培水氮高效利用研究[J]. 干旱區(qū)資源與環(huán)境,2011,25(6):127-131.

    Dong Wenjuan. Zhang Yongping, Xie Min, et al. Efficient utilization of water and nitrogen for spring wheat high-yield in Hetao Irrigation Region[J]. Journal of Arid Land Resources and Environment, 2011, 25(6): 127-131. (in Chinese with English abstract)

    [15] Pan B, Lam S K, Mosier A, et al. Ammonia volatilization from synthetic fertilizers and its mitigation strategies: A global synthesis[J]. Agriculture, Ecosystems & Environment, 2016, 232: 283-289.

    [16] Li Y, Huang G, Chen Z, et al. Effects of irrigation and fertilization on grain yield, water and nitrogen dynamics and their use efficiency of spring wheat farmland in an arid agricultural watershed of Northwest China[J]. Agricultural Water Management, 2022, 260: 107277.

    [17] 茍嘉超,李越,任東陽(yáng),等. 基于基因表達(dá)式編程的作物水分生產(chǎn)函數(shù)構(gòu)建[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(7):100-109.

    Gou Jiachao, Li Yue, Ren Dongyang, et al. A gene expression programming-based crop water production function[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(7): 100-109. (in Chinese with English abstract)

    [18] Li C, Xiong Y, Cui Z, et al. Effect of irrigation and fertilization regimes on grain yield, water and nitrogen productivity of mulching cultivated maize (L.) in the Hetao Irrigation District of China[J]. Agricultural Water Management, 2020, 232: 106065.

    [19] Li Y, Xu J, Liu S, et al. Salinity-induced concomitant increases in soil ammonia volatilization and nitrous oxide emission[J]. Geoderma, 2019, 361: 114053.

    [20] Howarth R W, Billen G, Swaney D, et al. Regional nitrogen budgets and riverine N & P fluxes for the drainages to the North Atlantic Ocean: Natural and human influences[J]. Biogeochemistry, 1996, 35(1): 75-139.

    [21] 馮兆忠,王效科,馮宗煒. 河套灌區(qū)地下水氮污染狀況[J]. 農(nóng)村生態(tài)環(huán)境,2005,21(4):74-76.

    Feng Zhaozhong, Wang Xiaoke, Feng Zongwei. N pollution of groundwater in Hettao Irrigation District[J]. Rural Eco-Environmental, 2005, 21(4): 74-76. (in Chinese with English abstract)

    [22] Anderson K A, Downing J A. Dry and wet atmospheric deposition of nitrogen, phosphorus and silicon in an agricultural region[J]. Water Air &Soil Pollution, 2006, 176: 351-374.

    [23] Xu W, Luo X S, Pan Y P, et al. Quantifying atmospheric nitrogen deposition through a nationwide monitoring network across China[J]. Atmospheric Chemistry and Physics, 2015, 15(21).

    [24] Liu X J, Xu W, Du E Z, et al. Environmental impacts of nitrogen emissions in China and the role of policies in emission reduction[J]. Philosophical Transactions Series A, Mathematical, Physical, and Engineering Sciences, 2020, 378: 20190324.

    [25] Bouwman A F, Drecht G V, Hoek K W V D. Surface N balances and reactive N loss to the environment from global intensive agricultural production systems for the period 1970-2030[J]. Science in China Series C-Life Science, 2005, 48: 767-779.

    [26] Liu C, Watanabe M, Wang Q. Changes in nitrogen budgets and nitrogen use efficiency in the agroecosystems of the Changjiang River basin between 1980 and 2000[J]. Nutrient Cycling Agroecosystems, 2007, 80(1): 19-37.

    [27] 劉健,姚寧,吝海霞,等. 冬小麥物候期對(duì)土壤水分脅迫的響應(yīng)機(jī)制與模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(21):115-124.

    Liu Jian, Yao Ning, Lin Haixia, et al. Response mechanism and simulation of winter wheat phenology to soil water stress[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(21): 115-124. (in Chinese with English abstract)

    [28] Butterbach-Bahl K, Baggs E M, Dannenmann M, et al. Nitrous oxide emissions from soils: How well do we understand the processes and their controls?[J]. Philological Translation Biology Science, 2013, 368: 1-13.

    [29] Xu J, Peng S, Yang S, et al. Ammonia volatilization losses from a rice paddy with different irrigation and nitrogen managements[J]. Agricultural Water Management, 2012, 104: 184-192.

    [30] Abera D, Kibret K, Beyene S, et al. Nitrate leaching under farmers’ fertilizer and irrigation water use in the Central Rift Valley of Ethiopia[J]. International Journal of Physical Science, 2018, 21: 1-17.

    [31] 閆雷,董天浩,喇樂(lè)鵬,等. 免耕和秸稈還田對(duì)東北黑土區(qū)土壤團(tuán)聚體組成及有機(jī)碳含量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(22):181-188.

    Yan Lei, Dong Tianhao, La Yuepeng, et al. Effects of no-tillage and straw returning on soil aggregates composition and organic carbon content in black soil areas of Northeast China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(22): 181-188. (in Chinese with English abstract)

    [32] 田慎重,張玉鳳,邊文范,等. 深松和秸稈還田對(duì)旋耕農(nóng)田土壤有機(jī)碳活性組分的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào),2020,36(2):185-192.

    Tian Shenzhong, Zhang Yufeng, Bian Wenfan, et al. Effects of subsoiling and straw return on soil labile organic carbon fractions in continuous rotary tillage cropland[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(2): 185-192. (in Chinese with English abstract)

    [33] Christensen B T, Olesen, J E. Nitrogen mineralization potential of organomineral size separates from soils with annual straw incorporation[J]. European Journal of Soil Science,1998, 49: 25-36.

    [34] 韓新忠,朱利群,楊敏芳,等. 不同小麥秸稈還田量對(duì)水稻生長(zhǎng)、微生物生物量及酶活性的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2021,31(11):2192-2199.

    Han Xinzhong, Zhu Liqun, Yang Minfang, et al. Effects of different amount of wheat straw returning on rice growth, soil microbial biomass and enzyme activity[J]. Journal of Agro-Environment Science, 2021, 31(11): 2192-2199. (in Chinese with English abstract)

    Irrigation and N application strategies for spring wheat fields in the Hetao irrigation district based on N balance

    Li Yue1,2, Li Gendong3, Chen Zhijun1,4, Zhang Xuechen1,2, Huang Guanhua1,2※

    (1,,100083,;2,,100083,; 3.,015000,; 4.,,110866,)

    Inappropriate irrigation and nitrogen (N) application have caused the severe N losses and budget imbalance in the spring wheat field in the Hetao Irrigation District (41°09′N, 107°39′E), Northwest China. This study aims to determine the rational irrigation and N application strategy for the spring wheat in the study area. A series of field experiments were carried out in the 2019-2021 seasons. A systematic investigation was made to clarify the effects of irrigation and N application on the soil N budget, NO3--N leaching, NH3volatilization, N2O emissions, crop N uptake, as well as soil N surplus and deficit. Three irrigation levels (i.e., 450, 315, and 180 mm (I1, I2 and I3)) and two N fertilization levels (i.e., 340 and 170 kg/hm2(N1 and N3)) were considered in the experiments. Furthermore, the control group was added in the medium N application level (i.e., 250 kg/hm2(N2)) under the condition of a medium irrigation level (i.e., 315 mm (I2)) in 2020 and 2021. The results showed that the fertilizer N was accounted for the largest proportion of soil N input, followed by the irrigation water, atmospheric deposition, and non-bio fixation. The wheat N uptake was accounted for the largest proportion of soil N output, followed by NO3--N leaching, NH3volatilization, and N2O emissions. The total N input in the soil showed a decreasing trend with the decrease of irrigation depth and N application. It infers that the reduced irrigation depth and N application significantly decreased the soil total N losses. But there was the lower threshold in the irrigation depth and N application treatment for the less N uptake of wheat. Among them, the control treatment (i.e., I1N1) resulted in the highest total N input, and the highest soil N losses (NO3--N leaching, NH3volatilization, and N2O emissions). Compared with the control treatment, the I2N2 and I3N3 treatments decreased the total N input by 25% and 48%, respectively, and the soil total N losses by 21%-29% and 16%-34%, respectively. The wheat N uptake in the I3N3 treatment was significantly lower than that in the control treatment. But there was no significant diffidence in the wheat N uptake between the I2N2 and control treatments. The N application presented a significant effect on the soil N surplus and deficit. The N1 treatment was resulted in the higher soil N surplus, while the N3 treatment with about 50% decrease of N application rate was resulted in the higher soil N deficit. However, a slight deficit was obtained in the N2 treatment under the I2 irrigation conditions, where the soil N deficit was only 28-50 kg/hm2. This N deficit value was also accounted for 11%-20% of the total N application amount. Once the 4-8 t/hm2per season wheat straw was returned to the field after harvest, the N balance could be maintained in the wheat field. Thus, the soil N losses can be effectively mitigated from the spring wheat farmlands in the Hetao Irrigation District, Northwest China. The soil N balance can also be expected to improve the irrigation and N application management, particularly with the appropriate straw returning.

    nitrogen; groundwater; greenhouse gas; water and fertilizer management; emissions; pollutants; environment

    10.11975/j.issn.1002-6819.2022.17.007

    S274

    A

    1002-6819(2022)-17-0061-12

    李越,李根東,陳志君,等. 基于氮收支平衡的河套灌區(qū)春小麥農(nóng)田灌溉和施氮策略[J]. 農(nóng)業(yè)工程學(xué)報(bào),2022,38(17):61-72.doi:10.11975/j.issn.1002-6819.2022.17.007 http://www.tcsae.org

    Li Yue, Li Gendong, Chen Zhijun, et al. Irrigation and N application strategies for spring wheat fields in the Hetao irrigation district based on N balance[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(17): 61-72. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.17.007 http://www.tcsae.org

    2022-05-08

    2022-08-10

    國(guó)家自然科學(xué)基金重點(diǎn)國(guó)際合作項(xiàng)目(52220105007);國(guó)家自然科學(xué)基金黃河水科學(xué)研究聯(lián)合基金重點(diǎn)項(xiàng)目(U2243217);內(nèi)蒙古自治區(qū)科技重大專項(xiàng)(2020,2021)

    李越,博士生,研究方向?yàn)檗r(nóng)業(yè)高效用水管理及其生態(tài)環(huán)境效應(yīng)。Email:yue_li@cau.edu.cn

    黃冠華,博士,教授,研究方向農(nóng)業(yè)水文模擬與農(nóng)業(yè)高效用水。Email:ghuang@cau.edu.cn

    tocl精华| 久久婷婷成人综合色麻豆| 亚洲国产精品sss在线观看| 亚洲av日韩精品久久久久久密| 中亚洲国语对白在线视频| e午夜精品久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利18| av有码第一页| 欧美人与性动交α欧美精品济南到| 久久久久国产一级毛片高清牌| 日本vs欧美在线观看视频| 国产免费av片在线观看野外av| 国产精品99久久99久久久不卡| 侵犯人妻中文字幕一二三四区| 亚洲一区中文字幕在线| 无限看片的www在线观看| 亚洲最大成人中文| 曰老女人黄片| 国产精品98久久久久久宅男小说| 97人妻天天添夜夜摸| 在线观看www视频免费| 色老头精品视频在线观看| 天天一区二区日本电影三级 | 99国产极品粉嫩在线观看| 亚洲伊人色综图| 亚洲 国产 在线| 搡老岳熟女国产| 国产免费av片在线观看野外av| 99在线人妻在线中文字幕| 欧美精品亚洲一区二区| 十分钟在线观看高清视频www| 一进一出抽搐动态| 操美女的视频在线观看| 国产区一区二久久| 亚洲成国产人片在线观看| 熟女少妇亚洲综合色aaa.| 国产高清videossex| xxx96com| 一区二区三区国产精品乱码| 国产成+人综合+亚洲专区| 九色亚洲精品在线播放| 老熟妇仑乱视频hdxx| www.999成人在线观看| 真人一进一出gif抽搐免费| 国产人伦9x9x在线观看| 亚洲av五月六月丁香网| 女人高潮潮喷娇喘18禁视频| 大陆偷拍与自拍| 亚洲人成伊人成综合网2020| 丝袜美足系列| 无遮挡黄片免费观看| 极品教师在线免费播放| 97超级碰碰碰精品色视频在线观看| 丝袜美足系列| 久久国产精品男人的天堂亚洲| 欧美激情高清一区二区三区| www.精华液| 亚洲午夜精品一区,二区,三区| 免费少妇av软件| 黄频高清免费视频| 国产精品 欧美亚洲| 一边摸一边抽搐一进一出视频| 国产又爽黄色视频| 麻豆av在线久日| 欧美国产日韩亚洲一区| 亚洲精品国产区一区二| 黑丝袜美女国产一区| 美女 人体艺术 gogo| 色播在线永久视频| av在线播放免费不卡| 亚洲国产精品成人综合色| 亚洲精品国产一区二区精华液| 又黄又爽又免费观看的视频| 国产精品久久久人人做人人爽| 亚洲五月天丁香| 黄色丝袜av网址大全| 国产亚洲av高清不卡| 欧美老熟妇乱子伦牲交| 国产成人一区二区三区免费视频网站| 亚洲午夜理论影院| 在线视频色国产色| 国产av精品麻豆| videosex国产| 老熟妇乱子伦视频在线观看| 看黄色毛片网站| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩瑟瑟在线播放| 岛国在线观看网站| 午夜福利高清视频| 久久人妻熟女aⅴ| 中文字幕色久视频| 手机成人av网站| 久久亚洲精品不卡| 久久久国产成人免费| 国产1区2区3区精品| 神马国产精品三级电影在线观看 | 91成人精品电影| 欧美成人性av电影在线观看| 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 欧美激情极品国产一区二区三区| 亚洲人成电影免费在线| 99国产精品一区二区三区| 中文字幕人成人乱码亚洲影| 久久天躁狠狠躁夜夜2o2o| 成人18禁高潮啪啪吃奶动态图| 不卡一级毛片| 欧美日韩亚洲综合一区二区三区_| 可以免费在线观看a视频的电影网站| 激情在线观看视频在线高清| 欧美老熟妇乱子伦牲交| e午夜精品久久久久久久| 在线十欧美十亚洲十日本专区| 男女做爰动态图高潮gif福利片 | 一边摸一边做爽爽视频免费| 久久久久久久午夜电影| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲av高清不卡| 国产欧美日韩精品亚洲av| 99久久国产精品久久久| 777久久人妻少妇嫩草av网站| 国产私拍福利视频在线观看| 99久久国产精品久久久| 一二三四在线观看免费中文在| 亚洲第一青青草原| aaaaa片日本免费| 99久久国产精品久久久| 欧美精品啪啪一区二区三区| 88av欧美| 亚洲欧美精品综合久久99| 国产精品1区2区在线观看.| 亚洲第一青青草原| 首页视频小说图片口味搜索| 9热在线视频观看99| 他把我摸到了高潮在线观看| 精品一区二区三区四区五区乱码| 男男h啪啪无遮挡| 天天躁狠狠躁夜夜躁狠狠躁| 欧美乱妇无乱码| 国产精品野战在线观看| 亚洲片人在线观看| 一级作爱视频免费观看| 老司机在亚洲福利影院| 色综合站精品国产| 亚洲第一欧美日韩一区二区三区| 色哟哟哟哟哟哟| 亚洲狠狠婷婷综合久久图片| 久久久久精品国产欧美久久久| 色哟哟哟哟哟哟| 黑丝袜美女国产一区| 国产亚洲欧美精品永久| 国产成+人综合+亚洲专区| 国语自产精品视频在线第100页| 久久久国产成人免费| 免费不卡黄色视频| 黄片播放在线免费| 热re99久久国产66热| 亚洲人成网站在线播放欧美日韩| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲伊人色综图| 在线免费观看的www视频| 免费女性裸体啪啪无遮挡网站| 这个男人来自地球电影免费观看| 97超级碰碰碰精品色视频在线观看| 国产在线精品亚洲第一网站| 国产精品亚洲美女久久久| 99久久久亚洲精品蜜臀av| av有码第一页| 一本久久中文字幕| 欧美国产日韩亚洲一区| 亚洲最大成人中文| 在线视频色国产色| 视频区欧美日本亚洲| 久久草成人影院| 精品不卡国产一区二区三区| 久久精品人人爽人人爽视色| 黄色丝袜av网址大全| ponron亚洲| 91av网站免费观看| 国产99久久九九免费精品| 亚洲国产精品999在线| www.www免费av| 精品久久久久久久人妻蜜臀av | 亚洲av第一区精品v没综合| svipshipincom国产片| 动漫黄色视频在线观看| 日本a在线网址| 国产精品一区二区免费欧美| 亚洲国产看品久久| 在线观看日韩欧美| 一二三四在线观看免费中文在| 亚洲中文av在线| 色av中文字幕| 嫩草影视91久久| 人成视频在线观看免费观看| 黑人操中国人逼视频| 美女午夜性视频免费| 欧美成狂野欧美在线观看| 精品久久久精品久久久| 欧美日本中文国产一区发布| 黄色丝袜av网址大全| 久久人人97超碰香蕉20202| 国产成人精品无人区| 人妻丰满熟妇av一区二区三区| 欧美av亚洲av综合av国产av| 国产亚洲av嫩草精品影院| 国产伦一二天堂av在线观看| 女生性感内裤真人,穿戴方法视频| 少妇的丰满在线观看| 亚洲成人免费电影在线观看| 亚洲aⅴ乱码一区二区在线播放 | 夜夜躁狠狠躁天天躁| 最好的美女福利视频网| 久9热在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕人妻丝袜一区二区| 国产精品 欧美亚洲| 日韩精品青青久久久久久| 啪啪无遮挡十八禁网站| 91在线观看av| 长腿黑丝高跟| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人不卡在线观看播放网| 在线观看免费日韩欧美大片| 亚洲av五月六月丁香网| 午夜亚洲福利在线播放| 欧美最黄视频在线播放免费| 91麻豆精品激情在线观看国产| 精品一区二区三区四区五区乱码| 欧美亚洲日本最大视频资源| 久久精品国产亚洲av高清一级| 午夜成年电影在线免费观看| 亚洲成av人片免费观看| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美| 精品国产超薄肉色丝袜足j| 精品人妻1区二区| 欧美黄色片欧美黄色片| 一个人观看的视频www高清免费观看 | 一区二区三区激情视频| 国产蜜桃级精品一区二区三区| 男女下面插进去视频免费观看| 精品无人区乱码1区二区| 国产亚洲精品综合一区在线观看 | 国产精品一区二区免费欧美| 成人特级黄色片久久久久久久| 精品国产美女av久久久久小说| 免费人成视频x8x8入口观看| 美女扒开内裤让男人捅视频| 亚洲欧美日韩另类电影网站| 婷婷六月久久综合丁香| 午夜福利视频1000在线观看 | 在线永久观看黄色视频| 一个人观看的视频www高清免费观看 | 女人精品久久久久毛片| 纯流量卡能插随身wifi吗| 亚洲av熟女| 夜夜躁狠狠躁天天躁| 99re在线观看精品视频| 国产成人欧美在线观看| 精品电影一区二区在线| 久久性视频一级片| 侵犯人妻中文字幕一二三四区| 国产精品av久久久久免费| 9色porny在线观看| 午夜a级毛片| 一级作爱视频免费观看| 在线观看www视频免费| av视频在线观看入口| 看免费av毛片| 深夜精品福利| 久久久久精品国产欧美久久久| 波多野结衣av一区二区av| 麻豆国产av国片精品| 精品福利观看| 大型av网站在线播放| 精品高清国产在线一区| 国产一区二区激情短视频| 成人18禁在线播放| 久久精品aⅴ一区二区三区四区| 99久久精品国产亚洲精品| 久久婷婷成人综合色麻豆| 久久精品91无色码中文字幕| 国产乱人伦免费视频| 日韩欧美免费精品| 在线观看舔阴道视频| 丰满的人妻完整版| 免费在线观看黄色视频的| 亚洲天堂国产精品一区在线| 亚洲成国产人片在线观看| 激情在线观看视频在线高清| 免费看a级黄色片| or卡值多少钱| 中出人妻视频一区二区| 亚洲国产精品合色在线| 亚洲精品久久成人aⅴ小说| 青草久久国产| 精品久久久久久久人妻蜜臀av | 人人妻人人爽人人添夜夜欢视频| 三级毛片av免费| 在线观看日韩欧美| 欧美日韩福利视频一区二区| 精品卡一卡二卡四卡免费| 婷婷精品国产亚洲av在线| 久久中文字幕人妻熟女| 亚洲精品国产色婷婷电影| 亚洲精品中文字幕在线视频| 午夜福利,免费看| 别揉我奶头~嗯~啊~动态视频| 无限看片的www在线观看| 午夜福利一区二区在线看| 两性夫妻黄色片| 免费观看精品视频网站| 男女午夜视频在线观看| 亚洲五月天丁香| 可以在线观看毛片的网站| 日韩成人在线观看一区二区三区| 久久天躁狠狠躁夜夜2o2o| 日日夜夜操网爽| 美女免费视频网站| 久久久久国内视频| 日日夜夜操网爽| 99re在线观看精品视频| 久久香蕉激情| 欧美激情久久久久久爽电影 | 亚洲成a人片在线一区二区| 亚洲专区国产一区二区| 黑丝袜美女国产一区| 久久人妻av系列| 久久久国产欧美日韩av| 后天国语完整版免费观看| 欧美国产精品va在线观看不卡| 久久香蕉国产精品| 国产精品 欧美亚洲| 最好的美女福利视频网| 真人一进一出gif抽搐免费| 久久久久国产精品人妻aⅴ院| 变态另类丝袜制服| 午夜久久久在线观看| 神马国产精品三级电影在线观看 | 日本欧美视频一区| 午夜精品久久久久久毛片777| 又黄又爽又免费观看的视频| 女性生殖器流出的白浆| 日本三级黄在线观看| 成人av一区二区三区在线看| 久久久久久久久久久久大奶| 看黄色毛片网站| 亚洲激情在线av| 免费不卡黄色视频| 激情视频va一区二区三区| 黄频高清免费视频| 中文字幕久久专区| 免费看a级黄色片| 99在线视频只有这里精品首页| 亚洲国产精品合色在线| 亚洲中文av在线| 亚洲天堂国产精品一区在线| 91成人精品电影| 香蕉国产在线看| 色老头精品视频在线观看| 91成年电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 人人妻人人澡欧美一区二区 | 国产99白浆流出| 国产精品 国内视频| 欧美日韩瑟瑟在线播放| 一区二区三区精品91| 亚洲国产高清在线一区二区三 | 女人高潮潮喷娇喘18禁视频| 黄片大片在线免费观看| 757午夜福利合集在线观看| 国产一区二区在线av高清观看| www.熟女人妻精品国产| 久久欧美精品欧美久久欧美| 欧美成人一区二区免费高清观看 | 12—13女人毛片做爰片一| 满18在线观看网站| 99国产精品99久久久久| 在线观看免费视频网站a站| 精品国产一区二区久久| 欧美精品亚洲一区二区| 午夜福利,免费看| 别揉我奶头~嗯~啊~动态视频| www日本在线高清视频| 亚洲久久久国产精品| 国产亚洲欧美精品永久| 亚洲伊人色综图| 亚洲欧美一区二区三区黑人| 成年女人毛片免费观看观看9| 亚洲 欧美 日韩 在线 免费| 男女之事视频高清在线观看| 怎么达到女性高潮| 免费少妇av软件| 国内精品久久久久久久电影| 黄色成人免费大全| 91成年电影在线观看| 婷婷丁香在线五月| 身体一侧抽搐| 欧美成人性av电影在线观看| 看片在线看免费视频| 大香蕉久久成人网| 亚洲中文日韩欧美视频| 少妇被粗大的猛进出69影院| 欧美人与性动交α欧美精品济南到| 日韩精品中文字幕看吧| 成人免费观看视频高清| 人人澡人人妻人| 午夜影院日韩av| 69精品国产乱码久久久| 9热在线视频观看99| 国内毛片毛片毛片毛片毛片| 看免费av毛片| 国产精品久久久人人做人人爽| 亚洲一区高清亚洲精品| 一区二区三区激情视频| 国产成人精品久久二区二区91| 精品卡一卡二卡四卡免费| 男女床上黄色一级片免费看| 亚洲国产中文字幕在线视频| 国产av在哪里看| 后天国语完整版免费观看| 九色亚洲精品在线播放| 欧美日韩乱码在线| 91精品国产国语对白视频| 欧美黄色片欧美黄色片| 欧美国产日韩亚洲一区| 99re在线观看精品视频| 国内精品久久久久久久电影| 精品人妻1区二区| 狂野欧美激情性xxxx| 国产一区在线观看成人免费| 亚洲在线自拍视频| 国产91精品成人一区二区三区| 欧美一级a爱片免费观看看 | 婷婷六月久久综合丁香| 日韩大尺度精品在线看网址 | 女性被躁到高潮视频| 日本 欧美在线| 动漫黄色视频在线观看| 久久久精品国产亚洲av高清涩受| 亚洲国产精品999在线| 国产亚洲欧美精品永久| 宅男免费午夜| 国产激情欧美一区二区| 久久人妻福利社区极品人妻图片| 欧美精品亚洲一区二区| 99在线视频只有这里精品首页| 神马国产精品三级电影在线观看 | 精品国产一区二区久久| 中文字幕人成人乱码亚洲影| xxx96com| 成人免费观看视频高清| 757午夜福利合集在线观看| 又大又爽又粗| 女人被躁到高潮嗷嗷叫费观| 国产成人欧美在线观看| 1024香蕉在线观看| 一区二区三区国产精品乱码| 亚洲一码二码三码区别大吗| 一区二区三区激情视频| 精品无人区乱码1区二区| 国产xxxxx性猛交| 欧美大码av| 久久香蕉国产精品| 少妇熟女aⅴ在线视频| 免费看十八禁软件| 天天躁夜夜躁狠狠躁躁| 人妻久久中文字幕网| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 男人舔女人下体高潮全视频| 精品日产1卡2卡| 亚洲美女黄片视频| 国内精品久久久久久久电影| 亚洲 国产 在线| 国产精品美女特级片免费视频播放器 | 日日爽夜夜爽网站| 黄片小视频在线播放| 高清在线国产一区| 精品人妻在线不人妻| 欧美黄色淫秽网站| 男女午夜视频在线观看| 中文字幕人妻熟女乱码| 亚洲电影在线观看av| 久久精品91无色码中文字幕| 亚洲国产精品合色在线| 国产亚洲精品久久久久5区| 亚洲av片天天在线观看| 久久久久久国产a免费观看| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| www.精华液| 欧美午夜高清在线| 中文字幕高清在线视频| 99国产精品一区二区蜜桃av| 国产成人欧美| 中文字幕色久视频| 国产成年人精品一区二区| 91九色精品人成在线观看| 久久精品91无色码中文字幕| 嫩草影视91久久| 精品久久久久久成人av| 天堂√8在线中文| 亚洲自偷自拍图片 自拍| 国产成人影院久久av| 国产精品久久久久久亚洲av鲁大| svipshipincom国产片| 美女高潮到喷水免费观看| 不卡av一区二区三区| 青草久久国产| 精品久久久久久,| 国产日韩一区二区三区精品不卡| 极品教师在线免费播放| 亚洲欧美激情在线| e午夜精品久久久久久久| 精品国产美女av久久久久小说| 精品高清国产在线一区| 色精品久久人妻99蜜桃| 日本免费一区二区三区高清不卡 | 岛国视频午夜一区免费看| 欧美乱码精品一区二区三区| 男人的好看免费观看在线视频 | 每晚都被弄得嗷嗷叫到高潮| 欧美一级毛片孕妇| 一二三四社区在线视频社区8| 欧美乱色亚洲激情| 成人手机av| 久久人人精品亚洲av| 老司机午夜十八禁免费视频| 久久精品aⅴ一区二区三区四区| 可以在线观看毛片的网站| av视频在线观看入口| 悠悠久久av| 人妻久久中文字幕网| 99在线人妻在线中文字幕| 亚洲国产精品合色在线| 在线观看日韩欧美| 欧美黑人精品巨大| 久久久久久久久中文| 视频在线观看一区二区三区| 日本五十路高清| 国产片内射在线| 嫩草影院精品99| 精品少妇一区二区三区视频日本电影| 非洲黑人性xxxx精品又粗又长| 日韩欧美一区二区三区在线观看| 看黄色毛片网站| 在线观看午夜福利视频| 国产亚洲精品第一综合不卡| 人人澡人人妻人| 日本a在线网址| 成人av一区二区三区在线看| 久久久国产成人精品二区| 深夜精品福利| 女性被躁到高潮视频| 在线观看日韩欧美| 久久精品国产亚洲av香蕉五月| 波多野结衣巨乳人妻| 视频在线观看一区二区三区| 男女午夜视频在线观看| 黄色a级毛片大全视频| 久久精品国产亚洲av高清一级| 女人被躁到高潮嗷嗷叫费观| 欧美中文综合在线视频| 一个人免费在线观看的高清视频| 亚洲人成网站在线播放欧美日韩| 性欧美人与动物交配| 中文字幕人妻丝袜一区二区| 亚洲一区中文字幕在线| 三级毛片av免费| 1024香蕉在线观看| 免费在线观看视频国产中文字幕亚洲| 精品欧美国产一区二区三| 中文字幕人成人乱码亚洲影| 一级毛片高清免费大全| 色精品久久人妻99蜜桃| 91成人精品电影| 日韩大尺度精品在线看网址 | 少妇粗大呻吟视频| 日本在线视频免费播放| av在线天堂中文字幕| 久久精品91蜜桃| 午夜福利成人在线免费观看| 90打野战视频偷拍视频| 久久精品91蜜桃| 亚洲av电影不卡..在线观看| 淫秽高清视频在线观看| 精品国产亚洲在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲无线在线观看| 精品国产亚洲在线| 美女扒开内裤让男人捅视频| 天天添夜夜摸| 麻豆久久精品国产亚洲av| 国产精品久久电影中文字幕| 黄色视频,在线免费观看| 999精品在线视频| 亚洲国产精品久久男人天堂| 久久中文字幕一级| 精品国产亚洲在线| 国产精品久久电影中文字幕| 99国产综合亚洲精品| 亚洲一区二区三区色噜噜| √禁漫天堂资源中文www| 免费高清在线观看日韩| 亚洲精品中文字幕在线视频| √禁漫天堂资源中文www| 国内毛片毛片毛片毛片毛片| 纯流量卡能插随身wifi吗| 国产aⅴ精品一区二区三区波| 狠狠狠狠99中文字幕| 乱人伦中国视频| 怎么达到女性高潮|