1961年以来火山喷发年份黑龙江省温水及主要灾害变化特征初探*

闫 平，王秋京，李永生，曹 蕾，初 征，刘 丹，翟 墨，姜丽霞

(1.黑龙江省气象科学研究所，哈尔滨 150030；2.黑龙江省气候中心，哈尔滨 150030；3.哈尔滨市气象台，哈尔滨 150028)

2022年1月15日，南太平洋岛国汤加发生的火山喷发引起了世界广泛关注。多年来，世界各国科学家持续研究火山喷发对气候变化的影响，研究表明，火山活动是气候变化的驱动因素之一[1-2]，利用火山活动观测资料分析指出，小火山喷发对全球温度影响很小，影响时间也较短，而大火山喷发后的1～2年全球年均温下降约0.3℃[3-4]，同时通过数值模拟研究发现，近百年强火山喷发可造成地表年平均温度下降约0.4℃[5-7]，可见观测与模拟研究结论较为接近。许多专家学者对火山活动与我国气候变化的关系开展了研究工作[8]，张先恭等[9]分析指出火山活动对我国气候确有一定影响，火山活动对气温的影响要比对降水的影响明显得多，于革等[10]利用三维大气环流模式进行气候模拟试验，结果显示火山灰阳伞效应造成我国长江流域以北的广大地区普遍出现了0.2℃的降温，李明启等[11]分析发现10°S～10°N的强火山喷发后第一年青藏高原东部出现降温；从近540年火山喷发对中国地区降水的影响看，由于火山纬度不同以及火山喷发强度不同，火山喷发对中国地区降水量影响不同，夏季降水量变化较全年更加显著，中国不同地区降水量对火山喷发的响应也存在差异，胡昌新[12]研究了火山喷发对长江流域干旱的影响，分析提出在太阳活动低值年(黑子数最小时) 背景下，当火山猛烈喷时，将形成灰尘幕—长江流域旱情—河口咸潮入侵。

总体来看，火山活动与气候变化关系较为复杂，火山喷发性质、季节、纬度以及区域不同，其产生影响也存在差异。目前关于火山喷发对黑龙江省气候影响的相关研究较少，在火山喷发年份，黑龙江省温水特征以及低温冷害、干旱特征如何变化仍不明确，基于此，本文利用1961—2021年气象资料，对1961年以来火山喷发年份黑龙江省温水条件状况进行分析，探讨温度、降水以及灾害特点，以期为气候预测评价、灾害监测评估及农业防灾减灾提供参考。

1 资料与方法

1.1 资料

气象资料包括1961—2021年黑龙江省80个气象观测站日平均气温、日降水资料。火山喷发强弱按火山喷发指数(VEI)区分，VEI共分为8级，即VEI-1级至VEI-8级，1级是最小级别，8级是最大级别[10]。5级和5级以上是比较强烈的火山喷发，在地球上并不多见，1961年以来5级及以上火山喷发6次，共计5年(表1)[13]。汤加火山在2009年、2014年、2015年、2019年和2021年曾发生多次喷发。多次均是小规模喷发，未达到5级。

表1 5级及以上强度火山喷发概况

1.2 数据处理方法

1.2.1 数据处理分析

利用80个气象观测站日平均气温、日降水资料分别计算黑龙江省1—12月各月、春季(3—5月)、夏季(6—8月)、秋季(9—11月)、冬季(12月—次年2月)、生长季(5—9月)平均气温、平均降水合量以及≥10℃积温、≥2℃无霜期、≥0℃初霜冻日期。研究表明，火山喷发及喷发后1～2年，其对气候产生较大影响[4,9,11,14]，因此本文重点分析火山喷发年及喷发后1～2年黑龙江省温水及灾害特征。

1.2.2 距平计算方法

考虑气候变暖影响，历年平均采取30年滑动平均方式计算，即滚动计算1961—1990年平均值、1962—1991年平均值、1963—1992年平均值，……，依次类推。在进行距平计算时，1961—1990年期间的各年份均与1961—1990年平均值进行比较，从1991年开始，则利用滚动计算的30年平均值，如1991年与1962—1991年平均值进行比较，1992年与1963—1992年平均值进行比较，依次类推。采用此种方式计算，可以减少气候变化的影响，确保数据比较结果的科学性，使计算结果准确性更高。

1.2.3 低温冷害年判识

利用日平均气温资料，计算1961—2021年黑龙江省80个气象观测站5—9月平均气温和。依据国家标准《北方水稻低温冷害等级：QX/T 34967—2017》[15]划分的水稻5—9月平均气温和距平等级(表2)，判定低温冷害年。

1.2.4 干旱年判识

采用水分亏缺距平指数(CWa)作为农业干旱指数计算基础[16]，CWa计算公式为：

(1)

CW=a×CWj+b×CWj-1+c×CWj-2+d×CWj-3+e×CWj-4。

(2)

式中：CWj为第j时间单位(取10 d)的水分亏缺指数(%)；CWj-1为第j-1时间单位的水分亏缺指数(%)；CWj-2为第j-2时间单位的水分亏缺指数(%)；CWj-3为第j-3时间单位的水分亏缺指数(%)；CWj-4为第j-4时间单位的水分亏缺指数(%)；a、b、c、d、e为权重系数，a取值为0.3；b取值为0.25；c取值为0.2；d取值为0.15；e取值为0.1。

CWj计算公式为：

(3)

式中：Pj为某10 d的降水量(mm)；Ij为某10 d的灌溉量(mm)；ETcj为作物某10 d实际蒸散量(mm),由ETcj=kc×ET0计算：ET0为某10 d的作物可能蒸散量，可采用联合国粮农组织(FAO 1998)推荐的Penman-Monteith公式计算[17]。

以作物水分亏缺距平指数为指标，根据表3确定干旱等级。

表2 东北地区不同热量区域的延迟型冷害指标(℃)

表3 作物水分亏缺距平指数的等级

2 结果与分析

2.1 火山喷发年份黑龙江省温度特征

2.1.1 强火山喷发年份温度特征

强火山喷发年，黑龙江省1—12月各月平均气温总体以偏高为主，其中正距平范围在0.1～3.1℃之间，负距平为-3.1～-0.1℃，正距平月份较负距平月份偏多28.0%。2月、8月未出现负距平，1月、3月、5月、9月、10月、12月出现3年负距平，其它月份出现1～2年负距平(图1)。

图1 强火山喷发年黑龙江省1—12月平均气温距平

从季节变化看，强火山喷发5年夏季、生长季均呈现正距平，以暖为主导，冬季总体偏暖，春季、秋季冷暖不一。无霜期有3年偏少1～3 d，初霜冻日期2年提前，≥10℃积温2年略偏少(表4)。

表4 强火山喷发年黑龙江省四季、生长季温度条件

进一步分析强火山喷发后1～2年温度变化情况，强火山喷发后1年黑龙江省1—12月各月正距平为0.1～3.2℃，负距平为-4.0～-0.1℃，月平均气温正距平月份较负距平月份偏少9.7%，6月、10月、11月出现4年负距平，2月、8月、9月出现3年负距平，其它月份出现1～2年负距平(图2a)；强火山喷发后2年1—12月各月正距平为0.1～4.4℃，负距平为-4.2～-0.1℃，正距平月份与负距平月份持平，4月出现4年负距平，1—3月出现3年负距平，其它月份出现1～2年负距平(图2b)。

图2 强火山喷发后1～2年黑龙江省1—12月平均气温距平

强火山喷发后1年，春季5年均为正距平，呈偏暖趋势，夏季、生长季出现4年负距平，总体偏冷，秋季和冬季冷暖并存，≥10℃积温以偏少为主，无霜期1年缩短，初霜冻日期3年提前；强火山喷发后2年，冷、暖并存在主导趋势，春季和冬季出现3年负距平，其它季节出现1～2年负距平，生长季平均气温和≥10℃积温1年呈现负距平，无霜期1年缩短，初霜冻日期2年提前(表5)。

2.1.2 汤加火山喷发年份温度分析

汤加火山喷发5年正距平为0.1～4.6℃，负距平为-2.8～-0.1℃，正距平月份多于负距平月份，偏多43.5%，2009年负距平出现最多为7个月，其它年份出现2～6个月负距平(图3a)；汤加火山喷发后1年1—12月各月正距平为0.1～4.0℃，负距平为-4.9～-0.1℃，正距平月份多于负距平月份，偏多53.8%，2020年负距平出现最少为2个月，2010年、2016年分别出现6和5个月负距平(图3b)。汤加火山喷发5年和后1年中2009年3个季节出现负距平，2014年未出现负距平，其他年份均是1～2个季节出现负距平，生长季平均气温2009年和2019年偏低0.1～0.2℃，≥10℃积温1年偏少，无霜期正常或偏长，初霜冻日期2年提前。汤加火山喷发后1年生长季平均气温和≥10℃积温均偏多，无霜期偏长，初霜冻日期2年提前(表6)。

表5 强火山喷发后1～2年黑龙江省四季、生长季温度条件

表6 汤加火山喷发年及后1年黑龙江省四季、生长季温度条件

图3 汤加火山喷发年及后1年黑龙江省1—12月平均气温距平

2.2 火山喷发年份黑龙江省降水量特征

2.2.1 强火山喷发年份降水量特征

强火山喷发年，黑龙江省1—12月各月降水量总体以偏少为主，其中正距平为2.8%～93.5%，负距平为-72.2%～-1.1%，负距平月份较正距平月份偏多18.2%。出现5年负距平的月份是11月，出现4年负距平的月份是2月和3月，其它月份出现1～3年负距平(图4)。

图4 强火山喷发年黑龙江省1—12月降水量距平

从季节变化看，强火山喷发年春季降水均为正距平，夏季、冬季、生长季出现3年负距平，秋季出现2年负距平(表7)。

表7 强火山喷发年黑龙江省四季、生长季降水距平

强火山喷发后1年正距平月份与负距平月份接近，但月份间降水距平差异较大，1983年4月降水距平高达226.2%，1964年2月降水距平少至-85.5%，二者相差311.7%，1992年负距平月份最多，占全年66.7%，其它年份33.3%～58.3%的月份出现负距平。6月未出现负距平，1月、2月出现4～5年负距平，其它月份出现1～3年负距平(图5a)；强火山喷发后2年正距平月份较负距平月份偏多13.0%，正距平月份与负距平月份差距较大，1965年负距平月份最多，占全年75.0%，其它年份25.0%～41.7%的月份出现负距平，1月、4月、5月、12月出现3年负距平，其它月份出现0～2年负距平(图5b)。

图5 强火山喷发后1～2年黑龙江省1—12月降水量距平

表8 强火山喷发后1～2年黑龙江省四季、生长季降水距平

表9 汤加火山喷发年及后1年黑龙江省四季、生长季降水距平

强火山喷发后1年，夏季、冬季、生长季降水以偏多为主导趋势，春季、秋季则出现2～3年负距平。强火山喷发后2年，夏季、秋季、生长季降水总体偏多，而春季和冬季出现2年负距平(表8)。

2.2.2 汤加火山喷发年份降水分析

汤加火山喷发5年降水量正距平月份多于负距平月份，偏多19.2%，2014年负距平月份最多，占全年66.7%，其它年份25.0%～50.0%的月份出现负距平，4月出现5年负距平，10月出现4年负距平，1月、3月出现3年负距平，其它月份出现0～2年负距平(图6a)。汤加火山喷发后1年正距平月份较负距平月份偏多109.1%，每年25.0%～33.3%的月份出现负距平，7月、12月出现2年负距平，其它月份出现0～1年负距平(图6b)。汤加火山喷发5年中0～2个季节出现负距平，生长季降水量均偏多。汤加火山喷发后1年中的1-2个季节出现负距平，生长季降水量为偏多趋势(表9)。

图6 汤加火山喷发年及后1年黑龙江省1—12月降水量距平

2.3 火山喷发年份灾害发生特征

2.3.1 低温冷害

强火山喷发年黑龙江省ΔT5-9>0，均未发生低温冷害；强火山喷发后1年中1964年、1981年、1983年发生严重低温冷害，1992年ΔT5-9为-1.0，未达到低温冷害程度，为低温年，2012年ΔT5-9>0，未发生低温冷害；强火山喷发后2年中仅1965年发生轻度低温冷害，其它年份ΔT5-9>0(图7)。由此可见，强火山喷发年黑龙江省低温冷害发生频率为0%，强火山喷发后1年低温冷害发生频率较大，为60%，强火山喷发后2年低温冷害发生频率为20%。

汤加火山喷发年均未达到低温冷害标准，其中2009年和2019年ΔT5-9<0，为低温年，2014年、2015年、2021年为正常或高温年；汤加火山喷发后1年均为正常或高温年(图7)。可见在汤加火山喷发年及次年，黑龙江省发生低温冷害频率为0%。

图7 黑龙江省5—9月平均气温和距平

由图8a至图8f可见，强火山喷发年低温冷害发生频率明显低于历年发生频率；强火山喷发后1年低温冷害发生频率明显高于历年；强火山喷发后2年低温冷害发生频率大部地区低于历年，西部和东南部少部分地区低温冷害发生频率略高于历年。汤加火山喷发年除了个别站点外，其它地区低温冷害发生频率低于历年平均状况，汤加火山喷发后1年全省低温冷害发生频率均低于历年平均状况。

图8 黑龙江省低温冷害发生频率分布图(审图号GS(2017)3320号，底图无修改)

图9 水分亏缺指数距平变化

2.3.2 干旱灾害

从月尺度来看，强火山喷发年、强火山喷发后1～2年黑龙江省发生干旱的月份较少，其中强火山喷发年1982年6月、7月出现轻旱，其发生频率为20.0%，分别较1961—2021年同期偏高15.0%、18.0%；强火山喷发后1年的1964年5月份发生轻旱，发生频率为20.0%，比1961—2021年同期偏高10.0%；强火山喷发后2年的1965年5月、6月达到轻旱等级，发生频率为20.0%，分别较1961—2021年同期偏高10.0%、15.0%。进一步分析干旱日变化发现，1982年主要为夏季干旱，6月12日—7月17日作物水分亏缺指数距平大于40.0%，持续干旱，其中6月30日—7月9日大部时段为重旱，最大值高达72.0%，旱情严重；1964年为春旱，5月4日—6月4日大部时段作物水分亏缺指数距平大于40.0%，但未达到重旱程度，干旱持续时间及影响程度明显弱于1982年；1965年为春旱和初夏旱，5月1—20日及5月31日—6月29日大部时段作物水分亏缺指数距平大于40.0%，干旱持续时间较长，其中5月3日达到重旱等级，最大值为70.0%，干旱影响程度弱于1982年，强于1964年(图9)。在生长季尺度上，强火山喷发年、强火山喷发后1年、强火山喷发后2年黑龙江省均未发生干旱。汤加火山喷发的5年及喷发后1年各月及生长季均未达到干旱程度。

3 讨论

本文研究结果表明强火山喷发后1年黑龙江省生长季温度呈现低温特征，夏季低温明显，低温冷害发生概率大，程度重，这与李平原等[4]、李明启等[11]、聂文政等[14]、WOOKDS等[18]研究结果一致，大量研究已明确造成低温的原因，即强火山喷发导致大量硫化物气体与水蒸汽进入平流层，半年至一年后逐渐转化形成气溶胶，大量的气溶胶扩散至周围区域甚至全球，反射和折射太阳辐射，进而使太阳对地表的直接辐射量减少，导致地表温度降低，低温多发生在强火山喷发后1年[4-5，19]。另外，本文分析表明强火山喷发后2年低温冷害发生概率不确定性较大，与李靖等[20]研究的强火山喷发后地面温度偏低持续1～3年略有差别，分析其原因，由于火山喷发对温度影响随时间、空间不同而存在差异，火山喷发对其附近的温度影响较明显[20]，而黑龙江省距离强火山较远，因此低温影响时间与全球平均状况存在差异。本研究还表明强火山喷发年、汤加火山喷发年和后1年，生长季温度呈现正常或偏高特征，低温冷害发生概率小，这与WIGNALL[21]研究结论相吻合，其分析指出火山喷发造成大气圈CO2浓度增加，引起温室效应，而强火山喷发短时间内以及低级别火山喷发并不能形成大范围气溶胶，因此不会造成低温，而温室效应却导致气温升高。

本研究中强火山喷发年、喷发后1～2年降水量及干旱趋势与张先恭等[9]研究结果相近，与张富国等[22]研究的强火山喷发对华北地区旱涝影响差异较大，已有研究结果显示世界不同火山区域强火山喷发对我国华南地区早涝影响亦不同，甚至相反[23]。从目前有限数据分析结果来看，强火山喷发对黑龙江省降水、旱涝影响及影响幅度存在不确定性。

4 结论

(1)强火山喷发年温度正常或偏高，作物生长季热量条件较好；强火山喷发后1年，生长季温度呈现低温特征，夏季低温明显，初霜冻提前，但无霜期有增加趋势；强火山喷发后2年，生长季温度呈现正常偏高特征，春季、冬季温度表现为低温特征，初霜冻正常或提前，无霜期有增加趋势。汤加火山喷发年和后1年，生长季温度正常或偏高，初霜冻正常或提前，无霜期有增加趋势。

(2)强火山喷发年春季降水略多，夏季、冬季、生长季降水呈现偏少趋势，作物生长期间水分条件略差；强火山喷发后1～2年，生长季降水充足，夏季有偏多趋势。汤加火山喷发年和后1年，月降水量分布不均，生长季总降水量呈现偏多趋势。

(3)强火山喷发年黑龙江省发生低温冷害概率小；强火山喷发后1年低温冷害发生概率大，程度重；强火山喷发后2年低温冷害发生存在不确定性，西北局部和东南局部发生低温冷害的概率略呈增加趋势，其它地区持平或呈减少趋势。在汤加火山喷发年及后1年，黑龙江省未发生低温冷害。

(4)强火山喷发年6月、7月发生干旱的概率增加，干旱发生程度较重；强火山发生后1年，5月发生干旱概率增加；强火山发生后2年，5月、6月发生干旱概率增加；汤加火山喷发年及后1年发生干旱概率低。

(5)受数据量限制，本文研究结果普适性仍需进一步验证，后期应继续收集火山喷发相关信息，结合气象数据、相关模拟模型进行深入研究。