An official website of the United States government
Abstract Over the last 7 years, geodetic data have detected periods of uplift and subsidence of the active volcano Ol Doinyo Lengai in Tanzania. Although numerous eruptions of the volcano have occurred historically, a systematic investigation of transient deformation using continuous Global Navigation Satellite System (GNSS) data has not been undertaken. We use the Targeted Projection Operator (TPO) to assess 7 years of continuous GNSS data from the TZVOLCANO network for transient signals and find rapid uplift spanning March 2022–December 2022 and then steady‐state uplift through August 2023. We conduct a nonlinear inversion of the GNSS velocities associated with the transient signal using dMODELS and find consistency with an inflating spheroidal source located 2.3 ± 0.6 km beneath the crater. Prior to March 2022, geodetic data indicated quiescence just below Ol Doinyo Lengai, thus detecting transient deformation with TPO allows for tracking changes in the magmatic system over time in the Natron Rift.
more »
« less
Versatile Modeling Of Deformation (VMOD) Inversion Framework: Application to 20 Years of Observations at Westdahl Volcano and Fisher Caldera, Alaska, US
Abstract We developed an open source, extensible Python‐based framework, that we call the Versatile Modeling Of Deformation (VMOD), for forward and inverse modeling of crustal deformation sources. VMOD abstracts from specific source model implementations, data types and inversion methods. We implement the most common geodetic source models which can be combined to model and analyze multi‐source deformation. VMOD supports Global Navigation Satellite System (GNSS), InSAR, electronic distance measurement, Leveling and tilt data. To infer source characteristics from observations, VMOD implements non‐linear least squares and Markov Chain Monte‐Carlo Bayesian inversions, including joint inversions using different sources of data. VMOD's structure allows for easy integration of new geodetic models, data types, and inversion strategies. We benchmark the forward models against other published results and the inversion approaches against other implementations. We apply VMOD to analyze deformation at Unimak Island, Alaska, observed with continuous and campaign GNSS, and ascending and descending InSAR time series generated from Sentinel‐1 satellite radar acquisitions. These data show an inflation pattern at Westdahl volcano and subsidence at Fisher Caldera. We use VMOD to test a range of source models by jointly inverting the GNSS and InSAR data sets. Our final model simultaneously constrains the parameters of two sources. Our results reveal a depressurizing spheroid under Fisher Caldera ∼4–6 km deep, contracting at a rate of ∼2–3 Mm3/yr, and a pressurizing spherical source underneath Westdahl volcano ∼6–8 km deep, inflating at ∼5 Mm3/yr. This and past applications of VMOD to volcanic unrest benefit from an extensible framework which supports jointly inversions of data sets for parameters of easily composable multi‐source models.
more »
« less
- PAR ID:
- 10499196
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 25
- Issue:
- 4
- ISSN:
- 1525-2027
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
{"Abstract":["Supplementary code and model files for the manuscript entitled "Elucidating the Magma Plumbing System of Ol Doinyo Lengai (Natron Rift, Tanzania) Using Satellite Geodesy and Numerical Modeling". OlDoinyoLengai_code_and_models.zip contains all necessary Matlab code, functions, input and output files for the GNSS, InSAR, and joint inversions presented in our manuscript necessary to reproduce the results. dMODELS is an open source code developed by the United States Geological Survey. The originally published program is available here: https://pubs.usgs.gov/tm/13/b1/ and the revised software archived here will also be available through the USGS website code.usgs.gov/vsc/publications/OlDoinyoLengai or by contacting Maurizio Battaglia. With this manuscript we are providing an update to dMODELS that includes improved graphics and joint inversion capabilities for both InSAR and GNSS data. <\/p>"],"Other":["This work was funded by the National Science Foundation (NSF) grant number EAR-1943681, Virginia Tech, Korean Institute of Geosciences and Minerals (KIGAM), and Ardhi University. Funding for this work also came from USAID via the Volcano Disaster Assistance Program and from the U.S. Geological Survey (USGS) Volcano Hazards Program.This material is based on services provided by the GAGE Facility, operated by UNAVCO, Inc., with support from the National Science Foundation, the National Aeronautics and Space Administration, and the U.S. Geological Survey under NSF Cooperative Agreement EAR-1724794. We acknowledge and thank Alaska Satellite Facility for making InSAR data freely available and TZVOLCANO GNSS data sets available through the UNAVCO data archive."]}more » « less
-
Abstract Uturuncu volcano in southern Bolivia last erupted around 250 ka but is exhibiting signs of recent activity, including over 50 yr of surface uplift, elevated seismic activity, and fumarolic activity. We studied the spatial and temporal scales of surface deformation from 1992 to 2021 to better understand subsurface activity. We tracked Uturuncu’s recent deformation using interferometric synthetic aperture radar (InSAR) data and the global navigation satellite system (GNSS) station UTUR, located near Uturuncu’s summit. We observed a spatially coherent signal of uplift from 2014 to 2021 from Sentinel-1 A/B satellites that indicates the Altiplano-Puna magma body, located 19–24 km below ground level, and previously noted as the source of the large region of deformation, is still active. The ground is now uplifting at a rate of ~3 mm/yr compared to prior rates of ~10 mm/yr. We corroborated this waning uplift with in situ data from station UTUR. We combined the Sentinel-1 data with TerraSAR-X interferograms to constrain an ~25 km2 region of subsidence located 11 km SSW of Uturuncu, with a source depth of 2.1 km below ground level to an active period of ~2.5 yr with ~5 mm/yr subsidence. We developed a conceptual model that relates these varying depths and time scales of activity in a transcrustal magmatic system. We associate the surface uplift with pressurization from ascending gases and brines from magmatic reservoirs in the midcrust. We infer the existence of brine lenses in the shallow hydrothermal system based on low subsurface resistivity correlated with surface subsidence.more » « less
-
Abstract Surface deformation and seismicity provide critical information to understand the dynamics of volcanic unrest. During 2006–2007, >80 mm/yr uplift was observed by interferometric synthetic aperture radar (InSAR) at the central Atka volcanic center, Alaska, coinciding with an increasing seismicity rate. On November 25, 2006, a phreatic eruption occurred at the Korovin volcanic vent, 5‐km north of the central Atka, following the drainage of its crater lake a month prior to the eruption. The InSAR data are assimilated into three‐dimensional finite element models using the Ensemble Kalman Filter to investigate: (1) the pressure source creating the surface deformation; (2) the triggering of the volcano‐tectonic (VT) earthquakes in the Atka volcanic center; and (3) the triggering of the phreatic eruption at Korovin. The models show that the pressure source required to create the surface deformation is a NE‐tilted, oblate ellipsoid, which rotated from steep to gentle dipping from June to November 2006 before the eruption. The modeled dilatancy in a pre‐existing weak zone, coinciding with the Amlia‐Amukta fault, driven by the pressure source has a spatial and temporal correlation with the evolution of the VT earthquakes during the unrest. The fault dilatancy may have increased the connected porosity and permeability of the fault zone allowing fluid injection which triggered the observed seismicity. In addition, the dilatated fault may have increased the fluid capacity of the fault zone by ∼105 m3, causing the discharge of the crater lake at Korovin. Consequently, the phreatic eruption of the Korovin volcano may have been triggered.more » « less
-
null (Ed.)Finite-fault models for the 2010 M w 8.8 Maule, Chile earthquake indicate bilateral rupture with large-slip patches located north and south of the epicenter. Previous studies also show that this event features significant slip in the shallow part of the megathrust, which is revealed through correction of the forward tsunami modeling scheme used in tsunami inversions. The presence of shallow slip is consistent with the coseismic seafloor deformation measured off the Maule region adjacent to the trench and confirms that tsunami observations are particularly important for constraining far-offshore slip. Here, we benchmark the method of Optimal Time Alignment (OTA) of the tsunami waveforms in the joint inversion of tsunami (DART and tide-gauges) and geodetic (GPS, InSAR, land-leveling) observations for this event. We test the application of OTA to the tsunami Green’s functions used in a previous inversion. Through a suite of synthetic tests we show that if the bias in the forward model is comprised only of delays in the tsunami signals, the OTA can correct them precisely, independently of the sensors (DART or coastal tide-gauges) and, to the first-order, of the bathymetric model used. The same suite of experiments is repeated for the real case of the 2010 Maule earthquake where, despite the results of the synthetic tests, DARTs are shown to outperform tide-gauges. This gives an indication of the relative weights to be assigned when jointly inverting the two types of data. Moreover, we show that using OTA is preferable to subjectively correcting possible time mismatch of the tsunami waveforms. The results for the source model of the Maule earthquake show that using just the first-order modeling correction introduced by OTA confirms the bilateral rupture pattern around the epicenter, and, most importantly, shifts the inferred northern patch of slip to a shallower position consistent with the slip models obtained by applying more complex physics-based corrections to the tsunami waveforms. This is confirmed by a slip model refined by inverting geodetic and tsunami data complemented with a denser distribution of GPS data nearby the source area. The models obtained with the OTA method are finally benchmarked against the observed seafloor deformation off the Maule region. We find that all of the models using the OTA well predict this offshore coseismic deformation, thus overall, this benchmarking of the OTA method can be considered successful.more » « less
