-
氮素是植物生长发育所必需的大量营养元素。研究表明,土壤氮素有效性是森林生产力的主要限制因子之一[1-2]。土壤中氮的有效性受控于决定其数量和形态的有机氮矿化速率、无机氮固定速率和硝化反硝化速率的大小[3]。近年来,随着大气环境污染研究的推进,大气氮沉降及其对森林生态系统的影响日益受到关注,森林土壤氮转化过程及其机理也成为当前生态环境领域的研究热点。
森林类型不同,输入到土壤中的林木凋落物和根系分泌物的种类、数量、分解速率等具有显著差异[4],进而导致森林土壤理化性质[5]以及微生物数量、活性和群体结构显著不同[6-7]。与阔叶林相比,针叶林凋落物中所含的营养物质较少,而木质素和纤维素等难降解物质较多,分解速率较慢,使得土壤有机碳总量和易分解有效碳数量相对较低[8],土壤C/N较高[9]。此外,针叶林凋落物分解过程中产生的有机酸和根系分泌的有机酸通常也会导致土壤pH显著低于阔叶林土壤[9]。这些性质上的差异可能会对微生物参与的土壤氮素转化过程产生显著影响,进而影响土壤中氮的有效性。目前,关于不同森林类型土壤净氮转化速率的研究较多[1,10-11],而针对森林土壤氮初级转化速率的研究只有零星报道[12]。研究发现,森林土壤中铵态氮和硝态氮的周转速率很快,有机质矿化产生的无机氮很快被微生物固定为有机氮,常规研究测定的土壤净氮矿化速率或净硝化速率为负值或者很低的情况实际上并不意味着森林土壤的氮矿化和硝化能力很低[13]。自然条件下,土壤净氮转化速率是控制氮转化的各个初级转化速率的综合结果[14],不能表征土壤氮实际转化速率的大小。因此,研究土壤氮初级转化速率可以从过程水平上深入认识控制森林土壤中无机氮形态和含量的各个氮转化过程的初始速率[15],对于了解森林生态系统生产力、氮素循环,评价森林土壤中氮素的可利用性以及森林经营管理具有重要意义。
针叶林和阔叶林是我国寒温带地区的两种主要林型,前期已有较多关于寒温带森林土壤净氮转化速率的研究[1,11,16],但尚未有从土壤氮初级转化速率角度开展的研究报道。此外,受进入土壤的林木凋落物数量、分解速率及林木根系活力差异的影响,不同深度土壤的理化性质和生物学性质存在差异,进而可能导致土壤氮初级转化速率的不同。因此,本研究以我国东北寒温带的人工红松针叶林(简称针叶林,下同)和天然次生阔叶混交林(简称阔叶林,下同)不同深度的表层土壤为研究对象,采用室内培养试验,通过15N同位素成对标记技术和FLUAZ数值优化模型研究土壤氮初级矿化速率、铵态氮固定速率、硝态氮固定速率、初级硝化速率特征对林型和土壤深度的响应,旨在揭示该地区森林土壤氮素转化特征,了解森林土壤氮供应能力及氮有效性,为森林生态系统的合理经营管理提供科学依据。
-
供试针叶林和阔叶林土壤采自黑龙江省北安市赵光镇,该地区位于47°35′~48°33′ N、126°16′~127°53′ E之间,海拔300~350 m,属寒温带大陆性季风气候,土壤类型为暗沃冷凉淋溶土,成土母质以花岗岩、砂页岩的风化残积物为主,有黏化层和暗沃表层[17-18]。北安市年平均降水量553 mm,年平均气温2 ℃。人工红松针叶林优势树种为红松(Pinus koraiensis Sieb. et Zucc.),其他树种为云杉(Picea asperata Mast.)和冷杉(Abies nephrolepis Maxim.),平均树龄超过50 a,树木间距约2 m。天然次生阔叶混交林主要树种包括桦树(Betula platyphylla Suk.)、杨树(Populus Bonatii Levt.)、糠椴(Tilia mandshurica Rmpr.et Maxim.)、榆树(Ulmus pumila Linn.)和蒙古栎(Quercus mongolica Fisch. ex Ledeb.),平均树龄超过50 a。灌木有毛榛子(Corylus mandshurica Maxim.)、暴马丁香(Syringa reticulate (Blume) H. Hara var. amurensis (Rupr.) J. S. Pringle)、东北鼠李(Rhamnus schneideri Levl.)、金银忍冬(Lonicera maackii (Rupr.) Maxim.)、东北溲疏(Deutzia parviflora Bunge var .amurensis Regel)等。
于2020年10月8日取样,每种林型土壤设置3个空间重复采样区域,空间重复间隔300 m左右,每个空间重复采样区域再设置5个采样点(1 m × 1 m),采样点间隔20 m左右,每个采样点去除表层枯枝落叶后用直径为5 cm的环刀分别采集0~10、10~20 cm土壤,将每个土层的10个样品混匀后作为此土层的一个空间重复样品[11,19]。此混合样品过2 mm筛后一部分新鲜土样置于4 ℃冰箱中储存用于培养试验,其余土样风干后用于土壤理化性质分析,2种林型土壤的基本理化性质见表1。
土壤性质
Soil propertiesCF土层 Soil layer of CF/cm BF土层 Soil layer of BF/cm 0~10 10~20 0~10 10~20 pH 4.34±0.11 b 4.29±0.08 b 5.13±0.19 a 5.06±0.14 a 有机碳 SOC/% 1.93±0.15 b 0.89±0.14 c 4.46±0.57 a 2.14±0.07 b 全氮 TN/% 0.105±0.014 c 0.045±0.015 d 0.357±0.050 a 0.162±0.024 b 碳氮比 C/N 18.6±1.56 a 20.4±3.60 a 12.5±0.91 b 13.4±2.16 b 水溶性有机碳 SOCw/(mg ·kg−1) 132±9.02 c 89.3±8.50 d 487±34.1 a 328±27.0 b 水溶性有机氮 SONw/(mg ·kg−1) 16.1±1.02 c 11.5±0.93 d 36.7±4.35 a 27.9±3.13 b 铵态氮 NH4+-N/(mg ·kg−1) 11.4±0.71 b 4.72±0.46 d 16.2±1.05 a 7.31±0.86 c 硝态氮 NO3−-N/(mg ·kg−1) 7.81±0.28 a 3.43±0.44 c 5.65±0.51 b 2.26±0.31 d 最大持水量 WHC/% 61.7±3.20 b 64.6±4.16 b 82.4±5.97 a 86.5±4.31 a 砂粒 Sand/% 43.2±2.43 a 34.5±2.10 b 32.4±1.76 b 23.8±1.29 c 粉粒 Silt/% 36.3±0.90 c 43.8±1.99 b 43.7±2.74 b 48.4±2.33 a 黏粒 Clay/% 20.5±1.54 c 21.7±0.51 c 23.9±3.35 bc 27.8±3.50 a 注:CF表示人工红松针叶林;BF表示天然次生阔叶混交林。同行不同小写字母表示土层间差异显著(p<0.05)。下同。
Notes: CF indicates Pinus koraienis planation coniferous forest; BF indicates natural secondary broad-leaved mixed forest. Different lowercase letters in the same row indicate significant difference among the four soil layers at 0.05 level. The same below.Table 1. Physico-chemical properties of the studied soils (mean ± SD)
-
本研究采用15N同位素成对标记室内培养试验,每个土壤设置2个15N标记处理:一个处理添加15NH4NO3,一个处理添加NH415NO3,每个处理3次重复。称取30 g(干基)新鲜土壤至250 mL三角瓶中,用移液管向三角瓶中均匀加入2 mL 15N标记的15NH4NO3或NH415NO3溶液(标记丰度为5% atom),使添加的铵态氮浓度达到20 mg·kg−1,同时将土壤水分含量调至最大持水量的60%。将三角瓶用锡箔封口,并在锡箔上扎4~5个洞以利于通气,然后置于25 ℃恒温培养箱中培养7 d,培养期间每隔1 d通过称质量法补充蒸发损失的水分。分别在培养后的0、1、3、5、7 d破坏性取样,向三角瓶中加入75 mL 2 mol·L−1的KCl溶液,250 rpm下振荡提取1 h后过滤,将滤液收集于100 mL塑料瓶中用于分析无机氮浓度及其15N丰度。过滤后的土壤用去离子水淋洗3次以去除残留的KCl,然后将土壤于60 ℃下烘干,球磨仪磨碎后用于分析土壤有机氮浓度及其15N丰度。
-
土壤pH用电位法测定(液土比为2.5:1);土壤有机碳和全氮用CN元素分析仪(vario MACRO cube, Germany)测定;土壤砂粒、粉粒、黏粒含量用比重计法测定;土壤铵态氮、硝态氮含量用MgO-定氮合金蒸馏法测定[20];土壤用水浸提后利用有机碳氮分析仪测定水溶性有机碳和水溶性有机氮含量[21];土壤最大持水量参照Fierer和Schimel的方法测定[22];土壤无机氮和有机氮的15N丰度采用同位素质谱仪(Thermo Fisher MAT 253, USA)测定。
-
土壤氮初级转化速率采用FLUAZ模型拟合得出[23]。将每一次取样获得的标记铵态氮库和硝态氮库的铵态氮、硝态氮、有机氮浓度及其15N丰度分别输入到模型中,通过数值优化拟合出土壤氮初级矿化速率(gm)、初级硝化速率(gn)、初级固定速率(gi)、铵态氮固定速率(ia)、硝态氮固定速率(in)。FLUAZ模型以最小平均加权误差(MWE, mean weighed error)为标准拟合数据,相对于最小平方和而言更加标准化,能够更加准确反映速率之间的变异,估算的土壤氮转化速率与实测值之间具有很好的一致性,而且运用简便、迅速。该模型可以拟合出不同培养时间段的氮初级转化速率,能够反映土壤氮转化的动态变化趋势,本文通过加权平均法计算出了整个培养期间的平均氮初级转化速率。文中数据均为3次重复的平均值,采用Origin软件做图,采用SPSS软件进行不同处理间的差异显著性检验。
-
2种林型土壤中无机氮含量的动态变化见图1。针叶林0~10 cm土壤铵态氮含量在0~1 d略有增加,之后没有显著变化;而10~20 cm土壤铵态氮含量在整个培养期间基本保持不变。阔叶林0~10 cm土壤铵态氮含量随着培养时间的增加显著增加,10~20 cm土壤铵态氮含量增加较缓(图1A)。说明阔叶林土壤有机氮矿化作用较强,有机氮转化为铵态氮的速率高于针叶林土壤。2种林型不同深度土壤中硝态氮含量在培养期间的动态变化很小(图1B),说明土壤硝化速率较低。无论是针叶林还是阔叶林,各个取样时间0~10 cm土壤无机氮含量均显著大于10~20 cm土壤(p<0.05)。无论0~10 cm还是10~20 cm土壤,阔叶林土壤铵态氮含量都显著高于针叶林土壤(p<0.05) (图1A),而阔叶林土壤硝态氮含量则显著低于针叶林土壤(p<0.05) (图1B)。
-
2种林型土壤中无机氮15N丰度的动态变化见图2。在15NH4NO3标记的土壤样品中,针叶林和阔叶林土壤中NH4 + 库的15N丰度随培养时间的增加有不同程度的下降,阔叶林土壤的15N丰度下降速率相对较快,说明有机氮矿化产生了自然丰度的铵态氮稀释了标记的NH4 + 库,阔叶林土壤的有机氮矿化速率比针叶林土壤更快。随着土壤深度的增加,2种林型土壤NH4 + 库的15N丰度下降幅度均变缓(图2A)。在15NH4NO3标记的土壤样品中,NO3−库的15N丰度在整个培养期间没有显著变化(图2B)。在NH415NO3标记的土壤样品中,无论哪种林型和土壤深度,培养期间各处理NH4 + 库的15N丰度均无明显变化(图2C),说明土壤中不存在标记的NO3−异化还原为NH4 + 或已同化的15NO3−再矿化。培养期间各处理土壤中NO3−库的15N丰度随着培养时间的延长逐渐下降,说明土壤中发生了硝化作用产生了自然丰度的NO3−稀释了标记丰度的NO3−库。同一深度下,针叶林土壤的下降幅度相对较大;随着土壤深度增加,2种林型土壤NO3−库的15N丰度下降幅度都有所降低(图2D)。
-
林型和土壤深度对土壤氮初级转化速率具有显著影响,二者的交互作用也显著影响除初级硝化速率之外的其余各个氮初级转化速率(图3)。无论是0~10 cm还是10~20 cm土壤,氮初级矿化速率、初级固定速率、铵态氮固定速率和硝态氮固定速率都表现为阔叶林土壤显著高于针叶林土壤,而初级硝化速率则是针叶林土壤显著高于阔叶林土壤(p<0.05)。无论是针叶林还是阔叶林,随着土壤深度的增加,各个氮初级转化速率都显著下降(p<0.05)。针叶林0~10 cm和10~20 cm土壤初级硝化速率与铵态氮固定速率比值(gn/ia)分别为1.56和1.92,显著高于阔叶林相同土层(0.27、0.15) (p<0.05)。
Figure 3. Average gross N transformation rates and gn/ia ratio in coniferous (CF) and broad-leaved (BF) forest soils at different soil depths
相关分析表明:土壤氮初级矿化速率、初级固定速率、铵态氮和硝态氮固定速率都与土壤pH、SOC、TN、SOCw、SONw和NH4 + -N含量呈显著正相关,初级矿化速率、初级固定速率和铵态氮固定速率都与C/N呈显著负相关(表2)。初级硝化速率与NO3−-N和砂粒含量呈显著正相关,与土壤pH、WHC、粉粒和黏粒含量呈显著负相关(表2)。
土壤性质 Soil properties gm gn gi ia in pH (H2O) 0.648* −0.602* 0.744** 0.751** 0.597* 有机碳 SOC 0.986** −0.027 0.970** 0.922** 0.956** 全氮 TN 0.968** −0.141 0.966** 0.930** 0.902** 碳氮比 C/N −0.641* 0.495 −0.691* −0.693* −0.561 水溶性有机碳 SOCw 0.860** −0.403 0.910** 0.897** 0.802** 水溶性有机氮 SONw 0.866** −0.397 0.879** 0.846** 0.787** 铵态氮 NH4+−N 0.926** 0.324 0.871** 0.824** 0.961** 硝态氮 NO3−-N 0.347 0.830** 0.194 0.175 0.448 最大持水量 WHC 0.493 −0.713** 0.593* 0.561 0.405 砂粒 Sand −0.046 0.936** −0.193 −0.213 0.041 粉粒 Silt 0.023 −0.870** 0.094 0.130 −0.107 黏粒 Clay 0.064 −0.741** 0.268 0.261 0.062 注:**表示p<0.01;*表示p<0.05。
Notes: ** indicates p<0.01; * indicates p<0.05.Table 2. The correlations between soil gross N transformation rates and soil physico-chemical properties
Characteristics of Gross Nitrogen Transformation Rates in the Soils of Pinus koraienis Coniferous Plantations and Natural Secondary Broad-leaved Mixed Forests in Northeast China
- Received Date: 2021-12-16
- Accepted Date: 2022-01-23
- Available Online: 2022-06-20
Abstract: