紅壤區(qū)旱作花生地剖面氮素遷移轉(zhuǎn)化模擬研究

摘要：為探究紅壤區(qū)花生種植地氮素的垂向遷移影響因素及其損失特征，本研究利用HYDRUS-ID模型構(gòu)建適用于紅壤地區(qū)的水氮耦合模型。該模型結(jié)合實地監(jiān)測數(shù)據(jù)和室內(nèi)淋溶實驗結(jié)果，對模型參數(shù)進行率定和驗證，并利用校驗過的模型模擬研究區(qū)不同環(huán)境條件下土壤剖面氮素的動態(tài)遷移與轉(zhuǎn)化過程。結(jié)果表明：HYDRUS-ID模型模擬值與實測值之間吻合度較高，模型可用于模擬花生地不同深度土壤水氮運移過程。土壤銨態(tài)氮和硝態(tài)氮在降雨施肥后迅速升高，隨后下降。施肥量變化主要影響0-40cm土層內(nèi)氮素濃度變化，而改變降雨量影響0-80 cm土層內(nèi)氮素的遷移轉(zhuǎn)化。此外，氮素濃度隨土層深度的增加而減小。模擬結(jié)果顯示，施氮量為225．0 kg·hm-2時，雨季發(fā)生總降雨量大于173 mm且降雨強度大于11 mm·d-1降雨事件時，或旱季發(fā)生強降雨事件，即總降雨量大于135 mm且降雨強度大于135 mm·d-1時，80 cm土層硝態(tài)氮濃度至少上升到1mg·L-1。受紅壤質(zhì)地黏重的影響，深層土壤氮素對降雨響應具有一定的滯后效應。隨著施肥量的增加，作物根系對氮素的吸收通量增大，但氮素吸收量與氮肥施用量的比例卻有所下降；降雨量對根系吸收氮素的影響相對較小。在花生生長季，土壤中的主要氮損失形式為硝態(tài)氮的淋失。模擬分析表明，降雨和施肥均會影響硝態(tài)氮的淋失量，其中降雨量的增加顯著提高了硝態(tài)氮的浸出量，硝態(tài)氮累積淋失量變化范圍在5.23-42.97 kg·hm-2之間，證明了花生地對紅壤區(qū)地下水硝酸鹽潛在污染不可忽視。

關(guān)鍵詞：紅壤；銨態(tài)氮；硝態(tài)氮；氮素運移；HYDRUS-ID

中圖分類號：S565.2 文獻標志碼：A 文章編號：1672-2043（2024）07-1557-11 doi：10.11654/jaes.2023-1032

南方紅壤區(qū)是我國重要的土地資源，占全國耕地的27.8%，在農(nóng)業(yè)生產(chǎn)中占有重要地位。為提高作物產(chǎn)量，農(nóng)民往往過量施用水肥。不合理的氮肥施用不僅加劇了土壤酸化問題，還影響了作物品質(zhì)。此外，該地區(qū)降雨季節(jié)分布不均，極端降雨事件頻發(fā)，導致土壤的保水保肥能力下降，從而增加了氮素的淋溶風險。研究發(fā)現(xiàn)，在雨季大量降水的人滲作用下，紅壤農(nóng)田淺層土壤中的氮素（尤其是硝態(tài)氮）容易滲入地下水，造成地下水中氮素含量過高。因此，深入研究紅壤區(qū)氮素在土壤中的垂向遷移轉(zhuǎn)化過程，定量評估氮素淋失量，對提高氮肥利用率和控制地下水污染具有重要意義。

當前研究農(nóng)田氮素遷移轉(zhuǎn)化的方法主要包括野外原位監(jiān)測和模型模擬。盡管田間試驗數(shù)據(jù)對于探究氮肥的轉(zhuǎn)化過程和損失途徑至關(guān)重要，但其全面性存在一定的限制。因此，利用數(shù)值模擬技術(shù)來研究土壤中水氮的動態(tài)變化已成為評估土壤生態(tài)健康的關(guān)鍵手段。在此背景下，許多研究人員對不同水肥條件下的土壤水分特性和氮素分布進行了研究，取得了豐富的成果。Yi等通過田間監(jiān)測和模擬，對黑河流域典型玉米田的水分平衡進行量化，并利用模型提出了優(yōu)化的灌溉方案。Mo等模擬了暴雨、大暴雨和極端暴雨三種降雨情景下氮素滲入地下水的變化情況，發(fā)現(xiàn)隨著降雨強度的增加，尿素態(tài)氮的浸出量大幅上升。Azad等探討氣候變化影響下玉米不同灌溉施肥管理下硝酸鹽向地下水的淋溶規(guī)律，并預測硝酸鹽向地下水轉(zhuǎn)移的量和速率。Krevh等研究表明土壤硝酸鹽的淋溶和積累主要受降水條件、土壤性質(zhì)和蒸散作用的影響，因此針對每種特定的環(huán)境條件優(yōu)化氮肥施用以減少硝酸鹽損失十分重要。

花生作為我國南方紅壤區(qū)的主要油料經(jīng)濟作物，種植面積和產(chǎn)量均位居全國前列。然而，由于耕作管理措施不當，花生種植期恰逢雨季，在土壤淋溶作用下氮素易隨降雨深層滲漏進入水體導致地下水污染。目前關(guān)于農(nóng)田氮素運移的相關(guān)研究主要集中于麥田、玉米地和稻田等土地利用方式上，對于旱作花生地的氮素運移和流失問題尚缺乏系統(tǒng)性研究。鑒于此，本研究選取南方紅壤區(qū)的旱作花生地為研究對象，通過室內(nèi)土柱試驗和HYDRUS-ID模擬，探討不同環(huán)境條件下旱作花生地土壤氮素遷移轉(zhuǎn)化的動態(tài)響應關(guān)系，闡明花生地土壤剖面氮輸移的主要過程及其影響因素，旨在為制定區(qū)域內(nèi)可持續(xù)的農(nóng)業(yè)氮肥利用管理策略提供參考價值，并為南方紅壤區(qū)有效控制非點源污染提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

研究區(qū)在江西省鷹潭市孫家農(nóng)業(yè)小流域，位于江西省鷹潭市余江縣劉家站（28°12＇N，116°55＇E）。該地區(qū)屬于中亞熱帶季風氣候，年平均氣溫17.8℃，年平均潛在蒸發(fā)量1378 mm。年平均降雨量1795 mm，且降水具有明顯的季節(jié)性，4-6月份降雨量約占全年降水量的50%。供試土壤為第四紀紅色黏土發(fā)育而來的典型紅壤。流域內(nèi)主要土地利用方式包括花生地、稻田和林地等。本研究所選研究對象為雨養(yǎng)農(nóng)用地，種植作物以花生為主，種植時間為4月中旬至8月中旬?；ㄉ诿磕?月進行一次性施肥，包括尿素、復合肥及少量K肥和P肥。

1.2 試驗設計

選取區(qū)域內(nèi)旱作紅壤作為土壤樣品，按照0-20、20-40、40-80 cm分層，挖出1 m×1 m×0.8 m的正方形深坑，每層取環(huán)刀原狀土樣以測量容重，并從每個土層平行采集10個樣品進行混合后于自封袋封存樣品，用于測定土壤基本理化性質(zhì)和制備土柱。

采集的土壤樣品在經(jīng)過自然風干處理后，清除植物殘體、根須等雜質(zhì)，并在混勻和研磨后過2 mm篩，用于測定土壤背景值和分析土壤基本理化性質(zhì)。本研究采用環(huán)刀法測定土壤容重；烘干法測定土壤含水率；甲種比重計法測量土壤機械組成；重鉻酸鉀容量法—外加熱法測定土壤有機質(zhì)；土壤中的銨態(tài)氮和硝態(tài)氮經(jīng)氯化鉀浸提后，分別采用納式試劑—分光光度法和紫外分光光度法進行測定。測定結(jié)果詳見表1。

模擬淋溶裝置為內(nèi)徑10 cm圓柱形亞克力柱，柱高90 cm。填裝土壤前，先鋪上200目濾布和2cm潔凈石英砂粒，然后在石英砂層上再鋪設200目尼龍濾布分隔石英砂和土壤界面。根據(jù)野外采集的環(huán)刀原狀土樣所測定的容重，依次將土壤分層填入亞克力管并壓實，模擬80 cm原狀土柱，土柱表面同樣鋪設濾布和石英砂。在土柱側(cè)面20、40、80 cm處設置3個取樣口，土柱底部設置排水口，用于承接淋溶液。

土柱制備完成后，為避免初始背景值的干擾，使用去離子水從上到下濕潤土柱，并充分淋洗2d。為模擬研究區(qū)常規(guī)水肥條件下的氮素淋溶情況，對當?shù)鼗ㄉ氐氖┓是闆r進行了調(diào)查。調(diào)查結(jié)果顯示，該地區(qū)花生尿素用量通常在150-300 kg·hm-2。基于年平均降雨量與實驗周期兩因素綜合考慮后設定土柱模擬條件：（1）率定條件：灌水量350 mm，尿素量200mg；（2）驗證條件：灌水量437.5 mm，尿素量250 mg。施肥時，將稱量好的尿素與土柱表層的石英砂混合均勻后回填至土柱內(nèi)，并通過馬氏瓶將水加至土壤表面。按一定的時間間隔（1、2、3、5、8、13、21、30 d）取土柱20、40、80 cm土樣和底部淋溶液，分別測定銨態(tài)氮和硝態(tài)氮的濃度。

1.3 HYDRUS-1D模型構(gòu)建與評價

1.3.1 基本方程

（1）非飽和水流動力學方程

HYDRUS-1D軟件水流模型的控制方程為改進后的Richards方程，其表達式如下：

式中：θ為土壤體積含水率，cm3·cm-3；t為時間，d；z為垂向坐標，cm，向下為正；h為壓力水頭，cm；K（h）為非飽和滲透系數(shù)，cm·d-1；s表示根系吸水項，cm3·cm-3·d-1。本次模擬運用van Genuchten模型來計算K（h）的值，表達式如下：

式中：θ為土壤體積含水率，cm3·3cm-3；θr為土壤殘余含水率，cm3·cm-3：θs為土壤飽和含水率，cm3·cm-3；α為進氣值倒數(shù)，cm-1；m為土壤水分特征曲線適線參數(shù)；n為孔徑分布系數(shù)，n＞1，m=1 - 1/n；Ks為飽和導水率，cm·d-1；Se為有效水分含量，無量綱；l表示孔徑連通系數(shù)，一般取經(jīng)驗值0.5。

作物根系吸水量采用Feddes模型來計算，表達式如下：

S（h）=α（h）Sp（5）

Sp=b（x）Tp（6）

式中：S（A）為根系實際吸水量；h為土壤基質(zhì)勢；Sp為潛在根系吸水量，d-1；α（A）為水分脅迫反應函數(shù)；b（x）為標準化根系吸水分布函數(shù)，cm-1；Tp為作物潛在蒸騰率，cm·d-1。

（2）氮素遷移與轉(zhuǎn)化方程

本實驗以施用尿素來討論氮素在土壤剖面的轉(zhuǎn)化和垂向運移過程。尿素的一階衰變鏈如圖2。

氮的轉(zhuǎn)化過程主要包括：尿素的水解作用；有機氮礦化成銨態(tài)氮；銨態(tài)氮的硝化作用；硝態(tài)氮的反硝化；土壤對銨態(tài)氮的吸附；作物對無機氮的吸收。礦化被認為是零級反應，其余反應均為一級反應。

尿素、銨態(tài)氮和硝態(tài)氮的對流—彌散運移方程為：式中：下標1，2和3分別表示尿素、銨態(tài)氮和硝態(tài)氮；θ為體積含水率，cm3·m-3；cw和cs分別為溶質(zhì)在液態(tài)和固態(tài)中的濃度，mg·L-1;ρ為土壤容重；D為溶質(zhì)在土壤層中的擴散系數(shù)，cm2·d-1；q為體積水通量，cm2·d-1；μ為一階氮素轉(zhuǎn)化常數(shù)，d-1；μ'為聯(lián)系尿素、銨態(tài)氮和硝態(tài)氮之間的一階氮素轉(zhuǎn)化常數(shù)，d-1；γ為零階氮素轉(zhuǎn)化常數(shù)，mg·L-1·d-1。銨態(tài)氮為吸附性物質(zhì)，假設其吸附過程符合線性等溫吸附，其中吸附方程為：cs，2=Kd·Cw，2，Kd為線性吸附平衡常數(shù)，cm3·g-1。ra，i=S·cw，i，ra，3、Ta，3分別表示作物根系吸收銨態(tài)氮和硝態(tài)氮的量，mg·L-1·d-1。

1.3.2 邊界條件與初始條件

水分運動的上邊界條件設置為有表面積水的大氣邊界條件，最大積水厚度（hmax）為6cm；下邊界條件則設置為自由排水邊界條件。水分初始條件采用土壤剖面的初始含水量；溶質(zhì)運移上邊界條件設置為濃度通量邊界，下邊界設置為零濃度梯度邊界。實驗主要考慮的溶質(zhì)為銨態(tài)氮和硝態(tài)氮，土柱實驗不考慮作物根系吸水作用。

1.3.3 參數(shù)確定

土壤水力參數(shù)初始值采用土壤基本物理性質(zhì)（表1）和神經(jīng)網(wǎng)絡Rosetta模型進行測算，再結(jié)合土柱淋溶實驗所測數(shù)據(jù)校驗，得出各土層土壤水力參數(shù)率定值，具體結(jié)果詳見表2。

對于溶質(zhì)運移參數(shù)，基于前人研究，銨態(tài)氮和硝態(tài)氮在土壤水中的擴散系數(shù)分別取經(jīng)驗值1.52cm2·d-1和1.64 cm2·d-1，銨態(tài)氮吸附系數(shù)取3.5 cm3·g-1。并利用銨態(tài)氮和硝態(tài)氮濃度對溶質(zhì)傳輸和轉(zhuǎn)化參數(shù)進行校準，結(jié)果見表3。

1.3.4 模型校驗

采用決定系數(shù)（R2）、均方根誤差（RMSE）和平均絕對誤差（MAE）對模擬精度進行評價。

2 結(jié)果與分析

2.1 模型率定與驗證

基于HYDRUS-1D模型對旱作花生地剖面0-80cm的土壤含水量、銨態(tài)氮和硝態(tài)氮濃度進行模擬計算，將其模擬值與花生地不同土層的實測值進行了比較（圖3），并對模型精度進行統(tǒng)計評價，結(jié)果如表4所示。0-80 cm土壤含水率R2、RMSE、MAE的變化范圍分別為0.882-0.985、0.005-0.014、0.004-0.011 cm3·cm-3。整體來看，不同土層土壤含水率模擬值與實測值吻合度較高。各層土壤銨態(tài)氮和硝態(tài)氮濃度模擬值與實測值之間的R2、RMSE、MAE分別為0.872-0.978、0.018-5.198、0.015-3.600 mg·L-1。除了0-20cm土層個別實測值與模擬值存在偏差外，整個土層土壤銨態(tài)氮和硝態(tài)氮濃度的模擬值與實測值變化趨勢基本接近。因此模型擬合效果較好，可應用于模擬紅壤區(qū)花生地土壤氮素遷移轉(zhuǎn)化過程。

2.2 模型應用

利用校驗好的模型參數(shù)對研究區(qū)2022年花生生長季水氮運移進行模擬與分析。上邊界條件設置為有地表徑流的大氣邊界條件，輸入日降水數(shù)據(jù)與蒸散發(fā)數(shù)據(jù)，并加入根系吸收模型，其余參數(shù)不作調(diào)整。利用Penman-Monteith公式和花生作物系數(shù)計算花生潛在蒸散量ETp：

ETp＝Kc·ETo（10）

式中：Ke為花生作物系數(shù)，苗期取0.58，花針期取0.8，結(jié)莢期取1.07，成熟期取0.65；ET0為參考作物潛在蒸散量。

潛在蒸散量ETp可分為潛在土壤蒸發(fā)量（Ep）和Tp兩部分：

Ep=ETp·eK·LAI（11）

Tp=ETp·（1-eK·LAI）（12）

式中：K為消光系數(shù)；LAI為葉面積指數(shù)。

模擬時間為2022年4月15日至8月12日共120d，0-90 d為雨季，91-120 d為旱季（圖4）。為研究花生生長期間施肥和降雨對土壤氮素運移的時空變化規(guī)律及作物氮素吸收情況的影響，設定了兩種試驗模式進行模擬分析：M模式，即氣象站實際降雨量（總降雨770 mm）下的不同施肥量（M1=112.5 kg·hm-2、M2=225.0 kg·hm-2、M3=337.5 kg·hm-2），氮肥一次性施入；N模式，即常規(guī)施肥條件（225.0 kg·hm-2）下的不同降雨量（N1=385 mm、N2=770 mm、N3=1 155 mm）模式，模擬0-80 cm土壤剖面氮素的變化過程及作物氮素吸收情況。

2.2.1 不同環(huán)境下土壤氮素濃度變化特征

花生地施用氮肥后，氮肥在表層土壤中經(jīng)水解產(chǎn)生大量銨態(tài)氮。隨后，這些銨態(tài)氮通過硝化作用轉(zhuǎn)化為硝態(tài)氮，導致硝態(tài)氮濃度迅速增加。在降雨作用下表層土壤的氮素通過水分傳輸遷移到深層土壤。但由于銨態(tài)氮具有較強的吸附性，大部分仍固定在表層土壤中，只有少量銨態(tài)氮能達到深層土壤。與此同時，由于硝態(tài)氮與土壤膠粒均帶負電荷，導致土壤中的硝態(tài)氮易隨水流運移，從而引起氮浸出。

圖5顯示了在M模式下花生生長季土壤剖面中氮素濃度的動態(tài)變化情況。結(jié)果表明，施肥量對0-20 cm土層土壤銨態(tài)氮量和0-40 cm土層土壤硝態(tài)氮量有顯著影響。隨著施肥量的增加，同一土層中氮素濃度的峰值變高，且峰值出現(xiàn)的時間更早。具體來看，在0-20 cm土層，M2和M3處理下銨態(tài)氮濃度峰值相較于M1分別提高了3.53%和24.31%，而硝態(tài)氮濃度峰值分別提高了44.51%和89.27%；在20-40 cm土層內(nèi)，M2和M3處理下硝態(tài)氮濃度峰值相較于M1分別提高了5.10%和23.92%。從土壤剖面來看，氮素隨土層深度的增加而減小，且M3降幅最大。受紅壤質(zhì)地黏重的影響，深層土壤中氮素濃度達到峰值的時間較晚。由此可知，0-40 cm土層氮素濃度對施肥量響應更加明顯，且主要表現(xiàn)在施肥初期。

圖6反映了花生生育期內(nèi)不同土層在N模式下土壤氮素變化情況。結(jié)果表明，氮素濃度隨降雨量的增加呈現(xiàn)出顯著的波動，0-80 cm不同土層氮素濃度均有不同程度增加。同一土層各處理土壤氮素濃度隨時間的波動趨勢相似，波動幅度隨土層深度增加而降低。在20 cm土層中，N2和N3處理下銨態(tài)氮濃度峰值相較于N1分別提高10.92%和92.44%，而硝態(tài)氮濃度峰值分別提高48.72%和100.20%；在40 cm土層內(nèi)，與N1相比，N2和N3處理下銨態(tài)氮濃度峰值平均分別增加32.77%和42.86%，硝態(tài)氮濃度峰值平均分別增加13.79%和75.10%。在80 cm土壤處，銨態(tài)氮濃度幾乎可以忽略，而硝態(tài)氮濃度受降雨影響仍變化較大，其在N1、N2和N3處理下的最高值分別為0.64、4.22 mg·L-1和7.28 mg·L-1?？傮w而言，降雨量對0-40 cm土層內(nèi)銨態(tài)氮濃度有一定影響，而0-80 cm土層內(nèi)硝態(tài)氮濃度受降雨量作用變化顯著。

2.2.2 作物氮素吸收規(guī)律

如圖7所示，在M模式下，隨著氮肥施用量的增加，作物根系氮素吸收通量增大，吸氮量與施氮量的比值卻呈下降趨勢。其中，M1、M2和M3處理下銨態(tài)氮的吸收量占施氮量的比例分別為16.44%、9.71%和7.54%；硝態(tài)氮的吸收量占施氮量的比例分別為64.27%、36.20%和25.77%。相比之下，不同降雨量對作物根系氮素吸收通量的影響較小。各處理銨態(tài)氮的吸收量占施氮量的比例范圍為8.45%-9.84%；硝態(tài)氮的吸收量占施氮量的比例范圍為33.98%-36.20%。

2.2.3 硝態(tài)氮淋失特征

圖8顯示了不同模式下硝態(tài)氮累積淋失通量隨時間變化圖。在M模式下，增加氮肥量后硝態(tài)氮的累積淋失量有升高趨勢但變化較小，相比于M1，M2和M3處理下的硝態(tài)氮淋失量僅分別增加了4.42%和8.93%。而在N模式下，硝態(tài)氮的累積淋失量隨降雨量的增加呈顯著上升趨勢。N2和N3處理下的淋失量相較于N1分別增加了148.18%和721.61%。這一發(fā)現(xiàn)表明，降雨量對硝態(tài)氮淋失的影響顯著大于氮肥施用量，降雨量的增加使得土壤硝態(tài)氮的淋失量呈非線性增長的趨勢。

3 討論

本研究基于土柱實驗數(shù)據(jù)，運用HYDRUS-1D模型進行了旱作花生地水氮運移模擬。結(jié)果表明，模型模擬得到的土壤含水率、銨態(tài)氮和硝態(tài)氮濃度與室內(nèi)試驗的實測值具有較好的一致性，且土壤含水率的模擬誤差小于土壤銨態(tài)氮和硝態(tài)氮的誤差，與汪順生等研究結(jié)果一致。這一現(xiàn)象可能源于環(huán)境條件對土壤氮素含量影響較大，土壤中的微生物活動、氧含量以及溫度等因素都會影響土壤氮素的遷移轉(zhuǎn)化過程。相關(guān)研究表明，深層土壤水氮模擬效果優(yōu)于淺層，這與本研究模擬結(jié)果類似。可能是由于淺層土壤易受灌溉和施肥等因素影響。總體上看，模型各評價指標均在可接受范圍內(nèi)，模型可適用于紅壤區(qū)花生地氮素遷移轉(zhuǎn)化規(guī)律的研究。

在本研究中，土壤中氮素濃度隨著施肥量的增加而升高，其中硝態(tài)氮濃度的變化更為顯著。有研究表明，施加氮肥顯著促進了硝化作用，且施氮量與土壤硝態(tài)氮含量顯著正相關(guān)，與本文研究結(jié)果一致。然而，施肥量的變化不僅導致土壤中底物質(zhì)量濃度的改變，還會影響土壤孔隙結(jié)構(gòu)和儲水能力，導致根系吸收氮素的差異。本研究表明，盡管作物根系氮素吸收通量增大，吸氮量與施氮量的比值卻呈下降趨勢。可能原因是花生以根瘤固氮為主，其對氮素積累的貢獻率最高可達50%，可有效供給花生生育所需要的大部分氮素營養(yǎng)。所以過量施用氮肥不僅導致花生產(chǎn)量增加不顯著，還降低了肥料利用率。

本研究結(jié)果顯示，雨養(yǎng)農(nóng)田土壤氮素的變化特征與降雨量具有明顯相關(guān)性。由于以灌溉為主的農(nóng)田在降雨事件發(fā)生時，受淹水層緩沖作用后土壤水氮對降雨并不敏感，所以大部分研究集中在對灌溉量的探討上，比如許雪婷等對平原河網(wǎng)區(qū)水稻田的模擬結(jié)果表明，土壤各形態(tài)氮素累積量隨著灌溉量的增加而增加，且硝態(tài)氮較氨氮更易隨土壤水分向土壤深層遷移。而降雨量作為旱作農(nóng)田主要的土壤水分補給，降雨強度影響著土壤表面水分再分配和土壤內(nèi)部向下的人滲水量，進而影響氮素的遷移轉(zhuǎn)化。進一步研究表明，施氮量為225.0 kg·hm-2時，在雨季當總降雨量大于173 mm且降雨強度大于11 mm·d-1時，80 cm土層硝態(tài)氮濃度明顯上升，至少達到1 mg·L-1；在花生生長季后期，天氣較為干旱，當發(fā)生單日強降雨事件，即總降雨量大于135 mm且降雨強度大于135 mm·d-1時，就會導致80 cm土層硝態(tài)氮濃度急速升高，硝態(tài)氮濃度也可達到1 mg·L-1。可能是低雨強長歷時降雨事件主要補充表層土壤水分，降水大部分被用于作物根系利用和蒸散發(fā)損失；而高雨強短歷時的降雨事件雨水勢能大，水壓增強，則會提高土壤入滲速率，雨水能夠攜帶硝態(tài)氮運移到深部土層。這就解釋了高雨強事件能在相對較少的降雨量下，土壤硝態(tài)氮就能運移到80 cm土層。因此，在滿足作物氮需求的前提下，避免在強降雨前施肥，同時做好農(nóng)田排水工作，降低降雨入滲率，有助于減少氮損失。

本研究發(fā)現(xiàn)，隨著降水量的增加，硝態(tài)氮浸出量顯著增加。先前的研究也指出，強降雨顯著影響農(nóng)業(yè)氮循環(huán)，并加速氮素的淋溶損失。雖然紅壤花生地土壤顆粒級配較差，孔隙率高，導致土壤氮素對降雨響應具有一定的滯后效應，但是本研究結(jié)果證明了花生地對紅壤區(qū)地下水硝酸鹽的潛在污染不容忽視。為了減少旱作土壤中的氮損失，并降低地下水面臨的潛在污染威脅，應在區(qū)域范圍內(nèi)優(yōu)化種植策略，有效控制紅壤丘陵區(qū)的非點源污染。

4 結(jié)論

（1）本研究采用HYDRUS-1D模型構(gòu)建了紅壤區(qū)花生地的水氮動態(tài)遷移模型，模型具有較高的模擬精度，模型表現(xiàn)出較高的R2（大于0．87），且均方根誤差和平均絕對誤差均在可接受范圍內(nèi)。因此，模型可用于模擬旱作花生地土壤剖面中氮素的動態(tài)遷移與轉(zhuǎn)化。

（2）土壤中氮素的運移受到施肥量和降水因素的共同影響。施肥量變化主要影響0-40 cm土層內(nèi)氮素濃度變化，而改變降雨量影響0-80 cm土層內(nèi)氮素的遷移轉(zhuǎn)化。銨態(tài)氮和硝態(tài)氮濃度與施肥量和降雨量均呈正相關(guān)，氮素濃度隨土層深度的增加而逐漸下降。

（3）隨著施肥量的增加，作物根系對氮素的吸收通量有所增大，但吸氮量與施氮量的比例卻呈下降趨勢。降雨條件對根系氮素吸收量的影響相對較小。

（4）硝態(tài)氮的淋失是花生種植地的主要氮損失方式。模擬分析表明，降雨和施肥都會影響硝態(tài)氮的淋失量。降雨量的增加加速了氮素向下運移，并顯著提高了硝態(tài)氮的浸出量，而施肥量的改變對硝態(tài)氮累積淋失量的影響較小。因此，在降水量較大且地下水埋藏層較淺的農(nóng)業(yè)種植區(qū)，需要特別關(guān)注硝態(tài)氮淋溶的風險。

（責任編輯：葉飛）

基金項目：第二次全國污染源普查項目（2110399）；國家自然科學基金項目（42107098）；江蘇省研究生科研與實踐創(chuàng)新計劃項目 （SJCX23_0180）