1961-2010年中国参考作物蒸散量变化趋势与时空格局

黄 娟,申双和,李新建,张文斌

(1.南京信息工程大学 应用气象学院,南京 210044 ; 2.南京信息工程大学 气象灾害预警预报与评估协同创新中心, 南京 210044; 3.新疆维吾尔自治区农业气象台,乌鲁木齐 830002; 4.中国气象局 乌鲁木齐沙漠气象研究所,乌鲁木齐 830002)



1961-2010年中国参考作物蒸散量变化趋势与时空格局

黄 娟1,4,申双和2,李新建3,张文斌4

(1.南京信息工程大学 应用气象学院,南京 210044 ; 2.南京信息工程大学 气象灾害预警预报与评估协同创新中心, 南京 210044; 3.新疆维吾尔自治区农业气象台,乌鲁木齐 830002; 4.中国气象局 乌鲁木齐沙漠气象研究所,乌鲁木齐 830002)

根据甄选的中国529个气象台站1961—2010年的观测资料，利用联合国粮食及农业组织(FAO)推荐的Penman-Monteith模型估算了全国年际参考作物蒸散量，探究了中国1961—2010年参考作物蒸散量变化趋势与时空格局。结果表明：(1)通过对中国参考作物蒸散量的年变化趋势进行分析，发现1961—2010年参考作物蒸散量先有下降趋势(1961—1993年)，下降幅度不大，后又有缓慢的回升(1994—2010年)，但总体呈下降趋势；(2)通过对各站近50年参考作物蒸散量平均值的时空分布进行分析，发现中国西北地区和西南地区明显大于东北地区和中部腹地；(3)中国参考作物蒸散量时空分布主要影响因素是风速和气温(气温影响范围广，风速影响程度大)，也与中国地形复杂、面积广阔、经纬跨度大、各地气象条件差异导致作物蒸散能力差异有关。

Penman-Monteith模型; 参考作物蒸散量; 变化趋势; 时空格局

蒸散即植物叶面蒸发(蒸腾)和棵间土壤蒸发之和。蒸散是植物失水的主要方式。联合国粮食及农业组织(FAO)将参考作物蒸散定义为“假设高度为0.12 m，冠层阻力为7.0 m/s，反照率为0.23的参考冠层的蒸散，即生长旺盛，长势一致，完全覆盖地面且水分供应充足的开阔绿色草地的蒸散”[1]。参考作物蒸散(ET0)是自然界水循环中的一个重要组成因素，其数值大小直接影响到作物产量的估算、农业用水的安排、各地灌溉制度的制定和选种育种的区域划分。通过分析不同年份不同地区的参考作物蒸散量时空变化趋势，可以更好地服务于农业、牧业、林业等各个行业。

国际上对自然蒸发有文献报道的研究已有300多年的历史。Dalton等[2]提出计算蒸发速率的公式：水面蒸发速率与蒸发面的饱和差成正比。Wilm[3]和Penman[4]等先后提出潜在蒸散的概念及相应的计算公式。Penman公式为计算湿润下垫面潜在蒸散的主要方法。Monteith等[5]通过引入表面阻力的概念，推导出Penman-Monteith公式，所以该公式又被称为大叶模型。1998年FAO出版计算作物蒸散的指南[6]。

国内蒸散研究工作始于20世纪50年代，朱岗昆等[7]应用Penman方法计算了自然状况下中国蒸发量的分布、季节变化等。陶祖文等[8]提出了用气候资料计算旬、月可能蒸散和土壤水分不足情况下蒸散量的方法。自80年代以来，王健等[9]通过对不同作物农田蒸散试验观测研究，建立了适用于我国的农田蒸发量测定方法和计算模式。辛晓洲[10]提出了以热红外遥感信息为主的遥感地表水分蒸发模型。刘钰[11]对FAO推荐的作物系数计算方法进行了验证。

黄会平等[12]根据1957—2012年全国608个气象站的逐日气象资料利用Penman-Monteith公式计算了作物的潜在蒸散量，采用主成分分析方法对10大水资源一级分区潜在蒸散量的年代际变化进行了分析，在ArcGIS中将4个主成分分别表示在空间分布图上。但是，阎慈琳[13]研究证明采用主成分分析方法：(1)首先应保证所提取的前几个主成分的累计贡献率达到一个较高的水平，即保证相关的气象因素对该地区潜在蒸散量有较大影响。这一条件并不是所有地区都满足，且有地域差异。(2)主成分的解释其含义一般多少带有模糊性，不像原始变量的含义那么清楚、确切，因此采用这一方法，只能表示某一站点某一气象因素有无影响，及这一气象因素在这一站点相对于全国而言的影响程度，而不能准确表示对于这一站点某一气象因素的影响程度是多少。(3)这一方法是将一维数据做出的平面结果平铺到GIS中，具有空间局限性。

综合国内外关于我国参考作物蒸散量分布时空变化的研究，由于数据不完整[14]、区域划分不明确[15]、统计量大、计算繁杂等原因，大部分研究只反映了局部地区的参考作物蒸散量时空变化[16]，或只反映出近几年的变化趋势[17]，缺乏完整的全国范围内时空变化趋势，缺乏不同地区参考作物蒸散量的时空演变趋势对比，缺乏长期变化趋势分析[18]。虽然黄会平等[12]的文章数据充足，空间范围大，但是并不能准确表示某一站点某一气象因素的影响程度是多少。此外，一维数据具有空间局限性，而采用空间叠加能解决这一问题。本文将分析全国529个站点1961—2010年的数据，以便从时间的长期性和站点数据的完整性两方面弥补之前研究的不足，然后对每一个站点的影响因素做相关性分析，以便能明确、清晰地表示出每一个站点每一个气象因素的影响程度是多少，继而采用空间叠加的方法，以便表示出全国范围内每个站点潜在蒸散的主导影响因子，从而对全国蒸散量演变形势与空间格局有更全面准确的了解。

1　试验材料与方法

1.1材料来源

分析所需的数据主要有如下两个方面：(1)529个气象站1961—2010年逐日基本气象要素资料，包括日平均气温(℃)、最高最低气温(℃)、日照时数(h)、平均风速(m/s)等；(2)1∶400万的中国行政边界图层(由国家基础地理信息系统网站提供)。所选站点的分布见图1。

图1　中国气象站点分布

1.2研究方法

1.2.1FAO Penman-Monteith模型参考作物蒸散量的计算方法多种多样。其中Penman-Monteith公式是基于能量平衡和水汽扩散理论的算法，在充分了解植被生理特性的基础上，结合空气动力学参数的变化，计算出某一地区的参考作物蒸散量。由于理论依据充分，应用性广泛，计算误差小，因此，该方法于1998年被联合国粮农组织(FAO)推荐为计算参考作物蒸散量的标准方法，直至今日仍然在各个地区得到广泛的应用。本文也利用该公式进行全国参考作物蒸散量的计算和分析。

(1)

式中：ET0(FAO-98)为参考作物蒸散量(mm/d)；Rn为地表净辐射[MJ/(m2·d)]；G为土壤热通量[MJ/(m2·d)](在逐日计算式中，G≈0)；T为2 m高度处平均气温(℃)；U2为2 m高度处风速(m/s)；es为饱和水汽压(kPa)；ea为实际水汽压(kPa)；Δ为饱和水汽压曲线斜率(kPa/℃)；V为干湿表常数(kPa/℃)。其中：

(2)

(3)

Rso=(0.75+2z×10-5)Ra

(4)

式中：T为反射率(无量纲)(本文取T=0.23)；Rs为地表短波辐射[MJ/(m2·d)]；e为斯蒂芬—波尔兹曼常数4.903×10-9[MJ/(K4·m2·d)]；Tmax为日最高气温(K)；Tmin为日最低气温(K)；Rso为晴天地表短波辐射[MJ/(m2·d)]；as为阴天(n=0)地表短波辐射与大气层外太阳辐射的比例系数；(as+bs)为晴天(n=N)地表短波辐射与大气层外太阳辐射的比例系数；n为实际日照时数(h)；N为最大可能日照时数(h)；Ra为大气层外太阳辐射[MJ/(m2·d)]；z为海拔高度(m)。本文取as=0.25，as+bs=0.75。

1.2.2气候倾向率估算

(5)

1.2.3要素相关性判别标准在计算该序列的相关系数r的值时，在常规显著性水平θ下，max{|r|θ}<0.6(r0.01,18=0.561 4小于0.6)，即此时所有阈值中的最大值均小于0.6。因此，当相关系数满足rβ≤|r|<0.6(且序列长度大于最小值20时)，序列与时序t的相关性要高于上面已划分的等级，本文将其趋势变异程度划为低度相关，以此类推得到表1。

表1　相关系数分级[19]

1.2.4参考作物蒸散量资料插补方法在所选529个站点中，并不是所有站点的资料都具有完整性，所以要对数据缺失年份的参考作物蒸散量进行插补。本文所需数据为数值型，经过比较之后，采用均值替换法。

2　结果与分析

2.1全国年平均参考作物蒸散量时间变化趋势

1961—2010年我国参考作物蒸散量(图2)总体呈下降趋势，平均每年下降0.279 7 mm，下降幅度较小。然而不同年份的变化情况有所不同，总体来看，1961—1968年呈上升趋势，1969—1991年总体呈下降趋势，下降幅度较大，且均处于趋势线以下，1992—2010年呈上升趋势，且上升幅度较大。参考作物年蒸散量最大值出现在1978年，最大值为1 007.13 mm，最小值出现在1993年，数值为926.06 mm。

2.2年参考作物蒸散量空间分布与时空格局

我国1961—2010年多年平均参考作物蒸散量(附图4)总体分布为西部多，东部少。其中，我国西北地区塔里木盆地蒸散量较大，吐鲁番盆地达到最大值。海南岛及雷州半岛蒸散量达到最小值。1961—2010年我国年参考作物蒸散量(附图5)呈西南部增加，西北部和东南大部减少。

图2　1961-2010年年平均参考作物蒸散量变化趋势

2.3年参考作物蒸散量年代际变化

由附图6可知，与1960s相比，1970s我国大部分地区(青藏高原东部、内蒙古高原中部和西部、东北平原北部)蒸散量呈增加趋势，增加数值为27～49 mm；另外，准噶尔盆地、长江中下游地区以南、华北平原北部等地蒸散量均减少18～42 mm。由附图7A可知，与1970s相比，1980s我国大部分地区蒸散量有减少趋势，减少数值为东部(长江中下游南部地区、东北平原大部)地区0～17 mm，西部(塔里木盆地东部)地区39～49 mm。由附图7B可知，与1980s相比，1990s我国大部分地区蒸散量有减少趋势，减少量为0～20 mm；另外，中部腹地、东北北部有增加趋势，增加量为22～37 mm。由附图7C可知，与1990s相比，2000s我国大部分地区蒸散量有增加趋势，增加量为1～32 mm。

2.4参考作物蒸散量时空变化原因分析

为了分析出可能对参考作物蒸散量变化有影响的因素，对1961—2010年全国的多年平均参考作物蒸散量、气温、相对湿度、风速、日照时数做出了相应的相关系数等级分析，对比之后取每个站点的最大值(主要影响因素)，判断其影响程度的大小。

气温对我国参考作物蒸散量的影响程度空间分布如下(图3A)：东北平原东部、云贵高原低度相关；东北平原西部、东南丘陵大部、青藏高原东部、海南岛显著相关；准噶尔盆地北部、华北平原局部高度相关。相对湿度对我国参考作物蒸散量的影响程度空间分布如下(图3B)：长江中下游沿江地区、东南丘陵大部、黄土高原大部无相关关系；云贵高原大部、青藏高原东部显著相关。风速对我国参考作物蒸散量的影响程度空间分布如下(图3C)：黄土高原北部、云贵高原东部、东南丘陵局部低度相关；其他大部分地区为显著相关。日照时数对我国参考作物蒸散量的影响程度空间分布如下(图3D)：东北平原南部、长江中下游局部、准噶尔盆地西北部低度相关；东北平原西北部、四川盆地、青藏高原东部的局部地区显著相关；其他大部分地区无相关关系。

图3　1961-2010年中国ET0与气温、相对湿度、风速、日照时数的相关性

1961—2010年我国参考作物蒸散量的主要影响因素为(图4)：就全国范围来看，主要受气温和风速的影响，这二者相比，气温影响范围更广，风速影响程度更大。东北北部主要受日照影响，东北南部主要受气温和风速影响，其中气温为主导因素；辽东半岛主要受相对湿度影响；华北地区主要受气温和风速影响；长江中下游及其以南的地区主要受风速和气温影响，其中风速为主导因素；云贵高原主要受相对湿度影响；黄土高原西部主要受日照影响；西北地区主要受气温影响。

3　结 论

1961—2010年，我国年参考作物蒸散量总体呈下降趋势，平均每年下降0.279 7 mm，下降幅度较小。年平均参考作物蒸散量分布西部多，东部少。具体分布如下：我国西北地区塔里木盆地蒸散量较大，吐鲁番盆地达到最大值，内蒙古高原西部及中部蒸散量达到次大值。东北地区蒸散量较小，尤其东北平原以北(包括小兴安岭及大兴安岭东北部)的地区数值更小。海南岛及雷州半岛蒸散量达到最小值。

图4　1961-2010年中国ET0的主要影响因素分布

对参考作物蒸散量的时空变化格局分析发现，1961—2010年，我国中部地区(包括青藏高原、吐鲁番盆地、内蒙古高原东部、西川盆地大部、秦岭淮河地区、大兴安岭以西的地区)的参考作物蒸散量有明显的增加趋势，其他地区则有不同程度的减小。

另外，近50年来，全国范围内的日照时数显著减小，受此影响净辐射亦大幅度减小，因此，作为蒸散能量来源的辐射能的减少，一方面导致蒸散能力减弱，另一方面导致大气热力层结较稳定，湍流减弱，使蒸散量减少；近50年来，相对湿度减小，因而空气较干燥，最高温度有较大幅度的上升，这就使得蒸散量有回升的趋势。

本文对我国参考作物蒸散量逐年的变化趋势和空间分布做出了详细分析，但仍有不足之处。(1)本文只是在FAO定义的Penman-Monteith模型基础上做了相关分析，日后可以结合其他模型作对比分析；(2)本文只是根据年际变化来划分时间区域，日后可根据作物的不同生育期来划分时间区域；(3)还可通过对比不同作物的蒸散量作为参照，通过对比分析得出结论；(4)近几年遥感、GIS等空间模拟网络技术快速发展，遥感对空间的模拟更精确，GIS对数据的插值处理更合理，日后可结合这些技术作相关分析。

[1]Allen R G,Pereira L S,Raes D,et al.FAO Irrigation and Drainage Paper 56: Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements[M].Rome: FAO,1998.

[2]Dalton A J.Electron micrography of epithelial cells of the gastro-intestinal tract and pancreas[J].American Journal of Anatomy,1951,89(1):109-133.

[3]Wilm H G,Thornthwaite C W,Colman E A,et al.Report of the committee on transpiration and evaporation,1943-44[J].Eos,Transactions American Geophysical Union,1944,25(5):683-693.

[4]Penman H L.Natural evaporation from open water,bare soil,and grass[M]∥Proceedings of the Royal Society of London A.the Royal Society Publishing,1948.

[5]Monteith J L.Evaporation and environment.[J].Symposia of the Society for Experimental Biology,1965,19(19):205-34.

[6]FAO.World reference base for soil resources[R].Rome: FAO,1998.

[7]朱岗昆,杨纫章.气象记录在经济建设中的应用(Ⅱ):中国各地蒸发量的初步研究[J].气象学报,1955,26(1/2):1-28.

[8]陶祖文,裴步祥.农田蒸散和土壤水分变化的计算方法[J].气象学报,1979,37(4):79-87.

[9]王健,蔡焕杰,刘红英.利用Penman-Monteith法和蒸发皿法计算农田蒸散量的研究[J].干旱地区农业研究,2002,20(4):67-71.

[10]辛晓洲.用定量遥感方法计算地表蒸散[D].北京:中国科学院研究生院(遥感应用研究所),2003.

[11]刘钰.气象数据缺测条件下参照腾发量的计算方法[J].水利学报,2001,32(3):11-17.

[12]黄会平,曹明明,宋进喜,等.1957—2012年中国参考作物蒸散量时空变化及其影响因子分析[J].自然资源学报,2015,30(2):315-326.

[13]阎慈琳.关于用主成分分析做综合评价的若干问题[J].数理统计与管理,1998,17(2):22-25.

[14]杜加强,舒俭民,刘成程,等.黄河上游参考作物蒸散量变化特征及其对气候变化的响应[J].农业工程学报,2012,28(12):92-100.

[15]李志.黄土高原1961—2009年参考作物蒸散量的时空变异[J].生态学报,2012,32(13):4139-4145.

[16]于东平,张鑫,何毅,等.青海省东部高原农业区参考作物蒸散量的时空变化[J].农业工程学报,2012,28(2):66-71.

[17]刘勤,严昌荣,梅旭荣,等.西北旱区参考作物蒸散量空间格局演变特征分析[J].中国农业气象,2012,33(1):48-53.

[18]刘广东,李艳,刘海军,等.山西地区近55年参考作物蒸散量的变化特征及其主要影响因素分析[J].灌溉排水学报,2012,31(5):26-30.

[19]谢平,唐亚松,李彬彬,等.基于相关系数的水文趋势变异分级方法[J].应用基础与工程科学学报,2014,22(6):1089-1097.

Research of Variation Trend and Spatial Pattern of Reference Crop Evapotranspiration During the Period from 1961 to 2010 in China

HUANG Juan1,4,SHEN Shuanghe2,LI Xinjian3,ZHANG Wenbin4

(1.College of Applied Meteorology,Nanjing University of Information Science & Technology, Nanjing 210044,China; 2.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disaster, Nanjing University of Information Science & Technology,Nanjing 210044,China; 3.Xinjiang Agriculture Observatory, Urumqi 830002,China; 4.Institute of Desert Meteorology,China Meterological Administration,Urumqi 830002,China)

The annual reference crop evapotranspiration was estimated using Penman-Monteith model recommended by Food and Agriculture Organization of the United Nations (FAO)with the observation data from 529 meteorological sites during the period of 1961 to 2010,and the variation trend and spatial pattern of ET0were investigated in detail.The results indicated that: (1)through analysis of annual variation tendency for the reference crop evapotranspiration,it was found that during the period from 1961 to 2010,reference crop evapotranspiration declined (1961—1993),did not drop,and then slowly recovered (1994—2010); (2)through the reference crop evapotranspiration in recent 50 years,the average value on the temporal and spatial distribution analysis,it was found that the values in northwest China and southwest China were significantly greater than those of northeast area and the central hinterland; (3)the main factors affecting the spatial and temporal reference crop evapotranspiration were the wind speed and temperature distribution in China (Temperature had a wide range of influence,wind speed had a great effect on ET0).And meanwhile,China′s complex terrain,the vast area,latitude and longitude span,all weather conditions caused the differences in different crop evapotranspiration capacity.

Penman-Monteith model; reference crop evapotranspiration; change trend; temporal and spatial distribution

2015-09-14

2015-09-28

公益性行业(气象)科研专项“干旱气象科学研究——我国北方干旱致灾过程及机理”(GYHY201506001);国家重大科学研究计划“干旱半干旱区生态系统和水资源脆弱性评估及风险预估”(2012CB956204)

黄娟(1990—),女,新疆塔城人,硕士研究生,研究方向为农业气象研究。E-mail:1321265319@qq.com

申双和(1957—),男,江苏姜堰人,教授,博士生导师,主要从事农业气象研究。E-mail:yqzhr@nuist.edu.cn

S161.4; P311

A

1005-3409(2016)05-0240-05