陆地蒸散发分割方法及影响因素研究进展
2021-07-06李新梦狄崇利王铁军
李新梦,狄崇利,2,刘 琴,韩 琼,王铁军,2
(1.天津大学 地球系统科学学院,表层地球系统科学研究院,天津 300072;2.天津市环渤海关键带科学与可持续发展重点实验室,天津 300072)
1 研究背景
陆地蒸散发(ET)是全球水文循环的重要组成部分,在全球尺度上约有73%的陆地降水通过蒸散发返回到大气中[1],在干旱地区甚至可以达到90%以上[2]。陆地蒸散发是土壤-植物-大气连续体系(soil-plant-atmosphere continuum, SPAC)能量传输的重要过程,蒸散发过程中水分汽化吸收的热量约为地表吸收太阳净辐射量的55%,同时陆地蒸散发过程也是联系生物圈、水圈和大气圈物质交换的主要环节[3]。陆地蒸散发研究对深入了解水循环过程以及准确评估气候变化对自然生态系统的影响等具有重要科学意义[4]。
陆地蒸散发由蒸腾(T)和蒸发(E)两部分组成。其中,蒸腾由植物蒸腾作用实现,而蒸发则主要包括土面和冠层蒸发。植物通过蒸腾作用向大气输送水分的同时从土壤中吸水,从而保持植物体内水分和温度的平衡,因此进行蒸散发分割对理解SPAC中的水分运动、揭示生态系统的生态水文过程具有重要作用[5]。此外,全球农业灌溉水的平均利用率仅为45%左右[6],其主要原因之一是对作物生育期蒸散发及各组分动态变化缺乏深入理解,无法准确估算蒸散发及各组分量。因此,研究陆地蒸散发分割对提高水资源利用效率[7],实现水资源优化管理具有重要的现实意义[8]。
蒸散发分割是围绕生态系统中蒸散发各组分组成及其动态变化而进行的研究。国内外学者就蒸散发分割提出了一系列的理论和方法,这些方法在不同生态系统中得到了广泛应用[9-10]。由于蒸腾和蒸发的影响因素存在差异[11-12],所以相较于蒸散发,蒸散发分割的影响因素更为复杂,现有研究对蒸散发分割影响因素的总结较少,且相关研究多集中于特定研究区[13-14]。本文在总结蒸散发分割主要方法的基础上,在不同时间尺度上对蒸散发分割的影响因素进行了系统地探讨,以期为蒸散发分割的进一步研究提供参考。
2 蒸散发分割方法
蒸散发分割方法可分为3类[15],即水文气象学方法[16-19]、稳定同位素方法[20-22]及数值模拟方法[23-25]。
2.1 水文气象学方法
水文气象学方法是利用涡度相关系统、波文比系统、蒸渗仪以及植物茎流计等仪器中的两种或两种以上相互结合,通过分别测量蒸散发量、蒸腾量和蒸发量中的两种或3种实现蒸散发分割。水文气象学方法在不同生态系统中均有应用[16-19]。祁秀娇等[16]利用植物茎流计和涡度相关系统在内蒙古小叶锦鸡儿灌丛中分别测量了灌丛的植物蒸腾量和蒸散发量,发现该灌丛中T/ET(蒸腾与蒸散发的比值)约为66%。Moran等[19]利用波文比系统和微型蒸渗仪分别测量了灌丛和草地的蒸散发量和蒸发量,结果表明灌丛和草地的T/ET分别为64%和79%。王渝淞等[17]以森林为研究对象,利用大型蒸渗仪和植物茎流计发现该林分(即内部特征大体一致,而与临近地段有明显区别的一片森林)蒸散发主要由蒸腾组成(约占89.6%)。
利用水文气象学方法可直接测量蒸散发及其各组分,物理过程明确,仪器测量精度高,所得蒸散发分割的结果准确,通常被作为其他分割方法的验证标准[23]。汉光昭等[23]利用基于涡度相关系统的水文气象学方法对Shuttleworth-Wallace(S-W)模型得到的结果进行检验,验证了S-W模型的准确性。但是,由于水文气象学方法受实验场地、观测时间、仪器测量精度以及仪器费用等因素限制,同时对下垫面稳定性和均一性要求高,因此通常用于小尺度的蒸散发分割研究[15]。
2.2 同位素方法
利用稳定同位素方法进行蒸散发分割,是根据蒸腾和蒸发产生水汽的氢氧稳定同位素组成差别实现的,但该方法对某些盐生、旱生植物并不适用[26]。目前,氢氧稳定同位素技术在灌木[27]、草地[27]、农田[28]以及森林[17]等生态系统中均有广泛应用。Wang等[27]在内蒙古灌丛-草原生态系统中利用同位素示踪方法量化了灌丛和草地的土壤水分利用情况,结果发现生长季灌丛蒸腾、草地蒸腾和土壤蒸发对总蒸散发量的贡献分别为24%、20%和56%。与大气水汽冷阱同位素质谱仪技术相比,大气水汽稳定同位素组成的原位连续监测技术已逐渐成为蒸散发分割研究中同位素方法研究的热点。娄源海[29]利用原位连续监测技术对侧柏林进行研究,结果表明在研究期间T/ET的平均值约为78%。
利用稳定同位素方法进行蒸散发分割的最大优势在于稳定同位素贯穿了生态系统水分循环的全过程[30],且自身的稳定性不易受外界因素干扰,数据测定精度高。但稳定同位素方法同样存在费用较高、仪器维修不便等问题,而且还会受到样品收集效率以及分析仪器精度的双重制约[30]。
2.3 数值模拟方法
数值模拟由于其计算简单、易操作的优点,在许多生态系统的蒸散发分割研究中均有应用,其中应用最为广泛的是S-W模型[31]和双作物系数模型[32]。Jiang等[33]利用S-W模型对黄河流域1981-2010年的蒸散发进行了划分,发现年均蒸散发量和蒸腾量分别为372.18和179.64 mm,T/ET约为48.3%。Bao等[34]利用S-W模型对科尔沁沙地典型地貌-土壤-植被组合(流动沙丘、半固定沙丘和固定沙丘)的蒸散发进行了模拟和划分,并用涡度相关系统对模型进行了验证,结果发现研究期内3种不同的地貌-土壤-植被组合E/ET(土壤蒸发与蒸散发的比值)分别为57.1%、25.2%和21.3%,随着植被覆盖度的增加,土壤蒸发所占的比重逐渐降低。于文颖等[35]利用FAO-56双作物系数模型、改进后的双作物系数模型以及基于涡度相关系统的水文气象学方法进行研究,发现利用这3种方法得到的趋势基本一样,研究期间蒸腾一直是蒸散发的主要组成部分。
相较于实测方法,数值模拟方法更为方便简单,并且利用数值模拟方法可以实现更大时空尺度上的蒸散发分割研究,Fatichi等[36]利用机理性生态水文模型发现在全球尺度上T/ET约为70%。但数值模拟也有很大的局限性,如S-W模型适用的时间尺度较短,通常少于一个生长期,当涉及更长的时间尺度时,模型的准确性降低[31];其次,模型参数的选择误差会造成模拟结果的准确度不高、模型普适性差,同时目前部分参数的获取也存在问题[9]。
综合分析比较上述3种蒸散发分割方法的优缺点(见表1),水文气象学方法和稳定同位素方法物理基础坚实且计算结果准确[37],但由于场地和观测时间限制,大多集中在田间尺度及植物生长季[9],同时大多数的测量仪器存在价格昂贵、维修不便等缺点;数值模拟方法操作简单且可实现区域尺度上的研究,但是由于部分模型参数地域性强、数量大,导致模型的求解难度大、可靠性低、适用性差,且模型在与遥感方法结合的过程中存在着参数不匹配等问题[9]。
表1 陆地蒸散发分割方法比较
3 蒸散发分割的环境影响因素
蒸散发及其组分的影响因素主要有气象、土壤和植被等因素[42-44]。由于环境因子对蒸发和蒸腾的控制不同(例如,蒸发和蒸腾受土壤含水量、土壤温度、风速等影响,此外蒸腾还受植被生理状况、固碳能力等因素影响[45]),所以同一因素对蒸散发和蒸散发分割的影响不同。
在不同时间尺度上,影响蒸散发分割的主导因素存在差异,以下拟从日尺度以及季节尺度两方面对蒸散发分割的环境影响因素(分别包括太阳辐射、土壤温度、植被类型及降水、土壤含水量、植被生长状况等)进行论述。
3.1 日尺度上蒸散发分割的环境影响因素
在日尺度上,对蒸散发分割影响最为明显的是气象因素,S-W模型[31]、TSEB(two source energy balance)模型[46]等均说明了气象因素对蒸散发分割的重要性;植被类型也是影响蒸散发分割的重要因素[47]。在日尺度上太阳辐射变化明显,是日尺度上影响蒸散发分割的重要环境因素,由此导致的温度、相对湿度变化对日尺度上蒸散发分割结果影响显著。同时,由植被类型和各类环境因素引起的气孔导度差异通过影响蒸腾作用进而影响蒸散发分割的结果。相对而言,植被生长状况、土壤含水量等在无水分输入情况下日变化不大,故在日尺度上对蒸散发分割的影响不明显。
3.1.1 太阳辐射和温度对蒸散发分割的影响 太阳辐射是影响蒸散发分割的重要因素[48],日尺度上太阳辐射与温度之间存在明显的正相关关系。在上午时分,太阳辐射增强,温度升高,蒸腾作用明显加强,T/ET呈现上升趋势;在正午时分,由于高温以及强辐射的影响,T/ET达到峰值后呈现下降趋势。在草地生态系统中利用氘同位素发现,日尺度上的T/ET随着太阳辐射的增强先呈现上升趋势,在正午时分达到最大值,随后由于太阳辐射增强导致植物叶片和土壤温度升高,植物气孔关闭,蒸腾作用减弱,蒸发增强,T/ET呈现下降趋势[49]。Wang等[48]对日本温带草原的研究以及孟飞[50]在疏勒河地区的蒸散发分割研究均发现在日尺度上T/ET与太阳辐射之间呈现正相关关系。在日尺度上,T/ET与温度之间也呈现正相关关系。Wang等[51]分别在与周围环境状况相同和额外提供热源两种不同的情况下进行实验,发现在加热条件下T/ET值明显增加。由此可见,在日尺度上,T/ET与太阳辐射以及温度之间均呈现正相关关系,但在正午时,由于太阳辐射过强,温度升高,植物水分流失严重,出现下降现象。
3.1.2 植被类型对蒸散发分割的影响 植被类型也是日尺度上影响蒸散发分割的重要因素,这主要体现在植被气孔对蒸腾作用的影响上[52],其影响程度一方面取决于植物表皮气孔的分布情况[53],另一方面则取决于气孔的张开程度[54],即气孔导度。其中植物表皮气孔的分布状况主要由植被类型决定[53],气孔导度则主要是植被与环境因素共同作用的结果[55-56]。关于气孔导度,一般而言,阔叶植物大于针叶植物,C3植物大于C4植物,热带植物大于温带植物,而在同一类甚至同一株植物中,新叶大于旧叶[57]。所以,日尺度上不同植被类型对蒸散发分割的影响体现在不同植被气孔导度的差异上,而同一植被类型蒸散发分割结果的差异主要受环境驱动的气孔导度影响。
Zhu等[47]分别在温带阔叶红松混交林(CBS)、亚热带针叶林(QYZ)和亚热带常绿混交林(DHS)中进行了蒸散发分割研究,结果发现3种森林中,气孔导度是日尺度上影响蒸散发分割的主要因素,CBS的气孔导度最大,相应的T/ET最高,QYZ的气孔导度最小,T/ET也最低,且在同一种植被类型中,T/ET随着气孔导度的增加呈现增大趋势。董军等[58]在玉米地进行蒸散发研究时,同样发现气孔导度是日尺度上影响蒸散发分割的因素,并认为这是由于太阳辐射的变化引起的,是环境因素变化导致的结果。由此可见,在日尺度上,植被类型差异和环境因素影响造成的气孔导度差异是影响蒸散发分割的重要因素,T/ET随气孔导度的增大而增大。
3.1.3 其他环境因素对蒸散发分割的影响 除了太阳辐射、温度和植被类型,蒸汽压差[59]、相对湿度[49]等同样会影响蒸散发分割的结果。Hu等[59]通过研究发现T/ET的日变化与蒸汽压差有明显的相关性,早晨蒸汽压差较低,蒸腾作用变化更强烈,而在下午,蒸汽压差增加,蒸发量的增加大于蒸腾量的增加,从而使得T/ET降低。在对肯尼亚中部某草原的研究过程中发现相对湿度也是日尺度上影响蒸散发分割的重要因素,相对湿度与T/ET之间呈现正相关关系,上午10:00之前,相对湿度逐渐增大,T/ET呈现上升趋势,下午(15:00-16:00)相对湿度小幅度上升的同时T/ET也出现小幅度上升[49]。
综上所述,在日尺度上,太阳辐射是影响蒸散发分割的最主要因素,伴随太阳辐射变化造成的温度、蒸汽压差、相对湿度的变化以及与植被相互作用导致的气孔导度的变化也是日尺度上影响蒸散发分割的重要因素。目前,日尺度上影响因素的研究大多侧重于太阳辐射以及植被两方面的研究,对其他因素的研究相对较少,且研究过程中通常研究单一气象因素对蒸散发分割的影响,缺少各气象和植被因素之间相互作用对蒸散发分割产生的影响。
3.2 季节尺度上蒸散发分割的环境影响因素
在较长的时间尺度上,蒸散发分割的影响因素更为多样。在季节尺度上,气象因素[60]、土壤因素[61]以及植被因素[62]对蒸散发分割均有明显影响。气象因素中降水的变化对蒸散发分割的影响最为突出,由降水导致的土壤含水量的变化对蒸散发分割的影响也较为强烈[61]。此外,在季节尺度上,植物的物候变化明显,植被对蒸散发分割的影响是不可忽略的,其中植被的生长状况,如叶面积指数(LAI)、归一化植被指数(NDVI)以及植被覆盖度[49]对蒸散发分割的影响尤为强烈。
3.2.1 降水对蒸散发分割的影响 水分是影响蒸散发各组分的重要因素,在干旱半干旱地区尤为明显。降水通过增加地表的湿润程度,进而影响蒸散发分割,使得T/ET降低。Jiao等[39]在中国黄土高原地区进行的蒸散发分割研究表明,在季节尺度上,蒸散发分割的结果与降水量关系密切,随着降水量的不断增加,T/ET呈现不断减小的趋势,两者之间呈对数关系。沈竞[63]通过实验发现,湿润区植物冠层截流作用产生的蒸发(Ec)是不能忽略的,Ec/ET与降水量之间呈正相关关系,与暴雨频次之间呈负相关关系。由此可见,不仅降水量会影响蒸散发分割的结果,降水频次或降水天数同样会影响蒸散发分割的结果。Scott等[64]通过对美国4个干旱半干旱地区的研究发现,地表湿润程度随降水天数的增加而增强,且土壤蒸发增量大于植物蒸腾,T/ET与降水天数之间存在负相关关系。但另有研究表明[65],在生长季尺度上,降水不是蒸散发分割的主要影响因素。根据以往研究结果,得到季节尺度上降水量与T/ET之间的关系如图1所示,分析其Pearson系数发现,降水量与T/ET之间的显著性较强,呈正相关关系(R=0.477,P<0.0001)。该结果与已有研究结果存在一定偏差[43],可能是由于数据中生态系统类型与植被类型多样,且农田生态系统中还存在灌溉、施肥以及覆盖地膜等人为活动的影响,增加了季节尺度上的水分输入,提高了水分利用效率,植物生长更加旺盛,T/ET在季节尺度上增加较为明显。
图1 季节尺度上降水量与T/ET的关系
3.2.2 土壤含水量对蒸散发分割的影响 在季节尺度上,土壤含水量的变化也是影响蒸散发分割的重要因素,表层土壤含水量越高,则T/ET越小。蒸发和蒸腾利用的分别是表层土壤水和较深层的土壤水,所以表层土壤水的增加会造成T/ET减小。Rigden等[61]对美国大陆进行研究发现,在土壤含水量降低过程中,T/ET逐渐减小,尤其是美国西部缺水地区,产生这种现象主要是由于表层土壤水用于蒸发作用,而较深层的土壤水则主要是由植物根部吸收并供给植物蒸腾[66],所以在表层土壤含水量减少时蒸发量的减小速度大于蒸腾量。Bao等[34]对科尔沁沙地典型地貌-土壤-植被组合的蒸散发进行了模拟和划分,并研究了土壤水分与蒸散发分割之间的关系,结果发现土壤水含量是影响蒸散发分割的重要因素,随着土壤水含量的增加土壤蒸发所占的比重逐渐升高。在农田生态系统中,人类活动(灌溉、耕作方式等)都会影响表层土壤含水量[67-68]。相较于滴灌,漫灌方式增加了表层土壤含水量,T/ET值更小[67]。耕作方式也影响表层土壤含水量进而影响蒸散发分割,在王月月[68]的研究中发现,与翻耕相比,免耕由于凋落物的影响,虽然增加了土壤水含量,但由于覆盖作用,使蒸发量减少,T/ET降低,其原理与地膜覆盖类似[69]。
3.2.3 植被生长状况对蒸散发分割的影响 植被在生长过程中会经历明显的季节变化,LAI[63]、NDVI[67]以及植被覆盖度[70]等都能够反映植被的生长状况,并影响蒸散发分割的结果。通常情况下LAI(或NDVI、植被覆盖度)的增加会导致T/ET的增大,这是由于LAI增大使植物蒸腾和叶片遮荫效果增强,地表温度降低,土壤蒸发减弱,最终导致T/ET增大。沈竞[63]通过对中亚热带人工针叶林的蒸散发及其组分拆分研究发现,在该生态系统中,T/ET约为85%,与LAI之间呈现正相关关系。此结果在草原生态系统[48]和灌丛生态系统[34]中有相似结论。Jiang等[33]通过研究发现NDVI是控制黄河流域植物蒸腾和蒸散发的重要因素,T/ET随着NDVI的增加呈上升趋势。Wang等[70]利用同位素方法发现植被覆盖度与T/ET同样存在上述关系。综上所述,植物生长状况是季节尺度上影响蒸散发分割的重要因素,LAI、NDVI以及植被覆盖度的增加都会导致T/ET的增大。
综上所述,在季节尺度上,蒸散发分割主要与植物的物候特征有关,而植物的物候特征又与降水、气温、土壤等因素有密切联系。由此可见,在季节尺度上,LAI、NDVI、植被覆盖度是影响蒸散发分割最重要的因素,降水、气温以及土壤等环境因素通过影响植物生长和蒸发耗水也对蒸散发分割造成影响。
不同时间尺度上T/ET的影响因素如图2所示。由图2可看出,在日尺度上,影响蒸散发分割的因素主要有太阳辐射、温度、相对湿度、蒸汽压差以及由植被类型和环境因素共同控制的气孔导度,其中与T/ET呈负相关关系的因素是蒸汽压差,其余因素均呈现正相关关系;在季节尺度上,蒸散发分割的控制因素主要有降水(降水量和降水频次)、土壤含水量以及植被生长状况(可用LAI、NDVI等指数表示),其中与T/ET呈正相关关系的因素是植被生长状况(LAI、NDVI等),其余因素则与T/ET呈负相关关系。
不同影响因素在蒸散发分割领域的研究已经较为普遍,但由于不同影响因素之间存在强烈的相互作用(见图2),如植被生长过程中通常伴随着降水以及温度的变化,所以在研究过程中通常忽略了各影响因素之间的相互作用。
注:⊕表示两因素间存在正相关关系;⊖表示两因素间存在负相关关系。
4 目前研究存在的问题与展望
4.1 目前研究存在的问题
随着研究的深入,关于蒸散发分割的方法已经越来越完善,在此基础上,不同时间尺度上蒸散发分割的影响因素也逐步清晰,但仍存在一些亟待解决的问题。
蒸散发分割的方法目前已经取得了长足的发展,但在一些方面仍需要进一步完善,这主要体现在:(1)多数实测方法(涡度相关法、波文比法等)以及数值模拟方法(双作物系数模型、互补相关模型等)要求下垫面均匀稳定,但是在实际情况下,这种要求通常难以满足,导致蒸散发分割的结果不够精确;(2)各种实测方法都有其物理基础及适用范围,例如,涡度相关系统基于雷诺分解可连续测得地表蒸散量,而蒸渗仪则是基于质量守恒以小时或天为单位进行测量,所以直接将各种实测方法进行结合、比较或验证仍存在问题;(3)目前已有的蒸散发分割模型在与遥感方法结合的过程中仍存在较大的误差,同时部分模型参数的确定需要借助经验公式,参数依赖于研究区,不具有普适性。
蒸散发分割影响因素的研究同样存在以下问题:(1)在探讨各因素对蒸散发分割的影响时,通常对单因素进行研究,忽略了各影响因子之间相互作用对蒸散发分割的影响,如太阳辐射与土壤温度之间的相互关系、植物的生长状况与各环境因素之间的关系等;(2)各因素对蒸散发分割的影响通常不止体现在一种时间尺度上,但目前研究通常集中于特定的时间尺度,如植被的生长状况既可以在日尺度上通过气孔导度的差异进行蒸散发分割,也可以在季节尺度上影响蒸散发分割的结果。
4.2 研究展望
蒸散发分割的实测方法是获取地表蒸散发分割数据的基本途径,也是检验模型可行性与准确性的有效手段。实测手段未来需要从两方面进行完善:(1)改进仪器,提高仪器的精确性以及适用性;(2)稳定同位素方法在蒸散发分割研究中的应用逐渐广泛,但由于设备价格昂贵、维修不便等问题,应寻找一种成本相对较低的同位素方法进行推广。
数值模拟可以向以下3个方向发展:(1)提出同时考虑时空延展性以及下垫面均匀性的动态模型;(2)优化模型模拟与遥感数据的结合以提高遥感数据在蒸散发分割研究中的应用;(3)采取多种手段进行研究,如目前已有的结合植物生物量[71]以及水分利用效率[14]等进行研究,以期能更方便地得到更准确的结果。
关于蒸散发分割的影响因素,未来可以从以下3个方面进行研究:(1)考虑各影响因素之间的相互作用,如人为活动、气候变化以及两者间的相互作用对蒸散发分割的影响;(2)对同一影响因素在不同时间尺度上影响蒸散发分割的机理研究是未来蒸散发分割影响因素研究的重点;(3)不同生态系统、不同气候区蒸散发分割对环境因子的响应也是本领域研究的重要内容。