scholarly journals Comparison of P-Wave Tomography Model in Sunda-Banda Arc by Using Data from BMKG and ISC

2021 ◽  
Vol 873 (1) ◽  
pp. 012058
Author(s):  
P T Brilianti ◽  
M S Haq ◽  
Haolia ◽  
M I Sulaiman ◽  
R P Nugroho ◽  
...  

Abstract The tectonic setting of our study area is located between the Island of Java and Timor Leste. The complexity of this region is started with two different plates, The Indo-Australian plate and the Eurasian plate that move with different orientations and convergence rate. This area also shows active seismic activity and has a series of active volcanoes as a product of subduction and collision. To deepen understand this area, we perform delay time tomography using FMTOMO package that includes 3-D finite-difference based ray tracing and sub-space inversion procedure. We used two different sets of data, the first one is 4 years data catalog from the Indonesian Agency of Meteorology, Climatology and Geophysics, and the second one is 47 years of data from the International Seismological Centre. Data from the local Indonesian show agency shows a fewer number of events but more focus clusters. Meanwhile, the data from ISC catalog has more events and evenly distributed data. However, we also noticed that data from ISC has cluster events located at the same depth that can be improved with events relocation for better depth estimation. The Checkerboard models from both data set show a comparable result, though data from ISC show a better recovered model at a deeper depth and shallow part in the eastern area. The checkerboard from the local Agency shows slightly better results in the shallow part. Next, we invert delay time for each data set using we optimized damping and smoothing parameters. Final tomogram models show that data from the local Agency show a more continuous fast velocity band representing a downgoing subducting slab and possible back-arc thrust while results from the ISC data show a more detached fast velocity band that could be contributed from fixed depth problem in the data set. However, we noticed that data from ISC show a higher amplitude low-velocity anomaly especially in the shallow depth

2021 ◽  
Vol 873 (1) ◽  
pp. 012056
Author(s):  
M F R Auly ◽  
A K Ilahi ◽  
I Madrinovella ◽  
S Widyanti ◽  
S K Suhardja ◽  
...  

Abstract The tectonic setting of Java island, located at southwestern edge of the Eurasia continent, is dominated by the subduction of Indo-Australia plate. One of the characteristics of active subduction is active seismicity, the generation of arc magmatism and volcanic activity. Mt. Merapi is one example of active volcano related with the subduction process. It is one of the most active volcanoes with location close to high population area. To better understand this area, we employed the Receiver Function technique, a method to image sub surface structure by removing the vertical component from horizontal component. First, we collected high magnitude events and processed RF with water level deconvolution method. Then, we constructed synthetic model with initial velocity input from previous tomography model. Note that we used reflectivity method in generating synthetic model with input parameters matched with parameters from real data processing. Next, we adjusted velocity inputs mainly on tops sediments (1-3 km) to include sediment layers and volcanic rocks, mid-depth low velocity zone that may be related with magma chamber and depth of crust-mantle boundary. Current forward velocity models show a relatively good agreement from 3 stations (ME25, ME32 and ME36). We estimate a thin layer of sediments followed a zone of velocity layer at a depth of 10-15 km and crust-mantle boundary ranging from 26-29 km. In this study, simulated that the signal of sediments layer and low velocity layers interfere main crust mantle boundary that supposed to be highest signal after the P wave in the typical receiver function study.


1996 ◽  
Vol 39 (6) ◽  
Author(s):  
C. Chiarabba ◽  
A. Amato

In this paper we provide P-wave velocity images of the crust underneath the Apennines (Italy), focusing on the lower crustal structure and the Moho topography. We inverted P-wave arrival times of earthquakes which occurred from 1986 to 1993 within the Apenninic area. To overcome inversion instabilities due to noisy data (we used bulletin data) we decided to resolve a minimum number of velocity parameters, inverting for only two layers in the crust and one in the uppermost mantle underneath the Moho. A partial inversion of only 55% of the overall dataset yields velocity images similar to those obtained with the whole data set, indicating that the depicted tomograms are stable and fairly insensitive to the number of data used. We find a low-velocity anomaly in the lower crust extending underneath the whole Apenninic belt. This feature is segmented by a relative high-velocity zone in correspondence with the Ortona-Roccamonfina line, that separates the northern from the southern Apenninic arcs. The Moho has a variable depth in the study area, and is deeper (more than 37 km) in the Adriatic side of the Northern Apennines with respect to the Tyrrhenian side, where it is found in the depth interval 22-34 km.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. P57-P70 ◽  
Author(s):  
Shaun Strong ◽  
Steve Hearn

Survey design for converted-wave (PS) reflection is more complicated than for standard P-wave surveys, due to raypath asymmetry and increased possibility of phase distortion. Coal-scale PS surveys (depth [Formula: see text]) require particular consideration, partly due to the particular physical properties of the target (low density and low velocity). Finite-difference modeling provides a pragmatic evaluation of the likely distortion due to inclusion of postcritical reflections. If the offset range is carefully chosen, then it may be possible to incorporate high-amplitude postcritical reflections without seriously degrading the resolution in the stack. Offsets of up to three times target depth may in some cases be usable, with appropriate quality control at the data-processing stage. This means that the PS survey design may need to handle raypaths that are highly asymmetrical and that are very sensitive to assumed velocities. A 3D-PS design was used for a particular coal survey with the target in the depth range of 85–140 m. The objectives were acceptable fold balance between bins and relatively smooth distribution of offset and azimuth within bins. These parameters are relatively robust for the P-wave design, but much more sensitive for the case of PS. Reduction of the source density is more acceptable than reduction of the receiver density, particularly in terms of the offset-azimuth distribution. This is a fortuitous observation in that it improves the economics of a dynamite source, which is desirable for high-resolution coal-mine planning. The final-survey design necessarily allows for logistical and economic considerations, which implies some technical compromise. However, good fold, offset, and azimuth distributions are achieved across the survey area, yielding a data set suitable for meaningful analysis of P and S azimuthal anisotropy.


1997 ◽  
Vol 40 (1) ◽  
Author(s):  
B. Alessandrini ◽  
L. Beranzoli ◽  
G. Drakatos ◽  
C. Falcone ◽  
G. Karantonis ◽  
...  

We present a tomographic view of the crust and uppermost mantle beneath the Central Mediterranean area obtained from P-wave arrival times of regional earthquakes selected from the ISC bulletin. The P-wave velocity anomalies are obtained using Thurber's algorithm that jointly relocates earthquakes and computes velocity adjustments with respect to a starting model. A specific algorithm has been applied to achieve a distribution of epicentres as even as possible. A data set of 1009 events and 49072 Pg and Pn phases was selected. We find a low velocity belt in the crust, evident in the map view at 25 km of depth, beneath the Hellenic arc. A low velocity anomaly extends at 40 km of depth under the Aegean back arc basin. High velocities are present at Moho depth beneath the Ionian sea close to the Calabrian and Aegean arcs. The tomographic images suggest a close relationship between P-wave velocity pattern and the subduction systems of the studied area.


2021 ◽  
Vol 873 (1) ◽  
pp. 012064
Author(s):  
M Yasir ◽  
P T Brilianti ◽  
S S Angkasa ◽  
S Widyanti ◽  
I Herawati ◽  
...  

Abstract The tectonic setting of Java Island is mainly controlled by the collision of Indo-Australian plate subducting the Eurasian plate. The high collision activity of Eurasian and Indo-Australian plates often causes megathrust earthquakes and the rise of arc magmatism that includes volcanic eruption. This study aims to determine the tectonic pattern beneath Central Java based on P-wave tomography inversion. We used the fast-marching method as ray tracing and subspace inversion to image subsurface velocity model to a depth of 150 km. The data used in this study are catalogue events data derived from a temporary seismometer network MERAMEX installed around central Java and DOMERAPI installed surround Mt. Merapi and Mt. Merbabu. We also include events collected from the International Seismological Centre. In total, we processed 563 earthquake events to illustrate velocity structures under central Java. The checker-board model shows that good resolutions can be identified at shallow depth, including offshore south Java contributed from Ocean Bottom Seismometer data. In vertical axis, good resolution models can be expected down to a depth 150 km following rich events from the Benioff zone. Current P wave model show a distinct low velocity zone under Mt Merapi that can be seen down to a depth of 40 km, suggesting a possible separated deep magma reservoir. To the south of Mt Merapi area also shows a low-velocity band that may be related with the southern mountain arc. Additionally, the northern part of Mt. Merapi displays a band of strong low-velocity anomaly to the East and West with the anomaly in the Eastern Part seems to have a deeper extension to a depth of ~50 km. We related this anomaly with Merapi Lawu Anomaly and Kendeng basin. Our results show a similar result with the previous tomography models in this region.


1990 ◽  
Vol 80 (1) ◽  
pp. 129-139 ◽  
Author(s):  
Jose Pujol ◽  
Richard Aster

Abstract Arrival time data from the Phlegraean Fields (Italy) earthquake swarm recorded by the University of Wisconsin array in 1983 to 1984 were reanalyzed using a joint hypocentral determination (JHD) technique. The P- and S-wave station corrections computed as part of the JHD analysis show a circular pattern of central positive values surrounded by negative values whose magnitudes increase with distance from the center of the pattern. This center roughly coincides with the point of the maximum uplift (almost 2 m) associated with the swarm. Corrections range from −0.85 to 0.10 sec for P-wave arrivals and from −1.09 to 0.70 sec for S-wave arrivals. We interpret these patterns of corrections as caused by a localized low-velocity anomaly in the epicentral area, which agrees with the results of a previous 3-D velocity inversion of the same data set. The relocated (JHD) epicenters show less scatter than the epicenters obtained in the velocity inversion, and move more of the seismic activity to the vicinity of the only presently active fumarolic feature. The capability of the JHD technique to detect low-velocity anomalies and at the same time to give reliable locations, particularly epicenters, was verified using synthetic data generated for a 3-D velocity model roughly resembling the model obtained by velocity inversion.


2021 ◽  
Vol 873 (1) ◽  
pp. 012067
Author(s):  
Haolia ◽  
M. I. Sulaiman ◽  
P. T. Brilianti ◽  
R. P. Nugroho ◽  
I. Madrinovella ◽  
...  

Abstract The Sunda-Arc transition to the Banda Arc is located on the south of the Flores Island, Indonesia, where the Australian lithosphere is moving to the north direction. On-going subduction process dictates the tectonic setting though some studies also suggest a collision and obduction may occur in the past due to of plate buoyancy. This area has active seismicity with frequent large magnitude events. To better understand the tectonic system in this region, we performed double-difference tomography inversion using regional events. We obtained the data catalog from the Indonesian Agency of Meteorology, Climatology, and Geophysics ranging from 116° to 125° east longitude and -6.5° to 12.5° latitude. We collected 4312 events data, detected from 15 stations from January 2015 to December 2019. Final relocated hypocenters showed a reduced fixed-depth problem and a more clustered event, although some deep events disappear. Most events are related to the subducting Benioff zone with some clustered events in the northern area may be related to back-arc thrust. We also observed clustered events near active volcano region and reduced shallow seismicity region to the west of the Timor Island. Resolution test using the checkerboard and Derivative weigh Sum (DWS) shows that fair P wave resolution can be achieved until 300 km, although a smearing start to show at a deeper depth. However, due to lack of arrival S wave data, the resolution test suggest good resolution can only be seen until a depth of 100 km. Tomogram P and S wave models show a clear dipping subducting slab from south to North down to a 250 km. We also spot a fast velocity band near the Timor Island area that similar to the previous tomography study, interpreted as sliver forearm. We spotted a band of lower Vp, lower Vs and higher Vp/Vs at shallow depth close to the volcanic line and we interpreted this as a zone of higher temperature, that may relate to magmatic activity in this region. We also noticed a zone of low velocity and higher Vp/Vs that may relate with dehydration and partial melting. However, we feel this still uncertain due to low Vs resolution.


2021 ◽  
Author(s):  
Michael J. Schmidtke ◽  
Ruth Keppler ◽  
Jacek Kossak-Glowczewski ◽  
Nikolaus Froitzheim ◽  
Michael Stipp

Abstract. Subduction and exhumation are key processes in the formation of orogenic systems across the world, for example, in the European Alps. For geophysical investigations of these orogens, it is essential to understand the petrophysical properties of the rocks involved. These are the result of a complex interaction of mineral composition and rock fabric including mineral textures (i.e. crystallographic preferred orientations). In this study we present texture-derived elastic anisotropy data for a representative set of different lithologies involved in the Alpine orogeny. Rock samples were collected in the Lago di Cignana area in Valtournenche, in the Italian Northwestern Alps. At this locality a wide range of units of continental and oceanic origin with varying paleogeographic affiliations and tectono-metamorphic histories are accessible. Their mineral textures were determined by time-of-flight neutron diffraction at the Frank Laboratory of Neutron Physics at the JINR in Dubna, Russia. From these data the elastic properties of the samples were calculated. The data set includes representative lithologies from a subduction-exhumation-setting. In subducted lithologies originating from the oceanic crust, the elastic anisotropies range from 1.4 to 5.0 % with average P-wave velocities of 7.01–8.24 km/s and VP / VS-ratios of 1.71–1.76. In the metasediments of the former accretionary prism the elastic anisotropies range from 4.7 to 8.2 %. This tectonic setting displays average P-wave velocities of 6.47–7.23 km/s and VP / VS-ratios of 1.60–1.76. Continental crust which is incorporated in the collisional orogen shows elastic anisotropies ranging from 1.8 to 2.8 % with average P-wave velocities of 6.42–6.51 km/s and VP / VS-ratios of 1.56–1.60. Our results suggest that mafic and felsic rocks in subduction zones at depth may be discriminated by a combination of seismic signatures: lower anisotropy and higher VP / VS ratio for mafic rocks, higher anisotropy and lower VP / VS ratio for felsic rocks and metasediments.


2021 ◽  
Author(s):  
Ricardo Correia ◽  
Manel Prada ◽  
Valenti Sallares ◽  
Irene Merino ◽  
Alcinoe Calahorrano ◽  
...  

<p>The SW Iberian Margin has a complex tectonic setting and crustal structure derived from a succession of rift events related to the opening of North Atlantic and Neotethys from the Mesozoic to the Lower Cretaceous, and to the subsequent convergence between Nubian and Eurasian plates from Lower Oligocene to present day. Neogene plate convergence led to the reactivation of pre-existing extensional faults originated during the Mesozoic rifting in a combination of thrust and strike-slip systems. Despite a slow convergence rate, these faults now represent a major regional seismological and tsunamigenic hazard, as demonstrated by the devastating 1755 Lisbon earthquake of M>8.5 and the ensuing tsunami. Thus, unveiling the lithospheric structure along the SW Iberian Margin is not only important to understand the different stages of rifting and compression, but also to characterize the distribution of major lithospheric-scale boundaries, currently active and potentially capable of generating great, destructive tsunamigenic earthquakes.</p><p>To this end, we use here a spatially coincident wide-angle seismic (WAS) and multichannel seismic (MCS) data set collected along a NW-SE ~320 km transect SW of São Vicente cape during the FRAME survey in 2018. WAS data were recorded by 24 ocean bottom seismometers and hydrophones (OBS/H) while the MCS data were recorded by a 6 km long, 480 channel streamer. From NW to SE, the transect runs from the Tagus Abyssal plain to the westernmost extension of the Gulf of Cadiz across four major thrust faults, namely, the Tagus Abyssal Plain fault, Marquês de Pombal fault, São Vicente fault, and Horseshoe fault.</p><p>We applied joint refraction and reflection travel-time tomography (TTT) using a combination of arrival times identified at both WAS and MCS recordings to invert for the 2D P-wave velocity (Vp) structure of the crust and uppermost mantle, as well as the geometry of the main structural boundaries identified as seismic reflectors: the top of the acoustic basement and the Moho. Combining WAS and MCS travel-times brings a remarkable increase in the resolution and accuracy of the structure of the upper layers (i.e. top of the basement) thanks to the huge increase of spatial sampling in the shallow parts of the crust provided  by MCS data as compared to WAS data alone.</p><p>The inverted model shows a Vp structure with abrupt lateral velocity and structural variations marked by a rough Top of Basement topography and sharp changes in crustal thickness. In the northernmost part of the model there is evidence of mantle exhumation. The Moho shallows beneath the NE continuation of the Horseshoe Basin and the Gorringe Bank, coinciding with the location of the Marquês de Pombal, São Vicente, and Horseshoe thrust faults. The inversion of deep seismic phases reveals the presence of four southwest dipping reflectors that sheds new light into the deep geometry these major regional thrust faults.</p>


2021 ◽  
Vol 873 (1) ◽  
pp. 012057
Author(s):  
P T Brilianti ◽  
Haolia ◽  
M I Sulaiman ◽  
S S Angkasa ◽  
S Widyanti ◽  
...  

Abstract Our study area is located near island Sumbawa, Sumba, Flores, West Timor, Indonesia and East Timor, popularly known as Sunda-Banda arc transition zone. The tectonic setting is mainly controlled by the movement of the oceanic lithosphere Indo-Australian plate subducting the Eurasian plate and Northward migration of Australian continental lithosphere into western Banda-arc in the region of Flores, Sumba and Timor island. We tried to image velocity structure beneath these regions using regional events and tomography inversion model. We collected 5 years of regional events from the Indonesian Agency of Meteorology, Climatology and Geophysics. In total, we reserved 3186 events recorded on 29 stations. For data processing, we used fast marching method as ray tracing between sources and receiver. We then employed subspace inversion as the tomography procedure to estimate the best velocity model representing the tectonic model in the region. Hypocenter data distribution is concentrated on shallow parts of the region and along the Benioff zone down to a maximum depth of 400 km. One of challenge of this study is that although events are abundance, the stations used are mostly located onshore and does not extend in the south-north direction that leads us to under determined problem in the inversion process. However, checker-board models show most our target area can be retrieved to its initial model with sign of smearing effects shown start from a depth of 50 km. After six iteration and optimized selection of damping and smoothing parameters, we observed low velocity anomaly under Bali, Lombok, Sumba, East Nusa Tenggara at shallow depth that may be related with volcanic activity. Deeper low anomaly can also be seen that may be related with partial melting process. A band of fast velocity is clearly seen that goes deepen to the north depicting subducting slabs own to a depth of 300 km. We also observed a possible of fast velocity in the northern part of our stations at shallow depth that we believe may represent the back arc thrust.


Sign in / Sign up

Export Citation Format

Share Document