鄧博文，許瑤瑤，陳逸飛，張 碩，鄭德祥＊，廖曉麗，林仁忠，靳少非,4

(1. 福建農(nóng)林大學(xué)林學(xué)院，福建 福州 350002；2. 閩江學(xué)院地理科學(xué)系，福建 福州 350108；3. 順昌縣國有林場，福建 順昌 353200；4. 華東師范大學(xué)崇明生態(tài)研究院，上海 202162)

生態(tài)化學(xué)計量學(xué)是研究生態(tài)系統(tǒng)中養(yǎng)分平衡的科學(xué)，為研究植物-環(huán)境相互作用提供了良好的指標(biāo)和研究方向[1-3]。探究植物葉片N、P、K與環(huán)境因子的相關(guān)性，有利于挖掘植物關(guān)鍵生物元素循環(huán)規(guī)律，為研究植物的元素平衡吸收利用提供理論依據(jù)[4-5]。在氣候變化、植物生理生化作用、土壤地球化學(xué)分布特征、植物群落組成等多重因素影響下[6-8]，森林生態(tài)化學(xué)計量存在著緯向變化特征，如對全球452個點(diǎn)、1 280種植物、5 087個觀測值的研究發(fā)現(xiàn)，隨著緯度的增加，平均溫度降低，植物葉片N、P含量增加，N/P比下降[7]。對中國區(qū)域及中國東部南北樣帶植物N、P含量的研究結(jié)果也顯示，葉片N、P含量呈現(xiàn)隨緯度降低而減少的趨勢，但N/P無顯著變化[9-11]。目前，對于另一關(guān)鍵化學(xué)計量元素K與其他營養(yǎng)元素的生態(tài)化學(xué)計量分析未得到廣泛研究[12]。

針葉林在我國各氣候帶內(nèi)均有分布，共6科30屬約200種，其廣闊的跨氣候帶的分布特征為研究針葉林生態(tài)化學(xué)計量空間大尺度變化特征提供了研究素材[13]，對探究針葉樹種葉片化學(xué)計量特征緯向變化具有重要意義。植物葉片生態(tài)化學(xué)計量的緯向變化特征是多環(huán)境影響因素的整體效應(yīng)，溫度、水分、土壤養(yǎng)分等環(huán)境因子也對植物生態(tài)化學(xué)計量特征產(chǎn)生一定的影響，N、P、K是影響植物生長發(fā)育的重要指標(biāo)，反映出土壤生產(chǎn)力的強(qiáng)弱[14-15]。隨著溫度升高，植物葉片中N/P呈增加趨勢[7,11]，閆曉俊等[16]在對杉木幼苗增溫實(shí)驗(yàn)中發(fā)現(xiàn)，大氣和土壤單獨(dú)增溫均可增加葉片N含量，大氣和土壤同時增溫顯著增加了葉片P含量，但對N/P影響不顯著；在干旱區(qū)域，葉片的化學(xué)計量與降水呈顯著相關(guān)，而在熱帶多雨地區(qū)，并未發(fā)現(xiàn)顯著關(guān)系[17]。自20世紀(jì)以來，人工合成氮肥技術(shù)的發(fā)展導(dǎo)致氮肥廣泛應(yīng)用、化石燃料大規(guī)模開采利用，增加了大氣N沉降量[18]。劉文飛等[19]在對杉木進(jìn)行大氣氮沉降模擬實(shí)驗(yàn)發(fā)現(xiàn)，增加N輸入可提高植物葉片N含量，降低植物P、K含量，進(jìn)而增加了N/P、N/K；過量的N輸入改變了植物生態(tài)化學(xué)計量特征，表現(xiàn)為葉片N含量增加，增加N/P比，導(dǎo)致植物生長更易受到P限制作用[20-22]。土壤生態(tài)化學(xué)計量特征可以直接影響葉片生態(tài)化學(xué)計量特征，并通過內(nèi)穩(wěn)態(tài)調(diào)節(jié)機(jī)制應(yīng)對土壤養(yǎng)分的變化[23-25]。目前，對于探究多環(huán)境因子對葉片化學(xué)計量變化的研究較少，因此，探究不同環(huán)境因子對植物葉片化學(xué)計量特征影響的貢獻(xiàn)及植物內(nèi)穩(wěn)態(tài)狀態(tài)對于綜合了解葉片生態(tài)化學(xué)計量特征具有重要的意義。

本研究選取中國區(qū)域針葉林優(yōu)勢樹種為對象，分析中國區(qū)域優(yōu)勢針葉林土壤－植物N、P、K含量以及化學(xué)計量特征，探究不同環(huán)境因子對葉片化學(xué)計量特征的影響，分析中國針葉樹種葉片化學(xué)計量穩(wěn)態(tài)性特征，該研究結(jié)果可為綜合研究針葉林生物地球化學(xué)循環(huán)特征與全球環(huán)境變化響應(yīng)提供數(shù)據(jù)基礎(chǔ)和理論支持。

1 材料和方法

1.1 數(shù)據(jù)來源與處理

本研究中所使用的數(shù)據(jù)均來自于已有文獻(xiàn)數(shù)據(jù)[26-33]，經(jīng)過匹配，共81個取樣點(diǎn)的土壤養(yǎng)分?jǐn)?shù)據(jù)被納入，計262個土壤剖面，每個剖面設(shè)置3個土壤深度：0～20、20～40、40～60 cm。

本研究中收集了包括西伯利亞落葉松(Larix sibiricaLedeb.)、黃花落葉松(L. olgensisHenry.)、魚鱗云杉(Picea jezoensisCheng.)、臭冷杉(AbiesnephrolepisMaxim.)、石松(Lycopodium clavammThunb.)、扁枝石松(L. complanatumHolub.)、玉柏石松(L. obscurumLinn.)、杉蔓石松(L. annotinumL.)、赤松(Pinus densifloraSieb.)、西伯利亞松(P.sylvestrisLitv.)、青海云杉(P. crassifoliaKom.)、油松(P. tabuliformisCarriere.)、白皮松(P. bungeanaZucc.)、馬尾松(P. massonianaLamb.)、云南松(P.yunnanensisFranch.)、 華 山 松 (P. armandiiFranch.)、杉木(Cunninghamia lanceolata(Lamb.)Hook.)、側(cè)柏(Platycladus orientalisFranco.)、柏木(Cupressus funebrisEndl.)、干香柏(C. duclouxianaHichel.)20種具代表性的針葉林樹種。

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

1.2.1 內(nèi)穩(wěn)態(tài)指數(shù) 內(nèi)穩(wěn)態(tài)指數(shù)（H）采用下式計算：

式中：x為土壤生態(tài)化學(xué)計量特征(如N、N/P等)，y為葉片生態(tài)化學(xué)計量，c為擬合常數(shù)。1/H表示log10(x)和log10(y)回歸斜率，其絕對值范圍為0～1，生物體內(nèi)穩(wěn)態(tài)越穩(wěn)定越高，其H值越大。

對內(nèi)穩(wěn)態(tài)按照穩(wěn)定性進(jìn)行分類級別劃分[34]，回歸分析顯著性檢驗(yàn)采用α= 0.1進(jìn)行單側(cè)檢驗(yàn)，若回歸關(guān)系不顯著(P＞ 0.1)，被定義為“絕對穩(wěn)態(tài)”；若回歸關(guān)系顯著，當(dāng)1/H= 1，則植物體不具有內(nèi)穩(wěn)態(tài)性質(zhì)；0 ＜ 1/H＜ 1的數(shù)據(jù)集分別定義為：0 ＜1/H＜ 0.25為“穩(wěn)態(tài)”，0.25 ＜ 1/H＜ 0.5為“弱穩(wěn)態(tài)”，0.5 ＜ 1/H＜ 0.75為“弱敏感態(tài)”，1/H＞ 0.75為“敏感態(tài)”。

1.2.2 環(huán)境因子相關(guān)性 為探究不同環(huán)境因子對葉片生態(tài)化學(xué)計量的影響，本研究利用結(jié)構(gòu)方程模型分析緯度、經(jīng)度、年降水量、年均氣溫、最暖月均氣溫、最冷月均氣溫、年均氣溫差、水熱指數(shù)、年蒸發(fā)量、水分虧欠指數(shù)、土壤不同層次生態(tài)化學(xué)計量等參數(shù)對中國優(yōu)勢針葉樹種葉片化學(xué)計量的影響。在分析過程中，設(shè)置4個潛變量，分別為溫度相關(guān)變量、降水相關(guān)變量、地理區(qū)域變量、土壤養(yǎng)分供應(yīng)；1個顯變量為水分虧欠指數(shù)，分析不同變量對葉片化學(xué)計量影響。具體的模型見圖1。

在結(jié)構(gòu)方程模型分析的基礎(chǔ)上，利用R軟件relaimpo包計算各個變量對葉片化學(xué)計量的相對貢獻(xiàn)；采用One-way ANOVA分析不同土壤層次中N、P、K、N/P、N/K、K/P等化學(xué)計量學(xué)特征；利用一元線性回歸分析針葉林葉片中N、P、K、N/P、N/K、K/P與土層次中的相關(guān)關(guān)系，并分析其穩(wěn)態(tài)性。研究中所采用的統(tǒng)計軟件為R (R version 3.5.3)。

圖 1 結(jié)構(gòu)方程模型結(jié)構(gòu)Fig. 1 The frame of structure equation model

2 結(jié)果與分析

2.1 中國針葉林優(yōu)勢樹種葉片及土壤生態(tài)化學(xué)計量分析

中國區(qū)域針葉林葉片N、P、K平均含量分別為13.16、1.02、6.82 g·kg-1，N/P、N/K、K/P平均比分別為15.48、2.17、7.77。土壤中N、N/P、N/K在不同土層內(nèi)存在差異（圖2），其中，土層0～20 cm N、P、K含量及N/P、N/K最高；0～20、20～40、40～60 cm土層的N平均含量分別為1.73、0.99、0.72 g·kg-1；0～20、20～40、40～60 cm土層的N/P分別為3.10、2.12、1.74，N/K分別為0.10、0.06、0.05。

2.2 不同環(huán)境要素對中國針葉林優(yōu)勢樹種葉片生態(tài)化學(xué)計量貢獻(xiàn)分析

結(jié)構(gòu)方程模型表明：多種環(huán)境因子與我國優(yōu)勢針葉林樹種葉片化學(xué)計量特征存在顯著關(guān)系（圖1）。由表1可知：緯度對葉片N含量貢獻(xiàn)率最高，為19.18%；水熱指數(shù)對葉片P含量貢獻(xiàn)率最高，為12.75%；緯度對葉片K含量貢獻(xiàn)率最高，為25.08%；對N/P、N/K、P/K生態(tài)化學(xué)計量貢獻(xiàn)率最高的分別為年降雨量（貢獻(xiàn)率15.65%）、土壤40～60 cm深度N/K值（貢獻(xiàn)率16.76%）及年降雨量（貢獻(xiàn)率18.21%）。此外，土壤因子中，N、P、K、N/P、N/K、K/P對葉片化學(xué)計量貢獻(xiàn)率分別為17.01%、11.20%、8.35%、16.27%、30.79%、27.52%。

2.3 中國針葉林優(yōu)勢樹種葉片化學(xué)計量穩(wěn)態(tài)性特征

對我國優(yōu)勢針葉樹葉片中化學(xué)計量學(xué)特征隨土壤環(huán)境的內(nèi)穩(wěn)態(tài)特征進(jìn)行分析，結(jié)果（表2）表明：葉片P、K、N/P、N/K、P/K利用穩(wěn)態(tài)性模型模擬結(jié)果不顯著，屬絕對穩(wěn)態(tài)型；而葉片N內(nèi)穩(wěn)態(tài)指數(shù)為0.261，屬于弱穩(wěn)態(tài)型。

3 討論

3.1 中國針葉林優(yōu)勢樹種葉片生態(tài)化學(xué)計量分析

本研究利用文獻(xiàn)資料探究了中國針葉林葉片N、P、K含量及其生態(tài)化學(xué)計量學(xué)特征，結(jié)果表明，中國針葉林優(yōu)勢樹種葉片N ＞ K ＞ P，這與之前對興安落葉松（Larix gmelinii Rupr.）及東部針葉樹、黃土高原地區(qū)葉片養(yǎng)分含量的研究結(jié)果相似（表3）[13，35-36]。由于數(shù)據(jù)來源具有差異性，本研究系統(tǒng)對比了其他學(xué)者研究結(jié)果發(fā)現(xiàn)，葉片K含量及與之相對應(yīng)的研究過少，因此，對N、P與之相對應(yīng)的化學(xué)計量重點(diǎn)研究。

植物葉片中N和P含量是研究關(guān)鍵生境要素相對限制性的重要指標(biāo)。依據(jù)已有的文獻(xiàn)指標(biāo)，當(dāng)N/P ＜ 14時，植物生長受到N限制；當(dāng)N/P ＞16時，則主要受到P限制；當(dāng)N/P介于14～16時，同時受到N、P限制[37]，這也是我國總體樣帶植物的葉片化學(xué)計量特征[9]。植物營養(yǎng)元素含量的多寡受到土壤中營養(yǎng)條件的影響，并表現(xiàn)為對土壤營養(yǎng)元素環(huán)境的適應(yīng)性[1-2,38]。

圖 2 中國優(yōu)勢針葉樹種不同土壤發(fā)生層中N、K、P含量以及化學(xué)計量分析Fig. 2 Box plots for the contents of N, P, and K, and their stoichiometry in different soil layers in China’s coniferous forests

表 1 不同環(huán)境因子對我國優(yōu)勢針葉樹種葉片化學(xué)計量特征貢獻(xiàn)率Table 1 Contribution of different environmental factors to leaf stoichiometry characteristics of dominant coniferous species in China %

表 2 中國優(yōu)勢針葉林葉片生態(tài)化學(xué)計量內(nèi)穩(wěn)態(tài)指數(shù)分析Table 2 Homeostasis analysis for the stoichiometry in China’s coniferous forests

表 3 中國針葉林N、P含量及N/P研究對比Table 3 Comparison of the N, P, and N/P among different studies

為研究中國區(qū)域針葉林土壤養(yǎng)分對植物養(yǎng)分的影響，選取3個土層的N、P、K數(shù)據(jù)用于研究植物養(yǎng)分供應(yīng)。除K、P外，土壤N含量以及對應(yīng)化學(xué)計量均隨土壤深度增加而降低，其中，N/P、N/K隨土壤深度的增加顯著降低，而P、K含量呈相對穩(wěn)定的分布特征。由此推斷，隨著中國N沉降的加劇，且伴隨著降雨導(dǎo)致土壤淋溶增加，中國針葉林土壤N/P、N/K有增大趨勢，進(jìn)而使中國針葉林土壤P限制潛勢繼續(xù)增加。

3.2 不同環(huán)境要素對中國優(yōu)勢針葉林優(yōu)勢樹種葉片生態(tài)化學(xué)計量貢獻(xiàn)分析

在本研究中，中國針葉林優(yōu)勢樹種葉片生態(tài)化學(xué)計量特征與環(huán)境因子之間具有顯著的關(guān)系。不同于前人的研究以緯度、溫度及降水之間的關(guān)系為主[7, 9, 11]，本研究中增加了更多的環(huán)境因子，如土壤養(yǎng)分化學(xué)計量特征、年平均溫度、年降雨量、水熱指數(shù)等。本研究發(fā)現(xiàn)，隨著緯度增加，植物葉片中P含量呈顯著增加，這表明我國針葉樹種葉片P含量在空間上呈顯著的緯向分布特征；此外，水熱指數(shù)和土壤能夠解釋12.75%和11.20%的P貢獻(xiàn)率，表明葉片中土壤P含量以及降水、溫度三者對葉片P含量有顯著影響。任書杰等[35]研究東北高緯度針葉林N、P含量隨緯度變化規(guī)律時發(fā)現(xiàn)，在一個相對小的高緯度地區(qū)，隨著緯度增加，葉片N含量會隨著緯度增加反而呈顯著減少的趨勢，而P含量未達(dá)到統(tǒng)計學(xué)顯著變化，這與中國區(qū)域其他功能群植物葉片變化特征不同[9-11]，本研究未發(fā)現(xiàn)葉片中N含量隨著緯度顯著增加，推測產(chǎn)生這種現(xiàn)象的原因，一方面是由于本研究所選取的天然林中的針葉樹相對于闊葉樹生長緩慢，且土壤因子對葉片N的貢獻(xiàn)合計達(dá)17.01%，表明中國優(yōu)勢針葉林葉片N含量主要受到了土壤N含量的影響；另一方面是由于針對興安落葉松的研究尺度較小，其他影響因素對葉片化學(xué)計量的影響超過了緯度的貢獻(xiàn)。本研究發(fā)現(xiàn)，針葉林葉片中K含量隨著緯度增加呈顯著增加，推測這與中國區(qū)域土壤K含量分布的規(guī)律一致，推測原因?yàn)樵谀戏礁邷囟嘤甑牡貐^(qū)，風(fēng)化造成K容易流失[39]，而北方溫度低且降雨少，土壤中K易于積累，不易成為限制元素。因此，中國南北氣候帶差異也造成了植物葉片內(nèi)K含量南低北高的空間分布規(guī)律。

3.3 中國針葉林優(yōu)勢樹種葉片化學(xué)計量穩(wěn)態(tài)性特征

植物養(yǎng)分含量和化學(xué)計量特征與土壤養(yǎng)分供應(yīng)之間存在著顯著關(guān)系[1,38]，在土壤養(yǎng)分環(huán)境發(fā)生變化時，植物可通過內(nèi)穩(wěn)態(tài)調(diào)節(jié)機(jī)制來保持體內(nèi)化學(xué)計量特征穩(wěn)定性。本研究通過分析中國針葉林優(yōu)勢樹種內(nèi)穩(wěn)態(tài)調(diào)節(jié)機(jī)制發(fā)現(xiàn)，我國針葉林優(yōu)勢樹種除葉片N含量易受到土壤N含量變化影響外，其他的葉片生態(tài)化學(xué)計量特征均呈絕對穩(wěn)態(tài)，這表明我國優(yōu)勢針葉樹種具有良好的土壤環(huán)境適應(yīng)特征，但目前人類工業(yè)活動釋放的過量N以及全球溫度增加對土壤中氮素含量變動產(chǎn)生了一定影響[16]。已有研究證明[19]，在模擬大氣氮沉降增加情景下，針葉葉片中N含量增加，能夠顯著影響與之相關(guān)的其他化學(xué)計量特征，因此，需要多加關(guān)注土壤N含量變動對葉片生態(tài)化學(xué)計量特征產(chǎn)生的影響。

4 結(jié)論

中國區(qū)域針葉林葉片N、P、K平均含量分別為13.16、1.02、6.82 g·kg-1，N/P、N/K、K/P平均比分別為15.48、2.17、7.77。緯度、土壤、降水和溫度因子共同對我國針葉林優(yōu)勢樹種葉片化學(xué)計量特征產(chǎn)生影響，緯度對葉片N和K貢獻(xiàn)率最高，水熱指數(shù)對葉片P含量貢獻(xiàn)率最高。對N/P、N/K、P/K影響最大的分別為年降雨量、土壤40～60 cm深度的N/K值以及年降雨量。除N含量外，中國針葉林優(yōu)勢樹種葉片生態(tài)化學(xué)計量特征均呈絕對穩(wěn)態(tài)，這表明我國優(yōu)勢針葉樹種具有良好的環(huán)境適應(yīng)特征。