-
黄土高原是干旱半干旱地区,降水量少,降水分布不均,生态环境脆弱。自20世纪50年代以来,该地区采取了一系列的植被恢复措施,并取得了显著成效,但由于不合理的植被建设,导致林水矛盾激化,部分林地出现“小老头树”和土壤水分亏缺等现象[1-3]。基于此,研究黄土高原地区林地树种耗水规律尤为重要。99.8%以上的植物蒸腾耗水来自树干液流,液流速率可以较为准确地描述单木蒸腾耗水过程[4-5]。早期对植物蒸腾研究多集中于昼间,忽略夜间液流的存在。随着观测技术的发展与研究时间尺度的精细化,越来越多的研究发现整日树干液流可划分为昼间与夜间,昼夜液流变化特征具有差异性[6]。昼间液流量占整日液流量的80%~96%,是植物生长的主要环节,研究昼间液流的变化特征有助于理解植物的生长策略[7];夜间液流量占比为4%~20%[8],占比较少,但目前研究普遍认为[9]:夜间液流有助于植物适应干旱环境,在补充水分、传输营养、供应氧气等方面发挥着重要的作用,因此在研究植物蒸腾耗水时,夜间液流不可被忽视。
已经有大量的研究表明:植物液流速率受自身生理结构和气象因子的影响。任启文等[10]研究发现:冀北山地油松Pinus tabulaeformis液流速率与气温、太阳辐射、水汽压亏缺、风速呈正相关,与空气湿度呈负相关;贾国栋等[11]通过研究北方土石山区的油松、刺槐Robinia pseudoacacia发现:影响植物蒸腾速率的主要气象因子为太阳辐射、水汽压亏缺和风速;也有研究认为[12]:不同时间尺度下树干液流与气象因子的响应不同。上述研究关注于整日树干液流与气象因子关系研究,然而昼夜树干液流对环境因子响应机制不同,若不考虑昼夜液流变化特征及其对环境因子响应差异性,可能会降低植物耗水估算精度,影响林地水资源管理。目前,关于晋西黄土区生长季刺槐昼夜树干液流速率的研究较少,不同月份气象因子对昼夜树干液流速率的影响有待进一步研究。基于此,本研究以山西省吉县蔡家川流域的主要造林树种刺槐作为研究对象,定位观测其树干液流和气象因子。主要研究目的是①分析刺槐昼夜树干液流变化特征;②分析各月夜间树干液流对总液流的贡献率;③探究各月昼夜树干液流变化对气象因子的响应;④建立各月昼夜树干液流与气象因子关系模型,估算样地林分耗水量。
-
选取生长良好,树干通直,无病虫害的刺槐标准木8株。样木基本情况如表1。各样树树干液流采用Granier热扩散探针(TDP)观测,观测时间为2021年5月1日至9月30日。在大约1.3 m胸高处选择观测点,并刮去树干死皮,用配套工具沿垂直方向钻2个小孔,将热源探针插入上方孔中,感应探针插在下方孔中,使用固体胶将树体与探针的空隙涂抹密封。将电缆线固定在树干上并与电源连接,用铝箔纸将安装部位整个树干包裹起来,防止阳光直射和雨水造成数据误差。将数据采集器设置为30 s采集1次数据,隔15 min计算平均值并记录。
样树编号 胸径/cm 树高/m 边材面积/cm2 样树编号 胸径/cm 树高/m 边材面积/cm2 1 22.5 10.2 142.2 5 16.8 9.5 85.2 2 21.9 11.2 141.5 6 18.0 11.0 93.4 3 20.6 9.8 138.7 7 17.4 10.3 90.2 4 21.3 11.2 140.2 8 18.2 10.8 92.6 Table 1. Main parameters of sample trees
树干液流速率采用Granier经验公式获得。计算公式为:
式(1)中:ΔTm为无液流时加热探针与参考探针的最大温差(℃);ΔT为瞬时温差值(℃);Fd为树干液流速率(cm3·h−1·cm−2)
-
研究区刺槐样地里安装有全自动气象站,可全天观测林冠上方(距地面20 m)气象因子。气温(T,℃)和空气相对湿度(HR,%)采用HMP155A空气温湿度传感器观测,太阳辐射(Rs,W·m−2)采用CMP-3总辐射传感器观测,风速(SW,m·s−1)采用010C-1风速传感器观测,林内20 cm土壤温度(Ts,℃) 采用109土壤温度传感器观测。设置数据采集时间间隔为30 min,并采用CR1000数据采集与记录器收集数据。选取水汽压亏缺(DVP)指标综合反映空气相对湿度与空气温度的协同效应。
-
采用baseliner软件将温差数据转化为液流数据,使用Excel整理和计算树干液流、气象因子数据,使用SPSS软件与R语言进行相关分析,使用Origin绘图。图表中所用数据均为平均值±标准差。采用配对样本t检验分析各月整日树干液流量与昼间树干液流量的差异性;采用Pearson相关分析法分析气象因子与树干液流的相关关系;采用随机森林回归模型中相对重要性得分指标[IncMSE,IncMSE指去除自变量后随机森林回归模型预测误差的百分数,自变量关于IncMSE得分越高,说明自变量越重要]与逐步回归方程中因子进入方程时拟合度(R2)的变化值分析气象因子对树干液流的影响程度(随机森林是基于分类树的算法[13],该算法预测精度较高,被多个领域用于评估指标重要性);采用随机森林回归法与逐步回归法构建气象因子与树干液流速率相关模型。
-
根据观测数据可知:研究期间太阳辐射、气温、土壤温度、空气相对湿度整体呈下降趋势,日太阳辐射为14.36~451.83 W·m−2,平均辐射为218.42 W·m−2;日气温为7.5~28.6 ℃,平均气温为20.0 ℃;日地温为12.0~28.6 ℃,平均地温为22.0 ℃;日相对湿度为23.22%~98.83%,平均相对湿度为66.58%。降水主要在8—10月,最大降水发生在2021年9月19日,降水量为61.7 mm。
-
由图1可知:各月刺槐树干液流日变化规律相对一致,呈先增后减变化趋势。晴天树干液流启动时间在6:00左右,启动后树干液流速率迅速上升,到12:00左右达到峰值,随后液流开始下降;夜间(20:00—6:00)有明显的树干液流,变化不稳定,但是呈现一定规律性,因此有必要进一步研究夜间液流的变化特征。
-
由图2可知:各月夜间刺槐树干液流速率均表现出递减趋势,前半夜(20:00—24:00)树干液流速率大于后半夜(24:00—6:00)树干液流速率,整个夜晚树干液流速率处于持续下降状态,凌晨5:00左右液流速率降到最低。
-
计算刺槐生长季各月昼间、夜间树干液流平均速率与它们对整日液流的贡献率(表2)。结果表明:各月昼夜刺槐树干液流速率及其贡献率差异明显。在5—9月,昼夜平均液流速率呈先增后减变化,昼间平均树干液流速率在7月最大,9月最小;夜间平均液流速率在6月最大,9月最小;各月昼间液流速率对整日液流速率贡献率为88%~93%,贡献率在9月最小,8月最大;各月夜间树干液流速率对整日液流速率贡献率为7%~12%,9月贡献率最大。
月份 昼间平均液流速率/
(cm3·h−1·cm−2)贡献
率/%夜间平均液流速率/
(cm3·h−1·cm−2)贡献
率/%5 2.71±1.07 91.26 0.29±0.23 8.74 6 2.94±1.05 91.16 0.43±0.21 8.84 7 3.24±0.85 92.25 0.35±0.17 7.75 8 3.18±1.61 92.52 0.33±0.18 7.48 9 1.74±1.26 88.39 0.27±0.12 11.61 Table 2. Diurnal and nocturnal characteristics of sap flow rate in growing season
对整日刺槐树干液流量与昼间树干液流量进行相关性分析和配对样本T检验发现(表3):各月整日刺槐树干液流量与昼间树干液流量呈极显著相关(P<0.01), 整日液流量与昼间液流量呈极显著差异(P<0.01),说明各月夜间树干液流对整日树干液流的影响不可忽视。
月份 相关分析 配对样本t检验 R2 P t P 5 0.991 0.000 12.582 0.000 6 0.993 0.000 9.261 0.000 7 0.998 0.000 9.054 0.000 8 0.976 0.000 8.234 0.000 9 0.993 0.000 7.982 0.000 Table 3. Paired samples t test between daily and diurnal sap flow of R. pseudoacacia
-
Pearson相关分析(表4)结果表明:各月昼间刺槐树干液流速率与太阳辐射、土壤温度、气温、水汽压亏缺呈极显著正相关(P<0.01),与空气相对湿度呈极显著负相关(P<0.01),5、6月风速与树干液流速率不存在相关关系(P>0.05)。各月夜间刺槐树干液流速率与土壤温度、气温、水汽压亏缺度呈极显著正相关(P<0.01),与空气相对湿度呈极显著负相关(P<0.01)。
昼夜 月份 太阳辐射 气温 土壤温度 水汽压亏缺 相对湿度 风速 昼间 5 0.575** 0.483** 0.429** 0.290** −0.278** 0.067 6 0.670** 0.539** 0.331** 0.389** −0.305** −0.031 7 0.762** 0.767** 0.497** 0.691** −0.722** 0.192** 8 0.689** 0.747** 0.546** 0.641** −0.638** 0.295** 9 0.674** 0.737** 0.608** 0.696** −0.692** 0.331* 夜间 5 0.582** 0.585** 0.556** −0.266** 0.351** 6 0.213** 0.224** 0.043** −0.049** 0.123** 7 0.512** 0.475** 0.483** −0.394** 0.338** 8 0.747** 0.546** 0.641** −0.638** 0.295** 9 0.141** 0.160** 0.107** −0.058** 0.215** 说明:**表示极显著相关(P<0.01); * 表示显著相关(P<0.05) Table 4. Correlation analysis between diurnal and nocturnal sap flow rate of R. pseudoacacia and meteorological factors
基于随机森林模型中IncMSE指标(IncMSE越大自变量越重要)与逐步线性回归R2变化值来判断各月气象因子对刺槐树干液流速率的影响程度。昼间,刺槐在不同月份产生树干液流速率的驱动因子差异不大,随机森林IncMSE指标与逐步线性回归均显示各月对树干液流速率影响较大的气象因子为太阳辐射和气温。夜间,刺槐在不同月份产生树干液流速率的驱动因子存在较大差异,随机森林IncMSE(图3)显示:5月气象因子对树干液流影响程度从大到小排序为水汽压亏缺、土壤温度、风速、相对湿度、气温;6月为土壤温度、水汽压亏缺、相对湿度、气温、风速;7月为气温、相对湿度、风速、土壤温度、水汽压亏缺;8月为水汽压亏缺、相对湿度、气温、土壤温度、风速;9月为风速、水汽压亏缺、土壤温度、相对湿度、气温。夜间,逐步线性回归表明(表5):5、6月土壤温度对刺槐树干液流速率影响最大;8月水汽压亏缺对树干液流速率影响最大;7月气温对树干液流速率影响最大;9月风速对树干液流速率影响最大。气象因子对昼夜树干液流速率的影响程度在2种分析方法上存在不同,但2种方法都表明昼间气象因子对刺槐树干液流速率的影响在各月差异较小,太阳辐射与气温是影响各月昼间树干液流速率的主导因子;夜间气象因子对树干液流速率的影响在不同月份差异较大。除7月外,各月均表现出水汽压亏缺对夜间树干液流速率影响的重要性。
因子 逐步回归R2变化量 昼间 夜间 5月 6月 7月 8月 9月 5月 6月 7月 8月 9月 太阳辐射 − 0.449 0.095 0.015 0.040 − − − − − 气温 0.334 0.026 0.589 0.559 0.543 0.127 − 0.362 − 0.008 水汽压亏缺 0.152 0.150 0.016 0.029 − 0.057 0.148 0.005 0.360 0.232 相对湿度 0.087 − 0.032 0.053 − 0.014 0.084 0.050 0.117 0.008 土壤温度 − 0.007 0.007 0.074 0.014 0.342 0.350 − − − 风速 − − 0.002 0.029 0.018 0.042 − 0.158 0.030 0.246 说明:−表示气象因子未进入到逐步回归方程,以上数据均具有显著性 Table 5. Changes of stepwise regression R2 of R. pseudoacacia sap flow rate and meteorological factors
-
采用随机森林回归法与逐步回归法构建各月昼夜树干液流与气象因子的相关模型。随机森林回归模型无解析表达式,逐步回归模型表达式如表6。根据拟合度来评估2个模型准确度,在6、7、8月2个模型拟合较好,5、9月2个模型拟合较差。总体而言,基于随机森林回归构建的模型拟合度优于逐步回归模型。
昼夜 月份 逐步回归模型 随机森林回归法拟合度/% 模型 拟合度/% 昼间 5 y=−1.29+0.579x2−3.427x4−0.063x5 57.3 67.3 6 y=−0.767+0.003x1+0.301x2−0.149x3−0.869x4 63.2 73.2 7 y=1.509+0.002x1+0.444x2−0.194x3−1.862x4−0.053x5+0.071x6 74.1 86.3 8 y=1.864+0.001x1+0.506x2−0.208x3−2.468x4−0.065x5+0.144x6 75.9 90.6 9 y=1.104+0.001x1+0.300x2−0.204x3+0.297x6 61.5 81.1 夜间 5 y=2.852+0.002x1+0.182x2−0.279x3+0.154x6 48.2 59.4 6 y=4.943+0.003x1+0.346x2−0.363x3−1.252x4−0.025x5−0.110x6 61.6 67.7 7 y=3.966+0.003x1+0.415x2−0.209x3−2.094x4−0.069x5 57.5 61.6 8 y=4.869+0.002x1+0.486x2−0.311x3−2.296x4−0.074x5+0.091x6 50.7 63.9 9 y=−5.886+0.502x2−3.033x4 49.4 54.6 说明:y表示液流速率,x1表示太阳辐射,x2表示气温,x3表示土壤温度,x4表示水汽压亏缺,x5表示相对湿度,x6表示风速 Table 6. Models of day and night sap flow and meteorological factors based on stepwise regression method and random forest algorithm
-
根据随机森林回归模型推算刺槐林分各月蒸腾量可知:7月林分蒸腾量最大,为37.5 mm,其他月份依次为8月33.5 mm,6月30.8 mm,5月30.6 mm,9月林分蒸腾量最小,为24.3 mm。推算结果符合刺槐生长规律。生长旺期的蒸腾最大,生长末期的蒸腾最小。
Diurnal and nocturnal changes in stem sap flow of Robinia pseudoacacia during growing season and its response to meteorological factors
doi: 10.11833/j.issn.2095-0756.20220159
- Received Date: 2022-02-13
- Accepted Date: 2022-06-19
- Rev Recd Date: 2022-06-12
- Available Online: 2022-11-21
- Publish Date: 2022-12-20
-
Key words:
- Robinia pseudoacacia /
- sap flow rate /
- growing season /
- diurnal variation /
- meteorological factors /
- loess region of western Shanxi Province
Abstract:
Citation: | ZHANG Rong, BI Huaxing, JIAO Zhenhuan, WANG Ning, ZHAO Danyang, YUN Huiya, HUANG Jinghan. Diurnal and nocturnal changes in stem sap flow of Robinia pseudoacacia during growing season and its response to meteorological factors[J]. Journal of Zhejiang A&F University, 2022, 39(6): 1238-1246. doi: 10.11833/j.issn.2095-0756.20220159 |