Distribution of C and N in black soil and its relationship with aggregate stability in sloping land
-
摘要:目的探究黑土区土壤团聚体碳氮分布规律及其与团聚体稳定性的关系，进一步对比耕地与林地团聚体破坏机制差异，为退耕还林以及其他植被恢复途径提高黑土区土地生产力等方面提供理论依据。方法以典型黑土区长期耕作的坡耕地和樟子松坡林地为研究对象，通过不同坡位(坡上、坡中上、坡中、坡中下、坡下)和不同粒级(2~5 mm、1~2 mm、0.5~1 mm、0.25~0.5 mm、 < 0.25 mm)水稳性团聚体总碳、总氮、铵态氮、硝态氮、碳氮比(C/N)及铵态氮与硝态氮比值(ANR)的测定与分析，探究团聚体碳氮分布特征及其与团聚体破坏率之间的相关关系。结果坡林地土壤团聚体有机碳和总氮在坡下沉积且富集在小粒径团聚体中；耕地土壤团聚体有机碳含量及C/N显著低于林地，氮对坡位和粒径变化的响应规律不明显。铵态氮含量在 < 1 mm土壤团聚体中含量较高，硝态氮含量则在>1 mm土壤团聚体中含量较高，且林地铵态氮含量显著高于耕地。耕地团聚体破坏率显著高于林地，以>1 mm粒级团聚体破坏率的表现最为显著，坡上、坡中和坡中下2~5 mm团聚体破坏率显著高于同坡位其他粒径团聚体。团聚体碳含量和团聚体氮含量分别与团聚体破坏率呈负相关和正相关趋势，但均未达显著水平，而综合指标C/N和ANR，以及单因素铵态氮均与黑土团聚体破坏率呈显著负相关关系；当耕地团聚体破坏率超过40%时，ANR、C/N均处于较低的水平。结论坡耕地与坡林地之间，由于养分归还方式与土壤侵蚀环境的不同，导致土壤团聚体碳氮分布规律与团聚体稳定性的差异。坡地黑土团聚体稳定性受团聚体内部有机碳氮的共同作用，与微生物代谢密切相关的综合性指标C/N和ANR，以及被微生物优先利用的铵态氮对黑土团聚体稳定性的影响较其他指标更为显著，可在微生物活性与团聚体稳定性的关系方面做进一步研究。
-
关键词:
- 黑土 /
- 团聚体破坏率 /
- 碳氮比(C/N) /
- 铵态氮与硝态氮比值(ANR) /
- 坡位
Abstract:ObjectiveThis study is to study the relationship between content and distribution of aggregate associated C and N and aggregate stability, and further analyze the difference of aggregate disruption mechanism between farmland and forested land, expectantly provide theoretical foundation for improving soil productivity by returning farmland to forestland and other vegetation rehabilitation measures in black soil region.MethodIn our study, a sloping farmland (preceding crop is corn) and an adjacent sloping forestland (artificial Pinus sylvestris forest) in typical black soil region of northeastern China were selected and each slope was divided into the upper slope position, upper to middle slope position, middle slope position, middle to lower slope position, and lower slope position along the three sampling lines which parallelled to the slope aspect. The soil samples were sieved individually and separated into five aggregate-size fractions (2-5 mm, 1-2 mm, 0.5-1 mm, 0.25-0.5 mm, < 0.25 mm) using the dry sieving and wet sieving method, then the percentage of aggregate disruption (PAD) of each fraction was calculated. Each fraction of soil water stable aggregates was passed through 0.25 mm sieve to determine aggregate-associated organic carbon, total nitrogen, ammonium-N, nitrate-N, and the ratio of ammonium-N to nitrate-N (ANR).ResultIn sloping forestland, aggregate-associated C and N deposited in the lower slope position and accumulated in smaller aggregates; aggregate-associated C and C/N in farmland was significantly lower comparing with which in forestland; but aggregate-associated N did not show a distinct trend with the change of slope position and aggregate size. Aggregate-associated ammonium-N predominated in < 1 mm aggregates and which of forestland was higher than farmland; however, aggregate-associated nitrate-N mainly concentrated in >1 mm aggregates. PAD was affected obviously by land use type, aggregate size and slope position: PAD of farmland was significantly lower than forestland especially for >1 mm; PAD of 2-5 mm in upper, middle, middle to lower slope positions was obviously higher than other aggregate sizes in the same slope position.Aggregate-associated C and N did not present a significantly negative and positive relationship with PAD at significance level of 0.05; but C/N, ANR and ammonium-N were negatively correlated with PAD, and both ANR and C/N were in a lower level when PAD was greater than 40%.ConclusionDistribution of aggregate-associated C and N and aggregate stability between farmland and forestland are caused by the difference of nutrient restitution and soil erosion environment. Aggregate stability of black soil is determined by the interaction of aggregate-associated C and N, not the individual concentration of aggregate-associated C or N. The integrated indicators C/N and ANR, which are closely related to soil microorganism, have significant influence on soil aggregate disruption; moreover, ammonium-N used by microorganism preferentially in the soil environment also shows a significant correlation with aggregate stability. In the future, we should pay more attention to the relationship between microorganism activity and soil aggregate stability. -
近年来，以全球变暖为主要标志的气候变化给生态系统带来了巨大的压力[1-4]，陆地生态系统的生产力以及水分利用也受到了严重的影响[5-7]。而树木水分利用效率（water use efficiency，WUE）可以反映生态系统对全球气候变化的响应和描述不同生态系统的碳水循环关系 [8]。因此研究WUE不仅可以揭示生态系统对气候变化响应机理，还可以为区域气候变化对生态系统碳水耦合关系的影响提供科学的评价依据。
因研究目的、对象不同，WUE一般可分为生态系统、群体、叶片、以及个体4个方面。生态系统WUE是指植物消耗单位质量水分所固定的二氧化碳或生产干物质的量，因研究尺度、研究领域及获取数据手段的差异，不同的研究对生态系统WUE的计算往往存在一定的区别[9]。群体水平WUE是植物累计一段时期的干物质积累量与作物耗水量的比值，在植物群体水平上研究WUE时，其测算需要在大田试验中进行，工作复杂繁琐，计算出来的WUE不够准确[10]。叶片尺度上的WUE定义为单位通过蒸腾作用消耗单位水量时光合作用形成的有机物量[11]，在叶片尺度上研究植物的WUE，只能测得叶片尺度瞬时值，缺乏长时间尺度WUE的研究。而基于树木年轮稳定碳同位素的水分利用效率研究可以反映植物个体的长期WUE[12]，已有研究表明这一方法的可行性 [13-15]，并且该方法采样破坏性小，测定简单且不受时间和季节的限制[16]。目前国内外学者也进行了许多关于树木年轮稳定碳同位素方法研究植物WUE对气候变化响应的工作，如路伟伟等[17]研究北京山区油松（Pinus tabuliformis）WUE发现，区域整体WUE年际变化与温度呈显著负相关关系，与降水量呈正相关关系。Kannenberg等[18]研究美国西部灌木丛和森林的WUE发现，随着气候变得越来越干燥，美国西南部植物的WUE快速增加。Gagen等[19]利用树轮δ13C研究北芬诺斯干地亚植物的WUE发现，樟子松（Pinus sylvestris）WUE随着CO2浓度的升高而升高，但随着CO2浓度逐渐升高，樟子松WUE的升高出现了阈值。由于生态系统植被类型的多样性以及水分利用的有效性，气候因子对不同生态系统WUE的影响存在较大的差距，不同生态系统WUE的变化特征以及对气候变化的响应也不同。为此针对不同区域生态系统树木WUE进行研究是非常有必要的。
日本柳杉（Cryptomeria japonica）是庐山植被类型的重要树种组成部分，主要分布在海拔700 m以上的常绿针叶人工林，在水源涵养、净化大气环境、森林游憩等生态功能上发挥了重要的作用[20]，并且有着较高的生态效益和经济效益。因此本文以庐山日本柳杉作为研究对象，基于树木年轮和稳定碳同位素方法，分析气候变化背景下庐山日本柳杉WUE的变化特征以及庐山日本柳杉WUE与主要气候因子之间的关系，有助于庐山日本柳杉的碳汇能力评估以及提升森林的经营管理水平，对评估庐山日本柳杉的碳水耦合关系提供科学的理论支持。
1. 研究方法与数据来源
1.1 研究区概况
庐山位于江西省九江市（115°51′ ~ 116°07′E，29°30′ ~ 29°41′N），东偎鄱阳湖，海拔25.0 ~ 1 473.8 m，是著名的避暑胜地，气温适度，年平均最高温度32 ℃，最低气温−16.8 ℃，全年平均温11.4 ℃，全年平均降雨量1 917 mm，属湿润气候区（图1）。充足的水热条件及海拔的差异，共同造就了庐山丰富的植被资源[21]。高达40种植物在庐山地区首次被发现或者以庐山（牯岭）来命名。在植被类型区域划分上虽然属于亚热带常绿阔叶林，由于海拔高度具有较大的差异，因此植被类型上又极具多样性。海拔在1 000 m以下多为常绿或者是常绿落叶混交林，海拔在1 000 ~ 1 300 m为落叶阔叶林带[22]。生长在海拔较低处的植被由于受到人为活动的影响，植被破坏较为严重，而海拔较高处植被生长相对较好，对生态系统水源涵养的调节具有一定的作用[23]。
1.2 树木年轮的采集与碳同位素的测定
本研究根据国际树木年轮数据库的采样准则[24]，于2019年7月份在庐山3个不同海拔点取样（采样点1、2、3，见图1），海拔分别为800、950和1 150 m，平均树高为18.6、20.0和19.5 m，平均胸径为29.2、30.5和30.1 cm，林分密度为780 株/hm 2，郁闭度 0.6。每个样点取40个样芯。3个采样点的庐山日本柳杉生长良好，受人类活动以及自然灾害影响较小。
通过对日本柳杉样品进行交叉定年[25]，选取其树轮无生长异常、年轮宽度变化趋势较为一致、早晚材年轮界限明显的树芯用于稳定同位素的δ13C测定。具体步骤：用不锈钢刀由外而内逐年剥取全木，某些极窄年份的剥取在显微镜的辅助下进行精确切割，将不同树芯的相同年份合并在一起，置于80 ℃的烘箱中烘干至恒质量，枝剪剪碎之后的木材碎屑用高通量组织研磨仪研磨粉碎，过100目筛后装于5 mL离心管中，用于碳同位素测定分析，然后将样品粉末置于V（苯）∶V（乙醇） = 2∶1混合溶液中抽提24 h，除去树脂、树蜡及单宁类等有机物[26]。天平（精度0.001 mg）称取（0.20 ± 0.05） mg抽提好的样品包裹于锡杯中，在赛默飞稳定同位素质谱仪（delta v advantange）中将样品高温燃烧转换为气体，测量精度0.01‰，试验误差小于0.2‰[27]，测得稳定碳同位素的比值，即13C丰度（δ13C），每测定10个样品插入1个标样来测定仪器的稳定性（处理的树轮样品在江西农业大学森林培育重点实验室内进行分析测定），根据国际标准进行校正后得到日本柳杉稳定碳同位素δ13C值。
1.3 数据来源及处理
本文所用的气象数据来源于中国气象数据网（http://data.cma.cn/site/index.html）的庐山站日气象数据（庐山站的位置见图1），主要包括1968—2018年的平均气温（T）、最高温度、最低温度、相对湿度、降水量（W）、日照时数、平均风速（Ws）等日值数据。
参考1998年联合国粮农组织推荐的计算公式，利用日照时数数据作为输入参数计算庐山多年的太阳辐射（Rs），详细计算步骤可具体参阅文献[28]。为更准确量化庐山干旱严重程度，使用标准化降水蒸散指数（standardized precipitation evapotranspiration index，SPEI）、标准化降水指数（standardized precipitation index，SPI）以及湿度指数（wet index，WI）来表示庐山地区的干旱情况。SPEI和SPI的计算步骤具体参阅文献[29]，WI的计算计算步骤具体参阅文献[30]。利用树木年轮稳定碳同位素值、大气CO2中碳同位素比值和WUE关系来估算庐山日本柳杉1969—2018年的水分利用效率变化，根据Farquhar[13]推荐的公式计算δ13C。
δ13C=(13C/12C)sample−(13C/12C)PDB(13C/12C)PDB×1000‰ (1) 式中：（13C/12C）sample是测试样品13C/12C的摩尔比率； （13C/12C）PDB是标样13C/12C的摩尔比率，δ13C是树木年轮稳定碳同位素比值。
WUE可通过Δ与Ca值之间的关系计算获得[31]。
Δ=δ13Cp−δ13C1+δ13Cp (2) WUE=Ca×(b−Δ)1.6×(b−a) (3) 式中：δ13Cp是大气CO2碳同位素比值，Δ是指树木叶片和空气在光合作用期间同位素水平的差异δ13C的判别值；a是CO2扩散过程中的同位素分馏，取4.4‰，b为CO2羧化过程中的同位素分馏，取27.4‰，此外，系数1.6表示空气中水蒸气与CO2的扩散率之比。Ca为大气二氧化碳浓度，大气二氧化碳浓度由全球系统研究实验室（https://ngdc.noaa.gov/ftp.html）提供。
1.4 统计方法
本文采用一元线性回归模型分析1969—2018年庐山日本柳杉WUE的动态变化趋势[32-33]。并采用曼−肯德尔（M-K）突变分析方法，研究庐山日本柳杉WUE的突变年份[34]。
为预测庐山日本柳杉水分利用效率未来变化趋势，利用Hurst指数（H）进行计算，H是定量描述时间序列的自相似性以及长程依赖性的方法，常采用R/S分析法[35]。本研究中R表示50年日本柳杉WUE序列的极差，S为该时间段WUE的标准差，R/S为极差与标准差的比值，ln为对数函数，τ为偶数序列。当0.5 < H < 1，表明时间序列变化具有持续性的序列，未来的变量与过去的变量是一致的关系；H = 0.5，说明时间序列为相互独立，不存在相互影响的随机序列；0 < H < 0.5，表明时间序列变化具有反持续性的序列，未来的变量与过去的变量是完全相反的序列。其中R2越接近1，显著性越强。植物在生长过程中，气候因子对其生长通常表现出一定的“滞后效应”[36]，树木生长不仅受当年生长季气候因子的影响，还会受到上一年生长季后期气候因子变化的影响。因此，本研究选取上一年7月至当年12月共18个月份的气象数据进行相关性分析和多元回归分析。
2. 结果与分析
2.1 庐山地区各气候因子动态变化
通过对庐山1968—2018年各月气候因子的统计分析得到如图2的结果。从图2可以看出，干旱指数SPEI、WI以及SPI月变化基本一致。大致从2月到6月表现出下降的趋势，呈干旱趋势发展，SPEI和SPI在7月到10月波动较平缓，WI波动较大；干旱指数SPEI、WI以及SPI在12月呈上升趋势（图2a）。月平均气温从1月到7月呈上升趋势，最高气温出现在7月，为22.5 ℃；从8月到12月呈下降趋势，最低气温出现在1月，为0.3 ℃（图2b）。庐山地区的月平均太阳辐射和月平均风速最高出现在7月，为590.37 MJ/m2和4.94 m/s，月平均太阳辐射最低在1月，为285.72 MJ/m2，月平均风速最低在12月，为 3.53 m/s（图2c）。月平均降水量从1月到6月呈上升趋势，降水量最高出现在6月，为299.54 mm，最低出现在12月，为58.33 mm（图2d）。
![]()
图 2 庐山地区1968—2018 年月均SPI、SPEI、WI、T、Rs、Ws和W的变化及趋势Figure 2. Average changes and trends of SPI, SPEI, WI, T, Rs, Ws and W in Lushan Mountain area from 1968 to 20182.2 庐山日本柳杉年轮稳定碳同位素值以及大气二氧化碳同位素值变化规律
如图3a所示，本研究中树轮δ13C序列没有因幼龄效应而表现出明显的趋势变化。日本柳杉δ13C序列变动范围在−23.09‰ ~ −25.67‰，年平均−24.55‰。这与O’Leary等[37]的研究得出自然条件下生长的陆地C3植物碳稳定同位素变化范围在−22‰ ~ −34‰的结论是相符的。庐山日本柳杉的δ13C总体呈递减趋势，平均每10年下降0.34‰。线性回归方程的相关系数R2 = 0.59，δ13C下降趋势显著。但在2001年后有小幅度上升趋势，最大值达到−24.35‰，在2010年呈下降趋势发展，最小值为−25.67‰。庐山δ13Cp序列（1969—2018）变动范围在−7.29‰ ~ −8.93‰（图3b），在1969—2018年总体上是逐年递减的趋势，线性回归方程的相关系数R2 = 0.98，下降趋势显著，平均每10年下降0.32‰。
![]()
图 3 庐山日本柳杉年轮稳定碳同位素值和大气CO2碳同位素值变化及趋势Figure 3. Changes and trends of stable carbon isotope values (δ13C) and atmospheric CO2 carbon isotope values (δ13Cp) of Cryptomeria japonica2.3 庐山日本柳杉水分利用效率的变化特征
如图4a所示，1969—2018年间，庐山日本柳杉的WUE变动范围在91.06 ~ 118.89 μmol/mol，年平均值102.37 μmol/mol，总体呈现上升趋势，平均每10年上升5.2 μmol/mol，2013年出现最大值，为118.89 μmol/mol，最小值出现在1985年，为91.06 μmol/mol。为明确庐山日本柳杉WUE的年际变化特征，本研究对WUE进行了M-K突变分析和Hurst指数分析。从图4b的M-K突变分析可知：WUE呈低−高的趋势走向。1969—1987年，UF呈下降趋势，表明在这一时期，庐山日本柳杉WUE呈下降趋势，1988—2018年，UF呈上升趋势，表明庐山日本柳杉WUE呈上升趋势。庐山日本柳杉WUE的UF、UB交点出现在2004年（图4b），但庐山日本柳杉的WUE交点并不在±1.96临界线值范围内，没有通过0.05的检验，因此该年份的WUE交点上升不具有突变性。在双对数坐标系可以拟合得到年际庐山日本柳杉WUE的时间序列ln（R/S）与ln（τ/2）线性关系图（图4c）。由图4c可知，R/S趋势预测点之间存在良好的线性关系，R2 = 0.91，接近1，模型可行度高，线性拟合效果好。线性拟合的斜率为0.816 8，即Hurst指数估计值为0.816 8，表明庐山日本柳杉的WUE时间序列具有持续性关系，在未来变化中与过去存在较强的相关性，因此未来庐山日本柳杉的WUE呈现上升趋势。
![]()
图 4 庐山日本柳杉水分利用效率的变化特征及趋势WUE. 水分利用效率； UB.逆序时间序列变化；UF.顺序时间序列变化；R.极差；S.标准差；τ.偶数序列；ln.对数函数。WUE, water use efficiency; UB, inverse-order time series variation; UF, sequential time series variation; R, range; S, standard error; τ, even sequences; ln, logarithmic function.Figure 4. Variation characteristics and trends of water use efficiency of Cryptomeria japonica2.4 庐山日本柳杉水分利用效率与气候变化的相关分析以及多元回归模型解释
从庐山日本柳杉WUE序列与月气候变化的相关性分析（图5）可知：庐山日本柳杉WUE与各月气温大多呈正相关关系，与上年的7月、9—11月和当年的2—6月、9—11月呈显著正相关系，其中与上年的9月和当年的3—5月以及9月温度达到极显著正相关系；庐山日本柳杉WUE与各月太阳辐射大多呈负相关，但相关性不显著；WUE与干旱指数上年7—12月份、当年1—12月份干旱指数SPEI、SPI以及WI均呈负相关，但相关性不显著；庐山日本柳杉WUE与上年7—12月、当年1—12月风速均呈极显著负相关关系；WUE与上年7—12月至当年1—12月降水均无显著关系。在庐山日本柳杉的生长过程中， T和Ws在很大程度上驱动了日本柳杉水分利用效率的变化。通过多元线性回归方法将相关气候因子数值进行标准化来描述庐山日本柳杉水分利用效率与月气候因子之间的关系，回归过程选用逐步回归。得到庐山日本柳杉水分利用效率的回归方程。
![]()
图 5 庐山日本柳杉WUE与月气候变化的相关性P 为上年月份，C 为当年月份。**表示在置信度（双侧）为0.01时相关性显著；*表示在置信度（双侧）为0.05时相关性显著。P represents the month in previous year, C represents the month in current year. ** indicates a significant correlation at a confidence level (bilateral) of 0.01, * indicates that the correlation is significant at a confidence level (bilateral) of 0.05.Figure 5. Correlations between WUE of Cryptomeria japonica and monthly climate change in Lushan MountainWUE=58.349−3.733Ws1+1.487T3+1.690T9+1.667T−9−3.377Ws12(P<0.05，R2=0.65) 式中：WUE表示庐山日本柳杉水分利用效率，庐山日本柳杉水分利用效率主要与Ws1（当年1月风速）、T3（当年3月均温）、T9 （当年9月均温）、T−9（上年9月均温）、Ws12（当年12月风速） 相关。
3. 讨 论
3.1 干旱对日本柳杉水分利用效率的影响
水分利用效率是用来衡量植物个体或生态系统水平上碳水耦合关系的重要指标，能够揭示植物个体或生态系统WUE的变化特征[38]。WUE的影响因子随地区和尺度的不同而变化[39]。在本次研究中，本文选3个干旱指数进行庐山日本柳杉水分利用效率对干旱的响应验证发现，无论是SPEI、SPI还是WI，对庐山日本柳杉长期WUE的影响较小或者没有影响，均呈负相关关系（P > 0.05）。Song等[40]对西南亚热带常绿原生林水分利用效率研究表明亚热带常绿森林在干旱年的WUE呈增加趋势。一般认为，不同的生态系统WUE对干旱的响应是不同的。Liu等[41]对中国陆地生态系统水分利用效率与干旱响应研究发现在中国南方，中度和极端干旱会导致WUE下降，严重干旱往往导致WUE略有上升，而在东北和内蒙古中部地区轻微干旱通常会导致WUE增加。杜晓铮等[7]在研究陆地生态系统水分利用效率对气候变化的响应研究综述表明，在不同的气候区域干旱对WUE的影响有所差异，不同生态系统WUE对干旱的响应不同。有研究表明，干旱对植物的水分利用效率存在一定的滞后影响[42]。在本研究中发现，干旱对庐山日本柳杉WUE没有产生滞后影响。在湿润区，干旱发生的同时由于太阳辐射的射入，可能会加快植物的生长速度，同时蒸腾作用也会损失一部分的水分，造成植物WUE呈现出下降的趋势。在本次研究我们也发现，降水与 WUE呈负相关关系（P > 0.05），但相关性不显著。这可能是由于庐山地处湿润区，降水丰富，导致植物叶片气孔导度增大、蒸腾速率增强，造成了庐山日本柳杉水分利用效率的降低。
3.2 气温对日本柳杉水分利用效率的影响
温度对于植物水分利用的影响较为复杂。我们发现，气温对庐山日本柳杉的水分利用效率存在一定的滞后影响，上年的9月气温与庐山日本柳杉WUE呈显著正相关关系，当年的9月气温与庐山日本柳杉WUE也呈显著正相关关系。可能是由于气温适中，植物酶活性增强，从而影响植物水分利用。周佳等[43]利用树木年轮研究河南民权与陕西白水刺槐（Robinia pseudoacacia）水分利用效率对气候的响应表明，两地刺槐WUE与平均温度存在显著正效应关系。一定范围内，温度升高，叶片的气孔导度增加，植物的净光合速率大于蒸腾速率，造成植物WUE升高。崔茜琳等[44]利用MODIS数据研究青藏高原植被水分利用效率，青藏高原植被的WUE 与气温呈显著正相关关系，WUE的敏感性随气温的升高而增加。仇宽彪等[45]对陕西省各植被类型WUE研究时发现，在温度低于11 ℃的地区WUE和温度呈现显著正相关（P < 0.01），而高于11 ℃的地区两者之间没有显著关系（P > 0.05）。在不同的气候区， WUE对温度的响应存在着较大的差异。
3.3 风速对日本柳杉水分利用效率的影响
风会影响植物周围的环境，如降低空气的相对湿度和温度，并通过植物叶片遮挡影响太阳辐射的射入[46]，同时会加速蒸腾速率而相应地降低植物的温度[47]。在不同的生态系统中，水分利用效率对Ws的响应幅度可能不同，甚至可能有完全不同的方向。王云霓[48]在研究宁夏六盘山典型树种水分利用效率中发现华北落叶松（Larix principis-rupprechtii）个体的WUE与Ws是正相关关系（P > 0.05）。在本次研究中发现，庐山日本柳杉WUE与Ws显著负相关关系，特别是当年的1月份和12月份的风速（P < 0.01），对庐山日本柳杉WUE产生了重要的影响。风速越大，植物的蒸腾作用会越强，植物体周围的空气湿度会降低，促使了植物体内的水分更快散失到大气中，造成植物WUE下降。在不同的地区，风速对于植物的WUE影响是有所差异的。庐山自然保护区风速呈现下降的趋势，平均每10年下降0.37 m/s，是日本柳杉WUE年际变化重要的影响因素之一。
3.4 太阳辐射对日本柳杉水分利用效率的影响
太阳辐射是植物光合作用所必须的重要元素之一，对植物WUE的变化也会产生重要影响。如徐晓桃[49]研究黄河源区植被的水分利用效率时发现，黄河源区的植被WUE与Rs是正相关关系。王云霓[48]研究宁夏六盘山典型树种华北落叶松时也发现WUE与Rs呈正相关关系，并且通过了显著性检验（P < 0.01）。一般认为，光照越强，植物的水分利用越高。但庐山日本柳杉WUE在本研究中与Rs未产生显著影响，这与他们的结果不同。一方面可能是由于太阳辐射增加时，植物的光合速率增大，但同时植物的蒸腾也在增大，WUE的变化难以确定。作为碳水循环的综合衡量指标，无论是植物个体的WUE或者是整个生态系统WUE，WUE的变化都是由多个相互作用的物理和生物过程决定，各个气候因子对WUE的变化具有相互协同作用。在评价树木WUE变化对气候因子的响应时，各个气候因子之间综合协同对WUE产生的影响还需进一步探究。
4. 结 论
本研究表明，1969—2018年期间，庐山日本柳杉树轮δ13C值呈下降趋势，WUE变化呈显著升高趋势。温度和风速是影响庐山日本柳杉WUE变化的主要气候因子，WUE与温度具有显著正相关关系，与风速具有显著负相关关系，并且温度对庐山日本柳杉WUE具有一定的滞后影响。而降水、太阳辐射以及干旱指数（SPI、SPEI、WI）对庐山日本柳杉年际WUE的变化未产生显著影响。本研究结果可为在全球气候变化下对庐山日本柳杉生长产生的影响提供理论依据，并且有助于庐山日本柳杉碳汇能力的评估以及提升杉木林的经营管理水平，同时为该地区杉木林生态系统的维护提供理论指导。
-
![]()
图 1 不同粒径和坡位处团聚体破坏率
不同大写字母表示相同坡位不同粒径之间差异显著；不同小写字母表示相同粒径不同坡位之间差异显著，P < 0.05。
Figure 1. PAD of water stable-soil aggregates for varied particle sizes and slope positions
Different capital letters represent significant difference among varied particle sizes for the same slope position; different lowercase letters represent significant difference among varied slope positions for the same particle size, P < 0.05.
表 1 不同粒径和坡位处水稳性团聚体碳氮含量
Table 1 Soil carbon and nitrogen content of water stable-soil aggregates for varied particle sizes and slope positions
项目Item 耕地Farmland 林地Forestland 总碳Total C/
(g·kg-1)总氮Total N/
(g·kg-1)C/N 总碳Total C/
(g·kg-1)总氮Total N/
(g·kg-1)C/N 2~5 mm 25.01Ba 1.95Aa 12.88Bab 40.25Ac 1.19Bb 39.02Aa 1~2 mm 24.30Bb 1.91Aa 12.79Bb 40.37Ac 1.67Ba 28.84Ab 0.5~1 mm 24.66Bb 1.93Aa 12.90Bab 46.82Ab 1.58Ba 26.99Ab 0.25~0.5 mm 25.28Ba 1.92Aa 13.51Ba 54.19Aa 2.01Aa 28.26Ab 坡上Upper slope 23.66Bc 1.87Ab 12.80Bb 30.46Ad 1.37Bc 23.95Ac 坡中上Upper to middle slope 23.50Bcd 1.92Aab 12.32Bb 35.94Ac 1.23Bc 32.42Ab 坡中Middle slope 23.09Bd 1.90Aab 12.18Bb 51.64Ab 2.16Ab 26.01Ac 坡中下Middle to lower slope 25.15Bb 2.05Aa 12.36Bb 37.28Ac 1.06Bc 42.03Aa 坡下Lower slope 28.66Ba 1.90Bab 15.44Ba 68.41Aa 2.69Aa 27.95Ac 注：不同大写字母表示相同粒径或坡位、不同土地利用类型之间差异显著，P < 0.05，n=20；不同小写字母表示相同土地利用类型、不同粒径或粒径之间差异显著，P < 0.05，n=16。表 2同此。Notes: different capital letters represent significant difference within the same row(i.e. the same in particle size or position, but different in land-use type), P < 0.05, n=20; different lowercase letters represent significant difference within the same column (i.e. the same in land-use type, but different in particle size or position), P < 0.05, n=16. The same in Tab. 2. 
下载: 导出CSV
表 2 不同粒径和坡位处水稳性团聚体铵态氮与硝态氮含量
Table 2 Soil ammonium-N and nitrate-N content of water stable-soil aggregates for varied particle sizes and slope positions
项目Item 耕地Farmland 林地Forestland 铵态氮
Ammonium-N/ (mg·kg-1)硝态氮
Nitrate-N/ (mg·kg-1)ANR 铵态氮
Ammonium-N/ (mg·kg-1)硝态氮
Nitrate-N/ (mg·kg-1)ANR 2~5 mm 27.54Bb 4.08Ba 8.20Ac 41.38Ac 8.84Aa 9.93Ad 1~2 mm 27.42Bb 3.68Bab 8.72Bbc 40.70Ac 4.77Ab 11.78Ac 0.5~1 mm 31.18Ba 3.86Aab 9.88Bab 51.25Ab 3.73Bc 17.09Ab 0.25~0.5 mm 32.19Ba 3.19Bb 11.48Ba 64.84Aa 3.65Ac 20.90Aa 坡上Upper slope 35.16Aa 4.92Aa 10.47Ba 35.07Ad 2.38Be 15.35Ab 坡中上Upper to middle slope 25.38Bc 2.92Bc 9.60Babc 65.08Ab 3.98Ab 18.76Aa 坡中Middle slope 29.95Bb 3.03Bc 10.35Bab 68.72Aa 3.55Ac 19.65Aa 坡中下Middle to lower slope 32.83Ba 4.21Aab 9.48Babc 52.68Ac 2.85Bd 19.39Aa 坡下Lower slope 27.41Bc 3.58Bbc 8.33Bc 30.36Ae 12.33Aa 3.49Bc 注：ANR表示铵态氮与硝态氮的比值。下同。Notes: ANR represents the ratio of ammonium-N to nitrate-N. The same as below.
下载: 导出CSV
表 3 团聚体破坏率与团聚体碳氮含量的相关系数(r)
Table 3 Correlation coefficient (r) for PAD with aggregate-associated carbon and/or nitrogen content
项目Item 总碳Total C 总氮Total N C/N 铵态氮Ammonium-N 硝态氮Nitrate-N ANR PAD -0.282 0.144 -0.357* -0.386* -0.108 -0.321* 注：*表示显著相关，P < 0.05，n=20。Notes:* represents significant correlation at P < 0.05 level, n=20.
下载: 导出CSV
表 4 不同处理及其交互作用对土壤团聚体碳氮含量和稳定性的影响(P值)
Table 4 Effects of different factors and their interactions on carbon and nitrogen content and stability of aggregates (P value)
影响因子
Influencing factor总碳
Total C总氮
Total N铵态氮
Ammonium-N硝态氮
Nitrate-NC/N ANR PAD J < 0.01 0.283 < 0.01 < 0.01 0.029 < 0.01 < 0.01 P < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 L < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 P×J < 0.01 0.140 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 P×L < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 < 0.01 L×J < 0.01 0.168 < 0.01 < 0.01 0.200 < 0.01 < 0.01 P×J×L < 0.01 0.229 0.042 < 0.01 < 0.01 < 0.01 < 0.01 注：J表示粒径，P表示坡位，L表示土地利用类型；×表示交互作用。Notes: J represents particle size, P represents slope position, L represents land-use type, × presents interaction of different factors.
下载: 导出CSV
-
[1] Tisdall J M, Oades J M. Organic matter and water-stable aggregates in soil[J]. Journal of Soil Science, 1982, 33(2): 141-163. doi: 10.1111/ejs.1982.33.issue-2
[2] 张明奎, 何振立, 陈国潮, 等.利用方式对红壤水稳定性团聚体形成的影响[J].土壤学报, 1997, 34(4): 359-365. doi: 10.3321/j.issn:0564-3929.1997.04.002 Zhang M K, He Z L, Chen G C, et al. Formation of water-stable aggregates in red soils as affected by land use[J]. Acta Pedologica Sinica, 1997, 34(4): 359-365. doi: 10.3321/j.issn:0564-3929.1997.04.002
[3] Srivastava P, Raghubanshi A S, Singh R, et al. Soil carbon efflux and sequestration as a function of relative availability of inorganic N pools in dry tropical agroecosystem[J]. Applied Soil Ecology, 2015, 96:1-6. doi: 10.1016/j.apsoil.2015.04.019
[4] Srivastava P, Singh P K, Singh R, et al. Relative availability of inorganic N-pools shifts under land use change: an unexplored variable in soil carbon dynamics[J]. Ecological Indicators, 2016, 64: 228-236. doi: 10.1016/j.ecolind.2015.12.043
[5] Chrenková K, Mataix-Solera J, Dlapa P, et al. Long-term changes in soil aggregation comparing forest and agricultural land use in different Mediterranean soil types[J]. Geoderma, 2014, 235-236(1): 290-299. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=e9d3f60ebad1ddff2a1cf727ca7d5096
[6] An S S, Mentler A, Mayer H, et al. Soil aggregation, aggregate stability, organic carbon and nitrogen in different soil aggregate fractions under forest and shrub vegetation on the Loess Plateau, China[J]. Catena, 2010, 81(3): 226-233. doi: 10.1016/j.catena.2010.04.002
[7] 彭畅.长期施肥条件下黑土有机碳库和氮库变化研究[D].北京: 中国农业科学院, 2006. http://cdmd.cnki.com.cn/article/cdmd-82101-2006110723.htm Peng C. Change of organic carbon and nitrogen pool of black soil under long-term fertilization[D]. Beijing: Chinese Academy of Agricultural Sciences, 2006. http://cdmd.cnki.com.cn/article/cdmd-82101-2006110723.htm
[8] 李鉴霖, 江长胜, 郝庆菊.土地利用方式对缙云山土壤团聚体稳定性及其有机碳的影响[J].环境科学, 2004, 35(12): 4695-4704. http://d.old.wanfangdata.com.cn/Periodical/hjkx201412037 Li J L, Jiang C S, Hao Q J. Impact of land use type on stability and organic carbon of soil aggregates in Jinyun Mountain[J]. Environment Science, 2004, 35(12): 4695-4704. http://d.old.wanfangdata.com.cn/Periodical/hjkx201412037
[9] 刘欣, 黄运湘, 袁红, 等.植被类型与坡位对喀斯特土壤氮转化速率的影响[J].生态学报, 2016, 36(9) : 2578-2587. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stxb201609016 Liu X, Huang Y X, Yuan H, et al. Effects of vegetation type and slope position on soil nitrogen transformation rate in karst regions[J]. Acta Ecologica Sinica, 2016, 36(9):2578-2587. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stxb201609016
[10] Bennie J, Mark O H, Baxter R, et al. Influence of slope and aspect on long-term vegetation change in British chalk grass[J]. Journal of Ecology, 2006, 94(2): 355-368. doi: 10.1111/jec.2006.94.issue-2
[11] Jackson L E, Schimel J P, Firestone M K. Short-term partitioning of ammonium and nitrate between plants and microbes in an annual grassland[J]. Soil Biology and Biochemistry, 1989, 21(3): 409-415. doi: 10.1016/0038-0717(89)90152-1
[12] 隋跃宇, 张兴义, 张少良, 等.黑龙江典型县域农田黑土土壤有机质现状分析[J].土壤通报, 2008, 39(1): 186-188. doi: 10.3321/j.issn:0564-3945.2008.01.035 Sui Y Y, Zhang X Y, Zhang S L, et al. Soil organic matter actuality of the black soil farmland in Heilongjiang counties[J]. Chinese Journal of Soil Science, 2008, 39(1): 186-188. doi: 10.3321/j.issn:0564-3945.2008.01.035
[13] 向艳文, 郑圣先, 廖育林, 等.长期施肥对红壤水稻土水稳性团聚体有机碳、氮分布与储量的影响[J].中国农业科学, 2009, 42(7): 2415-2424. http://d.old.wanfangdata.com.cn/Periodical/zgnykx200907021 Xiang Y W, Zheng S X, Liao Y L, et al. Effects of long-term fertilization on distribution and storage of organic carbon and nitrogen in water-stable aggregates of red paddy soil[J]. Scientia Agricultura Sinica, 2009, 42(7): 2415-2424. http://d.old.wanfangdata.com.cn/Periodical/zgnykx200907021
[14] Ashagrie Y, Zech W, Guggenberger G, et al. Soil aggregation, and total and particulate organic matter following conversion of native forests to continuous cultivation in Ethiopia[J]. Soil and Tillage Research, 2007, 94(1): 101-108. doi: 10.1016/j.still.2006.07.005
[15] 樊红柱, 张建辉, 王勇, 等.川北山区坡耕地侵蚀耕作对土壤团聚体碳的影响[J].农业机械学报, 2015, 46(11): 157-164. doi: 10.6041/j.issn.1000-1298.2015.11.022 Fan H Z, Zhang J H, Wang Y, et al. Tillage erosion impacts on soil aggregate associated carbon in mountainous region slope farmland of northern Sichuan[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(11): 157-164. doi: 10.6041/j.issn.1000-1298.2015.11.022
[16] 王勇, 张建辉, 李富程.耕作侵蚀对坡耕地土壤水稳定性团聚体和水分特征的影响[J].水土保持学报, 2015, 29(1): 180-185. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqsystbcxb201501036 Wang Y, Zhang J H, Li F C. Tillage erosion influence on soil water-stable aggregate and moisture properties on the hillslope[J]. Journal of Soil and Water Conservation, 2015, 29(1): 180-185. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trqsystbcxb201501036
[17] 陈晓芬, 李忠佩, 刘明, 等.不同施肥处理对红壤水稻土团聚体有机碳、氮分布和微生物生物量的影响[J].中国农业科学, 2013, 46(5): 950-960. doi: 10.3864/j.issn.0578-1752.2013.05.010 Chen X F, Li Z P, Liu M, et al. Effects of different fertilizations on organic carbon and nitrogen contents in water-stable aggregates and microbial biomass content in paddy soil of subtropical China[J]. Scientia Agricultura Sinica, 2013, 46(5): 950-960. doi: 10.3864/j.issn.0578-1752.2013.05.010
[18] 刘忠良, 宇万太, 周桦, 等.不同有机厩肥输入量对土壤团聚体有机碳组分的影响[J].土壤学报, 2011, 48(6): 1149-1157. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trxb201106006 Liu Z L, Yu W T, Zhou H, et al. Effect of application rate of barnyard manure on organic fraction of soil aggregates[J]. Acta Pedologica Sinica, 2011, 48(6): 1149-1157. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=trxb201106006
[19] 窦森, 郝翔翔.黑土团聚体与颗粒中碳、氮含量及腐殖质组成的比较[J].中国农业科学, 2013, 46(5): 970-977. doi: 10.3864/j.issn.0578-1752.2013.05.012 Dou S, Hao X X. Comparison of carbon, nitrogen contents and humus compositions in the aggregates and particles of black soil[J]. Scientia Agricultura Sinica, 2013, 46(5): 970-977. doi: 10.3864/j.issn.0578-1752.2013.05.012
[20] 徐文静, 丛耀辉, 张玉玲, 等.黑土区水稻土水稳性团聚体有机碳及其颗粒有机碳的分布特征[J].水土保持学报, 2016, 30(4): 210-215. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201604036 Xu W J, Cong Y H, Zhang Y L, et al. Distribution of organic carbon and particulate organic carbon in water-stable aggregates of paddy soil in black soil area[J]. Journal of Soil and Water Conservation, 2016, 30(4): 210-215. http://d.old.wanfangdata.com.cn/Periodical/trqsystbcxb201604036
[21] Li G Y, Fan H M. Effect of freeze-thaw on water stability of aggregates in a black soil of northeast China[J]. Pedosphere, 2016, 24(2): 285-290. http://cn.bing.com/academic/profile?id=72fa75e55251c5c61881871a72ac1af1&encoded=0&v=paper_preview&mkt=zh-cn
[22] 李露, 曹伟, 赵鹏志, 等.黑土区坡耕地与坡林地表土团聚体有机碳库差异[J].北京林业大学学报, 2016, 38(12): 46-53. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bjlydxxb201612006 Li L, Cao W, Zhao P Z, et al. Differences of soil organic carbon storage in surface soil aggregates between sloping farm land and sloping forest areas in black soil region[J]. Journal of Beijing Forestry University, 2016, 38(12): 46-53. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bjlydxxb201612006
[23] 赵鹏志, 陈祥伟, 王恩姮.东北黑土区典型坡面耕作侵蚀定量分析[J].农业工程学报, 2016, 32(12): 151-157. doi: 10.11975/j.issn.1002-6819.2016.12.022 Zhao P Z, Chen X W, Wang E H. Quantitative assessment of tillage erosion on typical sloping field in black soil area of northeast China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(12): 151-157. doi: 10.11975/j.issn.1002-6819.2016.12.022
[24] 刘孝义.土壤物理及土壤改良研究法[M].上海:上海科学技术出版社, 1981. Liu X Y. Soil physics and the methods of soil improvement[M]. Shanghai: Shanghai Science and Technology Press, 1981.
[25] 陈立新.土壤实验实习教程[M].哈尔滨:东北林业大学出版社, 2005. Chen L X. Soil experiment and practice course[M]. Harbin: Northeast Forestry University Press, 2005.
[26] 许中坚, 刘广深, 喻佳栋, 等.模拟酸雨对红壤结构体及其胶结物影响的实验研究[J].水土保持学报, 2002, 16(3): 9-11. doi: 10.3321/j.issn:1009-2242.2002.03.003 Xu Z J, Liu G S, Yu J D, et al. Effects of acid rain on aggregate and its cemedins of red soil by simulated experiments[J]. Journal of Soil and Water Conservation, 2002, 16(3): 9-11. doi: 10.3321/j.issn:1009-2242.2002.03.003
[27] 刘震, 徐明岗, 段英华, 等.长期不同施肥下黑土和红壤团聚体氮库分布特征[J].植物营养与肥料学报, 2013, 19(6) : 1386-1392. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zwyyyflxb201306013 Liu Z, Xu M G, Duan H Y, et al. Distribution of nitrogen in aggregates of black soil and red soil under long-term fertilization[J]. Journal of Plant Nutrition and Fertilizer, 2013, 19(6): 1386-1392. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zwyyyflxb201306013
[28] Six J, Elliott E T, Paustian K, et al. Aggregation and soil organic matter accumulation in cultivated and native grassland soils[J]. Soil Science Society of America Journal, 1998, 62:1367-1377. doi: 10.2136/sssaj1998.03615995006200050032x
[29] 葛翠萍, 赵军, 王秀峰, 等.东北黑土区坡耕地地形因子对土壤水分和容重的影响[J].水土保持学报, 2008, 28(6): 16-19. http://www.cnki.com.cn/Article/CJFDTOTAL-STTB200806005.htm Ge C P, Zhao J, Wang X F, et al. Influence of topographic factors on soil water and bulk density in a typical slope land in the black soil area of northeast China[J]. Journal of Soil and Water Conservation, 2008, 28(6): 16-19. http://www.cnki.com.cn/Article/CJFDTOTAL-STTB200806005.htm
[30] 李银坤, 陈敏鹏, 梅旭荣, 等.土壤水分和氮添加对华北平原高产农田有机碳矿化的影响[J].生态学报, 2014, 34(14) : 4037-4046. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stxb201414028 Li Y K, Chen M P, Mei X R, et al. Effects of soil moisture and nitrogen addition on organic carbon mineralization in a high-yield cropland soil of the North China Plain[J]. Acta Ecologica Sinica, 2014, 34(14): 4037-4046. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stxb201414028
[31] 刘晶鑫, 迟凤琴, 许修宏, 等.长期施肥对农田黑土微生物群落功能多样性的影响[J].应用生态学报, 2015, 26(10): 3066-3072. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb201510019 Liu J X, Chi F Q, Xu X H, et al. Effect of long-term fertilization on microbial community functional diversity in black soil[J]. Chinese Journal of Applied Ecology, 2015, 26(10): 3066-3072. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb201510019
-
期刊类型引用(3)
1. 陈彩云. 天然染料在鲜茧丝染色中的应用. 轻纺工业与技术. 2023(02): 145-147 . 
百度学术
2. 杨雨桐，张卿硕，符韵林，孙静. 巴里黄檀心材色素的提取及其抑菌活性. 东北林业大学学报. 2020(10): 100-103+119 . 百度学术
3. 王敬贤. 木材染色技术研究进展. 林业科技. 2020(06): 42-47 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 1922
- HTML全文浏览量: 512
- PDF下载量: 51
- 被引次数: 7
