宁夏河东沙地不同坡度柠条锦鸡儿（Caragana korshinskii）水分利用策略差异

高 阳， 韩 磊,2,3， 柳利利， 王娜娜， 彭 苓， 周 鹏， 展秀丽

（1.宁夏大学地理科学与规划学院，宁夏 银川 750021；2.中阿旱区特色资源与环境治理国际合作联合实验室，宁夏 银川 750021；3.宁夏旱区资源评价与环境调控重点实验室，宁夏 银川 750021；4.宁夏大学农学院，宁夏 银川 750021）

全球变暖加速了地球系统的水循环，使全球和区域降水格局发生变化［1］。宁夏河东沙地地处毛乌素沙地西南缘，在空间上是我国北方农牧交错区最具过渡性地域特征的沙地，也是黄土高原向鄂尔多斯台地过渡的地带，属于气候变化的敏感区［2］。研究干旱区典型固沙植物的水分利用策略，对评估降水格局变化背景下该地区的生态可持续发展具有重要意义［3］。稳定同位素技术在定量分析植物水分利用策略的相关研究中得到广泛应用［4-6］。研究表明，除极少数盐生、旱生和湿地植物外，陆地植物根系吸收的水分在到达叶片或未栓化的幼嫩枝条之前一般不发生同位素分馏［7］，且各潜在水源之间的氢氧稳定同位素组成存在显著差异。因此，通过比对植物木质部水与潜在水源的稳定同位素组成，可有效辨别和区分植物水分来源［8］。

多元线性混合模型（IsoSource）和贝叶斯混合模型（MixSIAR）被广泛应用于阐明干旱区植物的水分利用策略，但两种方法的量化结果可能存在一定差异［9-10］。曾祥明等［11］在西南喀斯特地区利用Mix⁃SIAR模型和IsoSource模型量化刺楸（Kalopanax sep⁃temlobus）、香椿（Toona sinensis）和化香（Platycarya strobilacea）的水分来源时，发现MixSIAR 模型计算植物对各水源的利用比例适用性高于IsoSource 模型；张宇等［12］利用IsoSource模型和MixSIR、MixSIAR模型研究了红瑞木（Swida alba）的水分来源，表明MixSIR 模型能够更精确地给出红瑞木主要的吸水层位以及对不同深度土层土壤水的吸收利用比例；Wang 等［13］利用IsoSource 模型和贝叶斯混合模型（SIAR、MixSIR、MixSIAR）研究了黄土高原半干旱地区黄荆（Vitex negundo）、细裂叶莲蒿（Artemisia gmeli⁃nii）、白刺花（Sophora davidii）的水分来源，发现SIAR、MixSIAR模型比MixSIR模型具有更好的植物水分溯源效果。部分学者利用IsoSource 模型对典型固沙植物的水分来源和水分利用特征进行了研究［14-15］，但该模型只是基于质量平衡进行迭代算法，存在未考虑随机测量误差、纳入潜在水源中稳定同位素组成的不确定性等缺点，使得IsoSource 模型的量化结果可靠性降低，因此有必要对不同模型的水源区分性能进行对比研究，分析其对解析植物水分利用策略的适用性。

本研究以宁夏河东沙地典型人工固沙植物柠条锦鸡儿为研究对象，利用氢氧稳定同位素技术、直接对比法、IsoSource 模型和MixSIAR 模型分析生长季不同时期、不同坡度柠条锦鸡儿的水分利用策略差异，评估两种模型计算植物水分来源的适用性，为阐释全球气候变化背景下干旱沙地人工固沙植物水分利用策略提供依据。

1 研究区概况

试验地位于宁夏回族自治区灵武市白芨滩国家级自然保护区（106°22′03″E，38°03′34″N，海拔1150 m）（图1），该地区地处毛乌素沙地西南缘，属温带大陆性气候。多年平均气温6.2~8.8 ℃，多年平均降水量198.9 mm，多年平均潜在蒸发量1928.4 mm，降水主要集中在7—9 月［16］。人工林树种主要有柠条锦鸡儿、丝棉木（Euonymus maackii）、沙枣（Elaeagnus angustifolia）、刺槐（Robinia pseudoaca⁃cia）、新疆杨（Populus alba）、沙拐枣（Calligonum mongolicum）和细枝山竹子（Corethrodendron scopari⁃um）等。地带性土壤以淡栗钙土、棕钙土和风沙土为主。

图1 研究区位置示意图Fig.1 Schematic diagram of the location of the study area

2 数据与方法

2.1 样品采集与测定

本研究在柠条锦鸡儿人工林内自然分布的4个坡度（6°、10°、16°和24°）进行研究，每个坡度分别布设3 块20 m×20 m 的样地，柠条锦鸡儿样地基本情况和土壤基本理化性质分别见表1和表2。2018年5—10 月每月分别采集一次植物木质部和土壤样品。在4个坡度样地内分别选取树龄、冠幅相似、地径相近、生长状况良好的3棵植株，用枝剪分别从植株4 个方向剪取中部直径约0.5 cm，长约3~5 cm 的非绿色栓化小枝，迅速剥除外皮与韧皮部，装入10 mL 的玻璃样品瓶中，并用Parafilm 封口膜密封后置于随身携带的保温箱（0~5 ℃）中冷藏，每个坡度样地均采集3 瓶样品作为重复。同时，在所选柠条锦鸡儿灌丛下方地势相对均一的地点，用土钻按照0~10 cm、10~20 cm、20~40 cm、40~60 cm、60~80 cm、80~120 cm、120~160 cm 和160~200 cm 分层采集土壤样品，每层采集2瓶土样，立即密封在已标记好的样品袋中，放入保温箱中。所有植物木质部和土壤样品带回实验室后放入冰柜（-20 ℃）中冷冻保存直到水分抽提。在无林空旷地内放置3个标准雨量筒（直径20 cm 的圆桶上放置漏斗，下端与雨量筒相连，漏斗口放置一个乒乓球，防止降水样品蒸发）收集降水量大于等于0.1 mm的降水样品，每次降水结束，立即将3 个标准雨量筒中的降水样品混合后转入10 mL玻璃样品瓶中，并用Parafilm封口膜密封后放入保温箱中带回实验室，置于冰箱（4 ℃）中冷藏保存待测。2018 年5—10 月，共收集到降水样品34个。利用Vantage Pro 2 小型气象站（Davis Inc，USA），每隔15 min 对大气温度、降水量等气象数据进行记录。

表1 不同坡度柠条锦鸡儿人工林基本情况Tab.1 Basic conditions of Caragana korshinskii plantations on different slopes

表2 土壤基本理化性质Tab.2 Basic physical and chemical properties of soil

利用全自动真空冷凝抽提系统（LI-2100，北京理加联合科技有限公司）抽提土壤水和植物木质部水。所有水样经0.45 μm 醋酸纤维素滤膜过滤后，利用基于稳定同位素红外光谱技术的液态水同位素分析仪（T-LWIA-45-EP，ABB-Los Gatos Research，CA）测定降水、土壤水、植物木质部水样品的δD 和δ18O 组成，δD 和δ18O 的测量误差分别低于±1‰和±0.3‰。测定结果以相对于维也纳标准平均海洋水（Vienna standard mean ocean water，VSMOW）的千分差值（‰）表示：

式中：Rsample和RV-SMOW分别为样品和维也纳标准平均海洋水中的稳定同位素比率（D/H和18O/16O）（‰）。

利用低温真空蒸馏法提取植物木质部水时，可能同时蒸馏出与水分子具有相似光谱吸收峰的甲醇和乙醇等有机物质，导致稳定同位素的测量误差超过仪器精度［17］。为减小测量误差，本研究利用Los Gatos 的光谱分析软件对植物木质部水的稳定同位素组成进行校正［18-19］。由于柠条锦鸡儿是干旱区常见的旱生植物，而旱生植物的根系在吸收水分的过程中会发生氢稳定同位素分馏［6］，因此，本研究采用氧稳定同位素计算植物的水分来源。

2.2 研究方法

试验地地下水埋深深度为18.0 m，人工固沙植物一般难以吸水和利用地下水，因此，本研究只将不同土层深度的土壤水作为潜在水源。

（1）多元线性混合模型（IsoSource）

将柠条锦鸡儿木质部水δ18O 值作为混合物，不同土层深度的土壤水及其δ18O 值作为来源代入模型，将来源增量设为1%，将质量平衡公差设为0.01‰，得到各潜在水源的平均水分贡献率。

（2）贝叶斯混合模型（MixSIAR）

将柠条锦鸡儿木质部水δ18O 值作为混合物数据，各层土壤水δ18O值作为源数据，并将分馏数据设为0，默认不发生氧稳定同位素分馏，Error structure选择Resid Only，马尔卡夫链蒙特卡罗（MCMC）的运行长度设为normal，模型结果使用Gelman-Rubin 和Geweke 诊断确定MixSIAR 模型是否收敛，如果Mix⁃SIAR 模型不收敛则增加MCMC 的运行时长，选择long 或者very long，模型输出结果为各潜在水源对应的平均贡献率［10］。

由于目前植物对不同水源的实际利用值无法直接观测［20］，因此，本研究通过评估植物木质部水氧稳定同位素的观测值和预测值之间的匹配程度来间接评估IsoSource 模型和MixSIAR 模型的性能［13］。将测定的柠条锦鸡儿木质部水δ18O 值作为观测值（Oi），预测值（Pi）的计算基于柠条锦鸡儿木质部水δ18O值是各潜在水源混合物的假设［21］。预测值Pi计算如下：

式中：n为植物木质部水的潜在水源个数（n=3）；fi为IsoSource 模型和MixSIAR 模型计算植物对第i个潜在水源的利用率（%）；δA为各潜在水源的氧稳定同位素值（‰）。

本研究根据不同土层深度土壤水δ18O值的相似性及变化趋势合并相邻土层［22］，将0~200 cm的土壤剖面自上而下依次划分为浅层（0~40 cm）、中层（40~120 cm）和深层（120~200 cm）土壤，各土层土壤水δ18O 值取对应深度土壤水δ18O 值的几何平均值。本研究采用Microsoft Excel 2019 软件进行数据整理，在SPSS 26.0中对不同坡度样地土壤水和柠条锦鸡儿木质部水δ18O值的差异显著性进行单因素方差分析。利用Origin 2018软件进行图像绘制，图表中数据为平均值±标准差。

2.3 模型评价

（1）均方根误差（RMSE）：

式中：n为验证样本的数量；Oi和Pi分别为植物木质部水δ18O 的观测值和预测值。RMSE 越低，表明该模型误差越小，预测结果越准确。

（2）模拟优度指数（G）：

式中：为观测值的平均值。G值越接近1 表示模型预测的可靠性越高；G为1 表示模型预测非常准确；G值为负值表示模型预测可靠性不高。

（3）平均绝对百分比误差（MAPE）：

MAPE<10%表示高度拟合；10%50%表示不可行拟合。

（4）最大预测误差（MaxE）：

（5）最小预测误差（MinE）：

3 结果与分析

3.1 降水分布及降水中δD和δ18O变化特征

试验地2018 年5—10 月植物生长季总降水量为172.1 mm，主要集中在7—9月（135.1 mm）（图2），占生长季总降水量的78.5%。降水中δD、δ18O 的变化范围分别介于-74.1‰~7.4‰、-10.5‰~0.7‰之间（图2）。对宁夏河东沙地事件尺度降水中δD、δ18Ο进行线性回归分析，拟合得到局地大气降水线（LM⁃WL：δD=5.12δ18O-12.53，R2=0.86，n=34）（图3）。LM⁃WL 的斜率和截距均小于全球大气降水线（GMWL：δD=8δ18O+10）［23-24］，说明雨滴在非饱和大气降落过程中受到了较强烈云下二次蒸发效应的影响，降水中稳定同位素发生不平衡分馏和重同位素富集，使得LMWL的斜率偏小［25］。

图2 不同时期气温、降水量及降水中δD、δ18O变化Fig.2 Variation of temperature,precipitation and δD and δ18O in precipitation during different periods

图3 局地大气降水线及土壤水、柠条锦鸡儿木质部水δD和δ18O变化Fig.3 Variation of local meteoric water line and δD and δ18O in soil water and xylem water of Caragana korshinskii

3.2 不同坡度柠条锦鸡儿木质部水δ18O变化特征

不同坡度样地柠条锦鸡儿木质部水δ18O差异显著（P<0.05），其中，坡度10°、16°和坡度24°样地柠条锦鸡儿木质部水δ18O随坡度的增大逐渐减小（图4），变化范围分别介于-6.9‰~-2.1‰、-7.7‰~-5.3‰和-5.9‰~-0.4‰，随时间的变化趋势基本一致，生长季初期（5月）值最大，生长季后期（9—10月）值最小；而坡度6°样地柠条锦鸡儿木质部水δ18O 的变化范围介于-7.8‰~-5.2‰，生长季后期值最大，生长季中期（6—8月）值最小。

图4 生长季不同坡度柠条锦鸡儿木质部水δ18O变化Fig.4 Variation of δ18O in xylem water of Caragana korshinskii at different slopes during the growing season

3.3 土壤含水量与土壤水δ18O变化特征

随着坡度的增大，柠条锦鸡儿样地土壤含水量逐渐减小，由坡度6°的（6.30±1.01）%减小到坡度24°的（1.33±0.58）%（图5）。坡度24°样地土壤水δ18O在0~40 cm和40~120 cm土层最高，120~200 cm土层则表现出坡度10°最高、坡度16°最低。不同坡度柠条锦鸡儿样地土壤水δ18O随土层深度的增加逐渐减小（图6）。

图5 生长季不同坡度土壤含水量变化Fig.5 Variation of soil water content at different slopes during the growing season

图6 生长季不同坡度土壤水δ18O变化Fig.6 Variation of soil water δ18O at different slopes during the growing season

3.4 柠条锦鸡儿的主要吸水层位

如图7所示，生长季初期，坡度6°样地柠条锦鸡儿主要利用中层（60~80 cm）、深层（160~180 cm）土壤水，坡度10°样地柠条锦鸡儿主要利用中层（40~60 cm、100~120 cm）、深层（140~160 cm）土壤水，而坡度16°样地柠条锦鸡儿主要利用中层（80~100 cm、120 cm）土壤水，坡度24°样地柠条锦鸡儿的主要吸水层位在60~80 cm土层。生长季中期，随降水量增加，坡度6°样地柠条锦鸡儿与土壤水存在多个交点，利用直接对比法不能判断出水分来源，坡度10°样地柠条锦鸡儿转为以利用浅层和中层土壤水为主，而坡度16°样地柠条锦鸡儿仍主要利用中层和深层土壤水，坡度24°样地柠条锦鸡儿增加了对浅层土壤水的利用，同时坡度24°样地柠条锦鸡儿木质部水δ18O与中层土壤水存在交点，表明坡度24°坡样地柠条锦鸡儿主要利用浅层和中层土壤水。生长季后期，坡度6°、10°样地柠条锦鸡儿转为以利用浅层（20~40 cm）土壤水为主，而坡度16°、24°样地柠条锦鸡儿转为利用更深层次的土壤水。

3.5 IsoSource模型和MixSIAR模型对植物水分来源定量分析

如图8所示，生长季初期，IsoSource模型和Mix⁃SIAR模型计算坡度6°、10°和16°样地柠条锦鸡儿对土壤水利用率的适用性均较高，随坡度增大，Iso⁃Source模型计算的柠条锦鸡儿对中层和深层土壤水的利用率分别高于MixSIAR 模型的计算结果23.5%、34.7%和24.4%，而MixSIAR 模型计算坡度24°样地柠条锦鸡儿土壤水利用率的适用性更高。生长季中期，IsoSource 模型和MixSIAR 模型计算坡度10°样地柠条锦鸡儿对浅层和中层土壤水的利用率较为一致（为91%），但IsoSource 模型计算坡度24°样地柠条锦鸡儿对浅层和中层土壤水的总利用率高于MixSIAR 模型的计算结果7.6%。生长季后期，坡度6°样地柠条锦鸡儿的主要水源由深层转移至浅层土壤，IsoSource 模型计算坡度6°、16°、24°样地柠条锦鸡儿对土壤水利用率的适用性高于Mix⁃SIAR 模型，坡度16°样地柠条锦鸡儿对深层土壤水的利用率增加至最大（52.4%），而坡度24°样地柠条锦鸡儿的主要水源由中层转移至深层土壤（73.9%）；坡度10°样地柠条锦鸡儿的主要水源由中层转移至浅层土壤，MixSIAR模型适用性更高，柠条锦鸡儿对浅层土壤水的利用率为38.3%。

3.6 IsoSource模型和MixSIAR模型的植物水分溯源效果

如表3 所示，MixSIAR 模型在计算坡度10°、16°、24°样地柠条锦鸡儿水分来源时RMSE 更小；IsoSource 模型计算坡度6°样地柠条锦鸡儿水分来源时RMSE 更小。MixSIAR 模型计算的坡度6°、16°、24°样地柠条锦鸡儿的G为负值，区分水源性能可靠性降低；当坡度为24°时，IsoSource 模型的MAPE>50%。IsoSource 模型计算坡度6°样地柠条锦鸡儿水分来源时具有更好的植物水分溯源效果；MixSIAR 模型计算坡度10°样地柠条锦鸡儿水分来源的溯源效果更好；计算较大坡度柠条锦鸡儿水分来源时，IsoSource模型的准确性低于MixSIAR模型。

表3 IsoSource模型和MixSIAR模型的预测性能Tab.3 Predictive performance of IsoSource model and MixSIAR model

4 讨论

坡度、坡向等地形条件对地表光、热、水资源具有再分配作用，进而影响林地根系的生长和分布［26-27］，而植物的根系分布特征决定植物的水分利用策略［28］。已有研究［14］表明，随着坡度的增大，细根根系参数的最大值由中深层土壤向浅层土壤转移，呈现出表聚现象，并且随土层深度的增加，细根的根长密度、根表面积密度等根系参数与土壤水分和养分的变化趋势基本一致，呈递减趋势。在整个生长季，柠条锦鸡儿根据土壤含水量的相对高低，在不同土层间灵活转换主要水分来源。生长季初期，柠条锦鸡儿对中层土壤水利用率随坡度增大呈先增后减趋势；生长季中期，柠条锦鸡儿的主要水分来源随坡度的增大向浅层土壤转移，主要利用浅层和中层土壤水；生长季后期，柠条锦鸡儿主要水分来源随坡度的增大向深层土壤转移，吸收利用深层土壤水。这可能与柠条锦鸡儿具有二态性根系结构有关［29］。朱雅娟等［30］对乌兰布和沙漠柠条锦鸡儿的研究也表明，柠条锦鸡儿在不同季节分别利用浅层和较深层的土壤水分，推测其可能具有双型根系。生长季初期，柠条锦鸡儿需水量少，水分来源单一；而生长季中期，一方面，降水补给土壤水作用明显，土壤含水量显著增加，植物利用浅层和中层土壤水消耗的能量相对较小［31］，另一方面，土壤水分和养分含量供给充足，细根数量显著增加，柠条锦鸡儿优先利用浅层土壤水［27］；生长季后期，柠条锦鸡儿吸收利用深层土壤水这一稳定性水源，说明柠条锦鸡儿采取了资源依赖型的水分利用策略［32］。

IsoSource模型和MixSIAR模型在量化植物水分来源上存在一定差异［11］。IsoSource模型和MixSIAR模型解析坡度6°、10°和16°样地柠条锦鸡儿主要水分来源利用比例的适用性较高。然而这两种模型在计算坡度24°样地柠条锦鸡儿主要水分来源利用比例时存在一定差异，因此借助直接对比法的定性判断结果发现，MixSIAR模型解析坡度24°样地柠条锦鸡儿主要水分来源利用比例的适用性高于Iso⁃Source模型。曾祥明等［11］通过直接相关法判断植物主要吸水源来衡量IsoSource 模型和MixSIAR 模型的适用性，结果也表明利用MixSIAR 模型计算植物对各水源的利用比例适用性高于IsoSource 模型。同时根据模型性能评价结果，IsoSource 模型和Mix⁃SIAR 模型分别计算坡度6°、10°样地柠条锦鸡儿水分来源时具有更好的植物水分溯源效果；而在计算坡度16°、24°样地柠条锦鸡儿水分来源时，IsoSource模型的准确性低于MixSIAR模型。因此本研究认为IsoSource 模型更适合解析较小坡度（6°和16°）样地柠条锦鸡儿的水分利用策略；而MixSIAR 模型解析较大坡度（10°和24°）样地柠条锦鸡儿水分利用策略的适用性更好。Wang 等［13］利用RMSE、G、R、MaxE和MinE 模型性能评价指标对IsoSource 模型和3 种贝叶斯混合模型（SIAR、MixSIR、MixSIAR）的水源分配性能进行了对比评价，同样发现IsoSource 模型的预测虽然准确，但该模型的性能不稳定，而SIAR 模型和MixSIAR 模型比MixSIR 模型具有更好的植物水分溯源效果，可能是由于后验分布仿真算法的不同造成了这种差异。

IsoSource 模型被广泛运用于解析植物水分来源，但该模型只是基于稳定同位素质量守恒定律，未明确纳入水源中同位素组成的不确定性［9］。而MixSIAR 模型是基于MixSIR 模型和SIAR 模型的优点所做的改进，通过考虑源值、分类和连续协变量和先验信息的不确定性来改进更简单的线性混合模型，以提高模型结果的准确性［11］。同时MixSIAR模型提供了一个最优解，而不是可行解的范围［10］。本研究中坡度6°、10°、16°样地柠条锦鸡儿木质部水及不同土层深度土壤水同位素组成的空间差异性较小，即使IsoSource模型计算时未考虑植物木质部水及不同土层深度土壤水同位素组成的空间差异性，该模型计算结果的可信度也较高，因此两种模型均适用且计算结果的不确定性较小；而坡度24°样地柠条锦鸡儿木质部水及不同土层深度土壤水同位素组成的空间差异性较大，使得IsoSource 模型在计算坡度24°样地柠条锦鸡儿主要水分来源利用比例时标准差较大，增加了该模型计算坡度24°样地柠条锦鸡儿主要水分来源利用比例的不确定性。MixSIAR模型解决了IsoSource模型计算时未考虑植物木质部水和各潜在水源稳定同位素组成空间差异性的问题［8］，其计算结果的标准差显著低于IsoSource 模型计算结果的标准差，有效降低了该模型计算结果的不确定性，即MixSIAR 模型的计算结果能更精确的给出坡度24°样地柠条锦鸡儿主要的水分来源以及对不同土层深度土壤水的利用比例，这与张宇等［12］的研究结果一致。李楠［33］研究也发现，SIAR模型计算的各水源贡献率的标准差明显低于IsoSource 模型计算结果的标准差。本研究进行了IsoSource 模型和MixSIAR 模型植物水分溯源效果的对比研究，并选择出了更适合解析干旱区不同坡度样地植物水分来源的溯源模型，解答了研究者在研究干旱区植物水分来源时，应该选择哪种溯源模型的困惑。

5 结论

（1）柠条锦鸡儿对不同土层深度土壤水的利用存在明显的季节性差异，生长季初期，随着坡度的增大，柠条锦鸡儿对主要水分来源的利用率呈现出先增后减趋势；生长季中期，柠条锦鸡儿主要水分来源随坡度的增大由深变浅；而生长季后期，柠条锦鸡儿主要水分来源随坡度增大由浅变深，呈相反的变化趋势。

（2）IsoSource 模型更适合解析较小坡度（6°和16°）样地柠条锦鸡儿的水分利用策略；而MixSIAR模型解析较大坡度（10°和24°）样地柠条锦鸡儿水分利用策略的适用性更好。