云南昆明地区年季降水量空间演变特征分析

胡嘉佳，曹壹梅

(云南省水文水资源局昆明分局，云南 昆明 650106)

降水空间变化有助于区域旱涝特征分析，是旱涝分析的重要数据基础[1]。近些年来，对于区域降水时空变化特征的研究有许多[2- 8]，这其中既有数理统计学方法[9- 10]，也有趋势诊断分析方法[11- 12]，但是针对其空间变化特征的研究还较少，正交经验函数近些年来在变量空间变化特征应用较为成熟[13- 15]，该方法结合方差分析原理，对变量空间变化的特征进行空间解译，全方位诠释变量空间变化的特性。但是区域降水空间变化特性研究还较少，特别是在云南地区还未得到相关应用，昆明地区属于低纬度热带区，一年四季如春，降水量充沛，但有时也出现旱情，当前，在气候变化情形下，对区域旱涝变化特征的研究有利于昆明地区的水资源优化配置。为此本文结合正交经验函数，以昆明某区域为研究范例，对其年季降水时空变化特征进行分析，成果对于研究昆明地区旱涝演变具有参考价值。

1 正交经验函数原理

正交经验函数又称作为EOF函数，其降水空间变化分析原理在于将降水站点数据通过内插方式得到全域降水分布，并进行标准化处理，处理后形成一个标准的降水矢量变化矩阵，方程为：

(1)

式中，C—降雨矩阵数据，mm；M—内插数据点，个；N—观测雨量站点数目，个。

在矩阵确定的基础上，结合各个向量之间的方差值，计算方程为：

(2)

式中，R—变量方差；T—方差向量个数，个。

在方差计算基础上，计算各特征雨量点的特征根和特征向量，计算方程为：

RN×N×VN×N=VN×N×ΛN×N

(3)

(4)

式中，ΛN×N—特征向量根；λ—特征根。

每个特征向量对应一个特征根，这就是正交函数，对其空间变化特征分析函数为：

(5)

式中，a—正交经验函数的求解系数；f(z)—正交经验求解函数值。

结合最小二乘方法原理计算各变化矩阵对应的方程求解，计算方程为：

A=(VTV)-1V·ΔC

(6)

式中，A—各阶数对应的正交函数求解值；ΔC—计算点之间的特征分析值。

2 昆明地区降水量时空变化特征分析

2.1 降水概况

昆明地区属于典型的低纬度温带区域，多年平均降水量约为900mm，年季降水量具有明显的变化特征，从1970—2018年降水统计其年降水量呈现较为明显的递减变幅，除冬季外，其他季节降水也呈现减少变化，冬季降水呈现略增多变化，86%的降水主要集中在夏季和秋季，主汛期降水占全年降水的50%～60%之间，春季和冬季降水量偏少，占全年降水量的20%左右，且由于降水量的偏少，在春季出现干旱的频率达到75%以上。降水空间分布十分不均匀，东南部、南部降水量总体低于800mm，西部山区、东部降水量偏多，总体高于1000mm。

2.2 年降水正交经验函数空间分析结果

结合正交经验函数对研究区域年降水量进行空间分析，其特征方差向量及三个主要特征向量的分析结果见表1和图1所示。

表1 研究区域前10特征向量年降水对应的方差

从表1中可以看出，在第四个主成分下其特征值下降较为明显，因此前三个主成分作为其主要分析的特征向量，其三个主成分特征值分别为26.0、3.2、2.4，前三个主成分其特征向量累积贡献方差达到85%，所有选用前三个主要特征向量进行区域年降水量空间分析。从三个特征向量空间分析结果可看出，第一特征向量代表年降水空间变化的一致性，即总体呈现一致偏小，第二特征向量表示其年降水量东南和西北方向空间变化差异，第二特征向量分析结果可看出，特征向量变化正值出现呈现西北方向，表明区域年降水量空间上呈现东南少，西北多的变化特征。第三个特征向量主要表征为东部和西部降雨的变化情况，从其分布情况可看出，特征向量正值主要出现在西部，因此空间表现为东西相反，东少西多的空间变化，从其空间变幅可分析出，各区域年降水空间均有所递减，递减幅度在3～5mm/a之间，且东南部区域降水递减空间变化及变幅最大，东北部变幅最小。

2.3 春季降水正交经验函数空间分析结果

结合正交经验函数对研究区域春季降水量的空间变化进行分析，其特征向量方差变化及主要特征向量空间分布情况见表2及图2所示。

表2 研究区域前10特征向量春季降水对应的方差

图1 研究区域年降水量不同特征向量空间分布结果

图2 研究区域春季降水量不同特征向量空间分布结果

图3 研究区域夏季降水量不同特征向量空间分布结果

与年降水空间变化特征向量方差变化较为一致，前三个特征向量为主要成分，从第四个特征向量开始，其特征值变化较为明显，三个主要成分特征向量其特征值分别为34.5、5.2、2.5，其特征向量的累积方差在52.2%～63.9%之间变化。因此选用前三个主成分特征向量作为其春季降水空间变化的主要向量进行分析。从第一特征向量可以看出，正值向量较多，表明区域春季降水总体呈现减少的空间变化一致性，且东南减少的趋势更为明显。第二特征向量表征为东南和西北相同或相反的表征，从图2中可看出，负值向量变化较多，表明东南和西北两个方向降水量变化具有相反特征，第三特征向量主要表征为东部和西部春季降水量空间变化的表征，分析结果可看出，其负值出现在东部区域，表明东部春季降水总体减少，而西部春季降水减少趋势不明显，从空间降水变幅可分析出，春季降水空间变异度较小，东南区域递减幅度最大，达1.3mm/a，北部降幅为0.8mm/a。

2.4 夏季降水正交经验函数空间分析结果

结合正交经验函数对研究区域夏季降水量的空间变化进行分析，其特征向量方差变化及主要特征向量空间分布情况见表3及图3所示。

表3 研究区域前10特征向量夏季降水对应的方差

从表3中可看出，研究区域夏季降水量也同样表征为前三个特征向量，前三个特征向量的主要累积贡献方差可达到89%，这主要是因为夏季降水量属于区域降水量较为集中的季节，其降水量空间变化特征更为明显。从其三个特征向量空间变化可看出，其夏季降水呈现较为明显的变化一致性，总体呈现减少变化，其中中部降水减少量最明显，正值分布较多，也较为集中。而对于第二特征向量而言，其特征向量的正值均分布在东部区域，表明东部夏季降水变化具有一致性，总体减少，而西北部特征向量分布较少，表明西北部夏季降水量变化较小。第三特征向量主要表征为东部和西部降水空间变化，从第三特征向量分布结果可看出，其在西部的特征向量正值变化较多，西部夏季降水量变化具有一致性，总体减少变化。从其降水空间变幅分析可看出，东南部夏季降水减少幅度均值为2.8mm/a，东北部有所递增，增幅均值达0.7mm/a。

2.5 秋季降水正交经验函数空间分析结果

结合正交经验函数对研究区域秋季降水量的空间变化进行分析，其特征向量方差变化及主要特征向量空间分布情况见表4及图4所示。

表4 研究区域前10特征向量秋季降水对应的方差

秋季降水量总体也表征在前三个主成分的特征向量，前三个特征向量累积方差达到77.6%，从三个特征向量降水空间变化分析结果可看出，第一特征向量下其正值向量总体分布在西部区域，东部区域正值向量分布较少，秋季降水整体空间变化具有一致性。而第二个特征向量表示为东部正值向量较多，西部正值向量分布较少，呈现较为明显的东南和西北降水空间相反变化。从第三个特征向量分布情况可看出，其特征向量正值主要分布在西部，东部降水特征向量正值分布较少，表征出较为明显的东西相反变化。

2.6 冬季降水正交经验函数空间分析结果

结合正交经验函数对研究区域冬季降水量的空间变化进行分析，其特征向量方差变化及主要特征向量空间分布情况见表5及图5所示。

表5 研究区域前10特征向量冬季降水对应的方差

图4 研究区域秋季降水量不同特征向量空间分布结果

图5 研究区域冬季降水量不同特征向量空间分布结果

区域冬季降水量较少，但其空间变化的整体性和其他季节变化具有相似性，也出现在前三个特征向量上，前三个主成分累积方差的贡献率可以达到69.5%，同样选用三个主成分的特征向量进行其降水空间变化。第一特征向量主要表征其冬季降水空间变化的趋势，从第一特征向量分析结果可看出，其冬季降水空间变化具有较大的差异性，总体变化不具有一致性。第二特征向量主要表征其东南和西北区域降水空间变化，从分析结果可看出，其特征向量正值较少，表现为其东南和西北降水的冬季降水变化具有相反特征。从第三特征向量的分布情况而言，其特征向量分布更少，这主要是冬季降水较少原因所致。从变化结果可看出，区域冬季降水量总体呈现增加变化，西部增加趋势高于东部。从其降水空间变幅可看出，冬季降水有所增加，但增幅较小，东南部有所增加，增幅低于0.2mm/a。

3 结语

(1)各区域年降水空间均有所递减，递减幅度在3～5mm/a之间，且东南部区域降水递减空间变化及变幅最大，东北部变幅最小，应调控降水递减显著区域的水资源配置。

(2)春季降水空间变异度较小，东南区域递减幅度最大，达1.3mm/a，北部降幅为0.8mm/a，春季是区域干旱频发季节，应针对降雨递减明显区域进行抗旱措施规划。

(3)东南部夏季降水减少幅度均值为2.8mm/a，东北部有所递增，增幅均值达0.7mm/a，总体呈现空间相反型变化。

(4)冬季降水有所增加，但增幅较小，东南部有所增加，增幅低于0.2mm/a，东、西部呈现空间相反变化。