潜在蒸散是指在一定的气象条件下, 水分充分供给时某一固定区域下垫面可能达到的最大蒸散量, 它是水文循环以及能量平衡的重要环节^[1]。潜在蒸散量在作物需水量、气候条件以及水资源可持续发展等方面均得到了广泛应用^[2], 因此, 研究潜在蒸散量也对作物需水量的计算、水资源的合理配置提供了重要依据。

近年来, 国内外学者对潜在蒸散量方面做了大量研究。国外倾向于对潜在蒸散量的估算方法研究, Valiantzas^[3]介绍了计算潜在蒸散量的彭曼公式, 认为将本地平均风速作为附加输入时, 公式的结果精度会进一步提高；Raziei和Pereira等^[4]通过用彭曼公式对伊朗的潜在蒸散量进行估算, 认为在伊朗东部、南部干旱和极度干旱的气候条件下, 一些站点的性能较差。国内许多专家学者倾向于研究潜在蒸散量的时空变化特征, 曹永强和李维佳^[5]研究发现1966—2015年辽宁省的潜在蒸散量时空差异性显著, 并且在近50年整体上呈现递减趋势；钟巧等^[6]通过对博斯腾湖流域潜在蒸散量的计算发现, 该流域山区与平原年均潜在蒸散量都呈现出先下降后上升的趋势；刘昌明和张丹^[7]认为我国各区域普遍存在“蒸发悖论”, 全国范围内, 影响潜在蒸散量最敏感的气象因子为水汽压。

辽宁省地形由山地、丘陵、平原构成, 属温带大陆性季风气候, 雨量不均, 东湿西干, 是气候变化敏感地区之一, 潜在蒸散量变化及气象因子响应对于气候变化十分重要。研究潜在蒸散量变化及气象因子响应, 对辽宁省生态保护和作物生长具有重要作用, 但目前对潜在蒸散量与气象因子敏感性分析较多, 对辽宁省长时间尺度下潜在蒸散量时空变异规律以及相关气象因子的分析研究较少。因此, 本文在已有研究成果基础上, 对辽宁省21个气象站点1965—2017年的潜在蒸散量时空变异规律以及相关气象因子进行分析, 为该区农业水资源优化管理、作物合理布局及灌溉的高效利用建立提供一定参考依据。

1 研究资料与方法 1.1 研究区概况与数据来源

辽宁省位于中国东北地区的南部, 118°53′—125°46′E, 38°43′—43°26′N之间。南濒黄、渤二海, 雨热同期, 四季分明, 夏季炎热多雨, 冬季寒冷干燥, 多年平均气温7—11℃, 年均降水量750 mm, 年日照时数2543 h。

本文所需气象数据来源于中国气象数据网(http://data.cma.cn/), 选取辽宁省21个气象站点1965—2017年逐日最高气温、最低气温、相对湿度、风速、日照时数等气象数据, 对缺测数据采用线性插补法进行处理, 具体气象站点的分布如图 1所示。

1.2 研究方法 1.2.1 Penman-Monteith(P-M)法

大量研究表明联合国粮农组织(FAO)在1998年修正的Penman-monteith(P-M)法是目前计算潜在蒸散量ET₀最为精确和最常用的估算方法之一^[8-9], 其计算公式如下：

(1)

式中, ET₀代表站点潜在蒸散量(mm)；Δ为饱和水汽压与温度关系曲线的斜率值(kPa/℃)；R_n代表地面净辐射量(MJ/m²)；G代表土壤热通量(MJ m^-2 d^-1)；γ代表湿度计常数(kPa/℃)；T代表日平均温度(℃)；U₂代表 2 m高度处风速(m/s)；e_s代表饱和水汽压(kPa)；e_α代表实际水汽压(kPa/℃)^[10]。各因子具体计算公式详见参考文献[11]。

1.2.2 小波分析法

Morlet小波分析目前多应用于数学与工程领域, 是一种强有力的处理非稳定信号的统计工具。该方法能够同时从时域和频域揭示时间序列的局部特征^[12], 常用来分析气候序列中所包含的周期性及其不均匀性。

1.2.3 敏感系数

McCuen^[13]于1973年提出敏感系数法, 该方法基于偏导数的形式计算蒸散发变化量与气象因子变化量之比, 表达简洁操作简单。本文利用无量纲相对敏感系数, 分析辽宁省潜在蒸散量对各气候因子变化的敏感程度。具体计算公式如下：

(2)

式中, S_v为潜在蒸散量ET₀对某气象因子x的敏感系数, 其特点在于通过变化率将各气象因子无量纲化, 便于不同气象因子之间的比较^[14], 这种方法广泛应用于潜在蒸散量研究中。敏感系数的绝对值越小, 气象因子对ET₀的影响越小^[7]；敏感系数为正表明气象因子与ET₀变化趋势一致, 敏感系数为负则表明两者变化趋势相反。

1.2.4 因子分析

因子分析是指从研究指标内部关系出发, 用少数或几个潜在的具有代表性的公共因子去解释原始数据的一种多元统计分析方法^[15]。本文采用SPSS软件对辽宁省1965—2017年的21个气象站点影响潜在蒸散量的5个气象因子进行因子分析。

2 结果与分析 2.1 时空变化特征 2.1.1 辽宁省ET₀年际变化特征

辽宁省ET₀年际变化如图 2所示。由图可知, 1965—2017年辽宁省年平均ET₀呈显著下降趋势, 变化率为2.89 mm/a。多年平均值955.82 mm, 2010年出现ET₀最小值, 为820.90 mm, 1974年出现ET₀最大值, 为1149.16mm；整体上来看, 辽宁省ET₀年际变化较大, 呈现阶段性起伏变化, 反映了气象因子对潜在蒸散量的影响。在M-K突变检测中(图 3), 由潜在蒸散量UF曲线可知, 1965—2017年辽宁省ET₀呈下降趋势。根据UF和UB交点的位置, 确定UF和UB曲线在1978年出现突变点, 表明辽宁省潜在蒸散量在1978年发生下降趋势突然增大的现象。1965—1978年ET₀呈减少趋势, 但并未超过0.05显著水平临界线(-1.96), 表明1978年以前辽宁省潜在蒸散量减少趋势不明显。1978—2017年ET₀减少的趋势超过0.05显著水平临界线, 表明1978年以后辽宁省潜在蒸散量下降趋势十分显著。1975—1980年辽宁省降水量大幅减少^[16], 也可能是潜在蒸散量下降的原因。

对辽宁省21个气象站点潜在蒸散量的趋势变化进行分析, 并进行显著性检验。如图 4所示, 辽宁省21个气象站点中有20个气象站点潜在蒸散量呈递减趋势, 其中14个气象站点通过了0.01显著性检验, 3个气象站点通过了0.05显著性检验, 仅有3个气象站点未通过显著性检验(瓦房店站、营口站、彰武站)。仅有1个气象站点潜在蒸散量呈递增趋势(阜新县站), 且未通过显著性检验。整体来看, 辽宁省21个气象站点潜在蒸散量呈下降趋势变化。

2.1.2 周期变化分析

将各年份ET₀通过小波分析得出小波变换等值线及小波方差图(图 5)。从图 5可知, 小尺度变化镶嵌在较大尺度下的复杂环境当中。1965—2017年辽宁省ET₀在20—35年尺度上振荡周期最明显, ET₀在周期震荡中存在“大—小—大”的变化特征, 此外, 在13—15年尺度上也存在若干小周期循环。为了进一步明确ET₀随时间变化的主周期, 绘制小波方差分析图, 分析得出：辽宁省ET₀在35年左右尺度下的小波方差极值最大, 是主周期；其次在20年左右还存在一个峰值, 峰值数值较小, 为次周期, 即辽宁省潜在蒸散量在35年左右会经历由多变少的变化过程。

2.1.3 辽宁省ET₀空间变化特征

对ET₀多年平均数据进行空间分析(图 6)发现, 辽宁省年均潜在蒸散量地域差异十分显著, 表现出由西向东逐渐减小的趋势, 最高值(1043.58 mm)出现在辽宁省西部的朝阳站, 朝阳地处大陆性季风气候区, 北部为内蒙古高原, 经常性风速较大的干燥冷空气侵入, 是造成该地区潜在蒸散量大的重要原因。最低值(497.52 mm)位于辽宁东部的清原站, 由于所处地区气候特点不同导致各区域潜在蒸散量存在明显空间差异。辽宁西部地区位于辽宁辽河以西与内蒙古河北接壤处, 气候类型为半干旱气候^[17], 该区域日照丰富, 所以潜在蒸散量较大。辽宁东部地区位于辽河以东, 该区濒临黄海、渤海, 气候湿润, 降水量较多, 相对湿度偏大, 是潜在蒸散量较小的主要原因。

2.2 敏感性分析

为了研究辽宁省ET₀对气候变化的敏感程度, 综合考虑辽宁省各地的气候特点以及省内空间分布, 选取叶柏寿站、清原站、大连站、开原站、营口站、锦州站、鞍山站、丹东站、抚顺站以及本溪站10个气象站作为研究对象, 计算10个气象站潜在蒸散量对气候因子的敏感系数(表 1)。10个气象站的潜在蒸散量对相对湿度的敏感系数皆为负值, 且敏感系数的绝对值也最大, 大连更达到了1.21, 表明如果相对湿度减少10%, 则潜在蒸散量增加12.1%。潜在蒸散量与最低气温、最高气温、平均风速以及日照时数均呈正敏感, 这4个影响因子中, 潜在蒸散量对最高气温、平均风速及日照时数的敏感系数十分接近。潜在蒸散量对最低气温的敏感系数最小, 几乎为零, 敏感性最低, 主要原因是夜晚温度一般为一天中的最低气温, 而夜晚是一天中蒸散发最弱的时间, 这与李耀军等^[18]研究甘肃省潜在蒸散量的研究结果一致。从空间尺度上来看, 在经度和纬度改变的情况下, 计算发现ET₀对最低气温、最高气温、平均风速、日照时数的敏感系数基本保持不变, 而相对湿度在经度、纬度变化的情况下呈现出无规律的起伏变化趋势, 并且潜在蒸散量对相对湿度的敏感系数在各种情况下绝对值都为最大, 这说明潜在蒸散量对相对湿度的敏感性最高, 不因地区的不同而改变。

2.3 因子分析

本文利用SPSS统计分析软件, 得出所选取的5个气象因子的共同度取值都在0.5以上, 表明适合采用因子分析。由表 2可知, 成分1是方差贡献最大的一个主成分, 方差贡献率为53.999%。另外, 在5个变量因子中, 其起始特征值大于1的有2个成分, 这2个成分的方差贡献率分别为53.999%、27.658%, 累计方差贡献率为81.658%, 这说明了起始特征值大于1的有2个成份可以解释原来5个变量81.658%的信息, 表示分析时只需提取前2个公共因子即可。

通过因子分析法提取两个主成分进行主成分分析, 得到各主成分的荷载矩阵(表 3)。

第一个主成分的方差贡献率(53.999%)最大, 主要在最高气温上的荷载较大, 表明气温上升是潜在蒸散量升高的主要影响因子, 但潜在蒸散量不仅只随温度改变而变化, 其他气象因子的作用也影响潜在蒸散量的变化。其中在相对湿度上荷载为-0.244, 表明相对湿度增加, 潜在蒸散量反而减少。

第二个主成分的方差贡献率(27.658%)明显低于第一个主成分。第二个主成分主要在平均风速和日照时数两个指标上的荷载较大, 说明平均风速和日照时数的增加在某种程度上成为潜在蒸散量上升的因素, 相对湿度在第二主成分的荷载为-0.262, 说明相对湿度对潜在蒸散量的上升有制约作用。综上所述, 最高气温、日照时数和平均风速对区域潜在蒸散量影响显著, 相对湿度对区域潜在蒸散量呈负相关影响。不同因子对潜在蒸散量影响程度不同, 其中平均风速是潜在蒸散量的主要影响因子, 辽宁省地处温带大陆性季风气候区, 受天气系统及海陆热力差异影响, 平均风速较大, 尤其东南部沿海地区^[19]。但近35年中国地面风速呈递减趋势, 东北地区平均风速下降趋势最为明显^[20]。辽宁省处在东亚季风区, 有研究表明, 温室气体的排放及气溶胶的变化, 导致东亚夏季风和冬季风正处在减弱趋势时期, 季风减弱则导致平均地面风速下降^[21], 该变化是造成辽宁省潜在蒸散量呈下降趋势的重要因素。由于人类破坏环境愈发严重, 太阳辐射量日渐减少, 日照时数也呈减少趋势, 导致地面接收的辐射能量减少, 潜在蒸散量下降。据研究表明, 辽宁省近些年及未来40年气候呈暖干趋势^[22], 相对湿度下降, 促进潜在蒸散量上升, 但由于最低气温、最高气温、日照时数、平均风速4个气象因子与潜在蒸散量呈正相关, 与相对湿度的负相关抵消, 从而最终导致潜在蒸散量整体呈下降趋势。

3 结论

(1) 从时空变异规律来看, 辽宁省1965—2017年潜在蒸散量多年平均值为955.82 mm, 变化速率为2.89 mm/a, 下降趋势显著。在空间上不同地区的潜在蒸散量差异明显, 表现出由西向东逐渐减少的趋势。

(2) 辽宁省潜在蒸散量对相对湿度的敏感性为负, 且最为敏感；对最低气温、最高气温、风速、日照时数均呈正敏感, 且对最低气温敏感性最小, 几乎为零。在经度和纬度改变的情况下, 最高气温、风速、最低气温、日照时数的敏感性几乎没有改变, 而相对湿度的敏感性呈现无规律的阶段性起伏变化。

(3) 平均风速是潜在蒸散量变化的主要影响因子, 相对湿度与潜在蒸散量呈负相关, 表明相对湿度增加, 则潜在蒸散量呈减少趋势。几个气象因子综合作用导致辽宁省潜在蒸散量呈显著下降趋势。

4 讨论

本文基于小波分析、敏感性分析和因子分析, 定量分析与定性分析相结合, 深入研究了辽宁省潜在蒸散量的时空变化特征及成因。结果表明, 辽宁省近53年年均潜在蒸散量呈下降趋势, 平均风速与日照时数的减少是辽宁省潜在蒸散量下降的主要原因, 气温上升并未导致潜在蒸散量上升, 表明辽宁省存在“蒸发悖论”的现象, 这与丛振涛等^[23]东北地区的研究结果一致。但是本文仅从潜在蒸散量传统公式方面分析了它的时空变化特征, 并未通过“蒸发悖论”成因与规律探究辽宁省潜在蒸散量的变化趋势, 结合下垫面要素、农业活动等人为因素对潜在蒸散量的影响, 潜在蒸散量的变化趋势特征还有待进一步探讨。另外, 影响因子分析中发现平均风速下降是辽宁省潜在蒸散量下降的主要影响因子, 但本文仅仅进行了初步探究, 其他气象要素以及植被、人为因素、地形等对潜在蒸散量的影响需要进一步研究。