Seismic Wave Velocities in Deep Sediments in Poland: Borehole and Refraction Data Compilation

Sedimentary cover has significant influence on seismic wave travel times and knowing its structure is of great importance for studying deeper structures of the Earth. Seismic tomography is one of the methods that require good knowledge of seismic velocities in sediments and unfortunately by itself cannot provide detailed information about distribution of seismic velocities in sedimentary cover. This paper presents results of P-wave velocity analysis in the old Paleozoic sediments in area of Polish Lowland, Folded Area, and all sediments in complicated area of the Carpathian Mountains in Poland. Due to location on conjunction of three major tectonic units — the Precambrian East European Craton, the Paleozoic Platform of Central and Western Europe, and the Alpine orogen represented by the Carpathian Mountains the maximum depth of these sediments reaches up to 25 000 m in the Carpathian Mountains. Seismic velocities based on 492 deep boreholes with vertical seismic profiling and a total of 741 vertical seismic profiles taken from 29 seismic refraction profiles are analyzed separately for 14 geologically different units. For each unit, velocity versus depth relations are approximated by second or third order polynomials.     

Improved streamlines and time-of-flight for streamline simulation on irregular grids

Streamline methods have shown to be effective for reservoir simulation. For a regular grid, it is common to use the semi-analytical Pollock’s method to obtain streamlines and time-of-flight coordinates (TOF). The usual way of handling irregular grids is by trilinear transformation of each grid cell to a unit cube together with a linear flux interpolation scaled by the Jacobian. The flux interpolation allows for fast integration of streamlines, but is inaccurate even for uniform flow. To improve the tracing accuracy, we introduce a new interpolation method, which we call corner-velocity interpolation. Instead of interpolating the velocity field based on discrete fluxes at cell edges, the new method interpolates directly from reconstructed point velocities given at the corner points in the grid. This allows for reproduction of uniform flow, and eliminates the influence of cell geometries on the velocity field. Using several numerical examples, we demonstrate that the new method is more accurate than the standard tracing methods.     

Structure of the upper mantle in Asia from phase and group velocities of Rayleigh waves

Dispersion curves of phase velocities of Rayleigh waves are determined by the method of frequency-time analysis in a range of periods of 10–200 s from data of 43 interstation traces in Central Asia. Because the joint use of phase and group velocities significantly decreases the uncertainty in the determination of S wave velocity structures, the same traces were used for calculating group velocities from tomographic reconstructions obtained in [Yanovskaya and Kozhevnikov, 2003, 2006] and determining average velocity structures along these traces. The velocity structures were calculated by the Monte Carlo and linear inversion methods, which gave consistent results. Using velocity values obtained at fixed depths by the 2-D tomography method, lateral variations in velocities at these depths were estimated, which allowed us to construct smoothed vertical velocity structures at some points in the region. The resulting structures were used as initial approximations for constructing local velocity structures solely from previously obtained local dispersion curves of group velocities in the area (32°–56°N, 80°–120°E). Based on these structures, we mapped the lateral distribution of velocity variations at upper mantle depths of 75–400 km and along three vertical profiles. The inferred velocity variations are in good agreement with data on the tectonics of the region.     

Layered Velocity Models of the Western Bohemia Region

A new robust and effective optimization algorithm – isometric algorithm – was used for the inversion of layered velocity models, with constant gradient in each layer, to find suitable 1-D models for the location of microearthquakes in the individual four subregions of the West Bohemian earthquake swarm region. Models which are considered as optimal yield the minimum sum of the absolute values of the travel-time residua in locating the whole group of earthquakes in the given subregion. The results obtained from the inversion of P and S waves and from P waves only are shown. For comparison, optimum homogeneous models derived by the grid search method, again using both P and S waves and P waves only, are given. The computations indicate that the models for the individual subregions differ from each other. For layered models the differences are more pronounced, as expected, in the upper parts, down to depths of about 5 km. In comparison with the subregions Nový Kostel and Plesná, the P and S wave velocities for subregion Lazy are relatively higher and the P and S velocities for subregion Klingenthal relatively lower. In the lower parts the differences are smaller and the velocities have practically identical gradients. The highest velocities were obtained for subregion Lazy and the lowest velocities for subregion Klingenthal, as well for the homogeneous models. The model that represents the whole swarm region was determined in a similar way. This model is compared with the previously published velocity-depth distribution, obtained from DSS profile VI/70 in the vicinity of the area under study.     

基于偏微分方程的定向插值在偏移成像中的应用

对于射线类偏移成像来说,求解射线追踪系统中所涉及的属性值不在网格节点上的插值计算问题是一个非常重要的环节,它影响到求解走时、路径和振幅信息的计算效率和精度,进而影响到整个偏移成像的质量和效率.本研究根据速度模型的空间梯度特点,考虑被插值点处速度的梯度在横向和纵向的分布特征,构建基于速度梯度空间变化的偏微分方程算法,将近几年发展起来的基于偏微分方程的定向插值算法引入到射线类偏移成像当中,实现射线追踪当中涉及的属性值不在网格节点上的插值计算.由于偏微分方程法本身固有的特性(局部特征不变性、解的唯一性和线性叠加性),因此,该算法可以实现不破坏原始速度模型空间梯度结构的非网格节点属性的插值计算.通过在常用的速度模型上的插值计算对比、不同速度模型上射线路径对比分析以及复杂介质模型上最后的偏移成像结果分析可以得出,应用基于速度梯度构建的偏微分方程插值算法在进行插值计算的过程当中可以实现不破坏原始速度模型空间速度梯度结构的属性计算,同时应用该算法可以最终提高射线类偏移成像的质量.     

Velocity structure of the lithosphere on the 2003 Mongolian-Baikal transect from <Emphasis Type="Italic">SV</Emphasis> waves

The S wave velocity distribution in the Earth’s crust and the first two hundred kilometers of the upper mantle is inferred from data of a seismological linear network including 18 broadband stations installed in the framework of the international teleseismic experiment carried out in 2003 in the south of Siberia and in Mongolia. Models were constructed by using P-to-S received function inversion beneath each station. Vertical cross sections of S wave velocities from the surface to depths of 65 and 270 km covering the entire 1000-km profile are constructed by the linear spline interpolation of individual velocity models. The vertical sections are also represented as anomalies relative to the standard velocity model. The most intense low velocity anomalies (from ?3 to ?6%) in the crust and upper mantle are identified beneath the Sayan, Khamar-Daban, and Khangai highlands and the Djida fold zone and agree both with other geophysical data and with the distribution of Late Cenozoic volcanic fields. The results of this work suggest that the activation of Mongolian-Siberian highlands is largely connected with uplift of the asthenosphere to the base of the crust.     

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 wave velocities in the sedimentary cover of Poland: Borehole data compilation

A knowledge of seismic wave velocities in the sedimentary cover is of great importance for interpreting reflection and refraction seismic data, deep seismic soundings and regional and global seismic tomography. This is particularly true for regions characterized by significant thicknesses and a complex sedimentary cover structure. This paper presents the results of an analysis of seismic P-wave velocities in the sedimentary cover of Poland, a complex area of juxtaposition of major tectonic units: the Precambrian East European Craton, the Palaeozoic Platform of Central and Western Europe, and the Alpine orogen represented by the Carpathian Mountains. Based on vertical seismic profiling data from 1188 boreholes, the dependence of velocity versus depth was determined for regional geological units and for successions from the Tertiary and Quaternary to the Cambrian. The data have been approximated by polynomials, and velocity-depth formulas are given down to 6000 m depth. The velocities in the sedimentary cover have been compared with those from other areas in Europe.     

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.     

Seismic Velocity and Q Structure of the Middle Eastern Crust and Upper Mantle from Surface-wave Dispersion and Attenuation

—Observed velocities and attenuation of fundamental-mode Rayleigh waves in the period range 7–82 sec were inverted for shear-wave velocity and shear-wave Q structure in the Middle East using a two-station method. Additional information on Q structure variation within each region was obtained by studying amplitude spectra of fundamental-mode and higher-mode Rayleigh waves. We obtained models for the Turkish and Iranian Plateaus (Region 1), areas surrounding and including the Black and Caspian Seas (Region 2), and the Arabian Peninsula (Region 3). The effect of continent-ocean boundaries and mixed paths in Region 2 may lead to unrealistic features in the models obtained there. At lower crustal and upper-mantle depths, shear velocities are similar in all three regions. Shear velocities vary significantly in the uppermost 10 km of the crust, being 3.21, 2.85, and 3.39 km/s for Regions 1, 2, and 3, respectively. Q models obtained from an inversion of interstation attenuation data show that crustal shear-wave Q is highest in Region 3 and lowest in Region 1. Q’s for the upper 10 km of the crust are 63, 71, and 201 for Regions 1, 2, and 3, respectively. Crustal Q’s at 30 km depth for the three regions are about 51, 71, and 134. The lower crustal Q values contrast sharply with results from stable continental regions where shear-wave Q may reach one thousand or more. These low values may indicate that fluids reside in faults, cracks, and permeable rock at lower crustal, as well as upper crustal depths due to convergence and intense deformation at all depths in the Middle Eastern crust.     

Shear Wave Velocity Estimates through Combined Use of Passive Techniques in a Tectonically Active Area

We made an attempt to assess the shear wave velocity values V _S and, to a lesser extent, the V _P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V _S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.     

Gradients and phase velocities of ULF geomagnetic disturbances used to determine the source of an impending strong earthquake

Results of studying the behavior of the vectors of gradients and phase velocities of ULF geomagnetic disturbances (F < 1 Hz) in the Japan seismic zone are presented. The gradient and phase velocity vectors along the Earth’s surface were determined using data of the group of three high-sensitivity three-component magnetovariation stations (MVC-3DS) located at triangle vertices at a small (～5 km) distance from one another (magnetic gradiometer). Two such groups of stations were installed in 1999 southwest and southeast of Tokyo at a distance of ～150 km from each other. It has been indicated that, several months before strong earthquakes (M > 5), the values of gradients and phase velocities start anomalously changing, and directions toward sources of impending strong earthquakes appear in the distribution of gradient vector directions. Directions from sources of impending earthquakes appear in the distribution of phase velocity vector directions. It is proposed to use gradients and phase velocities of ULF and ELF geomagnetic disturbances as one of the factors in a short-term prediction of strong earthquakes.     

The structure of the crust in TESZ area by kriging interpolation

A precise 3D model of the crust is necessary to start any tectonic or geodynamic interpretation. It is also essential for seismic interpretations of structures lying below as well as for correct analysis of shallow structures using reflection seismics. During the last decades, a number of wide-angle refraction experiments were performed on the territory of central and eastern Europe (POLONAISE’97, CELEBRATION 2000, SUDETES 2003), resulting in many high quality 2D models. It is an interesting and complicated transition zone between Precambrian and Palaeozoic Platforms. This paper presents 3D model of the velocity distribution in the crust and upper mantle interpolated from 2D models of the structure along 33 profiles. The obtained model extends to a depth of 50 km and accurately describes the main features of the crustal structures of Poland and surrounding areas. Different interpolation techniques (Kriging, linear) are compared to assure maximum precision. The final model with estimated uncertainty is an interesting reference of the area for other studies.     

Spatial relation between source regions of strong earthquakes (M ≥ 7) on Southern Kamchatka and the distribution of velocities and elastic moduli in the Benioff zone

The presence of a phenomenological relationship between high velocity regions in the Benioff zone and sources of relatively strong earthquakes (M ≥ 6) was established for the first time from the comparison of such earthquakes with the velocity structure of central Kamchatka in the early 1970s. It was found that, in the region with P wave velocities of 8.1–8.5 km/s, the number of M ≥ 6 earthquakes over 1926–1965 was 2.5 times greater than their number in the region with velocities of 7.5–8.0 km/s. Later (in 1979), within the southern Kurile area, Sakhalin seismologists established that regions with V _P = 7.3–7.7 km/s are associated with source zones of M = 7.0–7.6 earthquakes and regions with V _P = 8.1–8.4 km/s are associated with M = 7.9–8.4 earthquakes. In light of these facts, we compared the positions of M = 7.0–7.4 earthquake sources in the Benioff zone of southern Kamchatka over the period 1907–1993 with the distribution of regions of high P velocities (8.0–8.5 to 8.5–9.0 km/s) derived from the interpretation of arrival time residuals at the Shipunskii station from numerous weak earthquakes in this zone (more than 2200 events of M = 2.3–4.9) over the period 1983–1995. This comparison is possible only in the case of long-term stability of the velocity field within the Benioff zone. This stability is confirmed by the relationship between velocity parameters and tectonics in the southern part of the Kurile arc, where island blocks are confined to high velocity regions in the Benioff zone and the straits between islands are confined to low velocity regions. The sources of southern Kamchatka earthquakes with M = 7.0–7.4, which are not the strongest events, are located predominantly within high velocity regions and at their boundaries with low velocity regions; i.e., the tendency previously established for the strongest earthquakes of the southern Kuriles and central Kamchatka is confirmed. However, to demonstrate more definitely their association with regions of high P wave velocities, a larger statistics of such earthquakes is required. On the basis of a direct correlation between P wave velocities and densities, the distributions of density, bulk modulus K, and shear modulus μ in the upper mantle of the Benioff zone of southern Kamchatka are obtained for the first time. Estimated densities vary from 3.6–3.9 g/cm³ in regions of high V _P values to 3.0–3.2 g/cm³ for regions of low V _P values. The bulk modulus K in the same velocity regions varies from (1.4–1.8) × 10¹² to (0.8–1.1) × 10¹² dyn/cm², respectively, and the shear modulus μ varies from (0.8–1.0) × 10¹² to (0.5–0.7) × 10¹² dyn/cm², respectively. Examination of the spatial correlation of the source areas of southern Kamchatka M = 7.0–7.4 earthquakes with the distribution of elastic moduli in the Benioff zone failed to reveal any relationship between their magnitudes and the moduli because of the insufficient statistics of the earthquakes used.     

Performance of geostatistical interpolation methods for modeling sampled data with non-stationary mean

The measured ozone pollution peak in the atmosphere of Mexico City region was considered in order to study the effect of a non-stationary mean of the sampled data in geostatistics interpolation methods. With this objective the local mean value of the sampled data was estimated through a linear regression analysis of their values on the monitoring station’s coordinates. The residuals obtained by removing the data trend are considered as a set of stationary random variables. Several interpolation methods used in geostatistics, such as inverse distance weighted, kriging, and artificial neural networks techniques were considered. In an effort to optimize and evaluate its performance, we fit interpolated values to sampled data, obtaining optimal values for the parameters defining the used model, that means, the values of the parameters that give the lowest mean RMSE between the interpolated value and measured data at 20 stations at 1500 hours for a set of 21 days of December 2001, which was chosen as the training set. The training set is conformed by all the days in December 2001 excepting the days (3,6,9,12,...,27,30) which were considered as the testing set. Once the optimal model is obtained, it is used to interpolate the values at the stations at 1500 hours for the testing days. The RMSE between interpolated and measured values at monitoring stations was also evaluated for these testing values and is shown as a percentage in Table 2. These values and the defined generalization parameter G, can be used to evaluate the performance and the ability of the models to predict and reproduce the peak of ozone concentrations. Scatter plots for testing data are presented for each interpolation method. An interpretation of the ozone pollution levels obtained at 1500 hours at December 21 was given using the wind field that prevailed in the region 1 h before the same day.     

Upper mantle structure beneath central Western Europe from data on phase velocities of surface waves

The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.     

Shear wave velocity models retrieved using Rg wave dispersion data in shallow crust in some regions of southern Ontario, Canada

Many crucial tasks in seismology, such as locating seismic events and estimating focal mechanisms, need crustal velocity models. The velocity models of shallow structures are particularly important in the simulation of ground motions. In southern Ontario, Canada, many small shallow earthquakes occur, generating high-frequency Rayleigh (Rg) waves that are sensitive to shallow structures. In this research, the dispersion of Rg waves was used to obtain shear-wave velocities in the top few kilometers of the crust in the Georgian Bay, Sudbury, and Thunder Bay areas of southern Ontario. Several shallow velocity models were obtained based on the dispersion of recorded Rg waves. The Rg waves generated by an m _N 3.0 natural earthquake on the northern shore of Georgian Bay were used to obtain velocity models for the area of an earthquake swarm in 2007. The Rg waves generated by a mining induced event in the Sudbury area in 2005 were used to retrieve velocity models between Georgian Bay and the Ottawa River. The Rg waves generated by the largest event in a natural earthquake swarm near Thunder Bay in 2008 were used to obtain a velocity model in that swarm area. The basic feature of all the investigated models is that there is a top low-velocity layer with a thickness of about 0.5 km. The seismic velocities changed mainly within the top 2 km, where small earthquakes often occur.     

A Condition for Super-shear Rupture Propagation in a Heterogeneous Stress Field

— We have used numerical simulations with the boundary integral equation method to investigate a mechanism to excite super-shear rupture velocities in a homogeneous stress field including an asperity of increased initial stress. When the rupture, with the slip-weakening distance selected to generate sub-Rayleigh speed, encounters the asperity it either accelerates to super-shear velocities or maintains the sub-Rayleigh speed, dependent on the size and amplitude of the asperity. Three classes of rupture propagation are identified: the velocity (a) for the most narrow asperities increases slowly towards the Rayleigh wave speed, (b) for intermediate width of the asperities jumps to super-shear values for a short distance but then decreases to sub-Rayleigh wave speeds, and (c) for the widest asperities jumps to super-shear values and pertains to values between the S- and P-wave velocities. The transitions between the three classes of rupture propagation are characterized by very narrow (critical) ranges of rupture resistance. If the size of the initial asperity is smaller than critical, it becomes difficult for rupture to propagate with super-shear velocities even if the initial stress level is high. Our results suggest that stress variation along the rupture path helps homogenize the rupture velocity and propagate with sub-Rayleigh wave speeds.     

A Surface Wave Dispersion Study of the Middle East and North Africa for Monitoring the Comprehensive Nuclear-Test-Ban Treaty

—?We present results from a large-scale study of surface-wave group velocity dispersion across the Middle East, North Africa, southern Eurasia and the Mediterranean. Our database for the region is populated with seismic data from regional events recorded at permanent and portable broadband, three-component digital stations. We have measured the group velocity using a multiple narrow-band filter on deconvolved displacement data. Overall, we have examined more than 13,500 seismograms and made good quality dispersion measurements for 6817 Rayleigh- and 3806 Love-wave paths. We use a conjugate gradient method to perform a group-velocity tomography. Our current results include both Love- and Rayleigh-wave inversions across the region for periods from 10 to 60 seconds. Our findings indicate that short-period structure is sensitive to slow velocities associated with large sedimentary features such as the Mediterranean Sea and Persian Gulf. We find our long-period Rayleigh-wave inversion is sensitive to crustal thickness, such as fast velocities under the oceans and slow along the relatively thick Zagros Mts. and Turkish-Iranian Plateau. We also find slow upper mantle velocities along known rift systems. Accurate group velocity maps can be used to construct phase-matched filters along any given path. The filters can improve weak surface wave signals by compressing the dispersed signal. The signals can then be used to calculate regionally determined M _S measurements, which we hope can be used to extend the threshold of m _b :M _S discriminants down to lower magnitude levels. Other applications include using the group velocities in the creation of a suitable background model for forming station calibration maps, and using the group velocities to model the velocity structure of the crust and upper mantle.     

Correlated Crust and Mantle Strain Field in China

INTRODUCTION Moststudiesonactiveblockshavebeenfocusedonidentificationofblockboundariesandtheiractivity;inotherwords,mostoftheworkwaslimitedtothehorizontalmovementoftheblocks.Inreality,theblocksarenotonlysurroundedbyactivefaultsofhorizontalmotion,butalsoco…