A large hydrothermal reservoir beneath Taal Volcano (Philippines) revealed by magnetotelluric observations and its implications to the volcanic activity

Taal Volcano is one of the most active volcanoes in the Philippines. The magnetotelluric 3D forward analyses indicate the existence of a large high resistivity anomaly (∼100 Ω·m) with a volume of at least 3 km × 3 km × 3 km, which is capped by a conductive layer (∼10 Ω·m), beneath the Main Crater. This high resistivity anomaly is hypothesized to be a large hydrothermal reservoir, consisting of the aggregate of interconnected cracks in rigid and dense host rocks, which are filled with hydrothermal fluids coming from a magma batch below the reservoir. The hydrothermal fluids are considered partly in gas phase and liquid phase. The presence of such a large hydrothermal reservoir and the stagnant magma below may have influences on the volcano’s activity. Two possibilities are presented. First, the 30 January 1911 explosion event was a magmatic hydrothermal eruption rather than a base-surge associated with a phreato-magmatic eruption. Second, the earlier proposed four eruption series may be better interpreted by two cycles, each consisting of series of summit and flank eruptions.     

三维电磁扩散场数值模拟及磁化效应的影响

采用矢量有限元法实现了三维电磁扩散场数值模拟,并成功将其应用在大地电磁的正演研究中.为灵活精确地拟合起伏地形和地下不规则构造,采用由不规则四面体单元组成的非结构化网格,可根据模型设计的需要调整网格的大小.引入了基于二次场理论,将解析的一次场从总场中扣除,直接计算二次场,使得误差仅局限于相对较小的二次场,以提高总场计算精度.常规的节点有限元法不满足电性分界面上法向电场不连续和无源区单元内电流密度无散,违反麦克斯韦方程组.为克服节点有限元法的弊端,使用矢量有限元法求解基于二次电场的偏微分方程.另外,在算法设计中,考虑了磁导率参数的变化,可以模拟磁导率不均匀的模型.通过与COMMEMI模型已发表的结果对比,证明了本文算法的正确性和精确性.为突显非结构网格优势,计算了椭球异常体模型和任意地形模型的MT响应,并详细讨论了地形和磁化效应对三维数值模拟结果的影响.     

Global and domain-decomposed mixed methods for the solution of Maxwell's equations with application to magnetotellurics

A collection of global and domain decomposition mixed finite element schemes for the approximate solution of the harmonic Maxwell's equations on a bounded domain with absorbing boundary conditions at the artificial boundaries are presented. The numerical procedures allow us to solve efficiently the direct problem in magnetotellurics consisting of determining the electromagnetic scattered field in a two–dimensional earth model of arbitrary conductivity properties. Convergence results for the numerical procedures are derived. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 407–437, 1998