English   español  
Por favor, use este identificador para citar o enlazar a este item: http://hdl.handle.net/10261/133107
Compartir / Impacto:
Estadísticas
Add this article to your Mendeley library MendeleyBASE
Visualizar otros formatos: MARC | Dublin Core | RDF | ORE | MODS | METS | DIDL
Título

FEM-based linear inverse modeling using a 3D source array to image magma chambers with free geometry. Application to InSAR data from Rabaul Caldera (PNG).

AutorRonchin, Erika ; Masterlark, T.; Dawson, J.; Saunders, S.; Martí Molist, Joan
Fecha de publicación12-abr-2015
EditorEuropean Geosciences Union
CitaciónEuropean Geosciences Union General Assembly 2015 Vienna, Austria 12 – 17 April 2015 EGU2015-728, 2015
ResumenIn this study, we present a method to fully integrate a family of finite element models (FEMs) into the regularized linear inversion of InSAR data collected at Rabaul caldera (PNG) between February 2007 and December 2010. During this period the caldera experienced a long-term steady subsidence that characterized surface movement both inside the caldera and outside, on its western side. The inversion is based on an array of FEM sources in the sense that the Green’s function matrix is a library of forward numerical displacement solutions generated by the sources of an array common to all FEMs. Each entry of the library is the LOS surface displacement generated by injecting a unity mass of fluid, of known density and bulk modulus, into a different source cavity of the array for each FEM. By using FEMs, we are taking advantage of their capability of including topography and heterogeneous distribution of elastic material properties. All FEMs of the family share the same mesh in which only one source is activated at the time by removing the corresponding elements and applying the unity fluid flux. The domain therefore only needs to be discretized once. This precludes remeshing for each activated source, thus reducing computational requirements, often a downside of FEM-based inversions. Without imposing an a-priori source, the method allows us to identify, from a least-squares standpoint, a complex distribution of fluid flux (or change in pressure) with a 3D free geometry within the source array, as dictated by the data. The results of applying the proposed inversion to Rabaul InSAR data show a shallow magmatic system under the caldera made of two interconnected lobes located at the two opposite sides of the caldera. These lobes could be consistent with feeding reservoirs of the ongoing Tavuvur volcano eruption of andesitic products, on the eastern side, and of the past Vulcan volcano eruptions of more evolved materials, on the western side. The interconnection and spatial distribution of sources find correspondence in the petrography of the volcanic products described in literature and in the dynamics of the single and twin eruptions that characterize the caldera. As many other volcanoes, Rabaul caldera is an active and dangerous volcanic system whose dynamics still need to be understood to effectively predict the behavior of future eruptions. The good results obtained from the application of the method to Rabaul caldera show that the proposed linear inversion based on the FEM array of sources is suitable to generate models of magmatic systems. The method can image in space and time the complex free geometry of the source that generates the deformation, widening our understanding of deformational sources and their dynamics. This takes source modeling a step towards more realistic source models.
DescripciónThis article is subject to a CC Attribution 3.0 License.
URIhttp://hdl.handle.net/10261/133107
Aparece en las colecciones: (ICTJA) Comunicaciones congresos
Ficheros en este ítem:
Fichero Descripción Tamaño Formato  
Ronchin_EGU2015-728.pdf41,34 kBAdobe PDFVista previa
Visualizar/Abrir
Mostrar el registro completo
 


NOTA: Los ítems de Digital.CSIC están protegidos por copyright, con todos los derechos reservados, a menos que se indique lo contrario.