Subduction zone Calibration and teleseismic relocation of thrust zone events in the central Aleutian Islands

1981 ◽  
Vol 71 (6) ◽  
pp. 1805-1828
Author(s):  
Kazuya Fujita ◽  
E. R. Engdahl ◽  
Norman H. Sleep

Abstract Teleseismically determined epicenters for thrust zone earthquakes in the central Aleutian Islands are systematically mislocated by about 40 km to the north relative to local network solutions. We attribute the mislocation to travel-time anomalies produced by the higher P-wave velocity in the subducting slab. These travel-time anomalies can be estimated by seismic ray tracing through a thermally modeled velocity structure for the slab. Observed teleseismic station residuals and slab travel-time anomalies computed from local network hypocenters, plotted as a function of event distance to the center of curvature of the island arc, show good agreement. This distance versus residual relationship can be used to validate the local network hypocenters and to constrain the length of subducting lithosphere. Computed slab source corrections, when incorporated in a standard location program, reduce teleseismic mislocations to 20 km; the addition of empirical station corrections reduces the mislocation to about 10 to 15 km. Depth and origin time estimates can be further improved through use of better station corrections, a local station within 2°, and depth phases.

2021 ◽  
Vol 13 (13) ◽  
pp. 2449
Author(s):  
Huiyan Shi ◽  
Tonglin Li ◽  
Rui Sun ◽  
Gongbo Zhang ◽  
Rongzhe Zhang ◽  
...  

In this paper, we present a high resolution 3-D tomographic model of the upper mantle obtained from a large number of teleseismic travel time data from the ISC in the central Philippines. There are 2921 teleseismic events and 32,224 useful relative travel time residuals picked to compute the velocity structure in the upper mantle, which was recorded by 87 receivers and satisfied the requirements of teleseismic tomography. Crustal correction was conducted to these data before inversion. The fast-marching method (FMM) and a subspace method were adopted in the forward step and inversion step, respectively. The present tomographic model clearly images steeply subducting high velocity anomalies along the Manila trench in the South China Sea (SCS), which reveals a gradual changing of the subduction angle and a gradual shallowing of the subduction depth from the north to the south. It is speculated that the change in its subduction depth and angle indicates the cessation of the SCS spreading from the north to the south, which also implies that the northern part of the SCS opened earlier than the southern part. Subduction of the Philippine Sea (PS) plate is exhibited between 14° N and 9° N, with its subduction direction changing from westward to eastward near 13° N. In the range of 11° N–9° N, the subduction of the Sulu Sea (SS) lies on the west side of PS plate. It is notable that obvious high velocity anomalies are imaged in the mantle transition zone (MTZ) between 14° N and 9° N, which are identified as the proto-SCS (PSCS) slabs and paleo-Pacific (PP) plate. It extends the location of the paleo-suture of PSCS-PP eastward from Borneo to the Philippines, which should be considered in studying the mechanism of the SCS and the tectonic evolution in SE Asia.


1998 ◽  
Vol 41 (1) ◽  
Author(s):  
J. Plomerová ◽  
V. Babuska ◽  
R. Scarpa

Jeffreys-Bullen (absolute) and relative P-wave travel-time residuals were analyzed over Italy and its surrounding using P arrival times from the ISC bulletins supplemented by the data from local observatories. We analyzed the travel-time station corrections by two independent methods to obtain information on lateral variations of the velocity structure over the area and a view of possible upper mantle anisotropy. In the first method, the station corrections are computed as a constant and two cosine terms with appropriate phase shifts. Besides a static term, the second method allows us to study the relative residuals in dependence both on azimuths and incidence angles and thus to investigate their spatial variations and to map lateral variations of anisotropic structure of the subcrustal lithosphere. The high and low-velocity directions inferred from the spatial distribution of the relative residuals as well as the high- and low-velocity upper mantle heterogeneities reflect the geodynamic development of the region, governed by the collision between the African and Eurasian plates


Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. R45-R56 ◽  
Author(s):  
Lars Nielsen ◽  
Hans Thybo ◽  
Martin Glendrup

Seismic wide-angle data were recorded to more than 300-km offset from powerful airgun sources during the MONA LISA experiments in 1993 and 1995 to determine the seismic-velocity structure of the crust and uppermost mantle along three lines in the southeastern North Sea with a total length of 850 km. We use the first arrivals observed out to an offset of 90 km to obtain high-resolution models of the velocity structure of the sedimentary layers and the upper part of the crystalline crust. Seismic tomographic traveltime inversion reveals 2–8-km-thick Paleozoic sedimentary sequences with P-wave velocities of 4.5–5.2 km/s. These sedimentary rocks are situated below a Mesozoic-Cenozoic sequence with variable thickness: ∼2–3 km on the basement highs, ∼2–4 km in the Horn Graben and the North German Basin, and ∼6–7 km in the Central Graben. The thicknesses of the Paleozoic sedimentary sequences are ∼3–5 km in the Central Graben, more than 4 km in the Horn Graben, up to ∼4 km on the basement highs, and up to 8 km in the North German Basin. The Paleozoic strata are clearly separated from the shallower and younger sequences with velocities of ∼1.8–3.8 km/s and the deeper crystalline crust with velocities of more than 5.8–6.0 km/s in the tomographic P-wave velocity model. Resolution tests show that the existence of the Paleozoic sediments is well constrained by the data. Hence, our wide-angle seismic models document the presence of Paleozoic sediments throughout the southeastern North Sea, both in the graben structures and in deep basins on the basement highs.


2019 ◽  
Vol 110 (1) ◽  
pp. 26-37 ◽  
Author(s):  
Masumi Yamada ◽  
Thakur Kandel ◽  
Koji Tamaribuchi ◽  
Abhijit Ghosh

ABSTRACT In this article, we created a well-resolved aftershock catalog for the 2015 Gorkha earthquake in Nepal by processing 11 months of continuous data using an automatic onset and hypocenter determination procedure. Aftershocks were detected by the NAMASTE temporary seismic network that is densely distributed covering the rupture area and became fully operational about 50 days after the mainshock. The catalog was refined using a joint hypocenter determination technique and an optimal 1D velocity model with station correction factors determined simultaneously. We found around 15,000 aftershocks with the magnitude of completeness of ML 2. Our catalog shows that there are two large aftershock clusters along the north side of the Gorkha–Pokhara anticlinorium and smaller shallow aftershock clusters in the south. The patterns of aftershock distribution in the northern and southern clusters reflect the complex geometry of the Main Himalayan thrust. The aftershocks are located both on the slip surface and through the entire hanging wall. The 1D velocity structure obtained from this study is almost constant at a P-wave velocity (VP) of 6.0  km/s for a depth of 0–20 km, similar to VP of the shallow continental crust.


1967 ◽  
Vol 57 (5) ◽  
pp. 1017-1023 ◽  
Author(s):  
Max Wyss ◽  
James N. Brune

abstract The seismograms of the Alaskan earthquake of 28 March 1964 are characterized by multiple P-phases not predicted by the travel-time curves. Seismograms with low magnifications from 80 stations covering distances from 40° to 90° and a wide range of azimuths were analyzed. The character of the P-wave portion of the seismograms is interpreted in terms of an approximate multiple-event source mechanism where the propagating rupture triggers larger distinct events. Six events were located using the Gutenberg sine-curve method. The times after the initial origin time were 9, 19, 28, 29, 44 and 72 sec respectively, and the events were located 35, 66, 89, 93, 165 and 250 km away from the initial epicenter. Dividing the distance by the delay-time gives an average rupture velocity of 3.5 km/sec.


1990 ◽  
Vol 80 (2) ◽  
pp. 395-410 ◽  
Author(s):  
Glenn D. Nelson ◽  
John E. Vidale

Abstract We present a new method for locating earthquakes in a region with arbitrarily complex three-dimensional velocity structure, called QUAKE3D. Our method searches a gridded volume and finds the global minimum travel-time residual location within the volume. Any minimization criterion may be employed. The L1 criterion, which minimizes the sum of the absolute values of travel-time residuals, is especially useful when the station coverage is sparse and is more robust than the L2 criterion (which minimizes the RMS sum) employed by most earthquake location programs. On a UNIX workstation with 8 Mbytes memory, travel-time grids of size 150 by 150 by 50 are reasonably employed, with the actual geographic coverage dependent on the grid spacing. Location precision is finer than the grid spacing. Earthquake recordings at six stations in Bear Valley are located as an example, using various layered and laterally varying velocity models. Locations with QUAKE3D are nearly identical to HYPOINVERSE locations when the same flat-layered velocity model is used. For the examples presented, the computation time per event is approximately 4 times slower than HYPOINVERSE, but the computation time for QUAKE3D is dependent only on the grid size and number of stations, and independent of the velocity model complexity. Using QUAKE3D with a laterally varying velocity model results in locations that are physically more plausible and statistically more precise. Compared to flat-layered solutions, the earthquakes are more closely aligned with the surface fault trace, are more uniform in depth distribution, and the event and station travel-time residuals are much smaller. Hypocentral error bars computed by QUAKE3D are more realistic in that the trade-off of depth versus origin time is implicit in our error estimation, but ignored by HYPOINVERSE.


2020 ◽  
Author(s):  
Hanneke Paulssen ◽  
Wen Zhou

<p>Between 2013 and 2017, the Groningen gas field was monitored by several deployments of an array of geophones in a deep borehole at reservoir level (3 km). Zhou & Paulssen (2017) showed that the P- and S-velocity structure of the reservoir could be retrieved from noise interferometry by cross-correlation. Here we show that deconvolution interferometry of high-frequency train signals from a nearby railroad not only allows determination of the velocity structure with higher accuracy, but also enables time-lapse measurements. We found that the travel times within the reservoir decrease by a few tens of microseconds for two 5-month periods. The observed travel time decreases are associated to velocity increases caused by compaction of the reservoir. However, the uncertainties are relatively large. <br>Striking is the large P-wave travel time anomaly (-0.8 ms) during a distinct period of time (17 Jul - 2 Sep 2015). It is only observed for inter-geophone paths that cross the gas-water contact (GWC) of the reservoir. The anomaly started 4 days after drilling into the reservoir of a new well at 4.5 km distance and ended 4 days after the drilling operations stopped. We did not find an associated S-wave travel time anomaly. This suggests that the anomaly is caused by a temporary elevation of the GWC (water replacing gas) of approximately 20 m. We suggest that the GWC is elevated due to pore-pressure variations during drilling. The 4-day delay corresponds to a pore-pressure diffusivity of ~5m<sup>2</sup>/s, which is in good agreement with the value found from material parameters and the diffusivity of (induced) seismicity for various regions in the world. </p>


2021 ◽  
Author(s):  
Gregor Rajh ◽  
Josip Stipčević ◽  
Mladen Živčić ◽  
Marijan Herak ◽  
Andrej Gosar

<p>The investigated area of the NW Dinarides is bordered by the Adriatic foreland, the Southern Alps, and the Pannonian basin at the NE corner of the Adriatic Sea. Its complex crustal structure is the result of interactions among different tectonic units. Despite numerous seismic studies taking place in this region, there still exists a need for a detailed, smaller scale study focusing mainly on the brittle part of the Earth's crust. Therefore, we decided to investigate the velocity structure of the crust using concepts of local earthquake tomography (LET) and minimum 1-D velocity model. Here, we present the results of the 1-D velocity modeling and the catalogue of the relocated seismicity. A minimum 1-D velocity model is computed by simultaneous inversion for hypocentral and velocity parameters together with seismic station corrections and represents the best fit to the observed arrival times.</p><p>We used 15,579 routinely picked P wave arrival times from 631 well-located earthquakes that occurred in Slovenia and in its immediate surroundings (mainly NW Croatia). Various initial 1-D velocity models, differing in velocity and layering, were used as input for velocity inversion in the VELEST program. We also varied several inversion parameters during the inversion runs. Most of the computed 1-D velocity models converged to a stable solution in the depth range between 0 and 25 km. We evaluated the inversion results using rigorous testing procedures and selected two best performing velocity models. Each of these models will be used independently as the initial model in the simultaneous hypocenter-velocity inversion for a 3-D velocity structure in LET. Based on the results of the 1-D velocity modeling, seismicity distribution, and tectonics, we divided the study area into three parts, redefined the earthquake-station geometry, and performed the inversion for each part separately. This way, we gained a better insight into the shallow velocity structure of each subregion and were able to demonstrate the differences among them.</p><p>Besides general structural implications and a potential to improve the results of LET, the new 1-D velocity models along with station corrections can also be used in fast routine earthquake location and to detect systematic travel time errors in seismological bulletins, as already shown by some studies using similar methods.</p>


Sign in / Sign up

Export Citation Format

Share Document