免费文献传递   相关文献

Analysis of CO2 Efflux in Soil Profiles after Clear-cutting and Prescribed Burning of Natural Regenerated Castanopsis carlesis Secondary Forest in Subtropical China

米槠人促更新林皆伐火烧后初期土壤CO2通量动态



全 文 :林业科学研究 2016,29(1):53 60
ForestResearch
  文章编号:10011498(2016)01005308
米槠人促更新林皆伐火烧后初期
土壤 CO2通量动态
徐凯健,林伟盛,杨智杰,任 婕,刘小飞,熊德成,谢锦升,杨玉盛
(福建师范大学地理科学学院,湿润亚热带山地生态国家重点实验室培育基地,福建 福州 350007)
收稿日期:20150615
基金项目:国家自然科学基金重点项目“常绿阔叶林人促更新对红壤有机质维持的作用机制”(31130013)
作者简介:徐凯健(1991—),硕士生.主要研究方向:森林生态系统碳循环.Email:phoenixmax@qq.com
 通讯作者:硕士,助理实验师.主要研究方向:森林生态系统碳氮循环.电话:15659192612.Email:weilsnlin@fjnu.edu.cn
摘要:[目的]为科学评价皆伐、火烧对中亚热带常绿阔叶林不同深度土壤有机碳吸存的影响,[方法]以福建省中亚
热带36年生米槠人促更新林为研究对象,采用非散射红外 CO2浓度探测仪和 Licor8100土壤碳通量系统,并结合
Fick扩散法计算并分析0 80cm不同深度土层CO2通量的日动态特征。[结果]表明:(1)火烧地(RB)和皆伐地
(RR)不同土层CO2浓度均出现明显下降,其中,对照(CK)地土壤CO2浓度值(0 80cm)分别是 RB和 RR的1.
9、1.3倍;(2)各试验地土壤CO2通量(0 80cm)表现为RB(1.99μmol·m
2·s-1)>RR(0.99μmol·m2·s-1)
>CK(096μmol·m2·s-1),除20 40cm土层外,RB土壤各层CO2通量均显著高于RR和 CK(P<005);(3)
试验地不同土层CO2通量(0 80cm)日变化幅度表现为 RB>RR>CK,其中,RB土壤各层的变化幅度均显著大
于RR和CK(P<0.05),而RR与CK间的差异表现在0 5、10 20、20 40cm土层(P<0.05);(4)拟合分析表
明,各试验地不同深度土壤CO2通量与土壤温度呈显著相关,且 RB的决定系数(R
2)显著高于RR和CK;不同试验
地各土层温度、含水量的双因素模型拟合效果均优于单因素模型;Q10值显示,皆伐、火烧后初期土壤各层的温度敏
感性得到明显提高。
关键词:土壤CO2通量;土壤深度;皆伐;火烧;土壤CO2浓度
中图分类号:S714 文献标识码:A
AnalysisofCO2EfluxinSoilProfilesafterClearcuttingandPrescribed
BurningofNaturalRegeneratedCastanopsiscarlesisSecondaryForest
inSubtropicalChina
XUKaijian,LINWeisheng,YANGZhijie,RENjie,LIUXiaofei,XIONGDecheng,XIEJinsheng,YANGYusheng
(ColegeofGeographicalSciences,StateKeyLaboratoryBreedingBaseofHumidSubtropicalMountainEcology,
FujianNormalUniversity,Fuzhou 350007,Fujian,China)
Abstract:Clearcuting&prescribedburningisoneofthemostcommonmanagementmeasuresintheeconomicde
velopmentofChineseforestry,whichnotonlyhasastronginterferenceinwoodland,butalsoleadstothechangesin
soilphysicalandchemicalproperties.Inordertoevaluatetheimportanceofclearcutingandprescribedburningon
thesoilorganiccarbonuptakeandstorageofdiferentdepthsinmidsubtropicalevergreenbroadleavedforest,tak
ing36yearoldnaturalregeneratedCastanopsiscarlesissecondarystandinFujianProvinceastheobject,theFick’
sdifusionmethodwereusedtocalculatethedailydynamicchangesofsoilCO2fluxfrom0to80cmdepthofsoil.
AndtherelationshipbetweentheCO2fluxandsoiltemperature,moistureinsoilprofileswasalsoanalyzedbyusing
regressionmodels.Theresultsareasfolows:(1)TheCO2concentrationindiferentsoillayersdecreasedsignifi
林 业 科 学 研 究 第29卷
cantlyafterclearcutingandprescribedburning,andtheclosertothesurface,thelargerthedecline.From0to80
cmsoildepth,themeandailyCO2concentrationinthecontroled(CK)plotwas1.9timesinresidueburnt(RB)
plotand1.3timesinresidueretention(RR)plot.AndbothinRBandRRplot,diferentdepthofsoilCO2con
centrationdecreasedobviously.(2)ThedailydynamicchangesofCO2fluxinsoilprofileshowedamonopeak
trend,withitsmaximumvaluemainlyfrom12:30amto15:30pmaswelastheminimumfrom6:30pmto9:30
pm,themeandailysoilrespirationinRBplot(1.99μmol·m2·s-1)wassignificantlyhigherthanthatinRRplot
(0.99μmol·m2·s-1)andCKplot(0.96μmol·m2·s-1),thewholesoilCO2fluxinalthetestplotsde
creasedwiththeincreasingofsoildepth.ThedailyaverageCO2fluxinsoilprofileofRBplotwassignificantlyhigher
thantheothers,exceptforthesoildepthfrom20to40cm.AndtherewasnosignificantdiferenceondailyCO2
fluxinsoilprofilesofRRandCKplot.(3)Inaddition,thediurnalrangeofsoilCO2fluxinRBplot(147%)was
higherthanthatintheRR(99%)andtheCKplot(66%),andtherearesomesignificantdiferencesofRBvs.
RRandCK(P<0.05),whiletheRRvs.CKalsoexistedsignificantdiference(P<0.05)in0-5cm,10-20
cmand20-40cmsoillayer.(4)ByplotingCO2fluxvs.soiltemperature,itwasfoundthattheCO2efluxofdif
ferentdepthsoilcorelatedexponentialywithsoiltemperatureinalplots,andtheQ10valueshowedthatthesensi
tivityofsoiltemperatureincreasedsignificantlyafterclearcutingandprescribedburning.Butthesoilmoisturewas
notsignificantcorelatedwithCO2eflux.Andthehybridmodeloftemperatureandsoilmoisturewasprovedtobea
moreappropriatepredictorofCO2dynamicstateinaltestingplots.Consequently,clearcuting&prescribedburn
ingcoulddefinitelyafectCO2fluxindeepsoilaswelasitsdailydynamic,andwouldalsochangethesoilCO2flux
ofsoiltemperatureandmoisturecontentofdiferentsoillayerresponsemechanism.
Keywords:soilCO2eflux;soildepth;clearcuting;prescribedburning;soilCO2concentration
土壤呼吸是土壤向大气中释放 CO2的过程,是
全球陆地碳循环的重要环节[1]。全球每年由土壤呼
吸产生的碳通量远高于使用化石燃料产生的碳释放
总量[2]。土壤中的CO2来源可分为土壤动物、植物
根系和微生物的自养呼吸,以及土壤有机质和植物
凋落物分解产生的异氧呼吸两部分[3]。一般土壤
CO2浓度伴随土壤深度的增加而增加
[4],土壤内部
CO2浓度的垂直差异导致其不断向地表移动,因此,
土壤呼吸是不同深度土壤CO2浓度和扩散机制共同
作用的结果。深层土壤(20cm以下)的有机碳储量
通常占土壤总量的一半以上,是土壤碳循环中的关
键部分[5];同时,土壤有机碳的矿化过程在土壤底层
和表层(0 5cm)也同等重要[6];然而,当前大多数
的相关研究只关注土壤表层的呼吸动态,忽略对深
层部分的探索[7]。深层土壤在环境因子(温度、湿
度和氧气)、土壤因子(尤其是养分有效性)、微生物
以及土壤有机碳自身稳定性等方面均与表层土壤存
在明显的差异[5,8-10],因此,研究深层土壤呼吸对外
界环境变化的响应尤为重要。
皆伐、火烧是我国南方经济林造林最常用的营
林措施,其结果不仅对林地产生强烈干扰,同时导致
土壤的理化性质发生改变[11]。近年来,已有部分学
者针对皆伐、火烧开展相关研究[11-14],但对皆伐、火
烧后土壤呼吸变化,特别是针对不同深度土壤呼吸
的研究仍然缺乏。因此,本文使用SC8000便携式土
壤CO2测定仪,结合垂直扩散模型
[15-16]对不同土壤
深度(10、20、40、60、80cm)的 CO2通量进行估算,
探讨皆伐(保留采伐剩余物)和火烧(无采伐剩余
物)对中亚热带米槠(Castanopsiscarlesi(Hemsl.)
Hay.)次生林土壤呼吸日动态的影响,以期更科学
地评价皆伐、火烧等营林措施对亚热带常绿阔叶林
地土壤碳吸存动态的影响。
1 材料与方法
1.1 研究区概况
试验样地位于福建省三明市金丝湾的陈大林业
采育场(26°19′N,117°36′E),样地总面积达17.1
hm2。地处戴云山脉西北、武夷山脉东南面,属中亚
热带季风性气候,年均降水量1749mm,年均气温
19.1℃。地形以低山丘陵为主,平均海拔330m,平
均坡度25° 35°,该区的土壤厚度超过1m,且主要
是由黑云母花岗岩发育的红壤。土壤部分理化性质
见表1。试验地的拟采伐林分为1978年经强度择伐
后人促更新的米槠次生林,该林分密度为2650株
45
第1期 徐凯健,等:米槠人促更新林皆伐火烧后初期土壤CO2通量动态
·hm-2,平均树高19.7m,平均胸径13.5cm。主要
树种 为 米 槠、木 荷 (SchimasuperbaGardn.et
Champ.)和闽粤栲(CastanopsisfisaRehd.etWils)
等,其中,米槠为优势树种。林下植被主要是毛冬青
(IlexpubescensHook.etArn.)、狗骨柴(Diplospora
dubia(Lindl.)Masam.)、矩形叶鼠刺(Iteachinensis
Hook.etArn.)和狗脊蕨(Woodwardiajaponica(L.
f.)Sm.)等。
表1 试验地土壤本底属性(0 80cm)
土层/cm
有机C含量/
(g·kg-1)
C储量/
(t·hm-1)
全N/
(g·kg-1)
酸碱度
pH
土壤密度/
(g·cm-3)
微生物量碳/
(g·kg-1)
细根生物量/
(g·m-2)
0 10 28.32±2.51 25.65±2.37 1.74±0.17 4.4±0.05 0.95±0.03 476.03±33.20 141.16±7.13
10 20 15.32±4.10 14.57±2.15 1.24±0.50 4.4±0.11 1.04±0.03 377.39±35.20 59.90±5.28
20 40 11.50±0.98 22.16±4.63 0.61±0.01 4.3±0.27 1.20±0.06 317.34±43.10 70.80±4.99
40 60 8.21±0.61 19.23±6.40 0.50±0.07 4.4±0.20 1.31±0.07 220.86±32.10 39.36±5.11
60 80 6.25±0.07 17.19±3.22 0.380±0.11 4.4±0.06 1.37±0.07 133.80±28.40 41.89±2.87
  注:表中数据为均值±标准差。
1.2 研究方法
1.2.1 样地设置 2012年1—2月对拟采伐米槠人
促更新次生林进行皆伐,并将采伐剩余物设置2种
处理,即保留(RR)和火烧(RB),火烧时间为2012
年3月28日。每种处理面积约20m×20m,各处理
设置3个重复,6块样地的立地条件相似,按随机区
组排列,试验同时选择未皆伐样地为对照(CK)。
1.2.2 土壤剖面的 CO2浓度观测 皆伐、火烧后,
于2012年4月初在各试验样地内分别挖取深度1m
的标准剖面,使用环刀法取得土样以测量土壤含水
率、孔隙度、有机碳含量及微生物量等指标,结束后
修平剖面,并在10、20、40、60、80cm深度分别插入
PVC管,布设监测土壤 CO2浓度装置。实验装置布
设后稳定一段时间,于2012年4月21日至27日,
对每块试验地分别进行7d每天持续24h的观测,
观测期间天气条件稳定、无降雨干扰。选用 SC8000
便携式土壤CO2测定仪(GMT220系列,芬兰)连续
监测土壤各深度土层的 CO2浓度。该仪器由主机、
连接线和CO2检测器三部分组成,其中,检测器采用
新型硅基非散射红外CO2设备(NDIR),能够连续测
量土壤CO2浓度变化,同时避免灰尘、水汽等的影
响。在测定CO2浓度时,将检测器安置在不同的土
壤深度,然后连续测量不同土层CO2浓度,并自动保
存数据,仪器测定时间步长为30min。
1.2.3 土壤剖面 CO2通量计算 在测得土壤 CO2
浓度后,依据Fick气体扩散法计算土壤剖面 CO2通
量(FS):
Fs=-Ds
ΔC(z)
Δz
(1)
  式(1)中:Ds为土壤中 CO2的扩散系数(m
2·
s-1);C为深度z(m)土壤的CO2浓度(μmol·m
-3)。
Ds=εDα (2)
  式(2)中:ε为相对气体扩散系数;Da为自由大
气CO2扩散系数
[17]。关于ε的计算模型较多,课题
组的前期研究表明[18],使用 Moldrup扩散模型计算
的土壤表层CO2通量值与使用 Li8100腔室法的实
测结果值最接近,并达极显著相关(R2=0.64,p<
0.01)[19]。因此,本文选用该模型估算土壤剖面
CO2通量的变化。Moldrup模型的计算公式为:
ε=(φ-θ)
2.5
φ
(3)
  式(3)中:θ为土壤体积含水量(cm3·cm-3),φ
为土壤孔隙度[17]。
1.2.4 土壤表层的 CO2通量观测 采用 Li8100
腔室法同时段对各试验样地的土壤表层(0 5cm)
CO2通量进行日动态观测。
1.2.5 土壤剖面的温度与含水量观测 在土壤垂
直剖面各土层中分别插入 2支温度和水分探头
(H2OModelEC5,英国),进行实时持续观测,其时
间步长与土壤CO2浓度观测一致。
1.3 数据处理
使用SPSS19.0和 Origin8.0软件进行数据处理
及图形绘制。对所用数据均进行前处理,并通过正态
分布检验和F检验,使用剔除异常值后的均值进行计
算和分析。采用指数模型(R=a×ebT)、线性模型(R
=a×W+b)及双因素模型(R=a×ebT×Wc)(式中:R
为土壤呼吸通量;a、b、c为参数;T为土壤温度,W为
土壤含水量)拟合不同营林措施下土壤剖面CO2通量
与土壤温度、含水量的关系,并计算土壤各层的温度
敏感性指数Q10值(Q10=R
(T+10)/RT=e10b)[20]。
55
林 业 科 学 研 究 第29卷
2 结果与分析
2.1 不同更新方式下土壤剖面CO2浓度日变化
对RR、RB及 CK进行土壤 CO2浓度动态监测
结果(图1)显示:各试验地土壤 CO2浓度均随土壤
深度的增加而增大,其中,CK土壤CO2浓度(0 80
cm)日变化范围为3099 14250μmol·mol-1,RR
和RB土壤CO2浓度(0 80cm)日变化范围则分
别为1715 10400、798 9000μmol·mol-1,CK
的土壤CO2浓度(0 80cm)日均值分别是皆伐和
火烧后的1.3倍和1.9倍。CK各土层(10、20、40、
60、80cm,下同)的 CO2浓度日均值分别是3419、
5350、7374、9383、12427μmol·mol-1;与 CK相
比,RR各土层 CO2浓度值依次下降了 24.83%、
2094%、18.33%、20.01%、29.29%;而 RB各土层
CO2浓度值则分别下降了 66.78%、46.02%、
4752%、43.43%、50.53%。此外,各试验地不同深
度土层CO2浓度均大体呈单峰变化曲线,但峰值的
出现时间存在一定差异,表明皆伐和火烧处理初期
并未改变深层土壤CO2通量的日变化趋势。
图1 不同营林措施下土壤剖面的CO2浓度、土壤CO2通量、土壤温度及土壤含水量日动态
65
第1期 徐凯健,等:米槠人促更新林皆伐火烧后初期土壤CO2通量动态
2.2 不同试验地土壤CO2通量变化及垂直分布
利用Fick扩散法计算 RB、RR和 CK各层土壤
CO2通量,结果(图1)显示:CO2通量日变化均呈单
峰型,其峰值出现在12:30—15:30。各试验地土壤
CO2通量的差异体现在0 5、10 20、40 60、60
80cm土层中,RB各土层日均 CO2通量(分别为
5.67、2.86、0.82、0.39、0.23μmol·m-2·s-1)均显
著高于 RR(分别为 2.27、1.29、0.81、029、0.16
μmol·m-2·s-1)和 CK(分别为2.29、1.42、0.74、
0.33、0.2μmol·m-2·s-1;P<0.05),而RR与 CK
各层土壤日均CO2通量差异不显著(P>0.05),表
明火烧处理初期导致深层土壤CO2通量产生较大变
化,而皆伐处理初期并未产生类似的结果。此外,
RB、RR、CK土壤(0 80cm)的 CO2日通量分别为
1.99、0.99、0.96μmol·m-2·s-1,且火烧后土壤
CO2通量更集中于土壤表层(比例由 46%提升至
57%),这可能与不同处理初期各自土壤理化特性的
变化程度有关。
此外,RB、RR和CK土壤(0 80cm)CO2通量
的日变化幅度(以最大值与最小值之差除以平均值,
下同)分别为147%、99%和66%,且 RB在土壤各
层的CO2通量日变化幅度均显著高于 RR和 CK(P
<0.05;图2),RR与CK在0 5、10 20、20 40
cm土层的日变化幅度存在显著差异(P<0.05)。
由此可见,皆伐和火烧初期均使不同深度土壤的
CO2通量日动态过程产生明显扰动,并加剧其变化
过程的不稳定性。
图2 不同深度土壤CO2通量与日变化幅度分布(图中不同大写字母表示同一深度土层不同处理间差异显著(p<0.05))
2.3 土壤CO2通量与土壤温度、含水量的关系
2.3.1 土壤温度和含水量的日动态 RB、RR和
CK各土层之间土壤温度的日变化差异较大(图1),
其中,CK各土层温度的日变化相对平稳,无明显峰
值;而RB和RR与其自身 CO2通量变化模式类似,
呈明显的单峰型趋势(峰值时间为9:30—10:30);
各试验地各土层(0 10、10 20、20 40、40
60、60 80cm)土壤温度日变化幅度显示为RB(分
别为74.8%、52.0%、45.5%、45.0%、36.3%)>RR
(分别为 23.2%、26.5%、17.7%、18.2%、16.3%)
>CK(分别为 4.98%、5.21%、4.97%、4.52%、
4.11%)。不同试验地土壤含水量的日变化也存在
明显差异,其中,RR各土层含水量日变化幅度较大
(21% 32%),CK其次(5% 17%),而 RB则相
对平稳(4% 12%);在10、20、40cm土层,CK的平
均含水量显著高于RR与 RB(P<0.05),随着土壤
深度的下降,各试验地相同土层的平均含水量趋于
一致。
2.3.2 土壤温度、含水量对土壤CO2通量的影响 
运用指数模型(R=a×ebT)和线性模型(R=a×W+
b)分别对土壤各层 CO2通量和土壤温度、含水量进
行拟合,结果(表 2)显示:各试验地土壤(0 80
cm)CO2通量与土壤温度呈较好的指数关系(p<
001),其R2分别为0.540 0.797(RB)、0.423
0.608(RR)和0.128 0.380(CK),整体表现为RB
>RR>CK;随土壤深度的增加,CK土壤的 Q10值呈
上升趋势,但皆伐、火烧后部分土层的 Q10值垂直变
化趋势不明显(如 RB、RB的 60cm层),且 RR和
RB不同深度土壤的 Q10值均显著高于 CK(p<
005)。RR、RB和CK各土层 CO2通量与含水量线
性模型的R2普遍较低,但在相关性上,RB除80cm
75
林 业 科 学 研 究 第29卷
土层外均达极显著相关(P<0.01),CK除20、60cm
土层外,其余土层也达到极显著相关(P<0.01);火
烧后(RB)土壤(0 60cm)CO2通量与含水量的相
关系数(R2)均明显提高(P<0.01)。此外,采用土
壤温度与含水量的双因素模型(R=a×ebT×Wc),可
分别解释RB、RR和 CK各土壤深度土壤 CO2通量
变化的53.6% 96.1%、41.0% 72.9%和39.7%
65.4%,拟合结果明显优于以温度为单因子的指
数模型,表明土壤温度与含水量在较大程度上共同
影响土壤CO2通量的变化。
表2 各试验地不同深度土壤CO2通量与土壤温度、含水量的相关关系
试验地
土壤深度/
cm
土壤温度
参数
a b Q10 R2
土壤含水量
参数
a b R2
双因素
R2
RB 10 1.547 0.061 1.84 0.797 197.73 -29.47 0.285 0.961
20 0.117 0.113 3.10 0.718 -511.03 97.40 0.647 0.876
40 0.043 0.146 4.31 0.732 48.77 -8.81 0.621 0.760
60 0.070 0.084 2.32 0.540 29.71 -5.95 0.213 0.675
80 0.025 0.136 3.90 0.572 -1.10 0.48 0.010 0.536
RR 10 0.113 0.141 4.10 0.608 20.44 -1.76 0.188 0.729
20 0.106 0.151 4.53 0.564 -3.41 1.99 0.060 0.540
40 0.031 0.158 4.85 0.605 -0.66 0.96 0.030 0.597
60 0.027 0.116 3.19 0.596 1.81 -0.14 0.298 0.626
80 0.010 0.302 — 0.423 -0.70 0.34 0.022 0.410
CK 10 0.128 0.025 1.28 0.380 -11.91 5.59 0.132 0.513
20 0.720 0.071 2.03 0.340 17.63 -3.19 0.028 0.405
40 0.613 0.096 2.61 0.296 -10.04 3.05 0.177 0.654
60 0.151 0.107 2.92 0.337 -2.42 0.93 0.062 0.397
80 0.080 0.133 3.78 0.128 -14.59 3.94 0.354 0.466
  注:表示差异显著(P<0.05),表示差异极显著(P<0.01),—表示数值存在异常。
3 结论与讨论
3.1 土壤CO2浓度垂直梯度的日动态
本研究中,RR、RB和 CK的土壤 CO2浓度均随
土层深度的增加而上升,这与多数关于森林土壤剖
面CO2浓度的研究结论一致
[20-21]。在整个观测期
间,CK土壤(0 80cm)的CO2浓度日均值(0 80
cm)分别是RR和RB的1.3、1.9倍。同时与CK相
比,皆伐、火烧后土壤各层的CO2浓度均发生了不同
幅度的下降,其中,土壤表层(10cm)和底层(80
cm)CO2浓度的降幅最大,RB各土层 CO2浓度的下
降幅度也明显高于RR。这可能与皆伐、火烧处理有
效降低了土壤密度,使土壤内部的微生物活性降低,
加速CO2在土壤中的移动和扩散速率有关
[22-23];而
室内实验的结果表明,皆伐、火烧后初期,其土壤密
度(0 80cm)分别降低了10.65%、828%。对于
RB土壤各层 CO2浓度均显著低于 RR(P<0.05),
则可能是由于其表层土壤团聚体在高温火烧下遭到
破坏[24],导致表层土壤孔隙度进一步扩大,使聚集
在表层以下的土壤 CO2快速扩散至地表,最终造成
土壤深层 CO2浓度的大幅降低。此外,RR、RB和
CK的土壤 CO2日变化趋势均呈单峰型,且土壤表
层CO2浓度的变化幅度普遍高于土壤深层,与类似
的研究结果相近[25]。
3.2 皆伐、火烧对土壤剖面 CO2通量日变化的
影响
  皆伐和火烧作为我国常见的营林措施,会显著
改变林地土壤中植物根系、微生物的组成与活性
等[26],从而使土壤的 CO2通量发生变化。本研究
中,皆伐和火烧处理初期未改变土壤各层 CO2通量
的日变化趋势,RB、RR和 CK土壤的 CO2通量日变
化动态均呈单峰型,与刘乐中等[26]对杉木林的相关
研究结果相似。研究发现,各土层温度的日变化亦
呈单峰型,其变化趋势与对应土层的土壤 CO2通量
变化一致,因此,温度可能是主导土壤CO2通量变化
趋势的主要原因(回归分析显示,各土层土壤温度与
CO2通量普遍呈极显著相关,p<0.01)。此外,RB、
RR和CK土壤(0 80cm)的 CO2日通量表现为
RB(1.99μmol·m-2·s-1)>RR(0.99μmol·m-2
·s-1)>CK(0.96μmol·m-2·s-1)。各试验地之
间土壤CO2通量的差异主要体现在0 5、10 20、
40 60、60 80cm土层中,表现为 RB土壤日均
CO2通量显著高于RR和CK(图2;P<0.05),这与
国内外针对土壤 CO2 通量的相关观测结果相
85
第1期 徐凯健,等:米槠人促更新林皆伐火烧后初期土壤CO2通量动态
近[22,26-28]。火烧处理明显增加了深层土壤的 CO2
通量,这可能是由于火烧后不同土层密度、有机质含
量及微生物数量均低于未烧地[27-28],从而直接或间
接增加了土壤的CO2通量
[29]。Bauhus等[28]研究发
现,林木在火烧后会导致其土壤内微生物的呼吸量
明显上升,从而使 CO2通量大大超过原先水平。另
外,RR、RB和CK在20 40cm土层的CO2通量差
异不显著,这可能与该层通常作为土壤有机层和矿
质层的分界区,其 CO2传输作用强烈有关,王超
等[18]在对杉木林不同深度土壤 CO2通量的研究中
也发现类似结果。一般皆伐处理也会导致矿质土壤
中微生物活性及现存有机质分解速率的有效提
升[30],但本研究中,RR土壤剖面的 CO2通量与 CK
差异不显著,这可能与在皆伐处理初期,其土壤内部
环境的变化具有一定的时滞性有关[22]。
皆伐、火烧处理同样改变了土壤剖面 CO2通量
的日变化幅度,表现为RB高于RR和CK(图2,P<
0.05),而RR与CK在0 5、10 20和20 40cm
土层的变化幅度差异显著(P<0.05),这与目前针
对土壤表层的部分研究结果相近[22-23]。其差异的
产生可能是由于土壤中的细根受到火烧、皆伐的影
响而部分死亡,被微生物分解后释放出大量 CO2,而
不同土壤内分解速率的差异则会显著干扰其CO2通
量的变化[25];同时,森林经历皆伐、火烧失去茂密植
被覆盖物后,也会使其土壤温度和土壤内呼吸速率
变化过程的稳定性明显下降[11],如本研究中 RR、
RB不同土层温度的日变化幅度为 16.3%
748%,CK的变化幅度则不足5.2%。对于 RR与
RB之间土壤CO2通量的显著差异,则可能与RR内
保留了采伐物以及土壤中细根的存活量更多有
关[31],其结果导致了 RB土壤剖面 CO2通量的日变
化幅度增大(P<0.05)。此外,亚热带阔叶森林土
壤表层0 20cm内的CO2通量约占土壤呼吸总量
的74%,而火烧过后该比重会进一步提升至86%,
这可能与RB表层的土壤团聚体在高温火烧后产生
破坏、土壤孔隙度进一步增大,使聚集在土壤表层以
下部分的CO2更多地扩散到地表有关
[22,32]。
3.3 土壤CO2通量对土壤温度、含水量的响应
RR、RB和 CK不同深度土壤 CO2通量与土壤
温度普遍达到极显著相关(P<0.01),这与多数研
究土壤温度与土壤 CO2 通量关系的结论一
致[22-23,26-28]。其中,RB各土层 CO2通量与土壤温
度的R2最高,其能分别解释 RB、RR和 CK土壤各
层CO2通量变化的54% 79.7%、42.3% 60.8%
和12.8% 38%,这可能是由于森林经历火烧后完
全失去地表植被覆盖物,使土壤温度的日变化幅度
明显增大[11],且火烧后土壤孔隙度上升,可能导致
土壤CO2通量更易受外界温度变化的影响
[14],Rsik
等[33]对亚寒带森林土壤的研究也得到类似结果。
另外,林地在经过皆伐、火烧后,其植物根系逐渐死
亡[34-35],土壤CO2的来源逐渐以异氧呼吸方式为主;
而CK的土壤根系呼吸所占比重较大,其根系呼吸的
影响因素相对复杂[36],导致RB、RR的土壤CO2通量
与土壤温度具有更高的相关性。本研究中,RB、RR
和CK在10cm与20cm深度土壤CO2通量与温度的
相关性明显高于80cm土层,这可能是由于土壤剖面
理化性质不均一造成的,因为试验地本底0 20cm
的土壤密度明显偏低(表1)。
相比土壤温度,土壤含水量对CO2通量影响的贡
献率普遍较低,这可能与土壤含水量变化常处在稳定
的阈值范围内[37-38]及在短期内的变化不显著有
关[18],如本研究中土壤含水量的变化幅度明显小于
土壤温度;但土壤剖面CO2通量与土壤温度、含水量
双因素模型的拟合效果十分显著,使用该模型可以共
同解释 RB、RR和 CK的土壤呼吸变化的536%
96.1%、41.0% 72.9%和39.7% 654%,结果明
显优于土壤温度单因素模型,表明含水量可能通过影
响土壤温度,进一步对CO2通量产生影响。此外,我
国森林土壤的Q10值一般为1.33 5.53
[39],本研究中
土壤呼吸Q10值均处于该范围,同时RB、RR的土壤剖
面呼吸Q10值均明显高于CK,表明皆伐、火烧初期土
壤剖面呼吸速率的温度敏感性明显提升。
皆伐、火烧对森林生态系统的影响是一个复杂
的渐进过程,因此,未来还需要结合长期的定位观
测,同时进一步扩展对影响土壤呼吸的微生物群落
及酶等要素的研究,以便更深入地了解营林更新后
其整个生态系统碳吸存的变化过程。
参考文献:
[1]XuZF,TangSS,XiongL,etal.Temperaturesensitivityofsoil
respirationinChina’sforestecosystems:Paternsandcontrols[J].
AppliedSoilEcology,2015,93:105-110.
[2]BondLambertyB,ThomsonA.Temperatureassociatedincreasesin
theglobalsoilrespirationrecord[J].Nature,2010,464(7288):
579-583.
[3]PumpanenJ,IlvesniemiH,KulmalaL,etal.Respirationinborealfor
estsoilasdeterminedfromcarbondioxideconcentrationprofile[J].
SoilScienceSocietyofAmericaJournal,2008,72(5):1187-1196.
95
林 业 科 学 研 究 第29卷
[4]DeSuterTM,SauerTJ,ParkinTB,etal.Asubsurface,closed
loopsystemforsoilcarbondioxideanditsapplicationtothegradient
efluxapproach[J].SoilScienceSocietyofAmericaJournal,2008,
72(1):126-134.
[5]RumpelC,KgelKnabnerI.Deepsoilorganicmater—akeybut
poorlyunderstoodcomponentofterestrialCcycle[J].Plantand
Soil,2011,338(1-2):143-158.
[6]SaloméC,NunanN,PouteauV,etal.Carbondynamicsintopsoil
andinsubsoilmaybecontroledbydiferentregulatorymechanisms
[J].GlobalChangeBiology,2010,16:416-426.
[7]HashimotoS,TanakaN,KumeT,etal.Seasonalityofverticaly
partitionedsoilCO2productionintemperateandtropicalforest[J].
JournalofForestResearch,2007,12(3):209-221.
[8]盛 浩,罗 莎,周 萍,等.土壤 CO2浓度的动态观测模拟
和应用[J].应用生态学报,2012,23(10):2916-2922.
[9]Powlson1DS,Brookes1PC,WhitmoreAP,etal.SoilOrganic
Maters[J].EuropeanJournalofSoilScience,2011,62(1):1-4.
[10]TornMS,LapenisAG,TimofeevA,etal.Organiccarbonandcar
bonisotopesinmodernand100yearoldsoilarchivesoftheRussian
steppe[J].GlobalChangeBiology,2002,8(10):941-953.
[11]GuoJF,YangYS,ChenGS,etal.Efectsofclearcutingand
slashburningonsoilrespirationinChinesefirandevergreenbroad
leavedforestsinmidsubtropicalChina[J].PlantandSoil,2010,
333(1):249-261.
[12]安 静,王文杰,王洪岩,等.人工林皆伐对土壤碳及相关理
化性质的影响[J].东北林业大学学报,2012,40(9):57-63.
[13]MaranonjimenezS,CastroJ,KowalskiAS,etal.Postfiresoil
respirationinrelationtoburntwoodmanagementinaMediteranean
mountainecosystem[J].ForestEcologyandManagement,2011,
261(8):1436-1447.
[14]KulmalaL,AaltonenH,BerningerF,etal.Changesinbiogeo
chemistryandcarbonfluxesinaborealforestaftertheclearcuting
andpartialburningofslash[J].AgriculturalandForestMeteorolo
gy,2014,188:33-44.
[15]DeSuterTM,SauerTJ,ParkinTB,etal.Asubsurface,closed
loopsystemforsoilcarbondioxideanditsapplicationtothegradi
entefluxapproach[J].SoilScienceSocietyofAmericaJournal,
2008,72(1):126-134.
[16]LiangN,HiranoT,ZhengZM,etal.Continuousmeasurementof
soilCO2efluxinalarchforestbyautomatedchamberandconcen
trationgradienttechniques[J].BiogeosciencesDiscussions,2010,
7(1):1345-1375.
[17]JonesHG.PlantsandMicroclimate:AQuantitativeApproachto
EnvironmentalPlantPhysiology[M].2nded.NewYork:Cam
bridgeUniversityPress,1992:51-51.
[18]王 超,黄群斌,杨智杰,等.杉木人工林不同深度土壤 CO2
通量[J].生态学报,2011,31(19):5711-5719.
[19]MoldrupP,OlesenT,GamstJ,etal.Predictingthegasdifusioncoef
ficientinrepackedsoil:waterinducedlinearreductionmodel[J].
SoilScienceSocietyofAmericaJournal,2000,64(5):1588-1594.
[20]DavidsonEA,KathleenE,SusanE,etal.Verticalpartitioningof
CO2productionwithinatemperateforestsoil[J].GlobalChange
Biology,2006,12:944-956.
[21]HiranoT,KimH,TanakaY.Longtermhalfhourlymeasurementof
soilCO2concentrationandsoilrespirationinatemperatedeciduous
forest[J].JournalofGeophysicalResearch,2003,108:1-13.
[22]林伟盛,杨智杰,郭剑芬,等.米槠天然更新次生林皆伐火烧
后初期土壤呼吸日变化[J].福建师范大学学报:自然科学版,
2013,29(2):68-74.
[23]WuX,BrüggemannN.Longtermefectsofclearcutingandselec
tivecutingonsoilmethanefluxesinatemperatespruceforestin
southernGermany[J].EnvironmentalPolution,2011,159(10):
2467-2475.
[24]RyanMG,LawBE.Interpreting,measuring,andmodelingsoil
respiration[J].Biogeochemistry,2005,73(1):3-27.
[25]SchlesingerWH,AndrewsJA.Soilrespirationandtheglobalcar
boncycle[J].Biogeoehemistry,2000,48(1):7-20.
[26]刘乐中,杨玉盛,郭剑芬,等.杉木人工林皆伐火烧后土壤呼
吸研究[J].亚热带资源与环境学报,2008,3(1):8-14.
[27]HatenJA,ZabowskiD,OgdenA,etal.Soilorganicmaterina
ponderosapineforestwithvaryingseasonsandintervalsofpre
scribedburn[J].ForestEcologyandManagement,2008,255
(7):2555-2565.
[28]WanthongchaiK,BauhusJ,GoldammerJG.Nutrientlosses
throughprescribedburningofabovegroundliterandunderstoryin
drydipterocarpforestsofdiferentfirehistory[J].Catena,2008,
74(3):321-332.
[29]HartleyIP,HeinemeyerA,EvansSP,etal.Theefectofsoil
warmingonbulksoilvs.rhizosphererespiration[J].Global
ChangeBiology,2007,13(12):2654-2667.
[30]KimC.SoilCO2efluxinclearcutanduncutredpine(Pinusden
sifloraS.etZ.)standsinKorea[J].ForestEcologyandManage
ment,2008,255(8-9):3318-3321.
[31]杨玉盛,陈光水,王小国,等.皆伐对杉木人工林土壤呼吸的
影响[J].土壤学报,2005,42(4):584-590.
[32]GaudinskiJB,TrumboreSE,DavidsonEA,etal.Soilcarbon
cyclinginatemperateforest:Radiocarbonbasedestimatesofresi
dencetimes,sequestrationratesandpartitioningoffluxes[J].Bio
geochemistry,2000,51:33-69.
[33]RiskD,KelmanL,BeltramiH.Carbondioxideinsoilprofiles:
productionandtemperaturedependence[J].GeophysicalResearch
Leters,2002,29(6):11-14.
[34]杨玉盛,董 彬,谢锦升,等.林木根呼吸及测定方法进展
[J].植物生态学报,2004,28(3):426-434.
[35]HgbergP,NordgrenA,BuchmannN,etal.Largescaleforest
girdlingshowsthatcurentphotosynthesisdrivessoilrespiration
[J].Nature,2001,411(6839):789-792.
[36]TangJW,LaurentM,AlexanderG,etal.Continuousmeasure
mentsofsoilrespirationwithandwithoutrootsinaponderosapine
plantationintheSieraNevadaMountains[J].Agriculturaland
ForestMeteorology,2005,132(3/4):212-227.
[37]DavidsonEA,SavageKE,TrumboreSE,etal.Verticalpartitio
ningofCO2productionwithinatemperateforestsoil[J].Global
ChangeBiology,2006,12(6):944-956.
[38]HeinemeyerA,DiBC,LloydAR,etal.Soilrespiration:impli
cationsoftheplantsoilcontinuumandrespirationchambercolar
insertiondepthonmeasurementandmodelingofsoilCO2eflux
ratesinthreeecosystems[J].EuropeanJournalofSoilScience,
2011,62(1):82-94.
[39]陈光水,杨玉盛,吕萍萍,等.中国森林土壤呼吸模式[J].生
态学报,2008,28(4):1748-1761.
(责任编辑:徐玉秀)
06