蒸发是海洋和陆地水分进入大气的唯一途径,是地球物质、能量循环的主要环节之一,自然界中蒸发现象颇为复杂[1],水源充足的自然水面,水源和溶质浓度两个因素可以忽略不计,蒸发仅与气象因子有关。蒸发皿蒸发量虽然不能直接代表水面蒸发,但它与水面蒸发之间存在很好的折算系数,因此蒸发皿蒸发量的变化趋势可以代表水面蒸发的变化趋势[2]。东疆地区干旱少雨,同时也是大风、高温天气的高发地,蒸发旺盛,吐鲁番地区夏季蒸发量居全国之首,研究该地区的蒸发量变化趋势及其影响因子,对于东疆地区水资源评价、水利工程设计、气候区划、生态研究等具有重要的意义。
1 方法选取东疆地区吐鲁番、哈密、巴里坤3个气象观测站1961—2001年每日的蒸发量、日照时数及气温、相对湿度、风速、水汽压的日平均值等地面气象观测资料,其中巴里坤位于天山北部,吐鲁番和哈密位于天山南部。蒸发量用口径为20cm的小型蒸发皿测得,地表净辐射通量用经验公式获得,具体计算方法参见文献[3]。
1.1 趋势估测及突变检验蒸发量的变化趋势分析采用线性倾向估计和11年滑动平均两种方法,变化趋势的显著性检验和突变检验采用曼-肯德尔(Mann-Kendall)法,具体参见文献[4]。
1.2 通径分析通径分析(Path Analysis)是研究变量间相互关系、自变量对因变量作用方式、程度的多元统计分析技术。通径系数是变量标准化后的偏回归系数,它能有效地表示相关变量间单一原因因素对结果的直接影响效应,并兼顾两原因因素对结果的影响,估计出原因对响应因素的间接作用,从而直接比较各原因因素的相对重要性。
通径分析算法概要[5]:设有x1, x2, …,xn n个自变量,这些变量之间单相关系数为γij,各自变量对因变量y的直接通径系数为Pyxi,系数相互间组成正规方程组:
$ \left\{ \begin{array}{l} {P_{y{x_1}}} + {\gamma _{12}}{P_{y{x_2}}} + \cdots + {\gamma _{1n}}{P_{y{x_n}}} = {\gamma _{1y}}\\ {\gamma _{21}}{P_{y{x_1}}} + {P_{y{x_2}}} + \cdots + {\gamma _{2n}}{P_{y{x_n}}} = {\gamma _{2y}}\\ \cdots \;\;\;\; \cdots \;\;\;\; \cdots \;\;\;\; \cdots \;\;\;\; \cdots \;\;\;\; \cdots \\ {\gamma _{n1}}{P_{y{x_1}}} + {\gamma _{n2}}{P_{y{x_2}}} + \cdots + {P_{y{x_n}}} = {\gamma _{ny}} \end{array} \right. $ | (1) |
其矩阵式为:
$ \left[\begin{array}{l} 1\;\;\;{\gamma _{12}}\;\;\;{\gamma _{13}}\;\;\; \cdots \;\;\;{\gamma _{1n}}\\ {\gamma _{21}}\;\;1\;\;\;{\gamma _{23}}\;\;\; \cdots \;\;\;{\gamma _{2n}}\\ \cdots \;\;\; \cdots \;\;\; \cdots \;\;\; \cdots \;\;\; \cdots \\ {\gamma _{n1}}\;{\gamma _{n2}}\;\; \cdots \;\;{\gamma _{n\left( {n-1} \right)}}\;\;1 \end{array} \right]\left[\begin{array}{l} {P_{y{x_1}}}\\ {P_{y{x_2}}}\\ \cdots \\ {P_{y{x_n}}} \end{array} \right] = \left[\begin{array}{l} {\gamma _{1y}}\\ {\gamma _{2y}}\\ \cdots \\ {\gamma _{ny}} \end{array} \right] $ | (2) |
解式(2),可得到各直接通径系数Pyxi。xi变量通过xj变量对y的间接通径系数为:γijPyxj。xj变量通过xi变量对y的间接通径系数为:γijPyxi。某原因对结果的直接通径系数与间接通径系数总和恰好等于该原因和结果的相关系数:γiy=Pyxi+∑njPyxj。
各自变量对因变量的相对决定程度可用决定系数表示:
$ {D_y}{x_i} = P_{y{x_i}}^2 $ | (3) |
$ {D_y}{x_i}{x_j} = 2{\gamma _{ij}}{P_{y{x_i}}}{P_{y{x_j}}} $ | (4) |
$ {D_y}{x_i} + {D_y}{x_i}{x_j} + {D_y}E = 1 $ | (5) |
式(3)~(5)中Dyxi为原因xi对结果y的决定系数,Dyxixj为原因xi、xj协同作用对结果y的决定系数,DyE为剩余项的决定系数。
2 结果 2.1 蒸发量变化趋势及突变检验对东疆3站平均日蒸发量通过线性倾向估计和11年滑动平均分析得出(见图 1),东疆地区蒸发量呈逐年下降趋势。其线性倾向率为-0.28mm/10a,经检验,此变化趋势在0.05显著水平下达到显著。
对东疆3站1961—2001年平均日蒸发量用曼-肯德尔进行突变检验(见图 2),结果表明,东疆地区蒸发量共经历2次增多阶段和2次减少阶段,其中1960年代初至1960年代中期为增多趋势,1960年代后期为减少趋势,1970年代至1980年代初期为增多趋势,1980年代中期至21世纪初为减少趋势。其中1960年代初至1960年代中期和1970年代至1980年代初期的增多趋势及1980年代中期至21世纪初的减少趋势均超过显著水平0.05临界线,表明这几个阶段蒸发量增多或减少明显。东疆地区蒸发量1970年代至1980年代初期的增多和1980年代中期至21世纪初蒸发量的减少是突变现象,分别是从1972年和1990年开始的。
影响蒸发量的各个气象因子内部存在着相关性,为了清楚地了解各种气象因子在蒸发过程中所起的作用,采用通径分析进行分析。参考前人研究成果[1, 6-7]引入气温(x1)、相对湿度(x2)、风速(x3)、日照时数(x4)、水汽压(x5)、净地表辐射(x6)6种与蒸发相关的气象因子。
从表 1可以看出,6种气象因子与蒸发量的直接相关性均达到0.01极显著水平,除相对湿度与蒸发量呈负相关外,其他5种气象因子均与蒸发量呈正相关。表 2为计算出的通径系数,表示的是各种气象因子对蒸发量的直接影响及间接影响作用,通过对比表 1、表 2可以发现,相对湿度、风速、日照时数、水汽压的直接通径系数与其对应的直接相关系数正负符号相同,而气温、净地表辐射的直接通径系数与其对应的直接相关系数正负符号相反,这表明相对湿度、风速、日照时数、水汽压对蒸发量的影响作用是真实可信的,而气温、净地表辐射与蒸发量的显著相关可能存在着一定的虚假性。表 2中, 相对湿度与水汽压对蒸发量的直接影响作用符号相反, 这是由于水汽压表示的是大气中水汽压的分压强, 在大气压变化不定的自然条件下, 水汽压大并不一定抑制水汽的蒸发作用, 而相对湿度表示的是大气中水汽的饱和程度, 它的大小与蒸发量呈反比。
表 3为各气象因子对蒸发量的决定系数, 结果表明,6个气象因子单独对蒸发的决定程度按大小排序为:相对湿度、水汽压、气温、风速、日照时数、净地表辐射,分别位于第3、第4、第7、第9、第15、第21位,风速与日照时数、相对湿度与水汽压、气温与水汽压、气温与相对湿度的协同作用对蒸发量的决定作用很大,分别位于第1、第2、第5位和第6位。净地表辐射无论其单独作用或与其它气象因子的协同作用对蒸发量的决定作用都很小,排名全部在第13位之后。结合表 2的分析结果, 可以得出,气温对蒸发量的影响作用是显著的,排除其显著相关虚假的可能性,而净地表辐射对蒸发量的显著相关是虚假的。
根据前人研究[8-10],东疆近40年气候呈现降水显著增多、气温显著增暖的趋势,蒸发量呈显著减少趋势。依据气象因子与蒸发量的相关性可以看出,气温与蒸发量呈正相关,因而气温偏高对蒸发量减少不利,降水增多会增加空气的相对湿度,而相对湿度与蒸发量呈负相关,因而降水增多对蒸发量减少有利。在目前气候的变化趋势对蒸发量减少有利有弊的条件下,哪一因素起到了主导作用就成为一个亟待解决的问题。表 3表明,相对湿度对蒸发量的决定作用强于气温对蒸发量的决定作用。对1961—2001年东疆三地的相对湿度分析得出(见图 3),相对湿度呈增加趋势,其线性倾向率为0.09%/10a,其变化趋势在α=0.05显著水平下达到显著。这就表明近40年降水增多明显地增加了当地空气的相对湿度,相对湿度的显著增加是导致东疆的蒸发量逐年减少的重要原因。
东疆地区蒸发量呈显著下降趋势,1970年代至1980年代初期的增多和1980年代中期至21世纪初的减少是突变现象,分别是从1972年和1990年开始的。6种气象因子单独对蒸发的决定程度按大小依次为:相对湿度、水汽压、气温、风速、日照时数、净地表辐射,风速与日照时数、相对湿度与水汽压、气温与水汽压、气温与相对湿度的协同作用对蒸发量的决定作用很大。东疆三地近40年降水增多明显地增加了当地空气的相对湿度,相对湿度的显著增加是导致东疆的蒸发量逐年减少的重要原因。
净地表辐射一直都是研究蒸发量常用的气象因子之一,而该文却得出了净地表辐射对蒸发量的决定作用很小的结论,这可能是由于净地表辐射是由经验公式计算得出的,本身误差较大,同时这可能也是在分解日照时数及净地表辐射对蒸发量的影响作用过程中出现二者符号几乎完全相反现象的原因。
使用通径分析法分析各气象因素对蒸发量的影响作用只能在线性相关的基础上进行,而某些气象因子与蒸发量的关系可能并不是线性相关,因而在研究方法上存在着一定的局限性,有待于进一步完善。
段若溪, 姜会飞, 2002. 农业气象学[M]. 北京: 气象出版社, 94-95.
|
王艳君, 姜彤, 许崇育, 2005. 长江流域蒸发皿蒸发量及影响因素变化趋势[J]. 自然资源学报, 20(6): 864-870. DOI:10.11849/zrzyxb.2005.06.009 |
杜尧东, 刘作新, 张运福, 2001. 参考作物蒸散计算方法及其评价[J]. 河南农业大学学报, 35(1): 57-61. |
魏凤英, 2007. 现代气候统计诊断与预测技术[M]. 北京: 气象出版社, 36-66.
|
毛任钊, 张妙仙, 张玉铭, 2004. 海河低平原盐渍涝洼区表层土壤积盐影响因素通径分析[J]. 中国生态农业学报, 12(2): 50-53. |
闽骞, 1991. 建立计算水面蒸发量的实用性气候学模型[J]. 气象, 17(7): 43-45. |
张敏, 宛公展, 1995. 用不同气象要素估算不同蒸发力的模式研究[J]. 气象, 21(2): 36-37. DOI:10.7519/j.issn.1000-0526.1995.02.008 |
阿依夏木·尼亚孜, 周宁芳, 杨贵名, 2007. 近45年哈密地区温度变化特征[J]. 气象, 33(7): 89-97. DOI:10.7519/j.issn.1000-0526.2007.07.011 |
张学文, 张家宝, 2006. 新疆气象手册[M]. 北京: 气象出版社, 179-180.
|
买买提·阿不都拉, 玉苏甫·阿不都拉, 刘海涛, 等, 2006. 和田市近40年蒸发量的变化特征[J]. 气象, 32(8): 92-98. DOI:10.7519/j.issn.1000-0526.2006.08.016 |