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