利用SBAS—InSAR技术研究印尼Agung火山形变特征与岩浆房参数.docVIP

利用SBAS—InSAR技术研究印尼Agung火山形变特征与岩浆房参数 　　摘要：基于ALOS PALSAR影像，利用小基线集合成孔径雷达干涉测量技术，提取了位于印度尼西亚巴厘岛的Agung火山2007～2009年的地表形变时间序列，并基于Mogi点源模型和竖直椭球体模型反演了岩浆房参数。结果表明：Agung火山地区大气延迟相位干扰较严重，Agung火山在2007～2009年发生了较明显的隆升形变，且与时间呈正相关。竖直椭球体模型能够更好地拟合InSAR形变场，岩浆房位于火山体下方约5 km处。SBASInSAR结果表明，应加强跟踪监测Agung火山的潜在喷发危险性。 　　关键词：SBAS-InSAR；形变时间序列；Agung火山；岩浆房参数；大气延迟相位 　　中图分类号：P315727文献标识码：A文章编号：1000-0666（2013）03-0313-06 　　0引言 　　Agung火山位于印度尼西亚的巴厘岛东部，是该岛的最高峰，最高处海拔3 142 m，构造上位于爪哇海沟以北约300 km处，属于印度—澳大利亚板块向东南亚大陆下方俯冲的弧后火山。Agung火山体呈锥状，顶部存有520 m×375 m的漏斗状破火山口，无植被生长，火山体的侧面陡峭，植被较茂盛（Dilmy，1965）。历史上，Agung火山曾多次喷发，在过去的300年里，有记录的喷发事件分别发生在1808、1821（待考证）、1843、1963～1964年（Seach，2012）。其中，1963～1964年的布里尼式喷发夺去了至少1 000人的生命，爆炸性喷发所产生的巨大的火山灰云随后笼罩整个地球，造成了巨大的灾难（Frederic，1970）。喷发历史表明，Agung火山属于具有巨大喷发危险的当代活动火山，一旦喷发，将会给周围居民乃至全球气候带来巨大灾难。然而，由于当地社会经济条件落后等方面的原因，Agung火山目前缺乏地表形变监测资料。 　　近年发展起来的InSAR技术因其自身的优势已在火山区地表形变监测方面得到了广泛应用（Massonnet et al，1995；Hooper et al，2004，2007；Peltier et al，2010）。常规InSAR技术容易受到时间、空间失相干的影响，Berardino等（2002）提出了小基线算法（small baseline subsets，简称SBAS），该算法仅选取短基线干涉对进行时间序列形变求解，削弱了空间失相干的影响，同时大大降低了地形因素引入的误差和大气延迟相位。Jung等（2008）对小基线算法在4个方面进行了精化，提高了时间序列形变结果的可靠性：（1）选取高质量的干涉图（不含相位解缠误差，且相干性好）估计外部DEM误差和线性形变相位；（2）对形变时间序列和误差进行迭代计算；（3）利用有限差分平滑方法削弱时间域的噪声；（4）对参考点本身存在的误差进行了校正。目前小基线算法已在火山区缓慢地表形变监测方面表现出了极大的潜力和优势（Lundgren et al，2004；Pepe et al，2008；Casu et al，2008；季灵运等，2011）。 　　1InSAR数据处理与结果 　　Agung火山位于热带地区，植被覆盖茂盛（图1），本文选用了对植被具有较强穿透能力的日本ALOS PALSAR影像（波长236 cm），获取时间从2007年7月至2009年1月，选取的影像提供了15 a的地表形变观测资料。 　　11InSAR数据处理 　　常规的DInSAR技术容易受到大气延迟相位的干扰。Zebker等（1997）的研究表明，云的相对湿度在空间或时间域内变化20%将会导致10～14 cm的形变误差。对于Agung火山区，大气相位延迟在某些干涉图中表现的比较明显。图2列出了3幅干涉图（已解缠），可以看出图2a和图2b中的颜色分布相反，若差分干涉相位不包含误差，则两者反映了相反的形变态势，即在图2a干涉对的时间间隔内，火山区域发生了明显的下沉，最大下沉达75 cm；而在图2b干涉对的时间间隔内，火山区域则发生了明显的隆升，最大隆升达38 cm；图2c却未反映出较大的相位差信息。2007年7～10月火山并无喷发，如此短的时间内（92 d）形变态势发生转换，不符合火山区形变的物理解释。图2a，b中的相位分布与地形具有较强的相关性，即靠近火山体顶部，相位的绝对值较大，顶部和侧面的相位差达到3～4弧度，符合对流层水汽垂直分层效应对雷达信号延迟干扰的特征。由此可以得出，图2a，b中的较大相位差均由2007年8月24日这景影像携带的大气延迟相位所导致。1地震研究136卷第3期季灵运等：利用SBASInSAR技术研究印尼Agung火山的形变特征与岩浆房参数图1Agung火山周围SRTM D