2017年全球活动火山时空分布以及阿贡火山喷发前后的形变特征1

王佳龙 王家庆 邹 颖 孙嘉祥 马庆尊

2017年全球活动火山时空分布以及阿贡火山喷发前后的形变特征1

王佳龙1，2）王家庆1）邹 颖2）孙嘉祥2）马庆尊1）

1）中国地震局第一监测中心，天津 300180 2）中国地震局地质研究所，北京 100029

归纳总结2017年度全球81座活火山的活动情况，共计活动1058座次，平均每周记录20座活火山的活动信息。根据火山潜在喷发的危险性和火山活动的强弱程度对上述火山进行分级描述，火山活动主要反映了地球表层的构造活动，其中大角度俯冲带的弧后火山最为强烈，小角度的俯冲带、拉张裂谷和走滑为主的板块边界火山活动较为平静，火山活动频繁的印度尼西亚岛链是受灾最为严重的区域。预计全球火山活动将进一步加剧，印尼岛链受火山灾害威胁的程度依然较大。位于印尼岛链巴厘岛上的阿贡火山自2017年9月开始活动以来，整个喷发过程极具代表性，监测阿贡火山喷发过程可为全球典型火山喷发事件研究提供参考。

活火山 火山活动 构造活动 火山灾害 阿贡火山

引言

地质学家已经确定全球范围的全新世火山约有1500多座（许建东，2011）。火山类型主要为成层火山、复合式火山、熔岩锥、破火山口、盾状火山等（刘若新等，1999）。此外，还有众多火山位于深海海底，它们的喷发很少能到达水面，因此少有观察和记录。大部分火山的形成是板块运动的结果，是地球动力学过程的重要现象（Press等，1982）。全球现今活动构造划分为环太平洋、大洋中脊和大陆三大构造系统（马宗晋等，2003）。全球火山活动分西太平洋火山活动区、东太平洋火山活动区以及大西洋火山活动区（洪汉净等，2003，2009）。位于太平洋板块西南缘的印度尼西亚、菲律宾和日本等国家，拥有的火山数量为全球正在活动火山数量的1/3（Siebert等，2010），2017年该区域活动强度最大的火山为印度尼西亚巴厘岛的阿贡火山（王佳龙等，2018a）。

本文数据主要来自全球火山计划网站[1]，该网站由华盛顿国家自然历史博物馆矿物科学部门（Department of Mineral Sciences，National Museum of Natural History，Washington D.C）的史密森全球火山项目（Smithsonian Institution Global Volcanism Program）支撑，汇集了世界各地活动火山监测机构的监测信息及相应的火山监测网站链接，并于每周四发布监测周报。文中的活动火山指2016年12月30日—2017年12月31日发生活动并被该网站记录下来的火山。

通过整理该网站2017年发布的约858条全球火山监测活动信息，总结2017年度全球火山活动的时空关系，并针对阿贡火山的喷发过程进行了追踪，为监测火山喷发、研究火山活动机理、预测火山灾害等提供参考。

1 2017年全球活动火山空间分布特征

火山警戒等级是衡量火山活动危险程度的标识，USGS火山灾害项目[2]将火山预警等级划分为正常（Ⅰ级）、咨询（Ⅱ级）、注意（Ⅲ级）和警戒（Ⅳ级）4个等级，其代表的危险程度依次上升，火山活动越强烈，引发的火山灾害也越严重。具体划分标准为：Ⅰ级，火山警戒级别为正常，航空颜色代码为绿色，火山有活动迹象，但活动处于正常范围内，主要表现为火山喷气及火山灰柱飘散；Ⅱ级，火山警戒级别为低度警戒，航空颜色代码为黄色，火山活动信号上升，高于正常范围，主要表现为火山出现热异常，火山喷气柱升高，出现火山爆炸声，火山由正常活动向危险活动过渡，该类型火山一般不稳定；Ⅲ级，火山警戒级别为中度警戒，航空颜色代码为橙色，火山活动信号进一步上升，具有喷发的前兆，主要表现为火山地震频率和强度增大，可以监测到熔岩流溢出，热异常和火山灰柱活动都有所加强；Ⅳ级，火山警戒级别为高度警戒，航空颜色代码为红色，表现为火山即将喷发或正在喷发。一般情况下，火山的警戒等级随火山活动的强弱而改变，可能会在数月内出现连续变化（张传杰等，2016）。

图1 全球活动火山分布

全球活动火山约有80%分布在环太平洋地区（Siebert等，2010）。2017年度的活动火山约有90%分布在环太平洋地区，其中，太平洋北部阿留申群岛、堪察加半岛和北方四岛等分布12座活动火山，其警戒等级为Ⅱ级和Ⅲ级，警戒等级Ⅲ级的活动火山相对比较活跃；位于太平洋板块西缘的日本岛弧和东部海岛分布11座活动火山，警戒等级为Ⅰ级的火山相对比较活跃；位于澳大利亚板块北缘的印度尼西亚等一系列岛链，由澳大利亚板块向北俯冲到缅甸板块、欧亚板块和太平洋板块之下形成，构造复杂，火山活动能量大，2017年此区域有21座活动火山，包括3座警戒等级Ⅲ级的活动火山，可见其活动频度和强度都极高。

全球火山分布图（图1）的底图来源于USGS，对其做罗宾逊投影变换，中央子午线设为155.0°E。若某火山2017年每周均有持续活动，则将该火山的活动频率设为1；若火山全年无任何活动，则将火山的活动频率设为0。火山的活动频率一般介于0和1之间，它反映了火山的活跃程度。

2 2017年全球活动火山时间分布特征

2017年全球范围内，平均每月约有30座火山喷发的记录，最多达34座，最少为26座（图2）。年度内火山活动数量和活动强度有逐渐增多、增强的趋势，表明俯冲板块活动的加剧。可将活动的火山分为3种类型（图3）：①“冒泡”型，一般活动强度较弱，持续时间较短，警戒等级一般为I级或Ⅱ级，活动形式主要表现为火山喷气、火山震颤、火山地震等，通常不会造成地质灾害和影响周围居民的生活，典型代表有冰岛的卡特拉火山（Katla）、美国的希沙尔丁火山（Shishaldin）、印度尼西亚的赛梅鲁火山（Semeru）和俄罗斯的朱帕诺夫斯基火山（Zhupanovsky）等；②“持续”型，基本常年处于持续喷发状态，警戒等级一般为Ⅱ级或Ⅲ级，活动形式主要表现为火山喷气、火山地震、岩浆活动、火山碎屑流等，这类火山是全球火山监测的重点，对其进行的火山喷发监测和火山灾害预防都较为成熟，典型代表有美国夏威夷的基拉韦厄火山（Kilauea）（Fontijn等，2015；王佳龙，2018b）、印度尼西亚的锡纳朋火山（Sinabung）、秘鲁的萨班卡亚火山（Sabancaya）和俄罗斯勘察加半岛的希韦卢奇火山（Sheveluch），这几座火山几乎常年都在活动，火山区周围常引发地质灾害；③“剧烈”型，常年处于稳定状态，偶尔几周或几个月突然爆发，警戒等级一般为Ⅲ级或Ⅳ级，喷发的形式主要表现为喷发柱的突然上升，有时高达几千米至十几千米，随后火山碎屑物质喷出火口并造成坍塌，引发火山碎屑流，最后岩浆溢出形成熔岩流，典型代表有日本的新燃岳火山（Kirishimayama）、菲律宾的马荣火山（Mayon）以及印度尼西亚的阿贡火山（Agung）（王佳龙，2018a）。

图2 2017年度全球火山喷发数量

图3 活动火山分类

3 阿贡火山的喷发过程

全球的活火山主要分布在环太平洋地区，此区域也是著名的“火环”，而印度尼西亚岛链是西南太平洋最为活跃的1条“火链”（图4，图中火山名称下方数字，如2003/46/2代表的含义为2003年第1次监测到该火山的活动，目前共计活动了46周，2017年度从2月开始变活跃），其居民数量为世界之最，阿贡火山则是该“火链”中最为活跃的火山之一（Self等1996）。阿贡火山的喷发周期约50年，最近1次喷发于1963年（Zen等1964；Marinelli等，1968），距今56年。据印尼国家灾害管理局[3]（Badan Nasional Penanggulangan Bencana）的监测数据，自2017年8月开始，阿贡火山的地震明显增多，有复苏和再次喷发的前兆，9月份火山警戒等级升至Ⅲ级，10月份开始有喷发柱喷出，11月喷发柱急剧上升至4km并开始有岩浆活动，警戒等级也升至最高级Ⅳ级。印度尼西亚岛链位于澳大利亚板块与欧亚板块之间，由于澳大利亚板块向北的俯冲，在此区域形成了1套完整的沟-弧-盆体系。5—6级地震的震中主要位于岛链与巽他海沟之间，与岛链大致平行，距离岛链约100—150km，地震震源深度约100km；活火山主要位于岛链内部，岩浆来源大致在深度150—200km板块的交汇处。由此可见，由澳大利亚板块俯冲而下的矿物在岛链下方150—200km处发生脱水作用，脱水诱导下覆地震的部分熔融，产生低盐度、钙碱性岩浆；岩浆上升侵入上覆板块的岩石圈中。

图4 印度尼西亚岛链地震与火山分布

Fig. 4 Earthquakes and volcanoes in Indonesia island chain

4 阿贡火山喷发过程中的形变

在火山喷发过程中，前期地下岩浆活动，地表热异常，火山开始膨胀；当压力达到一定程度后，火山气体首先冲出，火山通道内部的压力减小，岩浆成分挥发，火山灰喷出，随后岩浆补给，岩浆爆炸喷出或溢出火口，后期岩浆冷却收缩，火口坍塌。目前，阿贡火山处在火山灰的喷发过程，未出现岩浆溢出。自2017年6月开始，阿贡火山有所膨胀，至2017年12月膨胀近15cm，随后开始逐渐收缩（图5）。膨胀的中心位于阿贡火山北侧附近，并未在火山的正下方，说明岩浆通道向北倾斜。

图5 阿贡火山喷发过程干涉条纹图及对应的地表形变

采用“二通”法对收集到的SAR复数影像进行干涉差分处理。在处理过程中，使用SRTM4 DEM数据消除地形相位影响，同时，为保证较高的配准精度，采用结合轨道和地形数据的图像配准技术，配准精度优于0.001个像元。由于巴厘岛大部分地区为山地，且气候温和多雨，岛上植被茂盛，易造成干涉像对的失相干，同时SAR影像在山区易出现叠掩、透视收缩、阴影等，也将加剧失相干的程度。失相干区或低相干区将导致相位解缠误差，为此，使用迭代自适应滤波算法进行降噪处理，并采用基于狄洛尼三角剖分的最小费用流算法（MCF）。首先，对高质量的相位区进行解缠获得可靠的参考相位模型，再利用参考相位实现对低相干区域的解缠，从而得到全局的最优结果，最后，经过地理编码获取了阿贡火山喷发过程的InSAR形变场。

5 阿贡火山喷发柱高度变化

阿贡火山自2017年9月开始活动以来，活动特征明显，浅源小震持续增多、热异常明显。2017年10月后，阿贡火山开始出现不同规模的爆炸，爆炸产生的火山灰柱上升至高空并随风向飘逸，爆炸的规模直接决定了喷发柱的高度。阿贡火山的整个喷发过程以及爆炸与喷发柱变化的关系，如图6所示。由图可见，阿贡火山在2017年10月份主要表现为火山地震和震颤，是地下岩浆通过管道向上运移的过程，此时山体表面也开始不断膨胀，该过程一直持续到11月27日，当日发生1次剧烈爆炸，造成火山剧烈喷发，喷发柱升至海拔4km，大量火山碎屑流顺火山翼冲下，当地居民被迫迁移，航班被迫取消；在随后的3个月内，不断有地下岩浆补给，经常性地发生爆炸并产生喷发柱；在2018年2月27日之后，岩浆补给减弱，火山警戒等级下降，爆炸规模减小，喷发柱高度降低，火山表面逐渐收缩；自2017年9月末至2018年3月末，阿贡火山历经了整个喷发过程，从开始的地震活动，到喷发时大规模爆炸，再到最后阶段的平静收缩，此过程对研究其它“剧烈”型火山具有借鉴意义。

图6 阿贡火山喷发过程喷发柱高度变化

图7 2015—2018年全球活动火山数量变化

6 结论

2017年全球火山喷发较往年强烈，年度内火山数量和强度均有增大的趋势（图7）。其中“剧烈”型火山造成的火山灾害较为严重，应对其加强监测和防灾。活火山主要分布在环太平洋地区，位于太平洋西南的印度尼西亚是2017年受灾最为严重的地区，据印尼国家灾害管理局的统计，阿贡火山的喷发造成10万余人无家可归。阿贡火山历时半年多的活动过程从膨胀到收缩，从平静到爆炸再到平静，使我们对火山喷发的过程有了更全面的了解。

洪汉净，于泳，郑秀珍等，2003. 全球火山活动分布特征. 地学前缘，10（S1）：11—16.

洪汉净，陈会仙，赵谊等，2009. 全球地震、火山分布及其变化特征. 地震地质，31（4）：573—583.

刘若新，樊祺诚，魏海泉等，1999. 中国活火山研究. 地质论评，45（S1）：3—15.

马宗晋，杜品仁，高祥林，2003. 全球构造研究的思考. 地学前沿，10（S1）：1—4.

王佳龙，张进奎，邹颖，2018a. 近期阿贡火山（巴厘岛，印尼）喷发警戒等级已升至最高级. 国际地震动态，（2）：2—8，36.

王佳龙，邹颖，张进奎等，2018b. 美国夏威夷基拉韦厄火山及其最新活动情况. 国际地震动态，（6）：2—14.

许建东，2011. 中国活动火山监测进展回顾. 矿物岩石地球化学通报，30（4）：390—392.

张传杰，李霓，龚丽文，2016. 2015年全球火山活动综述. 国际地震动态，（11）：4—11.

Fontijn K., Costa F., Sutawidjaja I., et al., 2015. A 5000-year record of multiple highly explosive mafic eruptions from Gunung Agung (Bali, Indonesia): implications for eruption frequency and volcanic hazards. Bulletin of Volcanology, 77(7): 59.

Marinelli G., Tazieff H., 1968. L'Ignimbrite et la caldera de Batur (Bali, Indonesie). Bulletin Volcanologique, 32(1): 89—120.

Press F., Siever R., 1982. Earth. 3rded. New York: W H Freeman and Company, 366—369.

Self S., King A. J., 1996. Petrology and sulfur and chlorine emissions of the 1963 eruption of Gunung Agung, Bali, Indonesia. Bulletin of Volcanology, 58(4): 263—285.

Siebert L., Simkin T., Kimberly P., 2010. Volcanoes of the world. Washington DC, Berkeley: Smithsonian Institution, University of California Press, 551—559.

Zen M. T., Hadikusumo D., 1964. Recent changes in the Anak-Krakatau volcano. Bulletin Volcanologique, 27(1): 259—268.

The Temporal and Spatial Distribution of Global Volcanism and the Deformation Characteristics Before and After Eruption of Mount Agung in 2017

Wang Jialong1, 2), Wang Jiaqing1), Zou Ying2), Sun Jiaxiang2)and Ma Qingzun1)

1) The First Monitoring Center, China Earthquake Administration, Tianjin 300180, China 2) Institute of Geology, China Earthquake Administration, Beijing 100029, China

In this article we summarize the activities of active volcanoes worldwide in 2017 (until December 31, 2017). There are a total of 81 volcanoes in the world with 1,058 eruptions. On average, 20 active volcanoes are recorded weekly. According to the potential eruption risk and the intensity of volcanic activity, the above volcanoes were classified. The volcanic activity mainly reflects the tectonic activities of the Earth's surface. Among them, the post-arc volcano of the high-angle subduction zone is relatively high, while the volcanic activity of the low-angle subduction zone, the stretched rift, and the plateau dominated by strike-slip is relatively low. The Indonesian island chain, where volcanic activity is high, is one of the regions most affected by the volcanic hazard. Hundreds of thousands of people have to move to shelters. It is expected that the global volcanic activity will further intensify in the next year, and the island chain is still threatened by volcanic hazard. The Agung volcano on Indonesian island chain has been representative of the entire eruption process since its inception in September 2017. Monitoring of the eruption process of the Agung volcano can provide reference for typical volcanic eruptions in the world.

Active volcano; Volcanic activity; Tectonic activity; Volcanic hazard; Agung volcano

10.11899/zzfy20190215

中国地震局地震监测预报司“全球典型火山事件跟踪”项目

2018-07-04

王佳龙，男，生于1985年。助理工程师。主要从事构造地质与火山地质研究。E-mail：85102605@qq.com

王佳龙，王家庆，邹颖，孙嘉祥，马庆尊，2019．2017年全球活动火山时空分布以及阿贡火山喷发前后的形变特征．震灾防御技术，14（1）：423—430．

[1] http://volcano.si.edu/reports_weekly.cfm

[2] https://volcanoes.usgs.gov/vhp/about_alerts.html

[3] https://www.bnpb.go.id