Evidence of a shallow persistent magmatic reservoir from joint inversion of gravity and ground deformation data: The 25–26 October 2013 Etna lava fountaining event

Abstract The 3 September 2017 Mw 5.2 North Korean underground nuclear test (DPRK2017) is the largest man-made explosion with surface displacements observed by Synthetic Aperture Radar (SAR) and showed as much as 3.5 m of horizontal permanent deformation. Although regional distance waveform-based seismic moment tensor (MT) inversion methods successfully identify this event as an explosion, the inverted solutions do not fit the SAR displacement field well. To better constrain the source, we developed an MT source-type inversion method that incorporates surface ground deformation (accounting for free-surface topography), regional seismic waveforms, and first-motion polarities. We applied the source-type inversion over a grid of possible source locations to find the best-fitting location, depth, and point-source MT for the event. Our best-fitting MT solution achieves ∼70% horizontal geodetic fit, ∼80% waveform fit, and 100% fit in the first-motion polarities. The joint inversion narrows the range of acceptable source types improving discrimination, and reduces the uncertainty in scalar moment and estimated yield. The method is transportable and can be applied to other types of events that may have measurable geodetic signals such as underground mine collapses and volcanic events.

Download Full-text

New insight into the December 2018 Etna eruption through the joint inversion of ground deformation and gravity data

10.5194/egusphere-egu2020-564 ◽

2020 ◽

Author(s):

Mahak Singh Chauhan ◽

Flavio Cannavò ◽

Daniele Carbone ◽

Filippo Greco

Keyword(s):

Search Algorithm ◽

Ground Deformation ◽

Joint Inversion ◽

Gravity Data ◽

Monitoring Network ◽

Pattern Search ◽

Gravity Change ◽

Model Parameters ◽

Geometrical Parameters ◽

Superconducting Gravimeters

<p>We focus on the eruption of Mt. Etna which took place on 24 December, 2018. The eruption occurred after a month of unrest and was accompanied by a seismic swarm that culminated in the M4.9 earthquake on the 26th, with epicentre on the eastern flank of the volcano. We jointly analyse ground deformation and gravity data to estimate the geometrical and kinematic parameters of the source structure, together with the density of the intruding material. The data used in this study were recorded by stations in the INGV-OE monitoring network (21 GPS stations and 2 gravity stations equipped with superconducting gravimeters), during the interval of 23 to 28 December (pre to post eruption). We assume a dike-type source for the forward calculation in the defined objective function. A pattern search algorithm (PSA) is used for the iterative minimization of the misfit error. In order to estimate the posterior probability density function (PDFs) of the model parameters, we also use a Markov Chain Monte Carlo (MCMC) approach. Indeed, the calculated PDFs provide more information about the uncertainties of the model parameters, which helps to understand overall tendencies of the solutions. We first test the constrained inversion of the gravity data, to calculate the density of eruptive magmatic body, by fixing the geometrical parameters of the dike, previously retrieved through inversion of the deformation data only. Using this approach, it is possible to suitable explain the deformation data and the gravity change observed at the station in the near field (MNT), while the gravity change at the other station (SLN) remain unexplained. We then invert jointly both deformation and gravity datasets, in order to adequately fit all the observations. The final model gives a density value of ~1.8-2.0 g/cm3. This value is significantly lower than the density of bubble-free magma and indicates either the involvement of gas in the intrusive process, or the formation of dry fissures during the emplacement of the dyke.</p>

Download Full-text

Failure Analysis Of Buried Gas Pipelines Crossing Seismic Faults

Academic Perspective Procedia ◽

10.33793/acperpro.03.02.2 ◽

2020 ◽

Vol 3 (2) ◽

pp. 781-790

Author(s):

M. Rizwan Akram ◽

Ali Yesilyurt ◽

A.Can. Zulfikar ◽

F. Göktepe

Keyword(s):

Ground Deformation ◽

Warning System ◽

Strike Slip ◽

Gas Pipelines ◽

Steel Pipes ◽

Seismic Fault ◽

Fault Parameters ◽

Dip Angle ◽

Permanent Ground Deformation ◽

Recent Earthquakes

Research on buried gas pipelines (BGPs) has taken an important consideration due to their failures in recent earthquakes. In permanent ground deformation (PGD) hazards, seismic faults are considered as one of the major causes of BGPs failure due to accumulation of impermissible tensile strains. In current research, four steel pipes such as X-42, X-52, X-60, and X-70 grades crossing through strike-slip, normal and reverse seismic faults have been investigated. Firstly, failure of BGPs due to change in soil-pipe parameters have been analyzed. Later, effects of seismic fault parameters such as change in dip angle and angle between pipe and fault plane are evaluated. Additionally, effects due to changing pipe class levels are also examined. The results of current study reveal that BGPs can resist until earthquake moment magnitude of 7.0 but fails above this limit under the assumed geotechnical properties of current study. In addition, strike-slip fault can trigger early damage in BGPs than normal and reverse faults. In the last stage, an early warning system is proposed based on the current procedure.&nbsp;

Download Full-text

Joint inversion of acoustic and resistivity data for the estimation of gas hydrate concentration

10.3133/b2190 ◽

2002 ◽

Keyword(s):

Gas Hydrate ◽

Joint Inversion ◽

Resistivity Data

Download Full-text

Maps showing landslide features and related ground deformation in the Woodlawn area of the Manoa Valley, City and County of Honolulu, Hawaii

Open-File Report ◽

10.3133/ofr89290 ◽

1989 ◽

Cited By ~ 1

Author(s):

Rex L. Baum ◽

Robert W. Fleming ◽

Stephen D. Ellen

Keyword(s):

Ground Deformation ◽

Valley City

Download Full-text

Joint Inversion of Hydrologic and Geophysical Data for Permeability Distribution of an Alluvial Aquifer

10.21236/ada358024 ◽

1998 ◽

Cited By ~ 2

Author(s):

Warren Barrash

Keyword(s):

Joint Inversion ◽

Alluvial Aquifer ◽

Geophysical Data ◽

Permeability Distribution

Download Full-text

Imaging the subsurface using induced seismicity and ambient noise: 3-D tomographic Monte Carlo joint inversion of earthquake body wave traveltimes and surface wave dispersion

Geophysical Journal International ◽

10.1093/gji/ggaa230 ◽

2020 ◽

Vol 222 (3) ◽

pp. 1639-1655

Author(s):

Xin Zhang ◽

Corinna Roy ◽

Andrew Curtis ◽

Andy Nowacki ◽

Brian Baptie

Keyword(s):

Surface Wave ◽

Ambient Noise ◽

Velocity Structure ◽

Joint Inversion ◽

Body Wave ◽

Source Parameters ◽

Wave Dispersion ◽

Inversion Method ◽

Surface Wave Dispersion ◽

Wave Data

SUMMARY Seismic body wave traveltime tomography and surface wave dispersion tomography have been used widely to characterize earthquakes and to study the subsurface structure of the Earth. Since these types of problem are often significantly non-linear and have non-unique solutions, Markov chain Monte Carlo methods have been used to find probabilistic solutions. Body and surface wave data are usually inverted separately to produce independent velocity models. However, body wave tomography is generally sensitive to structure around the subvolume in which earthquakes occur and produces limited resolution in the shallower Earth, whereas surface wave tomography is often sensitive to shallower structure. To better estimate subsurface properties, we therefore jointly invert for the seismic velocity structure and earthquake locations using body and surface wave data simultaneously. We apply the new joint inversion method to a mining site in the United Kingdom at which induced seismicity occurred and was recorded on a small local network of stations, and where ambient noise recordings are available from the same stations. The ambient noise is processed to obtain inter-receiver surface wave dispersion measurements which are inverted jointly with body wave arrival times from local earthquakes. The results show that by using both types of data, the earthquake source parameters and the velocity structure can be better constrained than in independent inversions. To further understand and interpret the results, we conduct synthetic tests to compare the results from body wave inversion and joint inversion. The results show that trade-offs between source parameters and velocities appear to bias results if only body wave data are used, but this issue is largely resolved by using the joint inversion method. Thus the use of ambient seismic noise and our fully non-linear inversion provides a valuable, improved method to image the subsurface velocity and seismicity.

Download Full-text