1D velocity structure beneath broadband seismic stations in the Cretaceous Gyeongsang Basin of Korea by receiver function analyses

The CIFALPS receiver-function (RF) profile in the southwestern Alps provided the first seismological evidence of continental subduction in the Alps, with the detection of waves converted on the European Moho at 75-80 km depth beneath the western edge of the Po basin (Zhao et al., 2015). To complement the CIFALPS profile and enhance our knowledge of the lithospheric structure of the Western Alps, we installed CIFALPS2, a temporary network of 55 broadband seismic stations that operated for ~14 months (2018-2019) across the North-Western Alps (Zhao et al., 2018). The CIFALPS2 line runs from the Eastern Massif Central to the Ligurian coast, across the Mont-Blanc and Gran Paradiso massifs and the Ligurian Alps. Seismic stations were installed along a quasi-linear profile with a spacing of 7-10 km.We will show 2 receiver-function CCP (common-conversion point) depth-migrated sections along the CIFALPS2 profile, the first one across the Alps, and the second one across the Ligurian Alps and the Po basin. The time-to-depth migration of RF data is based on the new 3-D Vs model of the Greater Alpine region derived by Nouibat et al. (2021) using transdimensional ambient noise tomography on a large dataset including the AlpArray seismic network. Depth sections across the Vs model are also useful for interpreting the RF CCP sections as they have striking similarities.The images of the lithospheric structure of the NW Alps along CIFALPS2 are surprisingly different from those of the SW Alps along CIFALPS. The deepest P-to-S converted phases on the European Moho are detected at 60-65 km depth beneath the Ivrea-Verbano zone, that is 15 km less than on CIFALPS. The negative polarity converted phase interpreted as the base of the Ivrea body mantle flake on the CIFALPS section is still visible on CIFALPS2, but with a lower amplitude. The RF section confirms the existence of a jump of the European Moho of ~10 km amplitude in less than 10 km distance, which is located within a few km from the western boundary of the Mont Blanc external crystalline massif. All these observations are confirmed by the Vs model that also displays a less deep continental subduction than on CIFALPS, weaker S-wave velocities in the Ivrea body wedge, and the jump of the European Moho.The Moho beneath the Ligurian Alps is detected at 25-30 km depth both on the RF and on the Vs depth sections. Moving northwards, this Ligurian Moho is separated from the Adriatic Moho by a puzzling S-dipping set of P-to-S converted waves with negative polarity. The crust of the Ligurian Alps is characterized by a set of north-dipping negative-polarity converted waves at 10 to 20 km depth beneath the Valosio massif, which is a small internal crystalline massif of (U)HP metamorphic rocks located north of Voltri. The similarity of this set of negative-polarity conversions to the one observed beneath the Dora Maira massif on the CIFALPS profile suggests that it may be a relic of the Alpine structure overprinted by the opening of the Ligurian sea.

Download Full-text

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).

10.5194/egusphere-egu21-9882 ◽

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

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.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.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 &#160;density and VS 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.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.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.

Download Full-text

Seismic Structure and Tectonic Evolution of Borneo and Sulawesi

10.5194/egusphere-egu21-16456 ◽

2021 ◽

Author(s):

Harry Telajan Linang ◽

Amy Gilligan ◽

Jennifer Jenkins ◽

Tim Greenfield ◽

Felix Tongkul ◽

...

Keyword(s):

Subduction Zones ◽

Receiver Function ◽

Leading Edge ◽

Sedimentary Basins ◽

Crustal Thickness ◽

The Philippines ◽

Seismic Structure ◽

Eurasian Plate ◽

Seismic Stations ◽

Preliminary Results

<div> <div> <div> Borneo is located at the centre of Southeast Asia, which is one of the most active tectonic regions on Earth due to the subduction of the Indo-Australian plate in the south and the Philippines Sea plate in the east. Borneo resides on the leading edge of the Sundaland block of the Eurasian plate and exhibits lower rates of seismicity when compared to the surrounding regions due to its intraplate setting. Sulawesi, an island which lies just southeast of Borneo, is characterised by intense seismicity due to multiple subduction zones in its vicinity. The tectonic relationship between the two islands is poorly understood, including the provenance of their respective lithospheres, which may have Eurasian and/or East Gondwana origin. Here, we present recent receiver function (RF) results from temporary and permanent broadband seismic stations in the region, which can be used to help improve our understanding of the crust and mantle lithosphere beneath Borneo and Sulawesi. We applied H-K stacking, receiver function migration and inversion to obtain reliable estimates of the crustal thickness beneath the seismic stations. Our preliminary results indicate that the crust beneath Sabah (in northern Borneo), which is a post-subduction setting, appears to be much more complex and is overall thicker (more than 35 km) than the rest of the island. In addition, we find that crustal thickness varies between different tectonic blocks defined from previous surface mapping, with the thinnest crust (23 to 25 km) occurring beneath Sarawak in the west-northwest as well as in the east of Kalimantan. We also present preliminary results from Virtual Deep Seismic Sounding (VDSS) in northern Borneo, where from the RF results we know that there is thick and complex crust. VDSS is able to produce well constrained crustal thickness results in regions where the RF analysis has difficulty recovering the Moho, likely due to complexities such as thick sedimentary basins and obducted ophiolite sequences. </div> </div> </div>

Download Full-text

A new inversion method to construct a 3-D crustal shear-wave velocity model from P-to-S converted waves and application to the Central Alps

10.5194/egusphere-egu2020-10386 ◽

2020 ◽

Author(s):

Leonardo Colavitti ◽

György Hetényi ◽

AlpArray Working Group

Keyword(s):

Simulated Annealing ◽

Wave Velocity ◽

Shear Wave Velocity ◽

Ambient Noise ◽

Velocity Structure ◽

Receiver Function ◽

Velocity Model ◽

Central Alps ◽

Travelling Salesman ◽

The Given

We developed a new method where teleseismic P-to-S converted waves are used to construct a fully 3-D shear-wave velocity model of the crust. The method differs from ambient noise and local earthquake tomography in its ray-paths being closer to vertical. Our approach requires a dense seismological network, and we first focus on the Central Alps considering the available permanent and temporary station datasets (e.g., Het&#233;nyi et al., 2018, Surv. Geophys.).We implemented an accurate ray-propagator which respects Snell&#8217;s law in 3-D at any interface geometry. Following a teleseismic P ray propagator (Knapmeyer, 2004) from event to station which uses a 1-D global velocity model (iasp91), P-to-S conversion at the Moho is calculated for the crustal S ray considering the true local dip. The corresponding arrival to the surface is typically several km away from the station, which we then adjust by changing the ray-parameter.&#160;In the Central Alps, using the 3-D P-velocity structure of Diehl et al. (2009) and the local Moho geometry of Spada et al. (2013), the mean distance between the arriving S-wave and the station is about 150 m (median ca. 40 m).For our approach we adopt a new model parameterization of velocities. It is rectangular in map view (nodes at 25x25 km in the Alps), while in depth we define a 2-layer model with separate velocities above and below each discontinuity. The introduction of this flexibility allows us to accommodate a velocity gradient within each layer and investigate velocity jumps across discontinuities.The inversion proceeds iteratively, by visiting every node of the map following a Travelling Salesman Path. At each node, receiver function rays in the surrounding volume are considered for inversion, and bundled into sub-blocks and ranges of back-azimuth (5x5 km size, 45&#176; or 60&#176; bins for the Central Alps). The velocity model at the given node is inverted using the technique of Simulated Annealing, followed by a pattern search algorithm to avoid falling in a local minimum. During iterations of the Simulated Annealing, individual velocity model corresponding to each receiver function is extracted from the 3-D model along its ray path.The inversion proceeds for 4 or 5 independent parameters: Moho and a hypothetical intra-crustal discontinuity depth, Vp/Vs ratio (either full crust, or separately for upper and lower crust) and the P-wave velocity jump at the intra-crustal discontinuity. Finally, the velocity structure is updated with the result obtained at the given node. We observe that a few rounds of Travelling Salesman Paths improve the overall misfit.First results on the Central Alps show that the Moho depth generally reflects well the roots of the Alpine orogen. Resolving crustal Vp/Vs ratio is more stable when considering the full crust, instead of two separate layers. The Conrad discontinuity remains difficult to resolve. The obtained velocity structure is compared along profiles to recent Vs results from 3-D ambient noise tomography (Lu et al., 2018).

Download Full-text