大兴安岭北部多年冻土区土壤碳氮含量及有机碳矿化特征

2021-10-09董星丰赵光影杨盛东臧淑英

生态学报 2021年17期

董星丰,赵光影,李苗,刘超,王迪,杨盛东,臧淑英

哈尔滨师范大学,寒区地理环境监测与空间信息服务黑龙江省重点实验室, 哈尔滨 150025

长期冷湿的环境使得北方多年冻土区积累了大量的碳,其储量约为当前大气中碳含量的二倍[1]。全球气候日益暖化情形下,多年冻土区生态环境稳定性和生物地球化学循环过程必将发生显著变化[2- 3]。众多学者认为,冻土退化将改变土壤水热状态,进而影响土壤碳氮的周转,导致储存在多年冻土中的有机碳(SOC)以温室气体的形式进入到大气中,并对气候变化起到正向反馈作用[4-6]。因此气候变暖情形下,多年冻土区土壤碳动态变化强度及其受控因素仍需进一步明确。

随着温度升高和冻土融化深度的增加,深层土壤碳输出贡献是评价土壤碳动态的关键[7]。Song等[8]研究指出长期冻结的深层泥炭在温度升高后生成的含碳温室气体释放量显著增加。Chen等[9]发现活动层中土壤碳的释放主要取决于可用基质,而多年冻土层中土壤碳的释放主要受控于微生物作用。Mu等[6]对青藏高原北缘黑河流域上流多年冻土有机碳分解的研究表明,深层土壤对水热变化的响应强度并不次于浅层土壤。可见多年冻土深层土壤有机碳已然暴露于微生物的分解。然而目前关于土壤有机碳储量及潜在排放量的研究主要集中在浅层土壤,对较深土层知之甚少,尤其是我国高纬度多年冻土区深层土壤碳库动态变化缺乏突破性研究进展。众所周知,不同深度土壤的物理化学性质不同,可能具有不同的微生物群落、丰度和活性[10-11]。森林生态系统碳储量巨大,其细微变化可能会显著影响地球系统温室效应,是碳循环相关研究的重点关注对象之一。因此,了解林区土壤中的碳储量、分布规律及其潜在排放量对于预测土壤碳库与全球变暖间的反馈关系意义重大。

大兴安岭作为我国第二大冻土区,现存多年冻土面积约3.9×105km2,冻土空间破碎化严重,活动层厚度逐年加深,热稳定性较差[13],是当前研究土壤生物地球化学过程的理想地区。然而该地区深层土壤的碳氮分布及其对水热变化的响应特征尚不清楚。本研究采集0—6 m土壤,分析了土壤碳氮随土壤深度的分布情况及其与土壤深度间的相互关系,通过室内模拟实验明确了多年冻土区土壤有机碳矿化对温度和水分变化的响应特征。该研究可为全面认识多年冻土退化背景下区域碳循环及预测气候变化提供基础支持。

1 材料与方法

1.1 研究区概况

大兴安岭属欧亚多年冻土带的南缘,主要分布有大面积的森林和湿地。本研究于大兴安岭漠河林场的黑龙江省漠河森林生态系统野外观测站试验区内选择落叶松-灌丛湿地(53°28′ N, 122°20′ E)为研究对象。平均海拔约298 m,优势植被为兴安落叶松(Larixgmelini),20 m×20 m 样方内平均树高约13.59 m,胸径为14.54 cm,郁闭度为0.6。灌木分布较少,主要由狭叶杜香(Ledumpalustre)、笃斯越橘(Vacciniumuliginosum)和杜鹃(Rhododendrondauricum)组成,草本以白毛羊胡子草(Eriophorumvaginatum)构成的直径不一的“塔头”为主。土壤以褐色针叶林土为主,同时存在泥炭及沼泽土。活动层土壤从4月中旬开始融化,至9月中旬融化深度达到最大(约为130 cm),活动层之下土壤富含冰晶,温度常年保持在0度以下。研究区受控于寒温带大陆性气候,监测数据显示2016—2017年均温为-3.3℃,年均降水量为442.95 mm,其中约66.47%的降水量以降雨形式集中在6—8月[14]。

1.2 样品采集

2017年3月于试验区内利用地质钻机采集一根长50 m,直径16 cm 的土柱(其中6.65 m以下为岩石层,土壤发育较差)。将土柱切割成30—40 cm长后按照常规取样方法装袋贴签,运送至实验室后放置于-80℃冷冻冰箱保存以备后续试验。因该土柱6 m以下砾石较多,选择0—6 m土壤样品进行试验,样品分层设置为:活动层土壤样品(0—20、20—40、40—60、60—80、80—100、100—120和120—150 cm),多年冻土层土壤样品(150—180、180—210、210—240、240—270、270—300、300—330、330—360、360—390、390—420、420—450、450—480、480—510、510—540、540—570和570—600 cm)[15]。各层样品手动去除植物根茎及碎石后用以测定土壤理化性质和培养实验。

1.3 样品测定

室内模拟实验的设计参考文献[16]。即采用温度与水分两种变量的交互实验(温度:T1,5℃、T2,10℃、T3,15℃；质量含水量:W1,30%、W2,45%、W3,60%)。单个样品用量为20 g,培养瓶为250 mL规格的玻璃瓶,称取样品后放入设定好温度的培养箱进行为期60 d的实验,期间用连接有三通阀的橡胶塞密封玻璃瓶防止气体漏出,共603个试验样品(22层×3温度×3水分×3重复+9空白)。培养过程中每3 d称重,以补偿水分损失。气体样品采集时间为实验的第1、3、6、9、15、21、28、36、45和60 d。采集之前打开封闭的培养瓶15 min并辅以注射器抽取以使培养瓶内外气体浓度相同。开始培养时关闭三通阀,将培养瓶放回恒温培养箱中,1 h后接通三通阀并使用注射器采集20 mL气体进行浓度分析(GC，Agilent 7890B，USA)。

1.4 数据分析

本研究中土壤有机碳矿化速率以CO2产生速率表示,具体计算过程见公式(1)[17]。土壤有机碳累积矿化量由公式(2)得到[18]:

F=ρ×ΔC×V×273/(273+T)×W

(1)

(2)

式中,F为土壤CO2产生速率(mg kg-1h-1)；ρ为标准状态下CO2的浓度(g/L)；△C为测试样品与对照样品间的浓度差值(mg/m3)；V为实验所用密闭玻璃瓶的容量(mL)；T为培养箱固定温度(℃)；W为实验用土质量(kg)。

土壤有机碳矿化对温度的响应程度用Q10表示[19-20]:

R=αeβΤ

(3)

Q10=e10β

(4)

式中,R为CO2释放速率(mg kg-1h-1)；T为培养温度(℃)；α和β为参数。

使用三因素方差分析法,明确温度、水分、土壤深度及三者的交互作用对有土壤有机碳累积矿化量的影响。利用Pearson相关分析法研究土壤有机碳累积矿化量和土壤理化性质之间的关系。所有统计均在SPSS 16.0中操作。图形处理采用Origin 8完成。

2 结果

2.1 土壤理化性质

图1 土壤理化因子的剖面分布特征Fig.1 Profile distribution characteristics of soil physicochemical factorsSWC: Soil water content，SOC: Soil organic carbon，DOC: Dissolved organic carbon，TN: Total nitrogen，TP: Total phosphorus

表1 土壤有机碳累积矿化量与理化因子间的Person相关性分析

2.2 土壤有机碳矿化特征

由图2可以看出,为期60 d的周期内,30%水分培养下,活动层土壤有机碳累积矿化量平均值在0.31(5℃)至0.77(15℃) mg C之间。多年冻土层土壤有机碳累积矿化量平均值在1.10(5℃)至1.99(15℃) mg C之间。45%水分培养下,活动层土壤有机碳累积矿化量平均值在0.30(5℃)至0.89(15℃) mg C之间。多年冻土层土壤有机碳累积矿化量平均值在0.75(5℃)至1.44(15℃) mg C之间。60%水分培养下,活动层土壤有机碳累积矿化量平均值在0.42(5℃)至0.92(15℃) mg C之间。多年冻土层土壤有机碳累积矿化量平均值在0.86(5℃)至1.58(15℃) mg C之间。

图2 60 d培养后的土壤有机碳累积矿化量Fig.2 Accumulated SOC mineralization after 60-days laboratory incubation experiment

整体来看,培养环境下温度增加有利于土壤有机碳的矿化。其对水分变化的响应比较复杂,呈现先减小(幅度较大)后增大(幅度较小)的特征。对比分析活动层和多年冻土层土壤有机碳累积矿化量,可见土壤有机碳累积矿化量在不同土壤深度上差异明显。各种处理下,有机碳矿化量的高值均出现于360—390 cm 之间,此深度土壤为深层泥炭,低值出现在150—180 cm 之间,此深度为多年冻土层的上界。方差分析表明温度、水分和深度及交互作用对累积矿化量影响显著(P<0.001)(表2)。

表2 土壤深度、温度、水分、及三者的交互作用对土壤有机碳累积碳矿化量的影响

2.3 土壤有机碳矿化的温度敏感性

Q10即温度升高10℃,土壤有机碳矿化程度随之变化的倍数。由图3可知,表层土壤和多年冻土层270 cm处的土壤具有较大的Q10值,分别在2.44—3.44和3.68—4.39之间。而在多年冻土层的顶部,土壤有机碳矿化对增温的响应并不符合指数形式,这使得计算得到的Q10值极低,在0.45—0.62之间波动。尽管土壤水分状况梯度分明,但Q10值在各水分处理之间没有显著差异(P>0.05)。具体来看,不同水分处理下,土壤Q10值在0.45—4.39之间波动,多年冻土层平均值低于活动层土壤,但随着水分的增加,两者的差值由0.85逐渐缩小至0.22。

图3 不同深度土壤有机碳矿化的温度敏感性Fig.3 Temperature sensitivity of the mineralization of SOC (Q10) at different depths

3 讨论

3.1 土壤碳氮分布特征

针对浅层土壤碳循环的相关研究已经取得了坚实的基础,但多年冻土层土壤碳储量及对气候变化的响应研究相对缺乏,将深层土壤碳周转机制与表层土壤混为一谈可能会忽略土壤性状和有机碳随土壤剖面分布的异质性[21-22],未来气候变化背景下深层土壤碳动态的潜在响应仍需进一步探讨。多年冻土层土壤碳库对水热变化的响应是冻土退化背景下区域碳循环及碳平衡预测的较大不确定性因素[23]。为明确土壤碳氮含量的剖面分布规律,本研究于大兴安岭连续多年冻土区钻取土柱进行测定,发现土壤碳氮含量的分布与其赋存机制有关。地表覆盖植被根茎的生长代谢、反复的冻融过程及伴生的复杂微生物活动使得活动层中SOC和TN含量随深度波动较大,而多年冻土层低温厌氧的条件使得SOC和TN较难分解,因此呈现出随深度的增加而增加的规律[24]。

3.2 土壤理化性质对有机碳矿化的影响

3.3 土壤温度、水分和深度对有机碳矿化的影响

一般认为,当土壤湿度适宜且其他因素不受限制时,土壤呼吸作用随温度的升高呈指数模式增长。本研究中,15℃条件下的土壤有机碳矿化量显著高于5℃培养条件(图2),与主流研究结果相符[35-39]。土壤含水量与微生物的氧利用率及土壤气体排放途径密切相关,而碳矿化又受到诸如微生物群落等生物因素影响[40],因此土壤水分状况对碳矿化的影响尚不明确。另外不同深度土壤持水能力的差异或是碳矿化作用与土壤水分关系复杂的原因之一。我们之前的研究表明土壤性状导致的持水差异使得两种林型土壤有机碳矿化随水分变化呈相反的趋势[16]。与前人研究结果一致,我们发现土壤温度、水分及二者交互作用对土壤有机碳累积矿化量影响显著[41-42]。此外,土壤深度是碳矿化的重要影响因素,大兴安岭泥炭地[43]及青藏高原高山草甸[23]的研究支持了这一结果。表2显示土壤有机碳累积矿化量、碳氮含量及深度彼此间都存在显著的相关关系,而有机质被认为是微生物活动的底物,因此可推测土壤碳氮含量随剖面的差异分布特征可能是土壤有机碳矿化量随深度变化的原因。另外土壤碳氮循环相关微生物的数量在不同土层中存在差异,这也是导致土壤有机碳矿化量随深度变化的主要原因[44]。

3.4 土壤有机碳矿化的温度敏感性

受生态系统类型和试验方案等因素综合影响,土壤有机碳矿化的温度敏感性存在较大的时空差异,但仍被视为理解土壤碳分解响应气候变暖的有效指标。试验表明0—6 m土柱各土层的Q10值变化范围为0.45—4.93之间。Q10值在各水分处理间没有显著性差异(P>0.05),指示了土壤水分对碳矿化作用的影响可能通过其他因素体现[40]。活动层中Q10值随深度增加而逐渐减小,平均值约为2.46,与同地区湿地土壤相似[43,45]。多年冻土层中Q10值随深度增加表现出先增后减的规律,与黑河流域上游多年冻土区土壤有机碳矿化温度敏感性的剖面规律类似[6]。增温对活动层底部(120—150 cm)和多年冻土层顶部(150—180 cm)土壤有机碳矿化的影响不同于其他土层,基于指数模型的拟合系数并不理想。这可能是因为多年冻土退化背景下,冻融作用使得活动层底部及多年冻土层顶部土壤水热条件和微生物多样性复杂。其次对应深度处土壤SOC、DOC含量较低(图1),而土壤活性碳是微生物作用的基质,模拟实验的高温环境促使微生物在较短时间内消耗土壤可利用性碳组分,减少了培养后期底物的供给[46],因此观测到的累积矿化量较低。此外,我们的实验观测到了较高的甲烷排放,多年冻土区过渡层 CO2通量的降低或与此深度较高的CH4排放有关。Jiang等[47]关于大兴安岭泥炭地土壤碳排放的研究证实了这一观点。孙晓新等[48]指出活动层深度增加必将影响CH4通量。这是由于水热的变化导致冻结在深层土壤中相对稳定的有机碳转变为易被微生物利用的有机碳,使得产甲烷菌可利用底物增加[49-50]。因此多年冻土区过渡层土壤有机碳矿化过程对温度增加的响应异于其他土层,呈现随温度升高反而减弱的现象。虽然整体来看深层多年冻土中Q10值小于活动层,但是随着水分的增加,两者的差值在逐渐缩小。而冻土层上限面被认为具有隔水层的作用,未来气候持续变暖导致的冻土退化会引起冻土层的消融,进而改变水分运移规律,可能使得多年冻土层土壤含水量增大。这意味着如果未来气候持续变暖,多年冻土层的融化将使得大兴安岭深层土壤有机碳分解潜力变大。