Leaf δ¹³C of the Main Afforestation Tree Species in Horqin Sandy Land

研究区位于科尔沁沙地南缘，辽宁省阜新市彰武县章古台镇辽宁省沙地治理与利用研究所试验基地（42°42′ N, 122°22′ E; 平均海拔226 m）。该区属于半干旱气候区。年均气温7.2 ℃，年降水量479 mm（1954—2019），年蒸发量1 300～1 800 mm^[4]。冬春两季风大而持续时间长，风速为4.15～5.10 m·s⁻¹ ，无霜期150 d。土壤为风沙土，pH 值6.7左右，养分元素N、P 等含量低^[22]。按地貌特征、地形起伏状况、植物覆盖和土壤性质等因素将章古台沙地划分为4种类型：流动与半流动沙丘、固定沙丘、平缓沙地和风蚀洼地^[24]。土壤盐分含量、土壤质地和土壤结构在研究区呈均匀分布。该地区主要树种有樟子松、杨树、油松和白榆等。林下植被以1年生草本植物为主，主要为禾本科（Graminaea）、藜科（Chenopodiaceae）、菊科（Compositae）、莎草科（Cyperaceae）^[22]。研究区平均地下水位为5.0 m 左右。2011年，研究期间（6—10月）降水量为361.1 mm，占常年同期降水量（384.7 mm）的94%；监测期间气温平均值为19.1 ℃，略高于同期平均气温（18.7 ℃），平均相对湿度为69%，平均饱和水汽压差为0.66 kPa，光合有效辐射为400.5 μmol·m⁻² ·s⁻¹，平均风速为2.0 m·s⁻¹（图1）

Figure 1.  Seasonal variation in the climatic variables for the study region in 2011

在生长季初期的5月，在对试验区进行全面踏查的基础上，选择31年生片状樟子松林分、12年生彰武小钻杨林带、40年生片状油松林分和42年生散生白榆（表1）。为了减小林分密度对叶片δ¹³C值的影响，在樟子松、油松和杨树林分内随机选择5株优势木，同时随机选取5株散生白榆，共计20株，挂牌作为研究样树。白榆、油松和樟子松林分属于中龄林，杨树林带接近中龄林（中华人民共和国林业行业标准：主要树种龄级与龄组划分）。于2011年生长季的6 — 10月每月中旬取1次叶片和土壤样品，整个生长季共计5次。

为了避免取样高度对叶片δ¹³C影响，不同树种每次取样高度保持基本一致（表1）。每次取样集中于上午8:00—10:00，用高枝剪采集每株样树树冠中部、向阳部位的叶片。樟子松和油松每株样树取当年生叶各100束；杨树和白榆每株样树取30片叶。将取得的叶立即装进信封内，迅速带回实验室。在实验室内用蒸馏水清洗叶片样品，在105 ℃烘箱中杀青20 min，然后在70 ℃烘箱中烘48 h，使样品完全干燥。利用粉碎机（FW100, 天津泰斯特）将叶样品粉碎，过100目筛使样品均一化，以备测定δ¹³C值。

叶片δ¹³C值在中国科学院沈阳应用生态研究所农产品安全与环境质量检测中心用稳定性同位素比率质谱仪（DELTA ^plus XP , America）进行测定，测试误差小于0.15 ‰，以PDB（Pee Dee Belemnite）为标准，根据下面公式（1）进行计算：

其中, δ¹³C表示样品¹³C/¹²C与标准样品偏离的千分率，(¹³C/¹²C)_PDB表示南卡罗来纳州白碚石(Pee DeeBelemnite)中的¹³C/¹²C。

每次叶片样品取样的当天，在每株样树树冠下选择1个样点，用土钻取土。取土深度为0～20，20～40 和40～60 cm。用烘干法测定土壤含水量。

用SPSS19.0软件进行双因素方差分析，新复极差法（Duncan）比较树种和取样日期对土壤含水量和当年生叶片δ¹³C值影响。采用单因素方差分析比较同一月份不同土壤深度土壤含水量和同一月份不同树种叶片δ¹³C值的差异。利用Pearson相关性分析方法分析不同树种叶δ¹³C逐月变化量（将当月与前月叶片δ¹³C值相减）与水热因子关系。显著性水平设定为α=0.05。

不同树种土壤含水量（0～60 cm）随着季节变化呈先增加后降低趋势（图2）。樟子松和彰武小钻杨土壤含水量按月份排序均为：8月和7月>10月、9月和6月（图2A 和2B, P<0.05）；油松土壤含水量按月份排序为：8月>10月和9月>7月和6月（图2C, P<0.05）；白榆土壤含水量按月份排序为：8月和7月>10月>6月和9月（图2D, P<0.05）。双因素方差分析结果表明（表2），不同树种和月份间土壤含水量存在显著性差异。多重比较结果表明（表3），樟子松和白榆的土壤含水量显著高于彰武小钻杨和油松；土壤含水量按月份排序为：8月>7月> 10月>6月和9月。

Figure 2.  Seasonal variations in soil water content at 0-60 cm depth for different tree species

不同月份樟子松叶片δ¹³C值介于−27.44‰～−26.36‰，平均值为−26.92‰；彰武小钻杨叶片δ¹³C值介于−27.87‰～−26.97‰，平均值为−27.40‰；油松叶片δ¹³C值介于−28.00‰～−26.72‰，平均值为−27.57‰；白榆叶片δ¹³C值介于−30.02‰～−28.45‰，平均值为−29.23‰（图3A）。樟子松、油松和白榆叶δ¹³C值月变化量（将当月与前月叶片δ¹³C值相减）生长季几乎都为负值，而彰武小钻杨叶δ¹³C值月变化7月和8月为正值，9月和10月为负值（图3B）。樟子松、油松和白榆叶δ¹³C值不同月份间差异都显著（P<0.05），排序都为：6月和7月> 8月、9月和10月；彰武小钻杨叶δ¹³C值不同月份间差异显著（P<0.05），排序为：8、7和9月>6月和10月。同一月份不同树种间叶δ¹³C值存在显著性显著（图3A, P<0.05），但是这种差异随月份变化而不同。双因素方差分析结果表明（表2），不同树种和月份间叶δ¹³C值（整个生长季平均值）存在显著性差异（P<0.05）。多重比较结果表明（表3），樟子松叶δ¹³C值显著高于彰武小钻杨和油松，而彰武小钻杨和油松显著高于白榆；6月和7月叶δ¹³C显著高于8月和9月，显著高于10月。

Figure 3.  Seasonal changes in the leaf δ¹³C values and leaf δ¹³C monthly change values for different tree species

将各月δ¹³C变化量与各月气象要素和土壤含水量进行相关分析。樟子松叶片δ¹³C变化量与降水量、相对湿度和土壤含水量呈负相关关系，与温度、饱和水汽压差和光合有效辐射呈正相关关系，相关关系都不显著（表3）。彰武小钻杨叶片δ¹³C变化量与气温呈显著正相关关系（P<0.05），与其他水热因子也呈正相关关系，相关关系不显著（表3）。油松叶片δ¹³C变化量与降水量、气温和相对湿度呈显著负相关关系（P<0.05），与其他水热因子也呈负相关关系，相关关系不显著（表3）。白榆叶片δ¹³C变化量与水热因子呈负相关关系，相关关系都不显著（表3）。

本研究发现，樟子松、油松和白榆叶片δ¹³C值随着季节性变化呈下降趋势，而彰武小钻杨叶片δ¹³C值随着季节变化呈先升高后降低趋势，表明樟子松、油松和白榆水分利用效率随着季节变化逐渐下降，而彰武小钻杨水分利用效率随着季节变化呈先增加后降低（图3）。这与其他研究结果基本一致。例如，郭树江等^[30]研究发现，民勤3种主要灌木植物水分利用效率随着季节变化呈降低趋势，主要与平均气温、相对湿度、蒸发量等环境因子季节性变化有关。曹生奎等^[31]研究表明，极端干旱区胡杨（Populus euphratica Oliv.）生长季水分利用效率逐渐下降，主要是因为气温升高、土壤含水量减小和地下水位埋深下降导致的气孔导度的变化引起的。方晓娟等^[32]研究发现随着季节变化，毛白杨（Populus tomentosa）杂种无性系水分利用效率呈先增加后降低趋势。本研究中不同树种水分利用效率随季节变化有所不同，这可能与不同树种对外界水热因子季节性变化响应程度不同有关^[33]。樟子松、油松和白榆叶片δ¹³C变化量与降水量、相对湿度和土壤含水量呈负相关关系，尤其是油松（表3），表明当土壤含水量降低, 空气湿度降低以及降水量不足时, 树木为了减少水分的蒸发, 降低气孔导度，提高水分利用效率^[33]。另外，油松叶片δ¹³C变化量与温度呈显著负相关关系，表明随着温度升高，油松水分利用效率降低。这与路伟伟等^[34]研究结果一致，即随着气温升高，北京市密云东部山区油松水分利用效率降低。一般情况下，干旱半干旱地区降水对树木水分利用效率的影响程度高于气温^[10]。6月份降水量和土壤含水量相对较低（图1），但是温度相对较高，土壤蒸发也强烈，因此，6月份，樟子松、油松和白榆为减少水分丧失，维持较低的气孔导度，从而具有较高的水分利用效率。6月份以后，由于降水量的增加，土壤含水量显著增加，能被树木利用的水分增加，树木气孔导度增加，蒸腾强度增加，导致树木在生长旺盛期的水分利用效率降低。8月之后，随着降水量减少和温度降低（图1），尽管土壤含水量降低，但是树木蒸腾需求量逐渐降低。因此，树木水分利用效率仍呈降低趋势^[35]。另外，彰武小钻杨叶片δ¹³C变化量与温度呈显著正相关关系，表明随着温度升高，促使彰武小钻杨气孔导度降低，水分利用效率增加。这与Sun 等^[36]研究结果一致，随着温度升高，张北地区小叶杨（Populus simonii Carr.）水分利用效率显著升高。6—8月，随着温度和光照强度增加，彰武小钻杨光合速率和蒸腾能力逐渐增加^[32]，尽管降水量增加，但是土壤含水量远远不能满足其蒸腾耗水的需要，彰武小钻杨通过降低气孔导度，从而提高水分利用效率，维持树木生长。8月份后，随着光照强度和温度逐渐降低，光合速率和蒸腾能力逐渐降低，但是光合速率的降低程度大于蒸腾速率^[37]，从而导致彰武小钻杨水分利用效率呈降低趋势。

另外，彰武小钻杨叶片δ¹³C变化量与相对湿度和饱和水汽压差呈正相关关系，而油松和白榆叶片δ¹³C变化量与相对湿度和饱和水汽压差呈负相关关系（表3），这表明不同树种叶片δ¹³C对空气相对湿度和饱和水汽压差变化具有不同的响应。这可能是因为随着空气相对湿度降低，彰武小钻杨叶片羧化效率降低，从而导致叶片δ¹³C降低；伴随着饱和水汽压差增加，气孔导度降低，引起彰武小钻杨叶片δ¹³C升高。然而，随着空气相对湿度降低，油松和白榆气孔导度和胞间CO₂浓度降低，导致其对¹³C分辨率下降，从而引起叶片δ¹³C值增加；同时，饱和水气压差升高，引起气孔导度下降，但是饱和水汽压差升高对光合速率限制作用可能高于气孔导度，从而导致油松和白榆叶片δ¹³C值降低。由于叶片δ¹³C是环境因子与生物因子综合作用的结果，有关树木叶片δ¹³C与水热因子关系有待于进一步研究。

本研究结果表明，樟子松叶δ¹³C值显著高于彰武小钻杨和油松，高于白榆，表明樟子松具有较高的水分利用效率，彰武小钻杨和油松次之，散生白榆最低。这与其他人的研究结果基本一致，例如，丁晓纲等^[38]比较了毛乌素沙地腹地樟子松和油松人工林生长季水分利用效率，发现生长季樟子松水分利用效率始终高于油松。高素华等^[39]研究表明，北方多粗沙区杨树水分利用效率显著高于白榆。本研究中，尽管樟子松土壤含水量较高，但是樟子松具有较高的水分利用效率，主要是因为樟子松叶片具有发达的角质层、气孔小而多、比叶面积小等特征^[40]，导致其扩散效应能力低，从而具有较高的水分利用效率。另外，与其他树种相比，樟子松根系分布较浅^[22,41]，容易受到干旱胁迫的威胁；因而采取更加保守的水分利用方式，从而具有较高的水分利用效率。尽管彰武小钻杨是落叶阔叶树种，但是其水分利用效率与针叶树种油松没有明显差异（图3），主要是因为彰武小钻杨的根系分布比油松深^[22]，可以利用深层土壤水和地下水；但是彰武小钻杨蒸腾耗水量高于油松，且油松和杨树0～60 cm土壤含水量并没有显著差异，因而彰武小钻杨和油松具有相似的水分利用效率。油松的气孔开度大于樟子松，而气孔密度小于樟子松^[40]，从而导致其扩散能力高于樟子松。因此，油松的水分利用效率低于樟子松。与其他树种相比，白榆具有最低的水分利用效率，可能是因为白榆叶片的气孔导度和比叶面积较大，其扩散能力高，导致白榆具有较低的水分利用效率。例如，陈仁升等^[42]研究表明，西北干旱区榆树的气孔导度显著高于樟子松和二白杨（Populus gansuensis C. Wang et H. L. Yang）。另外，白榆的土壤含水量高于杨树和油松，也可能是其具有较低水分利用效率原因之一。

樟子松的水分利用效率显著高于彰武小钻杨和油松，而彰武小钻杨和油松水分利用效率显著高于白榆。随着季节变化，樟子松、油松和白榆水分利用效率呈降低趋势，而彰武小钻杨水分利用效率呈先增加后降低趋势。叶片δ¹³C逐月变化量与水热因子的响应程度因树种而异，樟子松和白榆叶δ¹³C逐月变化量与水热因子相关关系不显著，而油松叶δ¹³C逐月变化量与降水量、气温和湿度呈显著负相关关系，而彰武小钻杨叶δ¹³C逐月变化量与气温呈显著正相关关系。不同树种水分利用效率差异不仅与树种生物学特性有关，而且也与受外界水热因子影响程度有关。因此，造林树种选择过程中，不仅要考虑树种生物学特性，还应该考虑树种对外界水热因子响应程度。