Inversion of seismic and gravity data for the composition and core sizes of the Moon

We model the internal structure of the Moon, initially homogeneous and later differentiated due to partial melting. The chemical composition and the internal structure of the Moon are retrieved by the Monte-Carlo inversion of the gravity (the mass and the moment of inertia), seismic (compressional and shear velocities), and petrological (balance equations) data. For the computation of phase equilibrium relations and physical properties, we have used a method of minimization of the Gibbs free energy combined with a Mie-Gr@uneisen equation of state within the CaO-FeO-MgO-Al₂O₃-SiO₂ system. The lunar models with a different degree of constraints on the solution are considered. For all models, the geophysically and geochemically permissible ranges of seismic velocities and concentrations in three mantle zones and the sizes of Fe-10%S core are estimated. The lunar mantle is chemically stratified; different mantle zones, where orthopyroxene is the dominant phase, have different concentrations of FeO, Al₂O₃, and CaO. The silicate portion of the Moon (crust + mantle) may contain 3.5–5.5% Al₂O₃ and 10.5–12.5% FeO. The chemical boundary between the middle and the lower mantle lies at a depth of 620–750 km. The lunar models with and without a chemical boundary at a depth of 250–300 km are both possible. The main parameters of the crust, the mantle, and the core of the Moon are estimated. At the depths of the lower mantle, the P and S velocities range from 7.88 to 8.10 km/s and from 4.40 to 4.55 km/s, respectively. The radius of a Fe-10%S core is 340 ± 30 km.     

Composition,temperature, and thickness of the lithosphere of the Archean Kaapvaal craton

A new model is proposed for the structure of the Kaapvaal craton lithosphere. Based on chemical thermodynamics methods, profiles of the chemical composition, temperature, density, and S wave velocities are constructed for depths of 100–300 km. A solid-state zone of lower velocities is discovered on the S velocity profile in the depth interval 150–260 km. The temperature profiles are obtained from absolute values of P and S velocities, taking into account phase transformations, anharmonicity, and anelastic effects. The examination of the sensitivity of seismic models to the chemical composition showed that relatively small variations in the composition of South African xenoliths result in lateral temperature variations of ～200°C. Inversion of some seismic profiles (including IASP91) with a fixed bulk composition of garnet peridotites (the primitive mantle material) leads to a temperature inversion at depths of 200–250 km, which is physically meaningless. It is supposed that the temperature inversion can be removed by gradual fertilization of the mantle with depth. In this case, the craton lithosphere should be stratified in chemical composition. The depleted lithosphere composed by garnet peridotites exists to depths of 175–200 km. The lithospheric material at depths of 200–250 km is enriched in basaltoid components (FeO, Al₂O₃, and CaO) as compared with the material of garnet peridotites but is depleted in the same components as compared with the fertile substance of the underlying primitive mantle. The material composing the craton root at a depth of ～275 km does not differ in its physical and chemical characteristics from the composition of the normal mantle, and this allows one to estimate the thickness of the lithosphere at 275 km. The results of this work are compared with data of seismology, thermal investigations, and thermobarometry.     

The nature of the boundary between the upper and the lower mantle and its influence on convection

In the PREM seismic model, the boundary between the upper and the lower mantle is accepted at a depth of 670 km, where seismic velocities and density increase. However, until recently there was an obvious inconsistency in this model. The density increases abruptly, and the velocities, in addition to the jumps, have also the subsequent zones of increased gradient. The discontinuity between the upper and the lower mantle is related to the transition of olivine from the ringwoodite phase into the mixture of perovskite and magnesiowustite. However, in the pyrolyte model, the transition zone of the upper mantle consists not wholly of olivine, but partly of olivine (60%) and partly of garnet (40%). The latest data of the garnet measurement at high pressures show that it also experiences phase transition, being converted into magnesium perovskite with the impurity of calcium perovskite. In contrast to the sharp transition in olivine (within a depth interval of only 5 km), the transition in garnet is spread over the interval of depths of 660–710 km. In the widely used PREM and AK135 models, this additional transition corresponds to the zone of the increased gradient in seismic velocities, while in the density distribution it is included in the sharp transition of ringwoodite. Thus, the mineralogy data indicate the need for correction of the PREM and AK135 seismic models: the density jump at a depth of 660 km should be reduced by approximately a factor of two, and a subjacent layer with the increased density gradient should be added at the depth interval of 660–710 km. The phase transition in olivine hampers the mantle flows, although in garnet it accelerates them. Therefore, with an allowance for the smaller jump in density, the decelerating effect of the subducting plates, caused by the phase transition in olivine, decreases, and, furthermore, the effect of their acceleration, caused by the phase transition in garnet, is added. The decrease in the density jump by almost a factor of two will lead to essential changes in the results of the majority of recent works addressing the assessment of the deceleration of convection at the upper/lower mantle discontinuity on the basis of the PREM model.     

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.     

High-pressure phase equilibria in a garnet lherzolite,with special reference to Mg2+Fe2+ partitioning among constituent minerals

Phase equilibria in a natural garnet lherzolite nodule (PHN 1611) from Lesotho kimberlite and its chemical analogue have been studied in the pressure range 45–205 kbar and in the temperature range 1050–1200°C. Partition of elements, particularly Mg²⁺Fe²⁺, among coexisting minerals at varying pressures has also been examined. High-pressure transformations of olivine(α) to spinel(γ) through modified spinel(β) were confirmed in the garnet lherzolite. The transformation behavior is quite consistent with the information previously accumulated for the simple system Mg₂SiO₄Fe₂SiO₄. At pressures of 50–150 kbar, a continuous increase in the solid solubility of the pyroxene component in garnet was demonstrated in the lherzolite system by means of microprobe analyses. At 45–75 kbar and 1200°C, the Fe²⁺/(Mg + Fe²⁺) value becomes greater in the ascending order orthopyroxene, Ca-rich clinopyroxene, olivine and garnet. At 144–146 kbar and 1200°C, garnet exhibits the highest Fe²⁺/(Mg + Fe²⁺) value; modified spinel(β) and Ca-poor clinopyroxene follow it. When the modified spinel(β)-spinel(γ) transformation occurred, a higher concentration of Fe²⁺ was found in spinel(γ) rather than in garnet. As a result of the change in the Mg²⁺Fe²⁺ partition relation among coexisting minerals, an increase of about 1% in the Fe₂SiO₄ component in (Mg,Fe)₂SiO₄ modified spinel and spinel was observed compared with olivine.These experimental results strongly suggest that the olivine(α)-modified spinel(β) transformation is responsible for the seismic discontinuity at depths of 380–410 km in the mantle. They also support the idea that the minor seismic discontinuity around 520 km is due to the superposition effect of two types of phase transformation, i.e. the modified spinel(β)-spinel(γ) transformation and the pyroxene-garnet transformation. Mineral assemblages in the upper mantle and the upper half of the transition zone are given as a function of depth for the following regions: 100–150, 150–380, 380–410, 410–500, 500–600 and 600–650 km.     

CPO-induced seismic anisotropy in UHP eclogites

Ultrahigh-pressure (UHP) eclogites often show strong plastic deformation and anisotropy of seismic properties. We report in this paper the seismic velocity and anisotropy of eclogite calculated from the crystallographic preferred orientations (CPOs) of constituent minerals (garnet, omphacite, quartz and rutile) and single crystal elastic properties. We also compared the calculated results with the measured results in similar eclogites. Our results suggest that (1) Except that garnet is a seismically quasi-isotropic mineral, omphacite, quartz, coesite and rutile all have strong seismic anisotropies (AVp = 23.0%―40.9%, Max. AVs = 18.5%―47.1%). They are the major sources for anisotropy in eclogite. The average seismic velocities are fast in garnet and rutile, moderate in omphacite and coesite, and slow in quartz. (2) The deformed eclogites have the maximum Vp (8.33―8.75 km/s) approximately parallel to foliation and lineation, the minimum Vp (8.25―8.62 km/s) approximately normal to foliation and lineation and the Vp anisotropies of 1.0―1.7%. Their Vs are 4.93―4.97 km/s. The corresponding maximum anisotropies (0.73%―1.78%) of Vs are at 45° to both foliation and lineation and the minimum anisotropies at positions normal to lineation on the foliation plane. The Vs1 polarization planes are approximately parallel to foliation. The mean Vp and Vs of eclogite under UHP peak metamorphism conditions (P = 3―5 GPa, T = 900―1100℃) are estimated to be 3.4%―7.2% and 6.3%―12.1% higher than those at ambient pressure and temperature conditions, respectively. (3) Omphacite component dominates the anisotropy of eclogite while garnet component reduces the anisotropy and increases the seismic velocities. Quartz component has a small effect on the anisotropy but reduces the seismic velocities of eclogite. The effect of rutile component is negligible on seismic properties of eclogite due to its trivial volume fraction. (4) The increase of volume fraction of omphacite in eclogite will reduce the seismic velocities and increase the anisotropy. Omphacitite has seismic velocities reduced by 6%―8% and anisotropies increased to 3%―4% compared to those of garnetite. Our results suggest that the seismic properties calculated with single crystal elastic properties and CPOs are equivalent to those measured in laboratory. Moreover, it provides insights into the mineral physical interpretations of eclogite seismic properties.     

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.     

Testing the reference Moon model in respect of the thermal regime and chemical composition of the mantle: Thermodynamics versus seismology

The VPREMOON seismic reference Moon model (Garcia et al., 2011) has been tested with respect to the thermal regime and chemical composition of the mantle. Based on a self-consistent thermodynamic approach and petrological models of the lunar mantle covering a wide range of concentrations of CaO, Al₂O₃, and FeO, we convert the P- and S-wave velocity profiles to the temperature–depth profiles. The solution procedure relies on the method of the Gibbs free energy minimization and the equations of state for the mantle material which take into account the effects of phase transformations, anharmonicity, and anelasticity. We find that regardless of the chemical composition, the positive P- and S-wave velocity gradient in the lunar mantle leads to a negative temperature gradient, which has no physical basis. For adequate mantle temperatures, the P- and S-wave velocities should remain almost constant or slightly decrease with depth (especially V_S) as a result of the effects of the temperature, which grows faster than pressure. These findings underscore the importance of the relationship of the thermodynamics and physics of minerals with seismology.     

Temperature dependence of the elastic moduli of ringwoodite

The elastic moduli of polycrystalline ringwoodite, (Mg_0.91Fe_0.09)₂SiO₄, were measured up to 470 K by means of the resonant sphere technique. The adiabatic bulk (K_S) and shear (μ) moduli were found to be 185.1(2) and 118.22(6) GPa at room temperature, and the average slopes of dK_S/dT and dμ/dT in the temperature range of the study were determined to be −0.0193(9) and −0.0148(3) GPa/K, respectively. Using these results, we estimate seismic wave velocity jumps for a pure olivine mantle model at 520 km depth. We find that the jump for the S-wave velocity is about 1.5 times larger than that for the P-wave velocity at this depth. This suggests that velocity jumps at the 520 km discontinuity are easier to detect using S-waves than P-waves.     

Source parameters of the earthquakes of the Baikal rift system

The dynamic parameters of the earthquake source—the seismic moment, the moment magnitude, the source radius, the stress drop, and the amplitude of displacement—are determined by the amplitude Fourier spectra of the body shear waves (S-waves) for 62 earthquakes of the Baikal rift system with the energy class of K _P = 9.1–15.7. In the calculations I used the classical Brune model. The seismic moment of the earthquakes being investigated changes from 3.65 × 10¹¹ N m to 1.35 × 10¹⁸ N m, and the radii of earthquake sources vary from 390 m to 1.84 km. The values of the drop in stress Δσ grow with an increase in the seismic moment up to 1.7 × 10⁸ Pa. For the group of weak earthquakes (M _w = 1.7–3.3), extremely low values of the drop in stress 10³–10⁴ Pa are observed. The maximum amplitude of displacement in the source amounts to 5.95 m. The empirical equations between the seismic moment and the other dynamic parameters of the source are determined. The regional dependence of the seismic moment and energy class is obtained: log M ₀ ± 0.60 = 1.03K _P + 3.17. The character of the relationship between the seismic moment and the corner frequency indicates that the classical scaling law of the seismic spectrum for the earthquakes in question is not fulfilled. The obtained estimates of the dynamic parameters are in satisfactory agreement with the published data concerning the analogous parameters of the other rift zones, which reflects the general regular patterns of the destruction of the lithosphere and the seismicity in the extension zones of the lithosphere.     

Generation of komatiite magma and gravitational differentiation in the deep upper mantle

The densities of silicate liquids with basic, picritic, and ultrabasic compositions have been estimated from the melting curves of minerals at high pressures. Silicate liquids generated by partial melting of the upper mantle are denser than olivine and pyroxenes at pressures higher than 70 kbar, and garnet is the only phase which is denser than the liquid at pressures from 70 kbar to at least 170 kbar. In this pressure range, garnet and some fraction of liquid separate from ascending partially molten diapirs. It is therefore suggested that aluminium-depleted komatiite with a high Ca/OAl₂O₃ ratio may be derived from diapirs which originated in the deep upper mantle at pressures from 70 kbar to at least 140 kbar (200–400 km in depth), where selective separation of pyropic garnet occurs effectively. On the other hand, aluminium-undepleted komatiite is probably derived from diapirs originating at shallower depths (< 200 km). Enrichment of pyropic garnet is expected at depths greater than 200 km by selective separation of garnet from ascending diapirs. The 200-km discontinuity in the seismic wave velocity profile may be explained by a relatively high concentration of pyropic garnet at depths greater than 200 km.     

Mineral chemistry,P-T-t paths and exhumation processes of mafic granulites in Dinggye,Southern Tibet

Lower crustal high grade metamorphic rocks have been successively found at Pamirs nearby the western Himalayan syntaxis, Namjagbarwa and Dinggye nearby the eastern Himalayan syntaxis and the central segment of the Himalayan Orogenic Belt, respec-tively[1―4]. In particular, some researchers deduced that there were probably eclogites at some locations[5]. Moreover, some geochronological data of these lower crustal granulites also have been accumulated. For example, the high-pressure granulit…     

Petrological and seismic studies of the lithosphere in the earthquake swarm region Vogtland/NW Bohemia,central Europe

New petrological and geochemical data of upper mantle and lower crustal xenoliths from a Quaternary tephra deposit in Mýtina, Czech Republic, are discussed in the frame of previous geophysical results (receiver functions, reflection seismology) of the western Eger/Ohře Rift area. The Vogtland/NW Bohemia region is well known for intraplate earthquake swarms, which are usually associated with volcanic activity. As previously reported, ³He/⁴He data of CO₂ emissions in mofettes and mineral-water springs point at ongoing magmatic processes in this area. Using teleseismic P receiver functions, an approximately 40-km-wide Moho updoming (from 31 to 27 km) and indications for a seismic discontinuity at 50 to 60 km depth were observed beneath the active CO₂-degassing field. The studied xenolith suite probes a lithospheric profile within the structural and gas geochemical anomaly field of the western Eger Rift.With regard to texture, composition, p–T estimates and origin, five xenolith groups can be discriminated. Upper crustal xenoliths (quartzites, phyllites, mica schists) resemble crystalline country rocks at surface. One noritic xenolith (6 kbar, 800 °C) could represent a sample of the lower crust. Clinopyroxenites and hornblendites probably represent cumulates of the nephelinitic magma or fragments of magmatic veins. Porous wehrlites and one hornblende peridotite xenolith reflect a metasomatied upper mantle. Megacrysts of Ti-rich amphibole, olivine, clinopyroxene, and phlogopite could be fragments of pegmatitic veins or high-pressure phenocrysts. Most of the ultramafic nodules (xenoliths and megacrysts) formed at pressures between 6 and 11 kbar (22 to 38 km depth), at temperatures well above regional geotherms of the Bohemian Massif calculated from surface heat flow studies. Orthopyroxene-bearing spinel-lherzolite xenoliths were not observed. Our petrographical, geochemical, and thermobarometric results indicate a lithospheric mantle strongly altered by magmatic processes. This metasomatism can cause slower than typical uppermost-mantle seismic velocities in a greater area and might help to explain observed seismic anomalies.     

Thermal structure and thickness of the lithospheric mantle underlying the Siberian Craton from the kraton and kimberlit superlong seismic profiles

A self-consistent approach is proposed for the investigation of the thermal conditions, chemical composition, and internal structure of the upper mantle of the Earth. Using this approach, the thermal state of the lithospheric mantle beneath the Siberian Craton (SC) is reconstructed from P velocities, taking into account the phase transitions, anharmonicity, and the effects of anelasticity. The velocities of seismic waves are more sensitive to temperature than to the composition of the mantle rocks, which allows the velocity models to be effectively used for reconstruction of the thermal regime of the mantle. The temperature at depths 100–300 km is reconstructed by inversion of the Kraton and Kimberlit superlong seismic profiles for compositions of the garnet harzburgite, lherzolite, and intermediate composition of garnet peridotite. The averaged temperature in the normal continental mantle is reconstructed by inversion of the IASP91 reference model for depleted and fertile substance. One-dimensional models and two-dimensional thermal fields undergo a substantial fall in temperature (～300–600°C) beneath the Siberian Craton as compared to the temperatures of the continental mantle and paleotemperatures inferred from the thermobarometry of xenoliths. Temperature profiles of the Siberian Craton deduced from seismic data lie between the conductive geotherms of 32.5–40.0 mW/m² and below the P(H)-T values obtained for low- and high-temperature xenoliths from the Mir, Udachnaya, and Obnazhennaya kimberlite pipes. The thickness of the thermal lithosphere estimated from the intersection with the potential adiabat is 300–320 km, which is consistent with the data on heat flows and seismotomographic observations. This provides grounds for the assumption that the low-temperature anomalies (thermal roots of continents) penetrate down to a depth of 300 km. The analysis of the sensitivity of seismic velocity and density to the variations in temperature, pressure, and chemical and phase composition of petrological models shows that recognition of fine differences in chemical composition of the lithospheric rocks by seismic methods is impossible.     

Phase transformations in a harzburgite composition to 26 GPa: implications for dynamical behaviour of the subducting slab

The mineralogy adopted by a depleted harzburgite composition has been studied over the pressure interval 5–26 GPa at temperatures of 1300–1400°C. The pyroxene-garnet component of the harzburgite composition (harzburgite minus 82 wt.% olivine) transforms to majorite garnet by 18–19 GPa, and further disproportionates to the assemblage of garnet + stishovite + Mg₂SiO₄ spinel above 20 GPa. At still higher pressures, first ilmenite (22–24 GPa) and then perovskite MgSiO₃ (24–26 GPa) are found to coexist with garnet. Garnet disappears at 26 GPa and almost complete transition to perovskite is achieved at this pressure. The mineral proportions and density profiles in the subducting oceanic lithosphere, modelled by a combination of 80% harzburgite + 20% primitive MORB compositions are calculated as a function of depth under conditions isothermal with surrounding pyrolite mantle, and also for a temperature distribution in which the slab is substantially cooler than surrounding mantle to below 700 km. Under isothermal conditions, the slab has a density similar to surrounding mantle to a depth of 600 km. However, between 600 and 700 km, the slab is up to 0.08 g/cm³ denser than surrounding mantle. This is caused primarily by the higher alumina content in pyrolite as compared to harzburgite, which causes the transition to perovskite in pyrolite to occur at substantially higher pressures than in harzburgite. The presence of alumina also smears out the garnet-perovskite transition in pyrolite over a depth interval of 50 km, whereas this transformation is much sharper in the harzburgite composition. Calculations based on the observed phase equilibria also show that a subducted cool slab remains much denser (by 0.1–0.3 g/cm³) than surrounding mantle to a depth of 700 km but possesses a density similar to surrounding mantle below this depth. These results have important implications for the dynamical behaviour of slabs possessing different thermal regimes when they encounter the 670 km discontinuity and also for the nature of this discontinuity.     

Titanium solubility in coexisting garnet and clinopyroxene at very high pressure: the significance of exsolved rutile in garnet

Exsolution microstructures including ilmenite±garnet in clinopyroxene and rutile in garnet are common in clinopyroxenite and eclogite from the Sulu ultrahigh-pressure (UHP) terrane. In order to understand the phase relations and Ti solubility in both garnet and clinopyroxene in a natural TiO₂-bearing system, several experiments at 5-15 GPa, 1000-1400°C were carried out using the multianvil high-pressure apparatus. The Hujianlin ilmenite-rich garnet clinopyroxenite showing exsolution microstructure was selected as starting material, because it closely approaches a composition lying in the TiO₂-CaO-MgO-FeO-Al₂O₃-SiO₂ system. Except for minor melt in one experiment at 1400°C and 5 GPa, other run products contain majoritic garnet+clinopyroxene±ilmenite (or rutile) and exhibit neoblastic texture. With increasing pressure, Ti and Ca, Mg and Si contents of neoblastic garnet increase with decreasing Al. The principal coupled substitutions are Ca²⁺Ti⁴⁺→2Al³⁺ and Si⁴⁺Mg²⁺→2Al³⁺ responding to majorite component increase. Titanium solubility (0.8-4.5 wt% as TiO₂) in garnet and Grt_Ti/Cpx_Ti ratio have a pronounced positive correlation with pressure between 5 and 15 GPa. On the other hand, the coexisting clinopyroxene contains low Ti (0.17-0.53 wt% as TiO₂), and shows no significant pressure effect. Rutile exsolution in garnet is coupled to that of pyroxene exsolution; both are exsolved from majoritic garnet on decompression. Therefore, the amount of such exsolved lamellae is a potential indicator of high-pressure metamorphism in exhumed rocks, whereas the TiO₂ content of clinopyroxene coexisting with garnet is not sensitive to pressure change.     

The model of the depth variations in the properties and state of the rocks in the upper,middle, and lower continental crust of the Kola-Norwegian block,Kola Peninsula

Based on the experimental and calculated data, the model of depth variations in density, as well as the velocities of longitudinal (P) and transverse (S) waves, within the upper, middle, and lower crusts of the Kola-Norwegian block down to a depth of about 40 km is suggested. The variations in density and the velocities of the P- and S-waves are primarily caused by the changes in the mineral composition of the rocks. The relative reduction in the velocities of the P- and S-waves under the action of the increasing pressure and temperature in the depth interval from 5 to 37 km is estimated at ～2%.     

Synthesis of pyrope-knorringite solid solution series

The garnet solid solution series between pyrope Mg₃Al₂Si₃O₁₂ and knorringite Mg₃Cr₂Si₃O₁₂ has been synthesized from oxide mixtures at pressures of 60–80 kbars and 1400–1500°C. Lattice parameters and refractive indices of solid solutions vary linearly with (molecular) composition within the limits of measurement. The lattice parameter of pure knorringite is 11.600Åand its refractive index is 1.83. The genetic significance of mineral inclusions in natural diamonds is discussed, particularly in the light of the very high knorringite contents often found in garnet inclusions. It is suggested that the most common mineral assemblage occurring as inclusions in diamonds (olivine + knorringite-rich garnet + enstatite) might be explained in terms of subduction into the mantle of olivine + chrome-spinel + enstatite cumulates originally formed by crystallization of mafic magmas within the oceanic crust. The cumulate assemblage experienced alteration by circulating hydrothermal solutions, resulting in the introduction of some carbonate and serpentine minerals. During subduction, this assemblage was partially melted at depth below 150 km, accompanied by reduction of carbonate, to form a reconstituted assemblage consisting of olivine + knorringite-rich garnet + enstatite ± diamond.