模擬氮沉降對油松土壤水溶性碳氮及其光譜特征的影響

翟珈瑩, 何俐蓉, 吳 陽, 王國梁, 薛 萐

(1.中國科學院 水利部 水土保持研究所 黃土高原土壤侵蝕與旱地農業國家重點實驗室,陜西 楊凌 712100; 2.中國科學院大學, 北京 100049; 3.陜西地建土地工程技術研究院, 西安 710016;4.西北農林科技大學 林學院, 陜西 楊凌 712100; 5.西北農林科技大學 水土保持研究所, 陜西 楊凌 712100)

黃土丘陵區由于其特殊的母質、氣候和地形特征，加上長期以來不合理的土地利用，水土流失嚴重，生態環境十分脆弱。當全球變化因子——氮沉降在該生態系統時，會對系統產生一系列影響，而其中的養分分配，尤其是對環境變化響應敏感的水溶性碳、氮組分對于氮沉降的響應目前還沒有一致的結論。因此，本文選擇黃土高原植被恢復中分布最廣泛的主要樹種之一——油松(PinustabuliformisCarr.)作為研究對象，通過氮添加的方式模擬氮沉降，分析土壤水溶性碳氮含量的變化特征，為評價未來全球變化背景下的土壤生態系統響應提供數據支持。

1 材料和方法

1.1 試驗區概況

試驗在中國陜西省北部，延安市宜川縣鐵龍灣林場(110°06′E，35°39′N)。試驗區屬大陸性氣候，年平均氣溫9.8℃，年平均降水量為584.4 mm，海拔860～1 200 m，坡度20°～25°，土壤類型為灰褐色森林土。地帶性植被為溫帶落葉闊葉林，遼東櫟(QuercuswutaishanseaMary)為頂極群落。主要喬木林為人工油松林(PinustabuliformisCarrière)和黃楊林(Buxussinica(Rehder & E.H.Wilson) M.Cheng)。灌叢中有興安胡枝子(Lespedezadaurica(Laxm.) Schindl.)、胡頹子(ElaeagnuspungensThunb.)、繡線菊(SpiraeasalicifoliaL.)和中間錦雞兒(CaraganaintermediaKuang et H.C.Fu)。草本植物包括異葉敗醬(PatriniaheterophyllaBunge)、茜草(RubiacordifoliaL.)、鐵桿蒿(Artemisiagmelinii)、野棉花(AnemonevitifoliaBuch-Ham.)等。試驗地選在氣候、位置、林形基本一致的人工油松林地，該人工油松林建于1966年，面積600 hm2。

1.2 試驗設計及樣品采集

在油松人工林地中設置16個小區，圍封面積為100 m2。參考全球氮沉降水平和國內外氮添加試驗的含量[16]，本研究設置4個處理CK,N3,N6,N9，依次為0,3,6,9 g/(m2·a)尿素的形式進行添加，以上每個處理設4個重復。自2014年4月起，每年4月、6月、8月、10月采用噴霧器對林下土壤分4次等量直接噴施，于2015年9月進行土壤樣品采集，方法為從每個地塊隨機選取6個樣點，取樣深度為0—10 cm和10—20 cm，相同層次樣品混勻，剔除根系后過2 mm的篩孔，用于水溶性有機碳、氮的測定。

1.3 測定方法

土壤水溶性有機碳、氮的提取及測定參考Jones等[17]的方法進行。

1.3.1 可溶性物質提取方法 稱取40.00 g鮮土(同時測定土壤含水量)，加蒸餾水120 ml(水土比為3∶1)浸提，在25℃下震蕩機振蕩0.5 h后，置高速離心機中(8 000 r/min)離心10 min，收集過濾液用0.45 μm濾膜進行抽濾，得到的浸提液置于4℃冰箱中保存備用。

1.3.3 三維熒光光譜測定 DOM的熒光光譜采用F-4600日立熒光光譜儀進行測定。PMT電壓設置為700 V，激發波長掃描范圍為200～450 nm，發射波長掃描范圍為250～450 nm，掃描速度為1 200 nm/min，狹縫寬度為5 nm。在分析之前，需要對水樣進行處理去除拉曼散射效應和瑞利散射效應。本研究中，一共有48個樣本，每個樣本包括了包含了51激發×41發射的矩陣。因此，總樣本構成了一個16×52×49的矩陣，通過PARAFAC法進行分析[18]。將三維熒光光譜分為2個熒光區域，識別2個熒光區域的熒光峰及對應的激發/發射光波長[19]，進行相應的區域積分獲得區域積分值[20]，計算出f450/f500指數，用以表征土壤中水溶性有機物的來源[21]。

1.3.4 土壤理化性質測定 土壤有機質(DOM)采用重鎘酸鉀外加熱法測定[22]。

1.4 統計分析

數據的統計采用Microsoft Office Excel完成，利用SPSS 22.0進行單因素方差分析(One-Way ANOVA)和LSD多重比較(顯著水平為p<0.05)。采用MATLAB的DOM Fluor工具盒[23]對熒光EEMs進行平行因子分析(PARAFAC)，并采用MATLAB R2016年(美國數學工程公司)對EEMs進行繪制。采用Origin 9.0進行制圖。

2 結果與分析

2.1 氮添加對不同土層水溶性氮及組分的影響

注：不同字母表示氮添加處理在各土層存在顯著差異(p<0.05)，下同。

氮添加改變了土壤水溶性氮組分占WSN的比例(圖2)，不同土層WSON/WSN比值隨氮添加水平先降低后升高，在N3顯著降低，隨后升高。0—10 cm土層WSIN/WSN含量隨氮添加水平先升高后降低，在N9達到了顯著水平，10—20 cm土層隨氮添加水平的增加而增加，在N3達到了顯著水平。

圖2 不同處理土壤的WON/WSN和WSIN/WSN

2.2 氮添加對不同土層水溶性碳的影響

氮添加改變了土壤水溶性碳及其占有機碳的比例(圖3)，0—10 cm土層WSOC隨氮添加水平先升高后降低，在N9顯著降低，10—20 cm土層各處理無顯著差異。0—10 cm土層WSOC/SOC隨氮添加水平的升高呈降低的趨勢，10—20 cm土層各處理無顯著差異。

圖3 不同處理土壤的WSOC含量以及WSOC/SO

2.3 氮添加對水溶性有機質三維熒光光譜特性的影響

樣品的三維熒光光譜經PARAFAC分析，得到水溶性有機質主要是由兩類物質組成。兩種組分的三維熒光光譜如圖4所示，組分1的Ex/Em波長以210～250/400～450 nm為中心，該組分被識別為類富里酸熒光基團，屬腐殖質物質。組分2中，有一個類似于峰狀的熒光峰，其Ex/Em集中于220～225/275～300 nm，該成分為類芳香族蛋白質熒光基團。分析得到不同土層中組分1和組分2在樣本中的含量(圖5—6)，不同土層類腐殖質隨著氮添加水平先降低后升高，組分2含量在各土層均隨氮添加水平的增加呈增長的趨勢。

圖4 組分1，2的三維熒光光譜

3 討論與結論

圖5 不同處理下組分1，2的含量

圖6 不同土層在不同處理下具有代表性的熒光發射激發基質(EEMs)

經對土壤中水溶性有機質(DOM)的三維熒光光譜PARAFAC分析后得到兩種存在土壤中的有機物，一種是類富里酸物質，一種是芳香族蛋白質物質。類富里酸物質在不同土層隨氮添加水平先降后升，這可能是因為氮添加促進了微生物對腐殖質的分解，高氮添使微生物降解中碳的限制加劇[4]，導致土壤DOM中類腐殖質物質有一個先降后升的趨勢，WSOC的含量也隨氮輸入先增后減。0—10 cm土層的這種物質含量高于10—20 cm，可能是因為隨著土層深度的增加，微生物數量減少，其活性下降，由微生物主導的土壤礦化速率降低[32]，10—20 cm土層腐殖化程度低。而芳香族蛋白質物質應該是微生物降解后的某種類蛋白質物質，在土壤表層呈增長趨勢，森林土壤中腐殖質的淋溶過程及微生物分解是DOM的主要形成途徑[33]，因淋溶作用10—20 cm的含量有所增加。因此在黃土丘陵區，氮輸入通過影響根際微生物的活動影響土壤的養分分配，從而影響水溶性有機碳、氮組分的含量變化。