劉 先,索沛蘅, 杜大俊, 劉雨暉, 胡亞林, 王玉哲,*

1 福建農(nóng)林大學林學院, 福州 350002 2 福建農(nóng)林大學莘口教學林場, 三明 365002

杉木[Cunninghamialanceolate(Lamb.) Hook]是我國南方重要的速生用材樹種,在我國種植面積廣大(1.21×107hm2),約占全國人工林總面積的25%[1],其在經(jīng)濟、生態(tài)價值方面都有不可取代的地位。杉木人工林取代天然林后,由于經(jīng)營方式單一,杉木純林和連栽導致了土壤質(zhì)量退化和生產(chǎn)力下降[2- 4]。研究表明,杉木多代連栽會改變土壤水分、pH、全碳、全氮含量、養(yǎng)分利用有效性以及林地生物地球化學循環(huán)過程[5- 7]。以往對杉木連栽的研究多集中于幼齡林和中齡林,發(fā)現(xiàn)杉木人工林林分生物量隨連栽代數(shù)遞減[8],土壤物理性質(zhì)惡化[9-10]且養(yǎng)分含量逐代下降(如土壤全氮、速效氮、速效鉀等)[11- 12],但對于連栽后近熟林-成熟林的研究還比較少。

土壤酶不僅參與調(diào)控土壤養(yǎng)分循環(huán)、有機質(zhì)的降解和轉(zhuǎn)化,還與土壤結(jié)構和有機質(zhì)形成有重要關系[13- 14]。因為土壤酶對環(huán)境變化的敏感性,所以可將其作為評估土壤系統(tǒng)健康狀況(如承受脅迫或處于恢復階段)的指示指標[15- 16]。研究表明,土壤酶是土壤結(jié)構和肥力的良好指示因子,如Dick等[17]發(fā)現(xiàn)土壤容重和脫氫酶、磷酸酶等幾種酶活性有顯著的負相關關系,Verstraete和Voets[18]發(fā)現(xiàn)土壤磷酸酶、土壤脲酶和β-葡萄糖苷酶與土壤肥力和作物產(chǎn)量顯著相關。在森林生態(tài)系統(tǒng)中,土壤酶的地位尤其重要,在不施肥的條件下森林土壤中的養(yǎng)分主要來源于凋落物的分解,其主要成分包括纖維素、木質(zhì)素、幾丁質(zhì)、蛋白、核酸、磷脂等,而目前研究的酶大都參與這些物質(zhì)的分解過程[19- 20]。土壤酶活性升高有利于凋落物和有機質(zhì)的分解,加快養(yǎng)分歸還的速度,對于提高森林土壤肥力和維持生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動具有重要意義。

森林生態(tài)系統(tǒng)土壤酶活性受多種因子影響,包括土壤因子(土壤pH、含水量、土壤溫度和有機質(zhì)含量等)、氣候因子(溫度和降雨量)和地理因子(林型、經(jīng)緯度以及海拔),這些因子相互作用共同調(diào)控土壤微生物群落結(jié)構,從而影響酶的釋放和活性[21]。常綠闊葉林轉(zhuǎn)換為杉木純林后,隨連栽代數(shù)增加林下微氣候環(huán)境發(fā)生改變,土壤保水能力下降,凋落物分解緩慢造成養(yǎng)分循環(huán)受阻且利用效率下降[22- 23]。大量研究表明連栽會改變杉木林的養(yǎng)分循環(huán)特征[6,23],而土壤微生物是土壤碳、氮循環(huán)過程的主要驅(qū)動力,直接調(diào)節(jié)土壤有機質(zhì)和養(yǎng)分的供給和轉(zhuǎn)化[24- 25]。土壤酶是表征微生物活性的重要指標,β-葡萄糖苷酶和N-乙酰氨基葡萄糖苷酶參與土壤碳循環(huán)過程,對土壤中不同穩(wěn)定性碳的分解和轉(zhuǎn)化起重要作用[26],土壤脲酶、蛋白酶、硝酸還原酶、亞硝酸還原酶和天冬酰胺酶參與土壤氮循環(huán)過程,能夠?qū)⑼寥乐胁煌螒B(tài)N分解為有效態(tài)N供植物吸收,對提高N素利用率有重要意義[27]。目前對多代連栽杉木人工林酶活性的研究主要關注脲酶、蛋白酶、蔗糖酶和過氧化氫酶等,且大部分結(jié)果表明隨杉木栽植代數(shù)增加酶活性逐代下降[28- 30],但是對于參與土壤碳氮循環(huán)過程的酶如β-葡萄糖苷酶、硝酸還原酶、亞硝酸還原酶、天冬酰胺酶則較少涉及。大多數(shù)研究主要在杉木幼齡林和中齡林中進行,缺乏對成熟林和過熟林的長期跟蹤研究。因此,研究不同連栽代數(shù)杉木人工林土壤碳(β-葡萄糖苷酶；N-乙酰氨基葡萄糖核苷酶)、氮循環(huán)(脲酶、蛋白酶、硝酸還原酶、亞硝酸還原酶、天冬酰胺酶)過程酶活性變化及其影響因子,對杉木人工林管理和土壤生態(tài)恢復有積極的指導意義。

1 材料與方法

1.1 試驗地概況與樣地設置

試驗地位于福建省三明市福建農(nóng)林大學莘口教學林場(東經(jīng)117°28′,北緯26°10′),屬武夷山支脈,海拔175—264 m,為低山丘陵地貌。土壤為粉砂巖發(fā)育的山地紅壤,土壤厚度為80—90 cm。該地區(qū)屬中亞熱帶季風氣候,年均氣溫19.1 ℃,年積溫5240.2—6677.9 ℃,年均降水量1749 mm,年平均蒸發(fā)量1585 mm,年平均相對濕度81%,全年日照時數(shù)1840 h,無霜期300 d左右。主要林下植被有狗脊蕨(Woodwardiajaponica(L.f.) Sm.)、江南卷柏(SelaginellamoellendorfiiHieron.)和三葉崖爬藤(TetrastigmahemsleyanumDiels et Gilg)等。天然林研究區(qū)喬木層主要由格氏栲、馬尾松和木荷等構成；草本層主要包括芒萁(Dicranopterispedata(Houtt.) Nakaike)、黑莎草(GahniatristisNees)和華山姜(AlpiniaoblongifoliaHayata)等種類。

在經(jīng)營方式、土壤類型相同以及海拔、坡向、坡度等自然狀況相似的杉木人工林建立12塊固定試驗樣地,每個處理4個重復,樣地面積均為20 m × 20 m。選擇不同連栽代數(shù)的杉木人工林(一代 FRP、二代 SRP、三代 TRP)和格氏栲天然林(NF)作為研究對象,不同連栽代數(shù)杉木人工林樣地的概況見表1。

表1 試驗地概況

FRP：First rotation plantation;SRP：Second rotation plantation;TRP：Third rotation plantation; NF：Natural forest

1.2 土壤樣品采集

1.3 土壤理化性質(zhì)和酶活性分析

測定土壤酶活性的鮮土用PBS緩沖液浸提(水土比1∶10)后,采用上海酶聯(lián)公司Elisa酶聯(lián)免疫分析試劑盒提取,利用多功能酶標儀(BioTek Synergy 2, Gene 5,America)在波長450 nm下測定其OD值,最后用標準物的濃度與OD值計算出標準曲線的回歸方程,將樣品的吸光值(OD)代入方程式,計算出樣品濃度,再乘以稀釋倍數(shù),即為樣品的實際濃度[31]。酶活性結(jié)果以單位質(zhì)量(g-1)干土單位時間(h-1)內(nèi)生成的酶的摩爾質(zhì)量(μmol)來表征。該試劑盒原理是應用雙抗體夾心法測定土壤樣本中的酶活性水平,以亞硝酸還原酶為例,用純化的亞硝酸還原酶抗體包被微孔板,制成固相抗體,往包被單抗的微孔中依次加入亞硝酸還原酶,再與辣根過氧化物酶(HRP)標記的亞硝酸還原酶抗體結(jié)合,形成抗體-抗原-酶標抗體復合物,經(jīng)過徹底洗滌后加底物四甲基聯(lián)苯胺(TMB)顯色,TMB在HRP酶的催化下轉(zhuǎn)化成藍色,并在酸的作用下轉(zhuǎn)化成最終的黃色,顏色的深淺和樣品中的亞硝酸還原酶呈正相關。其他酶均以類似方法測定。

表2 土壤酶名稱、縮寫、國際分類號及功能

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

首先對所有測得數(shù)據(jù)進行正態(tài)分布和方差齊性檢驗,然后采用一般線性模型中的重復測量方差(repeated measures ANOVA)分析連栽代數(shù)、采樣季節(jié)和二者交互作用對土壤理化性質(zhì)和酶活性的影響。每個采樣季節(jié)內(nèi)不同連栽代數(shù)土壤理化性質(zhì)及酶活性變化采用單因素方差分析(one-way ANOVA)進行檢驗,并用LSD法進行多重比較(α=0.05)。土壤酶活性和理化性質(zhì)的相關性采用Pearson相關分析。以上數(shù)據(jù)分析均在SPSS 22.0(SPSS Inc., USA)軟件中進行。為探討土壤酶活性的影響因素,利用Canoco 5.0(Microcomputer Power, Inc., Ithaca, NY)對酶活性和環(huán)境因子(土壤pH、含水量、全碳、全氮、碳氮比、銨態(tài)氮、硝態(tài)氮)進行冗余分析(RDA)。

2 結(jié)果與分析

2.1 連栽杉木林土壤理化性質(zhì)

表3 不同連栽代數(shù)和季節(jié)變化對土壤理化性質(zhì)的重復方差測量分析結(jié)果

表中數(shù)值是重復測量方差的P值,*P< 0.05,**P< 0.01,***P< 0.001

表4 不同連栽代數(shù)杉木人工林和天然林土壤理化性質(zhì)比較

Table 4 Comparisons of soil physiochemical properties in different generations of continuously-cropped Chinese fir plantation and natural forest

采樣季節(jié)Sampling season連栽代數(shù)Planting generationpH含水量/%Moisture全碳 Total carbon/(g/kg)全氮Total nitrogen/(g/kg)銨態(tài)氮Ammonium/(mg/kg)硝態(tài)氮Nitrate/(mg/kg)冬季FRP4.00±0.03c21.43±0.45b22.45±2.18ab1.54±0.09ab5.81±0.28b4.96±1.00aWinterSRP4.19±0.03ab27.46±1.23a25.58±1.79a1.77±0.05a6.38±0.68b4.58±0.48aTRP4.07±0.03bc20.51±0.13bc27.15±1.66a1.75±0.12a7.97±1.00b4.81±0.44aNF4.30±0.02a17.06±0.65c15.83±1.60b1.13±0.15b30.69±5.05a0.76±0.13b夏季FRP4.23±0.02b29.05±0.28a21.72±0.91b1.54±0.04a5.70±2.59a13.49±1.00aSummerSRP4.26±0.03b32.78±2.04a26.00±1.78ab1.79±0.09a7.40±0.18a2.97±0.57cTRP4.14±0.04b31.06±0.64a29.82±2.22a1.88±0.15a3.16±0.35a10.58±0.58abNF4.42±0.05a28.41±1.62a28.60±2.25ab1.81±0.07a6.65±1.47a8.61±1.40b

FRP,連栽第一代杉木林；SRP,連栽第二代杉木林；TRP,連栽第三代杉木林；NF,格氏栲天然闊葉林;表中數(shù)值為平均值±標準誤(n=4),不同字母代表同一季節(jié)不同連栽代數(shù)之間存在顯著性差異(P< 0.05)

2.2 參與土壤碳氮轉(zhuǎn)化過程酶活性

不同連栽代數(shù)對參與土壤碳氮循環(huán)酶活性無顯著影響(表5,P> 0.05),不同采樣季節(jié)對除亞硝酸還原酶之外的其他6種土壤酶活性均有顯著影響(表5,P< 0.01),而采樣季節(jié)和連栽代數(shù)的交互作用對土壤天冬酰胺酶和亞硝酸還原酶活性有顯著影響(表5,P≤ 0.05)。天然林和不同連栽代數(shù)杉木林土壤酶的測定結(jié)果表明,無論12月還是6月,連栽對大多數(shù)土壤酶的活性并無顯著影響(圖1,蛋白酶、β-葡萄糖苷酶、硝酸還原酶、亞硝酸還原酶、N-乙酰氨基葡萄糖苷酶)。 只有在12月,天然林的土壤脲酶活性(0.018—0.025 μmol h-1g-1)顯著高于連栽杉木人工林(0.005—0.014 μmol h-1g-1),約為后者的2—5倍,而不同代數(shù)連栽杉木林土壤脲酶活性并無顯著差異；天然林土壤天冬酰胺酶活性(22.09—30.92 μmol h-1g-1)顯著高于三代連栽杉木林(10.84—25.61 μmol h-1g-1),與一代、二代連栽杉木林并無顯著差異；在6月,4種林地的土壤脲酶和天冬酰胺酶活性均無顯著差異(圖1)。

圖1 不同連栽代數(shù)杉木人工林與天然林土壤酶活性冬夏兩季變化Fig.1 Changes of soil enzyme activities in different rotations of Chinese fir plantation and natural forest in December and June不同小寫字母表示同一季節(jié)土壤酶活性在不同連栽代數(shù)間存在顯著差異(P < 0.05)

表5 不同連栽代數(shù)和季節(jié)變化對土壤酶活性的重復方差測量分析結(jié)果

Table 5 Repeated measures ANOVA on the effects of planting generations and seasonal changes on soil enzyme activities

來源Source蛋白酶Proteaseβ-葡萄糖苷酶β-Glucosidase脲酶Urease天冬酰胺酶ASPG硝酸還原酶NAR亞硝酸還原酶NIRN-乙酰葡萄糖核苷酶NAG季節(jié) Time (T)<0.001???<0.001???<0.001???<0.001???0.008??0.253<0.001???栽植代數(shù) Generation (G)0.5130.556<0.001???0.4660.4620.5560.241季節(jié) × 栽植代數(shù) T × G0.0650.6640.1240.017?0.3340.036?0.075

表中數(shù)值是重復測量方差的P值,*P< 0.05,**P< 0.01,***P< 0.001

2.3 土壤酶活性與土壤理化性質(zhì)的相關性

表6 土壤酶活性與環(huán)境因子的Pearson相關分析(n=32)

表中數(shù)值是皮爾遜相關系數(shù)P值,所用數(shù)據(jù)為12月和6月兩次采樣的數(shù)據(jù);*P< 0.5;**P< 0.01;***P< 0.001

圖2 土壤酶活性與環(huán)境因子相關性的典范冗余分析(RDA),分析所用數(shù)據(jù)為12月和6月兩次采樣的數(shù)據(jù)Fig.2 Redundancy analysis (RDA) of relationships between soil enzyme activities and the environmental factors, using the bulked December and June dataASPG, Asparaginase， 天冬酰胺酶；BG，β-1,4-glucosidase, β-葡萄糖苷酶；NAG，N-acetyl-β-D-glucosaminidase，幾丁質(zhì)酶；urease，脲酶；protease，蛋白酶；NAR，Nitrate reductase，硝酸還原酶；NIR，Nitrite reductase，亞硝酸還原酶。TC，total carbon， 全碳；TN，total nitrogen， 全氮；C/N ratio，土壤碳氮比，moisture，土壤含水量

3 討論

本研究中除脲酶活性和天冬酰胺酶活性在天然林和杉木人工林中有顯著差異(圖1),連栽并未改變參與土壤碳、氮循環(huán)的幾種主要酶活性,這與前人報道中杉木連栽顯著降低土壤脲酶活性的結(jié)果并不一致[30, 43- 44]。這可能與杉木所處的發(fā)育階段不同有關,杉木林進入成熟林-過熟林階段后很多土壤理化指標已經(jīng)恢復,間伐等整地措施使林下植被環(huán)境變好,微生物活性增強,因此土壤狀況有所恢復[45]。前人研究表明,中度間伐(保留1200株/hm2)能顯著提高土壤過氧化氫酶、堿性磷酸酶、脲酶和蔗糖酶的活性[46],本研究中杉木人工林造林密度在2700—3000株/hm2,間伐后密度為1050—1845株/hm2,介于輕度間伐與中度間伐之間。間伐后,林下植被多樣性增加,促進土壤微生物種類和數(shù)量增加,加快土壤養(yǎng)分循環(huán)速度,另一方面使林分郁閉度降低,透光性增強,從而提高土壤酶活性[46]。本研究所涉及的參與土壤碳氮轉(zhuǎn)化過程酶活性在其他杉木連栽林中報道較少,雖然蛋白酶和β-葡萄糖苷酶在4種林地間無顯著差異,季節(jié)差異卻很明顯,這可能與微生物活躍程度或植物養(yǎng)分需求的季節(jié)性差異有關。

4 結(jié)論

由以上分析可知,杉木人工林進入成熟林-過熟林階段后,不同連栽代數(shù)林地許多土壤理化指標及酶活性均有不同程度恢復。2代和3代連栽杉木人工林的養(yǎng)分含量并不比1代杉木林低,有的(如全氮和硝態(tài)氮)甚至超過了天然林。參與土壤碳氮轉(zhuǎn)化過程酶活性并未呈現(xiàn)出隨連栽下降的趨勢,季節(jié)性變化大于連栽影響。土壤pH、含水量和養(yǎng)分含量是決定土壤碳氮循環(huán)過程酶活性的重要因子,但是由于不同酶對土壤條件響應機制的復雜性,想要精確建立土壤酶活性與土壤理化性質(zhì)和養(yǎng)分有效性的關系,還需進一步長期研究。如果能建立這種關系,對于預測森林生態(tài)系統(tǒng)土壤質(zhì)量的改變,探究人工林土壤退化及恢復過程中的酶學機理,有著重要的生態(tài)學意義。