Upwelling flow in Earth's mantle and sea-floor spreading

Upwelling flows in the Earth's mantle are accompanied by mass, momentum and energy transports from deep to upper layers. Those flows beneath the mid-ocean ridges give rise to sea-floor spreading. Mantle plumes, on the other hand, cause hot spots to be formed on the Earth's surface. Using the basic equations of fluid dynamics, temperature and velocity distributions in two-dimensional upwelling and cylindrical plumes can be obtained by an integral-relation method. Then the mass, momentum and energy transported to the lithosphere by these upwelling flows can readily be calculated. Based on those results we can more thoroughly discuss problems of plate dynamics, such as the driving mechanism of plate motion, the causes of formation of rift valleys over mid-ocean ridges, and the effect of mantle plumes on sea-floor spreading.     

On the activation volume for creep and its variation with depth in the Earth's lower mantle

This activation volume ΔV for creep may be derived from Keyes's elastic strain energy model or from Weertman's empirical relationship between viscosity and the melting temperature. These formulations are shown to be equivalent if the anharmonic Grüneisen parameters γ of all acoustic modes are equal and if the pressure dependence of the melting temperature follows Lindemann's law, both of which assumptions are valid for the close-packed mineral structure of the lower mantle. The pressure derivative of ΔV depends only on the bulk modulus and the acoustic γ, both of which are directly available from seismic models. Using the data of Brown and Shankland, we show that ΔV decreases by almost 50% between the top and the bottom of the lower mantle, which makes it easier to maintain a constant viscosity in this region. The isoviscous temperature profile can be adiabatic in the deep lower mantle only below 1700 km depth; it is super-adiabatic in the top 1000 km of the lower mantle.     

Thermal regime of the Earth's core

The outer core is assumed to consist of iron and sulfur, with a small amount of potassium that generates heat by radioactive decay of sim||pre|40 K. Two cases are considered, corresponding respectively to a high rate of heat production (Q = 2 · 10¹² cal./sec, about 0.1% K), and to a low rate (Q = 2 · 10¹¹ cal./sec). The temperature at a depth of 2800 km in the mantle is taken to be 3300°K (Wang, 1972). The temperature T_c at the core-mantle boundary depends on whether or not a density gradient in the lowermost layer D″ of the mantle prevents convection in that layer. In the first case, and for high Q, T_c = 4500–5000°K. In the second case, or for low Q, T_c ≈ 3500°K.The heat-conduction equation is used to calculate the temperature T_i at the inner-core boundary in the absence of convection. For high Q, T_i ? T_c ≈ 1600°K; for low Q, T_i ? T_c ≈ 160°K. Corresponding temperature gradients at r = r_c and r = r_i are listed in Table I.The adiabatic gradient at the top of the core is calculated by the method of Stewart (1970). It strongly depends on the parameters (ρ₀, c₀, γ₀, etc.) that characterize core material at low pressure. Stewart has drawn graphs that allow the selection of sets of parameters that are consistent with seismic velocities and a given density distribution in the core. Some acceptable sets of parameters are listed in Table II. Many sets yield temperatures T_c in the range 3500–5000°K; some give an adiabatic gradient steeper than the conductive gradient and are compatible with convection; others do not. Since properties of FeS melts remain unknown, there is at present no way of selecting any set in preference to another.Properties of the FeS system at low pressure suggest the possible appearance of immiscibility at high temperature in liquids of low sulfur content; accordingly, the inner-core boundary is thought to represent equilibrium between a solid (FeNi) inner core and a liquid layer containing only a small amount of sulfur; layer F in turn is in equilibrium with another liquid (forming layer E) containing more sulfur, and slightly less dense, than F. The temperature T_i at the inner-core boundary is about 6000–6500°K for high Q and T_c ≈ 4500–5000°K. It is consistent with Alder's (1966) and Leppaluoto's (1972) estimates of the melting point of iron at 3.3 Mbar, but not with that of Higgins and Kennedy (1971).     

Earth's magnetic field modelling and earth's structure and evolution

The advantages of the approximation of the Earth's magnetic field by means of the field of the so-called natural magnetic sources are discussed. The shifting of these natural magnetic sources, determined for different epochs, is used to forecast the Earth's magnetic field and to draw conclusions about the motion of the corresponding part of the Earth. On the basis of the representation of the Earth's magnetic field from several past geological epochs as a field of one optimum dipole a new theory about the Earth's evolution is proposed.     

An asymptotic solution of hydrodynamic equations at the mid-plane of a plume in the Earth's mantle

From the partial differential equations of hydrodynamics governing the movements in the Earth's mantle of a Newtonian fluid with a pressure- and temperature-dependent viscosity, considering the bilateral symmetry of velocity and temperature distributions at the mid-plane of the plume, an analytical solution of the governing equations near the mid-plane of the plume was found by the method of asymptotic analysis. The vertical distribution of the upward velocity, viscosity and temperature at the mid-plane, and the temperature excess at the centre of the plume above the ambient mantle temperature were then calculated for two sets of Newtonian rheological parameters. The results obtained show that the temperature at the mid-plane and the temperature excess are nearly independent of the rheological parameters. The upward velocity at the mid-plane, however, is strongly dependent on the rheological parameters.     

Array analysis of lateral inhomogeneities in the deep mantle

P-wave slowness and azimuth anomalies at LASA are critically dependent on the array configuration. This agrees with the interpretation that the anomalies arise by scattering at small-scale randomly distributed inhomogeneities in the crust and uppermost mantle beneath the array. Of particular importance is the result that numerous configurations can be chosen which yield ( dT/dΔ, ?) anomalies which are inconsistent with recent interpretations including a lateral inhomogeneity at the base of the mantle beneath Hawaii. Also configurations giving ( dT/dΔ, ?) anomalies inconsistent with the existence of mantle plumes under Iceland or the Galapagos Islands are found.     

Hafnium/rare earth element fractionation in the sedimentary system and crustal recycling into the Earth's mantle

Among long-lived radioactive parent-daughter element pairs, the ratio Lu/Hf is strongly fractionated relative to constant Sm/Nd in the Earth's sedimentary system. This is caused by high resistance to chemical weathering of the mineral zircon (Zr,Hf)SiO₄. Zircon-bearing sandy sediments on and near continents have very low Lu/Hf, while deep-sea clays have up to three times the chondritic Lu/Hf ratio. Turbidity currents mechanically carry the low-Lu/Hf sandy material onto the ocean floor. The results are important for the crust-to-mantle recycling discussion, where most recycled materials would be subducted oceanic sediments. Such sediment should be capable of explaining the HfNd mantle isotopic variation by mixing with peridotite, but in fact any average pelagic sediment has Nd/Hf and Lu/Hf too high to allow mixing curves to pass through the Hf/Nd isotopic array. The array could only be reproduced by subduction of turbidite sandstone with pelagic sediment in the approximate ratio 1.2 to 1, and by maintaining a good mixture between the two components. At least today, turbidites are available for subduction only at locations quite different and distant from those where pelagic sediments may be recycled; furthermore, mantle isotopic variation shows that the mantle often cannot mix itself well enough to homogenize these widely-separated sedimentary components to the degree required. The Lu/Hf fractionations place a severe restriction on the ability of recycled sediments to explain mantle isotopic patterns.     

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.     

The distribution of thermal diffusivity in the Earth's crust

Data in the literature and additional measurements on the thermal diffusivities of granites, granulites and ultrabasic rocks at temperatures up to 1000 K and pressures to 2 GPa, have been used to propose a new model for thermal diffusivity distribution in the crust and upper mantle.The laboratory measurements were made using a pulse method or the Angstroem method with cylindrical heat flow. After making particular assumptions about the pressure and temperature distribution within the top 60 km the pressure and temperature dependencies of diffusivity were transformed into a depth dependence.The model is characterised by a continuous decrease of diffusivity to a depth of ～30 km where there is a small but rapid increase to a nearly constant value of 7.3 × 10^?3 cm² s^?1.     

Noble gases in the earth's mantle

The abundances and isotopic compositions of noble gases in two samples from ultramafic xenoliths in alkali basalt, a young kaersutitic amphibole separated from a peridotite xenolith from Dish Hill, California and an ancient whole-rock lherzolite xenolith from Baja California, are reported and compared with the results of analyses on other mantle samples. In addition to previously recognized excesses of ³He and ¹²⁹Xe, our results indicate that ambient gases in the mantle show a general enrichment of the lighter-mass nonradiogenic isotopes of Ar, Kr and Xe, and Ar with ⁴⁰Ar/³⁶Ar = 3 · 10².     

The influence of extrinsic pressure changes on the Earth's dynamo

Reversals of the Earth's magnetic field have been claimed to correlate with ice ages, tectonic events and falls of tectites. A physical mechanism is needed to relate reversals with the other events before these correlations can be taken seriously. One possible connection lies through changes in pressure in the core. If events high up in the mantle were to lead to changes in core pressure, this would affect the rate of freezing of the liquid core and modify the power supplied to the dynamo. A sufficiently large modification could set off a reversal or perhaps change the mode of operation of the dynamo from a non-reversing to a reversing state.The model of Gubbins et al., allows a quantitative calculation to be made for the effect of a pressure change on the energy release. Any sufficiently sudden pressure change would change the power, but it seems unrealistic to consider less than a 1000 year time scale. Relaxation of shear forces in the mantle, overturning of core fluid, and changes in magnetic fields all take place on about this time scale. According to the model, a pressure change of 0.1 bar over a 1000 years could change the power supply drastically.A continuous process of mantle differentiation leading to the formation of the upper mantle from an initially homogeneous mantle can only provide 10% of the required pressure change, but the effect cannot be ruled out as a power source for the dynamo because uncertainties in the calculations can amount to at least an order of magnitude. The other effects produce changes of up to 1% in the power supply, which may be sufficient to alter the characteristics of the dynamo and produce reversals or a change in reversal behaviour. Further speculation must await a better understanding of the dynamics of reversals, and of mantle processes.     

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.     

The scattering of seismic waves through a crust and upper mantle with random lateral and vertical inhomogeneities

The propagation of seismic waves through Earth models with slightly random lateral and vertical inhomogeneities superimposed on one composed of layers with vertical velocity gradients was investigated. The maximum deviation of velocity from a mean value at a given depth and a correlation distance derived from a two-dimensional smoothing filter were two parameters used to vary the amplitude and size of the velocity anomalies. The resulting models show short discontinuous reflectors scattered about at various depths throughout the model, and are thus in agreement with many deep seismic reflection experiments. On the other hand numerical experiments using ray-tracing techniques showed that the effect of the lateral and vertical velocity anomalies is to scatter the energy, and break up the continuous travel-time lines from vertical gradient models into travel-time segments with different slopes similar to those observed in many long range seismic refraction experiments, and to those resulting from layering effects in the media. Many of the numerical experiments which modelled the random crust produced a Pg segment and a P* segment with an apparent Conrad discontinuity at a depth of 10–20 km, this apparent depth being related to the correlation distance.When a seismic wave propagates through a heterogeneous Earth the amount of its energy which is converted into scattered energy will be a function of the inhomogeneous characteristics of the medium through which it has passed. If a ray passes through a homogeneous Earth the energy arriving at an array station should be relatively coherent whereas if the ray encounters lateral and vertical inhomogeneities its energy will be incoherent and much more complex. A series of coherency measurements done on array recordings of earthquakes at various distances showed that large lateral and vertical variations in complexity exist for different ray paths through the Earth with the region below the 650 km discontinuity in the mantle tending to be much simpler than the region just below the lithosphere.     

Dissipation of the rare gases contained in the primordial Earth's atmosphere

If the Earth was formed by accumulation of rocky bodies in the presence of the gases of the primordial solar nebula, the Earth at this formation stage was surrounded by a massive primordial atmosphere (of about 1 × 10²⁶ g) composed mainly of H₂ and He. We suppose that the H₂ and He escaped from the Earth, owing to the effects of strong solar wind and EUV radiation, in stages after the solar nebula itself dissipated into the outer space.The primordial atmosphere also contained the rare gases Ne, Ar, Kr and Xe whose amounts were much greater than those contained in the present Earth's atmosphere. Thus, we have studied in this paper the dissipation of these rare gases due to the drag effect of outflowing hydrogen molecules. By means of the two-component gas kinetic theory and under the assumption of spherically symmetric flow, we have found that the outflow velocity of each rare gas relative to that of hydrogen is expressed in terms of only two parameters — the rate of hydrogen mass flow across the spherical surface under consideration and the temperature at this surface. According to this result, the rare gases were dissipated below the levels of their contents in the present atmosphere, when the mass loss rate of hydrogen was much greater than 1 × 10¹⁷ g/yr throughout the stages where the atmospheric mass decreased from 1 × 10²⁶ g to 4 × 10¹⁹ g.     

Dissolution of the primordial rare gases into the molten Earth's material

We have shown in a previous paper that, if the primordial solar nebula existed when the Earth was formed, the Earth was once surrounded by a dense and massive primordial atmosphere, whose temperature and pressure were about 4000 K and 900 atm, respectively, at the bottom. We suppose that this hydrogen-rich atmosphere escaped from the Earth after the solar nebula itself disappeared, both phenomena probably being due to the effect of strong solar wind and radiation.Using the results of our previous and new calculations on the structure of the primordial atmosphere, we have investigated the amount of dissolution of the rare gases, which were contained in the primordial atmosphere, into the molten Earth's material.The amount of the dissolved rare gases is found to be strongly dependent on the grain opacity of the atmosphere, i.e., on the amount of fine grains. However, their isotopic ratios and relative abundance are independent of the opacity and approximately equal to those in the primordial solar nebula, that is, to the present solar values. Especially, the dissolved neon is expected to have remained in the present mantle. Therefore, if a considerable amount of neon with nearly the solar isotopic ratio is discovered in present mantle material, this offers direct evidence for the proposition that the proto-Earth was once surrounded by the primordial atmosphere.     

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.     

Earth's melting due to the blanketing effect of the primordial dense atmosphere

When the proto-Earth was growing by the accretion of planetesimals and its mass became greater than about 0.1 M_E, where M_E is the present Earth's mass, an appreciable amount of gas of the surrounding solar nebula was attracted towards the proto-Earth to form an optically thick, dense atmosphere. We have studied the structure of this primordial atmosphere under the assumptions that (1) it is spherically symmetric and in hydrostatic equilibrium, and (2) the net energy outflow (i.e., the luminosity) is constant throughout the atmosphere and is given by GMM/R with M = M/10⁶yr or M/10⁷yr where M and R are the mass and the radius of the proto-Earth, respectively.The results of calculations show that the temperature at the bottom of the atmosphere, namely, at the surface of the proto-Earth increases greatly with the mass of the proto-Earth and it is about 1500°K for M = 0.25 M_E. This high temperature is due to the blanketing effect of the opaque atmosphere. Thus, as long as the primordial solar nebula was existing, the surface temperature of the proto-Earth was kept high enough to melt most of the materials and, hence, the melted iron sedimented towards the center to form the Earth's core.     

How strong is the magnetic field in the Earth's liquid core?

A crucial step in the investigation of the energetics of motions in the Earth's core and the generation of the geomagnetic field by the hydromagnetic dynamo process is the estimation of the average strength $\text{B}$ of the magnetic field B = B_p + B_T in the core. Owing to the probability that the toroidal field B_T in the core, which has no radial component, is a good deal stronger than the poloidal field B_p, direct downward extrapolation of the surface field to the core-mantle interface gives no more than an extreme lower limit to $\text{B}$. This paper outlines the indirect methods by which $\text{B}$ can be estimated, arguing that $\text{B}$ is probably about 10^?2 T (100 Γ) but might be as low as 10^?3 T (10 Γ) or as high as 5 × 10^?2 T (500 Γ).     

Evidence for long-term asymmetries in the Earth's magnetic field and possible implications for dynamo theories

Paleomagnetic data indicate that there is a north-south asymmetry in the time-averaged magnetic field and that there are small but significant differences between the normal and reverse polarity states. The geographical variation is most likely due to spatial variation in the boundary conditions at the core-mantle interface. The difference in the magnetic fields of the reverse and normal polarity states can be modeled in terms of a “standing field”. The paleomagnetic data are insufficient to determine whether or not this “standing field” is of core origin. However, consideration of mechanisms, including thermoelectric currents, indicates that there probably are important differences in core processes between the two polarity states. At first glance this interpretation is difficult to reconcile with the fact that the magnetic induction equation is antisymmetric with respect to the magnetic field. A way around this problem is the possibility that only certain transitions are allowed between acceptable eigenstates in dynamo models of the Earth's magnetic field.