固沙林恢复土壤酶化学计量特征与养分限制效应.pdf

第3 0卷第6期2 0 2 3年1 2月水土保持研究R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o nV o l.3 0,N o.6D e c.,2 0 2 3 收稿日期:2 0 2 2-1 0-2 8 修回日期:2 0 2 2-1 1-2 9 资助项目:陕西省自然科学基金重点项目(2 0 2 3-J C-Z D-1 0);陕西省农业专项资金项目NYK J-2 0 2 2-Y L(X N)1 6;国家大学生创新创业训练计划项目(2 0 2 2 1 0 7 1 2 1 5 8)第一作者:张延文(1 9 9 7),女,陕西咸阳人,硕士研究生,主要研究方向:人工林恢复碳循环。E-m a i l:1 5 9 9 1 0 8 2 4 8 41 6 3.c o m 通信作者:佟小刚(1 9 7 9),男(锡伯族),新疆伊犁人,副教授,博士,主要从事人工林恢复生态效应与机制研究。E-m a i l:x i a o g a n g t o n g 1 2 6.c o mh t t p:s t b c y j.p a p e r o n c e.o r gD O I:1 0.1 3 8 6 9/j.c n k i.r s w c.2 0 2 3.0 6.0 5 0.张延文,周雅洁,张不染,等.固沙林恢复土壤酶化学计量特征与养分限制效应J.水土保持研究,2 0 2 3,3 0(6):1 0 2-1 1 1.Z HAN GY a n w e n,Z HOUY a j i e,Z HAN GB u r a n,e t a l.S o l i E n z y m a t i cS t o i c h i o m e t r ya n dN u t r i e n tL i m i t a t i o nD u r i n gR e v e g e t a t i o n i nS a n d-F i x i n gF o r e s tJ.R e s e a r c ho fS o i l a n dW a t e rC o n s e r v a t i o n,2 0 2 3,3 0(6):1 0 2-1 1 1.固沙林恢复土壤酶化学计量特征与养分限制效应张延文1,周雅洁2,张不染2,拓卫卫1,杨 京1,佟小刚2(1.西北农林科技大学 水土保持研究所,陕西 杨凌7 1 2 1 0 0;2.西北农林科技大学 资源与环境学院,陕西 杨凌7 1 2 1 0 0)摘 要:目的 探究土壤胞外酶及酶化学计量比对沙漠化土地植被恢复过程的响应特征及其养分限制效应,可深入揭示人工固沙林重建植被过程养分转化能力与机理。方法 基于时空替代法,以榆林毛乌素沙区的半固定沙地为(0a)对照、选择恢复2 55 6a灌木与乔木固沙林时间序列样地,测定分析了01 0c m和1 02 0c m土层的胞外酶活性、酶生态化学计量的演变特征及其限制微生物养分需求效应。结果 随着植被恢复年限增加,两种林地土壤-葡萄糖苷酶(B G)、-1,4-乙酰氨基葡萄糖苷酶(NAG)、亮氨酸氨肽酶(L A P)和碱性磷酸酶(A P)活性在不同土层均呈显著升高趋势。但酶活性的不均衡变化也显著改变了土壤酶化学计量特征,5 6a植被恢复过程中,土壤酶CN,CP,NP呈增加趋势。土壤胞外酶化学计量比也表明植被恢复期间微生物生长受N限制显著增加,P限制并不明显,仅在乔木林恢复2 5a时1 02 0c m土层中出现P限制。结论 固沙林植被恢复显著提升了土壤酶活性,增强了土壤碳氮代谢的能力,且在恢复过程中土壤微生物活动主要受到N养分的限制。关键词:胞外酶活性;生态酶化学计量特征;养分限制;固沙林;毛乌素沙地中图分类号:S 7 1 4 文献标识码:A 文章编号:1 0 0 5-3 4 0 9(2 0 2 3)0 6-0 1 0 2-1 0S o l iE n z y m a t i cS t o i c h i o m e t r ya n dN u t r i e n tL i m i t a t i o nD u r i n gR e v e g e t a t i o ni nS a n d-F i x i n gF o r e s tZ h a n gY a n w e n1,Z h o uY a j i e2,Z h a n gB u r a n2,T u oW e i w e i1,Y a n gJ i n g1,T o n gX i a o g a n g2(1.I n s t i t u t eo fS o i la n dW a t e rC o n s e r v a t i o n,N o r t h w e s tA&FU n i v e r s i t y,Y a n g l i n g,S h a a n x i7 1 2 1 0 0,C h i n a;2.C o l l e g eo fN a t u r a lR e s o u r c e sa n dE n v i r o n m e n t,N o r t h w e s tA&FU n i v e r s i t y,Y a n g l i n g,S h a a n x i7 1 2 1 0 0,C h i n a)A b s t r a c t:O b j e c t i v eT h i ss t u d yw a sc o n d u c t e dt oe x p l o r et h er e s p o n s eo fs o i le x t r a c e l l u l a re n z y m e sa n de n z y m es t o i c h i o m e t r i cr a t i o st ot h er e v e g e t a t i o np r o c e s so fd e s e r t i f i e dl a n da n dt h e i rn u t r i e n tl i m i t a t i o ne f f e c t s,w h i c hc a np r o v i d ei n s i g h t i n t ot h en u t r i e n tc o n v e r s i o nc a p a c i t ya n dm e c h a n i s mo ft h er e v e g e t a t i o np r o c e s so f a r t i f i c i a l s a n d-f i x a t i o nf o r e s t.M e t h o d sB a s e do nt i m e-s p a c es u b s t i t u t i o n,t h es e m i-f i x e ds a n d yl a n dw a ss a m p l e da s 0y e a r s c o n t r o l s i t e,a n d t h e s h r u ba n da r b o r s a n d-f i x i n g f o r e s t sw i t hc h r o n o s e q u e n c e s o f2 55 6y e a r sr e v e g e t a t i o ni n M u U sS a n d yL a n do fY u l i n w e r ec h o s e n.T h ee v o l u t i o nc h a r a c t e r i s t i c so fe n z y m ea c t i v i t ya n d i t ss t o i c h i o m e t r y,a n dr e s o u r c e l i m i t a t i o no nm i c r o b i a ln u t r i e n t r e q u i r e m e n t i nb o t h01 0c ma n d1 02 0c ms o i l l a y e r sw e r e a n a l y z e d.R e s u l t sT h e a c t i v i t yo f-g l u c o s i d a s e(B G),-1,4-a c e t a m i-n o g l u c o s i d a s e(NAG),l e u c i n ea m i n o p e p t i d a s e(L A P)a n da l k a l i n ep h o s p h a t a s e(A P)i nb o t hs o i ll a y e r su n d e r t w of o r e s t l a n d s i n c r e a s e dc o n t i n u o u s l yw i t hy e a r i n c r e a s eo fv e g e t a t i o nr e s t o r a t i o n.H o w e v e r,u n b a l-a n c e dc h a n g e sa m o n gt h ee n z y m ea c t i v i t i e sr e s u l t e di ns i g n i f i c a n t i n c r e a s ei ns o l ie n z y m a t i cs t o i c h i o m e t r y.T h er a t i o so fCN,CPa n dNPc a l c u l a t e db a s e do ne n z y m ea c t i v i t ys h o w e dt h e i n c r e a s i n gt r e n dd u r i n g5 6y e a r so f r e v e g e t a t i o n i nb o t hf o r e s t l a n d s.A n d,s i g n i f i c a n t i n c r e a s e i ns o l i e n z y m a t i c s t o i c h i o m e t r yr e f l e c t e dt h a tm i c r o b i a l g r o w t hw a sl i m i t e db yn i t r o g e n(N)r a t h e rt h a np h o s p h o r u s(P).R e s o u r c el i m i t a t i o nb yPo n l yw a so b s e r v e d i n1 02 0c ms o i l l a y e r i na r b o r l a n dw i t h2 5y e a r so f r e v e g e t a t i o n.C o n c l u s i o nR e v e g e t a t i o n i ns a n d-f i x i n g f o r e s t i m p r o v e d t h e e n z y m e a c t i v i t i e s a n d e n h a n c e d t r a n s f o r m a t i o no f c a r b o na n dn i t r o g e n i n s o i l,a n d s o i lm i c r o b i a l a c t i v i t yw a sm a i n l y l i m i t e db yNn u t r i e n t d u r i n g t h ev e g e t a t i o nr e s t o r a t i o np r o c e s s.K e y w o r d s:e x t r a c e l l u l a re n z y m ea c t i v i t y;e c o l o g i c a l e n z y m es t o i c h i o m e t r i c;n u t r i e n t l i m i t a t i o n;s a n d-f i x a t i o nf o r e s t;M uU ss a n d y l a n d 沙漠化是全球生态环境防治的重大问题。过去七十年,我国通过大规模生态林工程的实施,实现了沙漠化的逆转,产生了巨大的生态恢复效益1-2。但限于沙漠化地区的脆弱生态环境,维持固沙林植被与土壤的自我恢复能力,成为沙区生态系统可持续恢复的关键3。大量研究表明林地系统养分转化与累积调控着植被与土壤的协同恢复4。以往人们多关注于固沙林植被恢复的生态功能,如土壤养分含量的增加、微生物活性的提高等5-7。现在越来越多的研究关注到植被和土壤之间的协同交互关系以及相关的生物学养分转化过程,这对认识植被恢复养分循环机理具有重要的科学意义。土壤胞外酶作为土壤有机质分解的驱动力,可以调节有机质分解转化和微生物吸收利用土壤C,N,P等养分元素8-9,有效催化生态系统中物质与能量的转化1 0。因此,探究生态酶活性可以表征微生物对能量与养分的利用效率1 1。特别是生态酶的化学计量特征与微生物的养分需求密切相关1 2,微生物活动受何种养分资源限制可以由土壤生态酶计量比表达1 3-1 4。如通过酶活性计算的矢量长度和矢量角度可反映微生物活动受限养分类型,也可帮助明确植被恢复过程中微生物营养限制的演变特征。S i n s a b a u g h等1 5研究发现酶的生态化学计量比可以更好地反映微生物的代谢和营养需求之间的 生化平衡。Z h a o等1 6将土壤碳、氮、磷的化学计量学与土壤微生物和-葡萄糖苷酶、亮氨酸氨肽酶、-1,4-乙酰氨基葡萄糖苷酶、碱性磷酸酶的化学计量学结合起来研究,结果表明土壤微生物和酶的化学计量学在很大程度上受土壤养分化学计量学的调控。因此,土壤胞外酶及其生态化学计量通过反映土壤微生物对碳氮磷养分的需求,从而可以显示出微生物的营养限制作用1 7。但目前生态酶化学计量特征及其养分限制效应对沙漠化土地植被恢复的响应研究还较为薄弱。毛乌素沙地是我国沙漠化治理的典型区域,通过6 0多年的治理,植被覆盖率已由5 0年代连片沙丘的1.8%1 8增加到目前的8 0%。鉴于此,本研究选择毛乌素沙地从半固定沙地到长期恢复的典型灌木和乔木固沙林为研究对象,探究沙漠化土地植被恢复过程中土壤生态酶活性及其化学计量比演变特征以及对微生物养分限制的影响,深入揭示人工固沙林重建植被过程养分转化能力与机理。1 材料与方法1.1 研究区概况项目研究区位于毛乌素沙地东缘陕西榆林市红石峡试验林场(3 8 2 0 N,1 0 9 4 4 E)(图1)。区域海拔10 5 015 0 0m,属典型的大陆性季风干旱和半干旱气候,平均温度为8.11 9,年平均降水量4 1 5mm,其中6 0%7 0%的降雨量集中在夏季。由于该地区是典型的农牧交错带,历史上经过过度的开垦、耕作和放牧导致了严重的沙漠化,风成砂质土中含沙量为9 2.6%2 0。在此之后该区域陆续开展了防风治沙的生态工程。先在流沙表面设置草方格机械沙障固定流沙,后有天然先锋植物侵入并将流沙变为半固沙地2 1。经过持续人工管理,每隔5a持续种植相同树种、相似密度的固沙灌木和乔木,逐渐形成了灌木草或乔木草的植被体系。目前该区域完整的植被面积扩增至3 3.3k m2,植被逐渐恢复,形成以樟子松(P i n u ss y l v e s t r i sv a r.m o n g h o l i c aL i t v.)为主要树种的乔木防护林2 2,以沙棘(H i p p o p h a er h a m n o i d e sL.)、花棒(H e d y s a r u ms c o p a r i u mF i s c h.e tB a s i n e r)和踏郎(H e d y s a r u m m o n g o l i c u mT u r e z)为主要树种的灌木防 护 林,林 下 形 成 狗 尾 草(S e t a r i av i r i d i sL.B e a u v.)、沙竹(P s a mm o c h l o av i l l o s aT r i n.B o r)、猫头刺(O x y t r o p i s a c i p h y l l aL e d e b.)、油蒿(A r t e m i s i ao r d o s i c aK r a s c h.)等植被,随后植被系统进入自然演替阶段。1.2 样地选取2 0 2 0年9月采用空间替代时间法,基于林场造林年限与树种情况的调研,选择半固定沙地为0a植被恢复对照样地,以及1 9 6 4年、1 9 7 4年、1 9 8 5年和301第6期 张延文等:固沙林恢复土壤酶化学计量特征与养分限制效应1 9 9 5年种植的柠条为主灌木林和以樟子松为主的乔木林,分别对应植被恢复5 6a,4 6a,3 5a和2 5a的林地。每个林地选择立地条件相似的3个样地作为研究重复,样地间的距离大于5 0m,在每个样地上设置2 0m2 0m的标准采样区,采用等距抽样法调查植物群落。所有样地都来自于流沙地,具有同样的原始条件。各样地分布及基本特征见表1,其中S F表示半固定沙地、S H表示灌木固沙林、A R表示乔木固沙林,下角标数字代表恢复年限。1.3 采样方法采用土钻(直径5c m 1 2 0c m长)在每个标准采样区上按照“S”型布 设1 6个点采集土样,混合成该小区代表土样,并同时采用环刀法测量土壤容重。所有土芯均被剖分成两层(01 0c m和1 02 0c m),来自同一层的土壤被混合为一个土壤样本。土壤样本通过2mm的筛,已清除植物残体与根系,同时将鲜土样同冰袋保存带回实验室冷藏,及时完成指标的测定。1.4 指标测定1.4.1 土壤理化性质的测定 土壤含水量是在1 0 5 的烘箱内加热2 4h后,采用称重的方法测定。土壤容重采用环刀法,土壤p H使用带有玻璃电极的p H计(P H S-3 C,上海仪电,中国)在水土比为12.5(干土质量:去离子水体积)下测得。图1 研究区域在毛乌素沙地中的位置F i g.1 L o c a t i o no f t h e s t u d ya r e a i nt h eM uU sS a n d yL a n d表1 样地基本特征及0-2 0c m土层深度的土壤理化性质T a b l e1 D e s c r i p t i o n so f t h e s a m p l i n gs i t e sa n ds o i l p h y s i c o c h e m i c a l p r o p e r t i e sa t 0-2 0c ml a y e r样地植被密度/(株h m-2)郁闭度/%含水率/%p H容重/(gc m-3)主要植被类型S F 0 a1 1 0 11 935.0 78.1 21.6 70.0 2 a籽蒿、沙竹S H 2 5 a1 1 3 57 01 05.0 27.9 61.6 10.0 2 a b柠条、油蒿、沙竹S H 3 5 a1 1 1 38 11 24.7 78.1 31.5 40.1 1 b c油蒿、柠条、沙竹、踏郎S H 4 6 a1 0 9 38 786.8 77.81.5 60.0 6 a b c柠条、油蒿、花棒、沙竹、踏郎、猫头刺S H 5 6 a1 1 1 89 01 25.7 17.9 81.5 30.1 5 b c柠条、油蒿、花棒、沙竹、踏郎、猫头刺A R 2 5 a1 0 7 76 795.1 47.8 81.4 50.0 5 c d乔木、狗尾草、紫穗槐A R 3 5 a1 1 1 38 01 45.6 17.7 71.5 40.0 7 b c乔木、狗尾草A R 4 6 a1 1 3 88 51 15.5 58.11.4 20.0 2 d乔木、狗尾草A R 5 6 a1 1 1 29 31 67.3 67.81.5 30.0 4 b c乔木注:每列中不同小写字母代表不同样地间指标在0.0 5水平下差异显著。1.4.2 酶活性测定 采用改进的微孔板荧光法9测定酶活性(表2)。具体为:称取1g鲜土置于2 5 0m l三角瓶中,加入1 2 0m l纯水和5m lp H为8.0的醋酸钠缓冲液,震荡2h。再用移液枪分别向酶标版的样品孔加入1 5 0l震荡悬浊液和5 0l酶底物;空白微孔加入1 5 0l震荡悬浊液+5 0l缓冲液;淬火标准孔加入1 5 0l震荡悬浊液和5 0l标准物质;阴性对照孔加入1 5 0l缓冲液和5 0l酶底物(合成的富含C-,N-或P-的底物);参考标准孔加入1 5 0l缓冲液和5 0l标准物质(标准物质为4-甲基伞形酮(4-M e t h y l u m b e l l i f e r o n e,MU B)。最后将酶标板放在2 5 下暗培养4h后加入1m o l/L的N a O H5 0l终止反应。使用多功能酶标仪在3 6 5n m激发光、4 5 0n m波长处读取吸光值1 1。酶活性表达为每克干土单位小时产生的底物物质的量n m o l/(gh)。土壤胞外酶活性计算公式9:Ab=F V/(e V1t m)(1)F=(f-fb)/q-fs(2)e=fr(CsV2)(3)q=(fq-fb)/fr(4)式中:Ab为土壤样品的酶活性n m o l/(gh);F为校正后的样品荧光值;V为土壤悬浊液的总体积(1 2 5401 水 土 保 持 研 究 第3 0卷m l);V1为微孔板每孔中加入样品悬浊液的体积(0.2m l);t为暗培养时间;m为干土样的质量(1g鲜土换算成干土样的质量);f为酶标仪读取样品微孔的荧光值;fb为空白微孔的荧光值;q为淬火系数;fs为阴性对照微孔的荧光值;e为荧光释放系数;fr为参考标准微孔的荧光值;Cs为参照标准微孔的浓度(1 0m o l/L);V2为加入参照标准物的体积(0.0 0 00 5L);fq为淬火标准微孔的荧光值。表2 土壤胞外酶种类T a b l e2 T y p e so f s o i l e x t r a c e l l u l a r e n z y m e s酶缩写底物E.C.-葡萄糖苷酶B G4-甲基伞形酮-D-葡萄糖苷3.2.1.7 4-N-乙酰氨基葡萄糖苷酶NA G4-甲基伞形酮-2-乙酰氨基-2-脱氧-D-吡喃葡萄糖苷3.2.1.5 2亮氨酸氨肽酶L A P1,2-二软酯酰基蛋黄素3.4.1 1.1碱性磷酸酶A P4-甲基伞形酮磷酸酯3.1.3.2注:E.C.为国际酶委员会编号。土壤胞外酶生态化学计量计算公式9:CNE E A=B G/(NAG+L A P)(5)CPE E A=B G/A P(6)NPE E A=(NAG+L A P)/A P(7)式中:E E A为土壤胞外酶活性;B G为-葡萄糖苷酶;NAG为N-乙酰氨基葡萄糖苷酶;L A P为亮氨酸氨肽酶;A P为磷酸酶。土壤胞外酶(L A P+N A G)或N A G的活性变化都可以反映土壤微生物对氮资源的需求状况,且B G(N A G+L A P)A P或B GN A GA P均可表示土壤胞外酶生态化学计量关系1 1,因此本试验测定了L A P和N A G的活性,并用(L A P+N A G)进行计算和分析。土壤微生物养分限制计算公式2 3:V e c t o r l e n g t h=l n B Gl n NAG+L A P()2+l n B Gl n A P2(8)V e c t o ra n g l e=D e g r e e sAT AN 2(l nB G/l nA P),l nB G/l n(NAG+L A P)(9)式中:V e c t o r l e n g t h为矢量长度;V e c t o ra n g l e为矢量角度;B G为-葡萄糖苷酶;NAG为-乙酰氨基葡萄糖苷酶;L A P为亮氨酸氨肽酶;A P为磷酸酶。矢量长度越长,C限制越大;矢量角度小于4 5 时,值越小表示N限制越大,而大于4 5 时,值越大表示P限制越大2 4-2 5。另一种表示方法中,评估微生物资源限制基于酶活性的比率,特别是B G/(L A P+N A G)较高为N限制,(N A G+L A P)/A P表示N和P限制2 6。1.5 数据统计分析数据结果用均值标准差(S D)来表示。采用方差分析(A N O V A)方法,分析年代序列对土壤酶活性差异的影响;比较最小显著差异(L S D)在p 0.0 5时的差异;使用皮尔森相关分析土壤养分特征与酶活性以及化学计量的相关关系;以上所有统计方法是用S P S S2 5.0版本完成,使用O r i g i n2 4.0绘制相关图。2 结果与分析2.1 土壤酶活性变化特征植被恢复后,土壤中-葡萄糖苷酶、-1,4-乙酰氨基葡萄糖苷酶、亮氨酸氨肽酶和碱性磷酸酶活性均随着恢复年限的增长呈显著增加趋势。在01 0c m土层,B G,NAG和L A P酶活性乔、灌木林随着植被年限的增加快速增加,在恢复5 6a后,B G活性分别是对照(0a)的1 0.0 8,1 4.9 9倍,NAG活性分别是对照(0a)的5.7 1,7.0 8倍,L A P活性分别是对照(0a)的8.3 8,6.9 7倍;A P活性在乔、灌木林随着植被年限的增加缓慢增加,且在恢复5 6a后,A P活性分别是对照(0a)的3.7 9,3.8 5倍。在1 02 0c m土层中,5 6a的B G活性乔、灌木林平缓增加分别比对照(0a)增加了3.9 4,1 6.5 8n m o l/(gh);L A P和A P活性在乔、灌木林随着植被年限的增加缓慢增加,L A P活性分别是对照(0a)的1 5.1,1 3.5倍,A P活性分别是对照(0a)的1 9.6,1 7.6倍;NAG乔木林随着植被年限的增加急剧增加是对照(0a)的5 8.1 3倍(图2)。B G活性的乔、灌木林地和NAG+L A P活性的乔木林地在两土层中有显著性差异;01 0c m和1 02 0c m的A P活性在乔木和灌木林地中存在显著差异(图3)。2.2 固沙林恢复过程酶活性比值演变特征植被重建后,灌木林地上,土壤酶CN,NP和CP在不同土层整体上随年限增加而增加;乔木林地上,1 02 0c m土层土酶CN和CP逐渐减少且小于对照(0a)。在01 0c m土层的乔、灌木林中,土壤酶CN,CP和NP逐渐增加均在第5 6501第6期 张延文等:固沙林恢复土壤酶化学计量特征与养分限制效应年达到最大,分别是对照(0a)的1.2倍和2.2倍,2.8倍和4.0倍,2.0倍和1.8倍;且在1 02 0c m土层的灌木林中,土壤酶CN,CP和NP也逐渐增加,均在第5 6年达到最大0.0 8,0.2 5,3.2 9分别是对照(0a)的1.2,3.2倍和2.7倍;在1 02 0c m土层的乔木林中,土壤酶CN和CP平缓减少,在恢复5 6a后分别比对照(0a)减少了0.0 3,0.0 1n m o l/(gh)。经过5 6a的植被恢复,灌木林地土壤酶CN和CP高于对应土层的乔木林地,两种林地土壤酶NP变化相似(图4)。图2 不同恢复年限下乔、灌木林地在不同土层的酶活性特征F i g.2 E n z y m ea c t i v i t yc h a r a c t e r i s t i c so f t r e e sa n ds h r u b s i nd i f f e r e n t s o i l l a y e r su n d e rd i f f e r e n t r e c o v e r yy e a r s2.3 酶化学计量指示微生物营养限制的作用随着植被恢复两类林地矢量长度均逐渐增加,随矢量长度越长,C限制增大,对比对照组在第5 6年达到最大值。乔、灌木林在01 0c m土层中矢量长度达到最大时分别是对照(0a)的1.8倍和2.0倍。在灌木林1 02 0c m土层中,矢量长度波动增加,第5 6601 水 土 保 持 研 究 第3 0卷年时最大为恢复6a的2.8倍(图5左);1 02 0c m土层乔木林中,矢量长度波动增加,在恢复3 6a时最大为对照组(0a)的1.9倍(图5左)。当矢量角度大于4 5 时为P限制,且随着其值的增加而增大,当矢量角度小于4 5 时为N限制,且随着其值的增加而减小。在植被恢复过程中,微生物N限制随着植被恢复显著增加,而P限制并不明显。在01 0c m土层灌木林矢量角度均小于4 5 且无显著变化,乔木林矢量角度逐渐减小且均小于4 5,在第5 6年达到最小;在1 02 0 c m土层灌木林中,矢量角度先增加后减小,且均小于4 5,即产生P限制,乔木林在植被恢复2 5 a时大于4 5,即无N限制产生P限制,之后在第5 6年时最小(图5右)。图3 不同林地与土层酶活性和酶化学计量特征F i g.3 E n z y m ea c t i v i t ya n ds t o i c h i o m e t r i cc h a r a c t e r i s t i c so fd i f f e r e n tw o o d l a n d sa n ds o i l l a y e r s 当(NAG+L A P)/A P大于1时,产生微生物N限制,小于1时产生微生物P限制,当B G/(NAG+L A P)大于1时产生微生物C和N限制以及C和P限制。生态酶化学计量比分析表明,在1 02 0c m土层中,乔木林仅在第2 5年存在微生物P限制,其他年份均为N限制,整个植被恢复阶段均无C和N限制以及C和P限制(图6)。3 讨 论3.1 固沙林恢复过程土壤胞外酶活性及其化学计量演变特征本研究两种固沙林恢复过程酶活性均呈现显著增加的趋势,其中以氮转化酶(NAG+L A P)活性对植被恢复变化响应最敏感。随着植被恢复,固沙林土壤有机碳、全氮和全磷含量不断增加,土壤性状、各理化指标显著提高2 7,引起微生物需求更多养分吸收和转化反应,分泌更多的酶,导致酶活性表现出增强趋势。其中碳转化酶(B G)活性和氮转化酶(NAG+L A P)活性为:整体上随年限增加而增加,磷转化酶(A P)活性为先增加后减小。恢复5 6a后,在不同土壤剖面中3类土壤酶活性都为表层(01 0c m)较高。这与李陆平等2 8的研究结果相同。由于植被恢复首先改善了表层土壤理化性质,使得表层土壤胞外酶活性也随之提高。同时,土壤表层水热条件和通气状况良好,又有微生物生长所需的营养物质,则微生物生长旺盛,代谢活跃使得表层的土壤酶活性较高2 9。701第6期 张延文等:固沙林恢复土壤酶化学计量特征与养分限制效应植被恢复4 4a后氮转化酶(NAG+L A P)的活性在01 0c m土层呈现逐渐稳定的趋势,在1 02 0c m的土层中继续增加,导致产生这种现象的原因可能是由于凋落物主要堆积在表层土壤,从而利于表层土壤肥力增加,深层土壤中酶含量具有滞后效应,当上层土壤中酶的含量达到平衡时,经过上层的淋溶与沉积作用,下层土壤的酶活性仍然继续上升。化学计量比代表了微生物需求的相对比例,当微生物需求的元素少时就会不断分泌相应的酶,以满足自身的需要,从而使化学计量比发生变化。酶的不同比例的增加也显著改变了酶的化学计量。不同土壤胞外酶活性的比值与土壤CNP值密切相关,可以用来评价土壤微生物C,N,P养分资源的需求状况3 0。本研究中,随着植被恢复年限增长,土壤CN,CP和NP逐渐增加,全球尺度上,土壤酶CN比为1.4 1,酶CP比为0.6 2,酶NP比为0.4 4。本研究中,灌、乔木林的酶CN最大值分别为0.3 0,0.1 9,远小于1.4 1的全球平均值;酶NP比的最小值分别为1.0 2,0.8 8,也显著高于全球均值(0.4 4),这意味着该区域土壤整体具有相对较高的N分解胞外酶活性,从而反映出土壤N元素的相对缺乏,使得微生物需要生产更多的N分解酶,以满足自身代谢的需要9。这是基于目前公认的个别碳氮转化酶之间的关系得出的生态酶化学计量关系,还不能完整反映对微生物养分需求的指示作用。可能需要考虑更多的与碳氮转化相关的酶进行化学计量分析。本研究中的灌木林地树种为豆科固氮树种,樟子松、油松为非固氮树种,然而随着植被恢复年代增加,土壤C,N限制仍然呈现增加趋势,建议建立更稳定的乔灌草相结合的植被体系3 1。图4 不同恢复年限下乔、灌木林地在不同土层的酶生态化学计量特征F i g.4 E c o s t o i c h i o m e t r i cc h a r a c t e r i s t i c so f e n z y m e s i nd i f f e r e n t s o i l l a y e r so fT r e e sa n ds h r u b s i nd i f f e r e n t r e s t o r a t i o ny e a r s801 水 土 保 持 研 究 第3 0卷3.2 固沙林恢复过程土壤酶活性变化的生态效应随着植被恢复,土壤中凋落物增加,促进微生物分泌更多胞外酶,使得土壤C,N,P的增加。榆林固沙林植被恢复过程中,土壤中微生物C限制在恢复5 6a达到最大,N限制随着植被恢复最终显著增加,而P限制并不明显。随着植被恢复,土壤C限制逐渐增强的原因可能是土壤N的相对过高,由于土壤中有效氮水平很高,导致C N没有达到有利于微生物生长的这两种元素的阈值比率,土壤C可能不足以支持微生物生长。土壤氮是研究生态系统中微生物生长的一个重要限制因素,而榆林毛乌素沙地固沙林土壤微生物基本上处于N限制,且微生物N限制相对于微生物P限制更为明显,灌木林土壤N限制普遍高于乔木林,这可能与物质的投入和枯落物的C/N密切相关,半固定沙地植被恢复后,以枯落的形式向土壤中输入了大量的C和N,这极大地促进了土壤有机质的积累。但是在不同林地中,有不同的枯落物进入土壤,因此导致土壤中有机质在长期的演变过程中表现出了不同的特征3 2。图5 不同恢复年限下矢量长度(左)与矢量角度(右)F i g.5 R e s p o n s e su n d e rv e c t o r l e n g t h(l e f t)a n dv e c t o ra n g l e(r i g h t)f o rd i f f e r e n t r e c o v e r yy e a r s901第6期 张延文等:固沙林恢复土壤酶化学计量特征与养分限制效应注:B G表示-1,4-葡萄糖苷酶;NA G表示-1,4-乙酰氨基葡萄糖苷酶;L A P表示亮氨酸氨肽酶;A P表示碱性磷酸酶。图6 土壤酶化学计量散点图显示微生物资源限制养分F i g.6 As c a t t e rp l o t o f s o i l e n z y m a t i c s t o i c h i o m e t r y s h o w i n g t h eg e n e r a l p a t t e r no fm i c r o b i a l r e s o u r c e l i m i t a t i o n4 结 论榆林毛乌素沙地随固沙林植被恢复,乔木和灌木林地土壤酶活性均随之显著提升,且土壤中氮转化酶(N A G+L A P)活性对土壤环境的变化反应最敏感,碳转化酶(B G)活性增加速度显著小于氮转化酶(N A G+L A P)和磷转化酶(A P)活性。微生物通过酶对C,N,P转化效应显著增强。同时,植被恢复