-
树木通过边材把根系吸收的水分传输到冠层,即树干液流过程(sap flow),是在个体尺度上分析树木的水分利用特征、诊断树木水分状况的有力工具[1-3],也是估算单株蒸腾量的有力手段之一[4-5]。目前,热示踪技术因能够连续、自动监测到边材液流速率而得到广泛运用[6-7]。然而,除了树木优势度[8]、不同分化等级[9]等树形因子可能引起液流速率在不同个体间的明显差异外,从某一方位测定的边材液流速率是否能够有效地代表整个树干的边材液流传输的平均值,从而为估算单株的液流通量提供可靠依据,仍然没有取得一致的结论,如报道认为西伯利亚落叶松(Larix sibirica Ledeb.)[10]、沙地樟子松(Pinus sylvestris var. mongolica Litv.)均表现出南侧液流速率高[11],而枣树(Ziziphus jujuba Mill.)西向液流速率较高[12],旱柳(Salix matsudana Koidz.)与小叶杨(Populus simonii Carr)方位差异的关系不固定,并随季节而变化等[13]。可见,明确一个树种边材液流速率的方位特征,对于更深入地认识其水分传输特征、提高单株蒸腾量的估算都具有重大意义。
赤松(Pinus densiflora Sieb. et Zucc.)是继樟子松沙地引种培育成功后的又一个优良固沙造林树种[14-15]。对于沙地樟子松,目前的研究已经对其树干液流特征有了较完整的认识[16-18],并发现了液流速率较明显的南、北方位差异[11],但是,对于赤松在这一地区表现的水分利用特征仍很少见报道。作为一个新的优良固沙树种,在大面积推广前迫切需要对其水分利用特征开展系统研究。为此,本研究利用热扩散技术(TDP)对赤松树干4个方位的液流速率开展连续监测,以期对沙地赤松的液流传输特征有较全面的认识,为沙地赤松的生产经营活动提供技术支撑。
-
试验地选择在科尔沁沙地东南缘的辽宁省固沙造林研究所实验林场三家子实验区。地理位置42°40′57" N,122°33′06" E,海拔207 m,年均气温6.3℃,全年无霜期150~160 d,年均降水量475.5 mm(1954—2010年),降水60%70%集中在6—8月,年蒸发量1 553.2 mm,约为降水量的3.27倍。该区土壤以风沙土为主,沙土颗粒均匀,沙层厚度126~128 m。林下草本主要为马唐(Digitaria Sanguinalis L.Scop.)、拂子茅(Calamagrostis epigeos L. Roth.)、碱草(Aneurolepidium chinense Trin. Kitag.)、中华隐子草(Cleistogenes chinensis (Maxim.) Keng)等18种草本,草本总盖度约93%。
-
试验林为辽宁省固沙造林研究所三家子实验区16林班9小班的赤松纯林,林龄40年,面积2.7 hm2。于1977年用2年生裸根苗造林,初植密度1 100株·hm-2,从1977年到2016年该林分经历2次间伐,现保存密度420株·hm-2,平均胸径20.15 cm,平均死枝下高4.2 m,平均活枝下高4.5 m,平均全树高10.9 m,东西冠幅平均5.3 m,南北冠幅平均5.9 m。根据林分调查结果,在样地内选取树干通直、生长良好、无病虫害的3株赤松作为树干液流测定样木,样木情况见表 1。
表 1 样木基本情况
Table 1. Information of sample trees for measuring sap flow density
编号
Tree number胸径
DBH/cm死枝下高
Dead branch
height /m活枝下高
Live branch
height/m全树高
Tree height/
m冠幅Crown diameter/m 东侧
East side南侧
South side西侧
West side北侧
North sideTree 1 21.3 4.5 5.4 11.2 3.0 3.2 2.6 2.6 Tree 2 20.9 3.9 5.2 9.5 2.9 3.0 2.5 3.1 Tree 3 25.1 - 3.6 10.8 3.1 3.6 3.4 2.7 -
应用FLGS-TDP植物液流计(Dynamax,U.S.A),于2017年3月分别在选定样株树干的东、南、西、北各安装1套TDP液流计30 mm长度探针传感器,具体安装过程及防护措施见文献[19-20]。数据采集器为CR1000(Campbell Scientific,U.S.A),采样间隔为10 min,每隔30 min记录1次数据的平均值。选取4—10月的数据代表 1个完整生长季进行数据分析。
液流速率(Js)的计算公式如下:
$ {J_s} = {\rm{ }}0.011{\rm{ }}9{\rm{ }}{K^{1.231}} \times 3{\rm{ }}{600^{[21]}} $
式中:Js为每小时的平均液流速率(cm·h-1);K = (dTm- dT)/dT,其中dTm为无液流时加热探针与参考探针的最大温差值,dT为瞬时温差值。dTm的确定以10天为1个周期,采用两次回归法确定[7, 20]。
-
利用便携式土壤水分测量仪(TRIME-Pico 64/32 TDR,IMKO,德国)监测土壤体积含水率。于2017年3月末在样地内随机埋设TDR测管3处,每处在0~180 cm层次内按10、20、50、100、150、180 cm共6个层次分层监测。每月10、20、30日为观测时间。
-
利用章古台生态站CAWS600-B型自动气象站(华云公司,中国)监测的太阳总辐射、降水量数据。
-
4—10月试验期间共观测到21次有效降水(单次降水>2.5 mm的次数),降水总量为335.8 mm,占该地区同期平均降水量(448.0 mm,1983—2017年4—10月间)的75%,其中约75%的降水集中在7、8月份,年内分布极不均匀(图 1a)。受降水的影响,10~20 cm土层的平均体积含水率(SWC, %)表现出与降水过程相一致的季节变化规律,而50 cm以下土层的SWC仅对于70.4 mm(8月3日)的强降水有响应,之后对降水不再有响应。从季节动态来看,4—6月份为水分衰减期,SWC变化在6.18%~9.85%范围内。7—8月份,随着降水的集中,土壤水分明显上升,SWC最高达17.9%,9月份之后又进入衰减期。试验期间(4—10月份)各月的SWC(10~180 cm平均值)分别为8.82%、7.03%、6.21%,7.56%,11.73%、8.92%、7.30%,8、9月份的值较高,6月份最小(图 1b)。受天气变化的影响,太阳总辐射(Rs,W·m-2)的日平均值波动较大,但自8月份之后Rs表现出较明显的递减趋势(图 1c)。
不同方位边材液流速率(Js,cm·d-1)总体上表现出较一致的季节变化动态。从整个生长季来看,Js的排序为南侧(97.7 ± 47.7 cm·d-1)>西侧(90.2 ± 49.6 cm·d-1)>北侧(85.7 ± 50.8 cm·d-1)>东侧(63.5 ± 32.1 cm·d-1),其中北侧液流速率是4个方位平均值的1.01倍,说明只测定北侧的液流速率就能较好地代表 4个方位的液流。
-
为了进一步比较不同方位Js的日变化进程及随季节的变化,对整个生长季逐月气象条件筛查,每月挑选出一典型晴天(确保该日数天前也为晴日),比较其不同方位Js的特征。就整个生长季节而言,不同方位Js的日变化进程有着不同的格局(图 2a),表现在南侧液流峰值出现早且峰值大(11.5 cm·h-1),其次为西侧液流峰值出现的时刻明显滞后。北侧液流表现出较典型的正态分布型,而东侧液流的峰值最低(9.68 cm·h-1)。就不同月份方位特征的比较来看,除了7、8月份西侧Js略高于南侧外,一般均以南侧Js最高。赤松液流传输主要来自于南、西、北3个方位,东侧Js相对较低(图 2a)。
图 2 不同方位液流速率的日变化过程,峰值及峰值时刻的季节性差异
Figure 2. Diurnal course, peak value and peak time of sap flow density from four directions of the trunk over a growing season.
一日内液流峰值出现的时刻均以东侧最早,其次为南侧或北侧,而以西侧为最晚。在生长季初期的4月份,东、南、西、北各方位液流速率的峰值时间均比较早(图 2b),在11:30—12:00之间,而在干旱月份的5—6月份,峰值出现的时刻明显滞后,如在6月份,东、南、西、北各方位液流速率的峰值时间在16:00—17:30(图 2d),在降雨充沛的7—8月,峰值时刻又趋于回归,如在8月份,东、南、西、北各方位液流速率的峰值时间在14:00—15:00(图 2f)。尽管如此,各方位液流速率峰值出现时刻差异不显著(P = 0.53)。
对不同方位液流速率以天为时间单位逐月进行统计结果表明,在绝大部分季节,不同方位液流速率以南侧最大,其次依次为西、北、东侧(表 2)。除10月份不同方位间差异不显著外,其余月份差异均显著(p = 0.05水平)。
表 2 不同方位液流速率差异统计
Table 2. Differences of sap flow density in different direction
月份
Month液流速率百分比Percentage/% 统计结果Statistical results 南侧
South
side西侧
West
side北侧
North
side东侧
East
side变异系数Cv 均方根
RMSP值
P value4 100 92 78 59 0.30 30.01 < 0.01 5 100 83 68 60 0.45 32.68 < 0.01 6 100 71 59 67 0.32 13.59 < 0.01 7 100 104 108 66 0.37 41.98 < 0.01 8 100 97 98 69 0.48 53.17 =0.01 9 100 96 99 69 0.33 35.28 < 0.01 10 100 86 81 66 0.58 24.96 =0.06 -
以旬为时间尺度的比较结果表明,4—6月,各方位Js的旬平均值遵循相似的动态变化进程,且具有较稳定的排序:南侧>西侧>北侧>东侧。其中4月份Js呈快速上升(图 3a),到下旬达到最高值:南侧(14.3 m·旬-1)>西侧(13.1 m·旬-1)>北侧(11.5 m·旬-1)>东侧(8.8 m·旬-1)。5—6月份Js快速下降,其间降水强度分别为3.8 mm(5月4日)、5.2 mm(5月5日)、6.2 mm(5月22)、7.4 mm(6月18日)、4.2 mm (6月19日)的次降雨事件只带来了液流的短暂回升(仅持续23天)(图 1d)。经过1个月的干旱期,北侧Js最低(图 3b)。7月份,随着几次强降雨事件(如30.7 mm,7月1日;23.8 mm,7月7日;17.9 mm,7月13日;26.4 mm,7月20日),Js快速恢复,其中北侧Js恢复强度最大(图 3c),到7月下旬北侧Js平均为15.4 m·旬-1。8月份,随着降雨更加频繁,Js能达到的峰值基本稳定,只是受天气的影响,各旬实际的液流速率间波动较大(图 3d)。可以看出,随着土壤水分状况的改善,南、北、西3个方位的Js值更高且更接近,与东侧Js间的差距变大。9月中旬后,随着气象蒸发条件的快速减弱,各方位Js均下降明显,进入消退期(图 3e),直至10月下旬的生长末季,各方位间Js比较接近(1.9~2.9 m·旬-1)。
-
将整个生长季划分为干旱季(5—6月份)与非干旱季,分别对不同方位液流速率与太阳辐射进行相关分析,结果表明:在非干旱期,各方位的Js对Rs的增加具有紧密的协同响应关系(Guass函数关系),表现出太阳辐射对树木蒸腾过程较强的驱动效应。但这一关系也具有阈值,当Rs升高到约287.4~313.9 W·m-2范围后,随着Rs的进一步增加,Js表现出下降趋势(图 4)。在干旱期,Js明显低于相同辐射强度下的非干旱期的Js,且与Rs间的关系不明确,反映了土壤水分亏缺对大气驱动蒸腾过程的抑制效应。
-
土壤水分对乔木树种液流过程的影响一般仅在较大的时间尺度上才能显现出来。试验期间的4—6月份,获得的降水量只占同期降水均值(1983—2017年)的23.0%,为历年该时段所占比例最低值(范围为23%~156%),赤松经历了较严重的水分胁迫,在6月末时,SWC低至6.2%,林下草本都已经处于萎蔫状态。尽管从图 1d和图 3b已经直观地展现了Js对水分胁迫的响应,进一步通过构造变量Js/Rs建立其与SWC之间的关系,结果表明:赤松4个方位的Js/Rs具有对SWC较一致的响应规律,当SWC下降到一定值后,Js/Rs快速下降,其中南、西、北3侧当SWC下降到约为7.5%8.0%时Js/Rs快速下降(图 5a~5c),而对于东侧,这一值约为8.0%~8.5%(图 5d)。
沙地赤松树干边材液流速率的方位特征研究
Azimuthal Variation of the Sap Flow Rate in the Trunk of Pinus densiflora Planted in Sandy Land
-
摘要:
目的 确定优良固沙树种沙地赤松(Pinus densiflora)树干边材液流速率(Js)方位变化规律,探讨环境因素对方位差异的影响。 方法 采用热扩散式树干液流计(TDP)连续监测树干液流并同步监测环境因素,比较分析不同方位树干液流速率的差异、季节动态及与降水、土壤体积含水率(θ)、太阳辐射强度(Rs)间的关系。 结果 表明:(1)在赤松树干东、南、西、北4个方位上,Js的日变化与季节变化均呈现出基本一致的格局,且均与Rs间表现出紧密的协同变化关系;(2)在中等强度干旱和极端干旱情况下,各方位边材Js均同步受到明显抑制,在土壤水分得到充分补充后又快速回升;(3)在典型晴日里,树干各方位边材Js其午间峰值出现的时刻有差异,呈从早到晚分别是东侧、南侧、西侧、北侧这种顺时针的方位规律;(4)在整个生长季,Js的日平均值大小一般为南侧>西侧>北侧>东侧,但方位间差异并不显著(P=0.35),北侧Js的平均值最接近4个方位的平均值(约为平均值的1.01倍)。 结论 为降低观测成本,通过北侧一个方位的测定来估算赤松单株的液流通量具有较好的可靠性;影响赤松各方位液流过程明显变化的θ的阈值为7.5%。 -
关键词:
- 液流速率
- / 赤松
- / 方位
- / 水分利用
- / 热扩散技术(TDP)
Abstract:Objective To determine the azimuthal variation of the sap flow rate of trunks (Js) and estimate the effects of environmental factors on azimuthal variation. Method The sap flow rate of Pinus densiflora, an important windbreak and sand fixation species, was monitored continuously during the growing season using a thermal dissipation probe (TDP) method, combined with simultaneous observations of precipitation, soil moisture, and solar radiation and other environmental factors. Result The results of comprehensive analysis showed that the daily and seasonal changes of the Js in the four directions of the trunk basically followed the same pattern, and all were in accordance with the solar radiation. In moderate-intensity drought and extreme drought conditions, the Js of each direction was significantly suppressed simultaneously and quickly recovered after the soil moisture was fully replenished. On a typical sunny day, there were differences in the time at which the peaks of the Js appeared at various times. From early to later, it ranked with the order of the east, south, west, and north sides clockwisely. During the whole growing season, the daily mean value of Js was generally ranked of south > west > north > east, but the difference between directions was not significant (P=0.35). Statistics show that the average of the Js in the north side was the closest to the average of four azimuths (approximately 1.01 times the average). Conclusion The sap flow estimated through a measurement on the north side in the actual measurement will get better reliability in order to reduce the cost of observation. Meanwhile, the threshold value of soil moisture affecting the sap flow of Pinus densiflora is about 7.5%. -
Key words:
- sap flow rate
- / Pinus densiflora
- / azimuth
- / water use
- / Thermal Dissipation Probe method (TDP)
-
表 1 样木基本情况
Table 1. Information of sample trees for measuring sap flow density
编号
Tree number胸径
DBH/cm死枝下高
Dead branch
height /m活枝下高
Live branch
height/m全树高
Tree height/
m冠幅Crown diameter/m 东侧
East side南侧
South side西侧
West side北侧
North sideTree 1 21.3 4.5 5.4 11.2 3.0 3.2 2.6 2.6 Tree 2 20.9 3.9 5.2 9.5 2.9 3.0 2.5 3.1 Tree 3 25.1 - 3.6 10.8 3.1 3.6 3.4 2.7 表 2 不同方位液流速率差异统计
Table 2. Differences of sap flow density in different direction
月份
Month液流速率百分比Percentage/% 统计结果Statistical results 南侧
South
side西侧
West
side北侧
North
side东侧
East
side变异系数Cv 均方根
RMSP值
P value4 100 92 78 59 0.30 30.01 < 0.01 5 100 83 68 60 0.45 32.68 < 0.01 6 100 71 59 67 0.32 13.59 < 0.01 7 100 104 108 66 0.37 41.98 < 0.01 8 100 97 98 69 0.48 53.17 =0.01 9 100 96 99 69 0.33 35.28 < 0.01 10 100 86 81 66 0.58 24.96 =0.06 -
[1] Chirino E, Bellot J, Sanchez J R. Daily sap flow rate as an indicator of drought avoidance mechanisms in five Mediterranean perennial species in semi-arid southeastern Spain[J]. Trees, 2011, 25(4):593-606. doi: 10.1007/s00468-010-0536-4 [2] Nadezhdina N. Sap flow index as an indicator of plant water status[J]. Tree Physiology, 1999, 19(13):885-891. doi: 10.1093/treephys/19.13.885 [3] Børja I, Světlík J, Nadezhdin V, et al. Sap flux-a real time assessment of health status in Norway spruce[J]. Scandinavian Journal of Forest Research, 2015, 31(5):450-457. [4] Gong D Z, Kang S Z, Yao L M, et al. Estimation of evapotranspiration and its components from an apple orchard in northwest China using sap flow and water balance methods[J]. Hydrological Processes, 2007, 21(7):931-938. doi: 10.1002/(ISSN)1099-1085 [5] Whitley R, Medlyn B, Zeppel M, et al. Comparing the Penman-Monteith equation and a modified Jarvis-Stewart model with an artificial neural network to estimate stand-scale transpiration and canopy conductance[J]. Journal of Hydrology, 2009, 373:256-266. doi: 10.1016/j.jhydrol.2009.04.036 [6] Lu P, Urban L, Zhao P. Granier's thermal dissipation probe (TDP) method for measuring sap flow in trees:theory and practice[J]. Acta Botanica Sinica, 2004, 46(6):631-646. [7] 党宏忠, 张劲松, 赵雨森.应用热扩散技术对柠条锦鸡儿主根液流速率的研究[J].林业科学, 2010, 46(3):29-36. [8] 李振华, 王彦辉, 于澎涛, 等.华北落叶松液流速率的优势度差异及其对林分蒸腾估计的影响[J].林业科学研究. 2015, 28(1):8-16. [9] 刘家霖, 满秀玲, 胡悦.兴安落叶松天然林不同分化等级林木树干液流对综合环境因子的响应[J].林业科学研究. 2016, 29(5):726-734. doi: 10.3969/j.issn.1001-1498.2016.05.015 [10] 白志强, 刘华, 佘春燕, 等.西伯利亚落叶松树干液流的动态变化[J].河北农业大学学报, 2016, 39(3):49-54. [11] 张友焱, 周泽福, 党宏忠, 等.利用TDP茎流计研究沙地樟子松的树干液流[J].水土保持研究, 2006, 13(4):78-80. doi: 10.3969/j.issn.1005-3409.2006.04.024 [12] 马建鹏, 汪有科, 陈滇豫, 等.不同时间尺度上枣树树干液流的变异特性[J].干旱地区农业研究, 2016, 34(3):95-101. [13] 徐丹丹, 尹立河, 侯光才, 等.毛乌素沙地旱柳和小叶杨树干液流密度及其与气象因子的关系[J].干旱区研究, 2017, 34(2):375-382. [14] 孟鹏.章古台沙地37年生赤松和樟子松生长特性研究[J].辽宁林业科技, 2013(5):20-23, 44. doi: 10.3969/j.issn.1001-1714.2013.05.005 [15] 张日升.章古台沙地3种针叶树生长规律对比研究[J].安徽农业科学, 2017, 45(18):140-142, 175. doi: 10.3969/j.issn.0517-6611.2017.18.045 [16] 韩辉, 张学利, 党宏忠, 等.基于树干液流通量的沙地樟子松合理林分密度的确定[J].林业科学研究, 2015, 28(6):797-803. doi: 10.3969/j.issn.1001-1498.2015.06.006 [17] 韩辉, 白雪峰, 徐贵军, 等.章古台樟子松树干液流的密度特征[J].东北林业大学学报, 2013, 41(4):27-31+82. [18] 卢志朋, 魏亚伟, 李志远, 等.辽西北沙地樟子松树干液流的变化特征及其影响因素[J].生态学杂志, 2017, 36(11):3182-3189. [19] Dang H, Zha T, Zhang J, et al. Radial profile of sap flow velocity in mature Xinjiang poplar (Populus alba L. var. pyramidalis) in Northwest China[J]. Journal of Arid Land, 2014, 6(5):612-627. doi: 10.1007/s40333-014-0007-7 [20] Lu P. A direct method for estimating the average sap flux density using a modified Granier measuring system[J]. Functional Plant Biology, 2001, 24(28):701-705. [21] Granier A. Evaluation of transpiration in a Douglas-fir stand by means of sap flow measurements[J]. Tree Physiology, 1987, 3(4):309-320. [22] Lu P, Woo K C, Liu Z T. Estimation of whole-plant transpiration of bananas using sap flow measurements[J]. Journal of Experimental Botany, 2002, 53(375):1771-1779. doi: 10.1093/jxb/erf019 [23] 孟秦倩, 王健, 张青峰, 等.黄土山地苹果树树体不同方位液流速率分析[J].生态学报, 2013, 33(11):3555-3561. [24] Naithani K J, Ewers B E, Pendall E. Sap flux-scaled transpiration and stomatal conductance response to soil and atmospheric drought in a semi-arid sagebrush ecosystem[J]. Journal of Hydrology, 2012, 464-465:176-185. doi: 10.1016/j.jhydrol.2012.07.008 [25] Clausnitzer F, Koestner B, Schwaerzel K, et al. Relationships between canopy transpiration, atmospheric conditions and soil water availability-Analyses of long-term sap-flow measurements in an old Norway spruce forest at the Ore Mountains/Germany[J]. Agricultural and Forest Meteorology, 2011, 151(8):1023-1034. doi: 10.1016/j.agrformet.2011.04.007 [26] 姚依强, 陈珂, 王彦辉, 等.华北落叶松树干液流速率主要影响因子及关系的时间尺度变化[J].干旱区资源与环境, 2017, 31(2):155-161. [27] 万艳芳, 于澎涛, 刘贤德, 等.祁连山青海云杉树干液流密度的优势度差异[J].生态学报, 2017, 37(9):3106-3114. [28] 李浩, 胡顺军, 朱海, 等.基于热扩散技术的梭梭树干液流特征研究[J].生态学报, 2017, 37(21):7187-7196. [29] 张晓艳, 褚建民, 孟平, 等.民勤绿洲荒漠过渡带梭梭(Haloxylon ammodendron(C.A.Mey) Bunge)树干液流特征及其对环境因子的响应[J].生态学报, 2017, 37(5):1525-1536. [30] 党宏忠, 杨文斌, 李卫, 等.新疆杨树干液流的径向变化及时滞特征[J].生态学报, 2015, 35(15):5110-5120. [31] 胡广录, 赵文智, 谢国勋.干旱区植被生态需水理论研究进展[J].地球科学进展, 2008, 23(2):87-94. [32] 焦树仁.辽宁省章古台樟子松固沙林提早衰弱的原因与防治措施[J].林业科学, 2001, 37(2):131-138. [33] 焦树仁.辽宁章古台樟子松人工林水分动态的研究[J].植物生态学与地植物学学报, 1987, 11(4):296-307.