氮輸入對中國東北地區土壤碳蓄積的影響

顧峰雪，黃 玫，張遠東，李 潔，閆慧敏，郭 瑞，鐘秀麗

1 中國農業科學院農業環境與可持續發展研究所，農業部旱作節水農業重點實驗室，北京 100081 2 中國科學院地理科學與資源研究所生態系統觀測與模擬重點實驗室，北京 100101 3 中國林業科學研究院森林生態環境與保護研究所，國家林業局森林生態環境重點實驗室，北京 100091

?

氮輸入對中國東北地區土壤碳蓄積的影響

顧峰雪1，黃 玫2，張遠東3,*，李 潔1，閆慧敏2，郭 瑞1，鐘秀麗1

1 中國農業科學院農業環境與可持續發展研究所，農業部旱作節水農業重點實驗室，北京 100081 2 中國科學院地理科學與資源研究所生態系統觀測與模擬重點實驗室，北京 100101 3 中國林業科學研究院森林生態環境與保護研究所，國家林業局森林生態環境重點實驗室，北京 100091

由于人類活動的干擾，通過沉降和施肥形式進入陸地生態系統的氮素持續增加，中國已經成為繼歐洲和北美之后的第三大氮沉降區，同時也是最大的化肥消費國。氮輸入與陸地生態系統生物地球化學循環的一系列過程都相互聯系，碳循環及其格局也受到氮輸入的影響。土壤有機碳庫在全球碳循環中具有重要作用，氮輸入能否或在多大程度上對土壤碳庫產生影響已經成為全球變化和氮沉降研究中不可回避的問題。東北地區是世界三大黑土帶之一，土壤碳的變化不僅對于土壤肥力維持具有重要意義，而且對區域碳收支具有重要影響。利用生態系統過程模型——CEVSA2模型，基于我國能源消費、施氮數據和降水數據生成了一套中國大氣氮沉降的時空網格數據，結合大氣CO2濃度、氣候、土地覆被、土壤類型和質地的時空數據，模擬評估了1961—2010年氮輸入對中國東北地區土壤碳蓄積的影響。結果表明：(1)1961—2010年東北地區的平均氮沉降速率為1.00 gN m-2a-1，年增長率為0.047 gN m-2a-1。東北農田總氮輸入速率達到5.78 gN m-2a-1，從20世紀80年代開始顯著增加。(2)氮輸入促進了東北地區土壤碳的蓄積，東北陸地生態系統的土壤碳密度平均增加了135 gC/m2，50a氮輸入共增加土壤碳蓄積0.16 PgC。(3)氮輸入引起的東北地區土壤碳蓄積量的變化呈現出東高西低、南高北低的空間格局，遼河平原、松嫩平原和三江平原的土壤碳密度增加量超過了300 gC/m2。(4)不同植被類型下的土壤碳密度對氮輸入的響應存在較大差異，農田土壤碳密度平均增加了230 gC/m2，森林、灌叢和草地則分別增加了76、169 gC/m2和89 gC/m2。氮輸入的空間差異和不同植被類型對氮輸入響應的差異共同決定了東北地區土壤碳增加量的空間格局。通過本研究闡明了氮輸入對東北農田土壤碳蓄積的影響，從而為農田生態系統的固碳減排和農田土壤碳氮管理提供了決策依據。

氮沉降；施肥；土壤碳蓄積；東北；CEVSA2模型

目前人類活動已經極大地改變了氮循環過程[1-2]。1860年以來，全球活性氮生產量逐年增加。據估計，全球人為活動導致的活性氮生產量由1860年的15 TgN/a提高到2000年的165 TgN/a[1]，其后這一數值仍在持續上升[3-4]。農業活動是重要的活性氮來源，工業固氮速率的增加使農業生產中日益增加的氮需求得以維持，從1910開始，工業固氮產生的活性氮由零的數量開始增加，到2000年工業固氮產生的活性氮已經超過的100 TgN/a，其中85%來自于化肥的生產[1,3]。我國作為世界上最大的發展中國家，每年消耗的肥料氮大于24 Tg，大約是全世界肥料氮施用量的30%左右[5]。這些活性氮通過沉降和施肥的方式進入陸地生態系統中，會引起土壤酸化，生物多樣性減少，改變生態系統生產力，引起一系列的生態問題[6]。氮輸入的增加引起了科學家和公眾的廣泛關注[1,2,7]。

土壤有機碳庫在全球碳循環中具有重要作用。據估計，土壤有機碳庫大約是大氣碳庫的3倍多，是植物碳庫的2.5—3倍[8]。土壤碳庫的動態平衡影響著作物產量和土壤肥力的高低，是土壤肥力保持和提高的重要研究內容[9]。氮輸入(氮沉降和施氮)是影響陸地生態系統土壤呼吸過程的重要因素[10-11]。氮沉降能否或在多大程度上對土壤碳庫產生影響已經成為全球變化和氮沉降研究中不可回避的問題[12]。國內外學者也開展了一些氮輸入與碳蓄積的模型研究,如，Wei等[13]基于FORECAST模型，研究了中國東南部杉木人工林碳蓄積對氮沉降的響應，模擬結果表明，在現有氮沉降水平(5—10 kg hm-2a-1)下，中國東南地區杉木林蓄積碳7.4×106MgC/a，其中土壤碳4.9×106MgC/a，如果氮沉降水平增加到7.5—15 kg hm-2a-1，蓄積的碳量將達到16×106MgC/a，土壤碳蓄積將達到9.9×106MgC/a；呂超群[14]應用動態陸地生態系統模型(DLEM)的研究表明，從1901—2005年，大氣氮沉降增加及其與其他環境變化的協同作用使得陸地碳庫增加6.43 PgC，其中植被、土壤和凋落物分別占48%，45%和7%，并且指出，長期持續的氮沉降增加必然會提高土壤碳庫在整個氮促碳庫中的重要性；Jain等[15]應用ISAM-NC和ISAM-C模型的研究表明，20世紀90年代，氮輸入(62.2 TgN/a)增加引起全球陸地碳匯增加0.26 PgC/a。盡管在全球和區域尺度上，開展了許多氮輸入與土壤碳蓄積響應的模擬研究，但關于氮輸入引起的土壤碳蓄積的估算仍存在很大的不確定性[16]，缺乏氮輸入對土壤碳蓄積影響的深入研究，因此，有必要對氮輸入對土壤碳蓄積的影響進行深入探討。

東北地區處于我國最高緯度地區，是中國最大的天然林區和重要的商品糧生產基地，也是世界三大黑土帶之一[17- 19]。農田作為東北地區重要的土地利用方式，其土壤碳的變化不僅對區域碳收支具有重要影響，而且對于土壤肥力維持具有重要意義。許多學者對我國東北地區的土壤碳的時空格局及其變化特征進行了研究[20- 25]，而且，研究表明氮輸入對該區土壤碳蓄積具有重要作用，如對東北黑土的長期定位試驗研究表明，施用化肥有利于黑土有機質的積累[26]，且土壤有機質含量對施肥的響應隨土層深度的加深呈下降趨勢[27]。但是氮輸入對東北地區土壤碳蓄積時空格局影響的研究較少。開展氮輸入對東北地區土壤碳蓄積的影響研究，尤其是農田的土壤碳蓄積，對于準確評估該區域的碳收支狀況，估算該區域土壤碳蓄積量以及農田土壤碳氮管理都具有重要意義。

1 研究方法

1.1 研究區域簡介

本研究中，東北地區(北緯 38°48′—53°33′N，東經 115°31′—135°05′E)主要包括黑龍江、吉林、遼寧以及內蒙古地區東北部，北部與俄羅斯接壤，東南部與朝鮮半島相接，南部濱臨中國渤海和黃海[28]。本研究中統計的總面積為11.73×105km2。東北地區屬寒溫帶大陸性季風氣候，四季分明，冷期較長，1961—2010東北地區年平均溫度在2.49—6.02 ℃之間，隨著緯度升高，年均溫逐漸減低，年降水量在427—680 mm，從東南向西北，降水量逐漸減少。

東北地區包括“三山兩平原一高原”，即西部為大興安嶺和呼倫貝爾高原，北部為小興安嶺，東部有長白山，中部為松遼平原，東北部為三江平原。東北平原的海拔在200 m左右，分布著大量黑土，面積7萬多km2，是我國重要的商品糧基地，水稻、玉米、大豆等作物播種面積均占全國首位[29]。黑土的土壤結構良好，土層疏松軟綿，呈中性或微酸性，有機質含量高，腐殖質較多，肥力較高，理化性能好，素有“土中之王”的美稱。東北陸地生態系統的土壤有機碳庫儲量為26.43 PgC，在全球碳庫預算中起著舉足輕重的作用[30]。

化肥是東北地區農業賴以發展的重要物質基礎，東北地區的化肥消費以氮肥為主，2010年全國單位面積施肥量為346.1 kg/hm2，其中東北地區的遼寧為343.9 kg/hm2，吉林為350.1 kg/hm2，黑龍江為176.3 kg/hm2，內蒙古為253.1 kg/hm2，近年來東北地區化肥施用量呈緩慢提升的趨勢[31-32]。

1.2 數據來源

1.2.1 氣候數據

模型中所使用的氣象數據來自于國家氣象信息中心，包括1961—2010年全國756個氣象臺站的每旬平均氣溫、降水量、云量和相對濕度。使用ANUSPLIN4.1插值軟件的樣條函數插值法對氣象數據進行內插到10km×10km的空間網格上[33]，得到全國1961—2010年的0.1°氣象柵格數據。

1.2.2 CO2,土壤類型和土壤質地,植被類型數據

大氣CO2濃度資料來源于美國夏威夷Mauna Loa觀測所(http://co2now.org/Current-CO2/CO2-Now/noaa-mauna-loa-co2-data.html)。土壤類型和質地資料取自1∶14000000土壤類型圖和第二次土壤普查數據，對其進行數字化并重采樣到0.1°。植被數據來自于Global Land Cover 2000數據庫(European Commission, Joint Research Centre, 2003. http://bioval.jrc.ec.europa.eu/products/glc2000/glc2000.php)共有22種土地覆被類型，在ArcGIS中進行插值重采樣，匹配于氣象數據的分辨率。

1.2.3 氮輸入相關數據：施肥、能源消費和氮沉降

本研究中的氮輸入主要包括通過沉降和施肥兩種途徑進入生態系統的氮。施肥數據是由統計年鑒的數據計算得到，氮沉降數據基于CEVSA2模型，根據柵格化的施肥和能源消費等數據在模型中運行得到。

施肥數據一方面是模型的氮輸入數據，同時也是計算氮沉降數據的關鍵變量。施肥數據來源于中國和各省的統計年鑒(NBS, http://www.stats.gov.cn/)，包括1978—2008年全國和各省的施肥量、施氮量，以及全國2000多個縣1992、1995、1999年的施肥數據，根據模型輸入數據的要求，對施肥數據做以下處理：(1)按照不同時期各縣所在省施肥總量的比例以及所在省1978—2008年的變化趨勢，計算得到全國各縣級單元施氮總量1978—2008年的變化趨勢。(2)基于1978—2008年各縣級單元施氮總量的線性變化趨勢，計算得到各縣級單元1961—1977年和2009、2010年的施氮總量數據。(3)將縣級單元施氮數據屬性表與中國縣市界先掛接，形成屬性空間化。(4)利用通過行政區劃求取的2000年各縣市面積的和，然后基于縣市名稱與各縣市耕地面積掛接，實現耕地面積總量的空間化。(5)最后利用施肥量/各省市縣耕地面積和的柵格數據，實現施肥量單元柵格化(gN m-2a-1)。

能源消費數據來源于中國統計局 (NBS, http://www.stats.gov.cn/)，查找自1978—2008年全國和各省的能源總量(萬噸標準煤/年)，全國和各省1961—1977年以及2009、2010年的能源消費數據根據1978—2008年的線性變換趨勢計算得到。根據省界圖，統計各省的面積，將各省的能源消費總量平均到單位面積上，在ArcGIS中提取模型輸入所需的0.1°格點的能源消費量(g標準煤m-2a-1)。

利用柵格化后的施肥、能源消費和降水對氮沉降的時空變化進行模擬，

Ndepo=Ammdepo+Nitdepo

(1)

Ammdepo=PRE×P1×fert+P2×fert

(2)

Nitdepo=PRE×P3×ener+P4×ener

(3)

式中，Ndepo為大氣氮沉降速率(gN m-2a-1)，Ammdepo(gN m-2a-1)為以NHx形式沉降的氮，Nitdepo(g N m-2a-1)為以NOy形式沉降的氮，PRE為降水量(mm)，fert為施氮量(gN m-2a-1)，ener為能源消費量(g標準煤m-2a-1)，P1、P2、P3和P4為擬合參數。通過文獻收集了中國區域內163個站點不同年份觀測的588個氮沉降數據(包含濕沉降、干沉降或干濕沉降均有的數據)，利用這些觀測數據，對P1—P4參數進行了擬合。該方法不僅具有較好的理論基礎，輸入數據和參數易于獲取，同時可以根據能源利用效率和氮肥利用率的變化來調整參數，從而實現對未來情景的預測。

1.2.4 CEVSA2模型簡介和運行

CEVSA2模型是一個基于生理生態過程模擬植物-土壤-大氣系統能量交換和水碳氮耦合循環的生物地球化學循環模型。CEVSA2模型基于目前已有的機理發現，包含了氮對光合、呼吸、分配和土壤碳分解等所有過程的影響模擬，從而能夠很好地表達氮輸入變化對于碳循環過程的影響。在空間模擬過程中輸入數據和參數易于獲取且空間分辨率較高。

本研究中，模型運行的時間分辨率為10d，空間分辨率為0.1°。首先應用1961—2010年的平均氣候數據運行模型至生態系統平衡態，即各個狀態變量如植被、土壤碳貯量以及土壤含水量等年際變化量小于0.1%，且凈初級生產力(NPP)、凋落物產生量 (LT)與異養呼吸(HR)相等，然后用1961—2010年每旬資料進行動態模擬，并反復運行模型以消除假定的生態系統狀態變量的初始值(即平衡態假設)對模擬結果的影響。從最終模擬結果中提取中國東北地區的數據進行分析。

2 模型驗證

本研究中，為了驗證CEVSA2模型在土壤碳儲量對氮輸入響應方面的模擬效果，利用從文獻中收集的公主嶺[34]和海倫[35]兩個站點施肥的長期試驗數據對其進行驗證。利用CEVSA2模型，分別計算兩個站點在16.5 gN/m2和12.8 gN/m2(公主嶺和海倫長期試驗中的施氮量)氮輸入水平下土壤有機碳密度，其與無氮條件下土壤碳密度的差值作為氮輸入引起的土壤碳儲量變化的模擬值。模擬值與站點實測值的比較如圖1。氮輸入引起的公主嶺地區土壤碳儲量變化13a的試驗值為 0.177 kgC/m2,模擬值為0.204 kgC/m2，相對誤差為 15.39%。氮輸入引起的海倫地區土壤碳儲量變化的19a平均值為0.300 kgC/m2，模擬值為0.280 kgC/m2，相對誤差為6.84%，因此，可以看出CEVSA2模型可以很好地模擬東北地區農田土壤碳對氮輸入的響應。模型模擬值與試驗值之間的誤差主要是由于：(1)土壤的空間異質性，土壤采樣點不一樣，會造成土壤實測值存在較大的變異性，通過統計分析發現，兩地的組內方差 (公主嶺為0.017，海倫為0.010) 均大于組間方差(公主嶺為0.003，海倫為0.009);(2)模型模擬的誤差，試驗中，氮肥非純氮，氮與其他肥料元素之間存在著復雜的交互作用，而模型則只考慮了氮的作用。

圖1 氮輸入對土壤碳蓄積影響的模擬值與實測值的比較 Fig.1 Comparison between the estimated C storage induced by N input with CEVSA2 and that from field measurements in Gongzhuling and Hailun station

3 結果與討論

3.1 東北地區氮輸入的時空變化

圖2 1961—2010年中國東北地區氮沉降的年際變化 Fig.2 The inter-annual variation of N deposition in Northeast China during 1961—2010

研究時段內，我國東北地區的大氣氮沉降速率平均為1.00 gN m-2a-1，高于全國的氮沉降速率0.81 gN m-2a-1，2010年，東北地區的氮沉降速率達到2.80 gN m-2a-1，是同時段全國平均氮沉降速率的1.64倍。東北地區氮沉降速率較高，與東北地區是我國的老工業基地和重要的農業基地有關，工業生產和氮肥施用引起高氮沉降。東北地區的氮沉降速率的年增長率為0.047 gN m-2a-1，高于全國氮沉降速率的增速0.036 gN m-2a-1(圖2)，20世紀60年代和70年代氮沉降速率基本保持不變，80年代開始增加，90年代以后氮沉降速率增加更加顯著。

從空間上看，黑龍江、吉林和遼寧地區的氮沉降速率高于內蒙古的東北部地區(圖3)。從1961至2010年，內蒙古東北部地區的氮沉降速率基本不變，始終在0—0.25 gN m-2a-1的范圍，而在黑龍江的西南部、三江平原、吉林省的中西部和遼寧地區的氮沉降速率發生了顯著變化，從20世紀60年代的0.25—0.75 gN m-2a-1，增加到80年代的0.75—2.0 gN m-2a-1，再到21世紀初的2.0—4.5 gN m-2a-1。

圖4 1961—2010年中國和東北地區施氮量的年際變化 Fig.4 The inter-annual variation of N fertilizer in China and Northeast China during 1961—2010

1961—2010年中國和東北地區的施氮量均呈增加趨勢(圖4)。我國東北地區的施氮量低于全國平均水平，變化范圍在1.72—13.82 gN m-2a-1，平均施氮量為4.78 gN m-2a-1。就變化趨勢來看，1961—2010年，東北地區每年增加0.21 gN m-2a-1。東北地區的施氮量20世紀70年代至21世紀初分別比60年代增加了15%、89%、250%和456%，由此可見，施氮量的增加主要在20世紀80年代以后，80年代以前我國的施氮量水平很低而且增量較小。

為了說明東北地區的施氮量在空間上的差異及其變化，本研究對1978—2008年施氮量統計資料進行空間化處理，得到了近31年施氮量的空間分布，并分析了20世紀80年代、90年代和21世紀初東北施氮量的空間分布(圖5)。近31年三江平原地區的農田施氮量普遍較低，在0—5 gN m-2a-1，松遼平原地區的施氮量普遍較高，在5.0 gN m-2a-1以上，部分地區的施氮量超過10 gN m-2a-1。近30年來，我國東北地區的農田施氮量持續增加。申建波等[36]對黑龍江稻作區的施肥情況研究發現，該區氮肥施用量從3.5—35.36 gN m-2a-1不等，平均為15 gN m-2a-1，與本研究結果非常相近。我國東北農田的施氮量低于全國平均水平，主要與東北地區農田施氮量分布差異較大有關，而且東北地區存在大面積的高施氮區，因此，有必要對東北地區農田施氮量的空間差異引起的生態效應進行進一步研究。

圖5 1981—2008年東北農田平均施氮量的空間分布和不同年代際的空間分布Fig.5 Spatial pattern of N fertilizer rate during 1981—2010, 1981—1990, 1991—2000 and 2001—2008

圖6 1961—2010年氮輸入引起的東北地區土壤碳密度的變化量 Fig.6 The inter-annual variation of soil C density induced by changes of N input during 1961—2010

氮肥施用量的增加，與氮肥工業的發展有直接關系。1965—1979年，碳酸氫銨的生產使中國氮肥產量迅速增加，是主要的氮肥消費形式。20世紀80年代以后，隨著無煙煤生產尿素技術的發展，氮肥產量迅速增加，1997年高濃度尿素超過碳酸氫銨成為中國第一大氮肥品種，氮肥消費迅速增長[37]。根據欒江等[31]的研究，化肥使用強度的增長是中國化肥施用總量增長的主因，但從2007年以后，使用強度的貢獻不斷下降，播種面積調整的貢獻有所提高。中國肥料的當季利用率低，氮肥利用率僅為30%—35%，氮肥的揮發、淋溶和徑流損失巨大，隨之引起土壤肥力下降，農作物品質降低和環境污染等[38]。從本研究的研究結果來看，我國農田施氮量持續增加，不同區域農田施氮量的差異很大，因此針對不同地區施氮量的差異，對我國農田施氮量進行管理，對于提高農作物產量和品質以及改善土壤肥力都將具有重要意義。

3.2 東北地區土壤碳儲量的時間變化對氮輸入的響應

本研究中，利用有氮輸入和無氮輸入兩種模擬情景中土壤碳密度的差值作為氮輸入增加引起的土壤碳密度變化量，并將2010年土壤碳蓄積的變化量作為1961—2010年由于氮輸入增加(氮沉降和施肥)引起的土壤碳蓄積增加量。由圖6可知，氮輸入速率的增加使得東北地區的土壤碳密度的增加量逐年增大。20世紀60年代，東北地區氮輸入增加使得東北地區的土壤密度增大了16 gC/m2，而在隨后的70年代，氮輸入增加引起的土壤碳密度增加量減小為10 gC/m2，從80年代開始，氮輸入引起的土壤碳密度的增加量逐漸增加，這與東北地區氮輸入的增長趨勢相一致，80年代以后，人類活動的干擾以及氮肥工業的發展，東北地區的氮沉降和施肥水平都大幅增加(表1)。研究時段內，東北地區的氮輸入量由1961年的0.97 gN m-2a-1增加到2010年的7.6 gN m-2a-1，使得東北陸地生態系統土壤碳密度增加135 gC/m2，土壤碳匯增加0.16 PgC。

表1 1961—2010年氮輸入引起的東北地區土壤碳密度增加量的年代際變化

氮輸入增加引起的土壤碳密度增加量在1961—1970年等于1970年的氮輸入引起的土壤碳密度的變化量，1971—1980等于1980年的氮輸入引起的土壤碳密度的變化量與1970年的差值，1981—1990等于1990年的氮輸入引起的土壤碳密度的變化量與1980年的差值，依次類推土壤碳庫的增減趨勢取決于碳輸入與輸出之間的平衡關系，一方面氮輸入增加導致葉氮濃度的升高[39- 41]，由此促進光合作用從而增加植物生物量。但氮輸入對凋落物的影響程度取決于樹種和養分條件[12]，一些研究發現氮輸入顯著增加了凋落量[42- 44]，但另外一些研究表明氮輸入對凋落物量影響不顯著[41,45]。在本研究中，氮輸入增加了植物的凋落物量碳，凋落物量平均每年增加0.21 gC m-2a-1，尤其是從20世紀80年代以后，氮輸入使得東北地區的凋落物量急劇增加(圖7)，從而導致進入土壤的碳素增加。另外，目前關于氮輸入對土壤呼吸的影響還沒有一致的結論，可能存在促進作用[46-49]、抑制作用[48,50-51]或無明顯影響[52]，同樣取決于樹種、氮沉降濃度和養分條件。本研究中氮輸入增加了土壤的異養呼吸，平均每年增加0.13 gC m-2a-1。植物根系的周轉、凋落物和土壤有機質的分解與土壤碳庫緊密相關，凋落物和土壤碳氮比與分解作用有著緊密聯系[53- 55]。總得來說，森林凋落物的分解速率因凋落物的種類和環境條件不同而不同，并且不同氣候條件和不同分解階段影響凋落物分解的各因子重要性可能存在差異。本研究結果表明氮輸入引起的植物凋落物碳輸入的增量大于土壤異養呼吸的增量，從而導致土壤碳儲量的上升。

圖7 1961—2010年氮輸入引起的東北地區土壤異養呼吸和凋落物量的變化Fig.7 Temporal variations of heterotrophical respiration and litter production change induced by N input

3.3 東北地區土壤碳蓄積的空間格局對氮輸入的響應

氮輸入引起的東北地區土壤碳密度的變化量呈現出東高西低，南高北低的空間格局(圖8)。研究時段內，氮輸入的增加使得土壤碳密度減少的面積僅占研究區域總面積的16%，主要分布在小興安嶺地區，大興安嶺和長白山的小部分地區。在絕大部分地區，氮輸入增加了東北地區的土壤碳密度。其中，土壤碳密度的增加量在0—100 gC/m2的地區占研究區域總面積的43%，主要分布在大興安嶺、呼倫貝爾高原、內蒙古高原等地區；土壤碳密度的增加量大于100 gC/m2的地區占研究區域總面積的41%，其中土壤碳密度增加量大于300 gC/m2的面積占研究區域總面積的19%，主要分布在松遼平原和三江平原地區，以及黑龍江的西南部地區。總體來看，東北地區的氮輸入使得土壤碳密度減少或者增加量較小的地區主要是在東北森林和草地的分布區，土壤碳蓄積量增加較多的地區主要在農田分布區。

圖8 1961—2010年氮輸入引起的東北地區土壤碳密度變化量的空間格局Fig.8 Spatial pattern of soil C density change induced by N input in Northeast China during 1961—2010

隨著氮輸入量的不斷增加，土壤碳蓄積的響應幅度不是維持不變的，其年代際的空間格局對于了解生態系統行為具有重要意義。圖9通過對比整個研究時段內不同年代的氮效應差值來揭示氮輸入增加引起的土壤碳密度變化量的年代際差異。從空間分布上來看，20世紀60年代和70年代氮輸入引起的土壤碳密度的變化不大，從80年代開始，氮輸入使得遼河平原、吉林中部和黑龍江西南部、三江平原等的部分地區土壤碳密度的增加量達到75—200 gC/m2，90年代達到75—200 gC/m2的區域面積進一步擴大，部分地區開始超過200 gC/m2，2000年后75—200 gC/m2的范圍在這些地區逐漸向外擴展，土壤碳密度的增加量進一步增多。而在大、小興安嶺、長白山以及內蒙古地區，氮輸入對這些地區的土壤碳密度的影響幾乎沒有變化，直到20世紀90年代，這些地區的土壤碳密度增加量才逐步達到10—50 gC/m2。

圖9 不同年代氮輸入引起的東北地區土壤碳密度變化量的空間格局Fig.9 Spatial pattern of soil C density change induced by N input in different decade in Northeast China

造成東北地區土壤碳密度量對氮輸入響應的空間格局存在差異的原因主要有兩方面：(1)東北地區氮輸入存在空間差異。研究發現，氮輸入引起的東北地區土壤碳密度變化量的空間格局與東北地區氮輸入的空間格局具有很好的一致性，如松遼平原和三江平原地區的氮沉降速率不論是50a平均值，還是不同年代的平均值始終都是東北地區氮沉降速率最高的地區，而且這些地區主要是農田分布區，有大量氮素以氮肥的形式施入生態系統中，進而影響植物生長，改變凋落量和土壤呼吸過程，因此這些地區由氮輸入引起的土壤碳密度增加量及其增速也同樣是最高的。而在氮沉降速率較小的大小興安嶺和內蒙地區的土壤碳密度的變化量同樣很小(圖9)。(2)不同植被類型對氮輸入的響應存在差異(表2)。從氮輸入對東北地區的土壤碳密度影響的空間格局及其不同年代土壤碳蓄積空間格局的影響來看，土壤碳密度變化較小的地區主要是在東北森林和草地的分布區，土壤碳密度增加較多的地區主要在農田分布區，由此可以初步認為氮輸入對土壤碳密度的影響與植被類型有很大的關系，尤其是農田的土壤碳密度受氮輸入影響很大。

表2 不同植被類型下的氮輸入使土壤有碳蓄積增加量

4 結論

(1)過去50a東北地區的氮輸入顯著增加，其中大氣氮沉降速率平均為1.00 gN m-2a-1， 2010年達到2.80 gN m-2a-1，年增長率為0.047 gN m-2a-1，高于全國的平均水平；而農田的施氮量則由1961年的1.721 gN m-2a-1增加到2010年3.82 gN m-2a-1的。東北平原地區是氮沉降和農田施氮量增加最為顯著的區域。

(2)氮輸入促進了東北地區土壤碳的吸收，東北陸地生態系統的土壤碳密度平均增加了135 gC/m2，1961—2010年氮輸入共增加土壤碳蓄積0.16 PgC。

(3)不同植被類型下的土壤碳密度對氮輸入的響應存在很大差異，氮輸入引起農田、灌叢、草地、森林的土壤碳密度分別增加230、169、89 gC/m2和76 gC/m2。

(4)氮輸入的空間格局和不同植被對氮輸入的響應差異決定了東北地區土壤碳蓄積變化的空間特征。氮輸入引起的東北地區土壤碳蓄積量的變化呈現出東高西低、南高北低的空間格局，東北平原和三江平原的土壤碳密度增加量超過了300 gC/m2。

[1] Galloway J N, Aber J D, Erisman J W, Seitzinger S P, Howarth R W, Cowling E B, Cosby B J. The nitrogen cascade. BioScience, 2003, 53(4): 341- 356.

[2] Vitousek P M, Aber J D, Howarth R W, Likens G E, Matson P A, Schindler D W, Schlesinger W H, Tilman D G. Human alteration of the global nitrogen cycle: sources and consequences. Ecological Applications, 1997, 7(3): 737- 750.

[3] Galloway J N, Townsend A R, Erisman J W, Bekunda M, Cai Z C, Freney J R, Martinelli L A, Seitzinger S P, Sutton M A. Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science, 2008, 320(5878): 889- 892.

[4] Gu B J, Chang J, Min Y, Ge Y, Zhu Q A, Galloway J N, Peng C H. The role of industrial nitrogen in the global nitrogen biogeochemical cycle. Scientific Reports, 2013, 3: 2579.

[5] 陳秋鳳. 杉木人工林林木養分和凋落物分解對模擬氮沉降的響應[D]. 福州: 福建農林大學, 2006.

[6] Galloway J N, Dentener F J, Capone D G, Boyer E W, Howarth R W, Seitzinger S P, Asner G P, Cleveland C C, Green P A, Holland E A, Karl D M, Michaels A F, Porter J H, Townsend A R, V?osmarty C J. Nitrogen cycles: past, present, and future. Biogeochemistry, 2004, 70(2): 153- 226.

[7] Lu C Q, Tian H Q, Liu M L, Ren W, Xu X F, Chen G S, Zhang C. Effect of nitrogen deposition on China′s terrestrial carbon uptake in the context of multifactor environmental changes. Ecological Applications, 2012, 22(1): 53- 75.

[8] Post W M, Emanuel W R, Zinke P J, Stangenberger A G. Soil carbon Pools and world life zones. Nature, 1982, 298(5870): 156- 159.

[9] 王紹強, 劉紀遠. 土壤碳蓄積量變化的影響因素研究現狀. 地球科學進展, 2002, 17(4): 528- 534.

[10] Cao M K, Woodward F I. Dynamic responses of terrestrial ecosystem carbon cycling to global climate change. Nature, 1998, 393(6682): 249- 252.

[11] 彭琴, 董云社, 齊玉春. 氮輸入對陸地生態系統碳循環關鍵過程的影響. 地球科學進展, 2008, 23(8): 874- 883.

[12] 鄧小文, 韓士杰. 氮沉降對森林生態系統土壤碳庫的影響. 生態學雜志, 2007, 26(10): 1622- 1627.

[13] Wei X H, Blanco J A, Jiang H, Kimmins J P H. Effects of nitrogen deposition on carbon sequestration in Chinese fir forest ecosystems. Science of the Total Environment, 2012, 416: 351- 361.

[14] 呂超群. 中國大氣氮沉降與陸地生態系統碳循環研究[D]. 北京: 中國科學院地理科學與資源研究所, 2009.

[15] Jain A, Yang X J, Kheshgi H, McGuire A D, Post W, Kicklighter D. Nitrogen attenuation of terrestrial carbon cycle response to global environmental factors. Global Biogeochemical Cycles, 2009, 23(4): 1- 13.

[16] de Vries W, Solberg S, Dobbertin M, Sterba H, Laubhann D, van Oijen M, Evans C, Gundersen P, Kros J, Wamelink G W W, Reinds G J, Sutton M A. The impact of nitrogen deposition on carbon sequestration by European forests and heathlands. Forest Ecology and Management, 2009, 258(8): 1814- 1823.

[17] 于磊, 張柏. 中國黑土退化現狀與防治對策. 干旱區資源與環境, 2004, 18(1): 99- 103.

[18] 朱鋒, 劉志明, 王宗明, 宋開山. 東北地區農田凈初級生產力時空特征及其影響因素分析. 資源科學, 2010, 32(11): 2079- 2084.

[19] 于穎, 范文義, 李明澤. 東北林區不同尺度森林的含碳率. 應用生態學報, 2012, 23(2): 341- 346.

[20] 奚小環, 楊忠芳, 崔玉軍, 孫淑梅, 于成廣, 李敏. 東北平原土壤有機碳分布與變化趨勢研究. 地學前緣(中國地質大學(北京); 北京大學), 2010, 17(3): 213- 221.

[21] 王麗麗, 宋長春, 葛瑞娟, 宋艷宇, 劉德燕. 三江平原濕地不同土地利用方式下土壤有機碳儲量研究. 中國環境科學, 2009, 29(6): 656- 660.

[22] 張春華, 王宗明, 任春穎, 宋開山, 張柏, 劉殿偉. 松嫩平原玉米帶農田土壤有機碳時空格局. 農業工程學報, 2010, 26(S1): 300- 307.

[23] 崔玉軍, 李延生, 劉立芬, 胡春見. 松嫩平原南部土壤近二十年有機碳含量變化特征. 黑龍江國土資源, 2010, (8): 41- 41.

[24] 王金達, 劉景雙, 劉淑霞, 于君寶. 松嫩平原黑土土壤有機碳庫的估算及其影響因素. 農業環境科學學報, 2004, 23(4): 687- 690.

[25] 方華軍, 楊學明, 張曉平. 東北黑土有機碳儲量及其對大氣CO2的貢獻. 水土保持學報, 2003, 17(3): 9- 12.

[26] 徐金忠, 孟凱, 崔曉陽, 馮軍, 宋金鳳. 不同施肥處理對黑土腐殖質組成的影響. 東北林業大學學報, 2009, 39(11): 84- 85.

[27] 徐寧, 周連仁, 苗淑杰. 長期施肥對黑土有機質及其組成的影響. 中國土壤與肥料, 2012, (6): 14- 16, 33- 33.

[28] 毛德華, 王宗明, 韓佶興, 任春穎. 1982—2010年中國東北地區植被NPP時空格局及驅動因子分析. 地理科學, 2012, 32(9): 1106- 1111.

[29] 王春春, 陳長青, 黃山, 鄧艾興, 張衛建, 卞新民. 東北氣候和土壤資源演變特征研究. 南京農業大學學報, 2010, 33(2): 19- 24.

[30] Wang S Q, Zhou C H, Liu J Y, Tian H Q, Li K R, Yang X M. Carbon storage in northeast China as estimated from vegetation and soil inventories. Environmental Pollution, 2002, 116(S1): S157-S165.

[31] 欒江, 仇煥廣, 井月, 廖紹攀, 韓煒. 我國化肥施用量持續增長的原因分解及趨勢預測. 自然資源學報, 2013, 28(11): 1869- 1878.

[32] 張四代, 王激清, 張衛峰, 王利, 馬文奇, 張福鎖. 我國東北地區化肥消費與生產現狀、問題及其調控策略. 磷肥與復肥, 2007, 22(5): 74- 78.

[33] Hutchinson M F. A new objective method for spatial interpolation of meteorological variables from irregular networks applied to the estimation of monthly mean solar radiation, temperature, precipitation and windrun // Fitzpatrick E A, Kalma J D, eds. Need for Climatic and Hydrological Data in Agriculture in South-East Asia. Proceedings of the United Nations University Workshop, December 1983, Division of Water Resources Technical Memorandum 89/5, CSIRO, Canberra, 1989: 95- 104.

[34] 彭暢, 朱平, 高洪軍, 劉淑環. 長期定位監測黑土土壤肥力的研究Ⅰ. 黑土耕層有機質與氮素轉化. 吉林農業科學, 2004, 29(5): 29- 33.

[35] 隋躍宇, 張興義, 焦曉光, 王其存, 趙軍. 長期不同施肥制度對農田黑土有機質和氮素的影響. 水土保持學報, 2005, 19(6): 190- 192.

[36] 申建波, 張福鎖. 水稻養分資源綜合管理理論與實踐[M]. 北京: 中國農業大學出版社, 2006: 310- 311.

[37] 張衛峰, 季玥秀, 馬驥, 王雁峰, 馬文奇, 張福鎖. 中國化肥消費需求影響因素及走勢分析——Ⅰ化肥供應. 資源科學, 2007, 29(6): 162- 169.

[38] 楊青林, 桑利民, 孫吉茹, 季志強, 袁文利, 郭玉煒, 蓋顏欣. 我國肥料利用現狀及提高化肥利用率的方法. 山西農業科學, 2011, 39(7): 690- 692.

[39] Erickson H, Keller M, Davidson E A. Nitrogen oxide fluxes and nitrogen cycling during postagricultural succession and forest fertilization in the humid tropics. Ecosystems, 2001, 4(1): 67- 84.

[40] Mirmanto E, Proctor J, Green J, Nagy L, Suriantata. Effects of nitrogen and phosphorus fertilization in a lowland evergreen rainforest. Philosophical Transactions of the Royal Society B: Biological Sciences, 1999, 354(1391): 1825- 1829.

[41] Pregitzer K S, Burton A J, Zak D R, Talhelm A F. Simulated chronic nitrogen deposition increases carbon storage in Northern Temperate forests. Global Change Biology, 2008, 14(1): 142- 153.

[42] Magill A H, Aber J D, Hendricks J J, Bowden R D, Melillo J M, Steudler P A. Biogeochemical response of forest ecosystems to simulated chronic nitrogen deposition. Ecological Applications, 1997, 7(2): 402- 415.

[43] Eriksson H M, Berdén M, Rosén K, Nilsson S I. Nutrient distribution in a Norway spruce stand after long-term application of ammonium nitrate and superphosphate. Water, Air, and Soil Pollution, 1996, 92(3/4): 451- 467.

[44] 樊后保, 劉文飛, 裘秀群, 徐雷, 王強, 陳秋鳳. 杉木人工林凋落物量對氮沉降增加的初期響應. 生態學雜志, 2007, 26(9): 1335- 1338.

[45] Deforest J L, Zak D R, Pregitzer K S, Burton A J. Atmospheric nitrate deposition and the microbial degradation of cellobiose and vanillin in a northern hardwood forest. Soil Biology and Biochemistry, 2004, 36(6): 965- 971.

[46] 涂利華, 胡庭興, 張健, 何遠洋, 田祥宇, 肖銀龍. 模擬氮沉降對華西雨屏區苦竹林細根特性和土壤呼吸的影響. 應用生態學報, 2010, 21(10): 2472- 2478.

[47] 李建敏, 丁維新, 蔡祖聰. 氮肥對玉米生長季土壤呼吸的影響. 應用生態學報, 2010, 21(8): 2025- 2030.

[48] 李凱, 江洪, 由美娜, 曾波. 模擬氮沉降對石櫟和苦櫧幼苗土壤呼吸的影響. 生態學報, 2011, 31(1): 82- 89.

[49] Berg B, Matzner E. Effect of N deposition on decomposition of plant litter and soil organic matter in forest systems. Environmental Reviews, 1997, 5(1): 1- 25.

[50] Franklin O, H?gberg P, Ekblad A, ?gren G I. Pine forest floor carbon accumulation in response to N and PK additions: bomb14C modelling and respiration studies. Ecosystems, 2003, 6(7): 644- 658.

[51] 賈淑霞, 王政權, 梅莉, 孫玥, 全先奎, 史建偉, 于水強, 孫海龍, 谷加存. 施肥對落葉松和水曲柳人工林土壤呼吸的影響. 植物生態學報, 2007, 31(3): 372- 379.

[52] Micks P, Aber J D, Boone R D, Davidson E A. Short-term soil respiration and nitrogen immobilization response to nitrogen applications in control and nitrogen-enriched temperate forests. Forest Ecology and Management, 2004, 196(1): 57- 70.

[53] Batjes N H, Bridges E M. A Review of Soil Factors and Processes that Control Fluxes of Heat, Moisture and Greenhouse Gases[M]. Wageningen, Netherlands: International Soil Reference and Information Centre, 1992.

[54] Enríquez S, Duarte C M, Sand-Jensen K. Patterns in decomposition rates among photosynthetic organisms: the importance of detritus C:N:P content. Oecologia, 1993, 94(4): 457- 471.

[55] Gholz H L, Wedin D A, Smitherman S M, Harmon M E, Parton W J. Long-term dynamics of pine and hardwood litter in contrasting environments: toward a global model of decomposition. Global Change Biology, 2000, 6(7): 751- 765.

Modeling the effect of nitrogen input on soil carbon storage in Northeast China

GU Fengxue1, HUANG Mei2, ZHANG Yuandong3,*, LI Jie1, YAN Huimin2, GUO Rui1, ZHONG Xiuli1

1 Key Laboratory of Dryland Agriculture, Ministry of Agriculture, Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences, Beijing 100081, China 2KeyLaboratoryofEcosystemNetworkObservationandModeling,InstituteofGeographicSciencesandNaturalResourcesResearch,ChineseAcademyofSciences,Beijing100101,China3KeyLaboratoryofForestEcologyandEnvironment,StateForestryAdministration,InstituteofForestEcology,EnvironmentandProtection,ChineseAcademyofForestry,Beijing100091,China

Anthropogenic activities have altered the global nitrogen (N) cycle, leading to increased N input into the ecosystem through N deposition and N fertilizer. China consumes the highest amount of fertilizer globally, and it has become the third largest N deposition region in the world. N input is an important factor in the terrestrial carbon (C) sink, and N input is implicated in a series of biogeochemical cycles that then influence the C cycle and its spatial pattern. The soil C pool plays an important role in the global C cycle; therefore, the question of whether and to what extent N input affects the soil C pool must be addressed in relation to global change and N deposition. Northeast China has the third-largest black belt in the world and it is an important commodity grain base of China. Cropland is the most prominent manner of land use in the Northeast, and the change in its soil C pool is an essential factor in explaining the regional C budget and maintaining soil fertility. Using a process-based carbon-water-nitrogen coupling model, CEVSA2, we simulated the spatial patterns of N deposition in China between 1961 and 2010. With the N deposition and environment data collected, we also used CEVSA2 to explore the effect of enhanced N input on soil C storage in Northeast China in this period. The results show that the N deposition in Northeast China (1.00 gN m-2a-1) is higher than the average for the whole country, and increases annually by 0.047 gN m-2a-1. From the 1980, the N input began to increase significantly. The model simulations indicated that enhanced N deposition has resulted in a net increase of soil C density by 135 gC/m2in Northeast, and the total soil C storage increased by 0.16 PgC over the past 50 years. The soil C sequestration induced by N input decreased from east to west and from south to north. In the Songliao Plain and the Sanjiang Plain, the soil C sequestration induced by N input was >300 gC/m2. The change in soil C density induced by N input varies significantly among the major biome types in the Northeast. Model simulations indicated that the increase in soil C density in cropland was 230 gC/m2. However, those in the forest, shrubland, and grassland were 76 gC/m2, 169 gC/m2and 89 gC/m2, respectively, and the spatial heterogeneity of N input as well as the different responses of vegetation types to it determined the spatial pattern of the increase in soil C storage. Elucidation of the effect of N input on soil C storage will provide a scientific foundation for C fixation as well as C and N management in cropland.

nitrogen deposition; fertilization; soil carbon storage; Northeast; CEVSA2 model

國家自然科學基金項目(41271118, 31370463, 31070398)；中國農業科學院科技創新工程項目;國家留學基金資助(留金發[2013]3018號)

2015- 02- 26;

日期:2015- 12- 14

10.5846/stxb201502260389

*通訊作者Corresponding author.E-mail: zyd@caf.ac.cn

顧峰雪，黃玫，張遠東，李潔，閆慧敏，郭瑞，鐘秀麗.氮輸入對中國東北地區土壤碳蓄積的影響.生態學報,2016,36(17):5379- 5390.

Gu F X, Huang M, Zhang Y D, Li J, Yan H M, Guo R, Zhong X L.Modeling the effect of nitrogen input on soil carbon storage in Northeast China.Acta Ecologica Sinica,2016,36(17):5379- 5390.