2022年1月15日世界标准时间4∶14，位于南太平洋汤加俯冲带火山链区域的Hunga Tonga-Hunga Ha'apai（HTHH，175.39°W，20.55°S）火山爆发并激发全球性海啸。聚焦火山喷发和海啸波的观测方法，介绍了海底传感器、大气波变化卫星图、电离层扰动对火山喷发和海啸的最新监测结果。结合该区域的地质构造与火山喷发历史资料，综述了2022年汤加HTHH火山喷发引起的海啸特征、最新观测手段，提出了对俯冲带火山海啸等进行详细的海洋地球物理调查、基于观测综合量化研究海底火山大规模喷发对全球气候的潜在影响等建议。
俯冲带是火山喷发、地震、大型滑坡等重大自然灾害的频发区域。这些海底构造事件的发生可能引发超大海啸，给临近沿岸地区造成大量人员伤亡与巨大经济损失。2022年1月14日，位于南太平洋汤加俯冲带弧后区域的Hunga Tonga-Hunga Ha'apai（以下简称HTHH）海底火山喷发（图1），释放里氏震级5.8级的能量，并引发了全球性海啸，这次火山喷发及其产生的罕见海啸在10h内迅速传播到了环太平洋沿岸地区，导致太平洋沿岸国家相继发布海啸预警；海啸信号随后到达了大西洋、印度洋、加勒比海和地中海等区域。此外，卫星观测数据显示此次火山喷发羽流高达58km，达到了大气中间层，造成了电离层扰动，并从火山喷发区域在大气层以波动的形式向全球范围传播，波动气压的变化被全球气压计网路成功记录（图1）。
图1 环太平洋海底压力计（蓝色、绿色倒三角表示）与气压计（红圈）台站分布图（a）与汤加火山区域地质图（b）
聚焦此次罕见火山喷发及其引发海啸的最新观测结果，综合整理与分析2022年HTHH火山喷发前地表的形变特征，喷发后引发的海啸以及大气波、电离层等多种新观测结果，对俯冲带火山海啸等海洋地质灾害的研究提出展望与建议。
汤加俯冲带区域构造及HTHH火山喷发前的地表形变特征
01 - 构造背景与历史活动
在西南太平洋区域，快速移动的太平洋板块向澳大利亚板块俯冲汇聚，形成了世界上仅次于马里亚纳海沟的第二深汤加-克马德克海沟（图1）。该俯冲带的汇聚速率约为20cm/a，具有非常活跃的弧后扩张和板片后撤等特点。在过去约45个百万年内，汤加俯冲带经历了多期次的弧后扩张，并在弧后区域形成汤加火山弧链（图1（b）），该火山弧跨度很长，从斐济的俯冲带一直延伸到新西兰东北边，其所覆盖的区域构造运动活跃，历史上火山、地震、海啸等地质灾害频发。
此次最新喷发的HTHH火山口位于汤加首都努库阿洛法（Nuku'alofa）以北约68km处（图1（b））。该火山口在2014—2015年连续喷发后连接了2座无人居住的Hunga-Ha'apai和Hunga-Tonga小岛。2022年喷发前的HTHH火山最高点海拔约100m，火山口处于海水下方，直径约5km，整个火山宽约20km。至20世纪有历史记录以来，HTHH火山共喷发6次，分别为1912年、1937年、1988年（持续3d）、2009年（持续1周）和2014—2015年（持续5周）。其中2009年的火山喷发的火山爆发指数为2，且被追踪到了火山的喷发阶段；Garvin等通过卫星数据分析2014—2015年火山喷发形成的火山锥侵蚀速率约为0.00256/a。以上历史数据显示HTHH火山的喷发间隔约为20~50a，最新的2022年喷发时间间隔仅为7a，表明该区域构造与火山活动性明显增加。
02 - 2022年HTHH火山喷发前的地表形变特征
对于HTHH火山喷发前的地表形变特征，胡羽丰等基于Sentinel-1卫星数据进行了分析，获取了该火山在喷发前2019-07-06至2021-12-10期间的卫星雷达视线（LOS）向累积形变时间序列（图2（a）~（h））。其中红色为负值，代表沿LOS向远离卫星；蓝色为正值，表示沿LOS向靠近卫星。该结果显示，在火山喷发前的观测时间内，火山体有明显的下沉现象且形变场呈类似于火山地表的半圆环形状，并在火山口向山底方向位移速率逐渐减小。进一步分析火山口形变区2个参考特征点和的时间序列显示（图2（i）和（j）），点在观测期间的累积位移下沉超过6cm，点的累积位移下沉约5cm，其下沉速率也呈阶段性变化，在2020年6月和2020年8月前基本稳定，随后、点分别以约3.1cm/a和2.2cm/a的速度下沉，定量揭示了火山喷发前的地表轻微形变的时空演化特征。
图2 HTHH火山喷发前地表形变记录
2022年HTHH火山喷发形变特征
01 - 引发全球性海啸
汤加HTHH火山喷发后引发的海啸在几到数十个小时内迅速扩张至全球。全球沿海测潮站网络，在太平洋、大西洋、印度洋以及在加勒比海和地中海等区域广泛观测到了明显的海啸波。Kubota等对环太平洋海底压力计和气压计记录的海啸信号进行了分析（图3），发现此次海啸的实际到时（图3（a）中的黑线）比理论上常规海底火山喷发引发的海啸到达时间（图3（a）中的绿线）早了2个多小时，传播时间更为迅速；汤加的最新报告显示当地的海啸可能至少达到15m，海啸波传到日本时高达2.7m。值得注意的是，太平洋周边沿岸海啸波的特征与世界其他地区明显不同。这些沿岸台站记录显示，在最大的海啸波（高达3m）达到之前先监测到了振幅较小的先导海啸波；在太平洋之外的海域则主要监测到了这些最初的小波幅海啸。
图3 HTHH火山喷发后环太平洋海底压力计信号记录（a）与气压记录（b）
汤加HTHH火山喷发引发的海啸与常见的俯冲带大地震引发的海啸相比有显著的不同。一般情况下地震引发的海啸会沿着长达数十或上百公里的断层线将海水从海底向上排开，而火山爆发引发的海啸更像是一个点源向外扩散。目前用于研究的海啸预报模型基本是基于地震触发进行模拟，不太符合火山海啸的情形。此外，这次喷发引起的海底压力变化持续时间更长，海啸的产生机制也更为复杂。最先到达沿岸的海啸波波幅较小且与火山爆发产生快速移动大气压力波的传播速度耦合（约300m/s），被认为是由大气波传播过程中与海水相互作用的结果。
02 - 激发大气波
从理论上讲，火山爆发和核试验等大型爆炸将产生不同长度和频率的大气波。在水平方向上的短波长尺度上，可以引发兰姆波（Lamb wave）、声波和重力波。此外，火山可以在最初喷发后通过上升气流与加热羽流对流形成持续的波源。HTHH火山上空羽流持续释放潜热超过12h，是此次大气波最重要的单个重力波源，激发了横跨整个太平洋的圆形波面，之前从未观测到单一震源产生如此量级的全球范围的响应（图4）。地球同步环境卫星GOES观测显示，在世界标准时间（UTC）04∶15左右，开始可以看到火山爆发，形成的羽流在30min内上升至30km高，覆盖200km的宽度。在羽流上升20~30min后，在10min分辨率的近红外地球静止图像中可以看到明显的大气波。Wright等通过地面压力数据的反向投影分析，判断触发源发生在UTC04∶28±0∶02，前波面以（318.2±6）m/s的近地面相速度在地表层传播，并根据信号的高相位速度、大振幅和非色散特性分析，将其识别为兰姆波。理论上兰姆波的能量密度随高度呈指数衰减，观测显示其以（308±5）至 （319±4）m/s的相速度在平流层中传播。兰姆波在全球范围内传播，在火山喷发后18.1h（±7.5min）经过阿尔及利亚对跖点，接下来的几天里绕地球传播了至少3次。基于大气红外测深仪（AIRS）、交叉航迹红外探测仪（CrIS）、红外大气探测干涉仪（IA⁃SI）和极轨热红外探测仪（IR）数据，初始重力波在兰姆波之后在全球范围内同样被广泛识别。重力波于1月16日UTC00∶30—02∶30（喷发后20~22h内）首次经过HTHH火山的对跖点，以（238±3）至（269±3）m/s的相速度在平流层中传播。综合观测分析显示，初始的兰姆波和重力波之间的传播距离差随着时间的推移而增大。
图4 HTHH火山爆发引发的大气波传播（地球同步环境卫星GOES观测）
03 - 引起电离层扰动
地球表面的突然运动引起的大气波会通过大气传播到上层大气，最终引起移动电离层的扰动（TIDs）。这种变化过去已经在地震、海啸、火山爆发、地下核爆炸等情景下被监测到。HTHH火山喷发后，火山上空电离层的垂直总电子含量（TEC）明显减少，全球导航卫星系统（GNSS）接收网络在日本上空监测到了电离层扰动的移动轨迹。Saito在TEC的扰动中观察到2种不同特征的电离层扰动信号。第一个扰动信号在火山喷发3h后抵达日本，距离火山口约7800km，振幅约为±0.5TECU。第二个扰动信号在火山喷发7h后抵达日本，振幅约为±1.0TECU。
除了日本上空，全球的其他区域也监测到了HTHH火山喷发引起的电离层扰动信号。Themens等使用了全球4735个全球导航卫星系统接收器的测量数据，第一时间跟踪了此次HTHH火山喷发造成相关的电离层扰动的移动信号（图5）。通过选取全球范围内不同区域的电离层异常信号进行分析（新西兰、澳大利亚、夏威夷、日本、北美东部、南非和北欧，图5（b）~（h）），发现2个不同的大尺度移动电离层扰动和几个随后的中尺度移动电离层扰动，这2个扰动信号都从喷发点向外径向传播。在距离喷发中心3000km的范围内，最初观测到波长为1600km的大尺度移动电离层扰动信号分别以约950m/s和约555m/s的速度传播，随后显著减缓至约600和约390m/s。该信号为全球电离层响应的主要信号，与火山喷发相关的大气表面压力扰动一致。中尺度移动电离层扰动信号在火山喷发后6h内观测到，其速度为200~400m/s。
图5 HTHH火山爆发在全球范围内引发的电离层效应观测记录
结论
总结了2022年汤加HTHH火山爆发以来关于喷发和海啸的最新观测与研究。HTHH火山海啸在世界各地引发的不同寻常地球物理场的变化，反映了海底火山爆发时，整个喷发过程与大气层、岩石圈和海洋相互作用的复杂性，及其可能造成的全球范围的海洋地质灾害影响。从这次事件获得的启示和对将来研究建议如下。
1）此次火山喷发引起的海啸特征复杂，尤其是火山体附近的地质结构变化及其对海啸的贡献尚不清楚，亟需详细的海洋地球物理调查来揭开这个谜团。现有的模型不足以提供准确的海啸模拟，因此需要建立一种新的非地震触发海啸模型，并纳入到全球海啸预警系统，以便及时为具有相似俯冲带地质构造特征的沿岸城市和居民提供相关紧急预警、疏散和救援等措施。
2）此次火山喷发向大气中释放了大量的能量和不同成分的气体，在全球范围内激发了重力波、大气波，造成电离层扰动等，迫切需要基于这些最新的观测综合量化研究海底火山大规模喷发对全球气候的潜在影响，且跨学科合作探讨海洋地质灾害与大气响应之间的关系，地球不同圈层之间的相互作用等。
3）利用新观测手段加强海底火山型海啸的观测，建立和完善全球俯冲带沿岸海啸预警与气候监测系统。
4）巨大型海洋地质灾害是全球范围内的，建议增强国际合作与救援，共同打造一个持续性宜居的家园。
本文作者：郑婷婷、邱强、林间
作者简介：郑婷婷，中国科学院边缘海与大洋地质重点实验室，中国科学院南海海洋研究所，南海生态环境工程创新研究院，南方海洋学与工程广东省实验室（广州），中国-巴基斯坦地球科学研究中心，中国科学院-巴基斯坦高等教育委员会，助理研究员，研究方向为海洋地质与地球物理；林间（通信作者），中国科学院边缘海与大洋地质重点实验室，中国科学院南海海洋研究所，南海生态环境工程创新研究院，南方海洋学与工程广东省实验室（广州），中国-巴基斯坦地球科学研究中心，中国科学院-巴基斯坦高等教育委员会，南方科技大学海洋科学与工程系，伍兹霍尔海洋研究所，教授，研究方向为海洋地质与地球物理、地球动力学。
原文发表于《科技导报》2023年第2期，欢迎订阅查看。