青海省祁連圓柏天然林林地土壤生態(tài)化學(xué)計(jì)量特征

井 賽，陳志林，侯曉巍，李智華，胡云云，向安民，侯 琳*

(1.西北農(nóng)林科技大學(xué) 林學(xué)院，陜西 楊陵 712100；2.國家林業(yè)和草原局西北調(diào)查規(guī)劃設(shè)計(jì)院 旱區(qū)生態(tài)水文與災(zāi)害防治國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710048)

生態(tài)化學(xué)計(jì)量學(xué)主要關(guān)注生態(tài)過程中多重化學(xué)元素(C、N、P)平衡，表征C、N、P的耦合關(guān)系，已被廣泛應(yīng)用于生態(tài)系統(tǒng)各層次的種群動(dòng)態(tài)、群落生產(chǎn)力限制因素、養(yǎng)分循環(huán)與分布等研究中[1-4]。土壤作為生態(tài)系統(tǒng)的重要組成部分，提供了植物生長所必須的養(yǎng)分元素，直接影響植被群落的結(jié)構(gòu)、林地生產(chǎn)力水平與生態(tài)系統(tǒng)的穩(wěn)定性[5-6]。探究土壤生態(tài)化學(xué)計(jì)量特征已逐步成為揭示林地養(yǎng)分狀況與營養(yǎng)循環(huán)，實(shí)現(xiàn)生態(tài)系統(tǒng)可持續(xù)經(jīng)營的新興工具[7-8]。

森林土壤C、N、P化學(xué)計(jì)量沿環(huán)境梯度變化的地理分布格局已成為共識[9-12]。海拔作為關(guān)鍵的地理因子，其綜合了氣候、土壤、植被因素的變化，圍繞海拔梯度森林土壤C、N、P生態(tài)化學(xué)計(jì)量研究已成為熱點(diǎn)[9-10]。有研究發(fā)現(xiàn)，土壤有機(jī)碳(SOC)、全N(TN)、C∶P、N∶P隨海拔升高呈現(xiàn)上升趨勢[8]，但其他研究卻發(fā)現(xiàn)相反的變化趨勢[5,11]。以往對青海林地土壤研究表明，SOC隨海拔呈單峰變化，TN隨海拔呈增加趨勢[12]。還有研究發(fā)現(xiàn)，土壤生態(tài)化學(xué)計(jì)量特征與海拔梯度間存在非線性關(guān)系[5]。已有的研究結(jié)果并不一致且多集中于不同植被類型間土壤生態(tài)化學(xué)計(jì)量特征的比較[5,8,13]，對祁連圓柏天然林林地土壤生態(tài)化學(xué)計(jì)量特征關(guān)注較少。

祁連圓柏(Juniperusprzewalskii)為柏科圓柏屬常綠喬木，因耐寒、耐旱、耐瘠薄成為青海省高寒生境中分布最廣的優(yōu)勢成林樹種之一，在水源涵養(yǎng)、防風(fēng)固沙、應(yīng)對氣候變化、維持生物多樣性等方面發(fā)揮著不可替代作用[14]。同時(shí)，祁連圓柏由于其廣泛的地理分布，生境水熱條件變化顯著，表現(xiàn)出明顯的垂直分異性[15]。目前，關(guān)于祁連圓柏的研究多聚焦于樹輪生長、造林、生理生態(tài)特性、葉片性狀變異[16-20]。本研究以青海省祁連圓柏7個(gè)典型分布區(qū)的天然林林地土壤為研究對象，探討祁連圓柏天然林土壤生態(tài)化學(xué)計(jì)量的垂直分布特征與影響因素，以期為林地土壤管理提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

選取祁連圓柏天然林分布的典型區(qū)域德令哈市、都蘭縣、互助縣、澤庫縣、祁連縣、烏蘭縣、興?？h7個(gè)地點(diǎn)為研究區(qū)(表1)。區(qū)域范圍97.33°-102.54° E，35.01°-38.22°N，海拔2 340～3 870 m，年平均氣溫-1.33～5.28℃，年降雨量51.60～502.50 mm，降水集中于5-9月，林下土壤類型多為山地褐色針葉林土和山地灰褐色森林土。7個(gè)研究區(qū)的年平均溫度和年平均降水等氣象數(shù)據(jù)從中國氣象數(shù)據(jù)網(wǎng)站(http://data.cma.cn/)獲取。

表1 研究區(qū)基本概況

1.2 研究方法

1.2.1 樣地設(shè)置與樣品采集 依據(jù)典型性原則，結(jié)合7個(gè)縣祁連圓柏林的海拔分布范圍(2 340～3 870 m)，以200 m為1個(gè)梯度進(jìn)行海拔劃分，即<2 900、2 900～3 100、3 100～3 300、3 300～3 500、3 500～3 700、>3 700 m共6個(gè)海拔梯度。在7個(gè)縣布設(shè)20 m×20 m的群落樣地35個(gè)。根據(jù)對角線法，在每個(gè)群落樣地內(nèi)設(shè)置5個(gè)5 m×5 m灌木和5個(gè)1 m×1 m草本調(diào)查樣方。利用GPS記錄每個(gè)樣地的經(jīng)度、緯度、海拔，對樣地內(nèi)的喬木進(jìn)行計(jì)數(shù)，測定胸徑與樹高；在各樣地的中心與四角位置設(shè)置5個(gè)5 m×5 m灌木和5個(gè)1 m×1 m草本調(diào)查樣方；灌木樣方中，記錄灌木的種類、株數(shù)，測定高度、蓋度、基徑，記錄草本植物種類、株數(shù)，測定高度、蓋度。

每個(gè)群落樣地內(nèi)按照“S”型取樣法選取土壤樣點(diǎn)5個(gè)，挖掘土壤剖面，因調(diào)查區(qū)土層厚度各異，為確保數(shù)據(jù)的可比性，在各土壤樣點(diǎn)按0～20、20～40、40～60 cm 3個(gè)土層用體積為100 cm3的環(huán)刀采集各層原狀土，同時(shí)采集同深度散土。同一樣地內(nèi)相同土層的散土混合均勻成一個(gè)土樣。將采集的環(huán)刀土和散土進(jìn)行現(xiàn)場稱重、編號并裝入土樣袋中帶回實(shí)驗(yàn)室。

1.2.2 樣品分析 將采集的土壤樣品風(fēng)干，去雜(撿去其中的植物殘?bào)w、礫石與侵入體)，磨細(xì)，過篩，用于測定土壤pH與SOC、TN、TP。

土壤有機(jī)碳(SOC)采用重鉻酸鉀-外加熱法測定；土壤全N(TN)采用凱氏定氮法測定；土壤全P(TP)采用氫氧化鈉熔融-鉬銻抗比色法測定；土壤pH利用pH計(jì)測定(水土比為2.5∶1)；土壤含水率采用經(jīng)典的烘干法測定；土壤容重利用環(huán)刀法測定[21]。

1.2.3 物種多樣性計(jì)算 本研究選取物種豐富度指數(shù)(S)、Shannon-wiener指數(shù)(H)來衡量群落多樣性特征，計(jì)算方法分別為：

物種豐富度指數(shù)(S)：S=樣方內(nèi)的物種數(shù)目

Shannon-wiener指數(shù)(H)：

H=-∑PilnPi

式中，S為樣方內(nèi)物種數(shù)目，Pi為物種i的相對重要值，即Pi=Ni/N，Ni為物種i的重要值，N為所有物種的重要值之和，重要值=(相對頻度+相對蓋度+相對高度)/3[22]。

1.3 數(shù)據(jù)處理

運(yùn)用Excel 2016和SPSS 22.0對數(shù)據(jù)進(jìn)行整理與分析。依據(jù)拉伊達(dá)準(zhǔn)則(3σ 準(zhǔn)則)排除異常值。采取單因素方差分析(ANOVA)檢驗(yàn)不同海拔與土層下土壤C、N、P化學(xué)計(jì)量特征的差異，若方差為齊性，用最小顯著差異法(LSD)進(jìn)行多重比較；若方差為非齊性，用Tamhane’s T2法進(jìn)行多重比較，所有檢驗(yàn)的顯著性水平為P<0.05或P<0.01。同時(shí)對海拔與土壤化學(xué)計(jì)量間進(jìn)行線性和二次回歸方程擬合，根據(jù)模型解釋度(R2)確定最優(yōu)回歸模型。運(yùn)用Pearson檢驗(yàn)衡量土壤C、N、P及其化學(xué)計(jì)量比間的相關(guān)性，利用RDA與Pearson相關(guān)性分析環(huán)境因素對土壤生態(tài)化學(xué)計(jì)量特征的影響。作圖采用Origin 2018。

2 結(jié)果與分析

2.1 祁連圓柏林土壤生態(tài)化學(xué)計(jì)量特征

土壤SOC、TN、TP含量平均值分別為36.49、3.43、0.62 g·kg-1，變化范圍為9.05～75.92、0.61～7.47、0.34～0.91 g·kg-1，變異系數(shù)為16.13%～45.77%(表3)。SOC與TN的變異系數(shù)相近，對外界環(huán)境變化響應(yīng)具有一致性，TP的變異較SOC和TN??；土壤化學(xué)計(jì)量比中，C∶N、C∶P、N∶P的變化范圍分別為7.13～17.17、15.60～105.44和1.24～9.64，變異系數(shù)在1.91%～37.02%，C∶P與N∶P變異系數(shù)相近(表3)。

表2 土壤C、N、P及其生態(tài)化學(xué)計(jì)量統(tǒng)計(jì)特征

2.2 祁連圓柏天然林土壤生態(tài)化學(xué)計(jì)量特征沿海拔梯度的變化

在同一土層不同海拔梯度下，土壤SOC隨海拔的增加呈降低趨勢，在<3 500 m時(shí)顯著高于其他海拔(3 500～3 700、>3 700 m)(P<0.05)(圖1)。土壤全N與有機(jī)碳變化規(guī)律相似，3個(gè)土層最高值均出現(xiàn)在<2 900 m，最低值出現(xiàn)在>3 700 m，隨海拔的增加而降低。在<2 900～3 700 m海拔范圍內(nèi)，土壤全P隨海拔上升小幅度下降，至海拔3 700 m略微上升，且在>3 700 m顯著低于其他5個(gè)海拔梯度(P<0.05)。土壤C∶N隨海拔梯度增加無明顯變化，且各海拔間無顯著差異(P>0.05)。C∶P與N∶P最大值出現(xiàn)在<2 900 m，最小值在最高海拔梯度(>3 700 m)，多重比較顯示，除海拔>3 700 m顯著低于其他海拔梯度外其余4個(gè)海拔梯度間差異并不顯著(P>0.05)(圖1)。

在同一海拔梯度不同土層內(nèi)，土壤SOC、TN含量均隨土壤深度的增加而逐漸降低，且0～20 cm土層的SOC、TN含量顯著高于20～40 cm與40～60 cm土層(P<0.05)(圖1)，表現(xiàn)為明顯的“表聚現(xiàn)象”。土壤TP雖隨土層增加表現(xiàn)為減小趨勢，但減小程度小于SOC與TN，同時(shí)0～20 cm土層磷含量顯著高于40～60 cm(P<0.05)(圖1)。在相同海拔梯度內(nèi)，土壤化學(xué)計(jì)量比隨土層加深呈現(xiàn)出各異的變化規(guī)律。其中，C∶P、N∶P表現(xiàn)為相同的變化趨勢，隨土層的加深而降低，而C∶N隨土層的加深無明顯變化，在各土層間的差異不顯著(P>0.05)(圖1)。

注：不同大寫字母表示同一土層不同海拔梯度下的差異；不同小寫字母表示同一海拔梯度不同土層間的差異。

進(jìn)一步對土壤SOC、TN、TP及其化學(xué)計(jì)量與各樣點(diǎn)的實(shí)際海拔進(jìn)行回歸分析，結(jié)果表明(圖2)，海拔對土壤3個(gè)土層SOC、TN、C∶P、N∶P具有顯著影響(P<0.05)，對土壤C∶N無顯著影響(P>0.05)。其中，土壤SOC、TN、C∶P、N∶P與海拔呈呈現(xiàn)二次曲線負(fù)相關(guān)關(guān)系，土壤TP與海拔間呈線性負(fù)相關(guān)關(guān)系，土壤C∶N與海拔間的回歸關(guān)系均未達(dá)到顯著水平(圖2)。

圖2 土壤化學(xué)計(jì)量與海拔間的回歸分析

相關(guān)分析表明，SOC、TN、TP含量兩兩之間均呈現(xiàn)極顯著正相關(guān)(P<0.01)，且SOC與TN之間相關(guān)系數(shù)最大(表3)。SOC與TP作為C∶P的分子和分母，理論上應(yīng)與C∶P分別呈正相關(guān)、負(fù)相關(guān)，而本研究結(jié)果SOC、TP與C∶P均呈極顯著正相關(guān)(P<0.01)(表3)，因此C∶P的變化主要由SOC引起。N∶P與TN的相關(guān)性強(qiáng)于TP，說明N∶P的變化主要有TN引起。

表3 土壤C、N、P及其化學(xué)計(jì)量的相關(guān)性

2.3 氣候因素沿海拔梯度的變化

對各研究區(qū)的年平均溫度與年平均降水與海拔梯度間進(jìn)行相關(guān)分析，結(jié)果表明(表4)，隨海拔的升高，年平均溫度與年平均降水呈降低的變化趨勢，且相關(guān)性均達(dá)到顯著水平(P<0.05)。

表4 研究區(qū)海拔與氣候因素的關(guān)系

2.4 群落特征與土壤指標(biāo)沿海拔梯度的變化

從青海祁連圓柏天然林群落多樣性指數(shù)可以看出，隨海拔的升高，總體灌木、草本的Shannon-Winner指數(shù)均呈減小趨勢，且在不同海拔梯度下差異性顯著(P<0.05)(表5)。然而從物種豐富度來看，灌木與草本物種豐富度隨海拔的升高無規(guī)律性變化特征，僅群落總體的物種豐富度隨海拔的升高而顯著降低(P<0.05)。林分密度隨海拔梯度增加雖呈減小趨勢，但在各海拔梯度間無顯著差異(P>0.05)(表5)。不同海拔梯度下土壤含水率、pH、容重差異顯著(P<0.05)，其中土壤含水率隨海拔梯度的升高呈減小趨勢，且在高海拔區(qū)域(3 500～3 700、>3 700 m)區(qū)域顯著低于低海拔區(qū)域(<3 500 m)；相反，土壤pH、容重最大值出現(xiàn)在海拔>3 700 m，最小值出現(xiàn)在海拔<2 900 m，表明土壤pH、容重隨海拔梯度的升高呈增加趨勢(表5)。

表5 不同海拔土壤理化性質(zhì)與群落特征

2.5 環(huán)境因素對土壤生態(tài)化學(xué)計(jì)量特征的影響

以土壤SOC、TN、TP及其化學(xué)計(jì)量比為響應(yīng)變量，氣候、土壤指標(biāo)、與群落特征為解釋變量，進(jìn)行冗余分析。結(jié)果表明，前2個(gè)排序軸累積解釋土壤SOC、TN、TP及其化學(xué)計(jì)量比變異的75.93%，其中第1排序軸解釋變異的75.28%，第2排序軸解釋變異的0.65%(圖2)。土壤SOC、TN與群落總體Shannon-Winner指數(shù)(H)、物種豐富度(R)、林分密度、降水、土壤含水率呈極顯著正相關(guān)(P<0.01)，與溫度呈顯著正相關(guān)(P<0.05)，與海拔、pH、容重呈極顯著負(fù)相關(guān)(P<0.01)。土壤TP與群落總體Shannon-Winner指數(shù)(H)、物種豐富度、林分密度、土壤含水率呈極顯著正相關(guān)，與海拔、pH、容重呈極顯著負(fù)相關(guān)(P<0.01)。土壤C∶N與各影響因素間相關(guān)性均未達(dá)到顯著水平(P>0.05)，土壤C∶P、N∶P與群落總體Shannon-Winner指數(shù)(H)、林分密度、降水、土壤含水率呈極顯著正相關(guān)，與海拔、pH、容重呈極顯著負(fù)相關(guān)(P<0.01)(圖3、表6)。冗余分析單獨(dú)效應(yīng)表明土壤含水率、海拔、土層是影響土壤生態(tài)化學(xué)計(jì)量的主要因素，分別解釋了變異的44.0%、15.1%、9.7%(圖3)，土壤含水率的解釋率遠(yuǎn)大于海拔與土層。

表6 土壤化學(xué)計(jì)量特征與環(huán)境因子的相關(guān)性

3 結(jié)論與討論

海拔對青海省祁連圓柏天然林土壤SOC、TN、TP及化學(xué)計(jì)量比具有顯著影響。隨海拔的升高與土層的加深，土壤SOC、TN、TP、C∶P、N∶P呈下降趨勢；C∶N無顯著的變化。不同海拔梯度上氣候、土壤理化性質(zhì)、植被因素的差異引起了土壤化學(xué)計(jì)量的變化，其中土壤含水率是影響土壤化學(xué)計(jì)量變化的主控因素，解釋了變異的44.0%。青海祁連圓柏天然林土壤P有效性低，N有效性高。

3.1 祁連圓柏天然林土壤C、N、P含量及化學(xué)計(jì)量比的海拔變化規(guī)律及影響因素

本研究中土壤SOC、TN變異系數(shù)相近，表現(xiàn)出良好的協(xié)同變化特征，均隨海拔的增加而減小，但與王艷麗等[12]調(diào)查的青海森林土壤SOC、TN垂直分布規(guī)律不一致?？赡艿脑蛟谟谘芯繉ο蟮牟町?。王艷麗等[12]研究了隨海拔梯度升高由落葉闊葉林過渡至針葉林的植被類型林地土壤生態(tài)化學(xué)計(jì)量特征的變化，而本研究的對象則為單一的祁連圓柏林。與闊葉林相比，針葉林土壤微生物胞外酶的活性顯著降低[23]，針葉難易分解[5,24]，因此2 900～3 100 m海拔區(qū)域土壤SOC、TN含量并未增加。盡管土壤C、N含量存在較大的空間變異性，但是土壤C、N比在不同海拔梯度間卻無顯著變化。這是因?yàn)镃、N作為土壤的結(jié)構(gòu)成分，在積累與消耗過程中，對外界環(huán)境變化的響應(yīng)是一致的[25]。

不同海拔梯度下土壤SOC、TN含量變化主要?dú)w因于氣候、植被因素與土壤理化性質(zhì)的變化(圖3，表6)。較大尺度范圍內(nèi)氣候是限制植被生長的重要因素[25]。自然條件下，土壤C、N含量受生物作用的影響，主要由土壤中有機(jī)質(zhì)和N輸入量與輸出量的相對大小決定[26]。研究區(qū)<2 900 m海拔區(qū)域降雨量高、溫度適宜，水溫條件好使得群落多樣性、物種豐度與林分密度增加，進(jìn)而使輸入至林下凋落物含量增加，凋落物為土壤輸送了更多的養(yǎng)分，從而促進(jìn)了土壤C、N含量的積累。不同海拔梯度下土壤理化性質(zhì)的變化也導(dǎo)致土壤C、N含量的變化。本研究表明，土壤SOC、TN含量與容重、pH呈極顯著負(fù)相關(guān)(P<0.01)，這是由于較高的容重與pH不利于土壤通氣，限制了土壤微生物活性和根系生長，不適宜有機(jī)質(zhì)分解與養(yǎng)分匯集[27]。土壤SOC、TN含量與土壤含水率呈極顯著正相關(guān)。已有研究表明，土壤濕度較高條件下，土壤微生物和酶的活性會(huì)提高[32]，有機(jī)質(zhì)礦化分解加快[28]。本研究中隨海拔的升高，土壤含水率降低，微生物與酶活性降低，有機(jī)質(zhì)分解與周轉(zhuǎn)周期加長，因而，土壤SOC、TN含量降低。

注:H.Shannon-Winner指數(shù); S.物種豐富度; SD.林分密度; Alt.海拔; T.溫度; P.降水P; SWC.土壤含水率; BD.容重; Lay.土層。

本研究中土壤TP變異系數(shù)為16.13%，顯著低于SOC與TN，空間變異性低。這是因?yàn)橥寥繮主要來源于長期而穩(wěn)定的巖石風(fēng)化過程，屬沉積性元素，遷移性差[29]。另外，土壤TP隨海拔升高而逐漸增加。土壤P含量決定于巖石分化形成的土壤母質(zhì)，研究區(qū)祁連圓柏林下土壤多為山地灰褐土，土壤母質(zhì)主要為黃土及殘積物[15]，P來源有限；同時(shí)，高海拔地區(qū)，pH增大，氧化還原電位升高，土壤中金屬離子對P的吸附作用增強(qiáng)[28]；加之容重增加與含水率的降低削弱了土壤風(fēng)化，土壤母質(zhì)中P的釋放量減少[30]，進(jìn)而造成土壤TP含量的降低。就土壤C/P與N/P比而言，隨海拔的升高，C/P與N/P均呈下降趨勢。本研究中土壤C∶P變化由SOC變化引起，TN變化導(dǎo)致N∶P的變化(表3)，海拔升高土壤C、N含量的下降解釋了這一變化。

3.2 土壤化學(xué)計(jì)量比的指示作用

土壤環(huán)境中C、N、P化學(xué)計(jì)量比值決定著植物和微生物對養(yǎng)分的可利用性與限制性。N、P通常是陸地生態(tài)系統(tǒng)常見的限制性元素。一般認(rèn)為土壤C∶P是土壤P有效性的標(biāo)志，低的C∶P有利于微生物通過分解過程釋放養(yǎng)分，促進(jìn)土壤有效P的增加。反之，較高的C∶P使土壤微生物對有機(jī)質(zhì)的分解受限，同時(shí)還會(huì)與植物競爭土壤無機(jī)磷[31]。土壤N∶P可用作N飽和的指標(biāo)，反映植被生長過程中土壤養(yǎng)分元素的有效性，并用來確定養(yǎng)分限制的閾值[4]。研究區(qū)土壤C∶P、N∶P平均值為57.96與5.43，高于中國土壤平均水平[29]，說明青海祁連圓柏天然林土壤P有效性低，N有效性高。同時(shí)，本研究土壤TP含量(0.62 g·kg)低于中國土壤平均水平(0.78 g·kg-1)[29]，從另一個(gè)方面驗(yàn)證了研究區(qū)土壤的P限制。冗余分析單獨(dú)效應(yīng)表明土壤含水率解釋了土壤化學(xué)計(jì)量比變化的44.0%，占所有環(huán)境因子總解釋率的58%。有研究指出隨著全球氣候變暖，土壤水分因蒸發(fā)量變多而變得比以前更加干旱[32]。本研究中土壤TP與含水率呈極顯著正相關(guān)關(guān)系，青海地處青藏高原東北部，對氣候變化敏感，未來在全球氣候變暖的背景下，土壤水分下降會(huì)進(jìn)一步加劇土壤P限制。因此，在天然林經(jīng)營與管理中應(yīng)注重緩解土壤P限制，如P肥的施用。