Seismic, Structural and Petrological Models of the Subcrustal Lithosphere in Southern Germany: A Quantitative Revaluation

—Anisotropy in the subcontinental lithosphere becomes increasingly important, because it is observed in many seismic studies especially for P _n -waves. Typical rocks of the uppermost mantle are peridotites, which predominantly exhibit a pronounced elastic anisotropy. This anisotropy is mainly caused by the anisotropic elastic properties and the lattice preferred orientation (here referred to as texture) of olivine. To evaluate the elastic anisotropy of peridotites from the subcontinental lithosphere, specimens of the Northern Hessian Depression (Germany) and the Balmuccia Ultramafic Massif (Northern Italy) have been used. They comprise four olivine texture types, which are characteristic for olivine textures observed worldwide. The bulk rock elastic properties have been calculated using olivine and orthopyroxene textures, their single-crystal elastic constants at ambient pressure/temperature conditions and their volume fraction. Clinopyroxene and spinel are assumed to be randomly distributed. The effect of four different orientations of the foliation within the uppermost mantle has been evaluated, since this orientation is usually unknown.¶Two of the olivine textures have a pronounced azimuthal dependence of compressional waves when a horizontal foliation within the uppermost mantle is presumed. These variations cause significant azimuthal variations of the P-wave reflections coefficients at the Moho. Primarily, we predict a significant azimuthal dependence of the critical points where the reflected amplitude increases from approximately 15% to 95%. Possibly, these azimuthal variations can be detected by seismic reflection measurements carried out at earth surface.¶The remaining two texture types only manifest a small directional dependence. When anisotropy of compressional waves is observed in seismic studies, these latter types can only be of subordinate importance. However, all of the peridotites investigated are able to explain the seismically observed azimuthal variations of compressional waves when a vertical foliation is proposed. This ambiguity can be substantially reduced when shear waves (S-waves) are considered. The directional distribution of S-wave velocities and of the S-wave splitting exhibits characteristic patterns for the different olivine texture types. This could be used to discriminate between different texture types and orientations of the foliation within the uppermost mantle. A fundamental requirement for a more comprehensive interpretation is the availability of detailed S-wave observations. The maximum S-wave splitting in the peridotites investigated coincides with the maximum of the faster (leading) S-wave. This may be of importance to detect S-wave splitting in future seismic studies.     

Shear velocity structure of crust and uppermost mantle in China from surface wave tomography using ambient noise and earthquake data

We present a 3D model of shear velocity of crust and upper mantle in China and surrounding regions from surface wave tomography.We combine dispersion measurements from ambient noise correlation and traditional earthquake data.The stations include the China National Seismic Network,global networks,and all the available PASSCAL stations in the region over the years.The combined data sets provide excellent data coverage of the region for surface wave measurements from 8 to 120 s,which are used to invert for 3D shear wave velocity structure of the crust and upper mantle down to about150 km.We also derive new models of the study region for crustal thickness and averaged S velocities for upper,mid,and lower crust and the uppermost mantle.The models provide a fundamental data set for understanding continental dynamics and evolution.The tomography results reveal significant features of crust and upper mantle structure,including major basins,Moho depth variation,mantle velocity contrast between eastern and western North China Craton,widespread low-velocity zone in midcrust in much of the Tibetan Plateau,and clear velocity contrasts of the mantle lithosphere between north and southern Tibet with significant E–W variations.The low velocity structure in the upper mantle under north and eastern TP correlates with surface geological boundaries.A patch of high velocity anomaly is found under the eastern part of the TP,which may indicate intact mantle lithosphere.Mantle lithosphere shows striking systematic change from the western to eastern North China Craton.The Tanlu Fault appears to be a major lithosphere boundary.     

Correlating P-wave Velocity with the Physico-Mechanical Properties of Different Rocks

In mining and civil engineering projects, physico-mechanical properties of the rock affect both the project design and the construction operation. Determination of various physico-mechanical properties of rocks is expensive and time consuming, and sometimes it is very difficult to get cores to perform direct tests to evaluate the rock mass. The purpose of this work is to investigate the relationships between the different physico-mechanical properties of the various rock types with the P-wave velocity. Measurement of P-wave velocity is relatively cheap, non-destructive and easy to carry out. In this study, representative rock mass samples of igneous, sedimentary, and metamorphic rocks were collected from the different locations of India to obtain an empirical relation between P-wave velocity and uniaxial compressive strength, tensile strength, punch shear, density, slake durability index, Young’s modulus, Poisson’s ratio, impact strength index and Schmidt hammer rebound number. A very strong correlation was found between the P-wave velocity and different physico-mechanical properties of various rock types with very high coefficients of determination. To check the sensitivity of the empirical equations, Students t test was also performed, which confirmed the validity of the proposed correlations.     

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

Ultrasonic elastic characteristics of five kinds of metamorphic deformed coals under room temperature and pressure conditions

The calibration of the elastic characteristics of deformed coals is essential for seismic inversion of such units, because the prediction of coal deformation is essential for both mining safety and methane production. Therefore, many samples of broken and mylonitic deformed coal were tested with ultrasonic waves in the laboratory. These samples came from four mining areas: the Huainan, Pingdingshan, Hebi and Jiaozuo coal mines, which present five different metamorphic ranks shown as cylinders striking across circular limits of steel. Under normal pressures and temperatures, ultrasonic P- and S-wave tests show that the velocities, quality factors, and elastic moduli of the deformed coals were greatly reduced compared with undeformed coals. Also, some correlation was found between the P- and S-wave velocities in the deformed coals. However, there is no evidence of linear correlations between velocity and density, velocity and quality factor, or the quality factors of P- and S-waves. Compared with the elastic characteristics of undeformed coals, such as P- and S-wave velocity ratios or Poisson’s ratio, those of deformed coals generally decrease and the P-wave quality factors are less than those of S-waves. Moreover, the analysis of the relationship between pore structure and elastic modulus shows a better correlation between the P- and S-wave velocities and effective porosity, pore volume and specific surface area. Also, there are similar relationships between the pore structure and the Young’s and shear moduli. However, there are no such correlations with other moduli. Correlations between these elastic moduli, pore structure, coal rank and density were not found for the various samples of deformed coals, which is consistent with only structural destruction occurring in the deformed coals with other physical properties remaining unchanged. The experimental results show that it is possible to predict the deformation of coals with multi-component seismic elastic inversion.     

Crust and mantle of the Tien Shan from data of the receiver function tomography

A 3-D velocity model of the Tien Shan crust and upper mantle is constructed through the inversion of the receiver functions of P and S waves together with teleseismic traveltime anomalies at nearly 40 local seismic stations. It is found that in the vast central region, where no strong earthquakes have been known over the past century, the S wave velocity at depths of 10–35 km is lower than in adjacent regions by up to 10%. These data are evidence for mechanical weakness of the crust preventing the accumulation of elastic energy. Apparently, the lower velocity and the weakness of the crust are due to the presence of water. The weakness of the crust is one of the possible reasons for the strain localization responsible for the formation of the present Tien Shan but can also be due in part to the young orogenesis. The crustal thickness is largest (about 60 km) in the Tarim-Tien Shan junction zone. The crust-mantle boundary in this region descends by a jump as a result of an increase in the lower crust thickness. This is probably due to the underthrusting of the Tien Shan by the Tarim lithosphere. This causes the mechanically weak lower crust of the Tarim to delaminate and accumulate in nearly the same way as an accretionary prism during the subduction of oceanic lithosphere. In the upper mantle, the analysis has revealed a low velocity anomaly, apparently related to basaltic outflows of the Upper Cretaceous-Early Paleogene. The Cenozoic Bachu uplift in the northern Tarim depression is also associated with the low velocity anomaly. The Naryn depression is characterized by a high velocity in the upper mantle and can be interpreted as a fragment of an ancient platform.     

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.     

柴达木盆地深部物理参数的反演

通过研究地球对表面载荷的响应来反演确定地球深部的物理参数.采用分为4层的地球模型,各层的介质被看作Maxwell粘弹性体.将柴达木盆地沉降发育过程中积累的沉积物作为随时间变化的载荷加于地表,通过拟合柴达木盆地沉降中心附近地区的沉降过程来研究该地区地壳的厚度及剪切模量、岩石层地幔的粘性及软流层的粘性等物理参数.计算结果表明,该地区地壳的剪切模量不会超过全球平均值的60％,岩石层地幔的粘性系数不大于1.5×10²³Pa·s.通过研究对该地区的深部物理参数得出了新的认识.     

Joint inversion of P- and S-receiver functions and dispersion curves of Rayleigh waves: The results for the Central Anatolian Plateau

The P- and S-wave receiver functions and dispersion curves of the fundamental Rayleigh wave are used to study the lithosphere within the Central Anatolian Plateau. The results for eight broadband seismic stations are presented. It is established that within the plateau, the crust with a thickness of about 35 km is underlain by the mantle lid with its bottom at a depth of about 60 km. The velocities of longitudinal (Vp) and shear (Vs) waves in this layer are at most 7.6 and 4.5 km/s, respectively, and the Vp/Vs ratio is close to 1.7 (i.e., by 6% lower than in the standard IASP91 and PREM models). Such a low velocity ratio is characteristic of rocks having high orthopyroxene content. Beneath the high-velocity mantle lid, the S-wave velocity decreases to 4.0–4.2 km/s and the Vp/Vs ratio is close to its standard value (1.8). At most stations, the P-wave receiver functions do not contain seismic phase P410s, which is formed at the global seismic boundary at a depth of 410 km. The seismic boundary at a depth of 410 km is related to the olivine-spinel phase transformation, and its absence can indicate the anomalously low olivine content and high basalt content. This anomaly is probably associated with the subduction of a large amount of oceanic crust during the closure of the Tethys. The results of the study overall indicate the high informativity of the used method.     

Scales of Heterogeneities in the Continental Crust and Upper Mantle

—A seismological characterization of crust and upper mantle can refer to large-scale averages of seismic velocities or to fluctuations of elastic parameters. Large is understood here relative to the wavelength used to probe the earth.¶In this paper we try to characterize crust and upper mantle by the fluctuations in media properties rather than by their average velocities. As such it becomes evident that different scales of heterogeneities prevail in different layers of crust and mantle. Although we cannot provide final models and an explanation of why these different scales exist, we believe that scales of inhomogeneities carry significant information regarding the tectonic processes that have affected the lower crust, the lithospheric and the sublithospheric upper mantle.¶We focus on four different types of small-scale inhomogeneities (1) the characteristics of the lower crust, (2) velocity fluctuations in the uppermost mantle, (3) scattering in the lowermost lithosphere and on (4) heterogeneities in the mantle transition zone.     

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.     

Velocity properties of the lithosphere in the ocean-continent transition zone in the Kamchatka region from seismic tomography data

Arrival times of P and S waves from local earthquakes in the Kamchatka area of the Kurile-Kamchatka Island Arc are used for calculating a spatial model of the elastic wave velocity distribution to a depth of 200 km. The lithosphere is shown to be strongly stratified in its velocity properties and laterally heterogeneous within the mantle wedge and seismic focal zone. A lower velocity layer (an asthenospheric wedge) is identified at depths of 70–130 km beneath the Eastern Kamchatka volcanic belt. The morphology of the Moho interface and the velocity properties of the crust are studied. The main tectonic structures of the region are shown to be closely interrelated with deep velocity heterogeneities. Regular patterns in the statistics of the earthquakes are analyzed in relation to variations in the elastic wave velocities in the focal layer. A mechanism of lithospheric block displacements along weakened zones of the lower crust and upper mantle is proposed.     

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.     

3-D P-wave Velocity Structure in Western Greece Determined from Tomography Using Earthquake Data Recorded at the University of Patras Seismic Network (PATNET)

—The 3-D P-wave velocity structure of the upper crust in the region of western Greece is investigated by inversion of about 1500 residuals of P-wave arrival times from local earthquake data recorded in the year 1996 by the newly established University of Patras Seismic Network (PATNET). The resulting velocity structure shows strong horizontal variations due to the complicated structure and the variation of crustal thickness. Relatively low-velocity contours are observed in the area defined by Cephallonia—Zakynthos Islands and northwestern Peloponnesos. This is in addition to some well localized peaks of relatively higher values of P-wave velocity may be related to the zone of Triassic evaporites in the region and correspond to diapirism that breaks through to the uppermost layer. Finally, a low P-velocity ‘deeping’ zone extending from Zakynthos to the Gulf of Patras is correlated with Bouguer anomaly map and onshore and offshore borehole drillings which indicate that thick sediments overly the evaporites which exist there at depth greater than 2 km.     

Siberian traps: Hypotheses and seismology data

Siberian traps are the result of huge basalt eruptions which took place about 250 Ma ago over a vast territory of Siberia. The genesis of Siberian traps is attributed to a mantle plume with a center in the region of Iceland or beneath the central Urals in terms of their present coordinates. The eruption mechanism is associated with delamination—replacement of the mantle lithosphere by the deep magma material. The receiver function analysis of the records from the Norilsk seismic station (NRIL) allows comparing these hypotheses with the factual data on the depth structure of the region of Siberian traps. The S-wave velocity section place the seismic lithosphere/asthenosphere boundary (LAB) at a depth of 155–190 km, commensurate with the data for the other cratons. The mantle lithosphere has a high S-wave velocity characteristic of cratons (4.6–4.8 km/s instead of the typical value 4.5 km/s). The seismic boundary, which is located at a depth around 410 km beneath the continents is depressed by ~10 km in the region of the NRIL station. The phase diagram of olivine/wadsleyite transformation accounts for this depression by a 50–100°С increase in temperature. At the depths of 350–400 km, the S-wave velocity drops due to partial melting. A new reduction in the S-wave velocities is observed at a depth of 460 km. The similar anomalies (deepening of the 410-km seismic boundary and low shear wave velocity at depths of 350–400 and 460–500 km, respectively) were previously revealed in the other regions of the Meso-Cenozoic volcanism. In the case of a differently directed drift of the Siberian lithosphere and underlying mantle at depths down to 500 km, these anomalies are barely accountable. In particular, if the mantle at a depth ranging from 200 to 500 km is fixed, the anomalies should be observed at the original locations where they emerged 250 Ma ago, i.e. thousands of km from the Siberian traps. Our seismic data suggest that despite the low viscosity of the asthenosphere, the mantle drift at depths ranging from 200 to 500 km is correlated with the drift of the Siberian lithospheric plate. Furthermore, the position of the mantle plume beneath the Urals is easier to reconcile with the seismic data than its position beneath Iceland because of the Siberian traps being less remote from the Urals.     

Mechanical behavior of magma reservoir envelopes: Elasticity of the olivine tholeiite solidus

A study of the elasticity of the solidus of igneous rocks provides a means of evaluating several high temperature mechanical characteristics of the zone that forms the transition region between solid rock and rock with a small fraction of partial melt — that is, the magma reservoirenvelope. This type of study has been performed for Hawaiian olivine tholeiite, by employing the theoretical treatment ofWalsh (1969). The effective aggregate bulk modulus (K*), aggregate rigidity (μ*), aggregate Young’s modulus (E*) and the aggregate Poisson’s ratio (n*) have been used to characterize the range in elastic behavior expected for this partially molten envelope. The details of mineralogy, mineral crystal chemistry and microstructure for a tholeiitic basalt from a deep-ponded flow unit (the Boiling Pots), as well as the geometry of melt pocket aspect ratios have been used to provide a basis for computing the elasticity for a shallow reservoir solidus. For a broad range in modal mineralogy, the internal 0–1% melt represents a transition region from an aggregate elasticity controlled by mineral composition, to an elastic behavior dominated by the melt phase (evaluated at a melt pocket aspect ratio ofε=0.001). For vesicular rocks in a low pressure environment, the internal 0–1% melt represents a transition from porosity-controlled elasticity to melt-controlled elasticity, evaluated at a prevalent aspect ratio ofε=0.001, and over the porosity range Φ=0.0?0.20 (20%). Thin melt films found to separate plagioclase microlites (0.0001=ε) are capable of effecting a 94.7% reduction in Young’s modulus, E*, over the 0–1% melt range, and suggest that the reservoir envelope may be critically weakened by rather small volume fractions of liquid, if dispersed in such narrow packets.     

Accumulation of elastic strains in the upper crust on the locked transform faults and the tectonomagnetic effect

Plate tectonics only allows small deformations in the lithospheric plates. The laboratory experiments with the rock specimens show that the creep is transient when the creep strain is at most 1%. Hence, if we assume that the creep strain in the lithospheric plates is below this threshold, the creep is transient. The present paper addresses the role of the elastic, brittle (pseudo-plastic), and creep rheology of the lithosphere during the accumulation of elastic shear strains on the locked faults in the Earth’s crust, i.e., during the process of preparation of the earthquakes. The effective viscosity characterizing the transient creep is lower than that under the steady-state creep and it depends on the characteristic time of a given process. The characteristic duration of the stress and strain accumulation process in the vicinity of the locked faults is a few dozen years. On these time intervals, the thin upper crustal layer behaves as brittle; the underlying layer behaves as elastic (it is just this layer which accommodates stress accumulation leading to the earthquake), whereas the transient creep is predominant in the lower crust and mantle lithosphere. Transient creep entails nonlinear time dependence of the strains arising in the vicinity of the locked fault in the elastic crust. The perturbations in the magnetic field induced by these strains can be treated as the magnetic precursor of the earthquake.     

Distributions of density and elastic parameters in the Earth

The velocity curve previously obtained for P waves in the Earth is used to determine the distributions of the density and elastic parameters. The density distribution in the new model differs from that in the standard PREM model only in the inner core. The distributions of the bulk and shear moduli can differ, depending on physical processes in the Earth. In particular, the bulk modulus can have a negative jump at the outer-inner core boundary, whereas the shear modulus can differ from zero in the lower part of the outer core.     

The evolution of the stress state in Southern California based on the geomechanical model and current seismicity

A three-dimensional geomechanical model of Southern California, which includes the mountain topography, fault tectonics, and main structural boundaries (the top of the lower crust and the Moho), is developed. The main stress state of the model is determined by the own weight of the rocks and by the horizontal tectonic motions identified from the GPS observations. The model enables tracking the changes which occur in the stress-strain state of the crust due to the evolution of the seismic process. As the input data, the model uses the current seismicity and treats each earthquake as a new defect in the Earth’s crust which brings about the redistribution of strains, elastic energy density, and yield stress of the crust. Monitoring the variations in the stress state of the crust and lithosphere arising in response to the seismic process shows that the model is suitable for forecasting the enhancement in seismic activity of the region and delineating the earthquake-prone areas with a reasonable probability on a given time interval.     

Ultrasonic Velocities, Acoustic Emission Characteristics and Crack Damage of Basalt and Granite

Acoustic emissions (AE), compressional (P), shear (S) wave velocities, and volumetric strain of Etna basalt and Aue granite were measured simultaneously during triaxial compression tests. Deformation-induced AE activity and velocity changes were monitored using twelve P-wave sensors and eight orthogonally polarized S-wave piezoelectric sensors; volumetric strain was measured using two pairs of orthogonal strain gages glued directly to the rock surface. P-wave velocity in basalt is about 3 km/s at atmospheric pressure, but increases by > 50% when the hydrostatic pressure is increased to 120 MPa. In granite samples initial P-wave velocity is 5 km/s and increases with pressure by < 20%. The pressure-induced changes of elastic wave speed indicate dominantly compliant low-aspect ratio pores in both materials, in addition Etna basalt also contains high-aspect ratio voids. In triaxial loading, stress-induced anisotropy of P-wave velocities was significantly higher for basalt than for granite, with vertical velocity components being faster than horizontal velocities. However, with increasing axial load, horizontal velocities show a small increase for basalt but a significant decrease for granite. Using first motion polarity we determined AE source types generated during triaxial loading of the samples. With increasing differential stress AE activity in granite and basalt increased with a significant contribution of tensile events. Close to failure the relative contribution of tensile events and horizontal wave velocities decreased significantly. A concomitant increase of double-couple events indicating shear, suggests shear cracks linking previously formed tensile cracks.