In situ measurements of sound velocities and densities across the orthopyroxene → high-pressure clinopyroxene transition in MgSiO3 at high pressure

Using acoustic measurement interfaced with a large volume multi-anvil apparatus in conjunction with in situ X-radiation techniques, we are able to measure the density and elastic wave velocities (V_P and V_S) for both ortho- and high-pressure clino-MgSiO₃ polymorphs in the same experimental run. The elastic bulk and shear moduli of the unquenchable high-pressure clinoenstatite phase were measured within its stability field for the first time. The measured density contrast associated with the phase transition OEN → HP-CEN is 2.6-2.9% in the pressure of 7-9 GPa, and the corresponding velocity jumps are 3-4% for P waves and 5-6% for S waves. The elastic moduli of the HP-CEN phase are K_S=156.7(8) GPa, G = 98.5(4) GPa and their pressure derivatives are K_S′=5.5(3) and G′ = 1.5(1) at a pressure of 6.5 GPa, room temperature. In addition, we observed anomalous elastic behavior in orthoenstatite at pressure above 9 GPa at room temperature. Both elastic wave velocities exhibited softening between 9 and 13-14 GPa, which we suggest is associated with a transition to a metastable phase intermediate between OEN and HP-CEN.     

Melting of forsterite Mg2SiO4 up to 15 GPa

The melting curve of forsterite has been studied by static experiment up to a pressure of 15 GPa. Forsterite melts congruently at least up to 12.7 GPa. The congruent melting temperature is expressed by the Kraut-Kennedy equation in the following form: T_m(K)=2163 (1+3.0(V₀ ? V)/V₀), where the volume change with pressure was calculated by the Birch-Managhan equation of state with the isothermal bulk modulus K₀ = 125.4 GPa and its pressure derivative K′ = 5.33. The triple point of forsterite-β-Mg₂SiO₄-liquid will be located at about 2600°C and 20 GPa, assuming that congruent melting persists up to the limit of the stability field of forsterite. The extrapolation of the previous melting data on enstatite and periclase indicates that the eutectic composition of the forsterite-enstatite system should shift toward the forsterite component with increasing pressure, and there is a possibility of incongruent melting of forsterite into periclase and liquid at higher pressure, although no evidence on incongruent melting has been obtained in the present experiment.     

A scaling law for K′0 for silicates with constant mean atomic mass

It is shown that Birch's formula for the isothermal pressure derivative of the isothermal bulk modulus, K′, can be used to generate reasonable values of K′₀ for a sequence of silicates with various ambient densities and constant mean atomic mass, m¯. The theory predicts values in fairly good agreement with experimental results, although there is a regrettable spread of experimental values of K′₀ for each solid. This first-order approximation theory for scaling between K′₀ and ?₀ is analogous to the law of corresponding states which scales K₀ and ?₀.     

Compression and phase behavior of solid CO2 to half a megabar

CO₂ has been investigated up to 514 kbar at23 ± 2°C by both optical and in situ X-ray diffraction studies using a diamond-anvil pressure cell. CO₂ solidifies in an unknown structure in the pressure range 5 to 23 kbar, and transforms to ordinary dry-ice structure above 23 kbar at room temperature. Isothermal compression data for dry ice have been obtained above about 24 kbar. These appear to be the first data at room temperature known in the literature. The data fitted to the Birch equation of state yieldK₀ = 29.3 ± 1.0kbar andK₀^′ = 7.8 assuming the volume of the hypothetical dry ice at zero-pressure and room temperature is 31.4 ± 0.2 cm³/mole. The isothermal bulk modulus(K₀) thus derived is consistent with the compression data and compressibilities for dry ice obtained at low temperatures using dilatometry and ultrasonic techniques, respectively, reported in the literature. By comparing shock-wave data for relevant materials, it is suggested that CO₂ is not likely to transform to one of the crystalline forms of SiO₂ which is otherwise expected from empirical grounds, but may instead decompose into C (diamond) + O₂, at high pressures.     

The thermodynamic properties of the earth's lower mantle

The thermodynamic properties of the lower mantle are determined from the seismic profile, where the primary thermodynamic variables are the bulk modulus K and density ρ. It is shown that the Bullen law (K ∝ P) holds in the lower mantle with a high correlation coefficient for the seismic parametric Earth model (PEM). Using this law produces no ambiguity or trade-off between ρ₀ and K₀, since both K₀ and K′₀ are exactly determined by applying a linear K?ρ relationship to the data. On the other hand, extrapolating the velocity data to zero pressure using a Birch-Murnaghan equation of state (EOS) results in an ambiguous answer because there are three unknown adjustable parameters (ρ₀, K₀, K′₀) in the EOS.From the PEM data, K = 232.4 + 3.19 P (GPa). The PEM yields a hot uncompressed density of 3.999 ± 0.0026 g cm^?3 for material decompressed from all parts of the lower mantle. Even if the hot uncompressed density were uniform for all depths in the lower mantle, the cold uncompressed mantle would be inhomogeneous because the decompression given by the Bullen law crosses isotherms; for example, the temperature is different at different depths. To calculate the density distribution correctly, an isothermal EOS must be used along an isotherm, and temperature corrections must be placed in the thermal pressure P_TH.The thermodynamic parameters of the lower mantle are found by iteration. Values of the three uncompressed anharmonic parameters are first arbitrarily selected: α₀ (hot), the coefficient of thermal expansion; γ₀, the Grüneisen parameter; and δ, the second Grüneisen parameter. Using γ₀ and the measured ρ₀ (hot) and K₀ (hot), the values of θ₀ (Debye temperature) and q = dlnγ/dlnρ are found from the measured seismic velocities. Then from (αK_T)₀ and q the thermal pressure P_TH at all high temperatures is found. Correlating P_TH against T to the geotherm for the lower mantle, P_TH is found at all depths Z. The isothermal pressure, along the 0 K isotherm, at every Z is found by subtracting P_TH from the measured P given by the seismic model. Using the isothermal pressure at depth Z, the solution for the cold uncompressed density ρ_0C and the cold uncompressed bulk modulus, K_T0 is found as a trace in the K_T0?_ρ0C plane. A narrow band of solutions is then found for ρ_0C and K_T0 at all depths.The thermal expansion at all T is found from [ρ_0C ? ρ₀ (hot)/ρ_0C. From Suzuki's formula, the best fit to the thermal expansion determines γ₀ and α₀ (hot). When the values of these two parameters do not agree with the original assumptions, the calculation is repeated until they do agree. In this way all the important thermodynamic parameters are found as a self-consistent set subject only to the assumptions behind the equations used.     

Compressibility of brownmillerite (Ca2Fe2O5): effect of vacancies on the elastic properties of perovskites

The evolution with pressure of the unit-cell parameters brownmillerite (Ca₂Fe₂O₅), a stoichiometric defect perovskite structure, has been determined to a maximum pressure of 9.46 GPa, by single-crystal X-ray diffraction measurements at room temperature. Brownmillerite does not exhibit any phase transitions in this pressure range. A fit of a third-order Birch–Murnaghan equation-of-state to the P–V data yields values of K_T0=127.0(5) GPa and K₀′=5.99(13). Analysis of the unit-cell parameter data shows that the structure compresses anisotropically. Compressional moduli for the axes are K_a0=141(1) GPa, K_b0=118(3) GPa and K_c0=122.2(2) GPa, with K_a0′=8.9(3), K_b0′=6.2(6) and K_c0′=4. The stiffest direction (i.e. along a) coincides with the direction of the FeO₄ tetrahedral chains. Comparison of these data with the elasticity systematics of Ca-perovskites shows that the presence of oxygen vacancies in the brownmillerite structure softens the structure by ∼25% and that the ordering of vacancies in the perovskite structure increases the anisotropy of compression.     

High-pressure phase transformations and isothermal compression in CaTiO3 (perovskite)

A polycrystalline CaTiO₃ (perovskite) was investigated under static pressures up to 38 GPa and temperatures up to 1000°C by using a diamond anvil pressure cell, a YAG laser, and the ruby fluorescence pressure calibration system. In situ x-ray diffraction data reveal that at room temperature, the orthorhombic CaTiO₃(I) transforms into a hexagonal CaTiO₃(II) at ∼ 10 GPa with a volume of change of 1.6%. At 1000°C, the orthorhombic CaTiO₃(I) first transforms into a tetragonal CaTiO₃(III) at 8.5 GPa and then transforms further into a hexagonal CaTiO₃(II′) at ∼ 15 GPa with molar volume changes of 0% and 1.6%, respectively. All three high-pressure polymorphs found in this study are nonquenchable.Isothermal compressibility of the orthorhombic CaTiO₃ was derived from measurements under truly hydrostatic environments (i.e., ⩽ 10.4 GPa). By assuming K^′₀ = 5.6 obtained ultrasonically on SrTiO₃ perovskite, the value of the bulk modulus (K₀) was calculated with the Birch-Murnaghan equation to be 210 ± 7 GPa.     

DHI evaluation by combining rock physics simulation and statistical techniques for fluid identification of Cambrian-to-Cretaceous clastic reservoirs in Pakistan

The use of seismic direct hydrocarbon indicators is very common in exploration and reservoir development to minimise exploration risk and to optimise the location of production wells. DHIs can be enhanced using AVO methods to calculate seismic attributes that approximate relative elastic properties. In this study, we analyse the sensitivity to pore fluid changes of a range of elastic properties by combining rock physics studies and statistical techniques and determine which provide the best basis for DHIs. Gassmann fluid substitution is applied to the well log data and various elastic properties are evaluated by measuring the degree of separation that they achieve between gas sands and wet sands. The method has been applied successfully to well log data from proven reservoirs in three different siliciclastic environments of Cambrian, Jurassic, and Cretaceous ages. We have quantified the sensitivity of various elastic properties such as acoustic and extended elastic (EEI) impedances, elastic moduli (K _sat and K _sat–μ), lambda–mu–rho method (λρ and μρ), P-to-S-wave velocity ratio (V _P/V _S), and Poisson’s ratio (σ) at fully gas/water saturation scenarios. The results are strongly dependent on the local geological settings and our modeling demonstrates that for Cambrian and Cretaceous reservoirs, K _sat–μ, EEI, V _P/V _S, and σ are more sensitive to pore fluids (gas/water). For the Jurassic reservoir, the sensitivity of all elastic and seismic properties to pore fluid reduces due to high overburden pressure and the resultant low porosity. Fluid indicators are evaluated using two metrics: a fluid indicator coefficient based on a Gaussian model and an overlap coefficient which makes no assumptions about a distribution model. This study will provide a potential way to identify gas sand zones in future exploration.     

Equation of state of gold and its application to the phase boundaries near 660 km depth in Earth’s mantle

The pressure-volume-temperature equation of state (EOS) of gold is fundamental to high-pressure science because of its widespread use as an internal pressure standard. In particular, the EOS of gold has been used in recent in situ multi-anvil press studies for determination of phase boundaries related to the 660-km seismic discontinuity. These studies show that the boundaries are lower by 2 GPa than expected from the depth of the 660-km discontinuity. Here we report a new P-V-T EOS of gold based on the inversion of quasi-hydrostatic compression and shock wave data using the Mie-Grüneisen relation and the Birch-Murnaghan-Debye equation. The previously poorly constrained pressure derivative of isothermal bulk modulus and the volume dependence of Grüneisen parameter (q=d lnγ/d ln V) are determined by including both phonon and electron effects implicitly: K′_0T=5.0±0.2 and q=1.0±0.1. This combined with other accurately measured parameters enables us to calculate pressure at a given volume and temperature. At 660-km depth conditions, this new EOS yields 1.0±0.2 GPa higher pressure than Anderson et al.’s EOS which has been used in the multi-anvil experiments. However, after the correction, there still exists a 1.5-GPa discrepancy between the post-spinel boundary measured by multi-anvil studies and the 660-km discontinuity. Other potential error sources, such as thermocouple emf dependence on pressure or systematic errors in spectroradiometry, should be investigated. Theoretical and experimental studies to better understand electronic and anharmonic effects in gold at high P-T are also needed.     

Thermal equation of state of majorite with MORB composition

In situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press at SPring-8 on majoritic garnet synthesized from natural mid-ocean ridge basalt (MORB), whose chemical composition is close to the average of oceanic crust, at 19 GPa and 2200 K. Pressure-volume-temperature data were collected using a newly developed high-pressure cell assembly to 21 GPa and 1273 K. Data were fit to the high-temperature Birch-Murnaghan equation of state, with fixed values for the ambient cell volume (V₀ = 1574.14(4) Å³) and the pressure derivative of the isothermal bulk modulus (K^′T = 4). This yielded an isothermal bulk modulus of K_T0 = 173(1) GPa, a temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.022(5) GPa K⁻¹, and a volumetric coefficient of thermal expansivity α = a + bT with values of a = 2.0(3) × 10⁻⁵ K⁻¹ and b = 1.0(5) × 10⁻⁸ K⁻². The derived thermoelastic parameters are very similar to those of pyrope. The density of subducted oceanic crust compared to pyrolitic mantle at the conditions in Earth's transition zone (410-660 km depth) was calculated using these results and previously reported thermoelastic parameters for MORB and pyrolite mineral assembledges. These calculations show that oceanic crust is denser than pyrolitic mantle throughout the mantle transition zone along a normal geotherm, and the density difference is insensitive to temperature at the pressures in lower part of the transition zone.     

High-temperature behaviour of the elastic moduli of LiF and NaF: Comparison with MgO and CaO

The elastic moduli of single-crystal LiF and NaF have been determined by the ultrasonic pulse superposition technique as a function of temperature from T = 298–650° K. These new data are consistent with low-temperature (T < 298° K) data obtained by other ultrasonic pulse techniques and are superior to previous high-temperature data from resonance experiments. The elastic moduli (c) are represented by quadratic functions in T over the experimental temperature range although the curvature is not in the same sense for all modes. For LiF, NaF, MgO and CaO, evaluation of the temperature derivatives of the elastic moduli at constant volume (V) indicates that the elastic moduli are only weakly dependent on T at constant volume. The fluoride—oxide analogue pair LiFMgO both exhibit high-temperature elastic behaviour at approximately the same absolute temperature. Mitskevich's theory and observed K_S-V systematics imply that (?c/?T)_P should be a function of the nearest neighbour distance for rocksalt fluorides and oxides; this result lends further support to a fluorideoxide modelling scheme based on similar ionic radii.     

The effect of ferromagnetism on the equation of state of Fe3C studied by first-principles calculations

First-principles calculations have been used to determine the equation of state of Fe₃C in both its low-pressure magnetically ordered and high-pressure non-magnetically ordered states; at 0 K the ferromagnetic transition was found to occur at about 60 GPa. In the high pressure, non-magnetically ordered regime at 0 K the material may be described by a Birch-Murnaghan third-order equation of state with V₀=8.968(7) ų per atom, K₀=316.62(2) GPa and K′=4.30(2). At atmospheric pressure the ferromagnetic phase transition in Fe₃C occurs at ∼483 K; preliminary measurements of the thermal expansion by powder neutron diffraction show that this transition produces a large effect on thermoelastic properties. The volumetric thermal expansion coefficient in the paramagnetic phase was found to be 4.34×10⁻⁵ K⁻¹ at T∼550 K. By applying a thermal expansion correction to the calculated equation of state at 0 K, predicted values for the density and adiabatic incompressibility of this material at core pressures and temperatures were obtained. These results appear to be sufficiently different from seismological data so as to preclude Fe₃C as the major inner core-forming phase.     

Shear and dilatational wave velocities for unsaturated soils

The primary objective of this study is for presenting some simple-to-use expressions relating the shear and dilatational wave velocities (V_S and V_P) to some physical and constitutive parameters of unsaturated soils. To this purpose, a simplified formulation is developed using the theory of linear poroelasticity in conjunction with some constitutive parameters widely used in geotechnical engineering. The derived expressions are of practical interest in view of the fact that they could be employed for evaluating the involved soil parameters from V_S and V_P measurements by in-situ or laboratory geophysical tests.     

Thermal equation of state of magnesiowüstite (Mg0.6Fe0.4)O

Volume measurements for magnesiowüstite (Mg_0.6Fe_0.4)O, were carried out up to pressures of 10.1 GPa in the temperature range 300–1273 K, using energy-dispersive synchrotron X-ray diffraction. These data allow reliable determination of the temperature dependence of the bulk modulus and good constraint on the thermal expansitivity at ambient pressure which was previously not known for magnesiowüstite. From these data, thermal and elastic parameters were derived from various approaches based on the Birch–Murnaghan equation of state (EOS) and on the relevant thermodynamic relations. The results from three different equations of state are remarkably consistent. With (∂K_T/∂P)_T fixed at 4, we obtained K₀=158(2) GPa, (∂K_T/∂T)_P=−0.029(3) GPa K⁻¹, (∂K_T/∂T)_V=−3.9(±2.3)×10⁻³ GPa K⁻¹, and α_T=3.45(18)×10⁻⁵+1.14(28)×10⁻⁸T. The K₀, (∂K_T/∂T)_P, and (∂K_T/∂T)_V values are in agreement with those of Fei et al. (1992) and are similar to previously determined values for MgO. The zero pressure thermal expansitivity of (Mg_0.6Fe_0.4)O is found to be similar to that for MgO (Suzuki, 1975). These results indicate that, for the compositional range x=0–0.4 in (Mg_1−xFe_x)O, the thermal and elastic properties of magnesiowüstite exhibit a dependence on the iron content that is negligibly small, within uncertainties of the experiments. They are consequently insensitive to the Fe–Mg partitioning between (Mg, Fe)SiO₃ perovskite and magnesiowüstite when applied to compositional models of the lower mantle. With the assumption that (Mg_0.6Fe_0.4)O is a Debye-like solid, a modified equation of heat capacity at constant pressure is proposed and thermodynamic properties of geophysically importance are calculated and tabulated at high temperatures.     

Elasticity of polycrystalline stishovite

Ultrasonic data for the velocities of SiO₂-stishovite have been determined as a function of pressure to 10 kbar at room temperature for polycrystalline specimens hot-pressed at pressures P = 120kbar and temperatures T = 900°C. These cylindrical specimens are 2 mm in diameter and 0.9–1.4 mm long and have a grain size less than 10 μm. Compressional and shear wave velocities were measured both parallel and perpendicular to the axis of pressing and were found to be isotropic at 10 kbar with ν_p = 11.0 ± 0.2km/sec andν_s = 6.9 ± 0.3km/sec; this shear velocit is substantially higher than that of Mizutani et al. (1972) perhaps due to the presence of crack orientations in their specimen which affected ν_s but not ν_p. The Murnaghan P-V trajectories calculated from the ultrasonic data [bulk modulus K_s = 2.5 ± 0.3Mbar and assuming (?K_s/?P)_T = 6 ± 2] are consistent with recent hydrostatic compression data and with the shock wave compression data above 600 kbar. The combined evidence from the data of the ultrasonic and hydrostatic compression techniques suggests that the most probable value of the bulk modulus of stishovite at zero pressure is close to the upper limit of the uncertainty of our ultrasonically determined value, K₀ = 2.7?2.8Mbar. Elasticity data for rutile-type oxides are not compatible with normal K_s-V₀ systematics perhaps due to the neglect of non-central forces in the lattice model. These new stishovite data would make it impossible to satisfy the elasticity-density data of the lower mantle using an oxide mixture with either olivine or pyroxene stoichiometry.     

High-pressure X-ray diffraction studies on β- and γ-Mg2SiO4

Pressure effects on the lattice parameters of β- and γ-Mg₂SiO₄ have been measured at room temperature and at pressures up to 100 kbar using a multi-anvil high-pressure X-ray diffraction apparatus. The volume changes (ΔV/V₀) at 90 kbar are 5.4 · 10^?2 and 4.2 · 10^?2 for β- and γ-Mg₂SiO₄, respectively. Isothermal bulk moduli at zero pressure have been calculated from least-square fits of the data to straight lines. They turn out to be 1.66 ± 0.4 and 2.13 ± 0.1 Mbar for β- and γ-Mg₂SiO₄, respectively. The α → γ transition obeys Wang's linear V_φ?ρ relation but the α → β transition does not.     

On the statistical distribution of seismic velocities in Earth's deep mantle

The existence of uncoupled shear (S) and compression (P) wave velocity variations in Earth's mantle is a characteristic that might only be explained by the presence of significant chemical and/or phase heterogeneity, with important implications for the dynamics and evolution of Earth's interior. While making a one-to-one comparison between tomographic models for P and S velocity (V_P and V_S) variations for a particular geographic region is ill-posed, their global statistical distributions reveal several robust characteristics indicative of the nature of uncoupled V_P and V_S in the deep mantle. We find that all of the V_P and V_S model distributions at a given depth are Gaussian-like throughout the lowermost mantle. However, a distinct low velocity feature is present in V_S distributions below ≈ 2200 km depth that is not present or is relatively weak in V_P models. The presence of anomalously low V_S material cannot be explained as an artifact, nor can the absence of a similarly strong feature in P models be ascribed to under-resolution. We propose that this feature can be partly explained by laterally variable occurrences of post-perovskite (pPv) lenses in the D″ layer, however, the persistence of significantly slow V_S regions at heights up to ≈ 700 km or more above the core–mantle boundary is likely to be incompatible with a pPv origin and might only be explained by the presence of a laterally discontinuous layer of chemically distinct material and/or some other kind of phase heterogeneity. There also exist significant discrepancies between tomographic models with respect to the width of the distributions as well as differences between the modeled peak values. We propose a scheme for comparison between different seismic models in which the widths of the dominant features in their statistical distributions is exploited.     

Applications of liquid state physics to the Earth's core

By use of the modern theory of liquids and some guidance from the hard-sphere model of liquid structure, the following new results have been derived for application to the Earth's outer core. (1) dK/$\text{d}\text{P ? 5 ? 5. 6P}$/K, where K is the incompressibility and P the pressure. This is valid for a high-pressure liquid near its melting point, provided that the pressure is derived primarily from a strongly repulsive pair potential φ. This result is consistent with seismic data, except possibly in the lowermost region of the outer core, and demonstrates the approximate universality of dK/dP proposed by Birch (1939) and Bullen (1949). (2) dlnT_M/dlnρ = (γC_V ? 1)/($\text{C}{}_{\mathrm{<mtext>V</mtext>}}\text{?}\text{3}\text{2}\text{),}\text{where}\text{T}{}_{\mathrm{<mtext>M</mtext>}}$ is the melting point, ρ the density, γ the atomic thermodynamic Grüneisen parameter and C_V the atomic contribution to the specific heat in units of Boltzmann's constant per atom. This reduces to Lindemann's law for C_V = 3 and provides further support for the approximate validity of this law. (3) It follows that the “core paradox” of Higgins and Kennedy can only occur if $\text{γ <}\text{2}\text{3}$. However, it is shown that $\text{γ <}\text{2}\text{3}\text{? ∫}{}^{\mathrm{\infty }}{}_0\text{(?g/?T)}{}_{\mathrm{\rho }}\text{r(}\text{d}\text{/}\text{d}\text{r)(r}{}^2\text{φ)}\text{d}\text{r > 0}$, which cannot be achieved for any strongly repulsive pair potential φ and the corresponding pair distribution function g. It is concluded that $\text{γ >}\text{2}\text{3}$ and that the core paradox is almost certainly impossible for any conceivable core composition. Approximate calculations suggest that γ ～ 1.3–1.5 in the core. Further work on the thermodynamics of the liquid core must await development of a physically realistic pair potential, since existing pair potentials may be unsatisfactory.     

Dynamic characteristics of soil-foundation interaction system detected from forced vibration test and earthquake observation

This paper presents some results on the following subjects obtained from in-situ forced vibration tests and earthquake observations. (1) The characteristics of the radiation damping of soil-foundation interaction systems vs. non-dimensional frequency a₀ (=ωr/V_s) were experimentally estimated by the equivalent damping ratios h_H ( = K′_H/2K_H) and h_R ( = K′_R/2K_R), which were defined by complex stiffnesses 1K_H (= K_H + iK_H) and 1K_R (= K_R + iK′_R) of soil. The results for h_H and h_R of base rock were compared with those of soft soil. (2) A comparative study of experimental and theoretical results was made. The theoretical results were obtained from elastic half-space theory. (3) A semi-empirical equation to estimate the equivalent S-wave velocity for the elastic half-space model is proposed here, considering the effects of layered media. (4) Various comparisons of the results of 1 K_H, 1 K_R, h_H and h_R of forced vibration tests and earthquake observations were made.     

Density profile and modern crustal and uppermost mantle geotherm of the Voronezh crystalline massif

A density profile and a modern temperature distribution in the lithosphere of the Voronezh crystalline massif (VCM) are derived through the use of the V_P(z), V_S(z) seismic velocity models, petrological data, measurements of V_P, V_S, density (ρ) and mean atomic weight ($\text{m}$) for several groups of rocks and minerals of different composition and genesis, as well as from pressure and temperature derivatives for different thermodynamic regimes.