Elastic Parameters of West Bohemian Granites under Hydrostatic Pressure

—The West Bohemian seismoactive region is situated near the contact of the Moldanu bian, Bohemian and Saxothuringian units in which a large volume is occupied by granitoid massifs. The spatial distribution of P-wave velocities and the rock fabric of five representative samples from these massifs were studied. The P-wave velocities were measured on spherical samples in 132 independent directions under hydrostatic pressure up to 400 MPa, using the pulse-transmission method. The pressure of 400 MPa corresponds to a depth of about 15 km in the area under study. The changes of P-wave velocity were correlated with the preferred orientations of the main rock fabric elements, i.e., rock forming minerals and microcracks. The values of the P-wave velocity from laboratory measurements on granite samples fit the velocity model used by seismologists in the West Bohemian seismoactive region.     

Seismic Anisotropy and Velocity Variations in the Mantle beneath the Saxothuringicum-Moldanubicum Contact in Central Europe

—We report on results of a passive seismic experiment undertaken to study the 3-D velocity structure and anisotropy of the upper mantle around the contact zone of the Saxothuringicum and Moldanubicum in the western margin of the Bohemian Massif in central Europe. Spatial variations of P-wave velocities and lateral variations of the particle motion of split shear waves over the region monitor changes of structure and anisotropy within the deep lithosphere and the asthenosphere. A joint interpretation of P-residual spheres and shear-wave splitting results in an anisotropic model of the lithosphere with high velocities plunging divergently from the contact of both tectonic units. Lateral variations of the mean residuals are related to a southward thickening of the lithosphere beneath the Moldanubicum.     

Three-dimensional P-wave Velocity Structure of the Bear Valley Region of Central California

—The three-dimensional P-wave velocity structure of the Bear Valley region of central California is determined by applying a circular ray-tracing technique to 1735 P-wave arrivals from 108 locally recorded earthquakes. Comparison of the results obtained from one-dimensional and laterally varying starting models shows that many of the features in the structure determined are fairly insensitive to the choice of the starting model. Velocities associated with the Gabilan granites southwest of the San Andreas Fault are slightly higher than those in the Franciscan formation to the northeast, and these two features are separated in the southern part of the region by a narrow fault zone with very low velocities. In the southeastern part of the region, where the Gabilan granites do not abut the San Andreas Fault, the low velocities of the fault zone cross over to the southwestern side of the fault. They also appear to extend to depths of at least 15km, thus locally reversing the contrast across the San Andreas Fault that prevails farther to the northwest. In the northwestern part of the region, the low velocities of the fault zone split and follow the surface traces of the San Andreas and Calaveras Faults, but do not appear to extend to depths much deeper than about 6km. There also appears to be a well-defined contrast in structure in the middle of the Santa Clara Valley, suggesting the existence of a fault in the basement of the valley that may be a southern extension of the Sargent Fault into this region. Relocated hypocenters beneath the San Andreas Fault cluster in a zone that dips about 80° southwest and intersects the surface trace of the fault in the southern part of region.     

Temporal Variations of Seismic Wave Velocity Associated with 1998 M6.1 Shizukuishi Earthquake

—?We attempt to detect temporal variations of seismic wave velocity before and after 1998 M6.1 Shizukuishi, northeastern Japan, earthquake by using waveform data from explosions and earthquake doublets spanning the period immediately before and after the earthquake. Direct P waves of the second explosion are delayed by ～20 ms at observation stations with epicentral distances less than 15 km. This tendency does not change if the analysis frequency band is changed. Our result suggests that the P-wave velocity decreased by at least 1% in the extremely shallow region of the hanging wall of the M6.1 thrust event after its occurrence. On the other hand, there was the frequency dependence of the coda wave delays for both artificial sources and for natural events. At 5–10 Hz, immediate sharp increases by more than 20 ms in time delays and lower coherency were observed at several stations. We estimated the region in which P-wave velocity might have decreased after the M6.1 earthquake. Maximum depth of the region is 13 km. The region includes the aftershock area of the M6.1 earthquake, but is eccentric to the earthquake in the west. Considering the frequency band analyzed (5–10 Hz), the scale of the spatial inhomogeneity which led to the coda wave delay is several hundreds meters. We investigated candidates for the cause of the direct P-wave and coda wave delay. Observed direct P-wave delay can be partly explained by the stress changes caused by coseismic fault slip. However, the coda wave delay cannot be explained by the stress changes that are limited to the superficial area. Crustal heterogeneity should have changed at coseismic time in the deeper area where aftershocks of the M6.1 earthquake occurred.     

Crustal structure of Tushka Region, Abu-Simbel, Egypt, inferred from spectral ratios of P waves of local earthquakes

The crustal structure of North Abu-Simbel area was studied using spectral ratios of short-period P waves. Three-component short period seismograms from the Masmas seismic station of the Egyptian National Seismic Network Stations were used. The Thomson-Haskell matrix formulation was applied for linearly elastic, homogeneous crustal layers. The obtained model suggests that the crust under the study region consists of a thin (0.8 km) superficial top layer with a P-wave velocity of 3.8±0.7 km/s and three distinct layers with a mean P-wave velocity of 6.6 km/s, overlaying the upper mantle with a P-wave velocity of 8.3 km/s (fixed). The results were obtained for 14 different earthquakes. The P-wave velocities of the three layers are: 5.8±0.6 km/s, 6.5±0.4 km/s and 7.2±0.3 km/s. The total depth to the Moho interface is 32±2 km. The crustal velocity model estimated using observations is relatively simple, being characterized by smooth velocity variations through the middle and lower crust and normal crustal thickness. The resultant crustal model is consistent with the model obtained from previous deep seismic soundings along the northern part of Aswan lake zone.     

Tomographic inversion for the three-dimensional seismic velocity structure of the Aswan region,Egypt

The Thurber iterative simultaneous inversion program is used to determine the three-dimensionalP-wave velocity structure in the Aswan seismic region of Egypt. The tomographic inversion presented in this study is based on 1131P-phase observations at 13 stations from 89 local earthquakes, all of which occurred within the Kalabsha fault zone. The assumed initial velocity model is that deduced from local explosion experiments. The results indicate that the Aswan region is characterized by a heterogeneous crust, consisting of a shallow, low-velocity zone and a deeper high-velocity anomaly. Seismic velocity structure within the shallow part demonstrates that the inferred change in velocity exists primarily across the east-west trending Kalabsha fault scarp, whereas the high-velocity zone is located south of this fault. Two well-resolved, low-velocity zones appear within the upper 6 km of the crust. The first coincides with a graben structure located between the Kalabsha and Seiyal faults and the second exists between the N-S Kurkur fault and the main axis of Lake Aswan. Both low-velocity zones occupy an area of approximately 30×40 km, located along the western bank of the lake. The most significant result of this study is that the location of the deeper, high-velocity anomaly coincides with the concentration of seismic activity in the lower crustal layer.     

A study of the velocity structure for the Benioff zone of Kamchatka: The Avacha Bay-Kronotskii Bay Segment

We present computations of 3D P-wave velocity field for the segment of the Benioff zone of Kamchatka between 51.5° N to Avacha Bay to Cape Shipunskii to Kronotskii Bay. P-wave travel times for regional earthquakes were used to compute and construct 3D models of the velocity field for this segment of the Benioff zone. We examined the velocity structure to interpret it in conjunction with neotectonic morphostructures. The computations showed a complex structure of the field and the presence of inhomogeneities, both along and across the Benioff zone. Of special interest are the results for the Avacha and Kronotskii basins, where high velocity masses are unusually shallow beneath the low velocity layers in the top of the earth section. This creates high-gradient zones that can possibly generate large earthquakes.     

Global Event Location with Full and Sparse Data Sets Using Three-dimensional Models of Mantle P-wave Velocity

—?In order to improve on the accuracy of event locations at teleseismic distances it is necessary to adequately correct for lateral variations in structure along the ray paths, either through deterministic model-based corrections, empirical path/station corrections, or a combination of both approaches. In this paper we investigate the ability of current three-dimensional models of mantle P-wave velocity to accurately locate teleseismic events. We test four recently published models; two are parameterized in terms of relatively long-wavelength spherical harmonic functions up to degree 12, and two are parameterized in terms of blocks of constant velocity which have a dimension of a few hundreds of km. These models, together with detailed crustal corrections, are used to locate a set of 112 global test events, consisting of both earthquakes and explosions with P-wave travel-time data compiled by the Internation al Seismological Centre (ISC). The results indicate that the supposedly higher resolution block models do not improve the accuracy of teleseismic event locations over the longer wavelength spherical harmonic models. For some source locations the block models do not predict the range of observed travel-time residuals as well as the longer wavelength models. The accuracy of the locations largely varies randomly with geographic position although events in central Asia are particularly well located. We also tested the effect of reduced data sets on the locations. Multiple location iterations using 30 P-wave travel times indicate that teleseismic events may be located within an area of 1000?km² of the true location 66% of the time with only the model-based corrections, and increasing to 75% if calibration information is available. If as few as 8 phases are available then this is possible only 50% of the time. Further refinement in models and/or procedure, such as the addition of P _n phases, azimuth data, and consideration of P-wave anisotropy may provide further improvement in the teleseismic location of small events.     

Simultaneous inversion of velocity structures and hypocentral locations: Application to the Friuli seismic area NE Italy

A layeredP- andS-wave velocity model is obtained for the Friuli seismic area using the arrival time data ofP- andS-waves from local earthquakes. A damped least-squares method is applied in the inversion.The data used are 994P-wave arrival times for 177 events which have epicenters in the region covered by the Friuli seismic network operated by Osservatorio Geofisico sperimentale (OGS) di Trieste, which are jointly inverted for the earthquake hypocenters andP-wave velocity model. TheS-wave velocity model is estimated on the basis of 978S-wave arrival times and the hypocenters obtained from theP-wave arrival time inversion. We also applied an approach thatP- andS-wave arrival time data are jointly used in the inversion (Roecker, 1982). The results show thatS-wave velocity structures obtained from the two methods are quite consistent, butP-wave velocity structures have obvious differences. This is apparent becauseP-waves are more sensitive to the hypocentral location thanS-waves, and the reading errors ofS-wave arrival times, which are much larger than those ofP-waves, bring large location errors in the joint inversion ofP- andS-wave arrival time. The synthetic data tests indicated that when the reading errors ofS-wave arrivals are larger than four times that ofP-wave arrivals, the method proposed in this paper seems more valid thanP- andS-wave data joint inversion. Most of the relocated events occurred in the depth range between 7 and 11 km, just above the biggest jump in velocity. This jump might be related to the detachment line hypothesized byCarulli et al. (1982). From the invertedP- andS-wave velocities, we obtain an average value 1.82 forV _p /V _s in the first 16 km depth.     

Velocity and Density Heterogeneities of the Tien-Shan Lithosphere

—The Tien-Shan orogene is a region in which the earth’s crust undergoes considerable thickening and tangential compression. Under these conditions the lithosphere heterogeneities (composi tion, rheological) create the prerequisites for the development of various phenomena of tectonic layering (lateral shearing, different deformation of layers). To study the distribution of velocity, density and other elastic parameters, the results from a seismic tomography study on P-wave as well as S-wave velocities were used. Using empirical as well as theoretical formulas on the relationship between velocity, density and silica content in rocks, their distribution in the Tien-Shan’s lithosphere has been calculated. In addition, other elastic parameters, such as Young’s modulus, shear modulus, Poisson’s ratio and coefficient of general compressions have been determined. Zoning of different types of crust was carried out for the region investigated. The characteristics of the "crust-mantle" transition have been investi gated. Large blocks with different types of the earth’s crust were distinguished. Layers with inverse values of velocity, density and shear and Young modulus are revealed in the Tien-Shan lithosphere. All of the above described features open new ways to solve geodynamics problems.     

Application of Multichannel Analysis of Surface Waves to S-Phase Wave Anisotropy Estimation

The Multichannel Analysis of Surface Waves (MASW) is an increasingly used technique for recognition of a shallow geological structure and estimation of geotechnical parameters, e.g., S-wave velocity, layer density, layer thickness, shear modulus, estimated P-wave velocity, and estimated Poisson ratio. MASW surveys were carried out in two limestone quarries in the southern part of Poland. The experimental areas are characterised by a simple geological structure: consolidated Triassic limestone. Measurement profiles were arranged as a shapely six-pointed star. For each survey line, 12 geophones with 2-meter (Deposit 1) and 3-meter (Deposit 2) spacing were applied. The research allowed to compare P- and S-wave velocity changes with the main crack systems in the studied rock masses.     

An attenuating structure beneath the Aso Caldera determined from the propagation of seismic waves

The attenuation of amplitude is seen in seismic waves which pass through the central region of the Aso caldera, in Kyushu, Japan. It is also recognized from spectral analysis of seismic waves that the higher frequencies of the P-wave are reduced in the waves which pass through the central region of the caldera. It is shown that the relative attenuation increases remarkably for the frequency range of 5 to 10 Hz. The specific attenuation factor Q of the P-wave train is about 100. From the surface projection of the ray paths with low Q values through the Aso caldera to each station, the attenuating region is located beneath the center of the caldera, extending to the north of the central cones. In conjunction with the low Q value of the P-wave and the decreases of S-wave amplitudes, the relative P-wave residual times have comparatively large values for seismic waves passing through the central region beneath the caldera. In order to attempt to provide additional information on the depth configuration of the attenuating material, the ray paths of P-wave's first arrivals are located in three-dimensional space. It indicates that the low-velocity material is located beneath the center of the caldera at depths of about 6 to 9 km. However, lowvelocity anomalies above the depth of 6 km and below the depth of 15 km were not able to be detected, because most of the available seismic ray paths had crossed the caldera at depths of about 6 to 15 km. Furthermore, the relative residual times have numerous errors resulting from incorrect hypocenter locations, origin times, inhomogeneities in the structure and uncertainty of the velocity structure. At shallow depths in the Aso caldera, refraction or reflection studies are required for an accurate estimate of the structure and more detailed properties of the attenuating material.     

Seismic Location Calibration of the European Arctic

—?A crustal velocity model has been developed for Fennoscandia, the Baltic shield and adjacent areas. This model represents a simplified average of various models developed for parts of this region. We show that P-wave travel times calculated with this model provide an excellent fit to observations at the Fennoscandian, KRSC and IRIS station networks for a set of seismic events with known or very well-constrained locations. The station-event paths cover large parts of Western Russia and the Barents Sea, thus indicating that this model, which we denote the Barents model, is appropriate for this entire region. We show by examples that significant improvements in event location precision can be achieved compared to using the IASPEI model. We finally use the Barents model to calculate locations of recent small seismic events in the Novaya Zemlya region of interest in a CTBT monitoring context.     

Constraining P-wave velocity variations in the upper mantle beneath Southeast Asia

We have produced a P-wave model of the upper mantle beneath Southeast (SE) Asia from reprocessed short period International Seismological Centre (ISC) P and pP data, short period P data of the Annual Bulletin of Chinese Earthquakes (ABCE), and long period PP-P data. We used 3D sensitivity kernels to combine the datasets, and mantle structure was parameterized with an irregular grid. In the best-sampled region our data resolve structure on scale lengths less than 150 km. The smearing of crustal anomalies to larger depths is reduced by a crustal correction using an a priori 3D model. Our tomographic inversions reveal high-velocity roots beneath the Archean Ordos Plateau, the Sichuan Basin, and other continental blocks in SE Asia. Beneath the Himalayan Block we detect high seismic velocities, which we associate with subduction of Indian lithospheric mantle. This structure is visible above the 410 km discontinuity and may not connect to the remnant of the Neo-Tethys oceanic slab in the lower mantle. Our images suggest that only the southwestern part of the Tibetan plateau is underlain by Indian lithosphere and, thus, that the upper mantle beneath northeastern Tibet is primarily of Asian origin. Our imaging also reveals a large-scale high-velocity structure in the transition zone beneath the Yangtze Craton, which could have been produced in multiple subduction episodes. The low P-wave velocities beneath the Hainan Island are most prominent in the upper mantle and transition zone; they may represent counter flow from the surrounding subduction zones, and may not be unrelated to processes beneath eastern Tibet.     

Shear wave velocity modelling in crustal rock for seismic hazard analysis

P-wave velocity data along with the thickness of sedimentary and crystalline layers within bedrock were collected from all global regions and presented in the Global Crustal Model CRUST2.0, published in 2001. This well-organised database provides invaluable potential contributions towards future seismic hazard modelling, particularly for stable continental regions (SCRs), where there is a scarcity of representative strong motion records for conventional modelling purposes. The P-wave velocity information presented in CRUST2.0 has been converted herein to S-wave velocity information. The latter is especially important for purposes of seismic hazard modelling. The value of the CRUST2.0 model has therefore been greatly enhanced by the important findings presented and further developed in this paper. By making the best use of available information on crustal conditions, the amplification behaviour of seismic waves affecting a region, an area or a site for any given earthquake scenario may be predicted. The developed methodology, which is intended for worldwide applications, has been illustrated by case studies in which model S-wave velocity profiles were developed for different geological regions within North America. The model profiles were found to be in excellent agreement with field measurements reported for each respective region.     

Application of first-arrival tomography to characterize a quick clay landslide site in Southwest Sweden

First-arrival traveltime tomography was applied to high-resolution seismic data acquired over a known quick-clay landslide scar near the Göta River in southwest Sweden in order to reveal the geometry and physical properties of clay-related normally consolidated sediments. Investigated area proved to be a challenging environment for tomographic imaging because of large P-wave velocity variations, ranging from 500 to 6000 m/s, and relatively steeply-dipping bedrock. Despite these challenges, P-wave velocity models were obtained down to ca. 150 m for two key 2D seismic profiles (each about 500-m long) intersecting over the landslide scar. The models portrait the sandwich-like structure of marine clays and coarse-grained consolidated sediments, but the estimated resolution (20 m) is too small to distinguish thin layers within this structure. Modelled velocity structures match well the results of reflection seismic processing and resistivity tomography available along the same profiles.     

Seismic velocity and Poisson’s ratio tomography of the crust beneath southwest Anatolia: an insight into the occurrence of large earthquakes

The western part of Anatolia is one of the most seismically and tectonically active continental regions in the world, and much of it has been undergoing NS-directed extensional deformation since the Early Miocene. In this study, we determine 3-D tomographic images of the crust under the southwestern part of the North Anatolian Fault Zone by inverting a large number of arrival time data of P and S waves. From the obtained P- and S-wave velocity models, we estimated the Poisson’s ratio structures for a more reliable interpretation of the obtained anomalies. Our tomographic results confirmed the major tectonic features detected by previous studies and revealed new structural heterogeneities related to the active seismotectonics of the studied area. High P-wave velocity anomalies are recognized near the surface, while at deeper crustal layers, low P-wave velocities are widely distributed. The crustal S-wave velocity and Poisson’s ratio exhibit more structural heterogeneities compared to the P-wave velocity structure. Microearthquake activity is intense along highly heterogeneous zones in the southwestern part, which is characterized by low to high P-wave velocity, low S-wave velocity, and high Poisson’s ratio anomalies. Large earthquakes are also concentrated in zones dominated by low velocities and low to high Poisson’s ratios. Results of the checkerboard and synthetic tests indicate that the imaged anomalies are reliable features down to a depth of 25 km. Moreover, they are consistent with many geological and geophysical results obtained by other researchers along the southwestern part of the North Anatolian Fault Zone. An erratum to this article can be found at     

S-wave Velocity Models in the Scotia Sea Region, Antarctica, from Nonlinear Inversion of Rayleigh Waves Dispersion

—More than 60 events recorded by four recently deployed seismic broadband stations around Scotia Sea, Antarctica, have been collected and processed to obtain a general overview of the crust and upper mantle seismic velocities.¶Group velocity of the fundamental mode of Rayleigh waves in the period between 10 s to 30–40 s is used to obtain the S-wave velocity versus depth along ten different paths crossing the Scotia Sea region. Data recorded by two IRIS (Incorporated Research Institutions for Seismology) stations (PMSA, EFI) and the two stations of the OGS-IAA (Osservatorio Geofisico Sperimentale—Instituto Antarctico Argentino) network (ESPZ, USHU) are used.¶The Frequency-Time Analysis (FTAN) technique is applied to the data set to measure the dispersion properties. A nonlinear inversion procedure, "Hedgehog," is performed to retrieve the S-wave velocity models consistent with the dispersion data.¶The average Moho depth variation on a section North to South is consistent with the topography, geological observations and Scotia Sea tectonic models.¶North Scotia Ridge and South Scotia Ridge models are characterised by similar S-wave velocities ranging between 2.0 km/s at the surface to 3.2 km/s to depths of 8 km/s. In the lower crust the S-wave velocity increases slowly to reach a value of 3.8 km/s. The average Moho depth is estimated between 17 km to 20 km and 16 km to 19 km, respectively, for the North Scotia Ridge and South Scotia Ridge, while the Scotia Sea, bounded by the two ridges, has a faster and thinner crust, with an average Moho depth between 9 km and 12 km.¶On other paths crossing from east to west the southern part of the Scotia plate and the Antarctic plate south of South Scotia Ridge, we observe an average Moho depth between 14 km and 18 km and a very fast upper crust, compared to that of the ridge. The S-wave velocity ranges between 3.0 and 3.6 km/s in the thin (9–13 km) and fast crust of the Drake Passage channel. In contrast the models for the tip of the Antarctic Peninsula consist of two layers with a large velocity gradient (2.3–3.0 km/s) in the upper crust (6-km thick) and a small velocity gradient (3.0–4.0) in the lower crust (14-km thick).     

Injection-induced Microseismicity in Colorado Shales

—Imperial Oil Resources Limited uses cyclic steam stimulation to recover oil from their Cold Lake oil field in Alberta. This operation, in particular situations, can be associated with the failure of well casings in the Colorado shales above the oil-bearing formation. A number of fluid injection operations was undertaken at this site and the associated microseismicity was detected using two three-component geophones and fifteen hydrophones. The purpose of this experiment was to simulate the occurrence of a casing failure, determine the feasibility of monitoring in a shallow environment, and characterize the microseismic activity. A calibration survey provided values of 1786 ± 108 m/s for P-wave velocity, 643 ± 56 m/s for S-wave velocity and 0.428 ± 0.017 for Poisson’s ratio in the shale formation. Estimates of the quality factor Q _P were 15 for the horizontal direction and 38 for the vertical direction, corroborating the evidence of velocity anisotropy. Calibration shots were located to within 10 m of the actual shot points using triangulation and polarization techniques. Several hundred microseis mic events were recorded and 135 events were located. The results showed that microseismic activity was confined to depths within 10 meters of the injection depth. The experiment clearly established the feasibility of detecting microseismicity induced by fluid injection rates typical of casing failures in shales at distances over 100 m.     

Subcrustal structure of the black sea basin from seismological data

The P-wave travel time data from the earthquakes offshore and onshore around the Black Sea are used for the tomographic reconstruction of the three-dimensional (3D) velocity distribution in the lithosphere of the region. The preliminary refinement of the foci parameters (the coordinates and origin time) has reduced the random errors in the travel-time data. The earthquake data were supplemented by the previous deep seismic sounding (DSS) data on the profiles in Crimea and offshore off the Black Sea. The dataset included more than 4000 travel times overall. In order to eliminate the crustal effect, the travel times were reduced to a surface at a depth of 35 km corresponding to the mean Moho depth in the region. The improved crustal model was used for removing the contribution of the crust from the initial data. The new tomography method, which was recently developed by one of the authors and which relies on the assumption of smoothness of the lateral velocity variations, was applied for reconstructing the velocity structure of the upper mantle beneath the Black Sea up to a depth of 95 km. The lateral velocity variation maps at different depths and the vertical velocity distributions along the meridional and sublatitudinal cross sections across the Black Sea were constructed. High velocities were revealed in the subcrustal lithosphere, and the structural difference below two subbasins—the West Black Sea (WBS) and the East Black Sea (EBS) ones—was established. It shows that the high-velocity body below the WBS is located deeper than below the EBS and is distinguished by higher velocities. Based on these results, it is concluded that the lithosphere beneath the Black Sea has a continental origin.