楊 霞， 陳麗華， 鄭學(xué)良

(北京林業(yè)大學(xué)水土保持學(xué)院，100083，北京)

生態(tài)化學(xué)計(jì)量學(xué)以生態(tài)學(xué)和化學(xué)計(jì)量學(xué)為基本原理，研究生物系統(tǒng)能量平衡和多種化學(xué)元素平衡，其為探索碳(C)、氮(N)和磷(P)等主要元素的地化循環(huán)和生態(tài)學(xué)過程提供了一種新思路[1-2]。

土壤C、N和P質(zhì)量分?jǐn)?shù)的變化可以反映土壤養(yǎng)分狀況，土壤C∶N、C∶P和N∶P計(jì)量特征不僅能體現(xiàn)土壤有機(jī)質(zhì)組成、土壤質(zhì)量狀況以及土壤養(yǎng)分供給能力，還可以指示土壤各養(yǎng)分的礦化能力和固持能力[3]。姜沛沛等[4]對(duì)子午嶺林區(qū)不同林齡(9、23、33和47年)油松(Pinustabuliformis)人工林土壤的研究表明，土壤有機(jī)C和全N含量隨林齡的變化規(guī)律為先增加后降低，P含量隨林齡先降低后增加。汪宗飛等[5]在子午嶺地區(qū)的研究結(jié)果為油松林土壤有機(jī)C和全N含量隨林齡增加而增加，P含量隨林齡逐漸減小。另外有學(xué)者研究了廣西憑祥不同林齡(6、16、23和35年)馬尾松(Pinusmassoniana)人工林土壤C、N和P變化特征，結(jié)果顯示：土壤有機(jī)C、全N和全P含量隨林齡增加先降低后逐漸增加，土壤有機(jī)C在不同林齡和土層之間的差異是影響土壤C∶N和C∶P變化的主要因素，N則是導(dǎo)致土壤N∶P變化的主要原因[6]；由此可見，人工林生態(tài)系統(tǒng)土壤肥力變化特征具有不一致性，這種不一致性是由區(qū)域氣候、植被生物學(xué)特性及人工林經(jīng)營(yíng)方式等造成的[7-8]。生態(tài)化學(xué)計(jì)量比研究可以為人工林經(jīng)營(yíng)過程中養(yǎng)分限制或養(yǎng)分循環(huán)提供重要的判斷作用，因此在特定地區(qū)和植被條件下開展土壤生態(tài)化學(xué)計(jì)量比研究是非常有必要的[9]。

油松作為我國的特有樹種，不僅是北方主要的造林用材樹種，也是遼東地區(qū)主要的人工水源涵養(yǎng)林樹種，在林業(yè)生產(chǎn)和森林生態(tài)系統(tǒng)中占有十分重要的地位。目前，國內(nèi)學(xué)者在黃土高原地區(qū)做了大量關(guān)于土壤養(yǎng)分特征的研究，這些研究對(duì)林分實(shí)際經(jīng)營(yíng)起到了一定的理論指導(dǎo)作用。本研究區(qū)位于遼寧省撫順縣溫道林場(chǎng)，北靠大伙房水庫，油松作為該區(qū)主要的人工水源涵養(yǎng)林之一，發(fā)揮著重要的生態(tài)作用。但近年全球氣溫回升、降雨量減少、蒸騰量增加，森林生態(tài)環(huán)境逐漸惡化，部分森林功能開始退化，林木生長(zhǎng)出現(xiàn)倒退現(xiàn)象[10]；因此，有必要對(duì)該區(qū)油松人工林土壤C、N和P質(zhì)量分?jǐn)?shù)及生態(tài)化學(xué)特征進(jìn)行研究。通過對(duì)不同林齡油松土壤C、N和P質(zhì)量分?jǐn)?shù)的測(cè)定和分析，探討不同林齡和土層下土壤C、N和P質(zhì)量分?jǐn)?shù)以及化學(xué)計(jì)量比的變化規(guī)律和異同，指導(dǎo)油松林木生長(zhǎng)環(huán)境的調(diào)節(jié)，為油松人工林土壤肥力的管理提供參考。

1 研究區(qū)概況

研究區(qū)位于遼寧省大伙房水庫上游社河流域的溫道林場(chǎng)，東臨新賓縣界，西接撫順市市郊，南連后安鎮(zhèn)，北靠大伙房水庫，地理坐標(biāo)：E 124°09′54.65″、N 41°49′23.50″，屬于長(zhǎng)白山脈龍崗支脈西麓的低山區(qū)，海拔200～500 m，坡度15～25°。氣候類型為溫帶季風(fēng)型大陸氣候，年均降雨量800～1 000 mm，多集中在6—8月，年均溫5～8 ℃，土壤主要為棕色森林土和暗棕色森林土，土質(zhì)較疏松，土壤有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)較高。該流域水源涵養(yǎng)林以人工林為主，此外有少量的天然次生林，森林覆蓋率為85%，主要樹種包括油松、落葉松(Larixgmelinii)、紅松(Pinuskoraiensis)、樟子松(Pinussylvestrisvar.mongolica)、刺槐(Robiniapseudoacacia)、山楊(Populusdavidiana)、蒙古櫟(Quercusmongolica)、紫穗槐(Amorphafruticosa)等[11-12]。

2 研究方法

2.1 樣地設(shè)置

于2018年6月在遼寧省撫順縣溫道林場(chǎng)進(jìn)行調(diào)查。經(jīng)實(shí)地調(diào)研并咨詢林場(chǎng)工作人員后，選取3個(gè)不同林齡油松人工純林(42、65和85年)作為研究對(duì)象，每個(gè)林齡布設(shè)3塊面積為20 m×20 m的樣地(共9塊樣地)，同齡林樣地之間的距離不超過1 km，保持立地條件基本一致。分別對(duì)各樣地進(jìn)行每木檢尺、生物多樣性調(diào)查，記錄海拔、坡度、坡向、郁閉度、樹高、胸徑、冠幅等。樣地基本情況見表1。

2.2 樣品采集及處理

按照坡上、坡中、坡下線路在上述布設(shè)的9個(gè)油松樣地中，各取3個(gè)土樣采集點(diǎn)(共27個(gè)采樣點(diǎn))，采樣點(diǎn)距邊界5 m以上。在采樣點(diǎn)挖土壤剖面，按照0～20、>20～40和>40～60 cm土層(根據(jù)試驗(yàn)區(qū)的立地條件，60 cm左右土層大石塊較多，故本次實(shí)驗(yàn)樣品均取到60 cm)分別取1 kg左右的土樣，最后將每個(gè)樣地3個(gè)采樣點(diǎn)相同土層的土樣混合均勻各取1 kg，風(fēng)干所取土樣，除去大的石塊和動(dòng)植物殘?bào)w，過2 mm土壤篩，并用四分法取部分風(fēng)干土樣[11]，研磨過0.149 mm土壤篩，以便測(cè)土壤養(yǎng)分。

表1 油松樣地基本特征Tab.1 Basic characteristics in the plots of Pinus tabuliformis

2.3 土壤養(yǎng)分測(cè)定

土壤有機(jī)C測(cè)定采用重鉻酸鉀氧化-濃硫酸外加熱法，土壤全N采用K1100全自動(dòng)凱氏定氮儀法，土壤全P采用氫氧化鈉熔融-鉬銻抗比色法[13-14]。本研究結(jié)果取各土樣3次重復(fù)測(cè)定的平均值。

2.4 數(shù)據(jù)分析

數(shù)據(jù)統(tǒng)計(jì)處理采用Microsoft Excel 2010。數(shù)據(jù)分析采用SPSS 22.0，用單因素方差分析(One-way ANOVA)和多重比較法分析油松人工林不同林齡或不同土層深度之間土壤有機(jī)C、全N和全P質(zhì)量分?jǐn)?shù)和化學(xué)計(jì)量比的差異，雙因素方差分析法(Two-way ANOVA)分析林齡和土層深度及其交互作用對(duì)土壤化學(xué)性質(zhì)和化學(xué)計(jì)量比的影響，用Pearson相關(guān)性分析化學(xué)計(jì)量比與土壤有機(jī)C、全N和全P的相關(guān)性。采用Microsoft Excel 2010作圖。

3 結(jié)果與分析

3.1 土壤碳、氮和磷質(zhì)量分?jǐn)?shù)特征

由方差分析可得(表2)，林齡和土層深度對(duì)土壤有機(jī)C質(zhì)量分?jǐn)?shù)均存在極顯著影響(P<0.01)，林齡和土層深度的交互作用對(duì)土壤有機(jī)C質(zhì)量分?jǐn)?shù)無顯著影響(P>0.05)。3個(gè)不同林齡油松林地0～60 cm層土壤有機(jī)C質(zhì)量分?jǐn)?shù)的均值分別為47.67 g/kg(65年)>22.74 g/kg(85年)>15.73 g/kg(42年)，即土壤有機(jī)C質(zhì)量分?jǐn)?shù)隨林齡增加先增加后逐漸減少。65年油松林地土壤有機(jī)C質(zhì)量分?jǐn)?shù)明顯高于42和85年，并且高于3個(gè)林地土壤有機(jī)C質(zhì)量分?jǐn)?shù)的平均值(28.71 g/kg)。另外，0～20和>40～60 cm土層內(nèi)有機(jī)C質(zhì)量分?jǐn)?shù)隨林齡的變化規(guī)律與上述規(guī)律相同，且各土層內(nèi)65年林地土壤有機(jī)C質(zhì)量分?jǐn)?shù)與其他2個(gè)林地差異顯著(P<0.05)。土壤有機(jī)C質(zhì)量分?jǐn)?shù)隨土層加深逐漸減少，其中0～20 cm層明顯高于其他土層且差異顯著(P<0.05)(圖1)。

表2 林齡和土層深度對(duì)油松人工林土壤碳、氮、和磷生態(tài)化學(xué)計(jì)量比的影響Tab.2 Effects of stand age and soil depth on soil C, N, and P ecological stoichiometry of Pinus tabuliformis plantations

同一林齡不同小寫字母表示不同土層間差異顯著(P<0.05)，同一土層不同大寫字母表示不同林齡間差異顯著(P<0.05)。下同。Different lowercase alphabets of the same stand age indicate significant differences in different soil layers (P<0.05), different capital alphabets of the same soil layer indicate significant differences among different stand ages (P<0.05). The same below. 圖1 不同林齡油松人工林土壤全C、全N、全P質(zhì)量分?jǐn)?shù)Fig.1 Soil total C, total N, and total P content in different stand age of Pinus tabuliformis plantations

林齡、土層深度以及林齡和土層深度的交互作用對(duì)土壤全N和全P質(zhì)量分?jǐn)?shù)均存在極顯著影響(P<0.01)(表2)。3個(gè)林齡油松林地0～60 cm土壤全N質(zhì)量分?jǐn)?shù)均值分別為0.40 g/kg(42年)、1.58 g/kg(65年)和0.57 g/kg(85年)，N質(zhì)量分?jǐn)?shù)隨林齡先增加后減少，這與土壤有機(jī)C質(zhì)量分?jǐn)?shù)變化規(guī)律是一致的。65年油松林地土壤全N質(zhì)量分?jǐn)?shù)均值明顯高于42和85年。0～20、>20～40和>40～60 cm土層內(nèi)土壤全N質(zhì)量分?jǐn)?shù)的變化與N質(zhì)量分?jǐn)?shù)在0～60 cm土層內(nèi)的變化是一致的，且在0～20和>20～40 cm土層內(nèi)，65年林地N質(zhì)量分?jǐn)?shù)明顯區(qū)別與其他2個(gè)林地(P<0.05)。另外，不同林地土壤全N質(zhì)量分?jǐn)?shù)均隨土層加深減少，0～20 cm土壤全N質(zhì)量分?jǐn)?shù)均顯著高于其他土層(P<0.05)。

不同林地0～60 cm土壤全P質(zhì)量分?jǐn)?shù)均值為2.27 g/kg(85年)>1.21 g/kg(65年)>1.11 g/kg(42年)，P質(zhì)量分?jǐn)?shù)隨林齡增加，這與上述土壤有機(jī)C和土壤全N質(zhì)量分?jǐn)?shù)隨林齡變化的規(guī)律稍有不同。且85年林地土壤全P質(zhì)量分?jǐn)?shù)均值明顯高于其他2個(gè)林地，約為其的1.80～2.20倍，其他林地則差別不明顯。在不同土層內(nèi)，85年林地土壤全P質(zhì)量分?jǐn)?shù)也均顯著高于42年、65年的(P<0.05)。另外，各林地土壤全P質(zhì)量分?jǐn)?shù)隨土層加深略有減少，除85年林地外，其他林地各土層間無顯著差異(P>0.05)。

圖2 不同林齡油松人工林土壤碳∶氮∶磷生態(tài)化學(xué)計(jì)量比特征Fig.2 Ecological stoichiometry of soil C∶N∶P in different stand age of Pinus tabuliformis plantations

3.2 土壤碳、氮和磷化學(xué)計(jì)量特征

根據(jù)方差分析結(jié)果，林齡和土層深度對(duì)土壤C∶N無顯著影響(P>0.05)，而林齡和土層深度的交互作用對(duì)土壤C∶N有極顯著影響(P<0.01)(表2)。不同林齡林地0～60 cm層土壤C∶N變化范圍為24.54～57.02，C∶N隨林齡先減小后增加，即42.43(85年)>37.30(42年)>36.19(65年)，各林地土壤C∶N均值差別不大。0～20 cm層土壤C∶N表現(xiàn)為42.67(42年)>37.67(85年)>24.54(65年)；>20～40 cm為53.41(85年)>42.56(42年)>27.02(65年)；>40～60 cm為57.02(65年)>36.22(85年)>26.67(42年)，且65年林地與其他2個(gè)林地差異顯著(P<0.05)。

由表2可知，林齡和土層深度對(duì)土壤C∶P和N∶P均有極顯著影響(P<0.01)。不同林地0～60 cm土層土壤C∶P范圍為4.27～55.26，其值隨林齡先增加后減小，表現(xiàn)為39.44(65年)>15.47(42年)>10.40(85年)，65年林地明顯高于其他林地。3個(gè)土層內(nèi)土壤C∶P的變化規(guī)律與上述規(guī)律一致，且在各土層中65年林地與42和85年林地均存在顯著差異(P<0.05)。另外，不同林地土壤C∶P隨土層加深逐漸降低，其中0～20 cm與其他2個(gè)土層均存在顯著差別(P<0.05)。3個(gè)油松林地0～60 cm層土壤N∶P為0.12～2.25，表現(xiàn)為隨林齡先增加后減小，即1.33(65年)>0.38(42年)>0.33(85年)，65年林地高于其他2個(gè)林地。從不同土層來看，土壤N∶P隨土層加深逐漸降低，各層土壤N∶P均值表現(xiàn)為1.16(0～20 cm)>0.54(>20～40 cm)>0.28(>40～60 cm)，且0～20 cm土層與>20～40 cm、>40～60 cm土層之間均有顯著差異(P<0.05)。土壤N∶P的變化規(guī)律與上述土壤C∶P的變化規(guī)律基本一致(圖2)。

3.3 土壤碳、氮和磷質(zhì)量分?jǐn)?shù)與生態(tài)化學(xué)計(jì)量比的相關(guān)性

根據(jù)Pearson雙變量相關(guān)性分析(表3)，土壤C∶N與土壤有機(jī)C無顯著相關(guān)性，與土壤全N質(zhì)量分?jǐn)?shù)極顯著相關(guān)。土壤C∶P與土壤有機(jī)C質(zhì)量分?jǐn)?shù)、土壤全N和土壤全P質(zhì)量分?jǐn)?shù)均存在極顯著的相關(guān)性。土壤N∶P與土壤全N和全P質(zhì)量分?jǐn)?shù)也均存在極顯著相關(guān)性。研究結(jié)果表明研究區(qū)油松人工林土壤C∶N、C∶P和N∶P與土壤有機(jī)C、全N和全P質(zhì)量分?jǐn)?shù)有非常密切的關(guān)系。

表3 油松人工林土壤碳、氮、和磷質(zhì)量分?jǐn)?shù)及生態(tài)化學(xué)計(jì)量比的相關(guān)性Tab.3 Correlation among soil C, N, P content and their ecological stoichiometry of Pinus tabuliformis plantations

4 討論

4.1 土壤碳、氮和磷質(zhì)量分?jǐn)?shù)特征

土壤有機(jī)C、全N和全P元素作為植物正常生長(zhǎng)發(fā)育必需的養(yǎng)分，在植物生長(zhǎng)過程中發(fā)揮著重要的作用[15]，其質(zhì)量分?jǐn)?shù)的多少以及組合成分，受植物種類、經(jīng)營(yíng)方式以及氣候和土壤等諸多因素的影響[16]。本研究結(jié)果顯示，油松人工林土壤有機(jī)C、全N質(zhì)量分?jǐn)?shù)隨林齡增加呈現(xiàn)先增加后減少的變化規(guī)律，這與姜沛沛等[4]、陳立新等[17]的研究結(jié)果一致。出現(xiàn)這種結(jié)果的原因可能是林下枯落物的累計(jì)以及林木對(duì)土壤養(yǎng)分的需求發(fā)生變化。土壤養(yǎng)分(C、N和P)主要來源于林下枯落物的分解和表層根系的分泌物[18]，隨著林齡的增加，林木生長(zhǎng)越來越旺盛，林下枯落物不斷積累，林木的根系也逐漸擴(kuò)大，這一時(shí)期大部分養(yǎng)分以枯落物等形式歸還土壤。在林齡達(dá)到一定階段之后，林木的生長(zhǎng)出現(xiàn)倒退現(xiàn)象，土壤中的養(yǎng)分也逐漸下降。一些學(xué)者的研究結(jié)果與上述研究稍有不同[19-20]，查閱文獻(xiàn)發(fā)現(xiàn)這些學(xué)者研究時(shí)所選取的油松林年限都較小，進(jìn)而猜想得出土壤有機(jī)C和全N質(zhì)量分?jǐn)?shù)隨林齡增加而增加的結(jié)論。另外，與以上2種結(jié)論不同的是，雷麗群等[6]對(duì)馬尾松人工林的研究表明，土壤有機(jī)C和全N質(zhì)量分?jǐn)?shù)隨林齡增加降低，至成熟林階段有所回升，曾凡鵬等[21]對(duì)遼東落葉松的研究結(jié)果為，土壤有機(jī)C和全N質(zhì)量分?jǐn)?shù)隨林齡增加降低。由此可猜想，林木在生長(zhǎng)過程中土壤C、N和P質(zhì)量分?jǐn)?shù)的變化可能出現(xiàn)多峰現(xiàn)象，而非僅單峰變化，這種現(xiàn)象可能是由林木在不同生長(zhǎng)時(shí)期對(duì)養(yǎng)分需求不同造成的。研究還發(fā)現(xiàn)土壤有機(jī)C、全N質(zhì)量分?jǐn)?shù)隨土層加深有明顯降低的趨勢(shì)，土壤C和N主要聚集在0～40 cm土層，這與以往學(xué)者的研究結(jié)果一致[22-25]。原因可能是枯落物聚集在土壤表層，微生物活性較高，枯落物分解速率較快，土壤較肥沃[5]。研究中土壤全P質(zhì)量分?jǐn)?shù)隨林齡增加，這與以往結(jié)果相似[18,26]。研究中42年油松林土壤全P質(zhì)量分?jǐn)?shù)為1.11 g/kg，與張?zhí)僮拥萚26]對(duì)遼西油松人工林(37 年)的研究結(jié)果(1.03 g/kg)基本相同。此外有學(xué)者在黃土高原的研究結(jié)果為0.13～0.61 g/kg，該研究顯示土壤P質(zhì)量分?jǐn)?shù)較低[27-28]。究其原因，主要是地域差異，土壤類型本身不同。其次是因?yàn)辄S土高原水土流失嚴(yán)重，土壤養(yǎng)分淋溶大.土壤P質(zhì)量分?jǐn)?shù)相對(duì)較。不同林齡油松林各土層土壤全P質(zhì)量分?jǐn)?shù)無顯著差異，造成這種現(xiàn)象的主要原因是土壤中P元素大多來自巖石風(fēng)化和淋溶，巖石風(fēng)化需要一個(gè)長(zhǎng)期的過程，在這個(gè)過程中土壤0～60 cm土層內(nèi)巖石風(fēng)化程度并無太大差別，并且P元素是沉積性的礦物，土壤中遷移相對(duì)較少，因此土壤P素質(zhì)量分?jǐn)?shù)在各土層的空間分布較穩(wěn)定[29]。

4.2 土壤碳、氮和磷生態(tài)化學(xué)計(jì)量特征

土壤生態(tài)化學(xué)計(jì)量比是評(píng)價(jià)土壤質(zhì)量的重要指標(biāo)，是確定土壤C、N和P平衡特征的重要參數(shù)[30]，對(duì)植物生長(zhǎng)有重要影響[28]。土壤C∶N可以反映土壤有機(jī)質(zhì)的分解狀況，是土壤N素礦化能力的標(biāo)志[31]。土壤C∶N小于25，不僅利于有機(jī)質(zhì)的轉(zhuǎn)化，而且超過生物生長(zhǎng)所需要的N素也會(huì)被釋放到土壤中，此時(shí)土壤N素充足；當(dāng)C∶N大于25，有機(jī)質(zhì)難以轉(zhuǎn)化，同時(shí)生物需要輸入N素來滿足自身的生長(zhǎng)需求，利于有機(jī)質(zhì)積累[32-33]。研究中3個(gè)林齡油松林土壤C∶N為36.19～42.43，C∶N隨林齡變化有所增加，但增加規(guī)律不明顯，林齡和土層深度對(duì)土壤C∶N無明顯影響。研究結(jié)果遠(yuǎn)大于最佳土壤C∶N值(25)，也高于中國陸地土壤平均水平(11.90)[34]，這說明研究區(qū)油松林土壤中有機(jī)質(zhì)難以轉(zhuǎn)化，土壤N素礦化能力較弱，土壤N素不足。王棣等[35]在秦嶺地區(qū)的研究結(jié)果為5～20，該結(jié)果與本研究差異顯著，主要是由土壤類型、氣候等條件不同引起的。相關(guān)性分析可以看出(表3)，土壤C∶N與土壤有機(jī)C無明顯相關(guān)性，與土壤全N顯著相關(guān)，說明N素在不同林齡和土層之間的變化是土壤C∶N發(fā)生變化的主要原因。土壤中的微生物可以從環(huán)境中吸收P元素或者將P從土壤有機(jī)物質(zhì)中釋放出來，這種吸附、釋放P素的能力可根據(jù)土壤C∶P進(jìn)行判斷[5]。研究中土壤C∶P隨林齡先增加后降低，隨土層深度逐漸降低，二者對(duì)其均有顯著影響。結(jié)果顯示研究區(qū)油松林土壤C∶P為10.40～39.44，該數(shù)據(jù)與遼西地區(qū)研究(14.39)[26]基本一致。C∶P越小，越利于土壤微生物礦化土壤有機(jī)質(zhì)釋放較多P素，補(bǔ)充土壤有效P庫[5]。研究中土壤C∶P低于我國陸地土壤平均水平(61)[34]，說明研究區(qū)土壤微生物礦化土壤有機(jī)物質(zhì)釋放P素的能力較強(qiáng)，土壤P素充足。分析可知，土壤C∶P與土壤有機(jī)C、全P均顯著相關(guān)，說明土壤C、P變化是C∶P變化的主要原因。土壤N∶P作為判斷生境中N或P不足的依據(jù)，可判別植物生長(zhǎng)過程中土壤養(yǎng)分的供給狀態(tài)。研究中，不同林齡油松林土壤N∶P變化范圍為0.26～1.33，小于全國平均水平(5.20)[34]。林齡和土層深度對(duì)土壤N∶P均有顯著影響，但土壤中N∶P變化相對(duì)穩(wěn)定。

5 結(jié)論

1)林齡和土層深度對(duì)土壤有機(jī)C、全N、全P、C∶P、N∶P均有顯著性影響。土壤有機(jī)C、全N質(zhì)量分?jǐn)?shù)隨林齡先增加后降低，土壤全P質(zhì)量分?jǐn)?shù)隨林齡逐漸增加，各土壤養(yǎng)分隨林齡的變化可能出現(xiàn)多峰現(xiàn)象，而非僅單峰變化。隨土層加深，以上三者的質(zhì)量分?jǐn)?shù)均表現(xiàn)為下降的趨勢(shì)。

2)不同林齡和土層之間土壤C∶N主要受N素影響，C∶P受C和P影響，N∶P受N和P影響。

3)土壤有機(jī)C和全N主要聚集在0～40 cm土層，且該區(qū)土壤N素礦化能力較弱，土壤N素不足。