Variation of the magnetic properties in a basaltic dyke with concentric cooling zones

The effects of the variation of magnetic grain size on the magnetic properties of rocks have been studied throughout a reversely magnetized basaltic dyke with concentric cooling zones.Except in a few tachylites in which the magnetic mineral is a Ti-rich titanomagnetite, in the bulk of the dyke the magnetization is carried by almost pure magnetite grains. Although the percentage p of these magnetic oxides varies slightly, the large changes in the various magnetic parameters observed across the dyke are essentially attributable to large variations in the grain size of the magnetic particles.From the outer scoria region, where the magnetic grains are a mixture of single-domain (SD) and superparamagnetic (SP) grains, to the tachylite zone with finely crystallized basaltic glass containing interacting elongated SD particles, one observes an increase of both the ratio of the saturation remanent magnetization and the saturation induced magnetization J_rs/J_is, the bulk coercive force H_c, the median destructive field MDF, the intensity of the remanent magnetization J_r, and the Koenigsberger ratio Q. In the tachylites these parameters reach unusually high values, for subaerial basalts:

$\text{〈}\text{J}{}_{\mathrm{rs}}\text{J}{}_{\mathrm{is}}\text{〉 = 0.3, 〈H}{}_{\mathrm{c}}\text{〉 = 460 O}{}_{\mathrm{e}}\text{, 〈MDF〉 = 620 O}{}_{\mathrm{e}}\text{r.m.s., 〈J}{}_{\mathrm{r}}\text{〉 = 2.7 · 10}{}^{\mathrm{?2}}\text{e.m.u. cm}{}^{\mathrm{?3}}\text{〈Q〉 = 24}$

These parameters decrease in the basalt toward the centre of the dyke where pseudo-single-domain (pseudo-SD) particles coexist together with multidomain (MD) grains. The susceptibility remains approximately constant from the inner basalt to the tachylite, but increases in the scoria up to values 10 times higher owing to the presence of SP particles. The magnetic viscosity increases also drastically toward the margin of the dyke due to an increase of the fraction of the SD particles just above the superparamagnetic threshold.     

Effect of loss-cone distribution on the generation of whistler waves in the magnetosphere

An anisotropic kappa velocity distribution with loss-cones is used to investigate whistler wave instability occurring in the magnetosphere. The elements of the dielectric tensor and dispersion relation using modified plasma dispersion function $\text{Z}{}_{\mathrm{\kappa }}{}^1\text{(ξ)}$ with loss-cone angle have been obtained for the linear waves propagating exactly parallel to a uniform local magnetic field in a homogeneous and hot plasma. The modified plasma dispersion function and integrals have been expressed in power-series form for argument of ξ≫1. Temporal/spatial growth rates for whistler wave in the magnetosphere have been evaluated by the method of numerical techniques. The results of such a kappa loss-cone distribution function on the generation of whistler waves are compared with those obtained by Maxwellian loss-cone distribution. Calculations show that either a loss-cone or a thermal anisotropy in the hot plasma component of the magnetosphere can lead to the generation of incoherent emission of low-frequency whistler waves. This methodology could be easily extended to the study of low frequency emissions from planetary magnetospheres under suitable choice of models of density and magnetic field and other plasma parameters.     

Elasticity of ilmenites

Ultrasonic data for the velocities of the ilmenites MgTiO₃ and CoTiO₃ have been determined as a function of pressure to 7.5 kbar at room temperature for polycrystalline specimens hot-pressed in a piston-cylinder apparatus at pressures up to 30 kbar. Titanate and germanate ilmenites define divergent isostructural trends on a Birch diagram of bulk sound velocity (υ_φ) vs. density (ρ). On a υ_φ vs. mean atomic weight ($\text{M}$) diagram, however, all of the ilmenite consistent with a single $\text{υ}{}_{\mathrm{\phi }}\text{M}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}\text{=}\text{constant trend}$. Elasticity systematics for isostructural sequences are used to e the bulk modulus (2.09 Mbar) and bulk sound velocity (7.4 km/sec) of MgSiO₃-elmenite.     

General relationships among sound speeds: II. Theory and discussion

The dependence of bulk sound speed V_φ upon mean atomic weight $\text{m}$ and density ρ can be expressed in a single equation:

$\text{V}{}_{\mathrm{\phi }}\text{=Bρ}{}^{\mathrm{\lambda }}\text{(}\text{m}{}_0\text{m}{}^{\mathrm{<mtext>[</mtext><mtext>1</mtext><mtext>2</mtext><mtext>+\lambda (1?c)]</mtext>}}\text{(}\text{km/sec}\text{)}$

Here B is an empirically determined “universal” parameter equal to 1.42, $\text{m}{}_0\text{= 20.2}$, a reference mean atomic weight for which well-determined elastic properties exist, and λ = 1.25 is a semi empirical parameter equal to $\text{γ ?}\text{1}\text{3}$ where γ is a Grüneisen parameter. The constant c = (? ln V_M/? ln $\text{m}\text{)}{}_{\mathrm{X}}$, where V_M is molar volume, is in general different for different crystal structure series and different cation substitutions. However, it is possible to use c_Fe = 0.14 for Fe²⁺Mg²⁺ and GeSi substitutions and c_Ca ? 1.3 for CaMg substitutional series. With these values it is pos to deduce from the above equation Birch's law, its modifications introduced by Simmons to account for Ca-bearing minerals, variations in the seismic equation of state observed by D.L. Anderson, and the apparent proportionality of bulk modulus K to V_M^?4.     

On the mechanism of the olivine-spinel phase transformation

Recent experiments on the olivine-spinel phase transformation in Ni₂SiO₄ have provided the following fundamental results: (1) the interphase structural orientation relationships with $\text{(100)}{}_{\mathrm{<mtext>ol</mtext>}}\text{|(111)}{}_{\mathrm{<mtext>sp</mtext>}}\text{and}\text{[001]}{}_{\mathrm{<mtext>ol</mtext>}}\text{|[}\text{1}\text{10]}{}_{\mathrm{<mtext>sp</mtext>}}$ between the olivine and spinel phases upon transformation; and (2) the growth of 300 μm spinel single crystal from single crystals of olivine. In the light of these observations, transformation mechanisms proposed have been compared. We conclude that neither a nucleation and growth model nor a martensite model is verified by the experimental results. The former model lacks consistency in terms of the interphase relationships, and the latter cannot offer a dislocation-multiplication mechanism for the growth process.     

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.     

Observational study of terrestrial eigenvibrations

A comprehensive analysis has been made of analog and digital recordings of eigenvibration ground motion obtained following four great earthquakes; August 1976 (Philippines), August 1977 (Indonesia), September 1979 (West Irian), and December 1979 (Colombia). The time series (ranging in length from ～28 to ～140 h) are assumed to be linear combinations of damped harmonics in the presence of noise. Tables are calculated from values of the four parameters: Θ, used in describing eigenvibrations, period of oscillation, amplitude, damping factor Q, and phase together with their statistical uncertainties (53 spheroidal modes, ${}_0\text{S}{}_4\text{to}{}_0\text{S}{}_{48}$, and 13 torsional modes, ${}_0\text{T}{}_8\text{to}{}_0\text{T}{}_{45}$). The estimation procedures are by the methods of complex demodulation and non-linear regression that specifically incorporate into the basic model the decaying aspect of the oscillations. These methods, extended to simultaneous estimations of groups of modes, help to eliminate measurement error and measurement bias from estimations of Θ. The result is that overtone modes very near in frequency to fundamental modes can, under certain conditions, be resolved through a non-linear regression technique, although parameter uncertainties are underestimated in general.Of the time series analyzed, 17 were from a northern California regional network of ultra-long period seismographs at Berkeley (three components), Jamestown (vertical component), and Whiskeytown (vertical component) following the four listed earthquakes. The other 7 time series were recorded digitally by the worldwide IDA network following the 1977 Indonesian earthquake. Weighted regional and worldwide averages were made for period and Q of each eigenvibration mode.From the theoretical viewpoint, comparisons of measured period, Q, amplitude, and phase for all modes analyzed led to five conclusions. First, there are no detectable systematic shifts in period, Q, or phase of eigenvibrations within a region whose dimensions are less than a wavelength. Second, though not conclusive, there may be slight systematic shifts in period (<0.65 s) and relative amplitudes within the California regional network due to different source positions and mechanisms. Differences in Q values are not statistically significant. Third, even though differences in period obtained worldwide were as great as 1.33 s (≈0.33%), differences between Q values (as great as 20%) for the same mode were not significant. The conclusion is that the damping characteristics of singlet eigenfunctions are not observed to be significantly different. Fourth, the assumption that a multiplet ${}_{\mathrm{n}}\text{S}{}_{\mathrm{l}}$ behaves as a single oscillation is valid from at least ${}_0\text{S}{}_7$ through ${}_0\text{S}{}_{30}$. Fifth, no systematic pattern emerged for the shift of eigenperiod as a function of order / or posit on the Earth.     

Spatial power spectra of the crustal geomagnetic field and core geomagnetic field

Lowes (1966, 1974) has introduced the function R_n defined by $\text{R}{}^{\mathrm{n}}\text{=(n + 1)}\text{∑}\text{m=0}\text{∞}\text{[(g}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{+ (h}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{]}\text{where}\text{g}{}_{\mathrm{n}}{}^{\mathrm{m}}\text{and}\text{h}{}_{\mathrm{n}}{}^{\mathrm{m}}$ are the coefficients of a spherical harmonic expansion of the scalar potential of the geomagnetic field at the Earth's surface. The mean squared value of the magnetic field B = ??V on a sphere of radius r > α is given by $\text{〈}\text{B}\text{·〉 =}\text{∑}\text{n=1}\text{∞}\text{R}{}^{\mathrm{n}}\text{(}\text{a/}\text{r}\text{)}{}^{\mathrm{2n=4}}\text{where}\text{a}$ is the Earth's radius. We refer to R_n as the spherical harmonic spatial power spectrum of the geomagnetic field.In this paper it is shown that Rⁿ = R^M_n = R^C_n where the components R_n^M due to the main (or core) field and R_n^C due to the crustal field are given approximately by $\text{R}{}^{\mathrm{M}}{}_{\mathrm{n}}\text{= [}\text{(n =1)/}\text{(n + 2)}\text{](1.142 × 10}{}^9\text{)(0.288}{}^{\mathrm{n}}\text{Λ}{}^2\text{R}{}^{\mathrm{C}}{}_{\mathrm{n}}\text{= [}\text{(n =1){[1 — exp(-n/290)]/}\text{(n/290)} 0.52 Λ}{}^2\text{where}\text{Iγ = 1 nT}$. The two components are approximately equal for n = 15.Lowes has given equations for the core and crustal field spectra. His equation for the crustal field spectrum is significantly different from the one given here. The equation given in this paper is in better agreement with data obtained on the POGO spacecraft and with data for the crustal field given by Alldredge et al. (1963).The equations for the main and crustal geomagnetic field spectra are consistent with data for the core field given by Peddie and Fabiano (1976) and data for the crustal field given by Alldredge et al. The equations are based on a statistical model that makes use of the principle of equipartition of energy and predicts the shape of both the crustal and core spectra. The model also predicts the core radius accurately. The numerical values given by the equations are not strongly dependent on the model.Equations relating average great circle power spectra of the geomagnetic field components to R_n are derived. The three field components are in the radial direction, along the great circle track, and perpendicular to the first two. These equations can, in principle, be inverted to compute the R_n for celestial bodies from average great circle power spectra of the magnetic field components.     

Evidences for a low-velocity core-mantle transition zone

Using the spectral ratios $\text{P}\text{PcP}\text{,}\text{ScS}{}_{\mathrm{n+1}}\text{ScS}{}_{\mathrm{n}}\text{,}\text{sScS}{}_{\mathrm{n+1}}\text{sScS}{}_{\mathrm{n}}\text{and}\text{SKS}\text{ScS}$, models for the core-mantle boundary are found. The models have close similarity with each other, implying an irregular surface with lateral variation in the core-mantle properties. The models are characterized by two to four low-velocity, high-density layers imbedded between the mantle and the core half space. The velocities of the imbedded layers decrease towards the core boundary with a lower bound of 9.3 km/sec for the compressional wave and 3.5 km/sec for the shear wave. The models fitted to the empirical data support the hypothesis of a finite rigid outer core with a higher bound for the shear velocity of 1.4 km/sec. Based on this finite rigidity in the outer core and a layered core-mantle transition zone, the value of Q for the whole mantle is 2,000. For the outer core Q ranges from 100–1,000 , which may indicate that it is chemically zoned.     

Synthesis of a new high-pressure ohase of tin dioxide and some geophysical implications

Tin dioxide (SnO₂) in the rutile structure as starting material has been found to transform to the orthorhombic α-PbO₂ structure (S.G. Pbcn) at about 155 kbar and 1000–1400°C when compressed in a diamond-anvil cell and heated by irradiation with a YAG laser. The lattice parameters at room temperature and 1 bar are $\text{a}{}_{\mathrm{<mtext>o</mtext>}}\text{= 4.719 ± 0.002, b}{}_{\mathrm{<mtext>o</mtext>}}\text{= 5.714 ± 0.002,}\text{and}\text{c}{}_{\mathrm{<mtext>o</mtext>}}\text{= 5.228 ± 0.002}\text{A}\text{?}\text{with}\text{Z = 4}$ for the orthorhombic form of SnO₂, which is 1.5% more dense than the rutile form. Crystal-chemical arguments suggest that stishovite (SiO₂) may also transform to the α-PbO₂ structure at elevated pressure and temperature with an increase in zero-pressure density of about 2–3%. Mineral assemblages containing the orthorhombic SiO₂ are unstable relative to those containing the perovskite MgSiO₃ under lower-mantle conditions.     

Physical constraints for the analysis of the geomagnetic secular variation

Slow changes in the magnetic field are believed to originate in the core of the Earth. Interpretation of these changes requires knowledge both of the vertical component of the field and of its rate of change at the core-mantle boundary (CMB). While various spherical harmonic models show some agreement for the field at the CMB, those for secular variation (SV) do not. SV models depend heavily on annual means at relatively few and poorly distributed magnetic observatories. In this paper, the SV at the CMB is modelled by fitting 15-year differences in the annual means of the X, Y and Z components (from 1959 to 1974). The model is made unique by imposing the constraint that $\text{?}{}_{\mathrm{CMB}}\text{B}\text{?}{}_{\mathrm{r}}{}^2\text{d}\text{S}$ be a minimum, using the method of Shure et al. (1982). If SV is attributed to motions of core fluid, then this model will yield, in some sense, the slowest core motions. The null space is determined by the distribution of observations, and therefore, to be consistent, only those observatories have been retained which recorded almost continuously throughout the interval 1959–1974.The method allows misfit between the model and the observations. The best value for the misfit can be derived from estimates of errors in the data, or alternatively, because larger misfit leads to smoother models (i.e., smaller $\text{?}\text{B}\text{?}{}_{\mathrm{r}}{}^2\text{d}\text{S}$), the best value can be estimated subjectively from the final appearance of the model. Both procedures have their counterparts in the conventional spherical harmonic expansion approach, when smoothing is achieved by lowering the truncation level. The new proposal made in this paper is to use objective criteria for determining the misfit, based on the assumption that diffusion is negligible, in which event all integrals $\text{∫}\text{B}\text{?}{}_{\mathrm{r}}{}^2\text{d}\text{S}$ will vanish when S_i is a region on the CMB bounded by a contour of zero vertical component of field. For the 1965 definitive model which is adopted here, and for most other contemporary models, there are six such areas, giving five independent integrals (the integrals over the six regions must sum to zero if ? · B = 0). Tabulating these integrals for various choices of the misfit gives minimum values near 2 nT y^?1. It is impossible to achieve this good a fit to the data using a reasonable model derived by truncating the spherical harmonic expansion. The value 2 nT y^?1 corresponds to errors of ～ 20 nT in individual annual means, which is rather larger than expected from the scatter in the data.     

Rate of seismicity of the dead sea region over the past 4000 years

The results of two millennia of earthquake documentation, a few decades of macroseismic and instrumental routine seismological observations and five months of microearthquake monitoring, are used to estimate the rate of seismic activity of the Dead Sea fault. It is found that these vastly diverse data which combine long- and short-term tectonic processes, are in good accord with the formula:

$\text{log}{}_{10}\text{N=2.54 ? 0.86M}{}_{\mathrm{<mtext>L</mtext>}}$

where N is the annual number of events of local magnitude M_L or greater. If this equation is extrapolated to ca. 2000 B.C., it yields a Richter magnitude M_s = 7 for the event of Sodom and Gomorrah which is believed to be associated with the strongest earthquake in the region during historical times.Comparing our findings with the results of other investigators in Turkey, Greece, Aegean Sea and Iran, we note that the b values along the Syrian-African rift zone (0.78–0.86) are smaller than those in Greece and its surrounding seas (0.94–1.16).     

Melting temperature distribution and fractionation in the lower mantle

The melting curve of perovskite MgSiO₃ and the liquidus and solidus curves of the lower mantle were estimated from thermodynamic data and the results of experiments on phase changes and melting in silicates.The initial slope of the melting curve of perovskite MgSiO₃ was obtained as $\text{d}\text{T}{}_{\mathrm{<mtext>m</mtext>}}\text{/}\text{d}\text{P?77}\text{K}\text{GPa}{}^{\mathrm{?1}}$ at 23 GPa. The melting curve of perovskite was expressed by the Kraut-Kennedy equation as $\text{T}{}_{\mathrm{<mtext>m</mtext>}}\text{(}\text{K}\text{)=917(}\text{1+29.6ΔV}\text{V}{}_0\text{)}$, where T_m?2900 K and P?23 GPa; and by the Simon equation, $\text{P(}\text{GPa}\text{)?23=21.2[}\text{(T}{}_{\mathrm{<mtext>m</mtext>}}\text{(}\text{K}\text{)}\text{2900)}{}^{1.75}\text{?1}\text{]}$.The liquidus curve of the lower mantle was estimated as $\text{T}{}_{\mathrm{<mtext>liq</mtext>}}\text{? 0.9 T}{}_{\mathrm{<mtext>m</mtext>}}$ (perovskite) and this gives the liquidus temperature T_liq=7000 ±500 K at the mantle-core boundary. The solidus curve of the lower mantle was also estimated by extrapolating the solidus curve of dry peridotite using the slope of the solidus curve of magnesiowüstite at high pressures. The solidus temperature is ～ 5000 K at the base of the lower mantle. If the temperature distribution of the mantle was 1.5 times higher than that given by the present geotherm in the early stage of the Earth's history, partial melting would have proceeded into the deep interior of the lower mantle.Estimation of the density of melts in the MgOFeOSiO₂ system for lower mantle conditions indicates that the initial melt formed by partial fusion of the lower mantle would be denser than the residual solid because of high concentration of iron into the melt. Thus, the melt generated in the lower mantle would tend to move downward toward the mantle-core boundary. This downward transportation of the melt in the lower mantle might have affected the chemistry of the lower mantle, such as in the D″ layer, and the distribution of the radioactive elements between mantle and core.     

Theoretical estimates of the westward drift

Virtually all dynamo models may be expected to give rise to a permanent differential rotation between mantle and core. Weak conductivity in the mantle permits small leakage currents which couple to the radial component of the magnetic field, producing a Lorentz torque. Mechanical equilibrium is achieved when a zero net torque is established at a critical rotation rate. An estimate of the drift is determined easily given the magnetic field structure predicted by any dynamo model. The result for the drift rate at the core-mantle interface along the equator is given by the product of three factors

$\text{U}{}_{\mathrm{\phi 1}}\text{=U}{}_{\mathrm{\phi }}\text{R}\text{λ}{}_{\mathrm{*}}\text{L}{}_{\mathrm{*}}$

The first of these is a geometrical factor which depends only on the structural character of the field. For a variety of model fields, this factor ranges from 16 to 35. The second factor is the ratio of r.m.s. toroidal to poloidal field. This ratio is an (implicitly) adjustable parameter of both α² and α-ω dynamos, and is a measure of the relative efficiency of the generation process for each component. The third (dimensional) term is the ratio of core magnetic diffusivity to core radius, 10^?4 cm s^?1.The result is essentially independent of the value of mantle diffusivity and its effective depth. The sign of the result may be positive or negative. For α² dynamos a westward drift is produced by choosing α > 0 in the Northern Hemisphere, which constitutes a dynamical assertion about the dynamo process. For an r.m.s. toroidal field of the order of 15 Gs, based on fairly general considerations, a drift rate comparable to observation is expected.     

Orthogonality of spherical harmonic coefficients

An essentially arbitrary function V(θ, λ) defined on the surface of a sphere can be expressed in terms of spherical harmonics $\text{V(θ, Λ) = a}\text{∑}\text{n=1}\text{∞}\text{∑}\text{m=0}\text{n}\text{p}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{(cos θ) (g}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{cos mΛ + h}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{sin mΛ)}$ where the P_n^m are the seminormalized associated Legendre polynomials used in geomagnetism, normalized so that $\text{〈[P}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{(cos θ) cos mΛ]〉}{}^2\text{=1/}\text{(2n+1)}$ The angular brackets denote an average over the sphere. The class of functions V(θ, λ) under consideration is that normally of interest in physics and engineering. If we consider an ensemble of all possible orientations of our coordinate system relative to the sphere, then the coefficients g_n^m and h_n^m will be functions of the particular coordinate system orientation, but $\text{〈:(g}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{〉) = 〈(h}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{=}\text{S}{}_{\mathrm{n}}\text{/}\text{(2n=1)}$ where $\text{S}{}_{\mathrm{n}}\text{=}\text{∑}\text{m=0}\text{n}\text{[(g}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{+ (h}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{]}$ for any orientation of the coordinate system (S_n is invariant under rotation of the coordinate system). The averages are over all orientations of the system relative to the sphere. It is also shown that 〈g^m_ng^l_p〉 = 〈h^m_nh^l_p〉 = 0 for l ≠ m or p ≠ n and 〈g^m_nh^l_p〉 = 0 fro all n, m, p, l.     

High-pressure phase transformations and compressions of ilmenite and rutile,I. Experimental results

Natural ilmenite (Fe,Mg)TiO₃ has been found to transform to the perovskite structure and then to disproportionate into its component oxides, (Fe,Mg)O plus a cubic phase of TiO₂, at loading pressures of 140 and 250 kbar respectively, and at temperatures of 1,400 to 1,800°C. Samples were compressed in a diamond-anvil press and heated by irradiation with a YAG laser. The lattice parameters of the perovskite phase of (Fe,Mg)TiO₃ at room temperature and 1 bar are $\text{a}{}_0\text{= 4.471 ± 0.004, b}{}_0\text{= 5.753 ± 0.005,}\text{and}\text{c}{}_0\text{= 7.429 ± 0.006}\text{A}\text{?}$ with 4 molecules per cell. The zero-pressure volume change is 8.0% for the ilmenite-perovskite transition, 13.3% for the perovskite-mixed-oxides transition, and 20.2% for the ilmenite-mixed-oxides transition. The cubic phase of TiO₂ can be indexed on the basis of space group Fm3m with $\text{Z = 4}\text{and}\text{a}{}_0\text{= 4.455 ± 0.008}\text{A}\text{?}$ at room temperature and 1 bar, which corresponds to a decrease in zero-pressure volume of 29.2% for the rutile-cubic-phase transition. An isentropic bulk modulus at zero pressure of 5.75 ± 0.30 Mbar and a pressure derivative greater than 8 were calculated for the high-pressure cubic phase. The calculated bulk modulus for the mixture of (Fe,Mg)O and cubic TiO₂ is 2.48 ± 0.25 Mbar. All the phase transformations, the calculated lattice parameters, and the bulk moduli observed in this study are in good agreement with published shock-Hugoniot data for ilmenite and rutile.     

Interionic repulsive force and compressibility of ions

The compressibility of an individual ion is examined, in comparison with a known set of data for the alkali halides. A simple extrapolation of ionic radius to high pressure is not acceptable, because the pressure derivative of ionic radius changes for different salts. According to the classical concept of an elastic ion, the repulsive potential energy between the ions i and j is specified by the nature of each ion as:

$\text{(ρ}{}_{\mathrm{i}}\text{+ ρ}{}_{\mathrm{j}}\text{) exp[}\text{(ρ}{}_{\mathrm{i}}\text{+ ρ}{}_{\mathrm{j}}\text{? r)}\text{(ρ}{}_{\mathrm{i}}\text{+ ρ}{}_{\mathrm{j}}\text{)}\text{]}$

as a function of the interionic distance r. In this expression, q_i and ρ_i are the ionic radius and ionic compressibility, respectively, in a suitably modified meaning. Such a form of the repulsive potential fits well to the data of lattice constants and bulk moduli. The parameters q_i and ρ_i are evaluated for alkali and halogen ions, and an anion turns out to be much more compressible than a cation. The present treatment may be usefully applied to the minerals in the Earth's mantle, which contain only a few major ions.     

High-pressure recovery of olivine: implications for creep mechanisms and creep activation volume

High-temperature and high-pressure recovery experiments were made on experimentally deformed olivines at temperatures of 1613–1788 K and pressures of 0.1 MPa to 2.0 GPa. In the high-pressure experiments, a piston cylinder apparatus was used with BN and NaCl powder as the pressure medium, and the hydrostatic condition of the pressure was checked by test runs with low dislocation density samples. No dislocation multiplication was observed. The kinetics of the dislocation annihilation process were examined by different initial dislocation density runs and shown to be of second order, i.e.

$\text{d}\text{ρ}\text{d}\text{t}\text{= ?p}{}^2\text{K}{}_0\text{exp}\text{[?(}\text{E}{}^1\text{+PV}{}^1\text{RT}\text{]}$

where ρ is the dislocation density, k₀ is a constant, $\text{E}{}^1\text{and}\text{V}{}^1$ are the activation energy and volume respectively, and P, R and T are pressure, gas constant and temperature, respectively. Activation energy and volume were estimated from the temperature and pressure dependence of the dislocation annihilation rate as $\text{E}{}^1\text{=389±59}\text{kJ mol}{}^{\mathrm{?1}}$ and $\text{V}{}^1\text{=14±2}\text{cm}{}^3\text{mol}{}^{\mathrm{?1}}$, respectively.The diffusion constants relevant to the dislocation annihilation process were estimated from a theoretical relation k=αD where $\text{k=k}{}_0\text{exp}\text{[?(E}{}^1\text{+ PV}{}^1\text{)/RT], D}$ is the diffusion constant and α is a non-dimensional constant of ca. 300. The results agree well with the self-diffusion constant of oxygen in olivine. This suggests that the dislocation annihilation is rate-controlled by the (oxygen) diffusion-controlled dislocation climb.The mechanisms of creep in olivine and dry dunite are examined by using the experimental data of static recovery. It is suggested that the creep of dry dunite is rate-controlled by recovery at cell walls or at grain boundaries which is rate-controlled by oxygen diffusion. Creep activation volume is estimated to be 16±3 cm³ mol^?1.