1972—2012年白洋淀湿地潜在蒸散量变化分析

贾秋兰，赵玉兵，王小娟，王 欢，吴云龙

（1邢台市气象局，河北邢台054000；2辛集市气象局，河北辛集052360）

0 引言

全球气候变暖已是不争的事实[1]，在气候变暖的大背景下，近100年来中国地表气温增加尤为明显，升温幅度约为0.5～0.8℃[2-3]。以增温为主要特征的气候变化将对中国的生态系统和社会经济产生重要影响[4-8]。

湿地作为一个独特的生态系统，在涵养水源，保护生物的多样性，调节区域气候，控制碳循环等方面具有重要的意义[9]。近些年来，一些学者对湿地气候变化开展了一些研究，取得了许多有益的成果，贾志军[10]对三江平原典型沼泽湿地蒸散量及其影响因子进行研究，认为沼泽湿地蒸散量时间变化特征明显。罗磊[11]分析认为全球气候变化是青藏高原湿地退化的重要气候背景。沃晓棠等[12]认为扎龙湿地的大气湿润度在减小，气候在向暖干方向发展。李风霞等[13]通过分析，揭示了气候变化对青海湖地区区域生态与环境的影响效应。王芳[14]研究认为七星河湿地对当地气候具有重要的调节作用。

白洋淀湿地位于河北省中部，在保定与沧州境内，是河北省第一大内陆湖，地理坐标在38°43′—39°02′N，115°38—116°26′E之间，面积366 km2，素来有“华北之肾”的称谓，但到目前，对白洋淀湿地的研究还不多，而潜在蒸散量是指充分供水条件下的区域蒸散发能力[15]。目前已被广泛应用于农田灌溉、作物需水量等水资源管理的研究中[16-17]。从气象学的角度来看，潜在蒸发量也是地表水量平衡和能量平衡的重要研究内容[18]。因此，笔者立足于白洋淀附近的安新气象站的多年气象资料，分析了近54年来白洋淀湿地潜在蒸散量的变化状况，此研究能够为白洋淀湿地的水平衡等生态保护问题提供积极地借鉴和参考。

1 资料与方法

1.1 资料来源

本研究数据来源于河北省气象局气候中心，用白洋淀湿地内的安新气象站1972—2012年逐日平均气温（含最高最低）、降水、日照、风速、相对湿度等观测数据资料，季节的划分按照3—5月为春季，6—8月为夏季，9—11月为秋季，12月至次年2月为冬季的划分方式生成逐季序列。

1.2 研究方法

1.2.1 潜在蒸散量计算公式的选取 潜在蒸散量又称参考作物蒸散量，是表征大气蒸散能力，评价气候干旱程度、植被耗水量、生产潜力以及水资源供需平衡的最重要的指标之一[19-24]，在计算该数值的过程中选取使用世界粮农组织(FAO)推荐使用的修正Penman-Monteinl方程[19]进行本区域潜在蒸散量的计算，方程式如(1)所示。

ET0是潜在蒸散量，Rn是冠层表面净辐射，冠层反射系数α取0.23；G是土壤热通量，其数值在本文忽略不计[25]；T是日平均气温；es是饱和水汽压；ea是实际水汽压；△是饱和水汽压－气温关系曲线在T处的切线斜率；γ是湿度计常数；U2是2 m高处的日平均风速。

1.2.2 采用方法 采用趋势方程、Mann-Kendall突变检验、标准差、Morlet小波分析对白洋淀湿地的地标干湿进行季节、年代的变化分析[26-27]。按春季（3—5月）、夏季（6—8月）、秋季（9—11月）、冬季（12月至2月）的季节划分生成逐季序列，采用月平均统计得到逐年平均气温。

2 结果与分析

2.1 趋势变化分析

由表1可以看到，1972—2012年白洋淀湿地年平均潜在蒸散量为993.0mm。其中春季蒸散量为330.3mm，占全年蒸散量的33%；夏季蒸散量402.8mm，占全年的41%；秋季蒸散量181.1mm，占全年的18%；冬季蒸散量79.0mm，占全年的8%。夏季潜在蒸散量最大，其次是春季和秋季，冬季最少。

表1 白洋淀湿地潜在蒸散量变化统计表（1972—2012年）

对白洋淀湿地的年、各季潜在蒸散量进行趋势分析，结合图1和表1可以看到，无论是年际还是各季节，白洋淀湿地的潜在蒸散量均呈波动下降趋势，年潜在蒸散量的趋势变化为-22.5mm/10 a，各季节中，春季的趋势变化最大，变化量为-9.5mm/10 a，冬季的趋势变化最小，变化量为-0.5mm/10 a；在41年的时间序列中，白洋淀湿地的年变化量、春季变化量、秋季变化量随时间的变化减少显著，趋势变化通过了99%的显著性检验，夏季和冬季的趋势变化没有通过任何的显著性检验；年潜在蒸散量的标准差最大，说明波动幅度最大，各季节的变化过程中，夏季的波动性最大，冬季的波动性最小。

2.2 时间序列的突变分析

Mann-Kendall突变检验简称M-K法，是一种非参数统计检验方法，可分析时间序列资料的变化趋势和突变点[26]。通过对白洋淀湿地年际与各季节潜在蒸散量的M-K检验分析图2发现，该湿地无论是年际还是各季均存在随着时间变化的突变点。年际潜在蒸散量在1983年发生突变，在1988年突破-1.96（95%显著性）的检验线，下降趋势显著；春季潜在蒸散量在1977年发生突变过程，在1986年突破-1.96的检验线，下降显著，但随着时间的变化于2009年又回归到检验线之内；夏季在1986年发生突变，在1994年和2001年2次突破-1.96的检验线，并与2001年后下降显著；秋季在1987年发生突变，在1991年下降显著；冬季在1984年发生突变，在1991—1997年间下降显著。

图1 白洋淀湿地潜在蒸散量变化趋势图（1972—2012年）

2.3 周期变化分析

在实际变化过程中，气象要素并非单调递增或者单调递减，常常在波动增减的基础上叠加了不同的变化周期，因此，采用Morlet小波方法研究了近41年来白洋淀湿地潜在蒸散量的周期性变化特征，正值代表增加，负值代表减少。在图3a可以看到，年潜在蒸散量存在8～10年的周期变化特征，在2010年左右潜在蒸散量进入减少阶段，该阶段将延伸至2014—2015年左右；春季潜在蒸散量存在5～7年左右的变化周期，潜在蒸散量在2012年后将进入增加阶段，此阶段的过程将维持3～4年；夏季存在准8年左右的周期变化特征，该周期在2011年左右进入减少阶段，预计在2015年后进入增加阶段；秋季存在特征在20世纪90年代中期以前不明显，在90年代后存在准8年左右的振荡周期，该周期在2008年左右进入减少阶段，在2012年后将进入增加的阶段；冬季存在准10年的变化周期，在2007年左右进入增加周期，预计在2012年后将进入减少的周期。

3 结论

（1）白洋淀湿地多年潜在蒸散量的年平均为993.0mm，各季节中，夏季潜在蒸散量最大，冬季最少，41年来该湿地的潜在蒸散量呈波动下降趋势，年潜在蒸散量的下降趋势为-22.5mm/10 a，各季节中，春季的趋势下降最大，变化量为-9.5mm/10 a，冬季的趋势下降最小，变化量为-0.5mm/10 a，在各季节的变化过程中，夏季的波动性最大，冬季的波动性最小。

（2）白洋淀湿地的年际、各季潜在蒸散量均存在时间序列的突变过程，其中年际变化在1988年下降趋势显著，夏季在2001年下降趋势显著，秋季在1991年下降趋势显著；春季在1986—2009年下降趋势显著。冬季在1991—1997年下降趋势显著。

（3）白洋淀湿地的年际、各季潜在蒸散量随时间的变化存在大小不同的周期变化特征，年潜在蒸散量存在8～10年的周期变化特征，春季潜在蒸散量存在5～7年左右的变化周期，夏季存在准8年左右的周期变化特征，秋季在20世纪90年代后存在准8年左右的周期特征，冬季存在准10年的变化周期。

1972—2012年，白洋淀湿地的潜在蒸散量呈减少趋势，尤其是在1988年后，其潜在蒸散量的年下降趋势呈进一步的减少状态，这对白洋淀湿地的水土保持是一个可喜的结果，对该地区的生态建设有促进的作用，本研究只分析了白洋淀湿地的潜在蒸散量的变化情况，对于湿地生态环境的其他要素的研究工作将在今后开展。

图2 白洋淀湿地Mann-Kendall突变检验分析图（1972—2012年）

图3 白洋淀湿地Morlet小波分析图（1972—2012年）

参考文献

[1]IPCC,IPCC ForthAssessmentReport(AR4)[M].Cambridge:Cambridge University press,2007:1-12.

[2]秦大河,陈振林,罗勇,等.气候变化科学的最新认知[J].气候变化研究进展,2007,3(2):63-73.

[3]丁一汇,任国玉,石广玉,等.气候变化国家评估报告(Ⅰ):中国气候变化的历史和未来趋势[J].气候变化研究进展,2006,2(1):3-08.

[4]蔡运龙,Smit B.全球气候变化下中国农业的脆弱性与适应对策[J].地理学报,1996,51(3):202-212.

[5]金之庆,方娟,葛道阔,等.全球气候变化影响我国冬小麦生产之前瞻[J].作物学报,1994,20(2):186-197.

[6]熊伟,林而达,居辉,等.气候变化的影响阈值与中国的粮食安全[J].气候变化研究进展,2005,1(2):84-87.

[7]郑景云,葛全胜,郝志新.气候增暖对我国近40年植物物候变化的影响[J].科学通报,2002,47(20):1582-1587.

[8]沙万英,邵雪梅,黄枚.20世纪80年代以来中国的气候变暖及其对自然区域界限的影响[J].中国科学,2002,32(4):317-326.

[9]张美玲,鹿翠华,张翠.近45年微山湖湿地周边大气湿度的变化特征[J].现代农业科技,2008.18:301-303.

[10]贾志军,宋长春,王跃思,等.三江平原典型沼泽湿地蒸散量研究[J].气候与环境研究,2007,12(4):496-502.

[11]罗磊.青藏高原湿地退化的气候背景分析[J].湿地科学,2005,3(3):190-199.

[12]沃晓棠,孙彦坤,玄明君,等.气候变暖背景下扎龙湿地气候变化特征[J].气象科技,2011,39(1):38-43.

[13]李凤霞,伏洋,杨琼,等.环青海湖地区气候变化及其环境效应[J].资源科学,2008,30(3):348-353.

[14]王芳,高永刚,姜春燕.黑龙江省七星河湿地近50年气候变化特征分析[J].安徽农业科学,2010,38(21):11563-11565.

[15]刘昌明,孙睿.土壤—植被—大气系统水分能量平衡研究进展[I].水科学进展,1999,10(3):23-28.

[16]Peterson T C,Golubev V S,Groisman P Y.Evaporation losing its strength[J].Nature,1995,377:687-688.

[17]Liu C M,Zeng Y.Changes of pan evaporation in the recent 40 years in the Yellow River Basin[J].Water International,2004,29(4):5l0-516.

[18]Bmtsaert W,Parlange M B.Hydrologic cycle explains the evaporation paradox[J].Nature,1998,396:30.

[19]Allen R G,Pereira L S,Raes D,et al.Crop Evapotranspiration and Drainage Paper No.56[M].Rome:FAO,1998.

[20]任国玉,郭军.中国水面蒸发量的变化[J].自然资源学报,2006,2l(1):31-44.

[21]谢贤群,王菱.中国北方近50年潜在蒸发的变化[J].自然资源学报,2007,22(5):683-691.

[22]高歌,陈德亮,任国玉,等.1956—2000年中国潜在蒸散量变化趋势[J].地理研究,2006,25(3):378-387.

[23]毛飞,张光智,徐祥德.参考作物蒸散量的多种计算方法及其结果的比较[J].应用气象学报,2000,11(增刊):128-136.

[24]刘钰,Pereira L S,Teixeira J L,等.参照腾发量的新定义及计算方法对比[J].水利学报,1997(6):27-33.

[25]吴绍洪,尹云鹤,郑度,等.近30年中国陆地表层干湿状况研[J].中国科学D辑地球科学,2005,35(3):276-283.

[26]魏凤英.现代气候统计诊断与预测技术[M].北京:气象出版社,1999:69-72.

[27]黄嘉佑.气象统计与预报方法[M].北京:气象出版社,2004:28-30.