Seismic Velocities and Anisotropy of the Lower Continental Crust: A Review

—Seismic anisotropy is often neglected in seismic studies of the earth’s crust. Since anisotropy is a common property of many typically deep crustal rocks, its potential contribution to solving questions of the deep crust is evaluated. The anisotropic seismic velocities obtained from laboratory measurements can be verified by computations based on the elastic constants and on numerical data pertaining to the texture of rock-forming minerals. For typical lower crustal rocks the influence of layering is significantly less important than the influence of rock texture. Surprisingly, most natural lower crustal rocks show a hexagonal type of anisotropy. Maximum anisotropy is observed for rocks with a high content of aligned mica. It seems possible to distinguish between layered intrusives and metasediments on the basis of in situ measurements of anisotropy, which can thus be used to validate different scenarios of crustal evolution.     

THEORETICAL AND EXPERIMENTAL INVESTIGATIONS FOR CAVITY RESEARCH WITH GEOELECTRICAL RESISTIVITY METHODS*

The exact localization of subterranean cavities and the determination of their dimensions is very important for the planning of geotechnical and mining activities. It is a complicated geophysical task often at the limit of detection. Nevertheless geophysical investigation is the only alternative to a dense and expensive grid of boreholes. This report tests the usefulness of geoelectrical resistivity methods for cavity detection under some new aspects. The basis for evaluation was a theoretical analysis of different conventional and focussing measuring arrays and of special arrays for a geoelectrical research between two boreholes. The limit of detectability of a cylindrical cavity of defined cross-section and depth was calculated for the different measuring arrays on the basis of computation of the apparent resistivity ρa. Furthermore, the influence of possible errors (current supply of the electrodes and the distance between the electrodes) is discussed for focussed systems. The second part of the article is directed at the behaviour of the apparent resistivity ρa, the disturbing potential δV_d caused by the cavity and the normal potential δV₀ of the measuring array all in relation to a homogeneous earth. Some new results are presented. In the last part of the article theoretical results are compared with some field measurements.     

Thermal effects of chemical reactions in the undifferentiated earth

Calculations of thermal effects of oxidation-reduction reactions in the open systems for the undifferentiated earth have been carried out up to 2000 kbar. It is shown that a change in the sign of the thermal effect occurs with increase of pressure in these open systems: endothermal reactions which occur at pressures P_s < P_c become exothermal at pressures P_s > P_c and vice versa. At a certain critical pressure P_s = P_c the enthalpy of reaction is equal to zero.The results of calculations show that the central part of the earth and the deep interior could become warm due to reduction processes (for example reduction of ferromagnesian silicates and oxides, and stishovite) during core formation. Reducing conditions in the undifferentiated earth during core formation imply that a large fraction of the total heat of chemical reactions was released in the deep interior and was absorbed in the upper parts of the planet.     

Phase transformations in MSiO4 compounds at high pressures and their geophysical implications

The phase behaviour of MSiO₄ compounds (MHf, Zr, U and Th0 has been investigated at high pressures and temperatures in a diamond-anvil press coupled with laser heating. All of these compounds have been found to undergo two or perhaps three phase transformations at pressures below 300 kbar. The high-pressure phase transformations of these compounds differ from one another, with the exception of HfSiO₄ and ZrSiO₄, which undergo identical phase transformations. The ultimate phase assemblages of these compounds are of dense component dioxides (although this is yet to be confirmed in ThSiO₄). It is suggested that the heat-producing elements U and Th would exist as dioxide solid solutions rather than silicates in the deep interior of the earth. Moreover, the densities of these dioxides are more than twice as great as mantle silicates and even slightly greater than pure iron under similar P, T conditions. Gravitational separation due to mandle convection may transport these dioxides to the deep interior of the earth to form deep heat sources. It is also suggested, however, that these deep heat sources are located in the inner-outer core boundary instead of in the lower mantle.     

Model studies of physical properties of mineral matter in high-pressure-temperature experiments

Data derived from experimental model studies of the elastic properties of rocks of various compositions under high pressures and temperatures inherent in the thermobaric regimes of the upper lithosphere are generalized. The geophysical implications of the data are discussed. New results of experimental studies of the electrical characteristics of metamorphic rocks of ultrabasic composition at P ? 20 kbar and T ? 1000°C are presented.     

Simultane interpretation geophysikalischer messdatensysteme

Summary After general definitions in1 and deducing a special method in2, it will be shown in3 on a simple concrete example the effectiveness of the simultaneous interpretation of two series of data belonging two different geophysical methods, namely, it can be determined in the simultaneous way deep values, too, although both methods, separately, are not capable under the given conditions even for estimation of the great order of the deep. In4 there are some conclusions; above all it is emphasized for any geophysical method the necessity and importance of the general applicable form of the interpretation procedures.     

指示地球深层动力过程的地震各向异性"化石"

地震各向异性是地球内部介质的一种重要属性。地震各向异性研究不仅有助于准确认识地球内部精细结构,还有助于研究地球活动过程。本文从海洋岩石圈与大陆岩石圈各向异性结构两个方面较详细地论述了地震各向异性与地球深层动力过程的相互关系,论证了地震各向异性是指示地球深层动力过程的地球物理“化石”。     

Study the elastic properties and the anisotropy of rocks using different machine learning methods

This paper aims to demonstrate that the elastic stiffnesses and the anisotropic parameters of rocks can be accurately predicted from geophysical features such as the porosity, the density, the compression stress, the pore pressure and the burial depth using relevant machine learning methods. It also suggests that the extreme gradient boosting method is the best method for this purpose. It is more accurate, extremely faster to train and more robust than the artificial neural networks and the support vector machine methods. Very high R-squared scores was obtained for the predicted elastic stiffnesses of a relevant dataset that is available in the literature. This dataset contains different types of rocks, and the values of the features are in large ranges. An optimal set of parameters was obtained by considering an appropriate sensitivity analysis. The optimized model is very easy to implement in Python for practical applications.     

Prediction of thermal conductivity in reservoir rocks using fabric theory

An accurate prediction of the thermal conductivity of reservoir rocks in the subsurface is extremely important for a quantitative analysis of basin thermal history and hydrocarbon maturation. A model for calculating the thermal conductivity of reservoir rocks as a function of mineral composition, porosity, fluid type, and temperature has been developed based on fabric theory and experimental data. The study indicates that thermal conductivities of reservoir rocks are dependent on the volume fraction of components (minerals, porosity, and fluids), the temperature, and the fraction of series elements (FSE) which represents the way that the mineral components aggregate. The sensitivity test of the fabric model shows that quartz is the most sensitive mineral for the thermal conductivity of clastic rocks. The study results indicate that the FSE value is very critical. Different lithologies have different optimum FSE values because of different textures and sedimentary structures. The optimum FSE values are defined as those which result in the least error in the model computation of the thermal conductivity of the rocks. These values are 0.444 for water-saturated clay rocks, 0.498 for water-saturated sandstones, and 0.337 for water-saturated carbonates. Compared with the geometric mean model, the fabric model yields better results for the thermal conductivity, largely because the model parameters can be adjusted to satisfy different lithologies and to minimize the mean errors. The fabric model provides a good approach for estimating paleothermal conductivity in complex rock systems based on the mineral composition and pore fluid saturation of the rocks.     

Inducing action of hydrogen migrating along faults on earthquakes

The reduction of the microhardness and the crystal constants of some non-metallic materials, such as calcite, dolomite, antigorite, etc., are observed after a short time of hydrogen permeating treatment at low pressure. It means that hydrogen diffusion can cause their strength dropping or weakening. The hydrogen, which is produced under the earth by various chemical reactions or accumulated when the earth formed, is migrating up continuously along faults, causing weakening of rocks and faults at the same time. So it is possible that rocks and faults break under lower tectonic stress condition. Hydrogen anomalies are passive reflection, precursor and accompaniment of fault activities or earthquakes on the face of it, but hydrogen migrating has active influence on faults and its moving. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 229–235, 1992.     

由地面、卫星重力资料研究岩石层密度

地球内部物质密度的分布决定了地球表面及外空的重力场,不过在特定的条件下,也可以根据地面或外空的重力场推求地球内部异常密度的分布。本文在单层位理论的基础上,利用上述资料确定了我国大陆下面地球深部的异常密度,试用分离重力场长、短波长的方法,求出了地球外层-岩石层内异常密度的侧向分布,对此,还初步作了地球物理的解释,其结果与地震测深的较为一致。     

Evaluation of approximations in modelling the cooling of magmatic bodies

The cooling of a magmatic intrusion is simulated by a simple model of a non-homogeneous earth, with thermal properties depending on temperature, in which heat transfer is assumed to take place by conduction only. The mathematical problem consists in solving a non-linear partial differential equation with continuity conditions on temperature and heat flux imposed at the contacts between different rocks. This has been done numerically by a finite difference method. The model is then adopted as “reality” against which a number of commonly used approximations are tested. It is found that the effect of latent heat liberation can be reasonably taken into account by attributing an effective initial temperature to the magma (errors within 20°C for t > 10⁵ years, when the temperature of the magma is still as high as 600°C); the effective specific heat approximation does not work as well. The dependence of thermal conductivity and specific heat on temperature may be eliminated by maintaining the errors within 30°C for t < 5 × 10⁵ years. The assumption that magma and country rocks have the same thermal properties allows an estimate of the temperature field in the host rocks with errors of 50°C at most. The assumption that all rocks have the same constant conductivity yields results that are far from “reality” (errors of 100–200°C even at shallow depth).     

Poisson's ratio behaviour in various crystalline rocks: application to the study of the Earth's interior

The study of Poisson's ratio (σ) behaviour in various crystalline rocks under different temperatures and pressures shows this parameter to depend upon the rock composition rather than upon P-T conditions. The results of this study are presented in the form of a comparison of σ(z) distributions within the consolidated crust and continental upper mantle and the specific variations of σ in crust and mantle rocks underlying the Voronezh crystalline massif (VCM). These investigations, which are based upon seismic and seismological data as well as high pressure experiments, should clarify in particular the composition and petrology of the Earth's interior.     

Elastic-anisotropic properties of garnet granulites from the lower crust of the belomorian mobile belt: Results of experimental study

The elastic-anisotropic properties of fifteen samples of garnet granulites from xenoliths in the diatreme on Elovyi Island, in Kandalaksha Bay in the White Sea are investigated. The PT-conditions of the formation of these rocks correspond to a depth of 25–40 km, which is the current depth of the lower crust in the region. According to acoustopolariscopy, the elastic anisotropy of these deep rocks is relatively low. The average seismic P and S velocities corresponding to the in-situ PT-conditions are 6.5–6.7 km/s and 3.7–3.8 km/s, respectively. It is shown that the garnet content in the deep-seated rocks can be estimated from the seismic data. According to the obtained results, the formation of the structural and textural pattern of the lower crustal rocks in the region was poorly affected by the paleotectonic factors.     

Modeling of mineralogical composition,density and elastic wave velocities in anhydrous magmatic rocks

We use the technique of direct minimization of the Gibbs free energy of the 8-component (K₂O-Na₂O-Fe₂O₃-FeO-CaO-MgO-Al₂O₃-SiO₂) multiphase system in order to determine the equilibrium mineral assemblages of rocks of different bulk chemical compositions equilibrated at various P-T conditions. The calculated modal compositions of rocks and experimental data on elastic moduli of single crystals are then used to calculate densities and isotropic elastic wave velocities of rocks together with their pressure and temperature derivatives. Sufficient accuracy of the calculations is confirmed by comparison with experimental data on the gabbro-eclogite transformation and precise ultrasonic measurements of elastic wave velocities in a number of magmatic and metamorphic rocks.We present calculated phase diagrams with isolines of density, elastic wave velocities, and their pressure and temperature derivatives for several anhydrous magmatic rocks, from granite to lherzolite. Density and elastic properties of rocks are controlled by their chemical compositions, especially the SiO₂ content, and by P-T of equilibration, and they increase with pressure due to mineral reactions changing mineral assemblages from plagioclase-bearing and garnet-free to garnetbearing and plagioclase-free. TheV _p -density correlation is high, and shows two clear trends: one for iron-poor ultramafic rocks and another for all the other rocks considered. Mineral reactions, which occur at high pressures, changeV _p and density of anhydrous magmatic rocks following the well-known Birch (or a similar) law.Felsic, intermediate and mafic rocks can be well distinguished in theV _p -V _p /V _s - diagram, although their values ofV _p can be close to one another. TheV _p -V _p /V _s -density diagrams together with calculated phase diagrams can serve as efficient instruments for petrologic interpretation of seismic velocities.     

Electrical Properties of Crustal and Mantle Rocks – A Review of Laboratory Measurements and their Explanation

The electrical properties of rocks and minerals are controlled by thermodynamic parameters like pressure and temperature and by the chemistry of the medium in which the charge carriers move. Four different charge transport processes can be distinguished. Electrolytic conduction in fluid saturated porous rocks depends on petrophysical properties, such as porosity, permeability and connectivity of the pore system, and on chemical parameters of the pore fluid like ion species, its concentration in the pore fluid and temperature. Additionally, electrochemical interactions between water dipoles or ions and the negatively charged mineral surface must be considered. In special geological settings electronic conduction can increase rock conductivities by several orders of magnitude if the highly conducting phases (graphite or ores) form an interconnected network. Electronic and electrolytic conduction depend moderately on pressure and temperature changes, while semiconduction in mineral phases forming the Earth’s mantle strongly depends on temperature and responds less significantly to pressure changes. Olivine exhibits thermally induced semiconduction under upper mantle conditions; if pressure and temperature exceed ~ 14 GPa and 1400 °C, the phase transition olivine into spinel will further enhance the conductivity due to structural changes from orthorhombic into cubic symmetry. The thermodynamic parameters (temperature, pressure) and oxygen fugacity control the formation, number and mobility of charge carriers. The conductivity temperature relation follows an Arrhenius behaviour, while oxygen fugacity controls the oxidation state of iron and thus the number of electrons acting as additional charge carriers. In volcanic areas rock conductivities may be enhanced by the formation of partial melts under the restriction that the molten phase is interconnected. These four charge transport mechanisms must be considered for the interpretation of geophysical field and borehole data. Laboratory data provide a reproducible and reliable database of electrical properties of homogenous mineral phases and heterogenous rock samples. The outcome of geoelectric models can thus be enhanced significantly. This review focuses on a compilation of fairly new advances in experimental laboratory work together with their explanation.     

Combined geothermal-geophysical models of the earth's crust and upper mantle for the European continent

An attempt is made to obtain a combined geophysical model along two regional profiles: Black Sea— White Sea and Russian Platform—French Central Massif. The process of the model construction had the following stages: 1. The relation between seismic velocity (V_p, km/s) and density (σ, g/cm³) in crustal rocks was determined from seismic profiles and observed gravity fields employing the trial and error method. 2. Relations between heat production HP (μW/m³), velocity and density were established from heat flow data and crustal models of old platforms where the mantle heat flow HF_M is supposed to be constant. The HF_M value was also determined to 11 ± 5 mW/m². 3. A petrological model of the old platform crust is proposed from the velocity-density models and the observed heat flow. It includes 10–12 km of acid rocks, 15–20 km of basic/metamorphic rocks and 7–10 km of basic ones. 4. Calculation of the crustal gravity effects; its substraction from the observed field gave the mantle gravity anomalies. Extensively negative anomalies have been found in the southern part of Eastern Europe (50–70 mgal) and in Western Europe (up to 200 mgal). They correlate with high heat flow and lower velocity in the uppermost mantle. 5. A polymorphic advection mechanism for deep tectonic processes was proposed as a thermal model of the upper mantle. Deep matter in active regions is assumed to be transported (advected) upwards under the crust and in its place the relatively cold material of the uppermost mantle descends. The resulting temperature distribution depends on the type of endogeneous regime, on the age and size of geostructure. Polymorphic transitions were also taken into account.     

The response of a horizontal conducting cylinder embedded in a uniform earth for uniform and line current sources

The response of a horizontal conducting cylinder embedded in a uniform conducting earth is studied using mathematical models of uniform and line current source excitation for the period range 10 to 10⁴ s. The line current source is located at heights ranging from 100–750 km above the surface of the earth. From the calculated results, it is shown that for periods greater than 10³ s the ratioE _x /H _y at the surface of the earth for localized fields, such as the auroral and equatorial electrojet normally situated at heights of about 100 km, is considerably different from that for a uniform source. The results presented also show that the magneto-telluric method of geophysical prospecting for ore bodies in regions of the electrojet may not be very practicable for periods exceeding 10³ s.     

龙门山构造带深部动力学过程与地表地质过程的耦合关系

发生在地球浅层的2008年汶川地震驱动了龙门山及前陆地区的地表同震垂向位移.根据冲断带-前陆盆地弹性挠曲模型理论,在进行弹性挠曲模拟反演的基础上,结合对深部地球物理特征(泊松比、电性结构)的分析,发现龙门山前陆盆地现今岩石圈有效弹性厚度(T_e)具有自东向西逐渐减薄的趋势,自川中地区的30~40 km减至龙门山地区的10~20 km.在对晚三叠世以来前陆盆地各阶段盆地结构进行刻画的基础上,进行弹性挠曲模拟反演,推断龙门山前陆盆地的前渊地区(四川盆地西部)岩石圈的T_e值自晚三叠世以来具有逐渐减薄的趋势.这可能与松潘一甘孜地块下方广泛存在的软流圈热物质对四川盆地西部岩石圈下部的长期加热而导致的熔融有关,反映了地球深部动力学过程与地球表层盆地演化之间的耦合关系.     

Continental basalts in the accretionary complexes of the South-west Japan Arc: Constraints from geochemical and Sr and Nd isotopic data of metadiabase

Abstract The Ryoke Belt is one of the important terranes in the South‐west Japan Arc (SJA). It consists mainly of late Cretaceous granitoid rocks, meta‐sedimentary rocks (Jurassic accretionary complexes) and mafic rocks (gabbros, metadiabases; late Permian–early Jurassic). Initial ?_Sr (+ 25– + 59) and ?_Nd (? 2.1–?5.9) values of the metadiabases cannot be explained by crustal contamination but reflect the values of the source material. These values coincide with those of island arc basalt (IAB), active continental margin basalt (ACMB) and continental flood basalt (CFB). Spiderdiagrams and trace element chemistries of the metadiabases have CFB‐signature, rather than those of either IAB or ACMB. The Sr–Nd isotope data, trace element and rare earth element chemistries of the metadiabases indicate that they result from partial melting of continental‐type lithospheric mantle. Mafic granulite xenoliths in middle Miocene volcanic rocks distributed throughout the Ryoke Belt were probably derived from relatively deep crust. Their geochemical and Sr–Nd isotopic characteristics are similar to the metadiabases. This suggests that rocks, equivalent geochemically to the metadiabases, must be widely distributed at relatively deep crustal levels beneath a part of the Ryoke Belt. The geochemical and isotopic features of the metadiabases and mafic granulites from the Ryoke Belt are quite different from those of mafic rocks from other terranes in the SJA. These results imply that the Ryoke mafic rocks (metadiabase, mafic granulite) were not transported from other terranes by crustal movement but formed in situ. Sr–Nd isotopic features of late Cretaceous granitoid rocks occurring in the western part of the Japanese Islands are coincident with those of the Ryoke mafic rocks. Such an isotopic relation between these two rocks suggests that a continental‐type lithosphere is widely represented beneath the western part of the Japanese Islands.