王喜娜 何嘉茜
(1.广东生态工程职业学院,广东 广州 510520;2.广东碧桂园学校,广东 佛山 528312)
蒸散发是全球水、能量和碳循环的重要组成部分和关键环节,与土壤-植物-大气系统密切相关,是地球表层中能量循环和物质转化最为强烈的活动层[1]。蒸散发过程涉及到土壤水文地质和近地大气紊流交换以及植物的生理活动,对局地微气候产生重要影响,是植物生长和生物产量的主要控制因素之一[2]。植物的蒸散发与园林、农业、林业和水资源管理以及气象监测预警等密切相关,因此,对植物蒸散发的研究和耗水监测分析有着重要的学术意义与社会经济价值。
植物蒸散是一个复杂的物理现象和生理过程,包括从地表和植物表面的水分蒸发,以及通过植物表面和植物体内的水分蒸散[1]。植物根系依靠渗透压从水中或土中吸收水分后输送到叶片, 然后经植物杆径、叶片以及花果表面蒸发到大气, 称为植物散发或蒸腾[3]。植物叶片的蒸腾作用可以帮助植物从土壤吸收水分,促进光合作用;同时也可以使植物排出多余水分,降低叶面温度,以适应其生长和生理需求。影响植物蒸散的因素很多,主要包括植物生理特性、气温、光照、土壤含水量、大气相对湿度和风速等[3]。在相同的外部条件下,不同的植物种类以及同种植物的不同生长期,其蒸散耗水量也各不相同[4-5]。比如,1995—1997 年间,鼎湖山人工松林的年均降雨量和蒸散量是1 878 mm 和938 mm,7 月蒸散量最大,高达149 mm,最小为1 月只有22 mm[6]。精确地监测和测算各个物种在不同生长期的蒸散量,是区域农业物种选择、城市生态景观规划和水资源精细化管理的基础,可以帮助决策者根据不同的水资源等外部环境条件选择种植不同的物种,或者根据不同的植物种类和生长周期辅以相应的浇灌策略。陈友媛等[7]于2011 年7—11 月份对对春羽Philodenron selloum、绿萝Scindapsus aureun和铜钱草Hydrocotyle vulgaris3 种水培盆栽景观植物的蒸散量开展了定量监测,发现铜钱草的耗水量高达1 550 g.m-2.d-1,是同期春羽和绿萝耗水量的2.1倍和3.6 倍。
盆栽花卉是城市生态景观与园林建设的主要内容之一,具有广阔的市场容量和巨大的社会经济价值与美学价值[8]。例如,长寿花Narcissus jonquilla长期占据丹麦盆栽花卉销售产值之冠,而荷兰的长寿花产值已达千亿美元[9]。城市家庭的盆栽花卉养殖已成为市民生活和城市景观的重要组成部分;然而,我国在盆栽花卉养殖技术研究中,除了如王梓等[10]和陈友媛等[7]研究外,对花卉的蒸散耗水和浇水量多数以定性描述为主,缺乏对花卉蒸散量和需水量的定量化监测[11-12]。
植物蒸散量的监测往往费时费力,成本较高,难以大规模的实施。多种方法已被开发出来估算区域或田块尺度的蒸散量,诸如空气饱和差法、气温积温法、Penman 公式法和布德科公式法[6]。这些常用的植被蒸散量测算模型多用于区域尺度和田块尺度的蒸散量测算[13-14],而对于单株城市盆栽花卉的蒸散量计算的研究和应用并不多见。本文的主要研究目的是选择几种常见的亚热带城市盆栽花卉,每天监测其蒸散量,分析典型盆栽花卉在不同季节的蒸散耗水量变化及其主要影响因素,并改进空气饱和差法,为快速测算盆栽花卉的潜在蒸散量和自动浇灌设计等提供参考。
本文研究对象为广州市城市家庭常见的4 种成年盆栽花卉:铜钱草(HV)、长寿花(NJ)、碧玉Peperomia tetraphylla(PT)和蝴蝶兰Phalaenopsis hybrid(PH)(图1)。其生境为广州市城区高层住宅阳台花架,阳台朝南,外有挡雨窗,免受雨水直淋,采光和透气性良好,栽培方式为盆栽土养,人工根部浇水。
图1 研究花卉与监测设备Figure 1 The study plants and instruments
铜钱草属伞形科和多年生挺水或湿生草本植物,株高可达15 cm,地下走茎发达,生长适合温度20~25 ℃[15-17]。长寿花属多年生草本植物,植株可高达40 cm,喜温暖湿润和阳光充足环境,生长适温为15~25℃[11]。碧玉属多年生簇生草本植物,植株高10~30 cm,喜湿润半阴环境,生长适温25 ℃左右,忌阳光直射[18-20]。蝴蝶兰属兰科多年附生草本植物,主要分布在热带及亚热带地区,喜暖畏寒,生长适温为16~27 ℃,10 ℃以下停止生长[12,21]。
常见的植物蒸散耗水量测定和估算方法有水文学法、微气象学法、植物生理学法、遥感方法以及综合模拟法,其中实测法有称重法、气孔计法、水文学法、能量平衡法等。针对家庭阳台盆栽花卉体积小和不受雨水直淋的特点,本文采用电子秤称重法直接测定盆栽花卉的蒸散量。电子秤由中山康诺电子有限公司生产,型号H-X,最大量程5 000 g,精度为1 g。同时,监测水杯内的水面蒸发量,用以表征整个大气环境和太阳辐射能量的实际蒸发潜力。水杯放在阳台外的窗台上,盆栽花卉放在通风良好的阳台内花架上,每个花盆下都配有接水托盘,水杯和花盆都不受雨水直淋;每天人工对花卉根部浇水,根据前1 d 的蒸散量控制浇水量,避免托盘积水,用称重法直接测量浇水前后的花盆质量,二者质量差即为当天的浇水量。在每次称重前,记录当时的温度和相对湿度,温度湿度计放在阳台内的花盆旁边,免受阳光直射;温度湿度计由温州盛测仪器仪表有限公司生产,型号为HTC-8(图1),温度测量范围:-50~70 ℃,测量精度:±1 ℃,分辨率:0.1 ℃;相对湿度测量范围:10%~99%,测量精度:±5%,分辨率:1%。
2019 年11 月1 日—2020 年10 月31 日,每天上午约7:30 时、晚上约19:30 时分别对花盆和托盘整体称重,计算出早、晚的质量差(ΔS),建立花盆(含托盘)的水量平衡公式(1):
公式(1)中,ΔS为花盆在早、晚两个时间段内的质量差,P为花盆接收到的降雨量,W为人工浇水量,ET 为花盆的蒸散量,Q为渗流量。
每次对花卉浇水都在晚上花盆称重后进行,先测得ΔS的值,此时浇水量(W)为零;花盆放在阳台内不受雨淋,降雨量(P)为零;每次浇水量都略等于前1 d 的蒸散量,避免过量浇水导致下渗,且花盆底放置有收集渗水的托盘,托盘和花盆一起称重,即使偶尔产生的少许渗流也被计入到花盆的蒸发量当中,这样渗流量(Q)也为零。根据花盆的水量平衡公式(1),当P、W和Q都为零时,早、晚两个时间段内花盆(包括托盘)的质量差值就代表每盆花卉的蒸散量(ET)。
根据表1 所列的各个花盆的直径(面积)和水的密度,可将每个花盆的蒸散质量转化为单位面积的蒸散量深度(mm)。花卉每日的蒸散量代表其每日所需的耗水量。由于花卉在12 h 内的生长量远小于花卉的蒸散量,在公式(1)以及本文的蒸散量和耗水量分析中,未单独考虑花卉的生长量。除了测量花盆的面积外,本文还测量了4盆花卉所有的叶片的总面积,计算出单位叶面积的年均蒸散量。本文所用风速和降雨量数据来自广州市气象局国家标准气象站(天河)。
表1 花盆直径、总叶面积及其监测期内(2019-11-1—2020-10-31)的蒸散量Table 1 The pot diameter, total leaf areas and the mean ET during 2019-11-1—2020-10-31
空气饱和差法根据月平均的气温(T)和相对湿度(RH)计算田块尺度的月蒸散量(ET),其通用公式为:
公式(2)中,k是转换系数,a和b表征植被ET 受气温和相对湿度的影响敏感度,T0为温度调节参数。
根据4 种盆栽花卉的实测蒸散量对气温和相对湿度影响的敏感性,将公式(2)改为公式(3):
根据植被蒸散量的连续监测显示,当气温低于植物的适宜生长温度时,植被的日蒸散量极低;而当气温高于某一温度时,植被的日蒸散量急剧上升。因此公式(3)中,增加一个温度参数Tg,可代表最适宜植被生长的气温下限,并将气温影响指数a移到内部作为实际气温与最佳生长气温比值的指数。每种花卉的T0、Tg、a和b4 个参数都用试差法,根据实测的 ET、T和RH 进行敏感性分析与最大相关性交叉拟合,找到一套最优的拟合参数范围,最后确定转换系数k。
将1 a 内每天连续监测的T、RH 和ET 数据按奇、偶天数分为两组。奇数组数据用于率定模型参数(T0、Tg、a和b),通过敏感性分析和试差法,计算ET和公式(3)右边除参数k以外的因子之间的Pearson 相关性系数,根据相关性最高确定最优参数组合。偶数组数据用于模型验证,并用估算的ET 与实测ET 之间的Pearson 相关性系数和相对偏差评估模型的性能。
水杯水面蒸发量表征花卉养殖环境的最大蒸发潜力。从2019 年11 月1 日至2020 年10 月31日,水杯水面累计蒸发7 243 g(或2 251 mm),日均蒸发6.2 mm,其中白天蒸发量是晚上的2.8倍(表1);日蒸发量最大值出现在2020 年7 月17 日—19 日(日最高温37 ℃,RH=60%,并伴随3~4 级东南风),高达16.2 mm。水面蒸发量与气温成正相关,与RH 成负相关;白天的蒸发量与RH 的相关性比晚上大(图2)。根据日均气温和RH,改进后的空气饱和差法能较好地估算出水面的潜在蒸发量(R2=0.50, RD = -3.6%)(图2)。
图2 水杯水面蒸发量与气温(a,b)和RH(c, d)间的散点图(2019/11/1-2020/10/31),以及模型估算蒸发量与实测蒸发量之间的散点图(e,f)Figure 2 The scatter plots of the half-day glass cup evaporation against air temperature (a,b)and RH(c, d),and the estimated evaporation against the measurements (e,f)
从月尺度看(表2),日均蒸发量的最大值(7月)和最小值(3 月)发生月份分别对应气温和RH 的最大值,说明在夏季气温对水面最大蒸发量起着主导作用,比如11 月和12 月份的日均RH 相同,但11 月的日均气温比12 月高4.6 ℃,11 月份的日均蒸发量也比12 月份的大1.6 mm。
表2 气温、相对湿度、水杯蒸发量和4 种花卉蒸散量每月的日均值,2019-11-1—2020-10-31Table 2 The daily mean air temperature, RH, glass evaporation and plant ET in each month during 2019-11-1—2020-10-31
在水面蒸发量中,由于住宅阳台环境条件限制,本研究采用直径为6.4 cm 的两个玻璃杯盛水后平行监测,取其平均值。但是,水杯的直径远小于常规的气象或者水文监测中用的D20(直径20 cm)和E601(61.8 cm)蒸发皿,所监测到的水面蒸发量大于自然开阔水面的蒸发量,如广州市地区2000—2012 年的年均水面蒸发量为1 150 mm[22],而本文在2019 年11 月1 日至2020 年10 月31 日期间用玻璃杯监测到的累计蒸发量高达2 251 mm,这主要是由于杯壁表面及其毛细作用对水面蒸发量的影响。本文增加了1个20.4 cm 口径的蒸发皿代表D20 蒸发皿的蒸发量,与水杯一起在2020 年5 月16 日至6 月17 日连续监测1 个月。玻璃杯的水面蒸发量与D20 水面蒸发量之间有很好的线性相关(CC=0.93),它与D20 水面蒸发量之间的折算系数为0.63。广东地区采用D20 蒸发皿测得的蒸发量与自然水面蒸发量之间的折算系数在水文行业的推荐值为0.69[22];因此,本文监测的玻璃杯水面蒸发量与开阔水面蒸发量的折算系数为0.63×0.69 = 0.44。在这365 d 的监测期,水面的自然蒸发量约等于2 251×0.44 =991 mm(2.7 mm/d),略低于2000—2012 年广州市的平均水面自然蒸发量1 150 mm,表明我们用水杯监测到的蒸发量及其与开阔水面自然蒸发量之间的折算系数是合理的。
本文所指的花卉蒸散量为花卉所在土壤(轻基质)蒸发量和花卉蒸散量之和。在这1 a 的监测期,玻璃杯的水面蒸发量明显高于铜钱草、长寿花、碧玉和蝴蝶兰的日均蒸散量(图3),分别为6.2、3.5、2.5、2.5 和1.8 mm/d,4 种盆栽花卉的蒸散量只相当于玻璃杯水面蒸发量的56%、41%、40%和30%(表1)。4 种花卉中,铜钱草蒸散耗水量最大,比开阔水面的自然蒸发量(2.7 mm/d)高28%,接近耗水量最小的蝴蝶兰的2 倍;在温度较高和相对湿度较低的9—11 月,铜钱草的日均蒸散量甚至达到蝴蝶兰日均蒸散量的3 倍(表2);这与陈友媛等[7]的研究结果一致,铜钱草在7 至11 月份的日均耗水量是春羽和绿萝耗水量的2.1 和3.6 倍;因此,铜钱草不是缺水城市适宜栽种的景观植物。
图3 日均气温、相对湿度、水杯日均蒸发量、4 种盆栽花卉的日均蒸散量Figure 3 The daily mean air T, RH and the daily ET for Glass cup, HV, NJ, PT, PH
除了计算单位花盆面积的蒸散量,本文还测量了4 种花卉的总叶面积,计算出单位叶面积的年蒸散量(表1)。与图1 所示的4 种花卉的生长形态一致,4 种花卉中长寿花总叶面积与花盆面积之比最高(52.8)、蝴蝶兰最低(5.5);长寿花的单位叶面积的年蒸散量最低,铜钱草的单位叶面积年蒸散量依旧最高、是长寿花的5 倍,而蝴蝶兰的单位叶面积年蒸散量排第二、是长寿花的4倍和碧玉的1.7 倍,这与根据单位花盆面积计算的年蒸散量有较大差异。由于各种花卉的生长形态和总的叶面积随季节略有不同,如何便捷的计算和比较不同物种的蒸散量还值得进一步探讨,本文余下内容主要使用单位花盆面积的蒸散量进行讨论。
在温度较低的2 月份以及RH 较高的3—4 月份,4 种盆栽花卉之间的蒸散量非常接近,而碧玉在此期间的蒸散量却相对较高(图3 和表2)。4 种花卉在白天的蒸散量都与RH 呈显著的负相关(95%置信度),而碧玉和蝴蝶兰在晚上的蒸散量与RH 的相关性不显著(表3)。令人意外的是,长寿花在白天和晚上的蒸散量都与气温没有显著的相关性,主要受到RH 的影响;而铜钱草、碧玉和蝴蝶兰在白天和晚上的蒸散量都与气温呈显著的正相关(表3)。
表3 气温、RH 以及水杯和4 种盆栽花卉在白天和晚上的蒸散量均值及其相互间的Pearson 线性相关系数矩阵Table 3 The Pearson Correlation Coefficient (CC) matrix during day time and night time
除了气温和RH 以外,风速对花卉的蒸散量影响也较大。如图3 所示,水杯、铜钱草和长寿花的蒸散量偶尔会出现较陡峭的峰值,这些峰值日期都有较大的风速。例如2019 年12 月27 日出现风速超过10 m/s 的阵风,日均气温和RH 为15.6 ℃和49%,水杯和4 种花的蒸散量分别是11.5、4.1、3.9、2.9 和2.9 mm/d;而无阵风的12月28 日,日均气温和RH 为16.7 ℃和46%,水杯和4 种花的蒸散量分别是8.2、2.2、2.8、1.7 和1.8 mm/d。由于风速的瞬时波动大,本文没有系统的同步监测风速,尚无法定量分析风速对花卉蒸散量的影响。
从2019 年10 月12 日广州市下了15 mm 的降雨后,直到12 月18 和19 日才下了5.4 mm 和4.3 mm 的降雨;连续67 d 无降雨,这是广州市比较罕见的干旱期。当每天晚上对4 种植物持续浇水并使花盆土壤接近饱和,使每种花卉都达到其近似的最大蒸散量,铜钱草的蒸散量是碧玉和蝴蝶兰的3 倍,并随RH 和温度大幅度波动。
从12 月3 日到12 日,连续10 d 的日均RH低于40%。在12 月6 日晚上对每盆花卉浇到饱和后,进行连续9 d 的水压力测试,直到12 月15 日晚上,当铜钱草和长寿花出现严重枯萎时,才恢复对各盆花浇水。铜钱草、长寿花、碧玉和蝴蝶兰在这9 d 的日均蒸散量只相当于12 月份其余日期的日均蒸散量的57%、74%、66%和39%;当水分受限后,各种花卉的蒸散量都趋于雷同,但铜钱草和长寿花的大部分叶片出现枯萎,而碧玉和蝴蝶兰的生长状态未发生明显枯萎,耐旱能力较强(图3)。
另一个特殊时期是在RH 很高的阴雨连绵日子,如2020 年3 月18 日至4 月7 日,广州市连续下雨,累计雨量218 mm,这20 d 的平均气温和RH 为20.3 ℃和85%,玻璃水杯、铜钱草、长寿花、碧玉和蝴蝶兰的日均蒸散量只有2.3、0.9、0.9、1.1 和0.7 mm/d,只相当于整个监测期内各自日均值的37%、26%、36%、45%和38%。这与干旱期水分受限时的蒸散量类似,当空气RH 较高时,各种花卉之间的蒸发量急剧降低且近乎相等。
根据上述4 种盆栽花卉的实测蒸散量受气温和RH 的不同影响,我们改进了空气饱和差法模型,并用这一年的奇数日记录的ET、T和RH 来确定改进模型的最优参数(图4a-d)。在模型最优参数率定过程中,估算ET 与实测ET 间的R2为59%~61%。其中参数b在0.5~5 之间的变化对R2的影响甚小,取为常数1。R2对参数a和Tg的变化较为敏感,a的取值大小与表3 中的相关性大小一致,比如长寿花的ET 与气温只有较弱的相关性,其a= 0.5;而碧玉和蝴蝶兰的ET 受气温的影响较大,其a= 2.4。Tg的取值与花卉的最低适宜生长温度有关,如铜钱草和碧玉为15 ℃,长寿花为10 ℃,蝴蝶兰为20 ℃。T0的取值上下浮动2 ℃,对R2的影响都小于0.01。
图4 4 种花卉的模型估算蒸散量与实测蒸散量之间的散点图Figure 4 The scatter plots of the estimated ET against the measured ET for the 4 potted plants at calibration (a~d) and validation (e~h)
另用偶数日记录的ET、T和RH 来验证改进模型及其最优参数。根据日均气温和RH,改进的空气饱和差法能较准确地测算4 种花卉的日蒸散量(R2=0.52~57,RD = -0.9%~4.4%),这对快速测算4 种花卉的潜在蒸散耗水量和自动浇灌设计有较重要的推广应用价值,基于该蒸散发估算模型,可以根据较为容易获得的气温和相对湿度数据,快速估算和预测植物的蒸散耗水量,从而做好浇水规划和每日自动浇水量设计(图4e-h)。
本文采用称重法监测了水杯的水面蒸发量和铜钱草、长寿花、碧玉和蝴蝶兰的蒸散量,并根据实测数据改进空气饱和差法测算4 种花卉的日均蒸散量。结果表明:
(1)水杯的水面蒸发量能表征环境的最大蒸发潜力,6.4 cm 口径水杯的蒸发量与D20 蒸发皿以及开阔水面的自然蒸发量的折算系数约为0.63和0.44。
(2)在这365 d 的监测期间,水杯、铜钱草、长寿花、碧玉和蝴蝶兰单位花盆面积的日均蒸散量6.2、3.5、2.5、2.4 和1.8 mm。铜钱草蒸散耗水量最大,比开阔水面的自然蒸发量(2.7 mm/d)高28%,是蝴蝶兰的2 倍。4 种花卉的单位叶面积年蒸散量分别为162、29、74 和123 mm;铜钱草依旧最高,长寿花最低,蝴蝶兰排第二、是长寿花的4 倍和碧玉的1.7 倍。
(3)4 种花卉在白天的蒸散量都与RH 呈显著的负相关(95%置信度),而碧玉和蝴蝶兰在晚上的蒸散量与RH 的相关性不显著。长寿花的蒸散量与气温没有显著的相关性。风速会提高4 种花卉的蒸散量。
(4)当土壤水分供给受限时,或大气相对湿度过高(>80%)时,4 种花卉的蒸散量急剧下降并趋于雷同。
(5)根据日均气温和RH,改进的空气饱和差法能较准确地测算4 种花卉的日蒸散量,这对快速测算4 种花卉的潜在蒸散量和自动浇灌设计有较重要的推广应用价值。