Shear modulus and Q of forsterite and dunite near partial melting from forced-oscillation experiments

An apparatus designed to determine the complex shear modulus of rock samples by forced torsion oscillations at high temperature and in the seismic frequency band 0.003–30 Hz is briefly described. Measurements were performed on natural dunite from Åheim (Norway) up to 1400°C and on polycrystalline forsterite up to 1500°C at 1 atm pressure. The two materials were chosen to study, by comparison, the effect of melt on the elasticity and anelasticity of mantle rocks.Between 1000 and 1200°C the absolute values of the shear modulus G are almost equal for both materials. Above 1200°C G for natural dunite decreases progressively with temperature and at 1400°C and 1 Hz reaches $\text{～}\text{1}\text{3}$ of its value at 1100°C. In contrast, G of pure forsterite depends little on temperature. For petrological reasons, supported by simultaneous measurements of the electric resistivity, there is strong evidence that the decrease of G in dunite above 1200°C is due to melt from the lower melting components of the dunite. Based on different models estimates of the melt fraction are made.At high temperature, in both materials Q^?1 is characterized by a monotonic decrease with frequency according to ω^?α, with α ≈ 0.25. An apparent activation energy of 38±5 kcal mol^?1 for forsterite and 48±8 kcal mol^?1 for dunite was found with no significant change in the regime of partial melting. From this it is concluded that Q^?1, even at partial melting, is dominated by solid state high temperature background absorption. There is no indication from these experiments for a constant-Q-band at low seismic frequencies or an increase of Q proportional to frequency as suggested by some seismologists. The present results are in good qualitative agreement with those for Young's modulus obtained previously by strain retardation experiments.     

Creep response of the lunar crust in mare regions from an analysis of crater deformation

The settling trends of 318 lunar mare craters are compared with predictions of numerical finite-element models in order to determine the creep response of the upper lunar mare crust. No settling is evident in craters smaller than 5 km in diameter. Settling rates of larger craters increase as function of crater size in a manner suggesting a non-linear lunar creep response corresponding to the power law $\text{ε}\text{?}\text{= 8.3 · 10}{}^{\mathrm{?34}}\text{σ}{}^2$ where έ is the strain rate and σ is the differential stress. However, the observed nonlinearity is probably an apparent nonlinearity resulting from the temperature induced viscosity decrease with depth due to a lunar crustal temperature gradient of 3° C/km and a creep activation energy of 20 kcal/mole. It is concluded that creep in the lunar medium is essentially Newtonian, and that the effective viscosity of the upper lunar mare crust is (1.6 ± 0.3) · 10²⁵ poise.     

The cooling Earth: A reappraisal

By treating the lithosphere as a diffusive boundary layer to mantle convection, the convective speed or mantle creep rate, $\text{?}\text{?}$, can be related to the mantle-derived heat flux, $\text{Q}\text{?}$. If cell size is independent of $\text{Q}\text{?}{}^2$ then $\text{?}\text{?}\text{∝}\text{Q}\text{?}$. (If cell size varies with $\text{Q}\text{?}$, then a different power law prevails, but the essential conclusions are unaffected.) Then the factthat for constant thermodynamic efficiency of mantle convection, the mechanical power dissipation is proportionalto $\text{Q}\text{?}$, gives convective stress $\text{σ ∝}\text{Q}\text{?}{}^{\mathrm{?1}}$, i.e. the stress increases as the convection slows. This means an increasing viscosityor stiffness of the mantle which can be identified with a cooling rate in terms of a temperature-dependent creep law. If we suppose that the mantle was at or close to its melting point within 1 or 2 × 10⁸ years of accretionof the Earth, the whole scale of cooling is fixed. The present rate of cooling is estimated to be about 4.6 × 10^?8 deg y^?1 for the average mantle temperature, assumed to be 2500 K, but this very slow cooling rate represents a loss ofresidual mantle heat of 7 × 10¹² W, about 30% of the total mantle-derived heat flux. This conclusion requires theEarth to be deficient in radioactive heat, relative to carbonaceous chondrites. A consideration of mantle outgassing and atmospheric argon leads to the conclusion that the deficiency is due to depletion of potassium, and that the K/U ratio of the mantle is only about 2500, much less than either the crustal or carbonaceous chondritic values. Thetotal terrestrial potassium is estimated to be about 6 × 10²⁰ kg. Acceptance of the cooling of the Earth removes the necessity for potassium in the core.     

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.     

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.     

The effect of interstitial gaseous pressure on the thermal conductivity of a simulated Apollo 12 lunar soil sample

The thermal conductivity of a simulated Apollo 12 lunar soil sample was measured with a needle probe under vacuum. The result showed that the sample, with bulk densities of 1.70–1.85 g cm^?3 held in a vertical cylinder (2.54 cm in diameter and 6.99 cm long) has a thermal conductivity ranging from 8.8 to 10.9 mW m^?1 K^?1. This is comparable to the lunar regolith's thermal conductivity as determined in situ. Besides the dense packing of the soil particles, an enhanced intergranular thermal contact, due to the self-compression of the sample, is necessary to raise the sample's thermal conductivity from the level of loose soil (< 5 mW m^?1 K^?1) to that of the lunar regolith deeper than 35 cm (～ 10 mW m^?1 K^?1). A model of the lunar regolith, a thin layer of loose soil resting on a compacted self-compressed substratum, is consistent with the lunar regolith's surface structure as deduced from an observation of the lunar surface's brightness temperature. Martian regolith surface structure is similar, except that its surface layer may be missing in places because of aeolian activity. Measurements of thermal conductivity under simulated martian surface conditions showed that the thermal properties of loose and compacted soils agreed with the two peak values of the martian surface's thermal inertia as observed from “Viking” orbiters, suggesting that drifted loose soil and exposed compacted soil are responsible for the bimodal distribution of the martian surface's thermal inertia near zero elevation. For compacted soil exposed to the martian surface to have the same thermal conductivity as that buried under the surface layer, a cohesion of the soil particles must be assumed.     

The electric field produced in the mantle by the dynamo in the core

A rigorous singular perturbation theory is developed to estimate the electric field E produced in the mantle M by the core dynamo when the electrical conductivity σ in M depends only on radius r, and when |r?_rln σ| ? 1 in most of M. It is assumed that σ has only one local minimum in M, either (a) at the Earth's surface ?V, or (b) at a radius b inside the mantle, or (c) at the core-mantle boundary ?K. In all three cases, the region where σ is no more than e times its minimum value constitutes a thin critical layer; in case (a), the radial electric field E_r ≈ 0 there, while in cases (b) and (c), E_r is very large there. Outside the critical layer, E_r ≈ 0 in all three cases. In no case is the tangential electric field E_S small, nearly toroidal, or nearly calculable from the magnetic vector potential A as ??_tA_S. The defects in Muth's (1979) argument which led him to contrary conclusions are identified. Benton (1979) cited Muth's work to argue that the core-fluid velocity u just below ?K can be estimated from measurements on ?V of the magnetic field B and its time derivative $\text{?}{}_{\mathrm{t}}\text{B}$. A simple model for westward drift is discussed which shows that Benton's conclusion is also wrong.In case (a), it is shown that knowledge of σ in M is unnecessary for estimating E_S on ?K with a relative error |r?_r 1nσ|^?1from measurements of E_S on ?V and knowledge of ?_tB in M (calculable from ?_tB on ?V if σ is small). Then, in case (a), u just below ?K can be estimated as $\text{?}\text{r}\text{×}\text{E}{}_{\mathrm{S}}\text{/B}{}_{\mathrm{r}}$. The method is impractical unless the contribution to E_S on ?V from ocean currents can be removed.The perturbation theory appropriate when σ in M is small is considered briefly; smallness of σ and of |r?_r ln σ|^?1 a independent questions. It is found that as σ → 0, B approaches the vacuum field in M but E does not; the explanation lies in the hydromagnetic approximation, which is certainly valid in M but fails as σ → 0. It is also found that the singular perturbation theory for |r?_r ln σ|^?1 is a useful tool in the perturbation calculations for σ when both σ and |r?_r ln σ|^?1 are small.     

Heat transport by variable viscosity convection and implications for the Earth's thermal evolution

The case is presented that the efficiency of variable viscosity convection in the Earth's mantle to remove heat may depend only very weakly on the internal viscosity or temperature. An extensive numerical study of the heat transport by 2-D steady state convection with free boundaries and temperature dependent viscosity was carried out. The range of Rayleigh numbers (Ra) is 10⁴?10⁷ and the viscosity contrast goes up to 250000. Although an absolute or relative maximum of the Nusselt number (Nu) is obtained at long wavelength in a certain parameter range, at sufficiently high Rayleigh number optimal heat transport is achieved by an aspect ratio close to or below one. The results for convection in a square box are presented in several ways. With the viscosity ratio fixed and the Rayleigh number defined with the viscosity at the mean of top and bottom temperature the increase of Nu with Ra is characterized by a logarithmic gradient β = ?ln(Nu)/? ln(Ra) in the range of 0.23–0.36, similar to constant viscosity convection. More appropriate for a cooling planetary body is a parameterization where the Rayleigh number is defined with the viscosity at the actual average temperature and the surface viscosity is fixed rather than the viscosity ratio. Now the logarithmic gradient β falls below 0.10 when the viscosity ratio exceeds 250, and the velocity of the surface layer becomes almost independent of Ra. In an end-member model for the Earth's thermal evolution it is assumed that the Nusselt number becomes virtually constant at high Rayleigh number. In the context of whole mantle convection this would imply that the present thermal state is still affected by the initial temperature, that only 25–50% of the present-day heat loss is balanced by radiogenic heat production, and the plate velocities were about the same during most of the Earth's history.     

Geothermal effects in the formation of electrically conducting zones and temperature distribution in the earth

Pollack and Chapman have shown that the surface heat flow in continental regions is dependent not only on the earth's crust below the observation site, but also on the upper mantle there. Therefore heat flow can be used to investigate the role of the thermal conditions in the creation of the electrically conductive zones in both the crust and mantle.Empirical exponential formulas describe the depth to the conductivity increase in the crust corresponding to granitization, the depth to the conductive zone at the top of the asthenosphere (SLVZ), as a function of heat flow. Comparing the latter with temperature estimations in the asthenosphere it is concluded that partial melting of the upper mantle occurs only where $\text{q ? 42}\text{m W m}{}^{\mathrm{?2}}\text{? 1}\text{HFU}$.The depth to the conductivity increase corresponding to the mineralogic phase transition in the upper mantle is increased with high temperatures. Such a conductive zone shows that the maximum temperature difference between stable platform areas and active zones is about 1000°C.     

Internal constitution and evolution of the moon

The composition, structure and evolution of the moon's interior are narrowly constrained by a large assortment of physical and chemical data. Models of the thermal evolution of the moon that fit the chronology of igneous activity on the lunar surface, the stress history of the lunar lithosphere implied by the presence of mascons, and the surface concentrations of radioactive elements, involve extensive differentiation early in lunar history. This differentiation may be the result of rapid accretion and large-scale melting or of primary chemical layering during accretion; differences in present-day temperatures for these two possibilities are significant only in the inner 1000 km of the moon and may not be resolvable. If the Apollo 15 heat-flow result is representative of the moon, the average uranium concentration in the moon is 0.05–0.08 p.p.m.Density models for the moon, including the effects of temperature and pressure, can be made to satisfy the mass and moment of inertia of the moon and the presence of a low-density crust inferred from seismic refraction studies only if the lunar mantle is chemically or mineralogically inhomogeneous. The upper mantle must exceed the density of the lower mantle at similar conditions by at least 5%. The average mantle density is that of a pyroxenite or olivine pyroxenite, though the density of the upper mantle may exceed 3.5 g/cm³. The density of the lower mantle is less than that of the combined crust and upper mantle at similar temperature and pressure, thus reinforcing arguments for early moon-wide differentiation of both major and minor elements. The suggested density inversion is gravitationally unstable and implies stresses in the mantle 2–5 times those associated with the lunar gravitational field, a difficulty that can be explained or avoided by: (1) adopting lower values for the moment of inertia and/or crustal thickness, or (2) postulating that the strength of the lower mantle increases with depth or with time, either of which is possible for certain combinations of composition and thermal evolution.A small iron-rich core in the moon cannot be excluded by the moon's mass and moment of inertia. If such a core were molten at the time lunar surface rocks acquired remanent magnetization, then thermal-history models with initially cold interiors strongly depleted in radioactive heat sources as a primary accretional feature must be excluded. Further, the presence of ～||pre|40 K in a FeFeS core could significantly alter the thermal evolution and estimated present-day temperatures of the deep lunar interior.     

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.     

Global electromagnetic induction in the Moon and planets

A summary of experiments and analyses concerning electromagnetic induction in the Moon and other extraterrestrial bodies is presented. Magnetic step-transient measurements made on the lunar dark side show the eddy current response to be the dominant induction mode of the Moon. Analysis of the poloidal field decay of the eddy currents has yielded a range of monotonic conductivity profiles for the lunar interior: the conductivity rises from 3·10^?4 mho/m at a depth of 170 km to 10^?2 mho/m at 1000 km depth. The static magnetization field induction has been measured and the whole-Moon relative magnetic permeability has been calculated to be $\text{μ}\text{μ}{}_0\text{= 1.01 ± 0.06}$. The remanent magnetic fields, measured at Apollo landing sites, range from 3 to 327 γ. Simultaneous magnetometer and solar wind spectrometer measurements show that the 38-γ remanent field at the Apollo 12 site is compressed to 54 γ by a solar wind pressure increase of 7·10^?8 dyn/cm². The solar wind confines the induced lunar poloidal field; the field is compressed to the surface on the lunar subsolar side and extends out into a cylindrical cavity on the lunar antisolar side. This solar wind confinement is modeled in the laboratory by a magnetic dipole enclosed in a superconducting lead cylinder; results show that the induced poloidal field geometry is modified in a manner similar to that measured on the Moon. Induction concepts developed for the Moon are extended to estimate the electromagnetic response of other bodies in the solar system.     

Observations of multiple core reflections of the PnKP and SnKP type and regional variations at the base of the mantle

Seismological results interpreted as evidence for large inhomogeneities near the base of the Earth's mantle below Hawaii have recently been published. It is possible to place constraints on the magnitude of such heterogeneities by identifying seismic phases multiply reflected within the Earth's core. The value of such a simple technique is illustrated by using array recordings of P and S5KP waves that have traversed the bottom of the mantle beneath Hawaii to show that there is no clear evidence for the unusual physical properties attributed to this region of the Earth. Identification of the phase S7KP is also reported.     

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.     

Core motions,electromagnetic core-mantle coupling and variations in the earth's rotation: New constraints from geomagnetic secular variation impulses

Recently observed secular acceleration impulses (SAI) of the geomagnetic field are interpreted in terms of organized motions of the outer core layers. Such motions have planetary dimensions (5000 km) and a large amplitude (3 × 10^?4 m s^?1) and are established in very short times (less than one year). The correlation of SAI observed in the Northern Hemisphere with minima in the Earth's rotation rate (around 1840, 1905 and 1970) is shown to be consistent with a simple model involving electromagnetic coupling of the weakly conducting (of the order of 100 ω^?1 m^?1) mantle, of a coherent outer core layer (thickness 100 to a few hundred kilometres) and of the rest of the core. The magnitude of the torque which acts suddenly on both parts of the core at the time of the impulses is estimated.     

Seismic codas on the Earth and the Moon: a comparison

The phenomenon of the seismic coda, which is composed of seismic energy delayed by scattering, is seen on both the Earth and the Moon. On the Moon the scattered coda is very large relative to body wave arrivals with a delay of the time of maximum energy, whereas on Earth scattered codas are relatively small and show no delay of the energy maximum. In both cases the form of the coda is controlled by three distance scales, the mean free path L, which is the average distance seismic energy travels before it is scattered, the attenuation distance $\text{x1}$, which is the average distance seismic energy travels before it is attenuated, and the source-receiver distance R. Two coda models are discussed based on these parameters; a strong scattering (diffusion) model, and a weak scattering (single scattering) model. A discussion of the diffusion scattering model indicates that if $\text{x1/L ? 1}$, diffusion scattering is an appropriate model, but if $\text{x1/L ? 1}$, single scattering is the appropriate model, within the appropriate range of R. A survey of the literature indicates that for the frequency range 0.5–10 Hz, diffusion scattering is important in lunar codas, but for the frequency range 1–25 Hz single scattering is important in terrestrial codas. Another important effect of attenuation is the elimination of scattering paths much longer than $\text{x1}$. On the Moon, this means that seismic energy in the coda can only propagate directly in the near-surface strong scattering zone between surface sources and the seismometer for source-seismometer separations of the order of $\text{(x1L)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$; otherwise, scattering is limited to regions near the source and the receiver. On Earth, this effect probably prevents multiple scattering.     

Inferences on the lunar temperature from gravity,stress state and flow laws

Inferences on the lunar temperature regime are made from the inversion of gravity for density anomalies and the stress-state of the Moon's interior, and by comparing these results with flow laws and estimates of likely strain-rates.The nature of the spectrum of the lunar gravitational potential indicates that the density anomalies giving rise to the potential are mainly of near-surface origon. The average stress-differences in the lunar mantle required to support these density anomalies are of the order of a few tens of bars and have persisted for more than 3 · 10⁹ years. If current flow laws for dry olivine can be extrapolated to the conditions of the lunar mantle, and the selenotherms based on electrical conductivity models are valid, the strain rates are too high to explain the preservation of the lateral near-surface density anomalies. We suggest that the present temperatures in the Moon are relatively low, of the order of 800°C or less, at a depth of about 300 km. This compares with 1100°C based on electrical conductivity models and is near the lower limit predicted by Keihm and Langseth (1977) from lunar heat-flow observations.     

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.