3-D Modeling of Irregular Volcanic Sources Using Sparsity-Promoting Inversions of Geodetic Data and Boundary Element Method

Research output: Contribution to journalArticle

Abstract

Geodetic observations of surface deformation associated with volcanic activities can be used to constrain volcanic source parameters and their kinematics. Simple analytical models, such as point and spherical sources, are widely used to model deformation data. The inherent nature of oversimplified model geometries makes them unable to explain fine details of surface deformation. Current nonparametric, geometry-free inversion approaches resolve the distributed volume change, assuming it varies smoothly in space, which may detect artificial volume change outside magmatic source regions. To obtain a physically meaningful representation of an irregular volcanic source, we devise a new sparsity-promoting modeling scheme assuming active magma bodies are well-localized melt accumulations, namely, outliers in the background crust. First, surface deformation data are inverted using a hybrid L1- and L2-norm regularization scheme to solve for sparse volume change distributions. Next, a boundary element method is implemented to solve for the displacement discontinuity distribution of the reservoir, which satisfies a uniform pressure boundary condition. The inversion approach is thoroughly validated using benchmark and synthetic tests, of which the results show that source dimension, depth, and shape can be recovered appropriately. We apply this modeling scheme to deformation observed at Kilauea summit for periods of uplift and subsidence leading to and following the 2007 Father's Day event. We find that the magmatic source geometries for these periods are statistically distinct, which may be an indicator that magma is released from isolated compartments due to large differential pressure leading to the rift intrusion.

Original languageEnglish (US)
JournalJournal of Geophysical Research: Solid Earth
DOIs
StateAccepted/In press - Jan 1 2017

Fingerprint

geodesy
datum (elevation)
boundary element method
geodetic datum
Boundary element method
volcanology
inversions
volume change
modeling
geometry
Geometry
magma
volcanic activity
methodology
subsidence
source parameters
Subsidence
fathers
differential pressure
outlier

Keywords

  • Geodetic inversion
  • InSAR deformation
  • Irregular source geometry
  • Kilauea volcano
  • Magma chamber
  • Sparsity regularization

ASJC Scopus subject areas

  • Geophysics
  • Oceanography
  • Forestry
  • Aquatic Science
  • Ecology
  • Condensed Matter Physics
  • Water Science and Technology
  • Soil Science
  • Geochemistry and Petrology
  • Earth-Surface Processes
  • Physical and Theoretical Chemistry
  • Polymers and Plastics
  • Atmospheric Science
  • Earth and Planetary Sciences (miscellaneous)
  • Space and Planetary Science
  • Materials Chemistry
  • Palaeontology

Cite this

@article{dae86b326c754934a414fd7e292bacf1,
title = "3-D Modeling of Irregular Volcanic Sources Using Sparsity-Promoting Inversions of Geodetic Data and Boundary Element Method",
abstract = "Geodetic observations of surface deformation associated with volcanic activities can be used to constrain volcanic source parameters and their kinematics. Simple analytical models, such as point and spherical sources, are widely used to model deformation data. The inherent nature of oversimplified model geometries makes them unable to explain fine details of surface deformation. Current nonparametric, geometry-free inversion approaches resolve the distributed volume change, assuming it varies smoothly in space, which may detect artificial volume change outside magmatic source regions. To obtain a physically meaningful representation of an irregular volcanic source, we devise a new sparsity-promoting modeling scheme assuming active magma bodies are well-localized melt accumulations, namely, outliers in the background crust. First, surface deformation data are inverted using a hybrid L1- and L2-norm regularization scheme to solve for sparse volume change distributions. Next, a boundary element method is implemented to solve for the displacement discontinuity distribution of the reservoir, which satisfies a uniform pressure boundary condition. The inversion approach is thoroughly validated using benchmark and synthetic tests, of which the results show that source dimension, depth, and shape can be recovered appropriately. We apply this modeling scheme to deformation observed at Kilauea summit for periods of uplift and subsidence leading to and following the 2007 Father's Day event. We find that the magmatic source geometries for these periods are statistically distinct, which may be an indicator that magma is released from isolated compartments due to large differential pressure leading to the rift intrusion.",
keywords = "Geodetic inversion, InSAR deformation, Irregular source geometry, Kilauea volcano, Magma chamber, Sparsity regularization",
author = "Guang Zhai and Manoochehr Shirzaei",
year = "2017",
month = "1",
day = "1",
doi = "10.1002/2017JB014991",
language = "English (US)",
journal = "Journal of Geophysical Research: Atmospheres",
issn = "2169-897X",
publisher = "Wiley-Blackwell",

}

TY - JOUR

T1 - 3-D Modeling of Irregular Volcanic Sources Using Sparsity-Promoting Inversions of Geodetic Data and Boundary Element Method

AU - Zhai, Guang

AU - Shirzaei, Manoochehr

PY - 2017/1/1

Y1 - 2017/1/1

N2 - Geodetic observations of surface deformation associated with volcanic activities can be used to constrain volcanic source parameters and their kinematics. Simple analytical models, such as point and spherical sources, are widely used to model deformation data. The inherent nature of oversimplified model geometries makes them unable to explain fine details of surface deformation. Current nonparametric, geometry-free inversion approaches resolve the distributed volume change, assuming it varies smoothly in space, which may detect artificial volume change outside magmatic source regions. To obtain a physically meaningful representation of an irregular volcanic source, we devise a new sparsity-promoting modeling scheme assuming active magma bodies are well-localized melt accumulations, namely, outliers in the background crust. First, surface deformation data are inverted using a hybrid L1- and L2-norm regularization scheme to solve for sparse volume change distributions. Next, a boundary element method is implemented to solve for the displacement discontinuity distribution of the reservoir, which satisfies a uniform pressure boundary condition. The inversion approach is thoroughly validated using benchmark and synthetic tests, of which the results show that source dimension, depth, and shape can be recovered appropriately. We apply this modeling scheme to deformation observed at Kilauea summit for periods of uplift and subsidence leading to and following the 2007 Father's Day event. We find that the magmatic source geometries for these periods are statistically distinct, which may be an indicator that magma is released from isolated compartments due to large differential pressure leading to the rift intrusion.

AB - Geodetic observations of surface deformation associated with volcanic activities can be used to constrain volcanic source parameters and their kinematics. Simple analytical models, such as point and spherical sources, are widely used to model deformation data. The inherent nature of oversimplified model geometries makes them unable to explain fine details of surface deformation. Current nonparametric, geometry-free inversion approaches resolve the distributed volume change, assuming it varies smoothly in space, which may detect artificial volume change outside magmatic source regions. To obtain a physically meaningful representation of an irregular volcanic source, we devise a new sparsity-promoting modeling scheme assuming active magma bodies are well-localized melt accumulations, namely, outliers in the background crust. First, surface deformation data are inverted using a hybrid L1- and L2-norm regularization scheme to solve for sparse volume change distributions. Next, a boundary element method is implemented to solve for the displacement discontinuity distribution of the reservoir, which satisfies a uniform pressure boundary condition. The inversion approach is thoroughly validated using benchmark and synthetic tests, of which the results show that source dimension, depth, and shape can be recovered appropriately. We apply this modeling scheme to deformation observed at Kilauea summit for periods of uplift and subsidence leading to and following the 2007 Father's Day event. We find that the magmatic source geometries for these periods are statistically distinct, which may be an indicator that magma is released from isolated compartments due to large differential pressure leading to the rift intrusion.

KW - Geodetic inversion

KW - InSAR deformation

KW - Irregular source geometry

KW - Kilauea volcano

KW - Magma chamber

KW - Sparsity regularization

UR - http://www.scopus.com/inward/record.url?scp=85038127664&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85038127664&partnerID=8YFLogxK

U2 - 10.1002/2017JB014991

DO - 10.1002/2017JB014991

M3 - Article

JO - Journal of Geophysical Research: Atmospheres

JF - Journal of Geophysical Research: Atmospheres

SN - 2169-897X

ER -