采煤沉陷区生态修复树种蒸腾特征及能量收支

徐立杰, 党晓宏,2, 高 永,3, 冯亚亚, 李传福

(1.内蒙古农业大学 沙漠治理学院, 内蒙古 呼和浩特 010018; 2.内蒙古杭锦荒漠生态系统国家定位观测研究站, 内蒙古 鄂尔多斯 017400; 3.中央与地方共建风沙物理重点实验室, 内蒙古 呼和浩特 010018)

柠条(Caraganakorshinskii)、沙柳(Salixpsammophila)、沙棘(Hippophaerhamnoides)是半干旱地区主要的乡土树种，因其具备良好的防风固沙能力，同时也具有较高的的经济价值，因此成为该地区主要的生态修复树种，研究其蒸腾耗水量及树种对环境的改善能力对于半干旱区水分调控及生态修复具有重要意义。目前对于3种植物的蒸腾特征的研究已有较多研究[1-2]。植物蒸腾是干旱、半干旱地区水分散失的主要途径，但常见的测量单株植物蒸腾的方法在测定沉陷区植物蒸腾速率时具有一定的局限性，称重法和光合仪容易对植物叶片造成损伤，影响植物正常生长；茎流计测得的茎流量不完全等于蒸腾耗水量；蒸渗仪测定结果最准确但昂贵的价格和笨重的体型导致其难以在野外中应用[3]。探索一种能够在野外条件下测定植株实际蒸腾耗水的方法对研究荒漠地区植物蒸腾具有重要意义。“三温模型”是邱国玉等人提出的测算蒸散发的方法，由于该模型的核心参数是植被表面温度、参考表面温度和气温，故被称为“三温模型”，具有参数少、计算简单、应用范围广、结果精准等特点[4-6]，红外热成像技术能观测人体肉眼不可见的红外波段的光谱，将其转变为可见的热图像[7]，具有高通量性、非接触性、高分辨等特征[8]。本试验以毛乌素沙地采煤沉陷区常见生态修复树种为研究对象，探究热红外遥感+“三温模型”的方法在半干旱地区应用过程中，适合采煤沉陷地区植物蒸腾测量的参考叶片，并确定3种主要树种的蒸腾特征以及在蒸腾过程中能量收支情况，为沙区植被恢复提供理论依据和技术支持。

1 材料和方法

1.1 研究区概况

研究区位于内蒙古自治区鄂尔多斯市伊金霍洛旗(38°50′—39°50′N，109°30′—110°30′E)。该区域为典型的大陆性季风气候，多年平均气温7.3 ℃，全年温差巨大，风沙天气频发，年均降水量362 mm，年均潜在蒸发量2 297.4～2 838.7 mm。土壤主要以风沙土为主，有机质含量低，氮、磷、钾元素十分缺乏。该区域由于采矿活动导致野生植物退化严重，植被以柠条、沙柳、沙棘等人工植被为主，零星生长有油蒿(Artemisiadesertorum)和沙米(Agriophyllumpungens)等灌草。

1.2 试验设计

1.2.1 研究对象选择 本试验于2019年9月25日(晴朗无云的天气)进行，选择毛乌素沙地东北缘采煤沉陷区生态恢复示范区作为试验区，面积约为1 762.85 m2，海拔1 228 m。依据树冠不受遮蔽原则，在试验区内选择7株不受其他植株遮蔽，生长健康的柠条作为试验对象，用于参考叶片的筛选。同时在该区域选择株高冠幅相似的柠条(高120 cm左右)、沙棘(高110 cm左右)、沙柳(高250 cm左右)各三株作为试验对象，用于植物蒸腾特征的研究(参考叶片为绿色卡纸)，沙柳为人工种植，沙棘和柠条为飞播种植。

1.2.2 参考叶片布设 卡纸和便签纸是研究区容易获得的绿色纸张，因此使用上述两种纸张和烘干植物叶片(校本部烘干带至研究区)作为参考叶片，两种纸张均购买于当地文具品店。将两种纸张裁剪为长约5 cm，宽约2 cm的矩形参考叶片。将参考叶片布设于无蒸腾冠层上方，参考叶片布设的倾斜角度、方位参数尽量和同时观测的植物叶片保持一致，为保证参考叶片和植物叶片受到相同的温度和光照，拍摄前半小时将参考叶片布设于试验样地。

1.2.3 表面温度拍摄及提取 使用Fluke TiR4 FT热成像仪获取植物叶片及各参考叶片的温度，该仪器配备320×240焦平面阵列(FPA)探测器，选用标准红外镜头，视场为23°×17°，热敏度为≤0.05 ℃，发射率设定为0.95。测量时段为9：00—17：00，步长2 h，相邻的两次拍摄顺序相反[9]，镜头垂直于冠层且冠层顶端距镜头1 m以上[10]。将拍摄的图像导入电脑端，在FLUKE中国官网下载Fluke Smart View 4.3软件对图像中冠层范围进行识别，并对其冠层温度进行提取(温度单位修改为摄氏度)，重复提取15次，以15次的均值作为该时刻柠条灌丛的叶温。

1.2.4 气象因子测定 使用Vaisala WXT520移动气象站同步测定研究区气温、太阳辐射、空气湿度、风速等气象因素，气象数据 10 min的平均值由 CR200X数据采集器计算和记录。

1.2.5 基于“三温模型”的植物蒸腾计算 当地表被植被覆盖时，可以忽略土壤热通量，植被表面的能量平衡方程可以表示为：

LEc=Rn,c-Hc

(1)

式中：LEc为植被的潜热通量(W/m2),即植物蒸腾；Rn,c是植被吸收的太阳净辐射(W/m2)；Hc是显热通量(W/m2)。Hc可表示为

HC=ρCp(Tc-Ta)/ra

(2)

式中：ρ是空气密度,为1.39 kg/m3；CP是空气定压比热(在平均空气状况下,其值约为1.013×10-3MJ/(kg·K)；TC植被冠层温度;Ta是空气温度,单位均为K；ra是空气动力学阻抗(s/m)。其值可通过参考叶片计算,参考叶片无蒸腾作用,且不会导致周围环境发生显著改变,ra计算公式为:

ra=ρCp(Tc,p-Ta)/Rn,cp

(3)

式中：Tc,p是参考植被的冠层温度(K)；Rn,cp是参考植被吸收的净辐射(W/m2)。

结合上述公式即可得到植被蒸腾的计算公式：

LEc=Rn,c-Rn,cp〔(Tc-Ta)/(Tc,p-Ta)〕

(4)

式中：(Tc-Ta)/(Tc,p-Ta)即为植物蒸腾扩散系数(hat),可用于反映植物水分状况。

净辐射可用以下公式推算：

Rn=(1-α)Rs+ΔR1

(5)

Rn=(1-α)Rs+5.675×10-8

(6)

ΔR1=(0.4+0.6Rs/Rs0)

(7)

(8)

式中：Rn太阳辐射(W/m2)；Rs是太阳辐射(W/m2)；α是地表反照率,(本试验中,α取值为0.22)； ΔR1是净长波辐射(W/m2)；Rs0是晴天太阳辐射量(W/m2),本研究在晴朗无云的典型晴天进行,因此认为Rs0与Rs相等；εa是大气发射率；εs是地表放射率(植被地表放射率取0.98)；σ是斯蒂芬—玻尔兹曼常数〔σ=5.675×10-8W/(m2·K4)〕。

2 结果与分析

2.1 研究区气象因子日变化

在“三温模型”中，气温和太阳辐射是测量植物蒸腾速率的关键因子。由图1可知，气温和太阳辐射均呈现先增加后减小的趋势；气温日变化类似“钟罩”形，波动范围较小，早上9：00气温最低，13：00时气温有最大值23.63 ℃；太阳辐射在13：00时达到峰值741.07 W/m2，13：00后太阳辐射变化幅度加快近于线性变化，17：00达到最小值123.16 W/m2。风速日变化同气温及太阳辐射日变化规律一致，13：00时风速最大为2.61 m/s，17：00时风速最小为0.81 m/s；空气湿度日变化趋势同气温及太阳辐射相反，9：00时太阳辐射和气温较低，空气中水分较多，湿度最高，为39.21%。随着气温和太阳辐射的上升，空气中的水分蒸发，空气湿度开始降低，在13：00处达到最低值22.94%，13：00—15：00间，相对湿度和气温变化幅度一致，15：00后湿度上升幅度加快，在17：00时湿度达到测定时间内第二峰值30.49%。

图1 研究区气象因子日变化

2.2 参考叶片的筛选

由图2可知，各参考叶片温度日变化为“单峰型”，且均在13：00时达到峰值，但便签纸温度在太阳辐射较高时远低于其他两种参考叶片，与植物叶片温度的差值较小，卡纸与干叶片温度变化在15：00之前有较高的拟合度，但是15：00以后随着太阳辐射的降低速度加快，卡纸与干叶片间的温差开始变大，17：00时各叶片温度大小关系为：干叶片>植物叶片>便签纸≈卡纸。经计算，当便签纸作为参考叶片时，柠条在9：00时的瞬时蒸腾速率为4.95 mm/h，远大于其余两种参考叶片同时刻的植物瞬时蒸腾速率，由图3可知，当以卡纸为参考叶片时，其计算结果与干叶片作为参考叶片的计算结果在9：00—15：00时有较高的吻合度，在13：00时达到最大值0.71 mm/h，在17：00时二者吻合程度较低。综上所述，以表面温度和计算结果作为选择参考叶片的指标发现绿色便签纸不适合作为采煤沉陷区柠条的参考叶片，绿色卡纸可以在9：00—15：00间替代干叶片作为植物蒸腾测量的参考叶片。

图2 各参考叶片温度日变化

图3 柠条蒸腾速率日变化

2.3 三种植物蒸腾特征

2.3.1 三种植物蒸腾速日变化规律及日蒸腾量 图4是3种植物11：00—15：00间瞬时蒸腾速率日内变化规律，3种植物蒸腾速率日内变化趋势均为“单峰型”，沙棘蒸腾速率在11：00时达到峰值0.61 mm/h，日蒸腾量为3.27 mm，柠条和沙柳瞬时蒸腾速率则在15：00时达到最大值，最大值分别为0.71和0.57 mm/h。日蒸腾量则为柠条3.80 mm，沙柳3.24 mm。和柠条相比，沙柳和沙棘蒸腾速率变化较为平缓，柠条在太阳辐射和温度高时蒸腾水平高于其他两种植物。

图4 3种植物瞬时蒸腾速率变化规律及日蒸腾量

2.3.2 三种植物蒸腾扩散系数变化规律 植被蒸腾扩散系数(hat)是衡量植被蒸散量并评价其水分利用状况的一种模型指标，可反映植物根系土壤水分状况，进而指示作物水分亏缺。hat取值范围为hat≤1，hat越小，说明植被无水分亏缺或不受环境胁迫;hat越大则表明植被受到最大水分亏缺或环境胁迫。由图5可知，3种植物的hat最大值在0.20左右，说明植物在测量时间内受到的水分胁迫程度轻微，主要原因是在测量前研究区经历了两天的降雨过程，土壤水分较为充足。其中，柠条和沙柳的hat值日变化规律为“单峰型”，沙柳则为“双峰型”。3种植物的hat日均值从大到小依次为：沙柳(0.18)>沙棘(0.11)>柠条(0.03)。说明柠条抵抗胁迫的能力最强，沙柳抵抗胁迫的能力最差。

图5 植物蒸腾扩散系数(hat)变化规律及日均值

2.4 常见荒漠植物能量收支特征

地表能量收支是一种陆面过程，包括太阳能和地球内能在地表物质之间的传递输送以及存在形式的转变等。在本研究中能量收支指的是针对植物的能量收支情况，净辐射通量(Rn)是指地表净得的短波辐射和长波辐射的和，是地表能量、动量、水分输送与交换过程中的主要能源。对于植物的蒸腾耗能需求具有重要意义。本研究中净辐射依据太阳辐射计算得到。显热通量(H)是物体或热系统之间的热量交换，不会发生相变，唯一改变的是温度。单位时间内，大气与地表间沿铅直向通过单位面积流过的热量即为显热通量，又称感热通量。对于本研究，显热通量指的是植株向周围环境释放的能量，正值为朝空气中释放热量，负值为吸收空气中热量。本研究显热通量通过能量平衡公式反推得到。由图6可知，柠条净辐射通量在13：00时最高，为463.62 W/m2；沙棘的净辐射通量在11：00，13：00时较高，在456～460 W/m2之间，沙柳的净辐射通量则在11：00最高，为450.24 W/m2。柠条和沙柳的显热通量存在小于0的情况，说明该时刻净辐射通量无法满足植物蒸腾的需要，植物开始从空气中吸收能量来保证蒸腾的正常进行，柠条从外界吸收的热量高于其他两种植物。

图6 3种植物能量收支状况

3 讨 论

不同地区的植物叶片在长期进化的过程中形成了不同的结构来适应当地环境，因此在选择参考叶片时应该因地制宜，选择适合该地区植物的参考叶片。在选取参考叶片时应同时考虑获取的简单性和结果的准确性两大特点，本着获取简单的原则，选择绿色卡纸和便签纸作为参考叶片进行植物蒸腾的计算，但和柠条叶片相比，便签纸厚度较薄，卡纸的厚度和柠条叶片类似。由于参考叶片不具备蒸腾降温作用，其温度变化主要受太阳辐射的影响，当辐射强度较高时，其温度远高于植物叶片。不同厚度的参考叶吸收太阳辐射的能力不同，便签纸温度在正午时分仅略高于叶片温度，和卡纸和干叶片温度相差极大，这一结果和Qiu、凌军等人[11-12]测得的在强太阳辐射下干叶温度高于植物叶的情况不符。而卡纸叶片与干叶片温度在太阳辐射较高时温度变化极为吻合，随着太阳辐射的减小吻合度略有降低但计算结果不存在显著性差异，结合Zhou等人[13]发现的三温模型在9：00—16：00时具有较高的精确度的研究结果，说明当太阳辐射强度较高时卡纸可以替代干叶片进行荒漠地区植物蒸腾速率的计算。

植物在绿化环境的同时也可以降低环境温度。苏从先的研究表明，西北干旱荒漠地区的植物相较于周围环境是一个冷源，形成冷岛效应[14]，植物的冷岛效应可以进一步改善城市热环境，提高生活质量[15]。本试验中3种植物的显热通量存在接近零甚至小于零的情况，表明植物为满足植物蒸腾的需要，从外界环境中吸收热量以满足植物蒸腾的消耗，从而达到降低外界环境温度的目的，与于小惠研究发现草坪净辐射通量低于潜热通量的结果一致[16]。但不同植物的蒸腾水平不同，对水分的消耗和抵御胁迫的程度也有所不同。高永等人研究发现植物蒸腾同株高冠幅间具有很大关系[17]。本研究中沙柳株高及冠幅大于柠条和沙棘，但日蒸腾量和沙棘相仿，低于柠条，表明同水平下沙柳的蒸腾耗水量低于其他两种植物。同时柠条累计蒸腾耗水量最高，说明柠条的吸水能力更强，其适应性也最强。邱国玉等人[18]在温室通过对番茄、甜瓜等作物进行的水分控制试验研究表明，hat值可以反映植物根系水分状况，hat值越接近1受到的水分胁迫越大，本试验中沙柳的hat值平均值高于沙棘和柠条，说明其抵御胁迫的能力最低，柠条抵抗胁迫的能力最强。

4 结 论

(1) 综合参考叶片表面温度变化规律和蒸腾计算结果表明绿色便签纸不适合作为参考叶片进行荒漠植物蒸腾的计算。绿色卡纸可以在9：00—15：00时作为参考叶片进行植物蒸腾的测定。

(2) 沙柳、沙棘、柠条的蒸腾速率日变化规律为单峰曲线，最高瞬时蒸腾速率分别为0.57，0.61和0.71 mm/h，日蒸腾量大小为：柠条>沙棘>沙柳。3种植物蒸腾扩散系数日均值从大到小依次为：沙柳(0.18)>沙棘(0.11)>柠条(0.03)。

(3) 植物可以从外界吸收热量用于蒸腾从而达到降低环境温度的目的，综合3种植物蒸腾量和hat值，在进行生态修复工作时，地势平坦，水分状况较好的地区应优先考虑沙柳；在地势险要、水分分状况较差的区域应优先考虑沙棘和柠条。