Fabric symmetry of low anisotropic rocks inferred from ultrasonic sounding: Implications for the geomechanical models

The geomechanical models were established based on the absence or presence of certain rock fabric elements — texture (crystallographic preferred orientation), microstructure (shape preferred orientation) and microcracks (flat voids). The proposed models include both (i) the ideal material showing random texture and structure but no microcracks, i.e. the material which is hardly to be found in nature, and (ii) the materials possessing various combinations of fabric elements that show different spatial arrangements. The mutual relationship between those parameters and seismic and geomechanical properties are discussed.Selected models were experimentally verified during laboratory experiments. These consist of measurement of P-wave velocities in 132 independent directions under several confining pressures in the range 0.1–400 MPa. From measured data 3D P-wave patterns can be constructed and the influence of microcracks and of texture and structure on the rock seismic anisotropy can be determined. The seismic anisotropy established at different levels of confining pressure can be used for the interpretation of rock fabric symmetry of rocks showing low anisotropy in macroscale and for the selection of directions in which the geomechanical test can be performed. The measured P-wave velocities were then mathematically processed by using a fitting function which reflects contribution of P-wave velocity in the mineral skeleton of an ideal sample without microcracks extrapolated to the atmospheric pressure level from high confining pressure interval (ca. 200–400 MPa) (v₀), linear compressibility of the samples (k_v), and confining pressure during which most of the cracks are closed (P₀). These parameters improve the understanding of the response of various rock fabric elements on increasing confinement and corresponding changes in elasticity.The observed seismic and geomechanical anisotropies reflect intensity of the fabric of rock-forming minerals and microcracks. The magnitude of seismic anisotropy measured at atmospheric pressure corresponds to the anisotropy of static elastic modulus and is governed by the spatial arrangement of microcracks. The magnitude of strength anisotropy (uniaxial compressive strength) correlates more likely to the seismic anisotropy determined at high confining pressure and is connected to the preferred orientations (either CPO or SPO or both) of rock-forming minerals.     

韧性剪切带中片麻岩和超高压榴辉岩变形特征及其与地震波速各向异性的关系：来自中国大陆科学钻探（CCSD）680～1200米岩心的证据

中国大陆科学钻探(CCSD)680-1200米区段发育了多个韧性剪切带,带中主要岩石类型包括片麻岩和超高压榴辉岩。片麻岩中的变形石英、面理化榴辉岩中的拉长石榴石和绿辉石的应变轴比都表现为X>Y>Z,Flinn系数分别为0.11-0.27、0.22-0.23和0.23-0.24。随着糜棱岩化作用的增强,变形石英的C轴组构由Z轴极密逐渐向Y轴极密和叶理面上的大圆环带转变。在常温常压下测试了样品的波速,计算出片麻岩Vp和Vs的各向异性分别为30.17％-60.97％和11.52％-35.79％,榴辉岩Vp和Vs的各向异性分别为0.17％-11.19％和2.41％-6.70％。影响各向异性的主要因素有岩石的结构构造、矿物的晶格优选方位(LPO)、形态优选方位(SPO)和定向微裂隙。随着糜棱岩化作用的增强,岩石的P波各向异性逐 渐升高。变形岩石中的黑云母、石英、绿辉石的LPO和SPO是地震波各向异性的主要控制因素。饱水后的片麻岩样品的P波各向异性明显低于干燥片麻岩样品。在东海钻井中的强反射带主要是由于不同岩层之间的波阻抗差异而造成的,榴辉岩/强退变榴辉岩和黑云斜长片麻岩之间的接触界面会产生较强的地震深反射。此外,与LPO相关的地震波各向异性会增强地震波的反射,所以韧性剪切带中的糜棱岩化片麻岩可能是地震反射的良好载体。韧性剪切带中岩石弹性波速度的强各向     

TRANSALP—Cross-line recording during the seismic reflection transect in the Eastern Alps

Cross-line recording formed a companion experiment of the TRANSALP seismic reflection transect through the Eastern Alps, conducted by partner institutions from Austria, Germany and Italy in three field campaigns in the period fall 1998 to fall 1999. Besides of the originally expected three-dimensional control for the north–south running main transect, additional information on seismic anisotropy and alternative images of crucial parts of the main transect could be gained.Conventionally processed sections along N–S running common-midpoint (CMP) binning lines confirm and strengthen the predominance of midcrustal reflective structures of the ‘Sub-Tauern-Ramp’ beneath and south of the Tauern Window. Velocity analysis of the first arrivals exhibit about 10% higher velocities in east–west propagating P-waves, compatible with texture-dominated rock anisotropy, recorded on cross-lines at the Tauern Window. Pre-stack depth migration of cross-line recordings shows dominant south dip of the Sub-Tauern-Ramp with easterly dip components and a sub-horizontal root zone of the Sub-Dolomites-Ramp.     

Exploring the relative contribution of mineralogy and CPO to the seismic velocity anisotropy of evaporites

We present the influence of mineralogy and microstructure on the seismic velocity anisotropy of evaporites. Bulk elastic properties and seismic velocities are calculated for a suite of 20 natural evaporite samples, which consist mainly of halite, anhydrite, and gypsum. They exhibit strong fabrics as a result of tectonic and diagenetic processes. Sample mineralogy and crystallographic preferred orientation (CPO) were obtained with the electron backscatter diffraction (EBSD) technique and the data used for seismic velocity calculations. Bulk seismic properties for polymineralic evaporites were evaluated with a rock recipe approach. Ultrasonic velocity measurements were also taken on cube shaped samples to assess the contribution of grain-scale shape preferred orientation (SPO) to the total seismic anisotropy. The sample results suggest that CPO is responsible for a significant fraction of the bulk seismic properties, in agreement with observations from previous studies. Results from the rock recipe indicate that increasing modal proportion of anhydrite grains can lead to a greater seismic anisotropy of a halite-dominated rock. Conversely, it can lead to a smaller seismic anisotropy degree of a gypsum-dominated rock until an estimated threshold proportion after which anisotropy increases again. The difference between the predicted anisotropy due to CPO and the anisotropy measured with ultrasonic velocities is attributed to the SPO and grain boundary effects in these evaporites.     

P-wave anisotropy in eclogites and relationship to the omphacite crystallographic fabric

Measurements of P-wave velocity at room temperature and confining pressures up to 500 MPa were carried out on three eclogite mylonites collected from a shear zone in the Monviso area (Western Alps). P-wave velocities at a pressure of 400 MPa range from 7.7 km/s to 7.9 km/s, yielding to a maximum anisotropy of 6%. From the CPO of omphacite we estimated a maximum contribution of omphacite to the P-wave anisotropy of only 1.3%. These results suggest that primarily the compositional layering and secondary the fabric of minor constituent minerals significantly contribute to the seismic anisotropy. Because of the anisotropy, the seismic reflectivity of subduction zones may vary with the direction of observation.     

Calculations of upper-mantle P-wave velocity anisotropies

Discovery of anisotropy for upper-mantle P-wave velocities has resulted recently in many laboratory determinations of velocities for natural mantle material, either by direct measurement on rock samples or through calculations based on single-crystal elastic constants, fabrics and modal compositions Comparison of field and laboratory data provides gross constraints for convection processes and flow mechanisms. As natural peridotites yeidl a wide range of maximum seismic anisotropies, their depth of origin and deformation history are inferred from pyroxene geobarometry and geothermometry, and from stress estimates derived from neoblast grain sizes. Large maximum anisotropies (ca. 1.0 km/s) are typical of samples with fabrics strengthened by low-temperature/highstress deformation related to emplacement. Other highly anisotropic mantle material includes peridotites beneath continental rifts and thinned continental crust. Maximum anisotropies for ordinary suboceanic and infracontinental peridotites are generally in the range of the field data (0.3–0.6 km/s) and support a model with [010]_ol nearly vertical (in accord with observations on ophiolites) and [100]_ol subparallel to the spreading direction. On the basis of fabric data, seismic anisotropy should decrease with depth.     

Experimental Study on Acoustoelastic Character of Rock Under Uniaxial Compression

This paper made a research about the change rule of elastic wave velocity with stress applied on rock from theoretical and experimental aspect. Firstly, a mathematical model of P-wave velocity and confining pressure of rock was set up from the point of acoustoelastic character. Effect of axial stress on P-wave velocity in granite and sandstone during uniaxial compression process was studied experimentally by using GAW-2000 rock mechanical testing system and RSM-SY5 ultrasonic wave testing system, and the relation curves of axial stress with P-wave velocity were obtained. Based on test data, acoustoelastic theoretical formulas of granite and sandstone were established and the best empirical formulas were fitted by using regression method. Meanwhile, a comparative analysis of the empirical and theoretical calculated values was carried out. Finally, the reliability of applying acoustoelastic theoretical formula in hard rock range was further verified based on the experimental data of granitic gneiss. The results show that the P-wave velocity experiences a rapid increase, gentle increase and then a sharp fall during the uniaxial compression process. The B-value in acoustoelastic theoretical formula (proportion coefficient determined by elastic modulus and third-order elastic constant) decreases exponentially with axial stress. The acoustoelastic theoretical formula can effectively reflect the relationship between rock acoustic velocity and stress within the allowable error, which can be the theoretical foundation of acoustoelastic geo-stress measurement of subsurface rock mass.     

Geochemistry, petrofabrics and seismic properties of eclogites from the Chinese Continental Scientific Drilling boreholes in the Sulu UHP terrane, eastern China

We present an integrated study of geochemistry, petrofabrics and seismic properties of strongly sheared eclogites from the Chinese Continental Scientific Drilling (CCSD) project in the Sulu ultrahigh-pressure (UHP) metamorphic terrane, eastern China. First, geochemical data characterize diverse protoliths of the studied eclogites. The positive Eu- and Sr-anomalies, negative Nb anomaly and flat portion of heavy rare earth elements in coarse-grained rutile eclogites (samples B270 and B295) suggest a cumulate origin in the continental crust, whereas the negative Nb anomaly and enrichment of light rare earth elements in retrograde eclogites (samples B504, B15 and B19) imply an origin of continental basalts or island arc basalts. Second, P-wave velocities (V_p) of three typical eclogite samples were measured under confining pressures up to 500 MPa and temperatures to 700 °C. At 500 MPa and room temperature, the mean V_p reaches 8.50-8.53 km/s in samples B270 and B295 but drops to 7.86 km/s in sample B504, and the P-wave anisotropy changes from 1.7-2.7% to 5.5%, respectively. The pressure and temperature derivatives of V_p are larger in the retrograde eclogite than in fresh ones. Third, the electron backscatter diffraction (EBSD) measurements of the eclogites reveal random crystal preferred orientation (CPO) of garnet and pronounced CPO of omphacite, which is characterized by a strong concentration of [001]-axes sub-parallel to the lineation and of (010)-poles perpendicular to the foliation. The asymmetric CPO of omphacite in sample B270 recorded a top-to-the-south shear event during subduction of the Yangtze plate. The calculated fastest V_p is generally sub-parallel to the lineation, but a different deformation environment during exhumation could form second-order variations in omphacite CPO and affect the V_p distribution in eclogites (e.g., the fastest V_p is at ~ 35° from the foliation in sample B295). Comparison between measured and calculated seismic properties indicates that the CPO of omphacite controls the seismic anisotropy of eclogites at high pressure, and compositional layering and retrograde minerals will increase the anisotropy. Calculated P-wave velocities agree well with velocities measured at 500 MPa and room temperature for fresh eclogites, but much higher than those of retrograde eclogite. As a case study, the laboratory-derived V_p-P and V_p-T relationships were used to estimate P-wave velocities of eclogites and peridotites beneath the Western Superior Province, Canada. The results indicate that besides the fabric-induced anisotropy, the direction dependence of pressure and temperature derivatives of V_p can significantly increase seismic anisotropy of eclogites with depth, which results in eclogites being an important candidate for the seismic anisotropy in the upper mantle. Due to their very high density and velocity, garnet-rich eclogites within peridotite could be detected in seismic reflections in subduction zones.     

Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.     

Orogenic structure of the Eastern Alps, Europe, from TRANSALP deep seismic reflection profiling

The TRANSALP Group, comprising of partner institutions from Italy, Austria and Germany, acquired data on a 340 km long deep seismic reflection line crossing the Eastern Alps between Munich and Venice. Although the field work was split into four campaigns, between fall 1998 and summer 2001, the project gathered for the first time a continuous profile across the Alps using consistent field acquisition and data processing parameters. These sections span the orogen itself, at its broadest width, as well as the editor Fred Davey and the two adjacent basins. Vibroseis and explosion data, complementary in their depth penetration and resolution characteristics, were obtained along with wide-angle and teleseismic data. The profile shows a bi-vergent asymmetric structure of the crust beneath the Alpine axis which reaches a maximum thickness of 55 km, and 80–100 km long transcrustal ramps, the southward dipping ‘Sub-Tauern-Ramp’ and the northward-dipping ‘Sub-Dolomites-Ramp’. Strongly reflective patterns of these ramps can be traced towards the north to the Inn Valley and towards the south to the Valsugana thrust belt, both of which show enhanced seismicity in the brittle upper crust. The seismic sections do not reveal any direct evidence for the presence of the Periadriatic Fault system, the presumed equivalent to the Insubric Line in the Western Alps. According to our new evolutionary model, the Sub-Tauern-Ramp is linked at depth with remnants of the subducted Penninic Ocean. The ‘crocodile’-type model describes an upper/lower crustal decoupling and wedging of both the European and the Adriatic–African continents.     

The effect of high temperature and high confining pressure on compressional wave velocities in quartz-bearing and quartz-free igneous and metamorphic rocks

Compressional wave velocities have been measured in granite, granulite, amphibolite and peridotite specimens under conditions of high temperature up to 700°C and confining pressures up to 6 kbar. In general, velocity increases with pressure and decreases with temperature.Quartz-bearing rocks show an anomalous behavior of their compressional wave velocities. The velocity—temperature relations exhibit a velocity-“deep” due to the high—low inversion of the constituent quartz crystals. The intrinsic effect of temperature on velocities is hard to determine due to thermal expansion and consequent loosening of the structure. The opening of new cracks and the widening of old cracks causes a large decrease in compressional wave velocities. The minimum pressure to prevent damage at a given temperature should, therefore, be about 1 kbar/100°C.The values obtained at these conditions are considered to be most nearly correct as intrinsic properties of the compact aggregates. Velocity anisotropies at high confining pressures and high temperatures correlate with preferred lattice orientation of the constituent minerals. The effect of dimensional orientation and microcracks on seismic anisotropy seems to be of minor importance in dry rocks. It is the more eliminated the higher the confining pressure. The data do not support the concept of a velocity maximum in depth of 10–20 km.     

Experimental study of the seismic properties of the Eastern Alps (Italy) along the Aurina–Tures–Badia Valleys transect

In order to better constrain the interpretation and the nature of the seismic reflectors, experimental measurements at high confining pressure (up to 300 MPa) and room temperature of the compressional wave velocity (Vp) on 10 samples representative of the most common lithologies along the Aurina (Ahrntal), Tures (Tauferer Tal), and Badia (Abtei Tal) Valleys profile (Eastern Alps, Italy) have been performed. For each sample, the speed of ultrasonic waves was measured in three mutually perpendicular directions, parallel and normal to the rock foliation and lineation.The main results are:(a) Good agreement between the calculated vs. measured modal compositions of the considered rocks, indicating that they were presumably equilibrated at the estimated P–T conditions; therefore, the seismic properties are representative of the crustal level indicated by the thermobarometry.(b) Measured and calculated average Vp are in good agreement, and are typical of mid-crustal level (6.0–6.5 km/s). Only the amphibolites show Vp typical of the lower crust (7.2 km/s).(c) The seismic anisotropy of metapelites is very high (12–27%), both with orthorhombic and transverse isotropy symmetry; amphibolites are transversely isotropic with an anisotropy of 8%; orthogneisses and granitoids are isotropic or weakly anisotropic.(d) The contacts between amphibolites and all other rock types may generate good reflections, provided they are not steeply inclined. Although the metamorphic foliation remains steeply inclined, discordant buried sub-horizontal igneous contacts may be detected.     

中国大陆科学钻探主孔100～2000m岩石弹性波速度：对地震深反射的约束

在常温常压条件对中国大陆科学钻CCSD主孔岩心的700样品进行了弹性波速度测量，并建立了主孔2000m的波速(Vp和Vs)连续剖面，为检验地球物理模型的合理解释提供了岩石物理学方面的宝贵资料。主孔中新鲜榴辉岩纵波速度(Vp)最大(7.86km／s)，正副片麻岩波速最小，又分别为5.53km／s和5．71km／s，榴辉岩的波速随着退变质作用的增强而明显减小。主孔2000m总平均Vp速度为6．2km／s，它与地球物理探测方法获得的大别-苏鲁造山带上地壳具有6.2-6．3km／s高速层结论是一致的。大部分岩石具有明显地震波各向异性。水饱和度使岩石纵波(Vp)速度和剪切波速度(Vs)分别增加19％和6％，而使Vp的各向异性降低3％～4％。不同岩性界面的反射系数(Rc)是产生地震反射的主要原因。金红石榴辉岩与片麻岩之间具有很高的反射系数(0.24-0．31)。韧性剪切带中糜棱岩化片麻岩和面理化榴辉岩使岩石各向异性和反射强度明显增加。岩石微裂隙与主孔原位波速变化有密切关系。饱水岩石速度(Vp和Vs)可以代表CCSD主孔原位状态的地震波速度。上述成果为本区地震反射体成因提供了重要的岩石物理性质约束。     

Crustal properties from S-wave and gravity data along a seismic refraction profile in Romania

VRANCEA'99 is a seismic refraction line that was carried out in 1999 to investigate the deep structure and physical properties of the upper lithosphere of the southeastern Carpathians and its foreland. It runs from the city of Bacau to the Danube River, traversing the Vrancea epicentral area of strong intermediate-depth seismicity and the city of Bucharest.

Interpretation of P-wave arrivals led to a velocity model that displays a multi-layered crust with velocities increasing with depth. The range of P wave velocities in the sedimentary cover increases from N to S and a structuring of the autochthonous basement of the Moesian Platform is observed. The crystalline crust displays thickness variations, but at the same time the lateral velocity structure along the seismic line remains almost constant. An intra-crustal boundary separates an upper crust from the lower crust. Within the upper mantle a low velocity zone is detected at a depth of about 55-km.

The interpretation of observable S-waves resulted in a velocity model that shows the same multi-layered crust, with S-velocities increasing similarly with depth as the P-waves. The corresponding Poisson's ratio is highly variable throughout the crust and ranges from 0.20–0.35 for the sedimentary cover to 0.22–0.25 for the crystalline crust. The interpretation of the V_p, V_s and Poisson's ratio in petrological terms suggests a large variety of rocks from sand and clay to sandstone, limestone and dolomite within the sedimentary cover. Within the crystalline crust the most probably rock types are granite, granodiorite, granite–gneiss and/or felsic amphibolite–gneiss in the upper part and gneiss and /or amphibolite in the lower part.

Based on the 2-D seismic velocity model, a density model is developed. Density values are assigned to each layer in agreement with the P-wave velocity model and with values accepted for the geological units in the area. After several iterations a good fit between the computed and observed Bouguer anomalies was obtained along the seismic line.     

Preliminary investigation results on fabrics and related physical properties of an anisotropic gneiss

In order to evaluate the fabric-dependent anisotropy of a particular gneiss type, we assessed the quartz lattice and shape preferred orientations as well as the microcrack pattern statistics. In the rock mechanics laboratory, several strength and deformability tests on drilled rock samples were run and the seismic wave propagation properties along the rock’s principal strain axes were determined. In the mechanical tests it turned out that despite a distinct stretching lineation with initially extreme grain elongation, the rock performs for the larger part as a transversely isotropic material with the schistosity as plane of isotropy. The anisotropy in seismic wave propagation within the plane of schistosity is attributed primarily to a highly anisotropic microcrack pattern, the distinct gneissose banding and to a lower degree to the LPO of quartz.     

苏鲁超高压变质岩中地震波速随围压的变化规律

我们在10~800MPa的静水压力下测量了中国大陆科学钻探(CCSD)主孔岩心和采自苏鲁地区地表露头上的共68块典型超高压变质岩试样的地震波速,并建立了波速和围压的定量关系,成功地解释了描述这一定量关系方程中4个参数的物理意义。笔者相信,本文将为超高压变质岩地区以及新老俯冲带地震波资料的解释提供必不可少的理论和实验支撑。     

NEHRP soil classification and estimation of 1-D site effect of Dehradun fan deposits using shear wave velocity

Shear wave velocity of the near surface soil at nearly 50 sites in the sub Himalayan mountain exit covering Doon fan deposits, was determined using Multi-channel Analysis of Surface Waves (MASW), a seismic reflection technique. Based on the average shear wave velocity of the upper 30 m soil column, sites in the Dehradun fan are predominantly classified as class ‘D’ (180–360 m/s). Similarly, sites located on the northwestern, eastern and southeastern sides of the fan deposit have shear wave velocities (in the upper 30 m soil) greater than 360 m/s, thereby classifying them as class ‘C’ (360–760 m/s) in accordance with NEHRP provisions. Some of the sites towards the southwestern side of the fan deposits had average shear wave velocities less than 180 m/s and could be classified as soil class ‘E’. One dimensional site effects, including amplification and dynamic period were calculated for the majority of the sites. However, some of the representative suite of sites across the north–south profile of Dehradun fan has been discussed here. Although the attenuation is greater on the southwestern side of the Dehradun fan deposits (i.e. thicker, low velocity sediments) and the sites had been classified as class ‘D’ and ‘E’ but the site amplification tends to be greater in the northern and northwestern part of the city due to large impedance contrast with in the near surface soils.     

Mineralogical and seismic properties of serpentinite of Ait Ahmane fault zone of Bou Azzer ophiolite,central Anti-Atlas of Morocco

This paper presents the first seismic measurements of serpentinite of Bou Azzer ophiolite, central Anti-Atlas of Morocco, including seismic velocities and anisotropy. Two serpentinite samples collected from the Ait Ahmane fault zone were analyzed in order to define the mineralogical and seismic features of the natural serpentinite of the Bou Azzer ophiolite. The mineralogical features were investigated using microscopic observation and Raman spectroscopy, while the seismic features were performed using an Electron Backscatter Diffraction (EBSD) instrument. Microscopy and spectroscopy analyses confirmed that the investigated serpentinite suffers from a variable degree of serpentinization, and the antigorite is the dominant variety of serpentine minerals in the study area. The crystal preferred orientation (CPO) results show that the axis [001] of the antigorite deformation is aligned subnormal to the foliation, while the axis [010] is mostly aligned subparallel to the lineation. The seismic anisotropy results are depending on serpentine amount in the rock samples. The sample with a low serpentine amount (30%) shows lowest P- and S-wave anisotropy (Vp = 7.2% and AVs = 6.55%), while the sample with a high amount of serpentine (85%) presents highest P-wave and S-wave anisotropy (Vp = 8.6% and AVs = 11.06%). Consequently, the results indicate that seismic anisotropy increases when increasing the antigorite amount.     

Fabric and seismic properties of carrara marble mylonite

Shear deformation in calcite-rich rocks can produce strong lattice preferred orientations (LPO), which result in a high anisotropy of bulk seismic properties because of the high elastic anisotropy of calcite (32% Vp anisotropy). Deformed rocks often show also strong shape preferred orientations (SPO). Theories for averaging the elastic properties have not yet satisfactorily predicted the contribution to the seismic anisotropy caused by the SPO alone.A calcite mylonite from Carrara (Italy) was investigated, which is characterised by a strong SPO and a weak LPO. It was composed of about 80% calcite, then white mica, quartz and hematite. Flattening of mica and of calcite grains defined the mylonitic foliation, and elongation of calcite grains defined the lineation. On average calcite grains have aspect ratios of about 2.5:1.6:1, and grain sizes of about 10 μm. At 400 MPa confining pressure, the measured Vp (km/s) parallel to the lineation (X direction) was highest (6.63), lower in the intermediate Y direction (6.47); the Vp normal to the foliation (Z direction) was lowest (6.30). This yielded a Vp anisotropy of 5%. The LPO, determined by automated electron backscatter diffraction (EBSD), was very weak (texture index 1.1), with intensities between 0.6 and 1.6 m.r.d. in the c-axis pole figure. Extrapolation of the texture index to an infinite number of orientation measurements indicated that the observed variations were mostly random noise in the orientation distributions and that the bulk rock texture was random. The Vp anisotropy of the Voigt, Reuss and Hill averages calculated from this calcite LPO is predicted to be close to zero. Adding 5% of muscovite with (001) perfectly aligned parallel to the foliation, we calculated a total anisotropy of 2.8%. The anisotropy calculated for the special directions X, Y and Z remained at 2.6% only.It was concluded that the measured seismic anisotropy cannot be explained by the LPO of calcite and by 5% of mica alone. It is also attributed to the strong SPO and to further grain boundary effects.     

Three-dimensional P-wave velocity image under the Carpathian arc

An inversion of P-wave travel time residuals from selected earthquakes in the distance range 30°–98° to two seismic station networks was used to model P-wave velocity anomalies down to 250 km depth. In the first inversion experiment a region between 43.5°–47.5°N and 21°–29°E was modelled, using 35 seismic stations, while in the second one a region between 44°–47°N and 25°–29°E was modelled, using 19 seismic stations. The 4-layer block model of the first inversion offers 19% reduction in residual variance, while the 5-layer block model of the second one offers 26% reduction, the rest being explained by noise and smaller scale heterogeneities. The obtained velocity anomalies correlate remarkably well with the gravity anomalies and with the tectonic model for the Vrancea region of Fuchs et al. (1979).