張亞捷，牛海山

（中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院，北京 100049）

氧化亞氮（N2O）作為重要的大氣溫室效應(yīng)氣體之一，其在大氣中濃度的不斷增加將會導(dǎo)致氣候變暖等全球性的環(huán)境問題［1-2］。與二氧化碳（CO2）和甲烷（CH4）等溫室氣體相比，N2O在大氣中存留時間長，達110～150 a，且具有較強的增溫潛勢，1分子N2O的潛在增溫效應(yīng)是1分子CO2的298倍［1-3］。N2O在大氣中雖以痕量存在，但在過去100 a中，其對溫室效應(yīng)的貢獻已達5%～10%。政府間氣候變化專門委員會（IPCC）在最新報告中指出，2011年人類活動造成的輻射強迫已達2.29 W·m-2，其中N2O排放造成的輻射強迫為0.17 W·m-2；大氣中N2O體積分?jǐn)?shù)也已升至0.324 μL·L-1，為80萬年以來最高值，且以每年近0.3%的速率增長，預(yù)計到2050年將達到0.35～0.40 μL·L-1［3-4］。由于N2O較強的增溫效應(yīng)在大氣中所起的重要作用，N2O濃度變化必將對全球生態(tài)環(huán)境產(chǎn)生極大影響。在造成N2O濃度增加的眾多來源中，農(nóng)田土壤是全球最主要的N2O排放源。當(dāng)前全球農(nóng)田土壤N2O的年均排放量（以氮計）高達6.4×1012g，約占全球N2O排放總量的 25%［5-6］。因此，在當(dāng)今氣候變化背景下，研究農(nóng)田土壤田土壤N2O排放已成為環(huán)境科學(xué)領(lǐng)域重要課題，引起學(xué)者廣泛關(guān)注。近年來國內(nèi)外已有大量關(guān)于農(nóng)田土壤N2O排放和相應(yīng)模型發(fā)展的研究報道，筆者將以此為基礎(chǔ)綜合評述農(nóng)田土壤N2O產(chǎn)生機制及相關(guān)模型的發(fā)展，以期深入了解農(nóng)田土壤N2O排放機制，并為發(fā)展相關(guān)模型提供理論依據(jù)。

1 農(nóng)田土壤N2O主要產(chǎn)生機制

農(nóng)田土壤是N2O的主要產(chǎn)生源。硝化作用（包括自養(yǎng)硝化和異養(yǎng)硝化）、反硝化作用（包括微生物自養(yǎng)反硝化、微生物異養(yǎng)反硝化和化學(xué)反硝化）、硝化細(xì)菌反硝化和硝態(tài)氮異化還原成銨作用等都能產(chǎn)生N2O。其中，土壤微生物主導(dǎo)的以銨態(tài)氮（-N）為反應(yīng)底物的硝化作用（主要指自養(yǎng)硝化）、以硝態(tài)氮（-N）為反應(yīng)底物的反硝化作用（主要指微生物異養(yǎng)反硝化）和硝化細(xì)菌反硝化作用是土壤N2O產(chǎn)生的最基本途徑（圖1）。

作為硝化作用的副產(chǎn)物和反硝化作用的中間產(chǎn)物，N2O在此過程中的產(chǎn)生量約占生物圈釋放到大氣中N2O總量的70%～90%［7-9］。硝酸鹽異化還原成銨作用也伴有N2O的釋放；但在自然條件下，此過程中N2O產(chǎn)生量基本可以忽略［10-11］。

圖1 N2O在土壤中的生成與轉(zhuǎn)化［10-11］Fig.1 Formation of N2O in soil as part of the soil nitrogen cycle

在硝態(tài)氮異化還原過程中，除以氣態(tài)氮化物為主要產(chǎn)物的反硝化作用外，還有以為終產(chǎn)物的硝態(tài)氮異化還原成銨（dissimilatory nitrate reduction to ammonium，DNRA）作用。DNRA作用生成N2O排放量占施氮量的1%～8%［10］。僅在某些特定條件下（如高pH值、高C/N比和厭氧環(huán)境），硝態(tài)氮異化還原成銨作用會在土壤氮素轉(zhuǎn)化過程中起較重要作用［11］。

硝化微生物反硝化作用同樣也是N2O產(chǎn)生的主要機制之一，指在低氧條件下僅在硝化微生物驅(qū)動下，NO2-被還原為N2O或N2的過程。該過程分為2個階段：第1階段是將NH4（+NH3）氧化成NO2-，第 2階段是將NO2-還原為NO、N2O或N2；整個過程中沒有NO3-生成，且僅由氨氧化細(xì)菌這類微生物參與［11-12］。高 NH4+-N（或NH3）含量、低有機碳和O2含量以及低pH值環(huán)境更有利于硝化細(xì)菌進行反硝化。

2 預(yù)測農(nóng)田土壤N2O排放的機制過程模型

在農(nóng)業(yè)生態(tài)系統(tǒng)中，土壤N2O等溫室氣體的產(chǎn)生是一個復(fù)雜的過程。在氣候、土壤、現(xiàn)行耕作制度和農(nóng)業(yè)管理措施等諸多條件下，任何因子的細(xì)小變化都可能改變土壤N2O產(chǎn)生及排放。由于這些因子存在地域差異，因此，建立一個受氣候、土壤和農(nóng)業(yè)管理措施綜合影響的農(nóng)田溫室氣體排放模型尤為重要。在明確農(nóng)田土壤N2O排放機制的基礎(chǔ)上，根據(jù)模型建立方法的不同，可將N2O排放模型分為經(jīng)驗統(tǒng)計模型和機制過程模型2類［15-16］。其中，機制過程模型是以土壤N2O產(chǎn)生機制為基礎(chǔ)，認(rèn)為N2O的排放主要來源于以土壤NO3--N和NH4+-N含量為基礎(chǔ)、受控于O2含量和水分狀況等因素的硝化和反硝化作用［15-16］。氣候、土壤、作物和農(nóng)業(yè)等因子通過相互作用影響土壤硝化和反硝化作用，繼而得到用數(shù)學(xué)公式表達的相關(guān)物理、化學(xué)和生物過程中的氮循環(huán)過程，以此模擬土壤N2O排放。機制過程模型的模擬過程復(fù)雜，所需參數(shù)較多，用于點位模擬時準(zhǔn)確度較高；用于區(qū)域模擬時，可探索不同土地利用方式和管理方式對溫室氣體排放的影響，繼而因地制宜地制定合理高效的減排措施［17-21］。

當(dāng)前，用來估算農(nóng)田土壤N2O排放并得到廣泛應(yīng)用的機制過程模型主要包括DNDC（denitrification-decomposition）、WNMM（water and nitrogen management model）、DAYCENT 和 Ecosys（ecosystem model）模型等。各模型的氮循環(huán)過程類似，但側(cè)重點不同，使得各模型N2O占硝化反應(yīng)的比例、反硝化N2O/N2或N2O/NO/N2的產(chǎn)生來源、硝化和反硝化作用臨界狀態(tài)的O2含量和水分狀況等一些基本參數(shù)設(shè)定有所不同，這也使得各模型N2O排放量模擬結(jié)果不同，應(yīng)用和適用范圍也各異［17］（表 1［15，17，22-31］）。當(dāng)前國外對于機制過程模型的研究較多，我國尚沒有自主開發(fā)的機制過程模型。

2.1 DAYCENT模型

DAYCENT模型是Century模型的衍生模型，所含子模塊包括作物生長模塊，作物殘留和土壤有機質(zhì)分解模塊，土壤含水量和溫度模塊，痕量氣體模塊等，主要輸入?yún)?shù)包括每日天氣變化、土壤屬性以及當(dāng)前和歷史土地利用方式等，可用于模擬土壤-大氣-植物碳氮循環(huán)相關(guān)的主要過程，包括作物生長、水分運移、熱量流動、有機碳分解、氮素礦化和固定作用、硝化作用和反硝化作用以及甲烷氧化作用等［22］（圖2）。

DAYCENT模型可用于模擬各種自然生態(tài)系統(tǒng)和人為管理下的農(nóng)業(yè)系統(tǒng)中N2O排放等，在世界范圍內(nèi)尤其是歐美一些國家得到廣泛應(yīng)用，被用于預(yù)測不同土地利用方式或管理措施下農(nóng)田N2O排放，取得了較好的模擬效果［32-34］。

我國學(xué)者應(yīng)用DAYCENT模型的研究相對較少。有研究探討了DAYCENT模型用于預(yù)測農(nóng)田系統(tǒng)不同管理模式下的溫室氣體減排潛力，并對不同管理模式下旱作農(nóng)業(yè)系統(tǒng)和稻作系統(tǒng)的產(chǎn)量和溫室氣體排放量進行模擬［35］。CHENG等［36］則采用國內(nèi)350處農(nóng)田實驗數(shù)據(jù)以驗證DAYCENT模型并預(yù)測中國農(nóng)田生態(tài)系統(tǒng)在不同土地管理模式下的溫室氣體減排潛力，結(jié)果表明，降低土壤有機質(zhì)含量、控制氮肥施用和減少翻耕次數(shù)可大幅減少雨養(yǎng)農(nóng)業(yè)生態(tài)系統(tǒng)下溫室氣體排放量。

2.2 Ecosys模型

Ecosys模型由加拿大阿爾伯達大學(xué)的GRANT等人開發(fā)，可用來模擬不同陸地生態(tài)系統(tǒng)下的生態(tài)行為，是一個較復(fù)雜的數(shù)學(xué)模型［37］。該模型包括7個子模型，主要用以模擬水、熱、碳、氧和氮等在土壤-大氣-植物界面下的遷移轉(zhuǎn)化過程，在模擬過程中同時兼顧土壤管理措施對這些過程的影響［38］。該模型需輸入?yún)?shù)較多，主要有氣候、土壤層、土壤屬性、詳細(xì)的作物信息和管理措施（施肥、耕作、灌溉、播種和收獲）等。

表1 預(yù)測農(nóng)田土壤N2O排放的機制過程模型簡介Table 1 Introduction of major N2O simulation models at the field scale

Ecosys模型結(jié)構(gòu)復(fù)雜，模擬過程繁瑣，模擬時間尺度從1 s到100 a，模擬空間尺度從1 mm到1 km，適用范圍廣，結(jié)果也較準(zhǔn)確［39-40］。但由于Ecosys模型主要為微生物數(shù)學(xué)模型，模擬過程所需參數(shù)較多，因此該模型應(yīng)用受到諸多條件限制，目前僅在加拿大等國家應(yīng)用較多，多用于模擬在不同空間尺度上施肥農(nóng)田N2O排放［41-42］。當(dāng)前，我國尚未開展有關(guān)Ecosys模型的研究。

圖2 DAYCENT模型中的碳氮收支平衡［23-24］Fig.2 The balance of carbon and nitrogen in DAYCENT model

2.3 WNMM模型

WNMM模型即農(nóng)業(yè)水氮管理模型，由澳大利亞墨爾本大學(xué)Li Yong與中國科學(xué)院南京土壤研究所張佳寶等人共同開發(fā)，用于模擬不同農(nóng)業(yè)管理模式下（灌溉、耕作和施肥等）水分遷移、溶質(zhì)運動、作物生長和碳氮循環(huán)等過程［28，43］（圖3）。WNMM模型適用于干旱半干旱氣候條件下的旱地作物開發(fā)，結(jié)構(gòu)較為簡單，所需輸入?yún)?shù)較少，著重于水肥管理措施對溫室氣體排放的影響，主要應(yīng)用于澳大利亞和中國等國家［15，44］。

WNMM模型在國際和國內(nèi)的應(yīng)用多集中在不同氣候條件及農(nóng)業(yè)管理措施下的農(nóng)田N2O等溫室氣體排放及作物產(chǎn)量的相關(guān)研究上，且均取得了較好的模型擬合結(jié)果［15］。LI等［45-46］在澳大利亞西部雨養(yǎng)麥田中應(yīng)用WNMM模型研究了不同氮肥施用所引發(fā)的N2O排放變化，并提出了適用于當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)的N2O減排措施。在國內(nèi)相關(guān)研究中，WNMM模型被用于模擬青藏高原陸地生態(tài)系統(tǒng)N2O排放總量，發(fā)現(xiàn)硝化作用造成了絕大部分N2O排放損失［47］。在發(fā)展了新的含氮氣體排放子模型后，WNMM還可用于模擬中國亞熱帶茶田N2O排放損失，發(fā)現(xiàn)反硝化作用為主要N2O排放來源［48］。

圖3 WNMM模型在農(nóng)田系統(tǒng)中的氮循環(huán)示意［44］Fig.3 Diagram of nitrogen cycling in soil-crop system in WNMM

2.4 DNDC模型

DNDC模型由美國新罕布什爾州大學(xué)陸地海洋空間研究中心李長生等人研發(fā)，最初被用來模擬美國農(nóng)業(yè)土壤N2O排放。模型輸入?yún)?shù)主要包括氣象條件、土壤屬性、植被和農(nóng)田管理措施等4個方面。在研發(fā)初期，DNDC模型包括土壤氣候過程子模型、有機質(zhì)分解過程子模型和反硝化作用子模型3個子模型，并對植物生長過程子模型和土地管理利用方式進行解釋［31］。隨后，DNDC模型增加了利用土壤氧化還原電位模擬發(fā)酵過程的子模型并進一步修正了N2O等氣體排放的模擬過程等，逐漸發(fā)展為現(xiàn)在的DNDC模型［49］（圖4）。

DNDC模型在模擬N2O排放方面的應(yīng)用非常廣泛。各國研究者應(yīng)用DNDC模型模擬不同作物管理模式下農(nóng)業(yè)土壤N2O等溫室氣體排放，取得了較好擬合效果，并發(fā)現(xiàn)N2O排放對施肥時機、施肥頻率、肥料類型、施肥方式和使用硝化抑制劑（或脲酶抑制劑）等因素都非常敏感［50-51］。我國學(xué)者也應(yīng)用DNDC模型做了大量相關(guān)工作。起初，DNDC模型被應(yīng)用于貴州省農(nóng)業(yè)土壤N2O排放研究中，初步論述了我國亞熱帶旱田生態(tài)系統(tǒng)N2O排放特征及影響因素，揭示并討論了施肥和耕作措施變化等對該地區(qū)農(nóng)業(yè)土壤N2O排放量的潛在影響［52］。近些年，關(guān)于DNDC模型的研究重點則逐漸集中在不同研究區(qū)農(nóng)業(yè)生態(tài)系統(tǒng)下N2O等溫室氣體的排放機制和減排措施，發(fā)現(xiàn)合理的水肥管理系統(tǒng)能有效減少N2O排放；DNDC模型還被廣泛應(yīng)用于碳氮循環(huán)模擬研究中，擬合效果均較理想［53-54］。

圖4 DNDC模型結(jié)構(gòu)［18］Fig.4 Structure of DNDC model

適用于中國特有農(nóng)業(yè)生態(tài)系統(tǒng)的China-DNDC也在不斷發(fā)展中。China-DNDC修改土壤水氮運移過程的缺省參數(shù)，引入地表徑流曲線和修正的通用土壤流失方程來控制和模擬地表徑流，還加入有關(guān)薄膜覆蓋管理模式的參數(shù)化模塊，并補充種養(yǎng)結(jié)合的相關(guān)子模型。這些改動提升了模型在中國主要農(nóng)業(yè)生態(tài)系統(tǒng)及生產(chǎn)管理模式下的應(yīng)用能力，并取得了理想的模擬效果［55-56］。與其他模型相比，DNDC模型適用環(huán)境范圍廣（旱地和水田等），模擬過程明確，模擬效果好，便于改造成適用于其他各種用途的綜合模型；且在利用GIS技術(shù)進行區(qū)域N2O排放模擬時，DNDC模型優(yōu)勢更明顯［18，57-58］。

3 結(jié)論與展望

近年來大氣中N2O濃度呈不斷上升趨勢，其引發(fā)的環(huán)境問題已成為當(dāng)今重要的全球性課題。N2O主要來源于土壤尤其是農(nóng)田土壤。隨著土壤N2O產(chǎn)生機制研究的不斷深入，影響農(nóng)田土壤N2O排放的相關(guān)模型研究也在發(fā)展，研究廣度和深度在不斷加強。機制過程模型可有效模擬陸地生化循環(huán)中的碳氮時空動態(tài)；但不同模型有著不同的物理生化過程，這也使各模型模擬N2O等溫室氣體排放過程時表現(xiàn)出一定差異。在實際應(yīng)用中，使用者應(yīng)根據(jù)研究目的和數(shù)據(jù)可用性來選擇模型。針對當(dāng)前研究現(xiàn)狀，還應(yīng)加強以下幾方面的探索：

在N2O產(chǎn)生機制研究方面，應(yīng)更多采用新的工具和技術(shù)（如同位素技術(shù)、宏基因組學(xué)等），加強土壤微環(huán)境碳氮循環(huán)機制與微生物群落參與的N2O產(chǎn)生與消耗過程的研究，尤其應(yīng)加強硝化細(xì)菌反硝化和硝態(tài)氮異化還原成銨過程的相關(guān)研究［10-11，59］。

在模型開發(fā)和應(yīng)用相關(guān)研究上，應(yīng)重視以下方面：（1）模型模擬N2O排放機制研究中，普遍只考慮硝化和反硝化作用對N2O排放的貢獻；對硝化細(xì)菌反硝化和硝態(tài)氮異化還原成銨過程的模擬研究工作應(yīng)繼續(xù)深入開展。（2）作物修剪、病蟲害發(fā)生和雜草生長等影響農(nóng)作物生長并作用于土壤N2O排放的相關(guān)模塊仍有待開發(fā)和改進。（3）模型中有關(guān)農(nóng)田管理措施的描述應(yīng)更詳盡，這將有助于探索更有效的土壤N2O減排措施（如施用控釋肥料，肥料深施或混施，應(yīng)用滴灌等節(jié)水灌溉技術(shù)，應(yīng)用秸稈還田和秸稈覆蓋技術(shù)以及使用硝化抑制劑等）［5，8，60-61］。（4）在中國，還應(yīng)加強其他N2O排放模擬模型（或模擬N2O排放的作物模型）相關(guān)研究，如NOE（the algorithm of nitrous oxide emission）、NGAS、Expert-N、NLOSS、STICS（simulateur multidiscplinaire pour les cultures standard）和 DSSAT（decision support system for agrotechnology transfer）等模型［62-67］。與N2O排放相關(guān)的以上模型的對比研究可參見文獻［68-69］。在對這些模型的研究基礎(chǔ)上，可更多采用多模型結(jié)合的方法估算農(nóng)田土壤N2O排放，以提高預(yù)測的可靠性與準(zhǔn)確性［70-71］。（5）在未來研究中，可用模型估算區(qū)域、國家乃至全球農(nóng)田土壤N2O排放總量，進一步研究不同氣候條件、土壤條件和農(nóng)田管理措施等對農(nóng)田土壤N2O排放的影響，這將有助于決策者制定更加合理有效的溫室氣體減排措施，維持農(nóng)業(yè)穩(wěn)定發(fā)展［1，72-73］。