Joint receiver function and gravity inversion for new constraints on the Ivrea body structure: a 2D high-resolution view along the Val Sesia profile (N. Italy).

2021 ◽  
Author(s):  
Matteo Scarponi ◽  
György Hetényi ◽  
Jaroslava Plomerová ◽  
Stefano Solarino
Keyword(s):  
High Resolution ◽  
Velocity Structure ◽  
Receiver Function ◽  
Joint Inversion ◽  
Gravity Data ◽  
Rock Physics ◽  
Receiver Functions ◽  
Gravity Inversion ◽  
European Alps ◽  
Horizontal Distance

<p>We present results from a joint inversion study of new seismic and gravity data to constrain a 2D high-resolution image of one of the most prominent geophysical anomalies of the European Alps: the Ivrea geophysical body (IGB). Our work exploits both new data and multidisciplinary a priori constraints, to better resolve the shallow crustal structure in the Ivrea-Verbano zone (IVZ), where the IGB is known to reach anomalously shallow depths and partially outcrop at the surface.</p><p>A variety of previous studies, ranging from gravity surveys to vintage refraction seismics and recent local earthquake tomographies (Solarino et al. 2018, Diehl et al. 2009), provide comprehensive but spatially sparse information on the IGB structure, which we aim at investigating at higher resolution, along a linear profile crossing the IVZ. To this purpose, we deployed 10 broadband seismic stations (MOBNET pool, IG CAS Prague), 5 km spaced along a linear West-East profile, along Val Sesia and crossing Lago Maggiore. This network operated for 27 months and allowed us to produce a new database of ca. 1000 seismic high-quality receiver functions (RFs). In addition, we collected new gravity data in the IVZ, with a data coverage of 1 gravity point every 1-2 km along the seismic profile. The newly collected data was used to set up an inversion scheme, in which RFs and gravity anomalies are jointly used to constrain the shape and the physical property contrasts across the IGB interface.</p><p>We model the IGB as a single interface between far-field constraints, whose geometry is defined by the coordinates of four nodes which may vary in space, and  density and V<sub>S</sub> shear-wave velocity contrasts associated with the interface itself, varying independently. A Markov chain Monte Carlo (MCMC) sampling method with Metropolis-Hastings selection rule was implemented to efficiently explore the model space, directing the search towards better fitting areas.</p><p>For each model, we perform ray-tracing and RFs migration using the actual velocity structure both for migration and computation of synthetic RFs, to be compared with the observations via cross-correlation of the migration images. Similarly, forward gravity modelling for a 2D density distribution is implemented and the synthetic gravity anomaly is compared with the observations along the profile. The joint inversion performance is the product of these two misfits.</p><p>The inversion results show that the IGB reaches the shallowest depths in the western part of the profile, preferentially locating the IGB interface between 3 and 7 km depth over a horizontal distance of ca. 20 km (between Boccioleto and Civiasco, longitudes 8.1 and 8.3). Within this segment, the shallowest point reaches up to 1 km below sea level. The found density and velocity contrasts are in agreement with rock physics properties of various units observed in the field and characterized in earlier studies.</p>

scholarly journals Joint Seismic and Gravity Data Inversion to Image Intra-Crustal Structures: The Ivrea Geophysical Body Along the Val Sesia Profile (Piedmont, Italy)

2021 ◽  
Vol 9 ◽  
Author(s):  
Matteo Scarponi ◽  
György Hetényi ◽  
Jaroslava Plomerová ◽  
Stefano Solarino ◽  
Ludovic Baron ◽  
...  
Keyword(s):  
Physical Properties ◽  
Joint Inversion ◽  
Gravity Data ◽  
Model Space ◽  
Vertical Gradient ◽  
Receiver Functions ◽  
Gravity Modeling ◽  
Thickness Variation ◽  
European Alps ◽  
Velocity Discontinuity

We present results from a joint inversion of new seismic and recently compiled gravity data to constrain the structure of a prominent geophysical anomaly in the European Alps: the Ivrea Geophysical Body (IGB). We investigate the IGB structure along the West-East oriented Val Sesia profile at higher resolution than previous studies. We deployed 10 broadband seismic stations at 5 km spacing for 27 months, producing a new database of ∼1000 high-quality seismic receiver functions (RFs). The compiled gravity data yields 1 gravity point every 1–2 km along the profile. We set up an inversion scheme, in which RFs and gravity anomalies jointly constrain the shape and the physical properties of the IGB. We model the IGB’s top surface as a single density and shear-wave velocity discontinuity, whose geometry is defined by four, spatially variable nodes between far-field constraints. An iterative algorithm was implemented to efficiently explore the model space, directing the search toward better fitting areas. For each new candidate model, we use the velocity-model structures for both ray-tracing and observed-RFs migration, and for computation and migration of synthetic RFs: the two migrated images are then compared via cross-correlation. Similarly, forward gravity modeling for a 2D density distribution is implemented. The joint inversion performance is the product of the seismic and gravity misfits. The inversion results show the IGB protruding at shallow depths with a horizontal width of ∼30 km in the western part of the profile. Its shallowest segment reaches either 3–7 or 1–3 km depth below sea-level. The latter location fits better the outcropping lower crustal rocks at the western edge of the Ivrea-Verbano Zone. A prominent, steep eastward-deepening feature near the middle of the profile, coincident with the Pogallo Fault Zone, is interpreted as inherited crustal thickness variation. The found density and velocity contrasts of the IGB agree with physical properties of the main rock units observed in the field. Finally, by frequency-dependent analysis of RFs, we constrain the sharpness of the shallowest portion of the IGB velocity discontinuity as a vertical gradient of thickness between 0.8 km and 0.4 km.

scholarly journals Reliable workflow for inversion of seismic receiver function and surface wave dispersion data: a “13 BB Star” case study

2019 ◽  
Vol 24 (1) ◽  
pp. 101-120
Author(s):  
Kajetan Chrapkiewicz ◽  
Monika Wilde-Piórko ◽  
Marcin Polkowski ◽  
Marek Grad
Keyword(s):  
Surface Wave ◽  
Inverse Problems ◽  
Receiver Function ◽  
Joint Inversion ◽  
Wave Dispersion ◽  
Receiver Functions ◽  
P Wave ◽  
Surface Wave Dispersion ◽  
Phase Velocity Dispersion ◽  
Wide Range

AbstractNon-linear inverse problems arising in seismology are usually addressed either by linearization or by Monte Carlo methods. Neither approach is flawless. The former needs an accurate starting model; the latter is computationally intensive. Both require careful tuning of inversion parameters. An additional challenge is posed by joint inversion of data of different sensitivities and noise levels such as receiver functions and surface wave dispersion curves. We propose a generic workflow that combines advantages of both methods by endowing the linearized approach with an ensemble of homogeneous starting models. It successfully addresses several fundamental issues inherent in a wide range of inverse problems, such as trapping by local minima, exploitation of a priori knowledge, choice of a model depth, proper weighting of data sets characterized by different uncertainties, and credibility of final models. Some of them are tackled with the aid of novel 1D checkerboard tests—an intuitive and feasible addition to the resolution matrix. We applied our workflow to study the south-western margin of the East European Craton. Rayleigh wave phase velocity dispersion and P-wave receiver function data were gathered in the passive seismic experiment “13 BB Star” (2013–2016) in the area of the crust recognized by previous borehole and refraction surveys. Final models of S-wave velocity down to 300 km depth beneath the array are characterized by proximity in the parameter space and very good data fit. The maximum value in the mantle is higher by 0.1–0.2 km/s than reported for other cratons.

scholarly journals P-wave velocity and density structure beneath Mt. Vesuvius: a magma body in the upper edifice?

2013 ◽  
Vol 56 (4) ◽  
Author(s):  
Paolo Capuano ◽  
Guido Russo ◽  
Roberto Scarpa
Keyword(s):  
High Resolution ◽  
Wave Velocity ◽  
Gravity Data ◽  
High Energy ◽  
P Wave ◽  
Gravity Inversion ◽  
Density Structure ◽  
Resolution Image ◽  
High Resolution Image ◽  
P Wave Velocity

<p>A high-resolution image of the compressional wave velocity and density structure in the shallow edifice of Mount Vesuvius has been derived from simultaneous inversion of travel times and hypocentral parameters of local earthquakes and from gravity inversion. The robustness of the tomography solution has been improved by adding to the earthquake data a set of land based shots, used for constraining the travel time residuals. The results give a high resolution image of the P-wave velocity structure with details down to 300-500 m. The relocated local seismicity appears to extend down to 5 km depth below the central crater, distributed into two clusters, and separated by an anomalously high Vp region positioned at around 1 km depth. A zone with high Vp/Vs ratio in the upper layers is interpreted as produced by the presence of intense fluid circulation alternatively to the interpretation in terms of a small magma chamber inferred by petrologic studies. In this shallower zone the seismicity has the minimum energy, whilst most of the high-energy quakes (up to Magnitude 3.6) occur in the cluster located at greater depth. The seismicity appears to be located along almost vertical cracks, delimited by a high velocity body located along past intrusive body, corresponding to remnants of Mt. Somma. In this framework a gravity data inversion has been performed to study the shallower part of the volcano. Gravity data have been inverted using a method suitable for the application to scattered data in presence of relevant topography based on a discretization of the investigated medium performed by establishing an approximation of the topography by a triangular mesh. The tomography results, the retrieved density distribution, and the pattern of relocated seismicity exclude the presence of significant shallow magma reservoirs close to the central conduit. These should be located at depth higher than that of the base of the hypocenter volume, as evidenced by previous studies.</p>

Crustal Thickness and Composition of the São Paulo Plateau and Florianópolis Ridge, SE Brazilian Margin

2020 ◽  
Author(s):  
Michelle Graça ◽  
Leanne Cowie ◽  
Nick Kusznir ◽  
Natasha Stanton
Keyword(s):  
Oceanic Crust ◽  
Seismic Reflection ◽  
Joint Inversion ◽  
Gravity Data ◽  
South Atlantic ◽  
Gravity Inversion ◽  
The South ◽  
Subsidence Analysis ◽  
Base Salt ◽  
Thermal Subsidence

&lt;p&gt;The S&amp;#227;o Paulo Plateau (SPP) and the Florian&amp;#243;polis Ridge (FR), located on the Santos segment of the SE Brazilian margin in the South Atlantic, are large positive bathymetric features with a combined lateral dimension of approximately 500 km. An important question is whether they are underlain by thinned continental crust or by anomalously thick magmatic crust. Each hypothesis has implications for the breakup of the South Atlantic and the evolution of the overlying saline Santos basin.&lt;/p&gt;&lt;p&gt;Integrated quantitative analysis consisting of gravity inversion, RDA (residual depth anomaly) analysis and flexural subsidence analysis has been applied to a deep long-offset seismic reflection line running NW-SE across the SPP and FR. Gravity inversion predicts crustal basement thicknesses in the range of 12 to 15 km for the SPP and FR, deceasing to 7-8 km thickness at the extreme SE end of the profile. The SPP and FR are separated by a region of thinner crust approximately 80 km wide. Thinning factors from subsidence analysis for SPP and FR are typically between 0.6 and 0.7.&lt;/p&gt;&lt;p&gt;RDA values close to zero and a thinning factor of 1 were obtained for the region with 7-8 km thick crust at the SE end of the profile which are all consistent with normal oceanic crust rather than previously interpreted exhumed mantle. This oceanic crust is highly tectonised and corresponds to the location of the Florian&amp;#243;polis Fracture Zone.&lt;/p&gt;&lt;p&gt;Flexural backstripping and reverse thermal subsidence modelling were performed to calculate palaeo-bathymetry at breakup and give 2.5 km below sea level at the SE end of the profile consistent with this region being oceanic crust. Flexural subsidence analysis applied to base salt shows that the observed base salt subsidence requires a component of syn-tectonic subsidence as well as post-rift thermal subsidence, and that the salt was deposited while the crust was still thinning.&lt;/p&gt;&lt;p&gt;Joint inversion of time seismic reflection and gravity data to determine the lateral variation in basement density by comparing seismic and gravity Moho in the time domain gives a basement density under the SPP and FR of between 2600 and 2700 kg/m&lt;sup&gt;3&lt;/sup&gt;. The same method gives a basement density of 900kg/m&lt;sup&gt;3&lt;/sup&gt; for the oceanic crust at the SE end of the profile. The FR basement in the NW shows a basement density similar to that of the SPP while in its SE the basement density is much higher approaching 2950 kg/m3.&amp;#160; We interpret the relatively low basement densities of the SPP with respect to that of oceanic crust as indicating a continental rather than magmatic composition. A similar analysis to determine basement density applied to the Evain et al. (2015) seismic refraction profile in the same location also gives a SPP basement density that supports a continental composition.&lt;/p&gt;

Receiver function analyses and joint inversion with gravity data: new constraints on the Ivrea geophysical body along a high-resolution profile

2020 ◽  
Author(s):  
Matteo Scarponi ◽  
György Hetényi ◽  
Jaroslava Plomerová ◽  
Stefano Solarino ◽  
Ludovic Baron
Keyword(s):  
Receiver Function ◽  
Joint Inversion ◽  
Gravity Data ◽  
Sampling Rate ◽  
Seismic Refraction ◽  
Western Alps ◽  
Transverse Component ◽  
Reference Points ◽  
Physical Parameters ◽  
S Wave

&lt;p&gt;We collected new seismological and gravity data in the Val Sesia and Lago Maggiore regions in NW Italy to constrain the geometry and properties of the Ivrea Geophysical Body. This piece of lower Adriatic lithosphere is known to be at anomalously shallow depth along the inner arc of the Western Alps, yet existing seismological constraints (vintage seismic refraction data, local earthquake tomography) are spatially sparse. With the aim to reach higher spatial resolution in imaging the structure of the IGB, we analyze the seismological data with various receiver function approaches to map the main velocity discontinuities, followed by joint inversion with gravity data to fill the bulk properties of bodies with densities.&lt;/p&gt;&lt;p&gt;The new data acquisition consisted of two type of campaigns. For seismology, we deployed 10 broadband seismic stations (MOBNET pool, IG CAS Prague) along a linear West-East profile at 5 km spacing along Val Sesia and across the Lago Maggiore. This network continuously recorded seismic data for 27 months at 100 Hz sampling rate. For gravimetry, we compiled existing datasets and then completed the spatial gaps by relative gravity surveys, tied to absolute reference points, to achieve 1 gravity point every 1-2 km along the profile.&lt;/p&gt;&lt;p&gt;The receiver function (RF) analyses aim at detecting velocity increases with depth: primarily the Moho and the shallow IGB interfaces and their crustal reverberations (multiples), together with their potential dip by analyzing the transverse component RFs. Furthermore, we aim at investigating the sharpness of the velocity gradient across the discontinuities by analyzing the frequency dependence of the corresponding RF peaks. We aim at reproducing the observations by simple synthetic models.&lt;/p&gt;&lt;p&gt;The 2D joint inversion combines S wave velocity V&lt;sub&gt;S&lt;/sub&gt; and bulk density as physical parameters to match both the seismological and gravimetry data. The relationship between the two parameters is initially chosen from the literature, but depending on the first results the relation itself may be inverted for, considering the various high-grade metamorphic rocks observed at the surface in the area, whose properties may not align with classical V&lt;sub&gt;S&lt;/sub&gt;&amp;#8211;density equations. In conclusion, we propose new constraints on the IGB, demonstrating the advantage of using multi-disciplinary geophysical observations and improved data coverage across the study area.&lt;/p&gt;

Lithospheric structure of the Arabian Shield from joint inversion of P- and S-wave receiver functions and dispersion velocities

2015 ◽  
Vol 65 (2) ◽  
pp. 239-255 ◽  
Author(s):  
Abdullah M. Al-Amri
Keyword(s):  
Green's Functions ◽  
Velocity Structure ◽  
Green’S Functions ◽  
Joint Inversion ◽  
Receiver Functions ◽  
Gulf Of Aqaba ◽  
Arabian Shield ◽  
S Wave ◽  
Velocity Models ◽  
Seismic Networks

Abstract New velocity models of lithospheric thickness and velocity structure have been developed for the Arabian Shield by three tasks: 1) Computing P-Wave Receiver Functions (PRFs) and S-Wave Receiver Functions (SRFs) for all the broadband stations within the Saudi seismic networks. The number of receiver function waveforms depends on the recording time window and quality of the broadband station. 2) Computing ambient noise correlation Green’s functions for all available station pairs within the Saudi seismic networks to image the shear velocity in the crust and uppermost mantle beneath the Arabian Peninsula. Together they provided hundreds of additional, unique paths exclusively sampling the region of interest. Both phase and group velocities for all the resulting empirical Green’s functions have been measured and to be used in the joint inversion. 3) Jointly inverted the PRFs and SRFs obtained in task 1 with dispersion velocities measured on the Green’s functions obtained in task 2 and with fundamental-mode, Rayleigh-wave, group and phase velocities borrowed from the tomographic studies to precisely determine 1D crustal velocity structure and upper mantle. The analysis of the PRFs revealed values of 25-45 km for crustal thickness, with the thin crust next to the Red Sea and Gulf of Aqaba and the thicker crust under the platform, and Vp/Vs ratios in the 1.70-1.80 range, suggesting a range of compositions (felsic to mafic) for the shield’s crust. The migrated SRFs suggest lithospheric thicknesses in the 80-100 km range for portions of the shield close to the Red Sea and Gulf of Aqaba and near the Arabian Gulf. Generally, the novelty of the velocity models developed under this paper has consisted in the addition of SRF data to extend the velocity models down to lithospheric and sub-lithospheric depths.

The Signature and Elimination of Sediment Reverberations on Submarine Receiver Functions: Imaging the Lithosphere of a Normal Ocean

2021 ◽  
Author(s):  
Ziqi Zhang ◽  
Tolulope Olugboji
Keyword(s):  
High Resolution ◽  
Receiver Function ◽  
Receiver Functions ◽  
Function Technique ◽  
Ocean Bottom ◽  
Sediment Layer ◽  
Waveform Modeling ◽  
Seismic Arrays ◽  
Filter Approach ◽  
Synthetic Waveform

&lt;p&gt;While the receiver function technique has been successfully applied to high-resolution imaging of sharp discontinuities within and across the lithosphere, it has been shown, however, that it suffers from severe limitations when applied to seafloor seismic recordings. This is because the water and sediment layer could strongly influence the receiver function traces, making detection and interpretation of crust and mantle layering difficult. This effect is often referred to as the singing phenomena in marine environments. Here, we show how one can silence this singing effect. We demonstrate, using analytical and synthetic waveform modeling, that this singing effect can be reversed using dereverberation filters tuned to match the elastic property of each layer. We apply the filter approach to high-quality earthquake records collected from the NoMelt seismic array deployed on normal, mature (~70 Ma) Pacific seafloor. An appropriate filter designed using the elastic properties of the underlying sediments, and obtained from prior studies, greatly improves the detection of Ps conversions generated from the moho (~8.6 km) and from a sharp discontinuity (&lt;~ 5 km) across the lithosphere-asthenosphere transition (~72 km). Sensitivity tests show that the filter is robust to small errors in the sediment properties. Our analysis suggests that appropriately filtering out the sediment reverberations from ocean seismic data could make inferences on subsurface structure more robust. We expect that this study will enable high-resolution receiver function imaging of the base of the oceanic plate across a growing fleet of ocean bottom seismic arrays being deployed in the global oceans.&lt;/p&gt;

Sign in / Sign up

Export Citation Format

Share Document